Symplectic partitioned Runge-Kutta method based onthe eighth-order nearly analytic discrete operator and its wavefield simulations

We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.

A spatiotemporal evolution model of a short-circuit arc to a secondary arc based on the improved charge simulation method

The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.

Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.

Risk assessment of coastal flooding disaster by storm surge based on Elevation-Area method and hydrodynamic models:Taking Bohai Bay as an example

The future inundation by storm surge on coastal areas are currently ill-defined.With increasing global sealevel due to climate change,the coastal flooding by storm surge is more and more frequently,especially in coastal lowland with land subsidence.Therefore,the risk assessment of such inundation for these areas is of great significance for the sustainable socio-economic development.In this paper,the authors use Elevation-Area method and Regional Ocean Model System(ROMS)model to assess the risk of the inundation of Bohai Bay by storm surge.The simulation results of Elevation-Area method show that either a 50-year or 100-year storm surge can inundate coastal areas exceeding 8000 km^(2);the numerical simulation results based on hydrodynamics,considering ground friction and duration of the storm surge high water,show that a 50-year or 100-year storm surge can only inundate an area of over 2000 km^(2),which is far less than 8000 km^(2);while,when taking into account the land subsidence and sea level rise,the very inundation range will rapidly increase by 2050 and 2100.The storm surge will greatly impact the coastal area within about 10-30 km of the Bohai Bay,in where almost all major coastal projects are located.The prompt response to flood disaster due to storm surge is urgently needed,for which five suggestions have been proposed based on the geological background of Bohai Bay.This study may offer insight into the development of the response and adaptive plans for flooding disasters caused by storm surge.

Research on Simulation Methods of Electric Field Intensity on Surface of 10 kV Cable Joint

The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field-circuit coupling method and equivalent circuit method. The average EFI of the inner surface of the outer semi-conducting shield can be calculated from the current in the measuring circuit. The relative error between these two methods is about 15%, which roughly proves the consistency of the two methods. Further practical application research enables online monitoring of cable joints.

Model of Universe as Described by Dynamic Universe Model

In this paper we will see the model of Universe according to Dynamic Universe Model of Cosmology by visualizing various processes that are happening in the Universe as per experimental evidences. For simplifying the matter here, we will see in part 1: about the Galaxy life cycle, where the birth and death of Galaxies discussed. Probably Universe gives guidance for the movement of Galaxies. We call this Part 1: Thinking and Reproducing Universe or Mindless Universe? (Galaxy life cycle). We see every day Sun, Stars, Galaxies etc., dissipating enormous energy in the form of radiation by the way of fusion of Hydrogen to helium. So after sometime all the Hydrogen is spent and Universe will die, is it not? … Dynamic Universe Model says that the energy in the form of electromagnetic radiation passing grazingly near any gravitating mass changes in frequency and finally will convert into neutrinos (mass). Hence Dynamic Universe Model proposes another process where energy will be converted back into matter and the cycle energy to mass to energy continues, sustaining the Universe to maintain this present status for ever in this form something like a Steady state model without any expansion. This we will see in Part 2: Energy - Mass - Energy Cycle. After converting energy into mass “how various elements are formed and where they are formed?” will be next logical question. Dynamic Universe Model says that these various particles change into higher massive particles or may get bombarded into stars or planets and various elements are formed. Here we bifurcate the formation of elements into 6 processes. They are for Elementary particles and elements generated in frequency changing process, By Cosmic rays, By Small stars, By Large Stars, By Super Novae and Manmade elements By Neutron Stars. This we will discuss in Part 3: Nucleosynthesis.

Numerical simulation on residual stress and deformation for WAAM parts of aluminum alloy based on temperature function method

Wire arc additive manufacture(WAAM) is a new technique to fabricate large-scale complex aluminum alloy components.However, the performance of the parts is critically influenced by residual stresses and deformation. A sequentially thermal-mechanical coupled model of residual stress and deformation for aluminum alloy WAAM parts was established based on commercial FE software ABAQUS. The temperature field was calculated by the moving heat source(MHS) method. The temperature function was obtained according to the distribution of the peak temperature. Furthermore, the MHS method and segmented temperature function(STF) method were used to calculate the residual stress and deformation. The results show that the STF method satisfies both the efficiency and accuracy requirements. 1-segment, 3-segment, and 5-segment STF methods can shorten the time for mechanical analysis by 91%, 79%, 63%, respectively.The error of the residual stress and deformation are all less than 20%. STF method provides an effective way to predict the residual stress and deformation of large-scale WAAM parts.

Modified slow-fast analysis method for slow-fast dynamical systems with two scales in frequency domain

