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.展开更多
This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been pu...This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been published,but very few on those generated from a mesh using triangular elements.The use of triangular elements is aimed at extending the application of the approach to any shape of modeled devices.Basic concepts of the approach are presented in the case of electromagnetic devices.The procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB environment are also included.展开更多
Polygonal finite elements remain an attractive option in finite element analysis due to their flexibility in modelingarbitrary shapes compared to triangles.In this study,a pentagonal membrane element was developed wit...Polygonal finite elements remain an attractive option in finite element analysis due to their flexibility in modelingarbitrary shapes compared to triangles.In this study,a pentagonal membrane element was developed with thestrain approach for the first time.The element possesses invariance,and the equilibrium constraint was appliedto the assumed strain field using corrective coefficients.Inspired by the advancing front technique,a pentagonalmesh was generated,and the mesh quality was enhanced with Laplacian smoothing.The performance of thedeveloped pentagonal element was assessed in a few numerical tests,and the results revealed its suitability inmodeling the bending of beams.Besides,the numerical results are enhanced when pentagonal elements are usedin mesh transitions along boundaries to smoothen curved edges and capture distributed loads.展开更多
A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is d...A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.展开更多
The Wilson finite element method is considered to solve a class of two- dimensional second order elliptic boundary value problems. By using of the particular structure of the element and some new techniques, we obtain...The Wilson finite element method is considered to solve a class of two- dimensional second order elliptic boundary value problems. By using of the particular structure of the element and some new techniques, we obtain the superclose and global superconvergence on anisotropic meshes. Numerical example is also given to confirm our theoretical analysis.展开更多
Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolat...Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method, The solution accuracy is further improved by implementing an adaptive meshing technique to generaie finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems tire: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to moving heat source.展开更多
The main aim of this paper is to study the approximation to viscoelasticity type equations with a Crouzeix-Raviart type nonconforming finite element on the anisotropic meshes. The superclose property of the exact solu...The main aim of this paper is to study the approximation to viscoelasticity type equations with a Crouzeix-Raviart type nonconforming finite element on the anisotropic meshes. The superclose property of the exact solution and the optimal error estimate of its derivative with respect to time are derived by using some novel techniques. Moreover, employing a postprocessing technique, the global superconvergence property for the discretization error of the postprocessed discrete solution to the solution itself is studied.展开更多
The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads...The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads to enhancement of accuracy at coarse meshes, and that the combination parameter plays an important role in the enhancement. As an improvement of conforming bilinear Q(4)-plane element, the combined hybrid method adopted the most convenient quadrilateral displacements-stress mode, i.e.,the mode of compatible isoparametric bilinear displacements and pure constant stresses. By adjusting the combined parameter, the optimized version of the combined hybrid element was obtained and numerical tests indicated that this parameter-adjusted version behaves much better than Q(4)-element and is of high accuracy at coarse meshes. Due to elimination of stress parameters at the elemental level, this combined hybrid version is of the same computational cost as that of Q(4)-element.展开更多
Mapping mesh generation is widely applied in pre-processes of Finite Element Method (FEM). In this study, the basic 3D mapping equations by Lagrange interpolating function are founded. Based these equations, a mappi...Mapping mesh generation is widely applied in pre-processes of Finite Element Method (FEM). In this study, the basic 3D mapping equations by Lagrange interpolating function are founded. Based these equations, a mapping pattern library, which maps essential configurations e.g. line, circle, rotary body, sphere etc. to hexahedral FEM mesh, has been built. Then available FEM mesh will be generated by clipping and assembling the mapped essential objects. Study case illustrates that the proposed method is simple and efficient to generate valid FEM mesh for complex 3D engineering structure.展开更多
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.展开更多
This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established with...This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established without any regularity assumption on the mesh.展开更多
Discretising a structure into elements is a key step in finite element(FE)analysis.The discretised geometry used to formulate an FE model can greatly affect accuracy and validity.This paper presents a unified dimensio...Discretising a structure into elements is a key step in finite element(FE)analysis.The discretised geometry used to formulate an FE model can greatly affect accuracy and validity.This paper presents a unified dimensionless parameter to generate a mesh of cubic FEs for the analysis of very long beams resting on an elastic foundation.A uniform beam resting on elastic foundation with various values of flexural stiffness and elastic supporting coefficients subject to static load and moving load is used to illustrate the application of the proposed parameter.The numerical results show that(a)Even if the values of the flexural stiffness of the beam and elastic supporting coefficient of the elastic foundation are different,the same proposed parameter“s”can ensure the same accuracy of the FE solution,but the accuracy may differ for use of the same element length;(b)The proposed dimensionless parameter“s”can indeed be used as a unified index to generate the mesh for a beam resting on elastic foundation,whereas the use of the same element length as a criterion may be misleading;(c)The errors between the FE and analytical solutions for the maximum vertical displacement,shear force and bending moment of the beam increase with the dimensionless parameter“s”;and(d)For the given allowable errors for the vertical displacement,shear force and bending moment of the beam under static load and moving load,the corresponding values of the proposed parameter are provided to guide the mesh generation.展开更多
Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation is proposed in this paper to improve the efficiency and accuracy of electromagnetic calculation. An image boundary extract...Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation is proposed in this paper to improve the efficiency and accuracy of electromagnetic calculation. An image boundary extraction algorithm is developed to map the image on the geometric domain. Identification algorithm for the location of nodes in polygon area is proposed to determine the state of the node. To promote the average quality of the mesh and the efficiency of mesh generation, a novel force-based mesh smoothing algorithm is proposed. One test case and a typical electromagnetic calculation are used to testify the effectiveness and efficiency of the proposed algorithm. The results demonstrate that the proposed algorithm can produce a high-quality mesh with less iteration.展开更多
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.展开更多
A modified weak Galerkin(MWG)finite element method is developed for solving the biharmonic equation.This method uses the same finite element space as that of the discontinuous Galerkin method,the space of discontinuou...A modified weak Galerkin(MWG)finite element method is developed for solving the biharmonic equation.This method uses the same finite element space as that of the discontinuous Galerkin method,the space of discontinuous polynomials on polytopal meshes.But its formulation is simple,symmetric,positive definite,and parameter independent,without any of six inter-element face-integral terms in the formulation of the discontinuous Galerkin method.Optimal order error estimates in a discrete H2 norm are established for the corresponding finite element solutions.Error estimates in the L^(2)norm are also derived with a sub-optimal order of convergence for the lowest-order element and an optimal order of convergence for all high-order of elements.The numerical results are presented to confirm the theory of convergence.展开更多
A mixed finite element solution of contact stresses in meshing gears is investigated with the consideration of coupled thermo-elastic deformation and impact behavior. A simulation procedure of finite element solution ...A mixed finite element solution of contact stresses in meshing gears is investigated with the consideration of coupled thermo-elastic deformation and impact behavior. A simulation procedure of finite element solution of meshing gears is developed. The versatility of the procedure for both numerical accuracy and computational efficiency is verified by numerical analysis of meshing gear teeth.展开更多
In this paper, firstly, a mathematical model for a specific kind of welted bifurcation is established, the parametric equation for the intersecting curve is resulted in. Secondly, a method for partitioning f...In this paper, firstly, a mathematical model for a specific kind of welted bifurcation is established, the parametric equation for the intersecting curve is resulted in. Secondly, a method for partitioning finite element meshes of the welted bifurcation is put forward, its main idea is that developing the main pipe surface and the branch pipe surface respectively, dividing meshes on each developing plane and obtaining meshes points, then transforming their plane coordinates into space coordinates. Finally, an applied program for finite element meshes auto-generation is simply introduced, which adopt ObjectARX technique and its running result can be shown in AutoCAD. The meshes generated in AutoCAD can be exported conveniently to most of finite element analysis soft wares, and the finite element computing result can satisfy the engineering precision requirement.展开更多
The cavitating flow around a Delft Twist-11 hydrofoil is simulated using the large eddy simulation approach.The volume-of-fluid method incorporated with the Schnerr-Sauer cavitation model is utilized to track the wate...The cavitating flow around a Delft Twist-11 hydrofoil is simulated using the large eddy simulation approach.The volume-of-fluid method incorporated with the Schnerr-Sauer cavitation model is utilized to track the water-vapor interface.Adaptive mesh refinement(AMR)is also applied to improve the simulation accuracy automatically.Two refinement levels are conducted to verify the dominance of AMR in predicting cavitating flows.Results show that cavitation features,including the U-type structure of shedding clouds,are consistent with experimental observations.Even a coarse mesh can precisely capture the phase field without increasing the total cell number significantly using mesh adaption.The predicted shedding frequency agrees fairly well with the experimental data under refinement level 2.This study illustrates that AMR is a promising approach to achieve accurate simulations for multiscale cavitating flows within limited computational costs.Finally,the force element method is currently adopted to investigate the lift and drag fluctuations during the evolution of cavitation structure.The mechanisms of lift and drag fluctuations due to cavitation and the interaction between vorticity forces and cavitation are explicitly revealed.展开更多
基金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.
文摘This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been published,but very few on those generated from a mesh using triangular elements.The use of triangular elements is aimed at extending the application of the approach to any shape of modeled devices.Basic concepts of the approach are presented in the case of electromagnetic devices.The procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB environment are also included.
基金supported by the Research Management Centre(RMC)of Multimedia University,Malaysia(Grant No.MMUI/220016).
文摘Polygonal finite elements remain an attractive option in finite element analysis due to their flexibility in modelingarbitrary shapes compared to triangles.In this study,a pentagonal membrane element was developed with thestrain approach for the first time.The element possesses invariance,and the equilibrium constraint was appliedto the assumed strain field using corrective coefficients.Inspired by the advancing front technique,a pentagonalmesh was generated,and the mesh quality was enhanced with Laplacian smoothing.The performance of thedeveloped pentagonal element was assessed in a few numerical tests,and the results revealed its suitability inmodeling the bending of beams.Besides,the numerical results are enhanced when pentagonal elements are usedin mesh transitions along boundaries to smoothen curved edges and capture distributed loads.
基金Supported by the National Natural Science Foundation of China (10671184)
文摘A lumped mass approximation scheme of a low order Crouzeix-Raviart type noncon- forming triangular finite element is proposed to a kind of nonlinear parabolic integro-differential equations. The L2 error estimate is derived on anisotropic meshes without referring to the traditional nonclassical elliptic projection.
基金Project supported by the National Natural Science Foundation of China (No. 10371113)
文摘The Wilson finite element method is considered to solve a class of two- dimensional second order elliptic boundary value problems. By using of the particular structure of the element and some new techniques, we obtain the superclose and global superconvergence on anisotropic meshes. Numerical example is also given to confirm our theoretical analysis.
文摘Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method, The solution accuracy is further improved by implementing an adaptive meshing technique to generaie finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems tire: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to moving heat source.
基金This research is supported by the NSF of China (10371113 10471133),SF of Henan ProvinceSF of Education Committee of Henan Province (2006110011)
文摘The main aim of this paper is to study the approximation to viscoelasticity type equations with a Crouzeix-Raviart type nonconforming finite element on the anisotropic meshes. The superclose property of the exact solution and the optimal error estimate of its derivative with respect to time are derived by using some novel techniques. Moreover, employing a postprocessing technique, the global superconvergence property for the discretization error of the postprocessed discrete solution to the solution itself is studied.
基金Foundation item: Supported by the NSF of China(10371113)Supported by the Foundation of Overseas Scholar of China(2001(119))Supported by the project of Creative Engineering of Province of China(2002(219))
文摘The combined hybrid finite element method is of an intrinsic mechanism of enhancing coarse-mesh-accuracy of lower order displacement schemes. It was confirmed that the combined hybrid scheme without energy error leads to enhancement of accuracy at coarse meshes, and that the combination parameter plays an important role in the enhancement. As an improvement of conforming bilinear Q(4)-plane element, the combined hybrid method adopted the most convenient quadrilateral displacements-stress mode, i.e.,the mode of compatible isoparametric bilinear displacements and pure constant stresses. By adjusting the combined parameter, the optimized version of the combined hybrid element was obtained and numerical tests indicated that this parameter-adjusted version behaves much better than Q(4)-element and is of high accuracy at coarse meshes. Due to elimination of stress parameters at the elemental level, this combined hybrid version is of the same computational cost as that of Q(4)-element.
基金Supported by the National Natural Science Foundation of China (A10102006)
文摘Mapping mesh generation is widely applied in pre-processes of Finite Element Method (FEM). In this study, the basic 3D mapping equations by Lagrange interpolating function are founded. Based these equations, a mapping pattern library, which maps essential configurations e.g. line, circle, rotary body, sphere etc. to hexahedral FEM mesh, has been built. Then available FEM mesh will be generated by clipping and assembling the mapped essential objects. Study case illustrates that the proposed method is simple and efficient to generate valid FEM mesh for complex 3D engineering structure.
基金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.
文摘This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established without any regularity assumption on the mesh.
基金the National Key Research and Development Program of China(Grant 2017YFB1201204)National Natural Science Foundation of China(Grants 51578552,U1334203).
文摘Discretising a structure into elements is a key step in finite element(FE)analysis.The discretised geometry used to formulate an FE model can greatly affect accuracy and validity.This paper presents a unified dimensionless parameter to generate a mesh of cubic FEs for the analysis of very long beams resting on an elastic foundation.A uniform beam resting on elastic foundation with various values of flexural stiffness and elastic supporting coefficients subject to static load and moving load is used to illustrate the application of the proposed parameter.The numerical results show that(a)Even if the values of the flexural stiffness of the beam and elastic supporting coefficient of the elastic foundation are different,the same proposed parameter“s”can ensure the same accuracy of the FE solution,but the accuracy may differ for use of the same element length;(b)The proposed dimensionless parameter“s”can indeed be used as a unified index to generate the mesh for a beam resting on elastic foundation,whereas the use of the same element length as a criterion may be misleading;(c)The errors between the FE and analytical solutions for the maximum vertical displacement,shear force and bending moment of the beam increase with the dimensionless parameter“s”;and(d)For the given allowable errors for the vertical displacement,shear force and bending moment of the beam under static load and moving load,the corresponding values of the proposed parameter are provided to guide the mesh generation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.52077203 and 61701467)the Natural Science Foundation of Zhejiang Province,China(Grant No.LY19E070003)。
文摘Two-dimensional finite element mesh generation algorithm for electromagnetic field calculation is proposed in this paper to improve the efficiency and accuracy of electromagnetic calculation. An image boundary extraction algorithm is developed to map the image on the geometric domain. Identification algorithm for the location of nodes in polygon area is proposed to determine the state of the node. To promote the average quality of the mesh and the efficiency of mesh generation, a novel force-based mesh smoothing algorithm is proposed. One test case and a typical electromagnetic calculation are used to testify the effectiveness and efficiency of the proposed algorithm. The results demonstrate that the proposed algorithm can produce a high-quality mesh with less iteration.
基金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.
基金M.Cui was supported in part by the National Natural Science Foundation of China(Grant No.11571026)the Beijing Municipal Natural Science Foundation of China(Grant No.1192003)Xiu Ye was supported in part by the National Science Foundation Grant DMS-1620016.
文摘A modified weak Galerkin(MWG)finite element method is developed for solving the biharmonic equation.This method uses the same finite element space as that of the discontinuous Galerkin method,the space of discontinuous polynomials on polytopal meshes.But its formulation is simple,symmetric,positive definite,and parameter independent,without any of six inter-element face-integral terms in the formulation of the discontinuous Galerkin method.Optimal order error estimates in a discrete H2 norm are established for the corresponding finite element solutions.Error estimates in the L^(2)norm are also derived with a sub-optimal order of convergence for the lowest-order element and an optimal order of convergence for all high-order of elements.The numerical results are presented to confirm the theory of convergence.
文摘A mixed finite element solution of contact stresses in meshing gears is investigated with the consideration of coupled thermo-elastic deformation and impact behavior. A simulation procedure of finite element solution of meshing gears is developed. The versatility of the procedure for both numerical accuracy and computational efficiency is verified by numerical analysis of meshing gear teeth.
文摘In this paper, firstly, a mathematical model for a specific kind of welted bifurcation is established, the parametric equation for the intersecting curve is resulted in. Secondly, a method for partitioning finite element meshes of the welted bifurcation is put forward, its main idea is that developing the main pipe surface and the branch pipe surface respectively, dividing meshes on each developing plane and obtaining meshes points, then transforming their plane coordinates into space coordinates. Finally, an applied program for finite element meshes auto-generation is simply introduced, which adopt ObjectARX technique and its running result can be shown in AutoCAD. The meshes generated in AutoCAD can be exported conveniently to most of finite element analysis soft wares, and the finite element computing result can satisfy the engineering precision requirement.
基金financially supported by the National Natural Science Foundation of China(Nos.U21A20126 and 52006197)the National Science Foundation of Zhejiang Province(Nos.LQ21E060012 and LR20E090001)the Key Research and Development Program of Zhejiang Province(No.2021C05006)。
文摘The cavitating flow around a Delft Twist-11 hydrofoil is simulated using the large eddy simulation approach.The volume-of-fluid method incorporated with the Schnerr-Sauer cavitation model is utilized to track the water-vapor interface.Adaptive mesh refinement(AMR)is also applied to improve the simulation accuracy automatically.Two refinement levels are conducted to verify the dominance of AMR in predicting cavitating flows.Results show that cavitation features,including the U-type structure of shedding clouds,are consistent with experimental observations.Even a coarse mesh can precisely capture the phase field without increasing the total cell number significantly using mesh adaption.The predicted shedding frequency agrees fairly well with the experimental data under refinement level 2.This study illustrates that AMR is a promising approach to achieve accurate simulations for multiscale cavitating flows within limited computational costs.Finally,the force element method is currently adopted to investigate the lift and drag fluctuations during the evolution of cavitation structure.The mechanisms of lift and drag fluctuations due to cavitation and the interaction between vorticity forces and cavitation are explicitly revealed.