This paper aims to investigate the tamed Euler method for the random periodic solution of semilinear SDEs with one-sided Lipschitz coefficient.We introduce a novel approach to analyze mean-square error bounds of the n...This paper aims to investigate the tamed Euler method for the random periodic solution of semilinear SDEs with one-sided Lipschitz coefficient.We introduce a novel approach to analyze mean-square error bounds of the novel schemes,without relying on a priori high-order moment bound of the numerical approximation.The expected order-one mean square convergence is attained for the proposed scheme.Moreover,a numerical example is presented to verify our theoretical analysis.展开更多
This paper presents a method for tracing a planar implicit curve f（x, y）=0 on a rectangular region based on continuation scheme. First, according to the starting track-point and the starting track-direction of the c...This paper presents a method for tracing a planar implicit curve f（x, y）=0 on a rectangular region based on continuation scheme. First, according to the starting track-point and the starting track-direction of the curve, make a new fimction F（x, y）=0 where the same curve withf（x, y）=0 is defined. Then we trace the curve between the two domains where F（x, y）〉0 and F（x, y）〈0 alternately, according to the two rules presented in this paper. Equal step size or adaptive step size can be used, when we trace the curve. An irregular planar implicit curve （such as the curve with large curvatures at some points on the curve）, can be plotted if an adaptive step size is used. Moreover, this paper presents a scheme to search for the multiple points on the curve. Our method has the following advantages： （1） it can plot Co planar implicit curves; （2） it can plot the planar implicit curves with multiple points; （3） by the help of using the two rules, our method does not need to compute the tangent vector at the points on the curve, and directly searches for the direction of the tracing curve; （4） the tracing procedure costs only one of two evaluations of function f（x, y）=0 per moving step, while most existing similar methods cost more evaluations of the function.展开更多
A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first ...A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.展开更多
In this paper,we present the semi-implicit Euler（SIE）numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square...In this paper,we present the semi-implicit Euler（SIE）numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square sense under the Local Lipschitz condition.展开更多
We deal with the numerical solution of the compressible Euler equations with the aid of the discontinuous Galerkin(DG)method with focus on the goal-oriented error estimates and adaptivity.We analyse the adjoint consis...We deal with the numerical solution of the compressible Euler equations with the aid of the discontinuous Galerkin(DG)method with focus on the goal-oriented error estimates and adaptivity.We analyse the adjoint consistency of the DG scheme where the adjoint problem is not formulated by the differentiation of the DG form and the target functional but using a suitable linearization of the nonlinear forms.Furthermore,we present the goal-oriented anisotropic hp-mesh adaptation method for the Euler equations.The theoretical results are supported by numerical experiments.展开更多
This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called...This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.展开更多
为高效模拟电磁散射问题,提出一种积分方程(integral equation,IE)辐射边界条件结合频域有限体积(finite volume frequency domain,FVFD)的电磁场计算方法。FVFD采用全隐格式时间推进和通量迎风分裂计算,IE算法在网格内电流等效面上积...为高效模拟电磁散射问题,提出一种积分方程(integral equation,IE)辐射边界条件结合频域有限体积(finite volume frequency domain,FVFD)的电磁场计算方法。FVFD采用全隐格式时间推进和通量迎风分裂计算,IE算法在网格内电流等效面上积分计算截断边界输入波通量,两者以松耦合方式结合。数值实验表明,IE算法对比常用的Mur算法,截断边界到物面距离可缩短到0.5倍波长,模拟多体散射问题可使用整体网格或独立的分离网格,也能模拟复杂外形目标电磁散射。所提方法在保证计算精度同时,可降低网格规模,有效提升效率。展开更多
Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix d...Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.展开更多
Natural gas hydrate, as a potential energy resource, deposits in permafrost and marine sediment with large quantities. The current exploitation methods include depressurization, thermal stimulation, and inhibitor inje...Natural gas hydrate, as a potential energy resource, deposits in permafrost and marine sediment with large quantities. The current exploitation methods include depressurization, thermal stimulation, and inhibitor injection. However, many issues have to be resolved before the commercial production. In the present study, a 2-D axisymmetric simulator for gas production from hydrate reservoirs is developed. The simulator includes equations of conductive and convective heat transfer, kinetic of hydrate decomposition, and multiphase flow. These equations are discretized based on the finite difference method and are solved with the fully implicit simultaneous solution method. The process of laboratory-scale hydrate decomposition by depressurization is simulated. For different surrounding temperatures and outlet pressures, time evolutions of gas and water generations during hydrate dissociation are evaluated, and variations of temperature, pressure, and multiphase fluid flow conditions are analyzed. The results suggest that the rate of heat transfer plays an important role in the process. Furthermore, high surrounding temperature and low outlet valve pressure may increase the rate of hydrate dissociation with insignificant impact on final cumulative gas volume.展开更多
This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using prec...This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.展开更多
In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian co...In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.展开更多
This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit m...This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.展开更多
This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditiona...This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.展开更多
This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods...This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods always lead to optimal convergence rates and have some other important features, especially the methods can be implemented parallelly.展开更多
Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the dr...Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the drilling boreholes. The purpose of this study is to use field geophysics to calculate the depth of mineral reserve. The study area is located 38 km from Zarand city called Jalalabad iron mine. In this study, gravimetric data were measured and mineral depth was calculated using the Euler method. 1314 readings have been performed in this area. The rocks of the region include volcanic and sedimentary. The source of the mineralization in the area is hydrothermal processes. After gravity measuring in the region, the data were corrected, then various methods such as anomalous map remaining in levels one and two, upward expansion, first and second-degree vertical derivatives, analytical method, and analytical signal were drawn, and finally, the depth of the deposit was estimated by Euler method. As a result, the depth of the mineral deposit was calculated to be between 20 and 30 meters on average.展开更多
A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme...A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.展开更多
Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme...Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.展开更多
The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF me...The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.展开更多
In this paper,two classes of Riesz space fractional partial differential equations including space-fractional and space-time-fractional ones are considered.These two models can be regarded as the generalization of the...In this paper,two classes of Riesz space fractional partial differential equations including space-fractional and space-time-fractional ones are considered.These two models can be regarded as the generalization of the classical wave equation in two space dimensions.Combining with the Crank-Nicolson method in temporal direction,efficient alternating direction implicit Galerkin finite element methods for solving these two fractional models are developed,respectively.The corresponding stability and convergence analysis of the numerical methods are discussed.Numerical results are provided to verify the theoretical analysis.展开更多
A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 a...A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.展开更多
基金supported by the National Natural Science Foundation of China(Nos.12471394,12371417)Natural Science Foundation of Changsha(No.kq2502101)。
文摘This paper aims to investigate the tamed Euler method for the random periodic solution of semilinear SDEs with one-sided Lipschitz coefficient.We introduce a novel approach to analyze mean-square error bounds of the novel schemes,without relying on a priori high-order moment bound of the numerical approximation.The expected order-one mean square convergence is attained for the proposed scheme.Moreover,a numerical example is presented to verify our theoretical analysis.
文摘This paper presents a method for tracing a planar implicit curve f（x, y）=0 on a rectangular region based on continuation scheme. First, according to the starting track-point and the starting track-direction of the curve, make a new fimction F（x, y）=0 where the same curve withf（x, y）=0 is defined. Then we trace the curve between the two domains where F（x, y）〉0 and F（x, y）〈0 alternately, according to the two rules presented in this paper. Equal step size or adaptive step size can be used, when we trace the curve. An irregular planar implicit curve （such as the curve with large curvatures at some points on the curve）, can be plotted if an adaptive step size is used. Moreover, this paper presents a scheme to search for the multiple points on the curve. Our method has the following advantages： （1） it can plot Co planar implicit curves; （2） it can plot the planar implicit curves with multiple points; （3） by the help of using the two rules, our method does not need to compute the tangent vector at the points on the curve, and directly searches for the direction of the tracing curve; （4） the tracing procedure costs only one of two evaluations of function f（x, y）=0 per moving step, while most existing similar methods cost more evaluations of the function.
文摘A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.
基金Supported by the NSF of the Higher Education Institutions of Jiangsu Province(10KJD110006)Supported by the grant of Jiangsu Institute of Education(Jsjy2009zd03)Supported by the Qing Lan Project of Jiangsu Province(2010)
文摘In this paper,we present the semi-implicit Euler（SIE）numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square sense under the Local Lipschitz condition.
基金Grant no.20-01074S of the Czech Science Foundation.
文摘We deal with the numerical solution of the compressible Euler equations with the aid of the discontinuous Galerkin(DG)method with focus on the goal-oriented error estimates and adaptivity.We analyse the adjoint consistency of the DG scheme where the adjoint problem is not formulated by the differentiation of the DG form and the target functional but using a suitable linearization of the nonlinear forms.Furthermore,we present the goal-oriented anisotropic hp-mesh adaptation method for the Euler equations.The theoretical results are supported by numerical experiments.
文摘This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.
文摘为高效模拟电磁散射问题,提出一种积分方程(integral equation,IE)辐射边界条件结合频域有限体积(finite volume frequency domain,FVFD)的电磁场计算方法。FVFD采用全隐格式时间推进和通量迎风分裂计算,IE算法在网格内电流等效面上积分计算截断边界输入波通量,两者以松耦合方式结合。数值实验表明,IE算法对比常用的Mur算法,截断边界到物面距离可缩短到0.5倍波长,模拟多体散射问题可使用整体网格或独立的分离网格,也能模拟复杂外形目标电磁散射。所提方法在保证计算精度同时,可降低网格规模,有效提升效率。
基金Project supported by the National Natural Science Foundation of China(Nos.60273048and60174023).
文摘Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.
基金supported by the National High Technology Research and Development Program of China(863 Program, Grant No.2006AA09A209-5)the National Natural Science Foundation of China (Key Program,Grant No.50736001)the Major Research Project of Ministry of Education of China (Grant No.306005)
文摘Natural gas hydrate, as a potential energy resource, deposits in permafrost and marine sediment with large quantities. The current exploitation methods include depressurization, thermal stimulation, and inhibitor injection. However, many issues have to be resolved before the commercial production. In the present study, a 2-D axisymmetric simulator for gas production from hydrate reservoirs is developed. The simulator includes equations of conductive and convective heat transfer, kinetic of hydrate decomposition, and multiphase flow. These equations are discretized based on the finite difference method and are solved with the fully implicit simultaneous solution method. The process of laboratory-scale hydrate decomposition by depressurization is simulated. For different surrounding temperatures and outlet pressures, time evolutions of gas and water generations during hydrate dissociation are evaluated, and variations of temperature, pressure, and multiphase fluid flow conditions are analyzed. The results suggest that the rate of heat transfer plays an important role in the process. Furthermore, high surrounding temperature and low outlet valve pressure may increase the rate of hydrate dissociation with insignificant impact on final cumulative gas volume.
文摘This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11261035 and 11171038)the Science Research Foundation of the Institute of Higher Education of Inner Mongolia Autonomous Region, China (Grant No. NJZZ12198)the Natural Science Foundation of Inner Mongolia Autonomous Region, China (Grant No. 2012MS0102)
文摘In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.
基金the National Natural Science Foundation of China(Nos.41807223 and 51908175)the Fundamental Research Funds for the Central Universities(No.B210202096)+1 种基金the Natural Science Foundation of Guangdong Province(No.2018A030310346)the Water Conservancy Science and Technology Innovation Project of Guangdong Province(No.2020-11),China。
文摘This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.
基金The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China （G2000048702, 2003CB716707）the National Science Fund for Distinguished Young Scholars （10025208）+1 种基金 the National Natural Science Foundation of China （Key Program） （10532040） the Research Fund for 0versea Chinese （10228028）.
文摘This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.
基金This work was supported by the National Natural Science Foundation of China
文摘This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods always lead to optimal convergence rates and have some other important features, especially the methods can be implemented parallelly.
文摘Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the drilling boreholes. The purpose of this study is to use field geophysics to calculate the depth of mineral reserve. The study area is located 38 km from Zarand city called Jalalabad iron mine. In this study, gravimetric data were measured and mineral depth was calculated using the Euler method. 1314 readings have been performed in this area. The rocks of the region include volcanic and sedimentary. The source of the mineralization in the area is hydrothermal processes. After gravity measuring in the region, the data were corrected, then various methods such as anomalous map remaining in levels one and two, upward expansion, first and second-degree vertical derivatives, analytical method, and analytical signal were drawn, and finally, the depth of the deposit was estimated by Euler method. As a result, the depth of the mineral deposit was calculated to be between 20 and 30 meters on average.
文摘A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using A USMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to rest this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.
文摘Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.
文摘The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.
基金supported by the Guangxi Natural Science Foundation[grant numbers 2018GXNSFBA281020,2018GXNSFAA138121]the Doctoral Starting up Foundation of Guilin University of Technology[grant number GLUTQD2016044].
文摘In this paper,two classes of Riesz space fractional partial differential equations including space-fractional and space-time-fractional ones are considered.These two models can be regarded as the generalization of the classical wave equation in two space dimensions.Combining with the Crank-Nicolson method in temporal direction,efficient alternating direction implicit Galerkin finite element methods for solving these two fractional models are developed,respectively.The corresponding stability and convergence analysis of the numerical methods are discussed.Numerical results are provided to verify the theoretical analysis.
基金the National Natural Science Foundation of China (Grant No.50579035)
文摘A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.
