Efficient Finite Difference/Spectral Method for the Time Fractional Ito Equation Using Fast Fourier Transform Technic

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.

High-Order Bound-Preserving Finite Difference Methods for Multispecies and Multireaction Detonations

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.

Numerical Study by Imposing the Finite Difference Method for Unsteady Casson Fluid Flow with Heat Flux

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 Application of the Generalized Finite Difference Method (GFDM) for Modelling Geophysical Test

The possibility of using a nodal method allowing irregular distribution of nodes in a natural way is one of the main advantages of the generalized finite difference method (GFDM) with regard to the classical finite difference method. Moreover, this feature has made it one of the most-promising meshless methods because it also allows us to reduce the time-consuming task of mesh generation and the numerical solution of integrals. This characteristic allows us to shape geological features easily whilst maintaining accuracy in the results, which can be a source of great interest when dealing with this kind of problems. Two widespread geophysical investigation methods in civil engineering are the cross-hole method and the seismic refraction method. This paper shows the use of the GFDM to model the aforementioned geophysical investigation tests showing precision in the obtained results when comparing them with experimental data.

An explicit finite element－finite difference method for analyzing the effect of visco－elastic local topography on the earthquake motion

An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.

COMPACT FINITE DIFFERENCE-FOURIER SPECTRAL METHOD FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

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.

A method of solving the stiffness problem in Biot＇s poroelastic equations using a staggered high-order finite-difference

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.

Implicit finite difference method for fractional percolation equation with Dirichlet and fractional boundary conditions

An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.

Numerical storm surge model with higher order finite difference method of lines for the coast of Bangladesh

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.

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

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.

Using finite difference method to simulate casting thermal stress

Thermal stress simulation can provide a scientific reference to eliminate defects such as crack,residual stress centralization and deformation etc.,caused by thermal stress during casting solidification.To study the thermal stress distribution during casting process,a unilateral thermal-stress coupling model was employed to simulate 3D casting stress using Finite Difference Method(FDM),namely all the traditional thermal-elastic-plastic equations are numerically and differentially discrete.A FDM/FDM numerical simulation system was developed to analyze temperature and stress fields during casting solidification process.Two practical verifications were carried out,and the results from simulation basically coincided with practical cases.The results indicated that the FDM/FDM stress simulation system can be used to simulate the formation of residual stress,and to predict the occurrence of hot tearing.Because heat transfer and stress analysis are all based on FDM,they can use the same FD model,which can avoid the matching process between different models,and hence reduce temperature-load transferring errors.This approach makes the simulation of fluid flow,heat transfer and stress analysis unify into one single model.

A Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations

In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Differential Equations (PDE) with mixed derivatives along with Dirichlet and Neumann boundary conditions. These two methods have almost the same accuracy from theoretical aspect with regular boundaries, but generally Finite Element Method produces better approximations when the boundaries are irregular. In order to investigate which method produces better results from numerical aspect, we apply these methods into specific examples with regular boundaries with constant step-size for both of them. The results which obtained confirm, in most of the cases, the theoretical results.

A Finite Difference Method for Determining Interdiffusivity of Aluminide Coating Formed on Superalloy

A numerical method has been developed to extract the composition-dependent interdiffusivity from the concentration profiles in the aluminide coating prepared by pack cementation. The procedure is based on the classic finite difference method (FDM). In order to simplify the model, effect of some alloying elements on interdiffusivity can be negligible. Calculated results indicate the interdiffusivity in aluminide coating strongly depends on the composition and give the formulas used to calculate interdiffusivity at 850, 950 and 1050癈. The effect on interdiffusivity is briefly discussed.

Full-vectorial finite-difference beam propagation method based on the modified alternating direction implicit method

A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.

A geometrically and locally adaptive remeshing method for finite difference modeling of mining-induced surface subsidence

