An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （...An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.展开更多
Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficien...Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.展开更多
A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved ...A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved boundaries, this method can handle the free-surface problems with complex geometries accurately and directly, without introducing any complicated boundary treatment or artificial diffusion. The method solves the volume transport equation geometrically through the Modified Lagrangian-Eulerian Re-map （MLER） method, which is applied to advective fluid volumes. Moreover, the PLIC method is adopted to give a second-order reconstructed interface approximation. To validate this method, two advection tests were performed for the establishment of the accuracy and convergence rate of the solutions. Numerical results for these complex tests provide convincing evidence for the excellent solution quality and fidelity of the method.展开更多
The paper presents a finite volume numerical method universally applicable for solving both linear and nonlinear aeroacoustics problems on arbitrary unstructured meshes. It is based on the vertexcentered multi-paramet...The paper presents a finite volume numerical method universally applicable for solving both linear and nonlinear aeroacoustics problems on arbitrary unstructured meshes. It is based on the vertexcentered multi-parameter scheme offering up to the 6th accuracy order achieved on the Cartesian meshes. An adaptive dissipation is added for the numerical treatment of possible discontinuities. The scheme properties are studied on a series of test cases, its efficiency is demonstrated at simulating the noise suppression in resonance-type liners.展开更多
A high-order Lagrangian cell-centered conservative gas dynamics scheme is presented on unstructured meshes. A high-order piecewise pressure of the cell is intro- duced. With the high-order piecewise pressure of the ce...A high-order Lagrangian cell-centered conservative gas dynamics scheme is presented on unstructured meshes. A high-order piecewise pressure of the cell is intro- duced. With the high-order piecewise pressure of the cell, the high-order spatial discretiza- tion fluxes are constructed. The time discretization of the spatial fluxes is performed by means of the Taylor expansions of the spatial discretization fluxes. The vertex velocities are evaluated in a consistent manner due to an original solver located at the nodes by means of momentum conservation. Many numerical tests are presented to demonstrate the robustness and the accuracy of the scheme.展开更多
Based on a new second-order neutron transport equation, self-adjoint angular flux (SAAF) equation, the spherical harmonics (PN) method for neutron transport equation on unstructured-meshes is derived. The spherical ha...Based on a new second-order neutron transport equation, self-adjoint angular flux (SAAF) equation, the spherical harmonics (PN) method for neutron transport equation on unstructured-meshes is derived. The spherical harmonics function is used to expand the angular flux. A set of differential equations about the spatial variable, which are coupled with each other, can be obtained. They are solved iteratively by using the finite element method on un- structured-meshes. A two-dimension transport calculation program is coded according to the model. The numerical results of some benchmark problems demonstrate that this method can give high precision results and avoid the ray effect very well.展开更多
To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any compl...To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any complicated boundary treatment or artificial diffusion, this method treated curved boundaries directly by utilizing the inherent merit of unstructured mesh in fitting curves. The PLIC (Piecewise Linear Interface Calculation) method was adopted to obtain a second-order accurate linearized reconstruction approximation and the MLER (Modified Lagrangian-Eulerian Re-map) method was introduced to advect fluid volumes on unstructured mesh. Moreover, an analytical relation for the interface’s line constant vs. the volume clipped by the interface was developed so as to improve the method’s efficiency. To validate this method, a comprehensive series of large straining advection tests were performed. Numerical results provide convincing evidences for the method’s high volume conservative accuracy and second-order shape error convergence rate. Also, a dramatic improvement on computational accuracy over its unstructured triangular mesh counterpart is checked.展开更多
A discontinuity-capturing scheme of finite element method(FEM)is proposed.The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unstea...A discontinuity-capturing scheme of finite element method(FEM)is proposed.The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unsteady flows,which exhibits the capability of capturing the shock waves and/or thin shear layers accurately in an unsteady viscous flow at high Reynolds number. In particular,a new testing variable,i.e.,the disturbed kinetic energy E,is suggested and used in the adaptive mesh computation,which is universally applicable to the capturing of both shock waves and shear layers in the inviscid flow and viscous flow at high Reynolds number.Based on several calculated examples,this approach has been proved to be effective and efficient for the calculations of compressible and incompressible flows.展开更多
Higher-order displacement-based finite element methods are useful for simulating bending problems and potentially addressing mesh-locking associated with nearly-incompressible elasticity,yet are computationally expens...Higher-order displacement-based finite element methods are useful for simulating bending problems and potentially addressing mesh-locking associated with nearly-incompressible elasticity,yet are computationally expensive.To address the computational expense,the paper presents a matrix-free,displacement-based,higher-order,hexahedral finite element implementation of compressible and nearly-compressible(ν→0.5)linear isotropic elasticity at small strain with p-multigrid preconditioning.The cost,solve time,and scalability of the implementation with respect to strain energy error are investigated for polynomial order p=1,2,3,4 for compressible elasticity,and p=2,3,4 for nearly-incompressible elasticity,on different number of CPU cores for a tube bending problem.In the context of this matrix-free implementation,higher-order polynomials(p=3,4)generally are faster in achieving better accuracy in the solution than lower-order polynomials(p=1,2).However,for a beam bending simulation with stress concentration(singularity),it is demonstrated that higher-order finite elements do not improve the spatial order of convergence,even though accuracy is improved.展开更多
Wavefield extrapolation is critical in reverse time migration(RTM).The finite diff erence method is primarily used to achieve wavefi eld extrapolation in case of the RTM imaging of tunnels.However,complex tunnel model...Wavefield extrapolation is critical in reverse time migration(RTM).The finite diff erence method is primarily used to achieve wavefi eld extrapolation in case of the RTM imaging of tunnels.However,complex tunnel models,including those for karsts and fault fracture zones,are constructed using regular grids with straight curves,which can cause numerical dispersion and reduce the imaging accuracy.In this study,wavefi eld extrapolation was conducted for tunnel RTM using the finite element method,wherein an unstructured mesh was considered to be the body-fi tted partition in a complex model.Further,a Poynting vector calculation equation suitable for the unstructured mesh considered in the fi nite element method was established to suppress the interference owing to low-frequency noise.The tunnel space was considered during wavefi eld extrapolation to suppress the mirror artifacts based on the fl exibility of mesh generation.Finally,the infl uence of the survey layouts(one and two sidewalls)on the tunnel imaging results was investigated.The RTM results obtained for a simple tunnel model with an inclined interface demonstrate that the method based on unstructured meshes can effectively suppress the low-frequency noise and mirror artifacts,obtaining clear imaging results.Furthermore,the two-sidewall tunnel survey layout can be used to accurately obtain the real position of the inclined interface ahead of the tunnel face.The complex tunnel numerical modeling and actual data migration results denote the eff ectiveness of the fi nite element method in which an unstructured mesh is used.展开更多
In this paper, the accuracy of Chang＇s unstructured space-time conservation element and solution element （CE/SE） scheme is analysed for the first time. Based on a redefinition of conservation elements and solution ...In this paper, the accuracy of Chang＇s unstructured space-time conservation element and solution element （CE/SE） scheme is analysed for the first time. Based on a redefinition of conservation elements and solution elements, an improved two-dimensional （2D） unstructured CE/SE scheme with an adjustable parameter β is proposed to accurately capture shock waves. The new scheme can be applied to any type of grid without special treatnmnt. Compared with Chang＇s original parameter a, larger/5 dose not cost extra computational resources. Numerical tests reveal that the new scheme is not only clear in physical concept, compact and highly accurate but also more capable of capturing shock waves than the popular fifth-order accurate weighted essentially non-oscillatory scheme.展开更多
Efficient and robust solution strategies are developed for discontinuous Galerkin （DG） discretization of the Navier-Stokes （NS） and Reynolds-averaged NS （RANS） equations on structured/unstructured hybrid meshes....Efficient and robust solution strategies are developed for discontinuous Galerkin （DG） discretization of the Navier-Stokes （NS） and Reynolds-averaged NS （RANS） equations on structured/unstructured hybrid meshes. A novel line-implicit scheme is devised and implemented to reduce the memory gain and improve the computational eificiency for highly anisotropic meshes. A simple and effective technique to use the mod- ified Baldwin-Lomax （BL） model on the unstructured meshes for the DC methods is proposed. The compact Hermite weighted essentially non-oscillatory （HWENO） limiters are also investigated for the hybrid meshes to treat solution discontinuities. A variety of compressible viscous flows are performed to examine the capability of the present high- order DG solver. Numerical results indicate that the designed line-implicit algorithms exhibit weak dependence on the cell aspect-ratio as well as the discretization order. The accuracy and robustness of the proposed approaches are demonstrated by capturing com- plex flow structures and giving reliable predictions of benchmark turbulent problems.展开更多
文摘An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.
基金National High Technology Research and Development Program(863 Program)(No.2006AA06Z105,2007AA06Z134)
文摘Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.
文摘A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved boundaries, this method can handle the free-surface problems with complex geometries accurately and directly, without introducing any complicated boundary treatment or artificial diffusion. The method solves the volume transport equation geometrically through the Modified Lagrangian-Eulerian Re-map （MLER） method, which is applied to advective fluid volumes. Moreover, the PLIC method is adopted to give a second-order reconstructed interface approximation. To validate this method, two advection tests were performed for the establishment of the accuracy and convergence rate of the solutions. Numerical results for these complex tests provide convincing evidence for the excellent solution quality and fidelity of the method.
基金Russian Foundation of Basic Research(No. 04-01-08034, 06-01-00293-a)
文摘The paper presents a finite volume numerical method universally applicable for solving both linear and nonlinear aeroacoustics problems on arbitrary unstructured meshes. It is based on the vertexcentered multi-parameter scheme offering up to the 6th accuracy order achieved on the Cartesian meshes. An adaptive dissipation is added for the numerical treatment of possible discontinuities. The scheme properties are studied on a series of test cases, its efficiency is demonstrated at simulating the noise suppression in resonance-type liners.
基金supported by the National Natural Science Foundation of China(Nos.11172050,11372051,and 11001027)
文摘A high-order Lagrangian cell-centered conservative gas dynamics scheme is presented on unstructured meshes. A high-order piecewise pressure of the cell is intro- duced. With the high-order piecewise pressure of the cell, the high-order spatial discretiza- tion fluxes are constructed. The time discretization of the spatial fluxes is performed by means of the Taylor expansions of the spatial discretization fluxes. The vertex velocities are evaluated in a consistent manner due to an original solver located at the nodes by means of momentum conservation. Many numerical tests are presented to demonstrate the robustness and the accuracy of the scheme.
基金Supported by pre-research fund of State Key Laboratory (51479080201 JW0802)
文摘Based on a new second-order neutron transport equation, self-adjoint angular flux (SAAF) equation, the spherical harmonics (PN) method for neutron transport equation on unstructured-meshes is derived. The spherical harmonics function is used to expand the angular flux. A set of differential equations about the spatial variable, which are coupled with each other, can be obtained. They are solved iteratively by using the finite element method on un- structured-meshes. A two-dimension transport calculation program is coded according to the model. The numerical results of some benchmark problems demonstrate that this method can give high precision results and avoid the ray effect very well.
基金the National Natural Science Foundation ofChina under Grant No. 50779043, 50779045
文摘To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any complicated boundary treatment or artificial diffusion, this method treated curved boundaries directly by utilizing the inherent merit of unstructured mesh in fitting curves. The PLIC (Piecewise Linear Interface Calculation) method was adopted to obtain a second-order accurate linearized reconstruction approximation and the MLER (Modified Lagrangian-Eulerian Re-map) method was introduced to advect fluid volumes on unstructured mesh. Moreover, an analytical relation for the interface’s line constant vs. the volume clipped by the interface was developed so as to improve the method’s efficiency. To validate this method, a comprehensive series of large straining advection tests were performed. Numerical results provide convincing evidences for the method’s high volume conservative accuracy and second-order shape error convergence rate. Also, a dramatic improvement on computational accuracy over its unstructured triangular mesh counterpart is checked.
基金The project supported by the National Natural Science Foundation of China (10125210),the Hundred-Talent Programme of the Chinese Academy of Sciences and the Innovation Project of the Chinese Academy of Sciences (KJCX-SW-L04,KJCX2-SW-L2)
文摘A discontinuity-capturing scheme of finite element method(FEM)is proposed.The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unsteady flows,which exhibits the capability of capturing the shock waves and/or thin shear layers accurately in an unsteady viscous flow at high Reynolds number. In particular,a new testing variable,i.e.,the disturbed kinetic energy E,is suggested and used in the adaptive mesh computation,which is universally applicable to the capturing of both shock waves and shear layers in the inviscid flow and viscous flow at high Reynolds number.Based on several calculated examples,this approach has been proved to be effective and efficient for the calculations of compressible and incompressible flows.
基金The research relied on computational resources[29]provided by the University of Colorado Boulder Research Computing Group,which is supported by the National1302 CMES,2021,vol.129,no.3 Science Foundation(Awards ACI-1532235 and ACI-1532236)University of Colorado Boulder,and Colorado State University.
文摘Higher-order displacement-based finite element methods are useful for simulating bending problems and potentially addressing mesh-locking associated with nearly-incompressible elasticity,yet are computationally expensive.To address the computational expense,the paper presents a matrix-free,displacement-based,higher-order,hexahedral finite element implementation of compressible and nearly-compressible(ν→0.5)linear isotropic elasticity at small strain with p-multigrid preconditioning.The cost,solve time,and scalability of the implementation with respect to strain energy error are investigated for polynomial order p=1,2,3,4 for compressible elasticity,and p=2,3,4 for nearly-incompressible elasticity,on different number of CPU cores for a tube bending problem.In the context of this matrix-free implementation,higher-order polynomials(p=3,4)generally are faster in achieving better accuracy in the solution than lower-order polynomials(p=1,2).However,for a beam bending simulation with stress concentration(singularity),it is demonstrated that higher-order finite elements do not improve the spatial order of convergence,even though accuracy is improved.
基金supported by the National Natural Science Foundation of China (Nos. 41804145, 41704146)Natural Science Foundation of Hebei Province (No. D2018210168)Project of Hebei Province Higher Educational Science and Technology Program (No.QN2019185)。
文摘Wavefield extrapolation is critical in reverse time migration(RTM).The finite diff erence method is primarily used to achieve wavefi eld extrapolation in case of the RTM imaging of tunnels.However,complex tunnel models,including those for karsts and fault fracture zones,are constructed using regular grids with straight curves,which can cause numerical dispersion and reduce the imaging accuracy.In this study,wavefi eld extrapolation was conducted for tunnel RTM using the finite element method,wherein an unstructured mesh was considered to be the body-fi tted partition in a complex model.Further,a Poynting vector calculation equation suitable for the unstructured mesh considered in the fi nite element method was established to suppress the interference owing to low-frequency noise.The tunnel space was considered during wavefi eld extrapolation to suppress the mirror artifacts based on the fl exibility of mesh generation.Finally,the infl uence of the survey layouts(one and two sidewalls)on the tunnel imaging results was investigated.The RTM results obtained for a simple tunnel model with an inclined interface demonstrate that the method based on unstructured meshes can effectively suppress the low-frequency noise and mirror artifacts,obtaining clear imaging results.Furthermore,the two-sidewall tunnel survey layout can be used to accurately obtain the real position of the inclined interface ahead of the tunnel face.The complex tunnel numerical modeling and actual data migration results denote the eff ectiveness of the fi nite element method in which an unstructured mesh is used.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10732010 and 10972010)
文摘In this paper, the accuracy of Chang＇s unstructured space-time conservation element and solution element （CE/SE） scheme is analysed for the first time. Based on a redefinition of conservation elements and solution elements, an improved two-dimensional （2D） unstructured CE/SE scheme with an adjustable parameter β is proposed to accurately capture shock waves. The new scheme can be applied to any type of grid without special treatnmnt. Compared with Chang＇s original parameter a, larger/5 dose not cost extra computational resources. Numerical tests reveal that the new scheme is not only clear in physical concept, compact and highly accurate but also more capable of capturing shock waves than the popular fifth-order accurate weighted essentially non-oscillatory scheme.
基金Project supported by the National Basic Research Program of China(No.2009CB724104)
文摘Efficient and robust solution strategies are developed for discontinuous Galerkin （DG） discretization of the Navier-Stokes （NS） and Reynolds-averaged NS （RANS） equations on structured/unstructured hybrid meshes. A novel line-implicit scheme is devised and implemented to reduce the memory gain and improve the computational eificiency for highly anisotropic meshes. A simple and effective technique to use the mod- ified Baldwin-Lomax （BL） model on the unstructured meshes for the DC methods is proposed. The compact Hermite weighted essentially non-oscillatory （HWENO） limiters are also investigated for the hybrid meshes to treat solution discontinuities. A variety of compressible viscous flows are performed to examine the capability of the present high- order DG solver. Numerical results indicate that the designed line-implicit algorithms exhibit weak dependence on the cell aspect-ratio as well as the discretization order. The accuracy and robustness of the proposed approaches are demonstrated by capturing com- plex flow structures and giving reliable predictions of benchmark turbulent problems.
