In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.T...Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.This study primarily presents a self-developed 2D ion cyclotron resonance antenna electromagnetic field solver(ICRAEMS)code implemented on the MATLAB platform,which solves the electric field wave equation by using the finite element method,establishing perfectly matched layer(PML)boundary conditions,and post-processing the electromagnetic field data.This code can be utilized to facilitate the design and optimization processes of antennas for ICRF heating technology.Furthermore,this study examines the electric field distribution and power spectrum associated with various antenna phases to investigate how different antenna configurations affect the electromagnetic field propagation and coupling characteristics.展开更多
Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis r...Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.展开更多
For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of ...For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.展开更多
A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines...A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.展开更多
In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique ...In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.展开更多
On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these meth...On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming （bi）linear macroelements or （bi）quadratic elements, and the rotation by conforming （bi）linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.展开更多
In this article,two kinds of expandable parallel finite element methods,based on two-grid discretizations,are given to solve the linear elliptic problems.Compared with the classical local and parallel finite element m...In this article,two kinds of expandable parallel finite element methods,based on two-grid discretizations,are given to solve the linear elliptic problems.Compared with the classical local and parallel finite element methods,there are two attractive features of the methods shown in this article:1)a partition of unity is used to generate a series of local and independent subproblems to guarantee the final approximation globally continuous;2)the computational domain of each local subproblem is contained in a ball with radius of O(H)(H is the coarse mesh parameter),which means methods in this article are more suitable for parallel computing in a large parallel computer system.Some a priori error estimation are obtained and optimal error bounds in both H^1-normal and L^2-normal are derived.Finally,numerical results are reported to test and verify the feasibility and validity of our methods.展开更多
In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
The transient behavior of a semiconductor device consists of a Poisson equation for the electric potential and of two nonlinear parabolic equations for the electron density and hole density. The electric potential equ...The transient behavior of a semiconductor device consists of a Poisson equation for the electric potential and of two nonlinear parabolic equations for the electron density and hole density. The electric potential equation is discretized by a mixed finite element method. The electron and hole density equations are treated by implicit-explicit multistep finite element methods. The schemes are very efficient. The optimal order error estimates both in time and space are derived.展开更多
This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their...This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
This paper is devoted to studying the superconvergence of streamline diffusion finite element methods for convection-diffusion problems. In [8], under the condition that ε ≤ h^2 the optimal finite element error esti...This paper is devoted to studying the superconvergence of streamline diffusion finite element methods for convection-diffusion problems. In [8], under the condition that ε ≤ h^2 the optimal finite element error estimate was obtained in L^2-norm. In the present paper, however, the same error estimate result is gained under the weaker condition that ε≤h.展开更多
The equilibrium equations and the functional for tubing buckling in arbitrary straight wells are derived. The entire buckling process of tubing in deviated wells is analyzed for the first time by utilizing the finite ...The equilibrium equations and the functional for tubing buckling in arbitrary straight wells are derived. The entire buckling process of tubing in deviated wells is analyzed for the first time by utilizing the finite element method. The effects of gravity and torques on the buckling are included in the analyses and the calculated results are well compared with existing solutions. It is shown that the buckling only occurs at the lower portion of the tubing where the axial load is the largest, and the contact force of the well, the bending moment of the tubing and the buckling displacement of this portion vary periodically. The buckling spreads upwards from the bit with the increase of axial load. There is no buckling at the upper portion of the tubing where the bending moment is zero. And the contact force of this section increases only slightly with the increase of the axial load. With the increase of the deviation angle, the length of buckling portion and buckling displacement amplitude decrease, the contact force increases with the increase of load at the upper portion and its amplitude decreases at the lower buckling section, the bending moment remains zero at the upper portion and its amplitude decreases at the lower buckling portion. The buckling displacement increases with the increase of the torque, but the increment is very small.展开更多
We derived and analyzed a new numerical scheme for the coupled Stokes and Darcy problems by using H（div） conforming elements in the entire domain. The approach employs the mixed finite element method for the Darcy e...We derived and analyzed a new numerical scheme for the coupled Stokes and Darcy problems by using H（div） conforming elements in the entire domain. The approach employs the mixed finite element method for the Darcy equations and a stabilized H（div） finite element method for the Stokes equations. Optimal error estimates for the fluid velocity and pressure are derived. The finite element solutions from the new scheme not only feature a full satisfaction of the continuity equation, which is highly demanded in scientific computing, but also satisfy the mass conservation.展开更多
A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations....A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.展开更多
A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface i...A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.展开更多
The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and ...The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass ...Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.展开更多
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
基金Project supported by the National MCF Energy R&D Program(Grant No.2022YFE03190100)the National Natural Science Foundation of China(Grant Nos.12422513,12105035,and U21A20438)the Xiaomi Young Talents Program.
文摘Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.This study primarily presents a self-developed 2D ion cyclotron resonance antenna electromagnetic field solver(ICRAEMS)code implemented on the MATLAB platform,which solves the electric field wave equation by using the finite element method,establishing perfectly matched layer(PML)boundary conditions,and post-processing the electromagnetic field data.This code can be utilized to facilitate the design and optimization processes of antennas for ICRF heating technology.Furthermore,this study examines the electric field distribution and power spectrum associated with various antenna phases to investigate how different antenna configurations affect the electromagnetic field propagation and coupling characteristics.
文摘Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.
基金supported by National Natural Science Foundation of China(11771257)the Shandong Provincial Natural Science Foundation of China(ZR2023YQ002,ZR2023MA007,ZR2021MA004)。
文摘For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.
基金supported by the National Natural Science Foundation of China(11871312,12131014)the Natural Science Foundation of Shandong Province,China(ZR2023MA086)。
文摘A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.
文摘In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.
基金supported by NSFC(11571266,91430106,11171168,11071132)NSFC-RGC(China-Hong Kong)(11661161017)
文摘On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming （bi）linear macroelements or （bi）quadratic elements, and the rotation by conforming （bi）linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.
基金Subsidized by NSFC (11701343)partially supported by NSFC (11571274,11401466)
文摘In this article,two kinds of expandable parallel finite element methods,based on two-grid discretizations,are given to solve the linear elliptic problems.Compared with the classical local and parallel finite element methods,there are two attractive features of the methods shown in this article:1)a partition of unity is used to generate a series of local and independent subproblems to guarantee the final approximation globally continuous;2)the computational domain of each local subproblem is contained in a ball with radius of O(H)(H is the coarse mesh parameter),which means methods in this article are more suitable for parallel computing in a large parallel computer system.Some a priori error estimation are obtained and optimal error bounds in both H^1-normal and L^2-normal are derived.Finally,numerical results are reported to test and verify the feasibility and validity of our methods.
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
文摘The transient behavior of a semiconductor device consists of a Poisson equation for the electric potential and of two nonlinear parabolic equations for the electron density and hole density. The electric potential equation is discretized by a mixed finite element method. The electron and hole density equations are treated by implicit-explicit multistep finite element methods. The schemes are very efficient. The optimal order error estimates both in time and space are derived.
文摘This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
基金Supported by the National Natural Science Foundation of China(10471103)
文摘This paper is devoted to studying the superconvergence of streamline diffusion finite element methods for convection-diffusion problems. In [8], under the condition that ε ≤ h^2 the optimal finite element error estimate was obtained in L^2-norm. In the present paper, however, the same error estimate result is gained under the weaker condition that ε≤h.
文摘The equilibrium equations and the functional for tubing buckling in arbitrary straight wells are derived. The entire buckling process of tubing in deviated wells is analyzed for the first time by utilizing the finite element method. The effects of gravity and torques on the buckling are included in the analyses and the calculated results are well compared with existing solutions. It is shown that the buckling only occurs at the lower portion of the tubing where the axial load is the largest, and the contact force of the well, the bending moment of the tubing and the buckling displacement of this portion vary periodically. The buckling spreads upwards from the bit with the increase of axial load. There is no buckling at the upper portion of the tubing where the bending moment is zero. And the contact force of this section increases only slightly with the increase of the axial load. With the increase of the deviation angle, the length of buckling portion and buckling displacement amplitude decrease, the contact force increases with the increase of load at the upper portion and its amplitude decreases at the lower buckling section, the bending moment remains zero at the upper portion and its amplitude decreases at the lower buckling portion. The buckling displacement increases with the increase of the torque, but the increment is very small.
基金The Key Technologies R&D Program ofSichuan Province (No.05GG006-0062)
文摘We derived and analyzed a new numerical scheme for the coupled Stokes and Darcy problems by using H（div） conforming elements in the entire domain. The approach employs the mixed finite element method for the Darcy equations and a stabilized H（div） finite element method for the Stokes equations. Optimal error estimates for the fluid velocity and pressure are derived. The finite element solutions from the new scheme not only feature a full satisfaction of the continuity equation, which is highly demanded in scientific computing, but also satisfy the mass conservation.
基金Supported by the National Natural Science Foundation of China(50976072,51106099,10902070)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)the Science Foundation for the Excellent Youth Scholar of Higher Education of Shanghai(slg09003)~~
文摘A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.
基金supported by the Open Project of Key Laboratory of Aerospace EDLA,CASC(No.EDL19092208)。
文摘A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.
文摘The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
基金supported by the National Natural Science Foundation of China(Grant No.11304160)the Natural Science Foundation of Jiangsu Provincial Higher Education Institutions,China(Grant No.13KJB140008)the Foundation of Nanjing University of Posts and Telecommunications,China(Grant No.NY213018)
文摘Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.
