A Novel Inverse-Free Neurodynamic Approach for Solving Absolute Value Equations

We propose a novel inverse-free neurodynamic approach (NIFNA) for solving absolute value equations (AVE). The NIFNA guarantees global convergence and notably improves convergence speed by achieving fixed-time convergence. To validate the theoretical findings, numerical simulations are conducted, demonstrating the effectiveness and efficiency of the proposed NIFNA.

Solution of Second-Order Ordinary Differential Equations via Simulated Annealing

In this paper, we approach the problem of obtaining approximate solution of second-order initial value problems by converting it to an optimization problem. It is assumed that the solution can be approximated by a polynomial. The coefficients of the polynomial are then optimized using simulated annealing technique. Numerical examples with good results show the accuracy of the proposed approach compared with some existing methods.

Modeling and simulation of JWL equation of state for reactive Al/PTFE mixture

An analytical method is presented to fit parameters of Jones-Wilkins-Lee （JWL） equation of state （EOS） for the chemical process of aluminum-polytetrafluoroethylene （ AI/PTFE ） mixture. Subroutine codes for both strength model and EOS were developed in explicit-FE code AUTODYN. Firstly, the shock Hugoniot data of reactive A1/PTFE mixture was analytically derived by implemen- ting this methodology. The JWL EOS was verified to fit shock Hugoniot data of both reacted and un- reacted A1/PTFE mixture, which gives reasonable results. Furthermore, to numerically ascertain the reaction phases of ignition and growth and quasi detonation of A1/PTFE mixture, characterized ex- periment was setup to validate the reaction phases and coefficients of JWL EOS for A1/PTFE mix- ture. From the test, a promising example of reactive mixture A1/PTFE is capable to enhance lethality of weapons, the status computation in clude quasi-detonation pressure and temperature of A1/PTFE mixture in different chemical reaction phases is validated.

An explicit method for numerical simulation of wave equations

In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.

New Numerical Scheme for Simulation of Hyperbolic Mild-Slope Equation

The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation, A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.

BGK-EQUATION SIMULATION OF THE EXHAUST PLUME FORMED BY A MANOEUVRE THRUSTER IN OUTER SPACE

The general principle of utilizing the BGK equation to simulate a macroscopic gas flow is illustrated. Two typical examples, i.e., a low-speed axisymmetric submerged jet and the Prandtl-Meyer expansion to a vacuum, are presented for validating the feasibility and accuracy of the BGK-equation simulation in continuum and non-continuum flow regimes. This approach is then used to simulate the exhaust plume formed by a small manoeuvre thruster of an artificial satellite in the outer space. The plume impingement on a flat surface perpendicular to the nozzle axis is also simulated by the same method. In the latter case the impingement force acting on the flat surface is calculated. When the flow reaches to the steady state the calculated impingement force is reasonably compared with the theoretical value of the nozzle thrust.

UNSTEADY/STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.

Approximate expression of Young's equation and molecular dynamics simulation for its applicability