Surface subsidence induced by underground mining is a typical serious geohazard.Numerical approaches such as the discrete element method(DEM)and finite difference method(FDM)have been widely used to model and analyze mining-induced surface subsidence.However,the DEM is typically computationally expensive,and is not capable of analyzing large-scale problems,while the mesh distortion may occur in the FDM modeling of largely deformed surface subsidence.To address the above problems,this paper presents a geometrically and locally adaptive remeshing method for the FDM modeling of largely deformed surface subsidence induced by underground mining.The essential ideas behind the proposed method are as follows:(i)Geometrical features of elements(i.e.the mesh quality),rather than the calculation errors,are employed as the indicator for determining whether to conduct the remeshing;and(ii)Distorted meshes with multiple attributes,rather than those with only a single attribute,are locally regenerated.In the proposed method,the distorted meshes are first adaptively determined based on the mesh quality,and then removed from the original mesh model.The tetrahedral mesh in the distorted area is first regenerated,and then the physical field variables of old mesh are transferred to the new mesh.The numerical calculation process recovers when finishing the regeneration and transformation.To verify the effectiveness of the proposed method,the surface deformation of the Yanqianshan iron mine,Liaoning Province,China,is numerically investigated by utilizing the proposed method,and compared with the numerical results of the DEM modeling.Moreover,the proposed method is applied to predicting the surface subsidence in Anjialing No.1 Underground Mine,Shanxi Province,China.

Solution precision of concrete temperature fields with finite difference method

With the finite difference method to calculate the temperature distribution in mass concrete structures, the solution precision will increase with a smaller step size, at the cost of computational time. In view of the basic characteristics of the finite difference method, a simple yet powerful improvement is introduced. By multiplying the adiabatic temperature function with a correction factor, the precision of the solution can be assured without an increase in the computation time. In addition, the correction rules for three types of commonly used concrete hydration formulas are investigated.

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.

Finite difference streamline diffusion method using nonconforming space for incompressible time-dependent Navier-Stokes equations

This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart （CR） element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.

Comparison of finite difference and pseudo-spectral methods in forward modelling based on metal ore model of random media

With more applications of seismic exploration in metal ore exploration,forward modelling of seismic wave has become more important in metal ore. Finite difference method and pseudo-spectral method are two important methods of wave-field simulation. Results of previous studies show that both methods have distinct advantages and disadvantages: Finite difference method has high precision but its dispersion is serious; pseudospectral method considers both computational efficiency and precision but has less precision than finite-difference. The authors consider the complex structural characteristics of the metal ore,furthermore add random media in order to simulate the complex effects produced by metal ore for wave field. First,the study introduced the theories of random media and two forward modelling methods. Second,it compared the simulation results of two methods on fault model. Then the authors established a complex metal ore model,added random media and compared computational efficiency and precision. As a result,it is found that finite difference method is better than pseudo-spectral method in precision and boundary treatment,but the computational efficiency of pseudospectral method is slightly higher than the finite difference method.

Finite difference method for dynamic response analysis of anchorage system

Based on some assumptions, the dynamic analysis model of anchorage system is established. The dynamic governing equation is expressed as finite difference format and programmed by using MATLAB language. Compared with theoretical method, the finite difference method has been verified to be feasible by a case study. It is found that under seismic loading, the dynamic response of anchorage system is synchronously fluctuated with the seismic vibration. The change of displacement amplitude of material points is slight, and comparatively speaking, the displacement amplitude of the outside point is a little larger than that of the inside point, which shows amplification effect of surface. While the axial force amplitude transforms considerably from the inside to the outside. It increases first and reaches the peak value in the intersection between the anchoring section and free section, then decreases slowly in the free section. When considering damping effect of anchorage system, the finite difference method can reflect the time attenuation characteristic better, and the calculating result would be safer and more reasonable than the dynamic steady-state theoretical method. What is more, the finite difference method can be applied to the dynamic response analysis of harmonic and seismic random vibration for all kinds of anchor, and hence has a broad application prospect.