基于DEM-FDM耦合的过渡段膨胀诱发钢轨上拱研究

研究目的:为分析涵洞过渡段地基膨胀引起的钢轨上拱响应,基于现场测试、室内膨胀试验数据,开展DEM-FDM耦合数值模拟,分析某涵洞附近路基土在膨胀范围为16 m,膨胀中心距离涵洞中心分别为0 m、5 m、10 m这三种工况下,不同膨胀率时基床填料的运动规律及钢轨的上拱响应。研究结论:(1)涵洞对于钢轨上拱位移的传递存在阻断作用,但会增大钢轨上拱的峰值,原位膨胀率下工况二的钢轨上拱峰值达到46 mm,当路基膨胀率为0.3%时钢轨上拱位移量达到无砟轨道钢轨可调节临界值(4mm);(2)过渡段钢轨上拱处同时产生轴向应力集中,其中原位膨胀率下工况二轴向应力峰值达到14.4 MPa;(3)对于膨胀区域位于涵洞下方的工况,钢轨轴向应力呈现出来的分布规律为钢轨上拱拱顶处为主拉应力状态,拱脚处为主压应力状态,因此一共包括三个压应力峰值点以及两个拉应力峰值点;(4)本文研究可为高铁涵洞过渡段路基膨胀病害解决方案的确定提供理论依据。

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.

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

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.

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.

高阶SF-SFDTD方法在含时薛定谔方程求解中的应用研究

时域有限差分方法(finite-difference time-domain,FDTD(2,2))被广泛用于量子力学中薛定谔方程的求解,然而受Courant-Friedrichs-Lewy(CFL)条件的影响,计算空间中的网格尺寸会限制时间步长的取值范围,极大降低了FDTD(2,2)方法的数值计算效率.另外,FDTD(2,2)方法在时间域和空间域只具有二阶数值精度,在计算中往往会导致较大的误差累计,影响仿真结果的正确性.为了克服这些问题,结合空间滤波方法(spatial filtering,SF)和高阶辛时域有限差分(symplectic finite-difference time-domain,SFDTD(3,4))方法(3和4分别表示时间和空间数值精度),提出了一种时间稳定性条件可扩展的SF-SFDTD(3,4)方法用于求解含时薛定谔方程.SF-SFDTD(3,4)方法无需对传统SFDTD(3,4)方法的迭代公式进行进一步的推导,只需要在每一次的数值迭代过程中加入空间滤波操作,滤除因采用不满足CFL条件的时间步长而产生的不稳定空间域高频分量,保证数值方法的稳定性,因此所提方法与传统SFDTD(3,4)方法具有较高的兼容性.同时,理论分析了SF-SFDTD(3,4)方法的数值色散误差.最后,通过数值算例验证了本文所提方法的正确性和有效性.

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.

基于ADI-FDTD法的探地雷达正演数值模拟及验证

探地雷达(GPR)是一种应用前景广泛、用于探测和定位地下物体的浅层地球物理方法。通过开展探地雷达正演模拟研究,可以获得复杂地下结构的探地雷达图像回波特征。时域有限差分(FDTD)法受到稳定性和收敛性条件的限制,导致效率和精度低。交替方向隐式-时域有限差分(ADI-FDTD)法克服FDTD法的稳定性限制,可以选择更大的时间步长来提高计算效率。从原理出发,对ADI-FDTD法进行公式推导,在基于ADI-FDTD法的基础上进行卷积完全匹配层(CPML)的结合。通过进行ADI-FDTD法的时间稳定性仿真和不同材质、不同填充介质的双管,沙槽正演模拟实验,结果表明,ADI-FDTD法可以保证时间无条件稳定性,相同情况下可采用更大的时间步长进行正演模拟,提高正演效率。对不同管线、不同材质情况下的正演剖面曲线特征进行解译分析,最后与沙槽实测进行对比,证明得到的探测解译结果与实际状况达到较好的吻合。

CFD-Based Method of Determining Form Factor k for Different Ship Types and Different Drafts

