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.展开更多
For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple boun...For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.展开更多
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.展开更多
In this paper,we propose a new conservative high-order semi-Lagrangian finite difference(SLFD)method to solve linear advection equation and the nonlinear Vlasov and BGK models.The finite difference scheme has better c...In this paper,we propose a new conservative high-order semi-Lagrangian finite difference(SLFD)method to solve linear advection equation and the nonlinear Vlasov and BGK models.The finite difference scheme has better computational flexibility by working with point values,especially when working with high-dimensional problems in an operator splitting setting.The reconstruction procedure in the proposed SLFD scheme is motivated from the SL finite volume scheme.In particular,we define a new sliding average function,whose cell averages agree with point values of the underlying function.By developing the SL finite volume scheme for the sliding average function,we derive the proposed SLFD scheme,which is high-order accurate,mass conservative and unconditionally stable for linear problems.The performance of the scheme is showcased by linear transport applications,as well as the nonlinear Vlasov-Poisson and BGK models.Furthermore,we apply the Fourier stability analysis to a fully discrete SLFD scheme coupled with diagonally implicit Runge-Kutta(DIRK)method when applied to a stiff two-velocity hyperbolic relaxation system.Numerical stability and asymptotic accuracy properties of DIRK methods are discussed in theoretical and computational aspects.展开更多
In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.Ho...In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.展开更多
We present a first-order finite difference scheme for approximating solutions of a mathematical model of cervical cancer induced by the human papillomavirus (HPV), which consists of four nonlinear partial differential...We present a first-order finite difference scheme for approximating solutions of a mathematical model of cervical cancer induced by the human papillomavirus (HPV), which consists of four nonlinear partial differential equations and a nonlinear first-order ordinary differential equation. The scheme is analyzed and used to provide an existence-uniqueness result. Numerical simulations are performed in order to demonstrate the first-order rate of convergence. A sensitivity analysis was done in order to compare the effects of two drug types, those that increase the death rate of HPV-infected cells, and those that increase the death rate of the precancerous cell population. The model predicts that treatments that affect the precancerous cell population by directly increasing the corresponding death rate are far more effective than those that increase the death rate of HPV-infected cells.展开更多
Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of t...Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference (FD) method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.展开更多
In this paper,a new type of finite difference mapped weighted essentially non-oscillatory(MWENO)schemes with unequal-sized stencils,such as the seventh-order and ninthorder versions,is constructed for solving hyperbol...In this paper,a new type of finite difference mapped weighted essentially non-oscillatory(MWENO)schemes with unequal-sized stencils,such as the seventh-order and ninthorder versions,is constructed for solving hyperbolic conservation laws.For the purpose of designing increasingly high-order finite difference WENO schemes,the equal-sized stencils are becoming more and more wider.The more we use wider candidate stencils,the bigger the probability of discontinuities lies in all stencils.Therefore,one innovation of these new WENO schemes is to introduce a new splitting stencil methodology to divide some fourpoint or five-point stencils into several smaller three-point stencils.By the usage of this new methodology in high-order spatial reconstruction procedure,we get different degree polynomials defined on these unequal-sized stencils,and calculate the linear weights,smoothness indicators,and nonlinear weights as specified in Jiang and Shu(J.Comput.Phys.126:202228,1996).Since the difference between the nonlinear weights and the linear weights is too big to keep the optimal order of accuracy in smooth regions,another crucial innovation is to present the new mapping functions which are used to obtain the mapped nonlinear weights and decrease the difference quantity between the mapped nonlinear weights and the linear weights,so as to keep the optimal order of accuracy in smooth regions.These new MWENO schemes can also be applied to compute some extreme examples,such as the double rarefaction wave problem,the Sedov blast wave problem,and the Leblanc problem with a normal CFL number.Extensive numerical results are provided to illustrate the good performance of the new finite difference MWENO schemes.展开更多
To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to...To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.展开更多
In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic fi...In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic finite difference scheme is presented for this equation.Some mathematical analyses of the scheme are then discussed.Lastly,to ascertain the efficacy and accuracy of the suggested technique,the numerical results are discussed and compared with the exact solution.展开更多
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.展开更多
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.展开更多
Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with ...Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.展开更多
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.展开更多
The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux ...The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.展开更多
In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is s...In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.展开更多
The dynamic calculations of slender marine risers, such as Finite Element Method (FEM) or Modal Expansion Solution Method (MESM), are mainly for the slender structures with their both ends hinged to the surface an...The dynamic calculations of slender marine risers, such as Finite Element Method (FEM) or Modal Expansion Solution Method (MESM), are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation (FDA) method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bemoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser's free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.展开更多
Basedonthe first-order Biot-equation with simplified coefficients,astaggered irregu- lar-grid finite difference method(FDM)is developed to simulate elastic wave propagation in 3-D heterogeneous anisotropic porous medi...Basedonthe first-order Biot-equation with simplified coefficients,astaggered irregu- lar-grid finite difference method(FDM)is developed to simulate elastic wave propagation in 3-D heterogeneous anisotropic porous media.The ‘slow’P wave in porous media wave simulation is highly dispersive.Finer grids are needed to get a precise wavefield calculation for models with curved interface and complex geometric structure.Fine grids in a global model greatly increase computation costs of regular grids scheme.Irregular fine or coarse grids in local fields not only cost less computing time than the conventional velocity-stress FDM,but also give a more accu- rate wavefield description.A dispersion analysis of the irregular-grid finite difference operator has confirmed the stability and high efficiency.The absorbing boundary condition is used to elimi- nate artificial reflections.Numerical examples show that this new irregular-grid finite difference method is of higher performance than conventional methods using regular rectangular grids in simulating elastic wave propagation in heterogeneous anisotropic porous media.展开更多
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.展开更多
基金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.
文摘For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.
基金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.
基金Research of Linjin Li and Jingmei Qiu is supported by the NSF grant NSF-DMS-1818924the Air Force Office of Scientific Computing FA9550-18-1-0257 and the University of Delawarethe Italian Ministry of Instruction,University and Research(MIUR)to support this research with funds coming from the PRIN Project 2017,No.2017KKJP4X and ITN-ETN Horizon 2020 Project,Project Reference 642768.
文摘In this paper,we propose a new conservative high-order semi-Lagrangian finite difference(SLFD)method to solve linear advection equation and the nonlinear Vlasov and BGK models.The finite difference scheme has better computational flexibility by working with point values,especially when working with high-dimensional problems in an operator splitting setting.The reconstruction procedure in the proposed SLFD scheme is motivated from the SL finite volume scheme.In particular,we define a new sliding average function,whose cell averages agree with point values of the underlying function.By developing the SL finite volume scheme for the sliding average function,we derive the proposed SLFD scheme,which is high-order accurate,mass conservative and unconditionally stable for linear problems.The performance of the scheme is showcased by linear transport applications,as well as the nonlinear Vlasov-Poisson and BGK models.Furthermore,we apply the Fourier stability analysis to a fully discrete SLFD scheme coupled with diagonally implicit Runge-Kutta(DIRK)method when applied to a stiff two-velocity hyperbolic relaxation system.Numerical stability and asymptotic accuracy properties of DIRK methods are discussed in theoretical and computational aspects.
文摘In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.
文摘We present a first-order finite difference scheme for approximating solutions of a mathematical model of cervical cancer induced by the human papillomavirus (HPV), which consists of four nonlinear partial differential equations and a nonlinear first-order ordinary differential equation. The scheme is analyzed and used to provide an existence-uniqueness result. Numerical simulations are performed in order to demonstrate the first-order rate of convergence. A sensitivity analysis was done in order to compare the effects of two drug types, those that increase the death rate of HPV-infected cells, and those that increase the death rate of the precancerous cell population. The model predicts that treatments that affect the precancerous cell population by directly increasing the corresponding death rate are far more effective than those that increase the death rate of HPV-infected cells.
基金This research was supported by the National Nature Science Foundation of China (No. 41074100) and the Program for NewCentury Excellent Talents in the University of the Ministry of Education of China (No. NCET- 10-0812).
文摘Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference (FD) method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.
基金the NSFC grant 11872210 and the Science Challenge Project,No.TZ2016002the NSFC Grant 11926103 when he visited Tianyuan Mathematical Center in Southeast China,Xiamen 361005,Fujian,Chinathe NSFC Grant 12071392 and the Science Challenge Project,No.TZ2016002.
文摘In this paper,a new type of finite difference mapped weighted essentially non-oscillatory(MWENO)schemes with unequal-sized stencils,such as the seventh-order and ninthorder versions,is constructed for solving hyperbolic conservation laws.For the purpose of designing increasingly high-order finite difference WENO schemes,the equal-sized stencils are becoming more and more wider.The more we use wider candidate stencils,the bigger the probability of discontinuities lies in all stencils.Therefore,one innovation of these new WENO schemes is to introduce a new splitting stencil methodology to divide some fourpoint or five-point stencils into several smaller three-point stencils.By the usage of this new methodology in high-order spatial reconstruction procedure,we get different degree polynomials defined on these unequal-sized stencils,and calculate the linear weights,smoothness indicators,and nonlinear weights as specified in Jiang and Shu(J.Comput.Phys.126:202228,1996).Since the difference between the nonlinear weights and the linear weights is too big to keep the optimal order of accuracy in smooth regions,another crucial innovation is to present the new mapping functions which are used to obtain the mapped nonlinear weights and decrease the difference quantity between the mapped nonlinear weights and the linear weights,so as to keep the optimal order of accuracy in smooth regions.These new MWENO schemes can also be applied to compute some extreme examples,such as the double rarefaction wave problem,the Sedov blast wave problem,and the Leblanc problem with a normal CFL number.Extensive numerical results are provided to illustrate the good performance of the new finite difference MWENO schemes.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foun- dation of China National Petroleum Corporation.
文摘To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.
文摘In the present paper,the numerical solution of It?type stochastic parabolic equation with a timewhite noise process is imparted based on a stochastic finite difference scheme.At the beginning,an implicit stochastic finite difference scheme is presented for this equation.Some mathematical analyses of the scheme are then discussed.Lastly,to ascertain the efficacy and accuracy of the suggested technique,the numerical results are discussed and compared with the exact solution.
基金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 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.
基金supported by China Scholarship Council and partially by the National "863" Program of China under contract No. 2007AA06Z218.
文摘Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.
基金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.
基金the National Key R&D Program of China(No.2018YFE0180900).
文摘The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.
基金The project supported by the National Natural Science Foundation of China(60073044)the State Key Development Programme for Basic Research of China(G1990022207).
文摘In this paper, we propose a novel incompressible finite-difference lattice Boltzmann Equation (FDLBE). Because source terms that reflect the interaction between phases can be accurately described, the new model is suitable for simulating two-way coupling incompressible multiphase flow The 2-D particle-laden flow over a backward-facing step is chosen as a test case to validate the present method. Favorable results are obtained and the present scheme is shown to have good prospects in practical applications.
基金supported and sponsored jointly by the National Natural Science Foundation of China(Grand Nos.51009092 and 50909061)Doctoral Foundation of the Ministry of Education of China(Grand No.20090073120013)the National High Technology Research and Development Program of China(863Program,Grand No.2008AA092301-1)
文摘The dynamic calculations of slender marine risers, such as Finite Element Method (FEM) or Modal Expansion Solution Method (MESM), are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation (FDA) method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bemoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser's free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.
文摘Basedonthe first-order Biot-equation with simplified coefficients,astaggered irregu- lar-grid finite difference method(FDM)is developed to simulate elastic wave propagation in 3-D heterogeneous anisotropic porous media.The ‘slow’P wave in porous media wave simulation is highly dispersive.Finer grids are needed to get a precise wavefield calculation for models with curved interface and complex geometric structure.Fine grids in a global model greatly increase computation costs of regular grids scheme.Irregular fine or coarse grids in local fields not only cost less computing time than the conventional velocity-stress FDM,but also give a more accu- rate wavefield description.A dispersion analysis of the irregular-grid finite difference operator has confirmed the stability and high efficiency.The absorbing boundary condition is used to elimi- nate artificial reflections.Numerical examples show that this new irregular-grid finite difference method is of higher performance than conventional methods using regular rectangular grids in simulating elastic wave propagation in heterogeneous anisotropic porous media.
基金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.