纺粘非织造喷射流场理论模型有限差分方程的建立

着重讨论了纺粘牵伸器喷射流场理论模型有限差分方程的建立,首先建立初步性的离散方程和通用化的离散方程,最后结合课题,以一维稳态有源项的对流—扩散问题,来分析如何建立纺粘二维喷射流场理论模型的有限差分方程,并推导二维稳态问题的离散化方程,同时简要介绍了三维问题的离散化方程。

A Hybrid Dung Beetle Optimization Algorithm with Simulated Annealing for the Numerical Modeling of Asymmetric Wave Equations

In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.

Fracture identification based on remote detection acoustic reflection logging

Fracture identification is important for the evaluation of carbonate reservoirs. However, conventional logging equipment has small depth of investigation and cannot detect rock fractures more than three meters away from the borehole. Remote acoustic logging uses phase-controlled array-transmitting and long sound probes that increase the depth of investigation. The interpretation of logging data with respect to fractures is typically guided by practical experience rather than theory and is often ambiguous. We use remote acoustic reflection logging data and high-order finite-difference approximations in the forward modeling and prestack reverse-time migration to image fractures. First, we perform forward modeling of the fracture responses as a function of the fracture-borehole wall distance, aperture, and dip angle. Second, we extract the energy intensity within the imaging area to determine whether the fracture can be identified as the formation velocity is varied. Finally, we evaluate the effect of the fracture-borehole distance, fracture aperture, and dip angle on fracture identification.

Three-dimensional acoustic wave equation modeling based on the optimal finite-difference scheme

Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.

Stability of Difference Systems with Finite Delay

In this paper, the authors establish some theorems that can ascertain the zero solutions of systemsx(n+1)=f(n,x n)(1)are uniformly stable,asymptotically stable or uniformly asymptotically stable. In the obtained theorems, ΔV is not required to be always negative, where ΔV(n,x n)≡V(n+1,x(n+1)) -V(n,x(n))=V(n+1,f(n,x n))-V(n,x(n)), especially, in Theorem 1, ΔV may be even positive, which greatly improve the known results and are more convenient to use.

New numerical algorithm of gas-liquid two-phase flow considering characteristics of liquid metal during mold filling

A new program is developed for gas-liquid two-phase mold filling simulation in casting. The gas fluid, the superheated liquid metal and the liquid metal containing solid grains are assumed to be governed by Navier-Stokes equations and solved through Projection method. The Level set method is used to track the gas-liquid interface boundary. In order to demonstrate the correctness of this new program for simulation of gas-liquid two-phase mold filling in casting, a benchmark filling experiment is simulated （this benchmark test is designed by XU and the filling process is recorded by a 16-mm film camera）. The simulated results agree very well with the experimental results, showing that this new program can be used to properly predicate the gas-liquid two-phase mold filling simulation in casting.

HIGH RESOLUTION POSITIVITY-PRESERVING DIFFERENCE SCHEMES FOR TWO DIMENSIONAL EULER EQUATIONS

A class of high resolution positivity preserving Boltzmann type difference schemes for one and two dimensional Euler equations is studied. First, the relation between Boltzmann and Euler equations is analyzed. By using a kind of special interpolation, the high resolution Boltzmann type difference scheme is constructed. Finally, numerical tests show that the schemes are effective and useful.

Pure quasi-P wave equation and numerical solution in 3D TTI media

Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis （VTI media）, a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis （TTI media）. This is achieved using projection transformation, which rotates the direction vector in the coordinate system of observation toward the direction vector for the coordinate system in which the z-component is parallel to the symmetry axis of the TTI media. The equation has a simple form, is easily calculated, is not influenced by the pseudo-shear wave, and can be calculated reliably when δ is greater than ε. The finite difference method is used to solve the equation. In addition, a perfectly matched layer （PML） absorbing boundary condition is obtained for the equation. Theoretical analysis and numerical simulation results with forward modeling prove that the equation can accurately simulate a quasi-P wave in TTI medium.

