The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted...The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.展开更多
An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
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.展开更多
In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite differ...In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.展开更多
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.展开更多
The numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional （2D） tank are studied. The irregular tank is mapped onto a fixed square domain th...The numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional （2D） tank are studied. The irregular tank is mapped onto a fixed square domain through proper mapping functions. A staggered mesh system is employed in a 2D tank to calculate the elevation of the transient fluid. A time-independent finite difference method, which is developed by Bang- fuh Chen, is used to solve the Euler equations for incompressible and inviscid fluids. The numerical results agree well with the analytic solutions and previously published results. The sloshing profiles of surge and heave motion with initial standing waves are presented. The results show very clear nonlinear and beating phenomena.展开更多
A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing e...A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer （HLL） Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws （MUSCL）. The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.展开更多
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.展开更多
We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is ...We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.展开更多
A class of finite difference Hermite radial basis functions weighted essentially non-oscillatory(HRWENO)methods for solving conservation laws was presented by Abedian(Int.J.Numer.Meth.Fluids,94(2022),pp.583-607).To re...A class of finite difference Hermite radial basis functions weighted essentially non-oscillatory(HRWENO)methods for solving conservation laws was presented by Abedian(Int.J.Numer.Meth.Fluids,94(2022),pp.583-607).To reconstruct the fluxes in HRWENO,the common practice of reconstructing the flux functions was employed.In this follow-up research work,an alternative formulation to reconstruct the numerical fluxes is considered.First,the solution and its derivatives are directly employed to interpolate point values at interfaces of computational cells.Afterwards,the point values at interface of cell in building block are considered to obtain numerical fluxes.In this framework,arbitrary monotone fluxes can be employed,while in HRWENO the classical practice of reconstructing flux functions can be considered only to smooth flux splitting.Also,in the process of reconstruction these type of schemes consider the effectively narrower stencil of HRWENO methods.Extensive test cases such as Euler equations of compressible gas dynamics are considered to show the good performance of the methods.展开更多
To solve the equation for gravity-gyroscopic waves in a rectangular domain, the distinguished algorithm for the solution of the Cauchy problem for a second-order transient equation is proposed. This algorithm is devel...To solve the equation for gravity-gyroscopic waves in a rectangular domain, the distinguished algorithm for the solution of the Cauchy problem for a second-order transient equation is proposed. This algorithm is developed by using the time-varying finite element method. The space derivatives in the gravity-gyroscopic wave equation are approximated with finite differences. The stability and accuracy of the method are estimated. The procedure for the implementation of the method is developed. The calculations were performed for determining the steady-state modes of fluctuations of the solutions of the gravity-gyroscopic wave equation depending on the problem parameters.展开更多
This paper deals with numerical methods for solving one-dimensional(1D)and twodimensional(2D)initial-boundary value problems(IBVPs)of space-fractional sine-Gordon equations(SGEs)with distributed delay.For 1D problems,...This paper deals with numerical methods for solving one-dimensional(1D)and twodimensional(2D)initial-boundary value problems(IBVPs)of space-fractional sine-Gordon equations(SGEs)with distributed delay.For 1D problems,we construct a kind of oneparameter finite difference(OPFD)method.It is shown that,under a suitable condition,the proposed method is convergent with second order accuracy both in time and space.In implementation,the preconditioned conjugate gradient(PCG)method with the Strang circulant preconditioner is carried out to improve the computational efficiency of the OPFD method.For 2D problems,we develop another kind of OPFD method.For such a method,two classes of accelerated schemes are suggested,one is alternative direction implicit(ADI)scheme and the other is ADI-PCG scheme.In particular,we prove that ADI scheme can arrive at second-order accuracy in time and space.With some numerical experiments,the computational effectiveness and accuracy of the methods are further verified.Moreover,for the suggested methods,a numerical comparison in computational efficiency is presented.展开更多
This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either...This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.展开更多
We propose an efficient numerical method for two population models, based on the nonstandard finite difference （NSFD） schemes and composition methods with complex time steps. The NSFD scheme is able to give positive...We propose an efficient numerical method for two population models, based on the nonstandard finite difference （NSFD） schemes and composition methods with complex time steps. The NSFD scheme is able to give positive numerical solutions that satisfy the conservation law, which is a key property for biological population models. The accuracy is improved by using the composition methods with complex time steps. Numerical tests on the plankton nutrient model and whooping cough model are presented to show the efficiency and advantage of the proposed numerical method.展开更多
In this paper, we consider the existence, uniqueness and convergence of weak and strongimplicit difference solution for the first boundary problem of quasilinear parabolic system: u_t=(-1)^(M+1)A(x,t,u,…,u_x^(M-1))u_...In this paper, we consider the existence, uniqueness and convergence of weak and strongimplicit difference solution for the first boundary problem of quasilinear parabolic system: u_t=(-1)^(M+1)A(x,t,u,…,u_x^(M-1))u_x^(2M)+f(x,t,u,…,u_x^(2M-1)), (x,t)∈Q_T={0<x<l, 0<t≤T}, (1) u_xk(0,t)=u_xk(l,t)=0,(k=0,1,…,M-1), 0<t≤T, (2) u(x,0) =ψ(x), 0≤x≤l, (3)where u, ψ and f are m-dimensional vector valued functions, A is an m×m positivelydefinite matrix and u_xk denotes ?~ku/?_xk. For this problem, the estimations of the differencesolution are obtained. As h→0, △t→0, the difference solution converges weakly in W_2^(2M,1) (QT)to the unique generalized solution u(x,t)∈W_2^(2M,1)(QT) of problems (1), (2), (3). Especially,a favorable restriction condition to the step lengths △t and h for explicit and weak implicitschemes is found.展开更多
A finite difference scheme for the generalized nonlinear Schr dinger equation with variable coefficients is developed. The scheme is shown to satisfy two conservation laws. Numerical results show that the scheme is a...A finite difference scheme for the generalized nonlinear Schr dinger equation with variable coefficients is developed. The scheme is shown to satisfy two conservation laws. Numerical results show that the scheme is accurate and efficient.展开更多
In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,whi...In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,which contain nonlinear high order derivatives,and possibly peakon solutions or shock waves.By introducing auxiliary variable(s),we rewrite the DP equation as a hyperbolic-elliptic system,and theµDP equation as afirst order system.Then we choose a linearfinite difference scheme with suitable order of accuracy for the auxiliary variable(s),and twofinite difference WENO schemes with unequal-sized sub-stencils for the primal variable.One WENO scheme uses one large stencil and several smaller stencils,and the other WENO scheme is based on the multi-resolution framework which uses a se-ries of unequal-sized hierarchical central stencils.Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil,both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights.Another advantage is that thefinal reconstructed polynomial on the target cell is a polynomial of the same de-gree as the polynomial over the big stencil,while the classicalfinite difference WENO reconstruction can only be obtained for specific points inside the target interval.Nu-merical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.展开更多
In this paper,we propose a new conservative semi-Lagrangian(SL)finite difference(FD)WENO scheme for linear advection equations,which can serve as a base scheme for the Vlasov equation by Strang splitting[4].The recons...In this paper,we propose a new conservative semi-Lagrangian(SL)finite difference(FD)WENO scheme for linear advection equations,which can serve as a base scheme for the Vlasov equation by Strang splitting[4].The reconstruction procedure in the proposed SL FD scheme is the same as the one used in the SL finite volume(FV)WENO scheme[3].However,instead of inputting cell averages and approximate the integral form of the equation in a FV scheme,we input point values and approximate the differential form of equation in a FD spirit,yet retaining very high order(fifth order in our experiment)spatial accuracy.The advantage of using point values,rather than cell averages,is to avoid the second order spatial error,due to the shearing in velocity(v)and electrical field(E)over a cell when performing the Strang splitting to the Vlasov equation.As a result,the proposed scheme has very high spatial accuracy,compared with second order spatial accuracy for Strang split SL FV scheme for solving the Vlasov-Poisson(VP)system.We perform numerical experiments on linear advection,rigid body rotation problem;and on the Landau damping and two-stream instabilities by solving the VP system.For comparison,we also apply(1)the conservative SL FD WENO scheme,proposed in[22]for incompressible advection problem,(2)the conservative SL FD WENO scheme proposed in[21]and(3)the non-conservative version of the SL FD WENO scheme in[3]to the same test problems.The performances of different schemes are compared by the error table,solution resolution of sharp interface,and by tracking the conservation of physical norms,energies and entropies,which should be physically preserved.展开更多
A compressible lattice Boltzmann-finite difference method is extended by the phase-field approach into a monolithic scheme to study fluid flow and heat transfer through regular arrangements of solid bodies of circular...A compressible lattice Boltzmann-finite difference method is extended by the phase-field approach into a monolithic scheme to study fluid flow and heat transfer through regular arrangements of solid bodies of circular,elliptical and irregular shapes.The advantage of using the phase-field method is demon-strated both in its simplicity of accounting for flow and thermal boundary conditions at solid surfaces with irregular shapes and in the capability of generating such complex-shaped objects.For an array of discs,numerical results for the overall solid-to-gas heat transfer rate are validated via experiments on flow through arrays of hot cylinders.The thus validated compressible LB-FD-PF hybrid scheme is used to study the dependence of heat transfer on flow and thermal boundary conditions(Reynolds number,temperature difference between the hot solid bodies and the inlet gas),porosity as well as on the shape of solid objects.Results are rationalized in terms of the residence time of the gas close to the solid body and downstream variations of gas velocity and temperature.Perspective for further applications of the proposed methodology are also discussed.展开更多
文摘The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
文摘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.
文摘In the article, the fully discrete finite difference scheme for a type of nonlinear reaction-diffusion equation is established. Then the new function space is introduced and the stability problem for the finite difference scheme is discussed by means of variational approximation method in this function space. The approach used is of a simple characteristic in gaining the stability condition of the scheme.
基金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 numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional （2D） tank are studied. The irregular tank is mapped onto a fixed square domain through proper mapping functions. A staggered mesh system is employed in a 2D tank to calculate the elevation of the transient fluid. A time-independent finite difference method, which is developed by Bang- fuh Chen, is used to solve the Euler equations for incompressible and inviscid fluids. The numerical results agree well with the analytic solutions and previously published results. The sloshing profiles of surge and heave motion with initial standing waves are presented. The results show very clear nonlinear and beating phenomena.
基金supported by the National Natural Science Foundation of China(Grant No.51579034)the Open Fund of the Key Laboratory of Ocean Circulation and Waves,Chinese Academy of Sciences(Grant No.KLOCW1502)
文摘A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer （HLL） Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws （MUSCL）. The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.
基金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.
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.
文摘A class of finite difference Hermite radial basis functions weighted essentially non-oscillatory(HRWENO)methods for solving conservation laws was presented by Abedian(Int.J.Numer.Meth.Fluids,94(2022),pp.583-607).To reconstruct the fluxes in HRWENO,the common practice of reconstructing the flux functions was employed.In this follow-up research work,an alternative formulation to reconstruct the numerical fluxes is considered.First,the solution and its derivatives are directly employed to interpolate point values at interfaces of computational cells.Afterwards,the point values at interface of cell in building block are considered to obtain numerical fluxes.In this framework,arbitrary monotone fluxes can be employed,while in HRWENO the classical practice of reconstructing flux functions can be considered only to smooth flux splitting.Also,in the process of reconstruction these type of schemes consider the effectively narrower stencil of HRWENO methods.Extensive test cases such as Euler equations of compressible gas dynamics are considered to show the good performance of the methods.
文摘To solve the equation for gravity-gyroscopic waves in a rectangular domain, the distinguished algorithm for the solution of the Cauchy problem for a second-order transient equation is proposed. This algorithm is developed by using the time-varying finite element method. The space derivatives in the gravity-gyroscopic wave equation are approximated with finite differences. The stability and accuracy of the method are estimated. The procedure for the implementation of the method is developed. The calculations were performed for determining the steady-state modes of fluctuations of the solutions of the gravity-gyroscopic wave equation depending on the problem parameters.
基金supported by the NSFC(Grant No.11971010)the Science and Technology Development Fund of Macao(Grant No.0122/2020/A3)MYRG2020-00224-FST from University of Macao,China.
文摘This paper deals with numerical methods for solving one-dimensional(1D)and twodimensional(2D)initial-boundary value problems(IBVPs)of space-fractional sine-Gordon equations(SGEs)with distributed delay.For 1D problems,we construct a kind of oneparameter finite difference(OPFD)method.It is shown that,under a suitable condition,the proposed method is convergent with second order accuracy both in time and space.In implementation,the preconditioned conjugate gradient(PCG)method with the Strang circulant preconditioner is carried out to improve the computational efficiency of the OPFD method.For 2D problems,we develop another kind of OPFD method.For such a method,two classes of accelerated schemes are suggested,one is alternative direction implicit(ADI)scheme and the other is ADI-PCG scheme.In particular,we prove that ADI scheme can arrive at second-order accuracy in time and space.With some numerical experiments,the computational effectiveness and accuracy of the methods are further verified.Moreover,for the suggested methods,a numerical comparison in computational efficiency is presented.
基金supported by the NSF under Grant DMS-2208391sponsored by the NSF under Grant DMS-1753581.
文摘This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.
文摘We propose an efficient numerical method for two population models, based on the nonstandard finite difference （NSFD） schemes and composition methods with complex time steps. The NSFD scheme is able to give positive numerical solutions that satisfy the conservation law, which is a key property for biological population models. The accuracy is improved by using the composition methods with complex time steps. Numerical tests on the plankton nutrient model and whooping cough model are presented to show the efficiency and advantage of the proposed numerical method.
文摘In this paper, we consider the existence, uniqueness and convergence of weak and strongimplicit difference solution for the first boundary problem of quasilinear parabolic system: u_t=(-1)^(M+1)A(x,t,u,…,u_x^(M-1))u_x^(2M)+f(x,t,u,…,u_x^(2M-1)), (x,t)∈Q_T={0<x<l, 0<t≤T}, (1) u_xk(0,t)=u_xk(l,t)=0,(k=0,1,…,M-1), 0<t≤T, (2) u(x,0) =ψ(x), 0≤x≤l, (3)where u, ψ and f are m-dimensional vector valued functions, A is an m×m positivelydefinite matrix and u_xk denotes ?~ku/?_xk. For this problem, the estimations of the differencesolution are obtained. As h→0, △t→0, the difference solution converges weakly in W_2^(2M,1) (QT)to the unique generalized solution u(x,t)∈W_2^(2M,1)(QT) of problems (1), (2), (3). Especially,a favorable restriction condition to the step lengths △t and h for explicit and weak implicitschemes is found.
文摘A finite difference scheme for the generalized nonlinear Schr dinger equation with variable coefficients is developed. The scheme is shown to satisfy two conservation laws. Numerical results show that the scheme is accurate and efficient.
基金supported by National Natural Science Foundation of China(Grant No.12071455)supported by National Natural Science Foundation of China(Grant No.11871428)。
文摘In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,which contain nonlinear high order derivatives,and possibly peakon solutions or shock waves.By introducing auxiliary variable(s),we rewrite the DP equation as a hyperbolic-elliptic system,and theµDP equation as afirst order system.Then we choose a linearfinite difference scheme with suitable order of accuracy for the auxiliary variable(s),and twofinite difference WENO schemes with unequal-sized sub-stencils for the primal variable.One WENO scheme uses one large stencil and several smaller stencils,and the other WENO scheme is based on the multi-resolution framework which uses a se-ries of unequal-sized hierarchical central stencils.Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil,both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights.Another advantage is that thefinal reconstructed polynomial on the target cell is a polynomial of the same de-gree as the polynomial over the big stencil,while the classicalfinite difference WENO reconstruction can only be obtained for specific points inside the target interval.Nu-merical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.
基金supported by AFOSR grant FA9550-09-1-0344 and NSF grant DMS-0914852supported by AFOSR grant FA9550-09-1-0126 and NSF grant DMS-0809086.
文摘In this paper,we propose a new conservative semi-Lagrangian(SL)finite difference(FD)WENO scheme for linear advection equations,which can serve as a base scheme for the Vlasov equation by Strang splitting[4].The reconstruction procedure in the proposed SL FD scheme is the same as the one used in the SL finite volume(FV)WENO scheme[3].However,instead of inputting cell averages and approximate the integral form of the equation in a FV scheme,we input point values and approximate the differential form of equation in a FD spirit,yet retaining very high order(fifth order in our experiment)spatial accuracy.The advantage of using point values,rather than cell averages,is to avoid the second order spatial error,due to the shearing in velocity(v)and electrical field(E)over a cell when performing the Strang splitting to the Vlasov equation.As a result,the proposed scheme has very high spatial accuracy,compared with second order spatial accuracy for Strang split SL FV scheme for solving the Vlasov-Poisson(VP)system.We perform numerical experiments on linear advection,rigid body rotation problem;and on the Landau damping and two-stream instabilities by solving the VP system.For comparison,we also apply(1)the conservative SL FD WENO scheme,proposed in[22]for incompressible advection problem,(2)the conservative SL FD WENO scheme proposed in[21]and(3)the non-conservative version of the SL FD WENO scheme in[3]to the same test problems.The performances of different schemes are compared by the error table,solution resolution of sharp interface,and by tracking the conservation of physical norms,energies and entropies,which should be physically preserved.
基金funded by the Deutsche For-schungsgemeinschaft(DFG,German Research Foundation)-422037413-CRC/TRR 287"BULK-REACTION".
文摘A compressible lattice Boltzmann-finite difference method is extended by the phase-field approach into a monolithic scheme to study fluid flow and heat transfer through regular arrangements of solid bodies of circular,elliptical and irregular shapes.The advantage of using the phase-field method is demon-strated both in its simplicity of accounting for flow and thermal boundary conditions at solid surfaces with irregular shapes and in the capability of generating such complex-shaped objects.For an array of discs,numerical results for the overall solid-to-gas heat transfer rate are validated via experiments on flow through arrays of hot cylinders.The thus validated compressible LB-FD-PF hybrid scheme is used to study the dependence of heat transfer on flow and thermal boundary conditions(Reynolds number,temperature difference between the hot solid bodies and the inlet gas),porosity as well as on the shape of solid objects.Results are rationalized in terms of the residence time of the gas close to the solid body and downstream variations of gas velocity and temperature.Perspective for further applications of the proposed methodology are also discussed.
