Development,hotspots and trend directions of groundwater numerical simulation:A bibliometric and visualization analysis

Groundwater is a vital component of the hydrological cycle and essential for the sustainable development of ecosystems.Numerical simulation methods are key tools for addressing scientific challenges in groundwater research.This study uses bibliometric visualization analysis to examine the progress and trends in groundwater numerical simulation methods.By analyzing literature indexed in the Web of Science database from January 1990 to February 2023,and employing tools such as Citespace and VOSviewer,we assessed publication volume,research institutions and their collaborations,prolific scholars,keyword clustering,and emerging trends.The findings indicate an overall upward trend in both the number of publications and citations concerning groundwater numerical simulations.Since 2010,the number of publications has tripled compared to the total before 2010,underscoring the increasing significance and potential of numerical simulation methods in groundwater science.China,in particular,has shown remarkable growth in this field over the past decade,surpassing the United States,Canada,and Germany.This progress is closely linked to strong national support and active participation from research institutions,especially the contributions from teams at Hohai University,China University of Geosciences,and the University of Science and Technology of China.Collaboration between research teams is primarily seen between China and the United States,with less noticeable cooperation among other countries,resulting in a diverse and dispersed development pattern.Keyword analysis highlights that international research hotspots include groundwater recharge,karst water,geothermal water migration,seawater intrusion,variable density flow,contaminant and solute transport,pollution remediation,and land subsidence.Looking ahead,groundwater numerical simulations are expected to play a more prominent role in areas such as climate change,surface water-groundwater interactions,the impact of groundwater nitrates on the environment and health,submarine groundwater discharge,ecological water use,groundwater management,and risk prevention.

Furnace structure analysis for copper flash continuous smelting based on numerical simulation

According to the innate characteristic of four types of furnace, the copper flash continuous smelting （CFCS） furnace can be considered a synthetic reactor of two relatively independent processes： flash matte smelting process （FMSP） and copper continuous converting process （CCCP）. Then, the CFCS thermodynamic model was proposed by establishing the multi-phase equilibrium model of FMSP and the local-equilibrium model of CCCP, respectively, and by combining them through the smelting intermediates. Subsequently, the influences of the furnace structures were investigated using the model on the formation of blister copper, the Fe3O4 behavior, the copper loss in slag and the copper recovery rate. The results show that the type D furnace, with double flues and a slag partition wall, is an ideal CFCS reactor compared with the other three types furnaces. For CFCS, it is effective to design a partition wall in the furnace to make FMSP and CCCP perform in two relatively independent zones, respectively, and to make smelting gas and converting gas discharge from respective flues.

STATISTICAL ANALYSIS AND NUMERICAL SIMULATION OF POYANG LAKE'S INFLUENCE ON TROPICAL CYCLONES

Based on the Tropical Cyclone(TC briefly thereafter)Yearbook 1980-2009,this paper first analyzes the number and intensity change of the TCs which passed directly over or by the side of Poyang Lake(the distance of TC center is less than 1°longitude or 1°latitude from the Lake)among all the landfalling TCs in China during the past 30 years.Two cases are examined in detail in this paper.One is severe typhoon Rananim with a speed of 3.26 m/s and a change of 1 hPa in intensity when it was passing the Lake.The other is super typhoon Saomai with a faster moving speed of 6.50 m/s and a larger change in intensity of 6 hPa.Through numerical simulation experiments,this paper analyzes how the change of underlying surface from water to land contributes to the differences in intensity,speed and mesoscale convection of the two TCs when they passed the Lake.Results show that the moisture and dynamic condition above the Lake were favorable for the maintenance of the intensity when Rananim was passing through Poyang Lake,despite the moisture supply from the ocean was cut off.As a result,there was strong convection around the lake which led to a rainfall spinning counter-clockwise as it was affected by the TC movement.However,little impact was seen in the Saomai case.These results indicate that for the TCs coming ashore on Poyang Lake with a slow speed,the large water body is conducive to the sustaining of the intensity and strengthening of the convection around the TC center and the subsequent heavy rainfall.On the contrary,a fast-moving TC is less likely to be influenced by the underlying surface in terms of intensity and speed.

Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

NUMERICAL SIMULATION AND MESOSCALE ANALYSIS OF A TORRENTIAL RAIN CAUSED BY TYPHOON

A heavy storm rainfall caused by Typhoon Aere (No.0418) when landing at Fujian has been successfully simulated by using AREM model. The simulation result is scale-separated by spatial band-pass filtering, which reveals the mesoscale low pressure and convergence line that has direct impact on this rainfall process. The physical characteristics of the two mesoscale systems and their relation with rainfall are also analyzed. Study shows that there exists a well corresponding relationship between the storm rainfall and mesoscale divergence and strong updraft arising from the convergence, which is caused by the interactions between the mesoscale systems and topographic features, and is directly responsible for the rainfall.