Improved finite difference method for pressure distribution of aerostatic bearing

An improved finite difference method （FDM）is described to solve existing problems such as low efficiency and poor convergence performance in the traditional method adopted to derive the pressure distribution of aerostatic bearings. A detailed theoretical analysis of the pressure distribution of the orifice-compensated aerostatic journal bearing is presented. The nonlinear dimensionless Reynolds equation of the aerostatic journal bearing is solved by the finite difference method. Based on the principle of flow equilibrium, a new iterative algorithm named the variable step size successive approximation method is presented to adjust the pressure at the orifice in the iterative process and enhance the efficiency and convergence performance of the algorithm. A general program is developed to analyze the pressure distribution of the aerostatic journal bearing by Matlab tool. The results show that the improved finite difference method is highly effective, reliable, stable, and convergent. Even when very thin gas film thicknesses （less than 2 Win）are considered, the improved calculation method still yields a result and converges fast.

Development and experimental validation of a one-dimensional dynamic hygrothermal modeling based on air humidity ratio

A modified one-dimensional transient hygrothermal model for multilayer wall was proposed using air humidity ratio and temperature as the driving potentials.The solution for the governing equations was obtained numerically by implementing the finite-difference scheme.To evaluate the accuracy of the model,a test system was built up to measure relative humidity and temperature within a porous wall and compare with the prediction of the model.The prediction results have good agreement with the experimental results.For the interface close to indoor side,the maximum deviation of temperature between calculated and test data is 1.87 K,and the average deviation is 0.95 K;the maximum deviation of relative humidity is 11.4%,and the average deviation is 5.7%.For the interface close to outdoor side,the maximum deviation of temperature between prediction and measurement is 1.78 K,and the average deviation is 1.1 K;the maximum deviation of relative humidity is 9.9%,and the average deviation is 4.2%.

The Stability Research for the Finite Difference Scheme of a Nonlinear Reaction-diffusion Equation

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.

PML and CFS-PML boundary conditions for a mesh-free finite difference solution of the elastic wave equation

Mesh-free finite difference(FD)methods can improve the geometric flexibility of modeling without the need for lattice mapping or complex meshing process.Radial-basisfunction-generated FD is among the most commonly used mesh-free FD methods and can accurately simulate seismic wave propagation in the non-rectangular computational domain.In this paper,we propose a perfectly matched layer(PML)boundary condition for a meshfree FD solution of the elastic wave equation,which can be applied to the boundaries of the non-rectangular velocity model.The performance of the PML is,however,severely reduced for near-grazing incident waves and low-frequency waves.We thus also propose the complexfrequency-shifted PML(CFS-PML)boundary condition for a mesh-free FD solution of the elastic wave equation.For two PML boundary conditions,we derive unsplit time-domain expressions by constructing auxiliary differential equations,both of which require less memory and are easy for programming.Numerical experiments demonstrate that these two PML boundary conditions effectively eliminate artificial boundary reflections in mesh-free FD simulations.When compared with the PML boundary condition,the CFS-PML boundary condition results in better absorption for near-grazing incident waves and evanescent waves.

Improving the Stability Problem of the Finite Difference Scheme for Reaction-diffusion Equation

This paper deals with the special nonlinear reaction-diffusion equation. The finite difference scheme with incremental unknowns approximating to the differential equation （2.1） is set up by means of introducing incremental unknowns methods. Through the stability analyzing for the scheme, it was shown that the stability conditions of the finite difference schemes with the incremental unknowns are greatly improved when compared with the stability conditions of the corresponding classic difference scheme.

3-D acoustic wave equation forward modeling with topography

In order to model the seismic wave field with surface topography, we present a method of transforming curved grids into rectangular grids in two different coordinate systems. Then the 3D wave equation in the transformed coordinate system is derived. The wave field is modeled using the finite-difference method in the transformed coordinate system. The model calculation shows that this method is able to model the seismic wave field with fluctuating surface topography and achieve good results. Finally, the energy curves of the direct and reflected waves are analyzed to show that surface topography has a great influence on the seismic wave＇s dynamic properties.

