Improved Model for Soil as a Two-Phase Mixture Based on Smoothed Particle Hydrodynamics (SPH)

It is desired to resolve soil contamination with reduced costs. “Insoluble treatment” is a soil improvement method for heavy metal containing soil, which uses soil mixers to mix soil and soil improvement liquid agents. To reduce the costs of this method, soil mixers have to be optimized. However, it is not achieved due to the lack of theoretical knowledge on mixing solid with liquid. Therefore, a numerical model to simulate the dynamic behavior of solid and liquid is on the development in this study using Smoothed Particle Hydrodynamics (SPH) method. To validate the numerical model, several experiments were carried out and numerically reproduced. The comparisons of the results showed that the numerical model replicated a liquid flow with an error rate of 2.1% and a seepage flow with an error rate up to 26.1%. Especially, the water distribution in the soil pores was highly improved with absolute gaps in volumetric water content up to 4.4% in the porosity range of 10% - 90%. For the water absorption into dry sand, the simulation result became more realistic by concerning soil suction.

A Hydrodynamic Modei for Slug Frequency in Horizontal Gas-Liquid Two-Phase Flow

The prediction of slug frequency has important significance on gas-liquid two-phase flow. A hydrody-namic modei was put forward to evaluate slug frequency for horizontal two-phase flow, based on the dependence of slug frequency on the frequency of unstable interfacial wave. Using air and water, experimental verification of the modei was carried out in a large range of flow parameters. Six electrical probes were installed at different positions of a horizontal plexiglass pipe to detect slug frequency development. The pipe is 30m long and its inner diameter is 24 mm. It is observed experimentally that the interfacial wave frequency at the inlet is about l to 3 times the frequency of stable slug. The slug frequencies predicted by the modei fit well with Tronconi (1990) modei and the experimental data. The combination of the hydrodynamic modei and the experimental data results in a conclusion that the frequency of equilibrium liquid slug is approximately half the miniraum frequency of interfacial wave.

Numerical Simulations of Wave Impact Forces on the Open-Type Sea Access Road Using A Two-Phase SPH Model

A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analyses presented by Chen et al.(2022),in which the sea access roads are semi-immersed.In this new configuration,the sea access roads are placed above the still water level,therefore the presence of the air phase becomes a relevant issue in the determination of the wave forces acting on the structures.Indeed,the comparison of wave forces on the open-type sea access roads obtained from the single and two-phase SPH models with the experimental results shows that the latter are in much better agreement.So in the numerical simulations,a two-phaseδ-SPH model is adopted to investigate the dynamical problems.Based on the numerical results,the maximum horizontal and uplifting wave forces acting on the sea access roads are analyzed by considering different wave conditions and geometries of the structures.In particular,the presence of the girder is analyzed and the differences in the wave forces due to the air cushion effects which are created below the structure are highlighted.

IMPROVED HYDRODYNAMIC MODEL OF TWO-PHASE SLUG FLOW IN VERTICAL TUBES

Merits of the Fernandes model(Fernandes et al.1983)for two-phase slug flow in verticaltubes are reviewed in this paper.While predicting many macroscopic parameters of slug flow in verti-cal tubes,it fails to present correctly the trend that the average voidage in liquid slugs increases asthe rising velocity of Taylor bubbles is increased.It is also desirable to extend its application toelectrolyte systems, and to churn flow conditions.Based on the diagnostic analysis,the model equa-tion for gas entrainment by falling liquid film is reformulated and the influence of surface tension isalso accounted for.Development of the falling liquid film is recognized in the revised model in or-der to suit the case of short Taylor bubbles as well.The modified model predicts the variation of av-erage voidage in liquid slugs in good agreement with available experimental data.

Numerical Investigation of Land Reclamation Effects on Hydrodynamics and Mangroves in Shacheng Bay for the Last 36 Years

