An Innovative Coupled Common-Node Discrete Element Method-Smoothed Particle Hydrodynamics Model Developed with LS-DYNA and Its Applications

In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.

A modified smoothed particle hydrodynamics method considering residual stress for simulating failure and its application in layered rock mass

Residual strength is an indispensable factor in evaluating rock fracture,yet the current Smoothed Particle Hydrodynamics(SPH)framework rarely considers its influence when simulating fracture.An improved cracking strategy considering residual stress in the base bond SPH method was proposed to simulate failures in layered rocks and slopes and verified by experimental results and other simulation methods(i.e.,the discrete element method).Modified Mohr–Coulomb failure criterion was applied to distinguish the mixed failure of tensile and shear.Bond fracture markψwas introduced to improve the kernel function after tensile damage,and the calculation of residual stress after the damage was derived after shear damage.Numerical simulations were carried out to evaluate its performance under different stress and scale conditions and to verify its effectiveness in realistically reproducing crack initiation and propagation and coalescence,even fracture and separation.The results indicate that the improved cracking strategy precisely captures the fracture and failure pattern in layered rocks and rock slopes.The residual stress of brittle tock is correctly captured by the improved SPH method.The improved SPH method that considers residual strength shows an approximately 13%improvement in accuracy for the safety factor of anti-dip layered slopes compared to the method that does not consider residual strength,as validated against analytical solutions.We infer that the improved SPH method is effective and shows promise for applications to continuous and discontinuous rock masses.

A multi-scale second-order autoregressive recursive filter approach for the sea ice concentration analysis

To effectively extract multi-scale information from observation data and improve computational efficiency,a multi-scale second-order autoregressive recursive filter(MSRF)method is designed.The second-order autoregressive filter used in this study has been attempted to replace the traditional first-order recursive filter used in spatial multi-scale recursive filter(SMRF)method.The experimental results indicate that the MSRF scheme successfully extracts various scale information resolved by observations.Moreover,compared with the SMRF scheme,the MSRF scheme improves computational accuracy and efficiency to some extent.The MSRF scheme can not only propagate to a longer distance without the attenuation of innovation,but also reduce the mean absolute deviation between the reconstructed sea ice concentration results and observations reduced by about 3.2%compared to the SMRF scheme.On the other hand,compared with traditional first-order recursive filters using in the SMRF scheme that multiple filters are executed,the MSRF scheme only needs to perform two filter processes in one iteration,greatly improving filtering efficiency.In the two-dimensional experiment of sea ice concentration,the calculation time of the MSRF scheme is only 1/7 of that of SMRF scheme.This means that the MSRF scheme can achieve better performance with less computational cost,which is of great significance for further application in real-time ocean or sea ice data assimilation systems in the future.

An Updated Lagrangian Particle Hydrodynamics (ULPH)-NOSBPD Coupling Approach forModeling Fluid-Structure Interaction Problem

A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.

Smoothed-Particle Hydrodynamics Simulation of Ship Motion and Tank Sloshing under the Effect of Regular Waves

Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to ship motion and liquid tank sloshing has been simulated by the Smoothed-Particle Hydrodynamics(SPH)method.Firstly,the sloshing flow in a rectangular tank was simulated and the related loads were analyzed to verify and validate the accuracy of the present SPH solver.Then,a three-dimensional simplified LNG carrier model,including two prismatic liquid tanks and a wave tank,was introduced.Different conditions were examined corresponding to different wave lengths,wave heights,wave heading angles,and tank loading rates.Finally,the effects of liquid tank loading rate on LNG ship motions and sloshing loading were analyzed,thereby showing that the SPH method can effectively provide useful indications for the design of liquid cargo ships.

Updated Lagrangian Particle Hydrodynamics (ULPH)Modeling of Natural Convection Problems

Natural convection is a heat transfer mechanism driven by temperature or density differences,leading to fluid motion without external influence.It occurs in various natural and engineering phenomena,influencing heat transfer,climate,and fluid mixing in industrial processes.This work aims to use the Updated Lagrangian Particle Hydrodynamics(ULPH)theory to address natural convection problems.The Navier-Stokes equation is discretized using second-order nonlocal differential operators,allowing a direct solution of the Laplace operator for temperature in the energy equation.Various numerical simulations,including cases such as natural convection in square cavities and two concentric cylinders,were conducted to validate the reliability of the model.The results demonstrate that the proposed model exhibits excellent accuracy and performance,providing a promising and effective numerical approach for natural convection problems.

On the Vacuum Hydrodynamics of Moving Bodies—The Theory of General Singularity