Three-dimensional discrete element numerical simulation of Paleogene salt structures in the western Kuqa foreland thrust belt

Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.

Reliability Analysis of Numerical Simulation of Sea-Crossing Bridge Deformation Under Wave Forces

To study on the numerical simulation calculation reliability of sea-crossing bridge under complex wave forces,the paper applied GPS deformation monitoring and numerical simulation calculation by researching Qingdao Jiaozhou Bay Sea-Crossing Bridge.The db3 wavelet three-layer decomposition was used on the horizontal movement of the sea-crossing bridge and the wind speed of the waves to analyze their correlation.The complex wave forces value of Qingdao Jiaozhou Bay Sea-Crossing Bridge was loaded on FLAC3D software successfully to make numerical simulation calculation of bridge deformation.Since the accuracy of the GPS deformation monitoring reaches millimeter level,it was used to monitor the exact value of the bridge deformation to judge the reliability of numerical simulation.The relative errors of displacement in X,Y and Z directions were between 33%and 41%through comparison.It could be seen that the numerical simulation error was relatively large,which was mainly due to various environmental factors and the deviation of applied wave forces.However,numerical simulation generally reflects the deformation law of the sea-crossing bridge under complex wave forces,providing an effectively technical support for the safe operation assessment of the sea-crossing bridge.

Numerical simulation of three-dimensional tidal current in the Bobal Sea

The model equations with tbree-dimensional, time-dependent, nonlinear Navier-Stokes equations are transformed by sigma-transformation.On the basis of the process splitting technique, the fluid flow problems are divided into two parts:the vertically-intopated equations (external mode) and the vertical structure equations(internal mace). The first set of equations being the propagation of the tidal weves and the ADI numerical scheme has ben chosen to solve them. Conerning the vertical structure equations, they are solved by means of leapfrog stepping procedure.The main features of the tide and associated tidal current in the Bohai Sea are examined with this 3-D model.To have a good reproduction of vertical structure, the column is divided into 10 layers and the M2 tidal current is computed in detail. The simulation reveal the spetial structure and some important characteristics of the tidal current of the Bohai Sea. The application of the 3-D madel to forecasting of the tidal current in the Bobal Sea has been Performed as an illustration.

CAREFUL NUMERICAL SIMULATION AND ANALYSIS OF MIGRATION-ACCUMULATION OF TANHAI REGION

Numerical simulation of careful parallel arithmetic of oil resources migration-accumulation of Tanhai Region ( three-layer) was done. Careful parallel operator splitting-up implicit iterative scheme, parallel arithmetic program, parallel arithmetic information and alternating-direction mesh subdivision were put forward. Parallel arithmetic and analysis of different CPU combinations were done. This numerical simulation test and the actual conditions are basically coincident. The convergence estimation of the model problem has successfully solved the difficult problem in the fields of permeation fluid mechanics, computational mathematics and petroleum geology.

Numerical Simulation of the Three-Dimensional Wave-Induced Currents on Unstructured Grid

By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.

Numerical analysis of hydraulic fracture propagation in deep shale reservoir with different injection strategies

Deep shale reservoirs are characterized by elevated breakdown pressures,diminished fracture complexity,and reduced modified volumes compared to medium and shallow reservoirs.Therefore,it is urgent to investigate particular injection strategies that can optimize breakdown pressure and fracturing efficiency to address the increasing demands for deep shale reservoir stimulation.In this study,the efficiency of various stimulation strategies,including multi-cluster simultaneous fracturing,modified alternating fracturing,alternating shut-in fracturing,and cyclic alternating fracturing,was evaluated.Subsequently,the sensitivity of factors such as the cycle index,shut-in time,cluster spacing,and horizontal permeability was investigated.Additionally,the flow distribution effect within the wellbore was discussed.The results indicate that relative to multi-cluster simultaneous fracturing,modified alternating fracturing exhibits reduced susceptibility to the stress shadow effect,which results in earlier breakdown,extended hydraulic fracture lengths,and more consistent propagation despite an increase in breakdown pressure.The alternating shut-in fracturing benefits the increase of fracture length,which is closely related to the shut-in time.Furthermore,cyclic alternating fracturing markedly lowers breakdown pressure and contributes to uniform fracture propagation,in which the cycle count plays an important role.Modified alternating fracturing demonstrates insensitivity to variations in cluster spacing,whereas horizontal permeability is a critical factor affecting fracture length.The wellbore effect restrains the accumulation of pressure and flow near the perforation,delaying the initiation of hydraulic fractures.The simulation results can provide valuable numerical insights for optimizing injection strategies for deep shale hydraulic fracturing.

A three-dimensional Eulerian method for the numerical simulation of high-velocity impact problems

