A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.

Techniques for improving computational speed in numerical simulation of casting thermal stress based on finite difference method

Finite difference method (FDM) was applied to simulate thermal stress recently, which normally needs a long computational time and big computer storage. This study presents two techniques for improving computational speed in numerical simulation of casting thermal stress based on FDM, one for handling of nonconstant material properties and the other for dealing with the various coefficients in discretization equations. The use of the two techniques has been discussed and an application in wave-guide casting is given. The results show that the computational speed is almost tripled and the computer storage needed is reduced nearly half compared with those of the original method without the new technologies. The stress results for the casting domain obtained by both methods that set the temperature steps to 0.1 ℃ and 10 ℃, respectively are nearly the same and in good agreement with actual casting situation. It can be concluded that both handling the material properties as an assumption of stepwise profile and eliminating the repeated calculation are reliable and effective to improve computational speed, and applicable in heat transfer and fluid flow simulation.

A simple method of depressing numerical dissipation effects during wave simulation within the Euler model

Numerical wave tanks are widely-acknowledged tools in studying waves and wave-structure interactions. They can generate waves under realistic scales and offers more information on the fluid field. However, most numerical wave tanks suffer from issues known as the numerical dissipation and numerical dispersion. The former causes wave energy to be slowly dissipated and the latter shifts wave frequencies during wave propagation. This paper proposes a simple method of depressing numerical dissipation effects on the basis of solving Euler equations using the finite difference method(FDM). The wave propagation solutions are solved analytically taking into account the influence of the damping terms. The main idea of the method is to append a source term to the momentum equation, whose strength is determined by how strong the numerical damping effect is. The method is verified by successfully depressing numerical effects during the simulation of regular linear waves, Stokes waves and irregular waves. By applying the method, wave energy is able to be close to its initial value after long distance of travel.

Numerical simulation of intermediate phase growth in Ti/Al alternate foils

To investigate the diffusion reaction between Ti/Al solid diffusion couple, Ti/Al alternate foils formed by hot pressing were annealed at 525, 550, 575 and 600 °C for time ranging from 1 to 40 h. The experimental results show that TiAl3 was the only observed phase at Ti/Al interface. The interface thermodynamics favored the preferential formation of TiAl3 in Ti/Al couple. The growth of TiAl3 layer occurred mainly towards Al foil side and exhibited a parabolic law. Using the interdiffusion coefficients calculated based on the contribution of grain boundary diffusion, the growth of TiAl3 was simulated numerically with the finite difference method, and the simulated results were in good agreement with the experimental ones.

Sloshing simulation of standing wave with time-independent finite difference method for Euler equations

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.

Numerical simulation of boundary heat flow effects on directional solidification microstructure of a binary alloy

The boundary heat flow has important significance for the microstructures of directional solidified binary alloy.Interface evolution of the directional solidified microstructure with different boundary heat flow was discussed.In this study, only one interface was allowed to have heat flow, and Neumann boundary conditions were imposed at the other three interfaces.From the calculated results, it was found that different boundary heat flows will result in different microstructures.When the boundary heat flow equals to 20 W·cm-2, the growth of longitudinal side branches is accelerated and the growth of transverse side branches is restrained, and meanwhile, there is dendritic remelting in the calculation domain.When the boundary heat flow equals to 40 W·cm-2, the growths of the transverse and longitudinal side branches compete with each other, and when the boundary heat flow equals to 100-200 W·cm-2, the growth of transverse side branches dominates absolutely.The temperature field of dendritic growth was analyzed and the relation between boundary heat flow and temperature field was also investigated.

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.

Numerical Simulation of 2D and 3D Sloshing Waves in a Regularly and Randomly Excited Container

In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.

Numerical Simulation of Groundwater Pollution Problems Based on Convection Diffusion Equation

The analytical solution of the convection diffusion equation is considered by two-dimensional Fourier transform and the inverse Fourier transform. To get the numerical solution, the Crank-Nicolson finite difference method is constructed, which is second-order accurate in time and space. Numerical simulation shows excellent agreement with the analytical solution. The dynamic visualization of the simulating results is realized on ArcGIS platform. This work provides a quick and intuitive decision-making basis for water resources protection, especially in dealing with water pollution emergencies.