The Theory of General Singularity is presented, unifying quantum field theory, general relativity, and the standard model. This theory posits phonons as fundamental excitations in a quantum vacuum, modeled as a Bose-Einstein condensate. Through key equations, the role of phonons as intermediaries between matter, energy, and spacetime geometry is demonstrated. The theory expands Einsteins field equations to differentiate between visible and dark matter, and revises the standard model by incorporating phonons. It addresses dark matter, dark energy, gravity, and phase transitions, while making testable predictions. The theory proposes that singularities, the essence of particles and black holes, are quantum entities ubiquitous in nature, constituting the very essence of elementary particles, seen as micro black holes or quantum fractal structures of spacetime. As the theory is refined with increasing mathematical rigor, it builds upon the foundation of initial physical intuition, connecting the spacetime continuum of general relativity with the hydrodynamics of the quantum vacuum. Inspired by the insights of Tesla and Majorana, who believed that physical intuition justifies the infringement of mathematical rigor in the early stages of theory development, this work aims to advance the understanding of the fundamental laws of the universe and the perception of reality.

The Convergence Rate of Fréchet Distribution under the Second-Order Regular Variation Condition

In this article we consider the asymptotic behavior of extreme distribution with the extreme value index γ>0 . The rates of uniform convergence for Fréchet distribution are constructed under the second-order regular variation condition.

Closed-form steady-state solutions for forced vibration of second-order axially moving systems

Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for the forced vibration when the damping effect and the coupling effect of multiple second-order models are considered.In this paper, Green's function method based on the Laplace transform is used to obtain closed-form solutions for the forced vibration of second-order axially moving systems. By taking the axially moving damping string system and multi-string system connected by springs as examples, the detailed solution methods and the analytical Green's functions of these second-order systems are given. The mode functions and frequency equations are also obtained by the obtained Green's functions. The reliability and convenience of the results are verified by several examples. This paper provides a systematic analytical method for the dynamic analysis of second-order axially moving systems, and the obtained Green's functions are applicable to different second-order systems rather than just string systems. In addition, the work of this paper also has positive significance for the study on the forced vibration of high-order systems.

Application of Smoothed Particle Hydrodynamics(SPH)for the Simulation of Flow-Like Landslides on 3D Terrains

Flow-type landslide is one type of landslide that generally exhibits characteristics of high flow velocities,long jump distances,and poor predictability.Simulation of its propagation process can provide solutions for risk assessment and mitigation design.The smoothed particle hydrodynamics(SPH)method has been successfully applied to the simulation of two-dimensional(2D)and three-dimensional(3D)flow-like landslides.However,the influence of boundary resistance on the whole process of landslide failure is rarely discussed.In this study,a boundary condition considering friction is proposed and integrated into the SPH method,and its accuracy is verified.Moreover,the Navier-Stokes equation combined with the non-Newtonian fluid rheologymodel was utilized to solve the dynamic behavior of the flow-like landslide.To verify its performance,the Shuicheng landslide event,which occurred in Guizhou,China,was taken as a case study.In the 2D simulation,a sensitivity analysis was conducted,and the results showed that the shearing strength parameters have more influence on the computation accuracy than the coefficient of viscosity.Afterwards,the dynamic characteristics of the landslide,such as the velocity and the impact area,were analyzed in the 3D simulation.The simulation results are in good agreement with the field investigations.The simulation results demonstrate that the SPH method performs well in reproducing the landslide process,and facilitates the analysis of landslide characteristics as well as the affected areas,which provides a scientific basis for conducting the risk assessment and disaster mitigation design.

Numerical and experimental study on the particle erosion and gas–particle hydrodynamics in an integral multi-jet swirling spout-fluidized bed

In this paper,using the computational fluid dynamics based on Euler Lagrange and the commercial software Barracuda VR,the gas-particle hydrodynamics and the erosion of particles on the inner wall and internal components of the spouted bed in the integrated multi-jet swirling spout-fluidized bed(IMSSFB)are studied.Erosion experiments have obtained the characterization of particle erosion on internal components and verified the relevant numerical models.The results show that:the particle distribution within the IMSSFB is uneven due to the cyclonic effect of the axial swirl vane(ASV),resulting in particle erosion for the ASV being concentrated on one side;when the gas reaches the top,too high an erosion gas velocity leads to gas backflow.As the filling height increases,there is a tendency for the erosion position of the particles on the ASV to expand upwards.However,the effect of increasing gas velocity on the erosion position is insignificant.

Three-dimensional mixed convection stagnation-point fow past a vertical surface with second-order slip velocity

This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.

Six New Sets of the Non-Elementary Jef-Family-Functions that Are Giving Solutions to Some Second-Order Nonlinear Autonomous ODEs