The value of form factor k at different drafts is important in predicting full-scale total resistance and speed for different types of ships. In the ITTC community, most organizations predict form factor k using a low-speed model test. However, this method is problematic for ships with bulbous bows and transom. In this article, a Computational Fluid Dynamics(CFD)-based method is introduced to obtain k for different type of ships at different drafts, and a comparison is made between the CFD method and the model test. The results show that the CFD method produces reasonable k values. A grid generating method and turbulence model are briefly discussed in the context of obtaining a consistent k using CFD.

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.

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.

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.

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.

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.

SUPERCONVERGENCE OF GENERALIZED DIFFERENCE METHOD FOR ELLIPTIC BOUNDARY VALUE PROBLEM

Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that for finite element method.

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.

A Comparative Study of the Difference Method and the Enzyme-Hydrolyzed Casein Method for Determining True Amino Acid Digestibilities and Endogenous Amino Acid Losses in Duck Feed

In this study, we examined the varia- tions between the difference method and the enzyme- hydrolyzed casein method for determining endogenous amino acid loss and the true amino acid digestibility in ducks fed normal protein-containing diets. These methods were compared to the nitrogen-free （N-free） diet method. The difference method was based on soy- bean meal as the only protein source, with the experi- mental diets containing crude protein levels at 15% and 20%. The enzyme-hydrolyzed casein method was based on enzyme-hydrolyzed casein meal as the pro- tein source, with the experimental diet containing a crude protein level of 17.5%. The N-free diet was prepared with starches and paper fibers. In each meth- od,64 Tianfu meat drakes （7-weeks-old） with an av- erage body weight of 2.77±0.16 kg were used and divided into four groups, and fed four different diets. Each group contained four replicates of four drakes and they were force fed trial diets according to the Sirbald method for detecting their apparent amino aciddigestibility, endogenous amino acid loss and true a- mino acid digestibility. The results demonstrated that using the difference, enzyme-hydrolyzed casein and N-free diet methods, endogenous amino acid losses were 0. 9946,1. 2243 and 0. 9297 mg/g dry matter in- take （ DMI）, respectively. The true amino acid digest- ibility measured by the difference method was 88.93 %±4.43 %. Using the enzyme-hydrolyzed ca- sein method with two dietary crude protein levels of 15% and 20%, the digestibility was 91.15%±4.33% and 91.97%±4. 16%, respectively, and by the N-free diet methods with two dietary crude protein levels of 15% and 20% ,it was 88.55%±4.29% and 88.82 %±4.61%, respectively. The results suggested that when the dietary protein level was 15% to 20 %, the true amino acid digestibility and endogenous ami- no acid loss as determined by the difference method was more accurate than the values determined by the enzyme-hydrolyzed casein method.

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 SIMPLE PROOF OF STABILITY AND CONVERGENCE IN L＿2 FOR SOME DIFFERENCE SCHEMES AND THEIR EXTRAPOLATION METHOD FOR PARABOLIC EQUATIONS

ＡＳＩＭＰＬＥＰＲＯＯＦＯＦＳＴＡＢＩＬＩＴＹＡＮＤＣＯＮＶＥＲＧＥＮＣＥＩＮＬ＿２ＦＯＲＳＯＭＥＤＩＦＦＥＲＥＮＣＥＳＣＨＥＭＥＳＡＮＤＴＨＥＩＲＥＸＴＲＡＰＯＬＡＴＩＯＮＭＥＴＨＯＤＦＯＲＰＡＲＡＢＯＬＩＣＥＱＵＡＴＩＯＮＳＳｕｎＺｈｉｚｈｏｎｇ...

Computing Bifurcation Diagrams of Steady State Kuramoto Sivashinsky Equation by Difference Method

Utilizing difference formulae, we obtained the discrete systems of steady state Kuramoto Sivashinsky (K S) equation. Applied Newton's method and continuation technology to the systems, the bifurcated solutions are derived, and the bifurcation diagrams are constructed. All the results are successful and satisfactory.