A modified slow-fast analysis method is presented for the periodically excited non-autonomous dynamical system with an order gap between the exciting frequency and the natural frequency.By regarding the exciting term as a slow-varying parameter,a generalized autonomous fast subsystem can be defined,the equilibrium branches as well as the bifurcations of which can be employed to account for the mechanism of the bursting oscillations by combining the transformed phase portrait introduced.As an example,a typical periodically excited Hartley model is used to demonstrate the validness of the method,in which the exciting frequency is far less than the natural frequency.The equilibrium branches and their bifurcations of the fast subsystem with the variation of the slow-varying parameter are presented.Bursting oscillations for two typical cases are considered,which reveals that,fold bifurcation may cause the the trajectory to jump between different equilibrium branches,while Hopf bifurcation may cause the trajectory to oscillate around the stable limit cycle.

基于遗传算法和最小二乘支持向量机的织物剪切性能预测

提出了一种基于最小二乘支持向量机的织物剪切性能预测模型,并且采用遗传算法进行最小二乘支持向量机的参数优化,将获得的样本进行归一化处理后,将其输入预测模型以得到预测结果.仿真结果表明,基于最小二乘支持向量机的预测模型比BP神经网络和线性回归方法具有更高的精度和范化能力. Abstract： A new method is proposed to predict the fabric shearing property with least square support vector machines ( LS-SVM ). The genetic algorithm is investigated to select the parameters of LS-SVM models as a means of improving the LS- SVM prediction. After normalizing the sampling data, the sampling data are inputted into the model to gain the prediction result. The simulation results show the prediction model gives better forecasting accuracy and generalization ability than BP neural network and linear regression method.

Numerical Simulation of Boundary Layer with Spectral Method

It is important to improve the speed of a ship,the friction resistance can be reduced by injection air at the bottom of a ship when the ship is running on the water.As the first part of the studying project,here numerical simulation study method,boundary condition and governing equations are presented.It is easy to study complicated problems from simple conditions,so the program concerning boundary layer condition is compiled to solve the problem.Here the spectral method is introduced,and the results are tested by Dorod’s results.

Simulation of Flow Field of Molten Salt in Neodymium Metal Electrolytic Cell Using Vortex-Flow Function Method

With the applications of Nd-Fe-B material extending in recent years, the materials of neodymium metal and other rare earth metal alloy confront the increased demand and the high quality request at the same time.These factors stimulated greatly to perfect the producing craft of RE metals and improve the equipments.The rare earth electrolysis cell is developing towards large-scale way.Notwithstanding the present electrolysis cell of Nd metal, include 6 kA and 10 kA cell, exists some insurmountable problems during operation and these problems lead to lower electric efficiency and higher operating costs.So it is significant to study the physical fields of rare earth electrolysis cell.In this paper,a numerical flow mode is established using vortex- flowing function method and the fluid flow field of 3000A Nd electrolysis cell is computed using MATLAB.The results of the study will be important reference in theory for improving and enlarging rare earth fluoride system cell.

An Isogeometric Cloth Simulation Based on Fast Projection Method

A novel continuum-based fast projection scheme is proposed for cloth simulation.Cloth geometry is described by NURBS,and the dynamic response is modeled by a displacement-only Kirchhoff-Love shell element formulated directly on NURBS geometry.The fast projection method,which solves strain limiting as a constrained Lagrange problem,is extended to the continuum version.Numerical examples are studied to demonstrate the performance of the current scheme.The proposed approach can be applied to grids of arbitrary topology and can eliminate unrealistic over-stretching efficiently if compared to spring-based methodologies.

Calculation Method of Passive Residual Heat Removal Heat Exchanger and Numerical Simulation

The tube inside and outside heat transfer mechanism of Passive Residual Heat Removal Heat Exchanger (PRHR HX) was analyzed. The calculation method of this special heat exchanger under natural convection condition in In-containment Refueling Water Storage Tank (IRWST) was carried out. The single-tube coupling model three-dimensional natural circulation in the IRWST was simulated numerically using Fluent. The heat transfer and flow characteristics of the fluid in IRWST were obtained. The comparison of the results between theoretical arithmetic and numerical simulation showed that the theoretical calculation method is suitable for the heat transfer calculation of PRHR HX.

Method“Monte Carlo”in healthcare