In this paper, we define some new sets of non-elementary functions in a group of solutions x(t) that are sine and cosine to the upper limit of integration in a non-elementary integral that can be arbitrary. We are using Abel’s methods, described by Armitage and Eberlein. The key is to start with a non-elementary integral function, differentiating and inverting, and then define a set of three functions that belong together. Differentiating these functions twice gives second-order nonlinear ODEs that have the defined set of functions as solutions. We will study some of the second-order nonlinear ODEs, especially those that exhibit limit cycles. Using the methods described in this paper, it is possible to define many other sets of non-elementary functions that are giving solutions to some second-order nonlinear autonomous ODEs.

Collisionless shock acceleration of protons in a plasma slab produced in a gas jet by the collision of two laser-driven hydrodynamic shockwaves

We have recently proposed a new technique of plasma tailoring by laser-driven hydrodynamic shockwaves generated on both sides of a gas jet[Marquès et al.,Phys.Plasmas 28,023103(2021)].In a continuation of this numerical work,we study experimentally the influence of the tailoring on proton acceleration driven by a high-intensity picosecond laser in three cases:without tailoring,by tailoring only the entrance side of the picosecond laser,and by tailoring both sides of the gas jet.Without tailoring,the acceleration is transverse to the laser axis,with a low-energy exponential spectrum,produced by Coulomb explosion.When the front side of the gas jet is tailored,a forward acceleration appears,which is significantly enhanced when both the front and back sides of the plasma are tailored.This forward acceleration produces higher-energy protons,with a peaked spectrum,and is in good agreement with the mechanism of collisionless shock acceleration(CSA).The spatiotemporal evolution of the plasma profile is characterized by optical shadowgraphy of a probe beam.The refraction and absorption of this beam are simulated by post-processing 3D hydrodynamic simulations of the plasma tailoring.Comparison with the experimental results allows estimation of the thickness and near-critical density of the plasma slab produced by tailoring both sides of the gas jet.These parameters are in good agreement with those required for CSA.

Hydrodynamic Cavitation Enhanced SR-Aops Degradation of Organic Pollutants in Water:A Review

SR-AOP(sulfate radical advanced oxidation process)is a novel water treatment method able to eliminate refractory organic pollutants.Hydrodynamic cavitation(HC)is a novel green technology,that can effectively produce strong oxidizing sulfate radicals.This paper presents a comprehensive review of the research advancements in these fields and a critical discussion of the principal factors influencing HC-enhanced SR-AOP and the mechanisms of synergistic degradation.Furthermore,some insights into the industrial application of HC/PS are also provided.Current research shows that this technology is feasible at the laboratory stage,but its application on larger scales requires further understanding and exploration.In this review,some attention is also paid to the design of the hydrodynamic cavitation reactor and the related operating parameters.

Hydrodynamic resistance of pore–throat structures and its effect on shale oil apparent permeability

Oil transport is greatly affected by heterogeneous pore–throat structures present in shale.It is therefore very important to accurately characterize pore–throat structures.Additionally,it remains unclear how pore–throat structures affect oil transport capacity.In this paper,using finite element(FE)simulation and mathematical modeling,we calculated the hydrodynamic resistance for four pore–throat structure.In addition,the influence of pore throat structure on shale oil permeability is analyzed.According to the results,the hydrodynamic resistance of different pore throat structures can vary by 300%.The contribution of additional resistance caused by streamline bending is also in excess of 40%,even without slip length.Fur-thermore,Pore–throat structures can affect apparent permeability by more than 60%on the REV scale,and this influence increases with heterogeneity of pore size distribution,organic matter content,and organic matter number.Clearly,modeling shale oil flow requires consideration of porous–throat structure and additional resistance,otherwise oil recovery and flow capacity may be overestimated.

Simulating the Escaping Atmosphere of GJ 436 b with Two-fluid Magnetohydrodynamic Models

Observations of transmission spectra reveal that hot Jupiters and Neptunes are likely to possess escaping atmospheres driven by stellar radiation.Numerous models predict that magnetic fields may exert significant influences on the atmospheres of hot planets.Generally,the escaping atmospheres are not entirely ionized,and magnetic fields only directly affect the escape of ionized components within them.Considering the chemical reactions between ionized components and neutral atoms,as well as collision processes,magnetic fields indirectly impact the escape of neutral atoms,thereby influencing the detection signals of planetary atmospheres in transmission spectra.In order to simulate this process,we developed a magnetohydrodynamic multi-fluid model based on MHD code PLUTO.As an initial exploration,we investigated the impact of magnetic fields on the decoupling of H^(+)and H in the escaping atmosphere of the hot Neptune GJ436b.Due to the strong resonant interactions between H and H^(+),the coupling between them is tight even if the magnetic field is strong.Of course,alternatively,our work also suggests that merging H and H^(+)into a single flow can be a reasonable assumption in MHD simulations of escaping atmospheres.However,our simulation results indicate that under the influence of magnetic fields,there are noticeable regional differences in the decoupling of H^(+)and H.With the increase of magnetic field strength,the degree of decoupling also increases.For heavier particles such as O,the decoupling between O and H^(+)is more pronounced.Our findings provide important insights for future studies on the decoupling processes of heavy atoms in the escaping atmospheres of hot Jupiters and hot Neptunes under the influence of magnetic fields.