In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.

Operation optimization of prefabricated light modular radiant heating system:Thermal resistance analysis and numerical study

The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.

Numerical simulation study on the mold strength of magnetic mold casting based on a coupled electromagnetic-structural method

The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.

Numerical Simulation and Parallel Computing of Acoustic Wave Equation in Isotropic-Heterogeneous Media

In this paper,we consider the numerical implementation of the 2D wave equation in isotropic-heterogeneous media.The stability analysis of the scheme using the von Neumann stability method has been studied.We conducted a study on modeling the propagation of acoustic waves in a heterogeneous medium and performed numerical simulations in various heterogeneous media at different time steps.Developed parallel code using Compute Unified Device Architecture(CUDA)technology and tested on domains of various sizes.Performance analysis showed that our parallel approach showed significant speedup compared to sequential code on the Central Processing Unit(CPU).The proposed parallel visualization simulator can be an important tool for numerous wave control systems in engineering practice.

High-Precision Flow Numerical Simulation and Productivity Evaluation of Shale Oil Considering Stress Sensitivity

Continental shale oil reservoirs,characterized by numerous bedding planes and micro-nano scale pores,feature significantly higher stress sensitivity compared to other types of reservoirs.However,research on suitable stress sensitivity characterization models is still limited.In this study,three commonly used stress sensitivity models for shale oil reservoirs were considered,and experiments on representative core samples were conducted.By fitting and comparing the data,the“exponential model”was identified as a characterization model that accurately represents stress sensitivity in continental shale oil reservoirs.To validate the accuracy of the model,a two-phase seepage mathematical model for shale oil reservoirs coupled with the exponential model was introduced.The model was discretely solved using the finite volume method,and its accuracy was verified through the commercial simulator CMG.The study evaluated the productivity of a typical horizontal well under different engineering,geological,and fracture conditions.The results indicate that considering stress sensitivity leads to a 13.57%reduction in production for the same matrix permeability.Additionally,as the fracture half-length and the number of fractures increase,and the bottomhole flowing pressure decreases,the reservoir stress sensitivity becomes higher.

Simulation Method and Feature Analysis of Shutdown Pressure Evolution During Multi-Cluster Fracturing Stimulation

Multistage multi-cluster hydraulic fracturing has enabled the economic exploitation of shale reservoirs,but the interpretation of hydraulic fracture parameters is challenging.The pressure signals after pump shutdown are influenced by hydraulic fractures,which can reflect the geometric features of hydraulic fracture.The shutdown pressure can be used to interpret the hydraulic fracture parameters in a real-time and cost-effective manner.In this paper,a mathematical model for shutdown pressure evolution is developed considering the effects of wellbore friction,perforation friction and fluid loss in fractures.An efficient numerical simulation method is established by using the method of characteristics.Based on this method,the impacts of fracture half-length,fracture height,opened cluster and perforation number,and filtration coefficient on the evolution of shutdown pressure are analyzed.The results indicate that a larger fracture half-length may hasten the decay of shutdown pressure,while a larger fracture height can slow down the decay of shutdown pressure.A smaller number of opened clusters and perforations can significantly increase the perforation friction and decrease the overall level of shutdown pressure.A larger filtration coefficient may accelerate the fluid filtration in the fracture and hasten the drop of the shutdown pressure.The simulation method of shutdown pressure,as well as the analysis results,has important implications for the interpretation of hydraulic fracture parameters.

Analysis of cranial aneurysm's 3D-computational hemodynamics numerical simulation and findings in operation

Objective With the cranial aneurysm’s 3D-computational hemodynamics numerical simulation,we can get the cranial aneurysms and its ’s mother arteries ’blood flow rate,path and wall shearing stress and so on. Taking use of the contrast analysis of the hemodynamics result and the findings in operation,we try to find hemodynamic character and mechanism of rupture

The numerical simulation for a 3D two-phase anisotropic medium based on BISQ model

Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.

Numerical Simulation of a Spark Ignited Two-Stroke Free-Piston Engine Generator

A numerical program is built to simulate the performance of a spark ignited two-stroke free-piston engine coupled with a linear generator. The computational model combines a series of dynamic and thermodynamic equations that are solved simultaneously to predict the performances of the engines. The dynamic analysis performed consists of an evaluation of the frictional force and load force introduced by the generator. The thermodynamic analysis used a single zone model to describe the engine＇ s working cycle which includes intake, scavenging, compression, combustion and expansion, and to evaluate the effect of heat transfer based on the first law of thermodynamics and the ideal gas state equation. Because there is no crankshaft, a time based Wiebe equation was used to express the fraction of fuel burned in the combustion. The calculated results were validated by using the experimental data from another research group. The results indicate that the free-piston generator has some advantages over conventional engines.