We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is ...We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations....A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.展开更多
An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same t...An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.展开更多
The Weak Galerkin (WG) finite element method for the unsteady Stokes equations in the primary velocity-pressure formulation is introduced in this paper. Optimal-order error estimates are established for the correspond...The Weak Galerkin (WG) finite element method for the unsteady Stokes equations in the primary velocity-pressure formulation is introduced in this paper. Optimal-order error estimates are established for the corresponding numerical approximation in an H1 norm for the velocity, and L2 norm for both the velocity and the pressure by use of the Stokes projection.展开更多
Through the construction of a new ramp function, the element-flee Galerkin method and finite element coupling method were applied to the whole field, and was made fit for the structure of element nodes within the inte...Through the construction of a new ramp function, the element-flee Galerkin method and finite element coupling method were applied to the whole field, and was made fit for the structure of element nodes within the interface regions, both satisfying the essential boundary conditions and deploying meshless nodes and finite elements in a convenient and flexible way, which can meet the requirements of computation for complicated field. The comparison between the results of the present study and the corresponding analytical solutions shows this method is feasible and effective.展开更多
An adaptive mixed least squares Galerkin/Petrov finite element method （FEM） is developed for stationary conduction convection problems. The mixed least squares Galerkin/Petrov FEM is consistent and stable for any co...An adaptive mixed least squares Galerkin/Petrov finite element method （FEM） is developed for stationary conduction convection problems. The mixed least squares Galerkin/Petrov FEM is consistent and stable for any combination of discrete velocity and pressure spaces without requiring the Babuska-Brezzi stability condition. Using the general theory of Verfiirth, the posteriori error estimates of the residual type are derived. Finally, numerical tests are presented to illustrate the effectiveness of the method.展开更多
A numerical method for coupled deformation between sheet metal and flexible-die was proposed. Based on the updated Lagrangian （UL） formulation, the elastoplastic deformation of sheet metal was analyzed with finite e...A numerical method for coupled deformation between sheet metal and flexible-die was proposed. Based on the updated Lagrangian （UL） formulation, the elastoplastic deformation of sheet metal was analyzed with finite element method （FEM） and the bulk deformation of flexible-die was analyzed with element-free Galerkin method （EFGM）. The frictional contact between sheet metal and flexible-die was treated by the penalty function method. The sheet elastic flexible-die bulging process was analyzed with the FEM-EFGM program for coupled deformation between sheet metal and bulk flexible-die, called CDSB-FEM-EFGM for short. Compared with finite element code DEFORM-2D and experiment results, the CDSB-FEM-EFGM program is feasible. This method provides a suitable numerical method to analyze sheet flexible-die forming.展开更多
In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy o...In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.展开更多
The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can ov...The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.展开更多
A new membrane finite element method for modeling fluid flow in a porous medium is presented in order to quickly and accurately simulate the geo-membrane fabric used in civil engineering. It is based on discontinuous ...A new membrane finite element method for modeling fluid flow in a porous medium is presented in order to quickly and accurately simulate the geo-membrane fabric used in civil engineering. It is based on discontinuous finite element theory, and can be easily coupled with the normal Galerkin finite element method. Based on the saturated seepage equation, the element coefficient matrix of the membrane element method is derived, and a geometric transform relation for the membrane element between a global coordinate system and a local coordinate system is obtained. A method for the determination of the fluid flux conductivity of the membrane element is presented. This method provides a basis for determining discontinuous parameters in discontinuous finite element theory. An anti-seepage problem regarding the foundation of a building is analyzed by coupling the membrane finite element method with the normal Galerkin finite element method. The analysis results demonstrate the utility and superiority of the membrane finite element method in fluid flow analysis of a porous medium.展开更多
In this article, we consider a two-dimensional symmetric space-fractional diffusion equation in which the space fractional derivatives are defined in Riesz potential sense. The well-posed feature is guaranteed by ener...In this article, we consider a two-dimensional symmetric space-fractional diffusion equation in which the space fractional derivatives are defined in Riesz potential sense. The well-posed feature is guaranteed by energy inequality. To solve the diffusion equation, a fully discrete form is established by employing Crank-Nicolson technique in time and Galerkin finite element method in space. The stability and convergence are proved and the stiffness matrix is given analytically. Three numerical examples are given to confirm our theoretical analysis in which we find that even with the same initial condition, the classical and fractional diffusion equations perform differently but tend to be uniform diffusion at last.展开更多
A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for-...A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.展开更多
The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite ...The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.展开更多
In this paper, the approximation of stationary equations of the semiconductor devices with mixed boundary conditions is considered. Two schemes are proposed for the system. One is Glerkin discrete scheme, the other is...In this paper, the approximation of stationary equations of the semiconductor devices with mixed boundary conditions is considered. Two schemes are proposed for the system. One is Glerkin discrete scheme, the other is hybrid variable discrete scheme. A convergence analysis is also given.展开更多
This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-d...This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.展开更多
The improved Boussinesq equation is solved with classical finite element method using the most basic Lagrange element k = 1, which leads us to a second order nonlinear ordinary differential equations system in time;th...The improved Boussinesq equation is solved with classical finite element method using the most basic Lagrange element k = 1, which leads us to a second order nonlinear ordinary differential equations system in time;this can be solved by any standard accurate numerical method for example Runge-Kutta-Fehlberg. The technique is validated with a typical example and a fourth order convergence in space is confirmed;the 1- and 2-soliton solutions are used to simulate wave travel, wave splitting and interaction;solution blow up is described graphically. The computer symbolic system MathLab is quite used for numerical simulation in this paper;the known results in the bibliography are confirmed.展开更多
Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space(Pk(K),P_(k−1)(∂K),[P_(k−1)(K)]^(2)).Optimal order a priori error estimates for both spac...Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space(Pk(K),P_(k−1)(∂K),[P_(k−1)(K)]^(2)).Optimal order a priori error estimates for both space-discrete scheme and implicit fully discrete scheme are derived in L1(L2)norm.This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes.Finite element algorithm presented here can contribute to a variety of hyperbolic problems where physical domain consists of heterogeneous media.展开更多
In this work,a modified weak Galerkin finite element method is proposed for solving second order linear parabolic singularly perturbed convection-diffusion equations.The key feature of the proposed method is to replac...In this work,a modified weak Galerkin finite element method is proposed for solving second order linear parabolic singularly perturbed convection-diffusion equations.The key feature of the proposed method is to replace the classical gradient and divergence operators by the modified weak gradient and modified divergence operators,respectively.We apply the backward finite difference method in time and the modified weak Galerkin finite element method in space on uniform mesh.The stability analyses are presented for both semi-discrete and fully-discrete modified weak Galerkin finite element methods.Optimal order of convergences are obtained in suitable norms.We have achieved the same accuracy with the weak Galerkin method while the degrees of freedom are reduced in our method.Various numerical examples are presented to support the theoretical results.It is theoretically and numerically shown that the method is quite stable.展开更多
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
基金Supported by the National Natural Science Foundation of China(50976072,51106099,10902070)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)the Science Foundation for the Excellent Youth Scholar of Higher Education of Shanghai(slg09003)~~
文摘A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.
基金Supported by the National Natural Science Foundation of China (10601022)Natural Science Foundation of Inner Mongolia Autonomous Region (200607010106)Youth Science Foundation of Inner Mongolia University(ND0702)
文摘An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.
文摘The Weak Galerkin (WG) finite element method for the unsteady Stokes equations in the primary velocity-pressure formulation is introduced in this paper. Optimal-order error estimates are established for the corresponding numerical approximation in an H1 norm for the velocity, and L2 norm for both the velocity and the pressure by use of the Stokes projection.
文摘Through the construction of a new ramp function, the element-flee Galerkin method and finite element coupling method were applied to the whole field, and was made fit for the structure of element nodes within the interface regions, both satisfying the essential boundary conditions and deploying meshless nodes and finite elements in a convenient and flexible way, which can meet the requirements of computation for complicated field. The comparison between the results of the present study and the corresponding analytical solutions shows this method is feasible and effective.
基金supported by the National Natural Science Foundation of China(Nos.10871156 and 11171269)the Fund of Xi'an Jiaotong University(No.2009xjtujc30)
文摘An adaptive mixed least squares Galerkin/Petrov finite element method （FEM） is developed for stationary conduction convection problems. The mixed least squares Galerkin/Petrov FEM is consistent and stable for any combination of discrete velocity and pressure spaces without requiring the Babuska-Brezzi stability condition. Using the general theory of Verfiirth, the posteriori error estimates of the residual type are derived. Finally, numerical tests are presented to illustrate the effectiveness of the method.
基金Project(51275130)supported by the National Natural Science Foundation of China
文摘A numerical method for coupled deformation between sheet metal and flexible-die was proposed. Based on the updated Lagrangian （UL） formulation, the elastoplastic deformation of sheet metal was analyzed with finite element method （FEM） and the bulk deformation of flexible-die was analyzed with element-free Galerkin method （EFGM）. The frictional contact between sheet metal and flexible-die was treated by the penalty function method. The sheet elastic flexible-die bulging process was analyzed with the FEM-EFGM program for coupled deformation between sheet metal and bulk flexible-die, called CDSB-FEM-EFGM for short. Compared with finite element code DEFORM-2D and experiment results, the CDSB-FEM-EFGM program is feasible. This method provides a suitable numerical method to analyze sheet flexible-die forming.
文摘In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.
文摘The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.
基金supported by the National Natural Science Foundation of China (Grant No. 50779012)
文摘A new membrane finite element method for modeling fluid flow in a porous medium is presented in order to quickly and accurately simulate the geo-membrane fabric used in civil engineering. It is based on discontinuous finite element theory, and can be easily coupled with the normal Galerkin finite element method. Based on the saturated seepage equation, the element coefficient matrix of the membrane element method is derived, and a geometric transform relation for the membrane element between a global coordinate system and a local coordinate system is obtained. A method for the determination of the fluid flux conductivity of the membrane element is presented. This method provides a basis for determining discontinuous parameters in discontinuous finite element theory. An anti-seepage problem regarding the foundation of a building is analyzed by coupling the membrane finite element method with the normal Galerkin finite element method. The analysis results demonstrate the utility and superiority of the membrane finite element method in fluid flow analysis of a porous medium.
文摘In this article, we consider a two-dimensional symmetric space-fractional diffusion equation in which the space fractional derivatives are defined in Riesz potential sense. The well-posed feature is guaranteed by energy inequality. To solve the diffusion equation, a fully discrete form is established by employing Crank-Nicolson technique in time and Galerkin finite element method in space. The stability and convergence are proved and the stiffness matrix is given analytically. Three numerical examples are given to confirm our theoretical analysis in which we find that even with the same initial condition, the classical and fractional diffusion equations perform differently but tend to be uniform diffusion at last.
基金the National Natural Science Foundation of China (Grants 41372301 and 51349011)the Preeminent Youth Talent Project of Southwest University of Science and Technology (Grant 13zx9109)
文摘A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.
基金This work was supported in part by the National Science Foundation under grant DMS-1620288。
文摘The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.
文摘In this paper, the approximation of stationary equations of the semiconductor devices with mixed boundary conditions is considered. Two schemes are proposed for the system. One is Glerkin discrete scheme, the other is hybrid variable discrete scheme. A convergence analysis is also given.
基金supported by the National Natural Science Foundation of China (Grants 11571223, 51404160)Shanxi Province Science Foundation for Youths (Grant 2014021025-1)
文摘This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
文摘The improved Boussinesq equation is solved with classical finite element method using the most basic Lagrange element k = 1, which leads us to a second order nonlinear ordinary differential equations system in time;this can be solved by any standard accurate numerical method for example Runge-Kutta-Fehlberg. The technique is validated with a typical example and a fourth order convergence in space is confirmed;the 1- and 2-soliton solutions are used to simulate wave travel, wave splitting and interaction;solution blow up is described graphically. The computer symbolic system MathLab is quite used for numerical simulation in this paper;the known results in the bibliography are confirmed.
文摘Weak Galerkin finite element method is introduced for solving wave equation with interface on weak Galerkin finite element space(Pk(K),P_(k−1)(∂K),[P_(k−1)(K)]^(2)).Optimal order a priori error estimates for both space-discrete scheme and implicit fully discrete scheme are derived in L1(L2)norm.This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes.Finite element algorithm presented here can contribute to a variety of hyperbolic problems where physical domain consists of heterogeneous media.
基金supported in part by National Natural Science Foundation of China (No.11871038).
文摘In this work,a modified weak Galerkin finite element method is proposed for solving second order linear parabolic singularly perturbed convection-diffusion equations.The key feature of the proposed method is to replace the classical gradient and divergence operators by the modified weak gradient and modified divergence operators,respectively.We apply the backward finite difference method in time and the modified weak Galerkin finite element method in space on uniform mesh.The stability analyses are presented for both semi-discrete and fully-discrete modified weak Galerkin finite element methods.Optimal order of convergences are obtained in suitable norms.We have achieved the same accuracy with the weak Galerkin method while the degrees of freedom are reduced in our method.Various numerical examples are presented to support the theoretical results.It is theoretically and numerically shown that the method is quite stable.