Second-Order MaxEnt Predictive Modelling Methodology. I: Deterministically Incorporated Computational Model (2nd-BERRU-PMD)

This work presents a comprehensive second-order predictive modeling (PM) methodology designated by the acronym 2nd-BERRU-PMD. The attribute “2nd” indicates that this methodology incorporates second-order uncertainties (means and covariances) and second-order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best- Estimate Results with Reduced Uncertainties” and the last letter (“D”) in the acronym indicates “deterministic,” referring to the deterministic inclusion of the computational model responses. The 2nd-BERRU-PMD methodology is fundamentally based on the maximum entropy (MaxEnt) principle. This principle is in contradistinction to the fundamental principle that underlies the extant data assimilation and/or adjustment procedures which minimize in a least-square sense a subjective user-defined functional which is meant to represent the discrepancies between measured and computed model responses. It is shown that the 2nd-BERRU-PMD methodology generalizes and extends current data assimilation and/or data adjustment procedures while overcoming the fundamental limitations of these procedures. In the accompanying work (Part II), the alternative framework for developing the “second- order MaxEnt predictive modelling methodology” is presented by incorporating probabilistically (as opposed to “deterministically”) the computed model responses.

Second-Order MaxEnt Predictive Modelling Methodology. II: Probabilistically Incorporated Computational Model (2nd-BERRU-PMP)

This work presents a comprehensive second-order predictive modeling (PM) methodology based on the maximum entropy (MaxEnt) principle for obtaining best-estimate mean values and correlations for model responses and parameters. This methodology is designated by the acronym 2nd-BERRU-PMP, where the attribute “2nd” indicates that this methodology incorporates second- order uncertainties (means and covariances) and second (and higher) order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best-Estimate Results with Reduced Uncertainties” and the last letter (“P”) in the acronym indicates “probabilistic,” referring to the MaxEnt probabilistic inclusion of the computational model responses. This is in contradistinction to the 2nd-BERRU-PMD methodology, which deterministically combines the computed model responses with the experimental information, as presented in the accompanying work (Part I). Although both the 2nd-BERRU-PMP and the 2nd-BERRU-PMD methodologies yield expressions that include second (and higher) order sensitivities of responses to model parameters, the respective expressions for the predicted responses, for the calibrated predicted parameters and for their predicted uncertainties (covariances), are not identical to each other. Nevertheless, the results predicted by both the 2nd-BERRU-PMP and the 2nd-BERRU-PMD methodologies encompass, as particular cases, the results produced by the extant data assimilation and data adjustment procedures, which rely on the minimization, in a least-square sense, of a user-defined functional meant to represent the discrepancies between measured and computed model responses.

Hydrodynamic characteristics and particle tracking of 90° lateral intakes at an inclined river slope

Lateral intakes are common in rivers.The pump effciency and sediment deposition are determined by the local hydrodynamic characteristics and mainstream division width.The hydraulic characteristics of lateral withdrawal from inclined river slopes at different intake elevations should be investigated.Meanwhile,the division width exhibits significant vertical non-uniformity at an inclined river slope,which should be clarified.Hence,a three-dimensional(3-D)hydrodynamic and particle-tracking model was developed with the Open Source Field Operation and Manipulation(Open FOAM),and the model was validated with physical model tests for 90°lateral withdrawal from an inclined side bank.The flow fields,withdrawal sources,and division widths were investigated with different intake bottom elevations,withdrawal discharges,and main channel velocities.This study showed that under inclined side bank conditions,water entered the intake at an oblique angle,causing significant 3-D spiral flows in the intake rather than two-dimensional closed recirculation.A lower withdrawal discharge,a lower bottom elevation of the intake,or a higher main channel velocity could further strengthen this phenomenon.The average division width and turbulent kinetic energy were smaller under inclined side bank conditions than under vertical bank conditions.With a low intake bottom elevation,a low withdrawal discharge,or a high main channel velocity,the sources of lateral withdrawal were in similar ranges near the local inclined bank in the vertical direction.Under inclined slope conditions,sediment deposition near the intake entrance could be reduced,compared to that under vertical slope conditions.The results provide hydrodynamic and sediment references for engineering designs for natural rivers with inclined terrains.