In 1805, Thomas Young was the first to propose an equation(Young's equation) to predict the value of the equilibrium contact angle of a liquid on a solid. On the basis of our predecessors, we further clarify that the contact angle in Young's equation refers to the super-nano contact angle. Whether the equation is applicable to nanoscale systems remains an open question. Zhu et al. [College Phys. 4 7(1985)] obtained the most simple and convenient approximate formula, known as the Zhu–Qian approximate formula of Young's equation. Here, using molecular dynamics simulation, we test its applicability for nanodrops. Molecular dynamics simulations are performed on argon liquid cylinders placed on a solid surface under a temperature of 90 K, using Lennard–Jones potentials for the interaction between liquid molecules and between a liquid molecule and a solid molecule with the variable coefficient of strength a. Eight values of a between 0.650 and 0.825 are used. By comparison of the super-nano contact angles obtained from molecular dynamics simulation and the Zhu–Qian approximate formula of Young's equation, we find that it is qualitatively applicable for nanoscale systems.

Modified Saint-Venant equations for flow simulation in tidal rivers

Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations （the SVN model） are extensively used to model tidal rivers. A force-corrected term expressed as the combination of flow velocity and the change rate of the tidal fevel was developed to represent tidal effects in the SVN model. A momentum equation incorporating with the corrected term was derived based on Newton＇s second law. By combing the modified momentum equation with the continuity equation, an improved SVN model for tidal rivers （the ISVN model） was constructed. The simulation of a tidal reach of the Qiantang River shows that the ISVN model performs better than the SVN model. It indicates that the corrected force derived for tidal effects is reasonable; the ISVN model provides an appropriate enhancement of the SVN model for flow simulation of tidal rivers.

Direct Numerical Simulation of Interaction Between Wave and PorousBreakwater Based on N-S Equation

In this paper, a numerical model is established. A modified N-S equation is used as a control equation for the wave field and porous flow area. The control equations are discreted and solved by the finite difference method. The free surface is tracked by the VOF method. The pressure field and velocity field of the whole flow area are solved by the reiterative iteration method. Finally, compared with the physical model test results of wave flume, the numerical model established in the present study is validated.

Changes of Stress-Strain Relationship of Rabbit Femoral Vein after Simulated Weightlessness

Objective: To observe the effect of simulated weightlessness on stress-strain relationship and the structural change of rabbit femoral vein. Methods: After seting up the Head-Down Tilt (-20°) (HDT) model to simulate weightlessness, 24 healthy male New-Zealand Rabbits were randomly divided into HDT-21d group, HDT-10d group and control group, with 8 in each. Femoral venous strips and rings were used to make uniaxial tensile test of the longitudinal and circumferential specimens of the vessels. At last we observed the microstructure of femoral vein wall in 3 groups. Results : With the increasing of load stress, both longitudinal and circumferential strains of vein samples from 3 groups increased significantly (P<0. 01). With the decrease of unload stress, strains decrease obviously (P<0. 01). The unloaded longitudinal and the circumferential strain from 3 groups increased much than those of the loaded. Under the same stress (longitudinal 0-2. 0 g, circumferential 0. 5-1. 0 g) , HDT-21d group and HDT-10d group increased obviously in tlie longitudinal or circumferential strain (load and unload) than control, and HDT-21d increased much than that of HDT-10d. The contents and structures of femoral vein walls of HDT-rabbits changed significantly. Some endotheli-um cells of femoral vein became short, columnar or cubic even fell off. Smooth muscle layers became thinner. Conclusion:The compliance of femoral venous increased significantly after weightlessness-simulation and increased much obviously after 21d-HDT than that of 10 d. The structure of femoral vein wall changed obviously. The changes may be one reason for the increase of femoral vein compliance.

LaNets:Hybrid Lagrange Neural Networks for Solving Partial Differential Equations

We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural network frameworks.Concretely,we first perform Lagrange interpolation in front of the deep feedforward neural network.The Lagrange basis function has a neat structure and a strong expression ability,which is suitable to be a preprocessing tool for pre-fitting and feature extraction.Second,we introduce small sample learning into training,which is beneficial to guide themodel to be corrected quickly.Taking advantages of the theoretical support of traditional numerical method and the efficient allocation of modern machine learning,LaNets achieve higher predictive accuracy compared to the state-of-the-artwork.The stability and accuracy of the proposed algorithmare demonstrated through a series of classical numerical examples,including one-dimensional Burgers equation,onedimensional carburizing diffusion equations,two-dimensional Helmholtz equation and two-dimensional Burgers equation.Experimental results validate the robustness,effectiveness and flexibility of the proposed algorithm.

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

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 GPU-based general numerical framework for plasma simulations in terms of microscopic kinetic equations with full collision terms

We have proposed a general numerical framework for plasma simulations on graphics processing unit clusters based on microscopic kinetic equations with full collision terms.Our numerical algorithm consistently deals with both long-range(classical forces in the Vlasov term)and short-range(quantum processes in the collision term)interactions.Providing the relevant particle masses,charges and types(classical,fermionic or bosonic),as well as the external forces and the matrix elements(in the collisional integral),the algorithm consistently solves the coupled multi-particle kinetic equations.Currently,the framework is being tested and applied in the field of relativistic heavy-ion collisions;extensions to other plasma systems are straightforward.Our framework is a potential and competitive numerical platform for consistent plasma simulations.

Simplifi ed traditional bubble-motion equation and air-gun wavelet simulation based on a Van der Waals gas model

An air-gun source is the most commonly used excitation method in off shore seismic exploration.The excitation characteristics of an air-gun source aff ect seismic data quality.Far-field wavelet simulation is an important approach to study these characteristics.Compared to the measured wavelet,far-field wavelet simulation based on a traditional bubble-motion equation and ideal gas wavelet model has some disadvantages,such as a greater amplitude and smaller pulse attenuation velocity.Here,we start from the linear acoustic wave equation in the spherical coordinate system to deduce an improved,simpler bubble-motion equation and develop a Van der Waals gas wavelet model based on this equation.Unlike the existing methods,our method considers the high-pressure environment during actual excitation,heat exchange between the bubble and outside water,and change in the air fl ow at the muzzle.The results show that the far-fi eld wavelet simulated using this model is closer to the measured wavelet than that of the ideal gas wavelet model.At the same time,our method has a more succinct equation and a higher calculation effi ciency.

Dirichlet-to-Neumann Map for a Hyperbolic Equation

In this paper, we provide an explicit expression for the full Dirichlet-to-Neumann map corresponding to a radial potential for a hyperbolic differential equation in 3-dimensional. We show that the Dirichlet-Neumann operators corresponding to a potential radial have the same properties for hyperbolic differential equations as for elliptic differential equations. We numerically implement the coefficients of the explicit formulas. Moreover, a Lipschitz type stability is established near the edge of the domain by an estimation constant. That is necessary for the reconstruction of the potential from Dirichlet-to-Neumann map in the inverse problem for a hyperbolic differential equation.

Numerical Simulation of Abnormal Wave Height Change Induced by Excavated Waterway with Boussinesq Equation

A nonlinear numerical model has been set up by use of Boussinesq Equation with finite difference method, and has been applied to the simulation of the abnormal change of wave height induced by excavated waterway. Numerical results demonstrate that the abnormal change of wave height is due to the adding of the reflected wave height induced by excavated waterway to the incident wave height. Because the angle between the incident wave and the axis of the waterway is smaller than the critical angle, the reflected wave produced by the waterway may propagate to the breakwater and may be added with the incident wave, then the abnormal change of wave height before the breakwater may be caused. So the wave reflection caused by the change of water depth cannot be neglected.

A modified equation of state for Xe at high pressures by molecular dynamics simulation

The exact equation of state （EOS） for the fission gas Xe is necessary for the accurate prediction of the fission gas behavior in uranium dioxide nuclear fuel, However, the comparison with the experimental data indicates that the applicable pressure ranges of existing EOS for Xe published in the literature cannot cover the overpressure of the rim fission gas bubble at the typical UO2 fuel pellet rim structure. Based on the interatomic potential of Xe, the pressure-volume-temperature data are calculated by the molecular dynamics （MD） simulation. The results indicate that the data of MD simulation with Ross and McMahan＇s potential [M. Ross and A. K. McMahan 1980 Phys. Rev. B 21 1658] are in good agreement with the experimental data. A preferable EOS for Xe is proposed based on the MD simulation. The comparison with the MD simulation data shows that the proposed EOS can be applied at pressures up to 550 MPa and 3 GPa and temperatures 900 K and 1373 K respectively. The applicable pressure range of this EOS is wider than those of the other existing EOS for Xe published in the literature.

Finite element equations and numerical simulation of elastic wave propagation in two-phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.