In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolso...In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.展开更多
In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical ...In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical reaction,say zi,is between 0 and 1,withΣz_(i)=1.Due to the lack of maximum-principle,most of the previous bound-preserving technique cannot be applied directly.To preserve those bounds,we will use the positivity-preserving technique to all the zi'is and enforceΣz_(i)=1 by constructing conservative schemes,thanks to conservative time integrations and consistent numerical fluxes in the system.Moreover,detonation is an extreme singular mode of flame propagation in premixed gas,and the model contains a significant stiff source.It is well known that for hyperbolic equations with stiff source,the transition points in the numerical approximations near the shocks may trigger spurious shock speed,leading to wrong shock position.Intuitively,the high-order weighted essentially non-oscillatory(WENO)scheme,which can suppress oscillations near the discontinuities,would be a good choice for spatial discretization.However,with the nonlinear weights,the numerical fluxes are no longer“consistent”,leading to nonconservative numerical schemes and the bound-preserving technique does not work.Numerical experiments demonstrate that,without further numerical techniques such as subcell resolutions,the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.展开更多
A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the c...A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.展开更多
This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip veloci...This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.展开更多
The possibility of using a nodal method allowing irregular distribution of nodes in a natural way is one of the main advantages of the generalized finite difference method (GFDM) with regard to the classical finite di...The possibility of using a nodal method allowing irregular distribution of nodes in a natural way is one of the main advantages of the generalized finite difference method (GFDM) with regard to the classical finite difference method. Moreover, this feature has made it one of the most-promising meshless methods because it also allows us to reduce the time-consuming task of mesh generation and the numerical solution of integrals. This characteristic allows us to shape geological features easily whilst maintaining accuracy in the results, which can be a source of great interest when dealing with this kind of problems. Two widespread geophysical investigation methods in civil engineering are the cross-hole method and the seismic refraction method. This paper shows the use of the GFDM to model the aforementioned geophysical investigation tests showing precision in the obtained results when comparing them with experimental data.展开更多
An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite...An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.展开更多
A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite differen...A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.展开更多
In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot's poroelastic equations will be e...In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot's poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.展开更多
An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation (FPE) with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for ...An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation (FPE) with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.展开更多
In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs...In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.展开更多
A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is ...A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.展开更多
Thermal stress simulation can provide a scientific reference to eliminate defects such as crack,residual stress centralization and deformation etc.,caused by thermal stress during casting solidification.To study the t...Thermal stress simulation can provide a scientific reference to eliminate defects such as crack,residual stress centralization and deformation etc.,caused by thermal stress during casting solidification.To study the thermal stress distribution during casting process,a unilateral thermal-stress coupling model was employed to simulate 3D casting stress using Finite Difference Method(FDM),namely all the traditional thermal-elastic-plastic equations are numerically and differentially discrete.A FDM/FDM numerical simulation system was developed to analyze temperature and stress fields during casting solidification process.Two practical verifications were carried out,and the results from simulation basically coincided with practical cases.The results indicated that the FDM/FDM stress simulation system can be used to simulate the formation of residual stress,and to predict the occurrence of hot tearing.Because heat transfer and stress analysis are all based on FDM,they can use the same FD model,which can avoid the matching process between different models,and hence reduce temperature-load transferring errors.This approach makes the simulation of fluid flow,heat transfer and stress analysis unify into one single model.展开更多
In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Di...In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Differential Equations (PDE) with mixed derivatives along with Dirichlet and Neumann boundary conditions. These two methods have almost the same accuracy from theoretical aspect with regular boundaries, but generally Finite Element Method produces better approximations when the boundaries are irregular. In order to investigate which method produces better results from numerical aspect, we apply these methods into specific examples with regular boundaries with constant step-size for both of them. The results which obtained confirm, in most of the cases, the theoretical results.展开更多
A numerical method has been developed to extract the composition-dependent interdiffusivity from the concentration profiles in the aluminide coating prepared by pack cementation. The procedure is based on the classic ...A numerical method has been developed to extract the composition-dependent interdiffusivity from the concentration profiles in the aluminide coating prepared by pack cementation. The procedure is based on the classic finite difference method (FDM). In order to simplify the model, effect of some alloying elements on interdiffusivity can be negligible. Calculated results indicate the interdiffusivity in aluminide coating strongly depends on the composition and give the formulas used to calculate interdiffusivity at 850, 950 and 1050癈. The effect on interdiffusivity is briefly discussed.展开更多
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.展开更多
Surface subsidence induced by underground mining is a typical serious geohazard.Numerical approaches such as the discrete element method(DEM)and finite difference method(FDM)have been widely used to model and analyze ...Surface subsidence induced by underground mining is a typical serious geohazard.Numerical approaches such as the discrete element method(DEM)and finite difference method(FDM)have been widely used to model and analyze mining-induced surface subsidence.However,the DEM is typically computationally expensive,and is not capable of analyzing large-scale problems,while the mesh distortion may occur in the FDM modeling of largely deformed surface subsidence.To address the above problems,this paper presents a geometrically and locally adaptive remeshing method for the FDM modeling of largely deformed surface subsidence induced by underground mining.The essential ideas behind the proposed method are as follows:(i)Geometrical features of elements(i.e.the mesh quality),rather than the calculation errors,are employed as the indicator for determining whether to conduct the remeshing;and(ii)Distorted meshes with multiple attributes,rather than those with only a single attribute,are locally regenerated.In the proposed method,the distorted meshes are first adaptively determined based on the mesh quality,and then removed from the original mesh model.The tetrahedral mesh in the distorted area is first regenerated,and then the physical field variables of old mesh are transferred to the new mesh.The numerical calculation process recovers when finishing the regeneration and transformation.To verify the effectiveness of the proposed method,the surface deformation of the Yanqianshan iron mine,Liaoning Province,China,is numerically investigated by utilizing the proposed method,and compared with the numerical results of the DEM modeling.Moreover,the proposed method is applied to predicting the surface subsidence in Anjialing No.1 Underground Mine,Shanxi Province,China.展开更多
With the finite difference method to calculate the temperature distribution in mass concrete structures, the solution precision will increase with a smaller step size, at the cost of computational time. In view of the...With the finite difference method to calculate the temperature distribution in mass concrete structures, the solution precision will increase with a smaller step size, at the cost of computational time. In view of the basic characteristics of the finite difference method, a simple yet powerful improvement is introduced. By multiplying the adiabatic temperature function with a correction factor, the precision of the solution can be assured without an increase in the computation time. In addition, the correction rules for three types of commonly used concrete hydration formulas are investigated.展开更多
A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. T...A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.展开更多
In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formul...In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formula,a one-to-one correspondence occurs between the physical and spurious states.Although these energies are exactly the same,the wave functions of the spurious states exhibit a much faster staggering than those of the physical states.Such spurious states can be eliminated when applying the finite difference method by introducing an extra Wilson term into the Hamiltonian.Furthermore,it was also found that the number of spurious states can be reduced if we improve the accuracy of the numerical differential formula.The Dirac equation is then solved in a momentum space in which there is no differential operator,and we found that the spurious states can be completely avoided in the momentum space,even without an extra Wilson term.展开更多
This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and th...This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart (CR) element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.展开更多
基金supported by the Key Laboratory of Road Construction Technology and Equipment(Chang’an University,No.300102253502)the Natural Science Foundation of Shandong Province of China(GrantNo.ZR2022YQ06)the Development Plan of Youth Innovation Team in Colleges and Universities of Shandong Province(Grant No.2022KJ140).
文摘In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.
基金the National Natural Science Foundation of China under Grant Number NSFC 11801302Tsinghua University Initiative Scientific Research Program.Yang Yang is supported by the NSF Grant DMS-1818467.
文摘In this paper,we apply high-order finite difference(FD)schemes for multispecies and multireaction detonations(MMD).In MMD,the density and pressure are positive and the mass fraction of the ith species in the chemical reaction,say zi,is between 0 and 1,withΣz_(i)=1.Due to the lack of maximum-principle,most of the previous bound-preserving technique cannot be applied directly.To preserve those bounds,we will use the positivity-preserving technique to all the zi'is and enforceΣz_(i)=1 by constructing conservative schemes,thanks to conservative time integrations and consistent numerical fluxes in the system.Moreover,detonation is an extreme singular mode of flame propagation in premixed gas,and the model contains a significant stiff source.It is well known that for hyperbolic equations with stiff source,the transition points in the numerical approximations near the shocks may trigger spurious shock speed,leading to wrong shock position.Intuitively,the high-order weighted essentially non-oscillatory(WENO)scheme,which can suppress oscillations near the discontinuities,would be a good choice for spatial discretization.However,with the nonlinear weights,the numerical fluxes are no longer“consistent”,leading to nonconservative numerical schemes and the bound-preserving technique does not work.Numerical experiments demonstrate that,without further numerical techniques such as subcell resolutions,the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.
基金the National Natural Science Foundation of China(No.11701103)the Young Top-notch Talent Program of Guangdong Province of China(No.2017GC010379)+4 种基金the Natural Science Foundation of Guangdong Province of China(No.2022A1515012147)the Project of Science and Technology of Guangzhou of China(No.202102020704)the Opening Project of Guangdong Province Key Laboratory of Computational Science at the Sun Yat-sen University of China(2021023)the Science and Technology Development Fund,Macao SAR(File No.0005/2019/A)the University of Macao of China(File Nos.MYRG2020-00035-FST,MYRG2018-00047-FST).
文摘A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.
文摘This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.
基金The authors acknowledge the support of the Escuela Tecnica Superior de Ingenieros Industriales(UNED)of Spain,project 2019-IFC02of the Universidad Politecnica de Madrid(UPM)(Research groups 2019).
文摘The possibility of using a nodal method allowing irregular distribution of nodes in a natural way is one of the main advantages of the generalized finite difference method (GFDM) with regard to the classical finite difference method. Moreover, this feature has made it one of the most-promising meshless methods because it also allows us to reduce the time-consuming task of mesh generation and the numerical solution of integrals. This characteristic allows us to shape geological features easily whilst maintaining accuracy in the results, which can be a source of great interest when dealing with this kind of problems. Two widespread geophysical investigation methods in civil engineering are the cross-hole method and the seismic refraction method. This paper shows the use of the GFDM to model the aforementioned geophysical investigation tests showing precision in the obtained results when comparing them with experimental data.
文摘An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.
基金the National Natural Science Foundation of China
文摘A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.
基金Project supported by the "100 Talents Project" of the Chinese Academy of Sciences and the Major Program of the National Natural Science Foundation of China (Grant No 10534040).
文摘In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot's poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.
基金supported by the National Natural Science Foundation of China(Nos.11171193 and11371229)the Natural Science Foundation of Shandong Province(No.ZR2014AM033)the Science and Technology Development Project of Shandong Province(No.2012GGB01198)
文摘An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation (FPE) with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.
文摘In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.
基金supported by the Yunnan Provincial Applied Basic Research Program of China(No. KKSY201207019)
文摘A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.
基金supported by the National Natural Science Foundation of China (No.50805056)New Century Excellent Talents in University (No.NCET-09-0396)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education (2009)
文摘Thermal stress simulation can provide a scientific reference to eliminate defects such as crack,residual stress centralization and deformation etc.,caused by thermal stress during casting solidification.To study the thermal stress distribution during casting process,a unilateral thermal-stress coupling model was employed to simulate 3D casting stress using Finite Difference Method(FDM),namely all the traditional thermal-elastic-plastic equations are numerically and differentially discrete.A FDM/FDM numerical simulation system was developed to analyze temperature and stress fields during casting solidification process.Two practical verifications were carried out,and the results from simulation basically coincided with practical cases.The results indicated that the FDM/FDM stress simulation system can be used to simulate the formation of residual stress,and to predict the occurrence of hot tearing.Because heat transfer and stress analysis are all based on FDM,they can use the same FD model,which can avoid the matching process between different models,and hence reduce temperature-load transferring errors.This approach makes the simulation of fluid flow,heat transfer and stress analysis unify into one single model.
文摘In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Differential Equations (PDE) with mixed derivatives along with Dirichlet and Neumann boundary conditions. These two methods have almost the same accuracy from theoretical aspect with regular boundaries, but generally Finite Element Method produces better approximations when the boundaries are irregular. In order to investigate which method produces better results from numerical aspect, we apply these methods into specific examples with regular boundaries with constant step-size for both of them. The results which obtained confirm, in most of the cases, the theoretical results.
文摘A numerical method has been developed to extract the composition-dependent interdiffusivity from the concentration profiles in the aluminide coating prepared by pack cementation. The procedure is based on the classic finite difference method (FDM). In order to simplify the model, effect of some alloying elements on interdiffusivity can be negligible. Calculated results indicate the interdiffusivity in aluminide coating strongly depends on the composition and give the formulas used to calculate interdiffusivity at 850, 950 and 1050癈. The effect on interdiffusivity is briefly discussed.
文摘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 National Natural Science Foundation of China(Grant Nos.11602235 and 41772326)the Fundamental Research Funds for the Central Universities of China(Grant No.2652018091)。
文摘Surface subsidence induced by underground mining is a typical serious geohazard.Numerical approaches such as the discrete element method(DEM)and finite difference method(FDM)have been widely used to model and analyze mining-induced surface subsidence.However,the DEM is typically computationally expensive,and is not capable of analyzing large-scale problems,while the mesh distortion may occur in the FDM modeling of largely deformed surface subsidence.To address the above problems,this paper presents a geometrically and locally adaptive remeshing method for the FDM modeling of largely deformed surface subsidence induced by underground mining.The essential ideas behind the proposed method are as follows:(i)Geometrical features of elements(i.e.the mesh quality),rather than the calculation errors,are employed as the indicator for determining whether to conduct the remeshing;and(ii)Distorted meshes with multiple attributes,rather than those with only a single attribute,are locally regenerated.In the proposed method,the distorted meshes are first adaptively determined based on the mesh quality,and then removed from the original mesh model.The tetrahedral mesh in the distorted area is first regenerated,and then the physical field variables of old mesh are transferred to the new mesh.The numerical calculation process recovers when finishing the regeneration and transformation.To verify the effectiveness of the proposed method,the surface deformation of the Yanqianshan iron mine,Liaoning Province,China,is numerically investigated by utilizing the proposed method,and compared with the numerical results of the DEM modeling.Moreover,the proposed method is applied to predicting the surface subsidence in Anjialing No.1 Underground Mine,Shanxi Province,China.
基金the Shanghai Excellent Young Professor Candidate Program (Grant No.04YQHB139)
文摘With the finite difference method to calculate the temperature distribution in mass concrete structures, the solution precision will increase with a smaller step size, at the cost of computational time. In view of the basic characteristics of the finite difference method, a simple yet powerful improvement is introduced. By multiplying the adiabatic temperature function with a correction factor, the precision of the solution can be assured without an increase in the computation time. In addition, the correction rules for three types of commonly used concrete hydration formulas are investigated.
文摘A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.
基金partly supported by the National Natural Science Foundation of China(No.11875070)the Natural Science Foundation of Anhui Province(No.1908085MA16)
文摘In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formula,a one-to-one correspondence occurs between the physical and spurious states.Although these energies are exactly the same,the wave functions of the spurious states exhibit a much faster staggering than those of the physical states.Such spurious states can be eliminated when applying the finite difference method by introducing an extra Wilson term into the Hamiltonian.Furthermore,it was also found that the number of spurious states can be reduced if we improve the accuracy of the numerical differential formula.The Dirac equation is then solved in a momentum space in which there is no differential operator,and we found that the spurious states can be completely avoided in the momentum space,even without an extra Wilson term.
基金supported by the National Natural Science Foundation of China(Nos.11271273 and 11271298)
文摘This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart (CR) element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.