Since the 1980s,the robust economic growth of China has prompted extensive land reclamation projects along its coastline,notably affecting local hydrodynamics and resulting in ecological repercussions.Using a nearshore finite volume ocean model,we constructed a hydrodynamic model for Shacheng Bay,a southeastern coastal region with a winding and narrow entrance.We examined the hydrodynamic changes and mangrove dynamics over the past 36 years and the relationship between hydrodynamic alterations and mangrove degradation.Simulation results reveal that extensive reclamation projects between 1984 and 2000 weakened the current in Shacheng Bay,leading to decreased water exchange capacity and a significant reduction in mangrove area from 0.3 to 0.06km^(2).During this period,over 37% of mangrove degradation was ascribed to time-changing hydrodynamic variables without the direct influence of land reclamation.The results also highlight the changes in local hydrodynamics and water exchange patterns that adversely influenced mangrove growth.From 2000 to 2020,there were minimal coastline changes in Shacheng Bay,demonstrating reduced land reclamation activities.This stopped the further weakening trend of the currents,with a slight increase during ebb tides,while the residual current continued to weaken due to the decreasing tidal prism and water exchange capacity.The mangrove area partially recovered during this period,expanding from 0.06 to 0.11 km^(2),predominantly in new areas instead of where mangroves disappeared from 1984 to 2000.This work underlines the intricate relationship between land reclamation,hydrodynamics,and mangrove ecosystems,underscoring the need for sustainable coastal development strategies.

Evaluation of frictional pressure drop correlations for air-water and air-oil two-phase flow in pipeline-riser system

Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.

Effect of length-width ratio of rounded rectangle aquaculture tank in dual-diagonal-inlet layout on hydrodynamics

To adapt to the change of aquaculture workshop site,optimize the shape of aquaculture tanks and improve the utilization rate of breeding space,it is necessary to determine the appropriate length width ratio parameters of aquaculture tanks.In this paper,computational fluid dynamics(CFD)technology is adopted to study the flow field performance of aquaculture tanks with different L/B ratios(L:the length;B:the width,of aquaculture tank)and different jet direction conditions(lengthways jet and widthways jet).A three-dimensional numerical calculation model of turbulence in rounded rectangle aquaculture tanks in dual-diagonal-inlet layout was established.Jet directions are arranged lengthways and widthways,and the water flow velocity,resistance coefficient change,vorticity,etc.are analyzed under two working conditions.Results show that the flow field performance in aquaculture tank decreases with the increase of the L/B ratio.The flow field performed well when L/B was 1.0-1.3,sharply dropped at 1.4-1.6,and poor at 1.7-1.9.The results provided a theoretical basis for the design and optimization in flow field performance of the industrialized circulating aquaculture tanks.

Insight into evolution of invasive patterns on fingering phenomenon during immiscible two-phase flow through pore structure

Understanding fingering, as a challenge to stable displacement during the immiscible flow, has become a crucial phenomenon for geological carbon sequestration, enhanced oil recovery, and groundwater protection. Typically governed by gravity, viscous and capillary forces, these factors lead invasive fluids to occupy pore space irregularly and incompletely. Previous studies have demonstrated capillary numbers,describing the viscous and capillary forces, to quantificationally induce evolution of invasion patterns.While the evolution mechanisms of invasive patterns have not been deeply elucidated under the constant capillary number and three variable parameters including velocity, viscosity, and interfacial tension.Our research employs two horizontal visualization systems and a two-phase laminar flow simulation to investigate the tendency of invasive pattern transition by various parameters at the pore scale. We showed that increasing invasive viscosity or reducing interfacial tension in a homogeneous pore space significantly enhanced sweep efficiency, under constant capillary number. Additionally, in the fingering crossover pattern, the region near the inlet was prone to capillary fingering with multi-directional invasion, while the viscous fingering with unidirectional invasion was more susceptible occurred in the region near the outlet. Furthermore, increasing invasive viscosity or decreasing invasive velocity and interfacial tension promoted the extension of viscous fingering from the outlet to the inlet, presenting that the subsequent invasive fluid flows toward the outlet. In the case of invasive trunk along a unidirectional path, the invasive flow increased exponentially closer to the outlet, resulting in a significant decrease in the width of the invasive interface. Our work holds promising applications for optimizing invasive patterns in heterogeneous porous media.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

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.

Simulation of nanofluid natural convection based on single-particle hydrodynamics in energy-conserving dissipative particle dynamics(eDPD)

