Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-s...As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.展开更多
After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir develo...After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.展开更多
At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures...At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.展开更多
Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce f...Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.展开更多
The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critica...The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.展开更多
The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improv...The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improved stochastic finite-fault method to systematically assess seismic impacts.Observed near-field recordings at MM.NGU station was used to determine the reliability of the theoretically derived stress drop as input for simulation.Far-field recordings constrained the frequency-dependent S-wave quality factors(Q(f)=283.305f^(0.588))for anelastic attenuation modeling.Comparisons of peak accelerations between simulation and empirical ground-motion models showed good agreement at moderate-to-large distances.However,lower near-fault simulations indicate a weaker-than-average source effect.Analysis of simulated instrumental seismic intensity revealed key patterns.Maximum intensity(Ⅹ)occurred in isolated patches within the ruptured fault projection,correlating with shallow high-slip areas.TheⅨ-intensity zone formed a north-south elongated band centered on fault projection.Significant asymmetry inⅧ-intensity distribution perpendicular to the fault strike was observed,with a wider western extension attributed to lower shear-wave velocities west of the fault.Supershear rupture behavior enhanced ground motions,expanding intensity ranges by~20%compared to sub-shear rupture.This study reveals the integrated effects of fault geometry,slip spatial distribution,rupture velocity,and site condition in governing ground motion patterns.展开更多
0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation ph...0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).展开更多
Iterative coupled methods are widely used in multi-fidelity simulation of rotating components due to the simple implementation,which iteratively eliminates the errors between the computational fluid dynamics models an...Iterative coupled methods are widely used in multi-fidelity simulation of rotating components due to the simple implementation,which iteratively eliminates the errors between the computational fluid dynamics models and approximate characteristic maps.However,the convergence and accuracy of the iterative coupled method are trapped in characteristic maps.In particular,iterative steps increase sharply as the operation point moves away from the design point.To address these problems,this paper developed an auxiliary iterative coupled method that introduces the static-pressure-auxiliary characteristic maps and modification factor of mass flow into the component-level model.The developed auxiliary method realized the direct transfer of static pressure between the high-fidelity models and the component-level model.Multi-fidelity simulations of the throttle characteristics were carried out using both the auxiliary and traditional iterative coupled methods,and the simulation results were verified using the experimental data.Additionally,the consistency between the auxiliary and traditional iterative coupled methods was confirmed.Subsequently,multi-fidelity simulations of the speed and altitude characteristics were also conducted.The auxiliary and traditional iterative coupled methods were evaluated in terms of convergence speed and accuracy.The evaluation indicated that the auxiliary iterative coupled method significantly reduces iterative steps by approximately 50%at the near-choked state.In general,the auxiliary iterative coupled method is preferred as a development of the traditional iterative coupled method in the near-choked state,and the combined auxiliary-traditional iterative coupled method provides support for successful multi-fidelity simulation in far-off-design conditions.展开更多
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 experi...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.展开更多
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 ...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.展开更多
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 comp...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.展开更多
The big-tapered profiled ring disk is a key component of engines for rockets and missiles.A new forming technology,as called spinning-rolling process,has been proposed previously for the high performance,high efficien...The big-tapered profiled ring disk is a key component of engines for rockets and missiles.A new forming technology,as called spinning-rolling process,has been proposed previously for the high performance,high efficiency and low-cost manufacturing of the component.Blank design is the key part of plastic forming process design.For spinning-rolling process,the shape and size of the blank play a crucial role in process stability,deformation behavior and dimensional accuracy.So this work proposes a blank design method to determine the geometry structure and sizes of the blank.The mathematical model for calculating the blank size has been deduced based on volume conservation and neutral layer length invariance principle.The FE simulation and corresponding trial production of an actual big-tapered profiled ring disk show that the proposed blank design method is applicative.In order to obtain a preferred blank,the influence rules of blank size determined by different deformation degrees(rolling ratio k)on the spinning-rolling process are revealed by comprehensive FE simulations.Overall considering the process stability,circularity of the deformed ring disk and forming forces,a reasonable range of deformation degree(rolling ratio k)is recommended for the blank design of the new spinning-rolling process.展开更多
This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydra...This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.展开更多
Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulati...Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.展开更多
Computer simulation permits answering theoretical and applied questions in animal and plant breeding.Blib is a novel multi-module simulation platform,which is able to handle more complicated genetic effects and models...Computer simulation permits answering theoretical and applied questions in animal and plant breeding.Blib is a novel multi-module simulation platform,which is able to handle more complicated genetic effects and models than most existing tools.In this study,we describe one major and unified application module of Blib,i.e.,ISB(abbreviated from in silico breeding),for simulating the three categories of breeding programs for developing clonal,pure-line and hybrid cultivars in plants.Genetic models on environments and breeding-targeted traits,one or several parental populations,and a number of breeding methods are key elements to run simulation experiments in ISB,which are arranged in three external input files by given formats.Applications of ISB are illustrated by three case studies,representing the three categories of plant breeding programs.Under the condition that 5000 F1 progenies were generated and tested from 50 heterozygous parents,Case study I showed that 50 crosses,each of 100 progenies,made the best balance between genetic achievement and field cost.In Case study II,one optimum breeding method was identified by which the pure lines with high yield and medium maturity could be developed.Case study III investigated the genetic consequence in hybrid breeding from five testers.One tester was identified for the simultaneous improvement in F1 hybrids and inbred lines.In summary,ISB identified a balanced crossing scheme,an optimum pure-line selection method,and an optimized tester in three case studies which are relevant to plant breeding.We believe the prediction by simulation would be highly required in front of the next generation of breeding to be driven by informatics and intelligence.展开更多
Without knowing the emittance value, it is difficult to optimize ion beam optics for minimum beam loss during transmission, especially considering the very high emittance values of electron cyclotron resonance(ECR) io...Without knowing the emittance value, it is difficult to optimize ion beam optics for minimum beam loss during transmission, especially considering the very high emittance values of electron cyclotron resonance(ECR) ion sources.With this in mind, to measure the emittance of the ion beams produced by the mVINIS ECR, which is part of the FAMA facility at the Vin?a Institute of Nuclear Sciences, we have developed a pepper-pot scintillator screen system combined with a CMOS camera. The application, developed on the Lab VIEW platform, allows us to control the camera's main attribute settings, such as the shutter speed and the gain, record the images in the region of interest, and process and filter the images in real time. To analyze the data from the obtained image, we have developed an algorithm called measurement and analysis of ion beam luminosity(MAIBL) to reconstruct the four-dimensional(4D) beam profile and calculate the root mean square(RMS) emittance. Before measuring emittance, we performed a simulated experiment using the pepper-pot simulation(PPS) program. An exported file(PPS) gives a numerically generated raw image(mock image) of a beam with a predefined emittance value after it has passed through a pepper-pot mask. By analyzing data from mock images instead of the image obtained by the camera and putting it into the MAIBL algorithm, we can compare the calculated emittance with PPS's initial emittance value. In this paper, we present our computational tools and explain the method for verifying the correctness of the calculated emittance values.展开更多
This study sheds light on how pore structure characteristics and varying dynamic pressure conditions influence the permeability of tight sandstone reservoirs,with a particular focus on the Paleozoic reservoirs in the ...This study sheds light on how pore structure characteristics and varying dynamic pressure conditions influence the permeability of tight sandstone reservoirs,with a particular focus on the Paleozoic reservoirs in the Qingshimao Gas Field.Using CT scans of natural core samples,a three-dimensional digital core was constructed.The maximum ball method was applied to extract a related pore network model,and the pore structure characteristics of the core samples,such as pore radius,throat radius,pore volume,and coordination number,were quantitatively evaluated.The analysis revealed a normally distributed pore radius,suggesting a high degree of reservoir homogeneity and favorable conditions for a connected pore system.However,it was found that the majority of throat radii measured less than 1μm,which severely restricted fluid flow and diminished permeability.Over 50%of the pores measured under 100μm^(3),further constraining fluid movement.Additionally,30%-50%of the pore network was composed of isolated and blind-end pores,which significantly impaired formation connectivity and reduced permeability.Based on this,the lattice Boltzmann method(LBM)was used for pore-scale flow simulation to investigate the influence mechanism of pore structure characteristics and dynamic-static parameters such as displacement pressure difference on the permeability performance of the considered tight sandstone reservoirs for various pressure gradients(0.1,1,and 10 MPa).The simulations revealed a strong relationship between pressure differential and both the number of streamlines and flow path tortuosity.When the pressure differential increased to 1 MPa,30 streamlines were observed,with a tortuosity factor of 1.5,indicating the opening of additional seepage channels and the creation of increasingly winding flow paths.展开更多
In response to the issue of determining the appropriate capacity when hybrid energy storage systems(HESS)collaborate with thermal power units(TPU)in the system’s secondary frequency regulation,a configuration method ...In response to the issue of determining the appropriate capacity when hybrid energy storage systems(HESS)collaborate with thermal power units(TPU)in the system’s secondary frequency regulation,a configuration method for HESS based on the analysis of frequency regulation demand analysis is proposed.And a corresponding simulation platform is developed.Firstly,a frequency modulation demand method for reducing the frequency modulation losses of TPU is proposed.Secondly,taking into comprehensive consideration that flywheel energy storage features rapid power response and battery energy storage has the characteristic of high energy density,a coordinated control strategy for HESS considering the self-recovery of state of charge(SOC)is put forward.Then,to measure the economic and technical performance of HESS in assisting the secondary frequency modulation of TPU,an optimized configurationmodel considering the full-life-cycle economy and frequency modulation performance of TPU and HESS system is constructed.Finally,a visual simulation platform for the combined frequency modulation of TPU and HESS is developed based on Matlab Appdesigner.Theresults of calculation examples indicate that the proposed configuration method can improve the overall economic efficiency and frequency modulation performance of TPU and HESS;The control strategy can not only prolong the service life of battery energy storage but also enhance the continuous response ability of HESS;The visual simulation platform is easy to use,and the simulation results are accurate and reliable.展开更多
The traditional detailed model of the dual active bridge(DAB)power electronic transformer is characterized by the high dimensionality of its nodal admittance matrix and the need for a small simulation step size,which ...The traditional detailed model of the dual active bridge(DAB)power electronic transformer is characterized by the high dimensionality of its nodal admittance matrix and the need for a small simulation step size,which limits the speed of electromagnetic transient(EMT)simulations.To overcome these limitations,a novel EMT equivalent model based on a generalized branch-cutting method is proposed to improve the simulation efficiency of the DAB model.The DAB topology is first decomposed into two subnetworks through branch-cutting and node-tearing methods without the introduction of a one-time-step delay.Sub-sequently,the internal nodes of each sub-network are eliminated through network simplification,and the equivalent circuit for the port cascade module is derived.The model is then validated through simulations across various operating conditions.The results demonstrate that the model avoids the loss of accuracy associated with one-time-step delay,the relative error across different conditions remains below 1%,and the simulation acceleration ratios improve as the number of modules increases.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金supported by the National Key Research and Development Program of China(2022YFC2904400)Guangxi Science and Technology Major Project(Gui Ke AA23023033)。
文摘As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.
基金funded by the 14th Five-Year Plan Major Science and Technology Project of CNOOC project number KJGG2021-0506.
文摘After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.
基金This work was supported by the National Key R&D Program of China(Grant No.2022YFE0207000)the National Natural Science Foundation of China(Grant Nos.12372289,11972250,and 12102298)+1 种基金the China Postdoctoral Science Foundation(Grant No.2021M702443)Tianjin Natural Science Foundation(Grant No.22JCZDJC00910).
文摘At low-Reynolds-number,the performance of airfoil is known to be greatly affected by the formation and burst of a laminar separation bubble(LSB),which requires a more precise simulation of the delicate flow structures.A framework based on the interior penalty discontinuous Galerkin method and large eddy simulation approach was adopted in the present study.The performances of various subgrid models,including the Smagorinsky(SM)model,the dynamic Smagorinsky(DSM)model,the wall-adapting local-eddy-viscosity(WALE)model,and the VREMAN model,have been analyzed through flow simulations of the SD7003 airfoil at a Reynolds number of 60000.It turns out that the SM model fails to predict the emergence of LSB,even modified by the Van-Driest damping function.On the contrary,the best agreement is generally achieved by the WALE model in terms of flow separation,reattachment,and transition locations,together with the aerodynamic loads.Furthermore,the influence of numerical dissipation has also been discussed through the comparison of skin friction and resolved Reynolds stresses.As numerical dissipation decreases,the prediction accuracy of the WALE model degrades.Meanwhile,nonlinear variation could be observed from the performances of the DSM model,which could be attributed to the interaction between the numerical dissipation and the subgrid model.
基金the financial support from National Key R&D Program of China(Grant number:2024YFC2815100)Natural Science Foundation of China(Grant number:52322110)Beijing Nova Program(Grant number:20230484341).
文摘Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.
基金National Key R&D Program of China(No.2017YFB1304000)Fundamental Research Funds for the Central Universities,China(No.2232023G-05-1)。
文摘The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.
基金National Key R&D Program of China under Grant No.2022YFC3003601。
文摘The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improved stochastic finite-fault method to systematically assess seismic impacts.Observed near-field recordings at MM.NGU station was used to determine the reliability of the theoretically derived stress drop as input for simulation.Far-field recordings constrained the frequency-dependent S-wave quality factors(Q(f)=283.305f^(0.588))for anelastic attenuation modeling.Comparisons of peak accelerations between simulation and empirical ground-motion models showed good agreement at moderate-to-large distances.However,lower near-fault simulations indicate a weaker-than-average source effect.Analysis of simulated instrumental seismic intensity revealed key patterns.Maximum intensity(Ⅹ)occurred in isolated patches within the ruptured fault projection,correlating with shallow high-slip areas.TheⅨ-intensity zone formed a north-south elongated band centered on fault projection.Significant asymmetry inⅧ-intensity distribution perpendicular to the fault strike was observed,with a wider western extension attributed to lower shear-wave velocities west of the fault.Supershear rupture behavior enhanced ground motions,expanding intensity ranges by~20%compared to sub-shear rupture.This study reveals the integrated effects of fault geometry,slip spatial distribution,rupture velocity,and site condition in governing ground motion patterns.
基金supported by the Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area(China Three Gorges University),Ministry of Education(No.2022KDZ03)the Science and Technology Projects of Yunnan Provincial Science and Technology Department(No.202401AT070328)+1 种基金the Young talents project of“Xingdian Talent Support Program”in Yunnan Province(No.YNWR-QNBJ-2020-019)the Fund Project of China Academy of Railway Sciences Co.,Ltd.(No.2021YJ178)。
文摘0 INTRODUCTION In recent years,modern railways have been actively under construction in the complex mountainous area of Southwest China.However,rockfall poses a significant threat to both construction and operation phases of railway projects(Yan et al.,2023;Chen et al.,2022;Fanos and Pradhan,2018).
基金funded by the Science and Technology Innovation Committee Foundation of Shenzhen,China(Nos.JCYJ20200109141403840 and ZDSYS20220527171405012)the National Natural Science Foundation of China(No.52106045)the Pearl River Talent Recruitment Program,China(No.2019CX01Z084)。
文摘Iterative coupled methods are widely used in multi-fidelity simulation of rotating components due to the simple implementation,which iteratively eliminates the errors between the computational fluid dynamics models and approximate characteristic maps.However,the convergence and accuracy of the iterative coupled method are trapped in characteristic maps.In particular,iterative steps increase sharply as the operation point moves away from the design point.To address these problems,this paper developed an auxiliary iterative coupled method that introduces the static-pressure-auxiliary characteristic maps and modification factor of mass flow into the component-level model.The developed auxiliary method realized the direct transfer of static pressure between the high-fidelity models and the component-level model.Multi-fidelity simulations of the throttle characteristics were carried out using both the auxiliary and traditional iterative coupled methods,and the simulation results were verified using the experimental data.Additionally,the consistency between the auxiliary and traditional iterative coupled methods was confirmed.Subsequently,multi-fidelity simulations of the speed and altitude characteristics were also conducted.The auxiliary and traditional iterative coupled methods were evaluated in terms of convergence speed and accuracy.The evaluation indicated that the auxiliary iterative coupled method significantly reduces iterative steps by approximately 50%at the near-choked state.In general,the auxiliary iterative coupled method is preferred as a development of the traditional iterative coupled method in the near-choked state,and the combined auxiliary-traditional iterative coupled method provides support for successful multi-fidelity simulation in far-off-design conditions.
基金supported by National Natural Science Foundation of China(Nos.92066108 and 51277061)。
文摘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.
基金the National Natural Science Foundation of China(No.51875062,No.52205336)the China Postdoctoral Science Foundation(No.2021M700567).
文摘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.
文摘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.
基金the National Natural Science Foundation of China(No.52275378)the National Key Laboratory for Precision Hot Processing of Metals(6142909200208)。
文摘The big-tapered profiled ring disk is a key component of engines for rockets and missiles.A new forming technology,as called spinning-rolling process,has been proposed previously for the high performance,high efficiency and low-cost manufacturing of the component.Blank design is the key part of plastic forming process design.For spinning-rolling process,the shape and size of the blank play a crucial role in process stability,deformation behavior and dimensional accuracy.So this work proposes a blank design method to determine the geometry structure and sizes of the blank.The mathematical model for calculating the blank size has been deduced based on volume conservation and neutral layer length invariance principle.The FE simulation and corresponding trial production of an actual big-tapered profiled ring disk show that the proposed blank design method is applicative.In order to obtain a preferred blank,the influence rules of blank size determined by different deformation degrees(rolling ratio k)on the spinning-rolling process are revealed by comprehensive FE simulations.Overall considering the process stability,circularity of the deformed ring disk and forming forces,a reasonable range of deformation degree(rolling ratio k)is recommended for the blank design of the new spinning-rolling process.
基金Funded by the Natural Science Foundation of Jiangsu Province(No.BK20241529)China Postdoctoral Science Foundation(No.2024M750736)。
文摘This study aims to develop a chloride diffusion simulation method that considers the hydration microstructure and pore solution properties during the hydration of tricalcium silicate(C3S).The method combines the hydration simulation,thermodynamic calculation,and finite element analysis to examine the effects of pore solution,including effect of electrochemical potential,effect of chemical activity,and effect of mechanical interactions between ions,on the chloride effective diffusion coefficient of hydrated C3S paste.The results indicate that the effect of electrochemical potential on chloride diffusion becomes stronger with increasing hydration age due to the increase in the content of hydrated calcium silicate;as the hydration age increases,the effect of chemical activity on chloride diffusion weakens when the number of diffusible elements decreases;the effect of mechanical interactions between ions on chloride diffusion decreases with the increase of hydration age.
基金financially supported by the Riset Pengembangan dan Penerapan(RPP),Diponegoro University 2023 research scheme with contract number 609-18/UN7.D2/PP/VIII/2023.
文摘Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.
基金supported by Biological Breeding-National Science and Technology Major Project(2023ZD0407501)National Natural Science Foundation of China(31861143003)Innovation Program of Chinese Academy of Agricultural Sciences.
文摘Computer simulation permits answering theoretical and applied questions in animal and plant breeding.Blib is a novel multi-module simulation platform,which is able to handle more complicated genetic effects and models than most existing tools.In this study,we describe one major and unified application module of Blib,i.e.,ISB(abbreviated from in silico breeding),for simulating the three categories of breeding programs for developing clonal,pure-line and hybrid cultivars in plants.Genetic models on environments and breeding-targeted traits,one or several parental populations,and a number of breeding methods are key elements to run simulation experiments in ISB,which are arranged in three external input files by given formats.Applications of ISB are illustrated by three case studies,representing the three categories of plant breeding programs.Under the condition that 5000 F1 progenies were generated and tested from 50 heterozygous parents,Case study I showed that 50 crosses,each of 100 progenies,made the best balance between genetic achievement and field cost.In Case study II,one optimum breeding method was identified by which the pure lines with high yield and medium maturity could be developed.Case study III investigated the genetic consequence in hybrid breeding from five testers.One tester was identified for the simultaneous improvement in F1 hybrids and inbred lines.In summary,ISB identified a balanced crossing scheme,an optimum pure-line selection method,and an optimized tester in three case studies which are relevant to plant breeding.We believe the prediction by simulation would be highly required in front of the next generation of breeding to be driven by informatics and intelligence.
基金funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia via the FAMA project (research topic “Physics and chemistry with ion beams”)。
文摘Without knowing the emittance value, it is difficult to optimize ion beam optics for minimum beam loss during transmission, especially considering the very high emittance values of electron cyclotron resonance(ECR) ion sources.With this in mind, to measure the emittance of the ion beams produced by the mVINIS ECR, which is part of the FAMA facility at the Vin?a Institute of Nuclear Sciences, we have developed a pepper-pot scintillator screen system combined with a CMOS camera. The application, developed on the Lab VIEW platform, allows us to control the camera's main attribute settings, such as the shutter speed and the gain, record the images in the region of interest, and process and filter the images in real time. To analyze the data from the obtained image, we have developed an algorithm called measurement and analysis of ion beam luminosity(MAIBL) to reconstruct the four-dimensional(4D) beam profile and calculate the root mean square(RMS) emittance. Before measuring emittance, we performed a simulated experiment using the pepper-pot simulation(PPS) program. An exported file(PPS) gives a numerically generated raw image(mock image) of a beam with a predefined emittance value after it has passed through a pepper-pot mask. By analyzing data from mock images instead of the image obtained by the camera and putting it into the MAIBL algorithm, we can compare the calculated emittance with PPS's initial emittance value. In this paper, we present our computational tools and explain the method for verifying the correctness of the calculated emittance values.
文摘This study sheds light on how pore structure characteristics and varying dynamic pressure conditions influence the permeability of tight sandstone reservoirs,with a particular focus on the Paleozoic reservoirs in the Qingshimao Gas Field.Using CT scans of natural core samples,a three-dimensional digital core was constructed.The maximum ball method was applied to extract a related pore network model,and the pore structure characteristics of the core samples,such as pore radius,throat radius,pore volume,and coordination number,were quantitatively evaluated.The analysis revealed a normally distributed pore radius,suggesting a high degree of reservoir homogeneity and favorable conditions for a connected pore system.However,it was found that the majority of throat radii measured less than 1μm,which severely restricted fluid flow and diminished permeability.Over 50%of the pores measured under 100μm^(3),further constraining fluid movement.Additionally,30%-50%of the pore network was composed of isolated and blind-end pores,which significantly impaired formation connectivity and reduced permeability.Based on this,the lattice Boltzmann method(LBM)was used for pore-scale flow simulation to investigate the influence mechanism of pore structure characteristics and dynamic-static parameters such as displacement pressure difference on the permeability performance of the considered tight sandstone reservoirs for various pressure gradients(0.1,1,and 10 MPa).The simulations revealed a strong relationship between pressure differential and both the number of streamlines and flow path tortuosity.When the pressure differential increased to 1 MPa,30 streamlines were observed,with a tortuosity factor of 1.5,indicating the opening of additional seepage channels and the creation of increasingly winding flow paths.
基金supported by a Key Project of the National Natural Science Foundation of China under Grant 52337004.
文摘In response to the issue of determining the appropriate capacity when hybrid energy storage systems(HESS)collaborate with thermal power units(TPU)in the system’s secondary frequency regulation,a configuration method for HESS based on the analysis of frequency regulation demand analysis is proposed.And a corresponding simulation platform is developed.Firstly,a frequency modulation demand method for reducing the frequency modulation losses of TPU is proposed.Secondly,taking into comprehensive consideration that flywheel energy storage features rapid power response and battery energy storage has the characteristic of high energy density,a coordinated control strategy for HESS considering the self-recovery of state of charge(SOC)is put forward.Then,to measure the economic and technical performance of HESS in assisting the secondary frequency modulation of TPU,an optimized configurationmodel considering the full-life-cycle economy and frequency modulation performance of TPU and HESS system is constructed.Finally,a visual simulation platform for the combined frequency modulation of TPU and HESS is developed based on Matlab Appdesigner.Theresults of calculation examples indicate that the proposed configuration method can improve the overall economic efficiency and frequency modulation performance of TPU and HESS;The control strategy can not only prolong the service life of battery energy storage but also enhance the continuous response ability of HESS;The visual simulation platform is easy to use,and the simulation results are accurate and reliable.
基金The Technology Project of State Grid Corporation of China Headquarters(No.5400-202318547A-3-2-ZN).
文摘The traditional detailed model of the dual active bridge(DAB)power electronic transformer is characterized by the high dimensionality of its nodal admittance matrix and the need for a small simulation step size,which limits the speed of electromagnetic transient(EMT)simulations.To overcome these limitations,a novel EMT equivalent model based on a generalized branch-cutting method is proposed to improve the simulation efficiency of the DAB model.The DAB topology is first decomposed into two subnetworks through branch-cutting and node-tearing methods without the introduction of a one-time-step delay.Sub-sequently,the internal nodes of each sub-network are eliminated through network simplification,and the equivalent circuit for the port cascade module is derived.The model is then validated through simulations across various operating conditions.The results demonstrate that the model avoids the loss of accuracy associated with one-time-step delay,the relative error across different conditions remains below 1%,and the simulation acceleration ratios improve as the number of modules increases.
