In this paper we study the higher accuracy methods - the extrapolation and defect correction for the semidiscrete Galerkin approximations to the solutions of Sobolev and viscoelasticity type equations. The global extr...In this paper we study the higher accuracy methods - the extrapolation and defect correction for the semidiscrete Galerkin approximations to the solutions of Sobolev and viscoelasticity type equations. The global extrapolation and the correction approximations of third order, rather than the pointwise extrapolation results are presented.展开更多
The shape and thickness qualities of strip are influenced by the vibration of rolling mill. At present, the researches on the vibration of rolling mill are mainly the vertical vibration and torsional vibration of sing...The shape and thickness qualities of strip are influenced by the vibration of rolling mill. At present, the researches on the vibration of rolling mill are mainly the vertical vibration and torsional vibration of single stand mill, the study on the vibration of tandem rolling mill is rare. For the vibration of tandem rolling mill, the key problem is the vibration of the moving strip between stands. In this paper, considering the dynamic of moving strip and rolling theory, the vertical vibration of moving strip in the rolling process was proposed. Take the moving strip between the two mills of tandem rolling mill in the rolling process as subject investigated, according to the theory of moving beam, the vertical vibration model of moving strip in the rolling process was established. The partial differential equation was discretized by Galerkin truncation method. The natural frequency and stability of the moving strip were investigated and the numerical simulation in time domain was made. Simulation results show that, the natural frequency was strongly influenced by the rolling velocity and tension. With increasing of the rolling velocity, the first three natural frequencies decrease, the fourth natural frequency increases; with increasing of the unit tension, when the rolling velocity is high and low, respectively, the low order dimensionless natural frequency gradually decreases and increases, respectively. According to the stability of moving strip, the critical speed was determined, and the matching relationship of the tension and rolling velocity was also determined. This model can be used to study the stability of moving strip, improve the quality of strip and develop new rolling technology from the aspect of dynamics.展开更多
The convergence and superconvergence properties of the discontinuous Galerkin （DG） method for a singularly perturbed model problem in one-dimensional setting are studied. By applying the DG method with appropriately...The convergence and superconvergence properties of the discontinuous Galerkin （DG） method for a singularly perturbed model problem in one-dimensional setting are studied. By applying the DG method with appropriately chosen numerical traces, the existence and uniqueness of the DG solution, the optimal order L2 error bounds, and 2p＋ 1-order superconvergence of the numerical traces are established. The numerical results indicate that the DG method does not produce any oscillation even under the uniform mesh. Numerical experiments demonstrate that, under the uniform mesh, it seems impossible to obtain the uniform superconvergence of the numerical traces. Nevertheless, thanks to the implementation of the so-called Shishkin-type mesh, the uniform 2p ＋ 1-order superconvergence is observed numerically.展开更多
With the improved moving least-squares (IMLS) approximation, an orthogonal function system with a weight function is used as the basis function. The combination of the element-free Galerkin (EFG) method and the IMLS a...With the improved moving least-squares (IMLS) approximation, an orthogonal function system with a weight function is used as the basis function. The combination of the element-free Galerkin (EFG) method and the IMLS approximation leads to the development of the improved element-free Galerkin (IEFG) method. In this paper, the IEFG method is applied to study the partial differential equations that control the heat flow in three-dimensional space. With the IEFG technique, the Galerkin weak form is employed to develop the discretized system equations, and the penalty method is applied to impose the essential boundary conditions. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the transient heat conduction equations and the boundary and initial conditions are time dependent, the scaling parameter, number of nodes and time step length are considered in a convergence study.展开更多
The paper presents the improved element-free Galerkin （IEFG） method for three-dimensional wave propa- gation. The improved moving least-squares （IMLS） approx- imation is employed to construct the shape function, w...The paper presents the improved element-free Galerkin （IEFG） method for three-dimensional wave propa- gation. The improved moving least-squares （IMLS） approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares （MLS） approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.展开更多
The aeroelastic analysis of high-altitude, long-endurance （HALE） aircraft that features high-aspect-ratio flexible wings needs take into account structural geometrical nonlinearities and dynamic stall. For a generic...The aeroelastic analysis of high-altitude, long-endurance （HALE） aircraft that features high-aspect-ratio flexible wings needs take into account structural geometrical nonlinearities and dynamic stall. For a generic nonlinear aeroelastic system, besides the stability boundary, the characteristics of the limit-cycle oscillation （LCO） should also be accurately predicted. In order to conduct nonlinear aeroelastic analysis of high-aspect-ratio flexible wings, a first-order, state-space model is developed by combining a geometrically exact, nonlinear anisotropic beam model with nonlinear ONERA （Edlin） dynamic stall model. The present investigations focus on the initiation and sustaining mechanism of the LCO and the effects of flight speed and drag on aeroelastic behaviors. Numerical results indicate that structural geometrical nonlinearities could lead to the LCO without stall occurring. As flight speed increases, dynamic stall becomes dominant and the LCO increasingly complicated. Drag could be negligible for LCO type, but should be considered to exactly predict the onset speed of flutter or LCO of high-aspect-ratio flexible wings.展开更多
In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the II...In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.展开更多
In this paper, we consider the local discontinuous Galerkin method （LDG） for solving singularly perturbed convection-diffusion problems in one- and two-dimensional settings. The existence and uniqueness of the LDG s...In this paper, we consider the local discontinuous Galerkin method （LDG） for solving singularly perturbed convection-diffusion problems in one- and two-dimensional settings. The existence and uniqueness of the LDG solutions are verified. Numerical experiments demonstrate that it seems impossible to obtain uniform superconvergence for numerical fluxes under uniform meshes. Thanks to the implementation of two-type different anisotropic meshes, i.e., the Shishkin and an improved grade meshes, the uniform 2p ＋ i-order superconvergence is observed numerically for both one-dimensional and twodimensional cases.展开更多
In this paper, a relaxed Hermitian and skew-Hermitian splitting （RHSS） preconditioner is proposed for saddle point problems from the element-free Galerkin （EFG） discretization method. The EFG method is one of the ...In this paper, a relaxed Hermitian and skew-Hermitian splitting （RHSS） preconditioner is proposed for saddle point problems from the element-free Galerkin （EFG） discretization method. The EFG method is one of the most widely used meshfree methods for solving partial differential equations. The RHSS preconditioner is constructed much closer to the coefficient matrix than the well-known HSS preconditioner, resulting in a RHSS fixed-point iteration. Convergence of the RHSS iteration is analyzed and an optimal parameter, which minimizes the spectral radius of the iteration matrix is described. Using the RHSS pre- conditioner to accelerate the convergence of some Krylov subspace methods （like GMRES） is also studied. Theoretical analyses show that the eigenvalues of the RHSS precondi- tioned matrix are real and located in a positive interval. Eigenvector distribution and an upper bound of the degree of the minimal polynomial of the preconditioned matrix are obtained. A practical parameter is suggested in implementing the RHSS preconditioner. Finally, some numerical experiments are illustrated to show the effectiveness of the new preconditioner.展开更多
Irregular plates are very common structures in engineering,such as ore structures in mining.In this work,the Galerkin solution to the problem of a Kirchhoff plate lying on the Winkler foundation with two edges simply ...Irregular plates are very common structures in engineering,such as ore structures in mining.In this work,the Galerkin solution to the problem of a Kirchhoff plate lying on the Winkler foundation with two edges simply supported and the other two clamped supported is derived.Coordinate transformation technique is used during the solving process so that the solution is suitable to irregular shaped plates.The mechanical model and the solution proposed are then used to model the crown pillars between two adjacent levels in Sanshandao gold mine,which uses backfill method for mining operation.After that,an objective function,which takes security,economic profits and filling effect into consideration,is built to evaluate design proposals.Thickness optimizations for crown pillars are finally conducted in both conditions that the vertical stiffness of the foundation is known and unknown.The procedure presented in the work provides the guidance in thickness designing of complex shaped crown pillars and the preparation of backfill materials,thus to achieve the best balance between security and profits.展开更多
This paper proposes and analyzes a class of robust globally divergence-free weak Galerkin （WG） finite element methods for Stokes equations. The new methods use the Pk/Pk-1 （k ≥ 1） discontinuous finite element com...This paper proposes and analyzes a class of robust globally divergence-free weak Galerkin （WG） finite element methods for Stokes equations. The new methods use the Pk/Pk-1 （k ≥ 1） discontinuous finite element combination for velocity and pressure in the interior of elements, and piecewise P1/Pk （l = k - 1, k） for the trace approximations of the ve- locity and pressure on the inter-element boundaries. Our methods not only yield globally divergence-free velocity solutions, but also have uniform error estimates with respect to the Reynolds number. Numerical experiments are provided to show the robustness of the proposed methods.展开更多
High-order Discontinuous Galerkin(DG) methods have been receiving more and more attentions in the area of Computational Fluid Dynamics(CFD) because of their high-accuracy property. However, it is still a challenge to ...High-order Discontinuous Galerkin(DG) methods have been receiving more and more attentions in the area of Computational Fluid Dynamics(CFD) because of their high-accuracy property. However, it is still a challenge to obtain converged solution rapidly when solving the Reynolds Averaged Navier–Stokes(RANS) equations since the turbulence models significantly increase the nonlinearity of discretization system. The overall goal of this research is to develop an Artificial Neural Networks(ANNs) model with low complexity acting as an algebraic turbulence model to estimate the turbulence eddy viscosity for RANS. The ANN turbulence model is off-line trained using the training data generated by the widely used Spalart–Allmaras(SA) turbulence model before the Optimal Brain Surgeon(OBS) is employed to determine the relevancy of input features.Using the selected relevant features, a fully connected ANN model is constructed. The performance of the developed ANN model is numerically tested in the framework of DG for RANS, where the‘‘DG+ANN' method provides robust and steady convergence compared to the ‘‘DG+SA' method. The results demonstrate the promising potential to develop a general turbulence model based on artificial intelligence in the future given the training data covering a large rang of flow conditions.展开更多
Discontinuous Galerkin (DG) methods are a class of finite element methodsusing discontinuous basis functions, which are usually chosen as piecewise polynomi-als. Since the basis functions can be discontinuous, these m...Discontinuous Galerkin (DG) methods are a class of finite element methodsusing discontinuous basis functions, which are usually chosen as piecewise polynomi-als. Since the basis functions can be discontinuous, these methods have the flexibilitywhich is not shared by typical finite element methods, such as the allowance of ar-bitrary triangulation with hanging nodes, less restriction in changing the polynomialdegrees in each element independent of that in the neighbors (p adaptivity), and localdata structure and the resulting high parallel efficiency. In this paper, we give a generalreview of the local DG (LDG) methods for solving high-order time-dependent partialdifferential equations (PDEs). The important ingredient of the design of LDG schemes,namely the adequate choice of numerical fluxes, is highlighted. Some of the applica-tions of the LDG methods for high-order time-dependent PDEs are also be discussed.展开更多
The initial value problems and the first boundary problems for the quasilinear wave equation u_(tt)-[a_0+na_1(u_x)^(n-1)]u_(xx)-a_2u_(xxtt)=0 are considered,where a_0,a_2>0 are constants,a_1 is an arbitrary real nu...The initial value problems and the first boundary problems for the quasilinear wave equation u_(tt)-[a_0+na_1(u_x)^(n-1)]u_(xx)-a_2u_(xxtt)=0 are considered,where a_0,a_2>0 are constants,a_1 is an arbitrary real number,n is a natural number.The existence and uniqueness of the classical solutions for the initial value problems and the first boundary problems of the equation (1) are proved by the Galerkin method.展开更多
In this paper a hybridized weak Galerkin(HWG) finite element method for solving the Stokes equations in the primary velocity-pressure formulation is introduced.The WG method uses weak functions and their weak derivati...In this paper a hybridized weak Galerkin(HWG) finite element method for solving the Stokes equations in the primary velocity-pressure formulation is introduced.The WG method uses weak functions and their weak derivatives which are defined as distributions.Weak functions and weak derivatives can be approximated by piecewise polynomials with various degrees.Different combination of polynomial spaces leads to different WG finite element methods,which makes WG methods highly flexible and efficient in practical computation.A Lagrange multiplier is introduced to provide a numerical approximation for certain derivatives of the exact solution.With this new feature,the HWG method can be used to deal with jumps of the functions and their flux easily.Optimal order error estimates are established for the corresponding HWG finite element approximations for both primal variables and the Lagrange multiplier.A Schur complement formulation of the HWG method is derived for implementation purpose.The validity of the theoretical results is demonstrated in numerical tests.展开更多
By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IE...By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IEFG) method is used for 2D transient heat conduction problems, and the finite difference method is applied in the splitting direction. The discretized system equation is obtained based on the Galerkin weak form of 2D problem; the essential boundary conditions are imposed with the penalty method; and the finite difference method is employed in the time domain. Four exemplary problems are chosen to verify the efficiency of the DSEFG method. The numerical solutions show that the efficiency and precision of the DSEFG method are greater than ones of the IEFG method for 3D problems.展开更多
文摘In this paper we study the higher accuracy methods - the extrapolation and defect correction for the semidiscrete Galerkin approximations to the solutions of Sobolev and viscoelasticity type equations. The global extrapolation and the correction approximations of third order, rather than the pointwise extrapolation results are presented.
基金supported by National Natural Science Foundation of China (Grant No. 50875231)Hebei Provincial Major Natural Science Foundation of China (Grant No. E2006001038)
文摘The shape and thickness qualities of strip are influenced by the vibration of rolling mill. At present, the researches on the vibration of rolling mill are mainly the vertical vibration and torsional vibration of single stand mill, the study on the vibration of tandem rolling mill is rare. For the vibration of tandem rolling mill, the key problem is the vibration of the moving strip between stands. In this paper, considering the dynamic of moving strip and rolling theory, the vertical vibration of moving strip in the rolling process was proposed. Take the moving strip between the two mills of tandem rolling mill in the rolling process as subject investigated, according to the theory of moving beam, the vertical vibration model of moving strip in the rolling process was established. The partial differential equation was discretized by Galerkin truncation method. The natural frequency and stability of the moving strip were investigated and the numerical simulation in time domain was made. Simulation results show that, the natural frequency was strongly influenced by the rolling velocity and tension. With increasing of the rolling velocity, the first three natural frequencies decrease, the fourth natural frequency increases; with increasing of the unit tension, when the rolling velocity is high and low, respectively, the low order dimensionless natural frequency gradually decreases and increases, respectively. According to the stability of moving strip, the critical speed was determined, and the matching relationship of the tension and rolling velocity was also determined. This model can be used to study the stability of moving strip, improve the quality of strip and develop new rolling technology from the aspect of dynamics.
基金This work is supported in part by the National Natural Science Foundation of China （10571053）Program for New Century Excellent Talents in University, and the Scientific Research Fund of Hunan Provincial Education Department （0513039） The second author is supported in part by the US National Science Foundation under grants DMS-0311807 and DMS-0612908
文摘The convergence and superconvergence properties of the discontinuous Galerkin （DG） method for a singularly perturbed model problem in one-dimensional setting are studied. By applying the DG method with appropriately chosen numerical traces, the existence and uniqueness of the DG solution, the optimal order L2 error bounds, and 2p＋ 1-order superconvergence of the numerical traces are established. The numerical results indicate that the DG method does not produce any oscillation even under the uniform mesh. Numerical experiments demonstrate that, under the uniform mesh, it seems impossible to obtain the uniform superconvergence of the numerical traces. Nevertheless, thanks to the implementation of the so-called Shishkin-type mesh, the uniform 2p ＋ 1-order superconvergence is observed numerically.
基金the National Natural Science Foundation of China (Grant No. 11171208)Shanghai Leading Academic Discipline Project (Grant No. S30106)
文摘With the improved moving least-squares (IMLS) approximation, an orthogonal function system with a weight function is used as the basis function. The combination of the element-free Galerkin (EFG) method and the IMLS approximation leads to the development of the improved element-free Galerkin (IEFG) method. In this paper, the IEFG method is applied to study the partial differential equations that control the heat flow in three-dimensional space. With the IEFG technique, the Galerkin weak form is employed to develop the discretized system equations, and the penalty method is applied to impose the essential boundary conditions. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the transient heat conduction equations and the boundary and initial conditions are time dependent, the scaling parameter, number of nodes and time step length are considered in a convergence study.
基金supported by the National Natural Science Foundation of China (11171208)Shanghai Leading Academic Discipline Project (S30106)
文摘The paper presents the improved element-free Galerkin （IEFG） method for three-dimensional wave propa- gation. The improved moving least-squares （IMLS） approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares （MLS） approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.
基金National Natural Science Foundation of China (10272012)
文摘The aeroelastic analysis of high-altitude, long-endurance （HALE） aircraft that features high-aspect-ratio flexible wings needs take into account structural geometrical nonlinearities and dynamic stall. For a generic nonlinear aeroelastic system, besides the stability boundary, the characteristics of the limit-cycle oscillation （LCO） should also be accurately predicted. In order to conduct nonlinear aeroelastic analysis of high-aspect-ratio flexible wings, a first-order, state-space model is developed by combining a geometrically exact, nonlinear anisotropic beam model with nonlinear ONERA （Edlin） dynamic stall model. The present investigations focus on the initiation and sustaining mechanism of the LCO and the effects of flight speed and drag on aeroelastic behaviors. Numerical results indicate that structural geometrical nonlinearities could lead to the LCO without stall occurring. As flight speed increases, dynamic stall becomes dominant and the LCO increasingly complicated. Drag could be negligible for LCO type, but should be considered to exactly predict the onset speed of flutter or LCO of high-aspect-ratio flexible wings.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.
基金Supported by the National Natural Science Foundation of China (10571053,10871066)the Programme for New Century Excellent Talents in University (NCET-06-0712).
文摘In this paper, we consider the local discontinuous Galerkin method （LDG） for solving singularly perturbed convection-diffusion problems in one- and two-dimensional settings. The existence and uniqueness of the LDG solutions are verified. Numerical experiments demonstrate that it seems impossible to obtain uniform superconvergence for numerical fluxes under uniform meshes. Thanks to the implementation of two-type different anisotropic meshes, i.e., the Shishkin and an improved grade meshes, the uniform 2p ＋ i-order superconvergence is observed numerically for both one-dimensional and twodimensional cases.
基金Acknowledgments. The authors express their thanks to the referees for the comments and constructive suggestions, which were valuable in improving the quality of the manuscript. This work is supported by the National Natural Science Foundation of China（11172192） and the National Natural Science Pre-Research Foundation of Soochow University （SDY2011B01）.
文摘In this paper, a relaxed Hermitian and skew-Hermitian splitting （RHSS） preconditioner is proposed for saddle point problems from the element-free Galerkin （EFG） discretization method. The EFG method is one of the most widely used meshfree methods for solving partial differential equations. The RHSS preconditioner is constructed much closer to the coefficient matrix than the well-known HSS preconditioner, resulting in a RHSS fixed-point iteration. Convergence of the RHSS iteration is analyzed and an optimal parameter, which minimizes the spectral radius of the iteration matrix is described. Using the RHSS pre- conditioner to accelerate the convergence of some Krylov subspace methods （like GMRES） is also studied. Theoretical analyses show that the eigenvalues of the RHSS precondi- tioned matrix are real and located in a positive interval. Eigenvector distribution and an upper bound of the degree of the minimal polynomial of the preconditioned matrix are obtained. A practical parameter is suggested in implementing the RHSS preconditioner. Finally, some numerical experiments are illustrated to show the effectiveness of the new preconditioner.
基金Project (51504044) supported by the National Natural Science Foundation of ChinaProject (14KF05) supported by the Research Fund of the State Key Laboratory of Coal Resources and Mine Safety(CUMT),China+2 种基金Project (2015CDJXY) supported by the Fundamental Research Funds for the Central Universities,ChinaProject (2015M570607) supported by Postdoctoral Science FoundationProject (2011DA105287-MS201503) supported by the Independent Subject of State Key Laboratory of Coal Mine Disaster Dynamics and Control,China
文摘Irregular plates are very common structures in engineering,such as ore structures in mining.In this work,the Galerkin solution to the problem of a Kirchhoff plate lying on the Winkler foundation with two edges simply supported and the other two clamped supported is derived.Coordinate transformation technique is used during the solving process so that the solution is suitable to irregular shaped plates.The mechanical model and the solution proposed are then used to model the crown pillars between two adjacent levels in Sanshandao gold mine,which uses backfill method for mining operation.After that,an objective function,which takes security,economic profits and filling effect into consideration,is built to evaluate design proposals.Thickness optimizations for crown pillars are finally conducted in both conditions that the vertical stiffness of the foundation is known and unknown.The procedure presented in the work provides the guidance in thickness designing of complex shaped crown pillars and the preparation of backfill materials,thus to achieve the best balance between security and profits.
文摘This paper proposes and analyzes a class of robust globally divergence-free weak Galerkin （WG） finite element methods for Stokes equations. The new methods use the Pk/Pk-1 （k ≥ 1） discontinuous finite element combination for velocity and pressure in the interior of elements, and piecewise P1/Pk （l = k - 1, k） for the trace approximations of the ve- locity and pressure on the inter-element boundaries. Our methods not only yield globally divergence-free velocity solutions, but also have uniform error estimates with respect to the Reynolds number. Numerical experiments are provided to show the robustness of the proposed methods.
基金co-supported by the Aeronautical Science Foundation of China (Nos. 20151452021and 20152752033)the National Natural Science Foundation of China (No. 61732006)
文摘High-order Discontinuous Galerkin(DG) methods have been receiving more and more attentions in the area of Computational Fluid Dynamics(CFD) because of their high-accuracy property. However, it is still a challenge to obtain converged solution rapidly when solving the Reynolds Averaged Navier–Stokes(RANS) equations since the turbulence models significantly increase the nonlinearity of discretization system. The overall goal of this research is to develop an Artificial Neural Networks(ANNs) model with low complexity acting as an algebraic turbulence model to estimate the turbulence eddy viscosity for RANS. The ANN turbulence model is off-line trained using the training data generated by the widely used Spalart–Allmaras(SA) turbulence model before the Optimal Brain Surgeon(OBS) is employed to determine the relevancy of input features.Using the selected relevant features, a fully connected ANN model is constructed. The performance of the developed ANN model is numerically tested in the framework of DG for RANS, where the‘‘DG+ANN' method provides robust and steady convergence compared to the ‘‘DG+SA' method. The results demonstrate the promising potential to develop a general turbulence model based on artificial intelligence in the future given the training data covering a large rang of flow conditions.
基金The research of the first author is support by NSFC grant 10601055,FANEDD of CAS and SRF for ROCS SEMThe research of the second author is supported by NSF grant DMS-0809086 and DOE grant DE-FG02-08ER25863.
文摘Discontinuous Galerkin (DG) methods are a class of finite element methodsusing discontinuous basis functions, which are usually chosen as piecewise polynomi-als. Since the basis functions can be discontinuous, these methods have the flexibilitywhich is not shared by typical finite element methods, such as the allowance of ar-bitrary triangulation with hanging nodes, less restriction in changing the polynomialdegrees in each element independent of that in the neighbors (p adaptivity), and localdata structure and the resulting high parallel efficiency. In this paper, we give a generalreview of the local DG (LDG) methods for solving high-order time-dependent partialdifferential equations (PDEs). The important ingredient of the design of LDG schemes,namely the adequate choice of numerical fluxes, is highlighted. Some of the applica-tions of the LDG methods for high-order time-dependent PDEs are also be discussed.
文摘The initial value problems and the first boundary problems for the quasilinear wave equation u_(tt)-[a_0+na_1(u_x)^(n-1)]u_(xx)-a_2u_(xxtt)=0 are considered,where a_0,a_2>0 are constants,a_1 is an arbitrary real number,n is a natural number.The existence and uniqueness of the classical solutions for the initial value problems and the first boundary problems of the equation (1) are proved by the Galerkin method.
基金supported by National Natural Science Foundation of China(Grant Nos.11271157,11371171 and 11471141)the Program for New Century Excellent Talents in University of Ministry of Education of China
文摘In this paper a hybridized weak Galerkin(HWG) finite element method for solving the Stokes equations in the primary velocity-pressure formulation is introduced.The WG method uses weak functions and their weak derivatives which are defined as distributions.Weak functions and weak derivatives can be approximated by piecewise polynomials with various degrees.Different combination of polynomial spaces leads to different WG finite element methods,which makes WG methods highly flexible and efficient in practical computation.A Lagrange multiplier is introduced to provide a numerical approximation for certain derivatives of the exact solution.With this new feature,the HWG method can be used to deal with jumps of the functions and their flux easily.Optimal order error estimates are established for the corresponding HWG finite element approximations for both primal variables and the Lagrange multiplier.A Schur complement formulation of the HWG method is derived for implementation purpose.The validity of the theoretical results is demonstrated in numerical tests.
基金supported by the National Natural Science Foundation of China(Grant Nos.11571223,and 51404160)the Science and Technology Innovation Foundation of Higher Education of Shanxi Province(Grant No.2016163)
文摘By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IEFG) method is used for 2D transient heat conduction problems, and the finite difference method is applied in the splitting direction. The discretized system equation is obtained based on the Galerkin weak form of 2D problem; the essential boundary conditions are imposed with the penalty method; and the finite difference method is employed in the time domain. Four exemplary problems are chosen to verify the efficiency of the DSEFG method. The numerical solutions show that the efficiency and precision of the DSEFG method are greater than ones of the IEFG method for 3D problems.