In the present study,the nanofliud natural convection is investigated by the energy-conserving dissipative particle dynamics(eDPD)method,where the nanoparticles are considered at the single-particle level.The thermal expansion coefficientβand the viscosityμof the simulated system containing nanoparticles are calculated and found to be in close alignment with the previous simulation results.The single-particle hydrodynamics in e DPD enables simulations of nanofluid natural convection with higher Rayleigh numbers and greater nanoparticle volume fractions.Additionally,this approach is utilized to simulate the nanoparticle distribution during the enhanced heat transfer process in the nanofluid natural convection.The localized aggregation of nanoparticles enhances the heat transfer performance of the nanofluid under specific Rayleigh numbers and nanoparticles volume fractions.

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.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

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.

Simulation of Two-Phase Flowback Phenomena in Shale Gas Wells

The gas-water two-phaseflow occurring as a result of fracturingfluidflowback phenomena is known to impact significantly the productivity of shale gas well.In this work,this two-phaseflow has been simulated in the framework of a hybrid approach partially relying on the embedded discrete fracture model(EDFM).This model assumes the region outside the stimulated reservoir volume(SRV)as a single-medium while the SRV region itself is described using a double-medium strategy which can account for thefluid exchange between the matrix and the micro-fractures.The shale gas adsorption,desorption,diffusion,gas slippage effect,fracture stress sensitivity,and capillary imbibition have been considered.The shale gas production,pore pressure distribution and water saturation distribution in the reservoir have been simulated.The influences of hydraulic fracture geometry and nonorthogonal hydraulic fractures on gas production have been determined and discussed accordingly.The simulation results show that the daily gas production has an upward and downward trend due to the presence of a large amount of fracturingfluid in the reservoir around the hydraulic fracture.The smaller the angle between the hydraulic fracture and the wellbore,the faster the daily production of shale gas wells decreases,and the lower the cumulative production.Nonplanar fractures can increase the control volume of hydraulic fractures and improve the production of shale gas wells.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

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.

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

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

Assessment of Low Global Warming Potential Refrigerants for Waste Heat Recovery in Data Center with On-Chip Two-Phase Cooling Loop

Data centers(DCs)are highly energy-intensive facilities,where about 30%–50%of the power consumed is attributable to the cooling of information technology equipment.This makes liquid cooling,especially in twophase mode,as an alternative to air cooling for the microprocessors in servers of interest.The need to meet the increased power density of server racks in high-performance DCs,along with the push towards lower global warming potential(GWP)refrigerants due to environmental concerns,has motivated research on the selection of two-phase heat transfer fluids for cooling servers while simultaneously recovering waste heat.With this regard,a heat pump-assisted absorption chiller(HPAAC)system for recovering waste heat in DCs with an on-chip twophase cooling loop driven by the compressor is proposed in the present paper and the low GWP hydrofluoroolefin refrigerants,including R1224yd(Z),R1233zd(E),R1234yf,R1234ze(E),R1234ze(Z),R1243zf and R1336mzz(Z),are evaluated and compared against R245fa as server coolant.For theHPAAC system,beginning with the development of energy and economic models,the performance is analyzed through both a parametric study and optimization using the coefficient of performance(COP),energy saving ratio(ESR),payback period(PBP)and net present value(NPV)as thermo-economic indicators.Using a standard vapor compression cooling system as a benchmark,the results indicate that with the evaporation temperature between 50℃and 70℃and the subcooling degree ranging from5℃to 15°C,R1233zd(E)with moderate compressor suction pressure and pressure ratio is the best refrigerant for the HPAAC systemwhile R1234yf performs the worst.More importantly,R1233zd(E)is also superior to R245fa based on thermo-economic performance,especially under work conditions with relatively lower evaporation temperature as well as subcooling degree.Under the given working conditions,the overall COP,ESR,NPV,and PBP of R1233zd(E)HPAAC with optimum subcooling degree range from4.99 to 11.27,25.53 to 64.59,1.13 to 4.10×10^(7) CNY and 5.77 to 2.22 years,respectively.Besides,the thermo-economic performance of R1233zd(E)HPAAC under optimum working conditions in terms of subcooling degree varying with the evaporation temperature is also investigated.