A numerical simulation of surface wave excitation in a rectangular planar-type plasma source

The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.

Numerical Simulation of Microstructure Formation and Forecast of Mechanical Properties for Spheroidal Graphite Iron

The numerical simulation can overcome the hardship of mathematical analysis and experimental research, explicate the mechanism of microstructure shaping, predict mechanical properties and operating life of castings and then optimize technology and control microstructure formation to obtain the qualified castings. The finite difference method (FDM) is applied to the simulation of temperature field based on all kinds of nucleation and growth models on all stages of solidification of spheroidal graphite cast iron. Visual C++ is used to program the numerical simulation software, QTstructure-1 to simulate the solidification process of spheroidal graphite cast iron and the formation of all phases in solidification process. Finally, the result of simulation is well agreed with the experimental result.

Numerical simulation of tungsten alloy in powder injection molding process

The flow behavior of feedstock for the tungsten alloy powder in the mold cavity was approximately described using Hele-Shaw flow model. The math model consisting of momentum equation, consecutive equation and thermo-conduction equation for describing the injection process was established. The equations are solved by the finite element/finite difference hybrid method that means dispersing the feedstock model with finite element method, resolving the model along the depth with finite difference method, and tracking the movable boundary with control volume method, then the pressure equation and energy equation can be resolved in turn. The numerical simulation of the injection process and the identification of the process parameters were realized by the Moldflow software. The results indicate that there is low temperature gradient in the cavity while the pressure and shear rate gradient are high at high flow rate. The selection of the flow rate is affected by the structure of the gate. The shear rate and the pressure near the gate can be decreased by properly widening the dimension of the gate. There is a good agreement between the process parameters obtained by the numerical simulation and the actual ones.

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.

Failure mechanism of directional roof cutting and design method optimization

Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock mass,the implementation of this technology often encounters design challenges,leading to suboptimal results and increased costs.This paper establishes a structural analysis model of the goaf working face roof,revealing the failure mechanism of DRC,and clarifies the positive role of DRC in improving the stress of the roadway surrounding rock and reducing the subsidence of the roof through numerical simulation experiments.On this basis,the paper further analyses the roadway pressure and roof settlement under different DRC design heights,and ultimately proposes an optimized design method for the DRC height.The results indicate that the implementation of DRC can significantly optimize the stress environment of the working face roadway surrounding rock.At the same time,during the application of DRC,three scenarios may arise:insufficient,reasonable,and excessive DRC height.Insufficient height will significantly reduce the effectiveness of the technology,while excessive height has little impact on the implementation effect but will greatly increase construction costs and difficulty.Engineering verification shows that the optimized DRC design method proposed in this paper reduces the peak stress of the protective coal pillar in the roadway by 27.2%and the central subsidence of the roof by 41.8%,demonstrating excellent application results.This method provides technical support for the further promotion of NCMSE mining method.

Parallelization strategies for resolved simulations of fluid-structure-particle interactions

Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.

NUMERICAL METHOD FOR THREE-DIMENSIONAL NONLINEAR CONVECTION-DOMINATED PROBLEM OF DYNAMICS OF FLUIDS IN POROUS MEDIA

For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fractional steps techniques are needed to convert a multi-dimensional problem into a series of successive one-dimensional problems. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators, and the theory of prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources and predicting the consequences of seawater intrusion and protection projects.

High accuracy compact finite difference methods and their applications

Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.

THE UPWIND FINITE DIFFERENCE METHOD FOR MOVING BOUNDARY VALUE PROBLEM OF COUPLED SYSTEM

Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.

Characteristic Finite Difference Methods for Moving Boundary Value Problem of Numerical Simulation of Oil Deposit

The software for oil-gas transport and accumulation is to describe the history of oil-gas transport and accumulation in basin evolution. It is of great value in rational evaluation of prospecting and exploiting oil-gas resources. The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary value problem. This paper puts forward a kind of characteristic finite difference schemes, and derives from them optimal order estimates in l^2 norm for the error in the approximate solutions. The research is important both theoretically and practically for the model analysis in the field, for model numerical method and for software development.