Uniform Convergence Analysis of Finite Difference Scheme for Singularly Perturbed Delay Differential Equation on an Adaptively Generated Grid

Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.

Analyzing Parabolic Profile Path for Underwater Towed-Cable

This article discusses the dynamic state analysis of underwater towed-cable when tow-ship changes its speed in a direction making parabolic profile path. A three-dimensional model of underwater towed system is studied. The established governing equations for the system have been solved using the central implicit finite-difference method. The obtained difference non-linear coupled equations are solved by Newton＇s method and satisfactory results were achieved. The solution of this problem has practical importance in the estimation of dynamic loading and motion, and hence it is directly applicable to the enhancement of safety and the effectiveness of the offshore activities.

Continuum and Discrete Element Coupling Approach to Analyzing Seismic Responses of a Slope Covered by Deposits

Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problematic, however, to extend these approaches to the simulation of the slide, collapse and disintegration of the deposits under seismic loading. Contrary to this, a discrete element method （DEM） provides a means to consider large displacement and rotation of the non-continuum. To take the advantages of both methods of continuum and non- continuum analyses, seismic responses of a slope covered by deposits are studied by coupling a twodimensional （a-D） finite difference method and a 2-D DEM, with the bedrock being modelled by the finite difference grids and the deposits being represented by disks. A smooth transition across the boundaries of the continuous/discontinuous domains is obtained by imposing the compatibility condition and equilibrium condition along their interfaces. In the course of computation, the same time-step value is chosen for both continuous and discontinuous domains. The free-field boundaries are adopted for lateral grids of bedrock domain to eliminate the radiation damping effect. When the static equilibrium under gravity load is obtained, dynamic calculation begins under excitation of the seismic wave input from the continuum model bottom. In this way, responses to the earthquake of a slope covered by deposits are analyzed dynamically. Combined with field monitoring data, deformation and stability of the slope are discussed. The effects of the relevant parameters of spectrum characteristic, duration, andpeak acceleration of seismic waves are further investigated and explained from the simulations.

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.

Transient analysis of Casson fluid thermo-convection from a vertical cylinder embedded in a porous medium: Entropy generation and thermal energy transfer visualization

Thermal transport in porous media has stimulated substantial interest in engineering sciences due to increasing applications in filtration systems,porous bearings,porous layer insulation,biomechanics,geomechanics etc.Motivated by such applications,in this article,a numerical study of entropy generation impacts on the heat and momentum transfer in time-dependent laminar incompressible boundary layer flow of a Casson viscoplastic fluid over a uniformly heated vertical cylinder embedded in a porous medium is presented.Darcy’s law is used to simulate bulk drag effects at low Reynolds number for an isotropic,homogenous porous medium.Heat line visualization is also included.The mathematical model is derived and normalized using appropriate transformation variables.The resulting non-linear time-dependent coupled governing equations with associated boundary conditions are solved via an implicit finite difference method which is efficient and unconditionally stable.The outcomes show that entropy generation and Bejan number are both elevated with increasing values of Darcy number,Casson fluid parameter,group parameter and Grashof number.To analyze the heat transfer process in a two-dimensional domain,plotting heat lines provides an excellent approach in addition to streamlines and isotherms.It is remarked that as the Darcy number increases,the deviations of heat lines from the hot wall are reduced.

An Investigation of Buckled of Column on Liner Elastic foundation by Using Finite-Difference Mode

The use of columns on elastic foundation is very common in Civil Engineering, like bridge pier, the foundation of the buildings etc. So, it will be useful to find the critical load for the structure, the problem in this paper will be solved by Finite-Difference Mode, that＇ s simple and has an extensive use. The way it works is that by dividing the component into many units. Finite-difference methods （FDM） are numerical methods for anoroximating, the solutions to differential eauations usine finite difference equations to approximate derivatives.