In public health,simulation modeling stands as an invaluable asset,enabling the evaluation of new systems without their physical implementation,experimentation with existing systems without operational adjustments,and testing system limits without real-world repercussions.In simulation modeling,the Monte Carlo method emerges as a powerful yet underutilized tool.Although the Monte Carlo method has not yet gained widespread prominence in healthcare,its technological capabilities hold promise for substantial cost reduction and risk mitigation.In this review article,we aimed to explore the transformative potential of the Monte Carlo method in healthcare contexts.We underscore the significance of experiential insights derived from simulated experimentation,especially in resource-constrained scenarios where time,financial constraints,and limited resources necessitate innovative and efficient approaches.As public health faces increasing challenges,incorporating the Monte Carlo method presents an opportunity for enhanced system construction,analysis,and evaluation.

Review of the algebraic linear methods and parallel implementation in numerical simulation of groundwater flow

The desire to increase spatial and temporal resolution in modeling groundwater system has led to the requirement for intensive computational ability and large memory space. In the course of satisfying such requirement, parallel computing has played a core role over the past several decades. This paper reviews the parallel algebraic linear solution methods and the parallel implementation technologies for groundwater simulation. This work is carried out to provide guidance to enable modelers of groundwater systems to make sensible choices when developing solution methods based upon the current state of knowledge in parallel computing.

COMPUTER SIMULATION OF STRUCTURE OF MOLTEN DyF_3－BaF_2－LiF SYSTEM BY MONTE CARLO METHOD

ＣＯＭＰＵＴＥＲＳＩＭＵＬＡＴＩＯＮＯＦＳＴＲＵＣＴＵＲＥＯＦＭＯＬＴＥＮＤｙＦ_３－ＢａＦ_２－ＬｉＦＳＹＳＴＥＭＢＹＭＯＮＴＥＣＡＲＬＯＭＥＴＨＯＤ￥ＸｕＣｈｉ；ＳｕＨａｎｇ；ＣｈｅｎＮｉａｎｙｉ（ＳｈａｎｇｈａｉＩｎｓｔｉｔｕｔｅｏｆＭｅｔａｌｌ...

Modeling and Simulation of Particle-Packing Structures and Their Stability Using the Distinct Element Method

A numerical method for simulating the stability of particle-packing structures is presented. The packing structures were modeled on the basis of face-centered cubic (fcc) and body-centered cubic (bcc) structures, and the stability of these structures was investigated using the distinct element method. The interaction between the particles was simplified by considering repulsive, adhesive, and damping forces, and the stability against the gravitational force was simulated. The results under a certain set of parameters showed characteristic deformation when the particles were arranged in an fcc array. Focusing on the local structure, the resulting model was divided into several domains: The bottom base, four top corners, and intermediate domains. The bottom base notably became a body-centered tetragonal (bct) structure, which corresponds to a uniaxially compressed bcc structure. Conversely, the models based on the bcc arrangement were structurally stable, as no specific deformation was observed, and a monotonously compressed bct structure was obtained. Consequently, the bcc arrangement is concluded to be more stable against uniaxial compression, such as the gravitational force, in a particle-packing system.

基于FAST和Sobol指数法的雷达系统效能敏感性分析

为了以较小的代价快速提升对空探测雷达系统的综合效能,首先搭建了对空探测雷达系统综合效能的评估指标体系,在此基础上完成了综合效能评估模型的建立。选择两种基于方差原理的全局敏感性分析方法,即傅里叶幅度敏感性测试(Fourier amplitudes sensitivity test,FAST)和Sobol指数法,分别从频域和时域对建立的评估模型进行敏感性分析。在基于FAST进行敏感性分析时,提出了21维输入指标所对应的特征频率集的选取方法;在基于Sobol指数法进行敏感性分析时,考虑了二阶及高阶敏感性指数的影响。对比FAST和Sobol指数法的分析结果,两种方法获得的敏感性分析结果具有较好的一致性。其中,功率和频率为对空探测雷达系统的主要影响指标,为雷达系统的研发和升级提供了可靠性依据。

Numerical Simulation of Molten Metal Droplet Behavior in Gas Metal Arc Welding by Three-Dimensional Incompressible Smoothed Particle Hydrodynamics Method

The numerical model was developed using a SPH (Smoothed Particle Hydrodynamics) method and the projected transfer phenomena during a GMA (Gas Metal Arc) welding were simulated by the model to clarify mechanisms of the phenomena. As a result, the droplet transfer mode obtained from this calculation was regarded as a projected transfer mode in which the liquid column grew about 1 mm and a droplet grew up until its diameter became large the same as a wire diameter,?after that it was detached from the tip of the column. In addition, 10 droplets were formed for 0.1 s through these growth and detachment processes at the tip of a wire. To compare with the numerical results, actual GMA welding was carried out and molten metal droplet transfers were taken by high speed camera. The diameter of a wire, the length of a liquid column, the velocity of a droplet right before it reached a weld pool obtained by simulation showed good agreement with experiment.