In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and ...In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and co-state variables, and piecewise constant function is used to approximate control variables. Generally, the optimal conditions for the problem are solved iteratively until the control variable reaches error tolerance. In order to calculate all the variables individually and parallelly, we introduce a gradient recovery based two-grid method. First, we solve the small scaled optimal control problem on coarse grids. Next, we use the gradient recovery technique to recover the gradients of state and co-state variables. Finally, using the recovered variables, we solve the large scaled optimal control problem for all variables independently. Moreover, we estimate priori error for the proposed scheme, and use an example to validate the theoretical results.展开更多
In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element...In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.展开更多
The present work aims to evaluate the increase in the number of spot welds in the 16 × 16 type fuel assembly structure that connects guide thimbles and spacer grids, in order to provide a proper joint for this co...The present work aims to evaluate the increase in the number of spot welds in the 16 × 16 type fuel assembly structure that connects guide thimbles and spacer grids, in order to provide a proper joint for this connection. This new and improved process can provide more stiffness to the whole structure, since the number of spots raised from four to eight. A 3-D geometric model of a guide thimble section was generated in a CAD （computer aided design） program （SolidWorks）. After that, the geometric model was imported to a CAE （computer aided engineering） program （ANSYS Mechanical APDL, Release 14.0）, where the finite element model was built, considering the guide thimble geometry assembled with the spacer grid through the welded connections. Boundaries conditions were implemented in the model in order to simulate the correct physical behavior due to the operation of the fuel assembly inside the reactor. The analysis covered specific loads and displacements acting on the entire structure. The method used to solve this finite element analysis was a linear static simulation in order to perform the connection between a spacer grid cell and a guide thimble section. Hence, four models was evaluated, differing on the spot weld number in the spacer grid and guide thimble connection. The rotational stiffness results of each model were compared. The results acquired from four and eight spot weld were validated with physical test results. The behavior of the structure under the acting force/displacement and the related results of the analysis, mainly the stiffness, were satisfied. The results of this analysis were used to prove that the increasing spot welds number is an improvement in the dimensional stability when submitted to loads and displacements required on the fuel assembly design. This analysis aid to get more information of extreme importance such as, the pursuance to develop better manufacturing process and to improve the fuel assembly performance due to the increasing of the bum-up.展开更多
In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the m...In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.展开更多
In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of fini...In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.展开更多
Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow pheno...Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.展开更多
In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG) stabilized method and three-step finite element method(FEM),for large eddy simulation(LES) is developed to pred...In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG) stabilized method and three-step finite element method(FEM),for large eddy simulation(LES) is developed to predict the wind flow with high Reynolds numbers.Weak form of LES motion equation is combined with the SUPG stabilized term for the spatial finite element discretization.An explicit three-step scheme is implemented for the temporal discretization.For the numerical example of 2D wind flow over a square rib at Re=4.2×105,the Smagorinsky's subgrid-scale(SSGS) model,the DSGS model,and the DSGS model with Cabot near-wall model are applied,and their results are analyzed and compared with experimental results.Furthermore,numerical examples of 3D wind flow around a surface-mounted cube with different Reynolds numbers are performed using DSGS model with Cabot near-wall model based on the present stabilized method to study the wind field and compared with experimental and numerical results.Finally,vortex structures for wind flow around a surface-mounted cube are studied by present numerical method.Stable and satisfactory results are obtained,which are consistent with most of the measurements even under coarse mesh.展开更多
How to automatically generate three-dimensional finite element Delaunay mesh by a peifected node connection method is introduced, where nodes are generated based on existing elements, instead of independence of node c...How to automatically generate three-dimensional finite element Delaunay mesh by a peifected node connection method is introduced, where nodes are generated based on existing elements, instead of independence of node creation and elements generation in traditional node connection method. Therefore, Ihe the difficulty about how to automatically create nodes in the traditional method is overcome.展开更多
For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid an...For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.展开更多
Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine fini...Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.展开更多
How the outer substance could penetrate through the skin lies in the stratum corneum, because it is the main barrier in the multi-layers of the skin. Supposing the keratin cell with a special geometry as tetrakaidecah...How the outer substance could penetrate through the skin lies in the stratum corneum, because it is the main barrier in the multi-layers of the skin. Supposing the keratin cell with a special geometry as tetrakaidecahedron, the penetration property of stratum corneum was the key problem which was numerically simulated with finite element method. At first the discretization of the stratum corneum region was given in two steps： first, the discretization of the keratin cell; second, the discretization of fattiness that surrounds the keratin. Then there was the work of numerical simulation. In this procedure, the finite element method and the multi-grid method were used. The former was to obtain the discretization of basic elements; the latter was to decrease the high frequency error. At last the visualization of the numerical simulation was shown.展开更多
In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed...In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.展开更多
This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and...This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and temperature loads. In the calculation mesh, the contact surface of pair nodes is located at places on the arch dam where cracking is possible. A new effective iterative method, the mixed finite element method for friction-contact problems, is improved and used for nonlinear simulation of the cracking process. The forces acting on the structure are divided into two parts： external forces and contact forces. The displacement of the structure is chosen as the basic variable and the nodal contact force in the possible contact region of the local coordinate system is chosen as the iterative variable, so that the nonlinear iterative process is only limited within the possible contact surface and is much more economical. This method was used to simulate the cracking process of the Shuanghe Arch Dam in Southwest China. In order to prove the validity and accuracy of this method and to study the effect of thermal stress on arch dam cracking, three schemes were designed for calculation. Numerical results agree with actual measured data, proving that it is feasible to use this method to simulate the entire process of nonlinear arch dam cracking.展开更多
A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiat...A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases.展开更多
This paper presents a hybrid finite volume/finite element method for the incompressible generalized Newtonian fluid flow （Power-Law model）. The collocated （i.e. non-staggered） arrangement of variables is used on t...This paper presents a hybrid finite volume/finite element method for the incompressible generalized Newtonian fluid flow （Power-Law model）. The collocated （i.e. non-staggered） arrangement of variables is used on the unstructured triangular grids, and a fractional step projection method is applied for the velocity-pressure coupling. The cell-centered finite volume method is employed to discretize the momentum equation and the vertex-based finite element for the pressure Poisson equation. The momentum interpolation method is used to suppress unphysical pressure wiggles. Numerical experiments demonstrate that the current hybrid scheme has second order accuracy in both space and time. Results on flows in the lid-driven cavity and between parallel walls for Newtonian and Power-Law models are also in good agreement with the published solutions.展开更多
In this paper we study the convergence of adaptive finite element methods for the gen- eral non-attine equivalent quadrilateral and hexahedral elements on 1-irregular meshes with hanging nodes. Based on several basic ...In this paper we study the convergence of adaptive finite element methods for the gen- eral non-attine equivalent quadrilateral and hexahedral elements on 1-irregular meshes with hanging nodes. Based on several basic ingredients, such as quasi-orthogonality, estimator reduction and D6fler marking strategy, convergence of the adaptive finite element methods for the general second-order elliptic partial equations is proved. Our analysis is effective for all conforming Qm elements which covers both the two- and three-dimensional cases in a unified fashion.展开更多
A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid m...A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.展开更多
The method of determining grid sizes of finite element method of water quality models is presented on the basis of characteristics of pollutant distributed in the just downstream grid of discharge. A given length of t...The method of determining grid sizes of finite element method of water quality models is presented on the basis of characteristics of pollutant distributed in the just downstream grid of discharge. A given length of the grid is applied to calculate the minimum distance for the pollutant to complete mixing in the lateral direction in the just downstream grid of the discharge, and the minimum distance is regarded as the maximum width of the grid. Application of the maximum width of the grid and flow rate per unit width of the river results in the number of flow zones in the lateral direction of the river. Consequently division of grids and numerical calculation of water quality with finite element method can be carried out on the basis of available method and topographic map of the river course. Analysis of sensitivity shows that the effect of the width of the river and the effect of the length of the grids on the minimum number of flow zones in the lateral direction is relatively small, and that the effect of the given lateral mixing degree of the pollutants on the minimum number of flow zones is acceptable, which indicates that the method presented is relatively stable. The method can avoid the blindness and randomness in determining the number of flow zones, and can decrease the amount of computation under the condition that the accuracy requirement is met. The problem that determining the sizes of the grids must rely on empirical knowledge has been solved.展开更多
For the path dependency and nonlinearity introduced by incremental construction, numerical method has been widely used in deformation analysis of geo engineering.In the numerical simulation scheme commonly used in the...For the path dependency and nonlinearity introduced by incremental construction, numerical method has been widely used in deformation analysis of geo engineering.In the numerical simulation scheme commonly used in the past, the excavating loads are extracted from nodal stresses, which are deduced linearly from the stresses at Gauss point in finite element method.The unneglectable calculation error is contained in this process when elastic plastic constitutive model is employed.The error mentioned above is analyzed in detail.Based on the analysis of excavation process and the principle of finite element theory, a new simulation scheme for excavation is proposed.At the end of this paper, an application in rock engineering is given out.展开更多
INTERNODES is a general purpose method to deal with non-conforming discretizations of partial differential equations on 2D and 3D regions partitioned into two or several disjoint subdomains. It exploits two intergrid ...INTERNODES is a general purpose method to deal with non-conforming discretizations of partial differential equations on 2D and 3D regions partitioned into two or several disjoint subdomains. It exploits two intergrid interpolation operators, one for transfering the Dirichlet trace across the interfaces, and the other for the Neumann trace. In this paper, in every subdomain the original problem is discretized by either the finite element method (FEM) or the spectral element method (SEM or hp-FEM), using a priori non-matching grids and piecewise polynomials of different degrees. Other discretization methods, however, can be used. INTERNODES can also be applied to heterogeneous or multiphysics problems, that is, problems that feature different differential operators inside adjacent subdomains. For instance, in this paper we apply the INTERNODES method to a Stokes- Darcy coupled problem that models the filtration of fluids in porous media. Our results highlight the flexibility of the method as well as its optimal rate of convergence with respect to the grid size and the polynomial degree.展开更多
文摘In this paper, the optimal control problem of parabolic integro-differential equations is solved by gradient recovery based two-grid finite element method. Piecewise linear functions are used to approximate state and co-state variables, and piecewise constant function is used to approximate control variables. Generally, the optimal conditions for the problem are solved iteratively until the control variable reaches error tolerance. In order to calculate all the variables individually and parallelly, we introduce a gradient recovery based two-grid method. First, we solve the small scaled optimal control problem on coarse grids. Next, we use the gradient recovery technique to recover the gradients of state and co-state variables. Finally, using the recovered variables, we solve the large scaled optimal control problem for all variables independently. Moreover, we estimate priori error for the proposed scheme, and use an example to validate the theoretical results.
文摘In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.
文摘The present work aims to evaluate the increase in the number of spot welds in the 16 × 16 type fuel assembly structure that connects guide thimbles and spacer grids, in order to provide a proper joint for this connection. This new and improved process can provide more stiffness to the whole structure, since the number of spots raised from four to eight. A 3-D geometric model of a guide thimble section was generated in a CAD （computer aided design） program （SolidWorks）. After that, the geometric model was imported to a CAE （computer aided engineering） program （ANSYS Mechanical APDL, Release 14.0）, where the finite element model was built, considering the guide thimble geometry assembled with the spacer grid through the welded connections. Boundaries conditions were implemented in the model in order to simulate the correct physical behavior due to the operation of the fuel assembly inside the reactor. The analysis covered specific loads and displacements acting on the entire structure. The method used to solve this finite element analysis was a linear static simulation in order to perform the connection between a spacer grid cell and a guide thimble section. Hence, four models was evaluated, differing on the spot weld number in the spacer grid and guide thimble connection. The rotational stiffness results of each model were compared. The results acquired from four and eight spot weld were validated with physical test results. The behavior of the structure under the acting force/displacement and the related results of the analysis, mainly the stiffness, were satisfied. The results of this analysis were used to prove that the increasing spot welds number is an improvement in the dimensional stability when submitted to loads and displacements required on the fuel assembly design. This analysis aid to get more information of extreme importance such as, the pursuance to develop better manufacturing process and to improve the fuel assembly performance due to the increasing of the bum-up.
文摘In this paper, we present the a posteriori error estimate of two-grid mixed finite element methods by averaging techniques for semilinear elliptic equations. We first propose the two-grid algorithms to linearize the mixed method equations. Then, the averaging technique is used to construct the a posteriori error estimates of the two-grid mixed finite element method and theoretical analysis are given for the error estimators. Finally, we give some numerical examples to verify the reliability and efficiency of the a posteriori error estimator.
基金Project supported by the National Natural Science Foundation of China(Nos.11671157 and11826212)
文摘In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.
基金King Mongkut’s University of Technology North Bangkok (KMUTNB)the Office of the Higher Education Commission (OHEC)the National Metal and Materials Technology Center (MTEC) for supporting this research work
文摘Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.
基金Project supported by the National Natural Science Foundation of China(No.51078230)the Research Fund for the Doctoral Program of Higher Education of China(No.200802480056)the Key Project of Fund of Science and Technology Development of Shanghai(No.10JC1407900),China
文摘In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG) stabilized method and three-step finite element method(FEM),for large eddy simulation(LES) is developed to predict the wind flow with high Reynolds numbers.Weak form of LES motion equation is combined with the SUPG stabilized term for the spatial finite element discretization.An explicit three-step scheme is implemented for the temporal discretization.For the numerical example of 2D wind flow over a square rib at Re=4.2×105,the Smagorinsky's subgrid-scale(SSGS) model,the DSGS model,and the DSGS model with Cabot near-wall model are applied,and their results are analyzed and compared with experimental results.Furthermore,numerical examples of 3D wind flow around a surface-mounted cube with different Reynolds numbers are performed using DSGS model with Cabot near-wall model based on the present stabilized method to study the wind field and compared with experimental and numerical results.Finally,vortex structures for wind flow around a surface-mounted cube are studied by present numerical method.Stable and satisfactory results are obtained,which are consistent with most of the measurements even under coarse mesh.
基金This project is supported by Provincial Natural Science foundation of Guangdong!(970516)
文摘How to automatically generate three-dimensional finite element Delaunay mesh by a peifected node connection method is introduced, where nodes are generated based on existing elements, instead of independence of node creation and elements generation in traditional node connection method. Therefore, Ihe the difficulty about how to automatically create nodes in the traditional method is overcome.
文摘For two-dimension nonlinear convection diffusion equation, a two-grid method of characteristics finite-element solution was constructed. In this method the nonlinear iterations is only to execute on the coarse grid and the fine-grid solution can be obtained in a single linear step. For the nonlinear convection-dominated diffusion equation, this method can not only stabilize the numerical oscillation but also accelerate the convergence and improve the computational efficiency. The error analysis demonstrates if the mesh sizes between coarse-grid and fine-grid satisfy the certain relationship, the two-grid solution and the characteristics finite-element solution have the same order of accuracy. The numerical is more efficient than that of characteristics example confirms that the two-grid method finite-element method.
基金Project supported by the National Natural Science Foundation of China (No. 10371096)
文摘Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.
文摘How the outer substance could penetrate through the skin lies in the stratum corneum, because it is the main barrier in the multi-layers of the skin. Supposing the keratin cell with a special geometry as tetrakaidecahedron, the penetration property of stratum corneum was the key problem which was numerically simulated with finite element method. At first the discretization of the stratum corneum region was given in two steps： first, the discretization of the keratin cell; second, the discretization of fattiness that surrounds the keratin. Then there was the work of numerical simulation. In this procedure, the finite element method and the multi-grid method were used. The former was to obtain the discretization of basic elements; the latter was to decrease the high frequency error. At last the visualization of the numerical simulation was shown.
文摘In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.
基金supported by the National Nature Science Foundation of China (Grant No 90510017)
文摘This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and temperature loads. In the calculation mesh, the contact surface of pair nodes is located at places on the arch dam where cracking is possible. A new effective iterative method, the mixed finite element method for friction-contact problems, is improved and used for nonlinear simulation of the cracking process. The forces acting on the structure are divided into two parts： external forces and contact forces. The displacement of the structure is chosen as the basic variable and the nodal contact force in the possible contact region of the local coordinate system is chosen as the iterative variable, so that the nonlinear iterative process is only limited within the possible contact surface and is much more economical. This method was used to simulate the cracking process of the Shuanghe Arch Dam in Southwest China. In order to prove the validity and accuracy of this method and to study the effect of thermal stress on arch dam cracking, three schemes were designed for calculation. Numerical results agree with actual measured data, proving that it is feasible to use this method to simulate the entire process of nonlinear arch dam cracking.
基金The project supported by the Research Grant Council of Hong Kong (HKUST 722196E, 6039197E)the National Natural Science Foundation of China(59809003)the Foundation of University Key Teacher by the Chinese Ministry of Education
文摘A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases.
基金supported by the National Natural Science Foundation of China (10771134).
文摘This paper presents a hybrid finite volume/finite element method for the incompressible generalized Newtonian fluid flow （Power-Law model）. The collocated （i.e. non-staggered） arrangement of variables is used on the unstructured triangular grids, and a fractional step projection method is applied for the velocity-pressure coupling. The cell-centered finite volume method is employed to discretize the momentum equation and the vertex-based finite element for the pressure Poisson equation. The momentum interpolation method is used to suppress unphysical pressure wiggles. Numerical experiments demonstrate that the current hybrid scheme has second order accuracy in both space and time. Results on flows in the lid-driven cavity and between parallel walls for Newtonian and Power-Law models are also in good agreement with the published solutions.
基金supported by the Special Funds for Major State Basic Research Project (No. 2005CB321701)
文摘In this paper we study the convergence of adaptive finite element methods for the gen- eral non-attine equivalent quadrilateral and hexahedral elements on 1-irregular meshes with hanging nodes. Based on several basic ingredients, such as quasi-orthogonality, estimator reduction and D6fler marking strategy, convergence of the adaptive finite element methods for the general second-order elliptic partial equations is proved. Our analysis is effective for all conforming Qm elements which covers both the two- and three-dimensional cases in a unified fashion.
文摘A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.
文摘The method of determining grid sizes of finite element method of water quality models is presented on the basis of characteristics of pollutant distributed in the just downstream grid of discharge. A given length of the grid is applied to calculate the minimum distance for the pollutant to complete mixing in the lateral direction in the just downstream grid of the discharge, and the minimum distance is regarded as the maximum width of the grid. Application of the maximum width of the grid and flow rate per unit width of the river results in the number of flow zones in the lateral direction of the river. Consequently division of grids and numerical calculation of water quality with finite element method can be carried out on the basis of available method and topographic map of the river course. Analysis of sensitivity shows that the effect of the width of the river and the effect of the length of the grids on the minimum number of flow zones in the lateral direction is relatively small, and that the effect of the given lateral mixing degree of the pollutants on the minimum number of flow zones is acceptable, which indicates that the method presented is relatively stable. The method can avoid the blindness and randomness in determining the number of flow zones, and can decrease the amount of computation under the condition that the accuracy requirement is met. The problem that determining the sizes of the grids must rely on empirical knowledge has been solved.
文摘For the path dependency and nonlinearity introduced by incremental construction, numerical method has been widely used in deformation analysis of geo engineering.In the numerical simulation scheme commonly used in the past, the excavating loads are extracted from nodal stresses, which are deduced linearly from the stresses at Gauss point in finite element method.The unneglectable calculation error is contained in this process when elastic plastic constitutive model is employed.The error mentioned above is analyzed in detail.Based on the analysis of excavation process and the principle of finite element theory, a new simulation scheme for excavation is proposed.At the end of this paper, an application in rock engineering is given out.
文摘INTERNODES is a general purpose method to deal with non-conforming discretizations of partial differential equations on 2D and 3D regions partitioned into two or several disjoint subdomains. It exploits two intergrid interpolation operators, one for transfering the Dirichlet trace across the interfaces, and the other for the Neumann trace. In this paper, in every subdomain the original problem is discretized by either the finite element method (FEM) or the spectral element method (SEM or hp-FEM), using a priori non-matching grids and piecewise polynomials of different degrees. Other discretization methods, however, can be used. INTERNODES can also be applied to heterogeneous or multiphysics problems, that is, problems that feature different differential operators inside adjacent subdomains. For instance, in this paper we apply the INTERNODES method to a Stokes- Darcy coupled problem that models the filtration of fluids in porous media. Our results highlight the flexibility of the method as well as its optimal rate of convergence with respect to the grid size and the polynomial degree.
