A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of ...A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensiti...We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.展开更多
Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive...Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.展开更多
In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two th...In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.展开更多
Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a ...Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.展开更多
This study introduces a coupled electromagnetic–thermal–mechanical model to reveal the mechanisms of microcracking and mineral melting of polymineralic rocks under microwave radiation.Experimental tests validate the...This study introduces a coupled electromagnetic–thermal–mechanical model to reveal the mechanisms of microcracking and mineral melting of polymineralic rocks under microwave radiation.Experimental tests validate the rationality of the proposed model.Embedding microscopic mineral sections into the granite model for simulation shows that uneven temperature gradients create distinct molten,porous,and nonmolten zones on the fracture surface.Moreover,the varying thermal expansion coefficients and Young's moduli among the minerals induce significant thermal stress at the mineral boundaries.Quartz and biotite with higher thermal expansion coefficients are subjected to compression,whereas plagioclase with smaller coefficients experiences tensile stress.In the molten zone,quartz undergoes transgranular cracking due to theα–βphase transition.The local high temperatures also induce melting phase transitions in biotite and feldspar.This numerical study provides new insights into the distribution of thermal stress and mineral phase changes in rocks under microwave irradiation.展开更多
The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research ...The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research on the unique characteristics of staged filling.This study aims to investigate the vacuum consolidation process of staged-filled soil slurry through laboratory model tests and numerical simulations,also assessing the impact of anionic polyacrylamide.Comparative analyses are conducted between vacuum consolidation with and without anionic polyacrylamide,as well as self-weight consolidation without anionic polyacrylamide.Results reveal contour lines of excess pore pressure,water content,and soil strength forming an ellipse around the prefabricated horizontal drain board.During the consolidation process,a higher degree of consolidation,lower water content,and higher soil strength were observed closer to the prefabricated horizontal drain board.After treatment,the uppermost filling layer exhibits an average water content that was approximately 40%higher than the lower filling layer,and its average strength was about 60%lower.This discrepancy is primarily due to the absence of sealing on the top surface and the relatively short vacuum consolidation time caused by staged filling.The introduction of anionic polyacrylamide-induced flocculation significantly improves the initial consolidation rate but minimally affects the dewatering capacity of vacuum preloading.Using flocculant can enhance both the staged filling rate and soil strength(by 1e2 times).Additionally,employing a staggered arrangement between different prefabricated horizontal drain layers is advisable to prevent top-down penetration in areas with low soil strength.展开更多
Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important me...Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.展开更多
Twisted and coiled polymer actuator(TCPA)is a type of artificial muscle that can be driven by heating due to its structure.A key issue with TCPA performance is the low driven frequency due to slow heat transfer in hea...Twisted and coiled polymer actuator(TCPA)is a type of artificial muscle that can be driven by heating due to its structure.A key issue with TCPA performance is the low driven frequency due to slow heat transfer in heating and cooling cycles,especially during cooling.We developed a numerical model of coating heating and nitrogen gas cooling that can effectively improve the driven forces and frequencies of the TCPA.Results indicate that natural cooling and electric fan cooling modes used in many experiments cannot restore the TCPA to its initial configuration when driven frequencies are high.Nitrogen gas cooling,at high driven frequencies,can fully restore the TCPA to its initial configuration,which is crucial for maintaining artificial muscle flexibility.In addition,as driven frequency increases,the corresponding driven force decreases.Systematic parametric studies were carried out to provide inspirations for optimizing TCPA design.The integrative computational study presented here provides a fundamental mechanistic understanding of the driven response in TCPA and sheds light on the rational design of TCPA through changing cooling modes.展开更多
The Pearl River Estuary(PRE)is one of China’s busiest shipping hubs and fishery production centers,as well as a region with abundant island tourism and wind energy resources,which calls for accurate short-term wind f...The Pearl River Estuary(PRE)is one of China’s busiest shipping hubs and fishery production centers,as well as a region with abundant island tourism and wind energy resources,which calls for accurate short-term wind forecasts.First,this study evaluated three operational numerical models,i.e.,ECMWF-EC,NCEP-GFS,and CMA-GD,for their ability to predict short-term wind speed over the PRE against in-situ observations during 2018-2021.Overall,ECMWF-EC out-performs other models with an average RMSE of 2.24 m s^(-1)and R of 0.57,but the NCEP-GFS performs better in the case of strong winds.Then,various bias correction and multi-model ensemble(MME)methods are used to perform the deterministic post-processing using a local and lead-specific scheme.Two-factor model output statistics(MOS2)is the optimal bias correction method for reducing(increasing)the overall RMSE(R)to 1.62(0.70)m s^(-1),demonstrating the benefits of considering both initial and lead-specific information.Intercomparison of MME results reveals that Multiple linear regression(MLR)presents superior skills,followed by random forest(RF),but it is slightly inferior to MOS2,particularly for the first few forecasting hours.Furthermore,the incorporation of additional features in MLR reduces the overall RMSE to 1.53 m s^(-1)and increases R to 0.74.Similarly,RF presents comparable results,and both outperform MOS2 in terms of correcting their deficiencies at the first few lead hours and limiting the error growth rate.Despite the satisfactory skill of deterministic post-processing techniques,they are unable to achieve a balanced performance between mean and extreme statistics.This highlights the necessity for further development of probabilistic forecasts.展开更多
Sporadic E(Es)layers in the ionosphere are characterized by intense plasma irregularities in the E region at altitudes of 90-130 km.Because they can significantly influence radio communications and navigation systems,...Sporadic E(Es)layers in the ionosphere are characterized by intense plasma irregularities in the E region at altitudes of 90-130 km.Because they can significantly influence radio communications and navigation systems,accurate forecasting of Es layers is crucial for ensuring the precision and dependability of navigation satellite systems.In this study,we present Es predictions made by an empirical model and by a deep learning model,and analyze their differences comprehensively by comparing the model predictions to satellite RO measurements and ground-based ionosonde observations.The deep learning model exhibited significantly better performance,as indicated by its high coefficient of correlation(r=0.87)with RO observations and predictions,than did the empirical model(r=0.53).This study highlights the importance of integrating artificial intelligence technology into ionosphere modelling generally,and into predicting Es layer occurrences and characteristics,in particular.展开更多
In this work, we present numerical modelling of coupled heat and mass transfer within porous materials. Our study focuses on cinder block bricks generally used in building construction. The material is assumed to be p...In this work, we present numerical modelling of coupled heat and mass transfer within porous materials. Our study focuses on cinder block bricks generally used in building construction. The material is assumed to be placed in air. Moisture content and temperature have been chosen as the main transfer drivers and the equations governing these transfer drivers are based on the Luikov model. These equations are solved by an implicit finite difference scheme. A Fortran code associated with the Thomas algorithm was used to solve the equations. The results show that heat and mass transfer depend on the temperature of the air in contact with the material. As this air temperature rises, the temperature within the material increases, and more rapidly at the material surface. Also, thermal conductivity plays a very important role in the thermal conduction of building materials and influences heat and mass transfer in these materials. Materials with higher thermal conductivity diffuse more heat.展开更多
Deep coalbed methane(CBM)has become one of the most significant potential sources of natural gas in China.However,the exploration and development of deep CBM in China is still in an initial stage,and its accumulation-...Deep coalbed methane(CBM)has become one of the most significant potential sources of natural gas in China.However,the exploration and development of deep CBM in China is still in an initial stage,and its accumulation-forming characteristics require further study.Therefore,taking the No.8 deep coal seam in the central-eastern region of Ordos Basin as an example,this study investigated the geologic characteristics of CBM accumulations to establish a numerical model.The evolution of the burial and accumulation of CBM in the area was reconstructed.The modeling results suggest that the No.8 coal seam experienced continuous subsidence from the Late Cretaceous to the Triassic,alternating subsidence and uplift during the Jurassic,rapid burial throughout the Early Cretaceous,and continuous uplift since the Late Cretaceous.The coal reached its maximum maturity at the end of the Early Cretaceous.Furthermore,CBM generation in the region was divided into four stages of thermal eventsdbiogenic and early thermogenic gas,cracking of light oil into gas,cracking of the remaining kerogen into gas,and hydrocarbon generation ceasingdwhich accelerated coal maturity and generation.The adsorption capacity presented an overall declining trend prior to the end of the Cretaceous,followed by a rapid increase since the Late Cretaceous.As for adsorption mass evolution,the CBM successively un-derwent unsaturated minor adsorption,unsaturated rapid-rising adsorption,saturated decreasing adsorption,and saturated rising adsorption.The in-situ gas mass was found to be controlled by a combination of generation,adsorption,and expulsion of hydrocarbons,with its present-day value being 9-29×10^(4)t/km^(2)and the corresponding gas volume per ton of coal being 12-28 m^(3)/t.Moreover,free gas evolution initially showed an increasing trend,followed by a decline,ultimately accounting for 11%-28%of the total gas content.展开更多
Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model o...Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in<sup>2</sup>/s for humidity, 2.993 in<sup>2</sup>/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .展开更多
Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical e...Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical experimental measurement and numerical simulation pose research challenges.This study focuses on the ice load of a cylinder structure breaking upward through the ice sheet form underneath in the Small Ice Model Basin of China Ship Scientific Research Center(CSSRC SIMB).A high-speed camera system was employed to observe the ice sheet failure during the tests,in which,with the loading position as center,local radial cracks and circumferential cracks were generated.A load sensor was used to measure the overall ice load during this process.Meanwhile,a numerical model was developed using LS-DYNA for validation and comparison.With this model,numerical simulation was conducted under various ice thicknesses and upgoing speeds to analyze the instantaneous curves of ice load.The calculation results were statistically analyzed under different working conditions to determine the influence of the factors on the ice load of the cylinder.The study explores the measurement method about ice load of objects vertically breaking through model ice sheet and is expected to provide some fundamental insights into the safety design of underwater structures operating in ice waters.展开更多
Numerical modelling is an effective technique to improve the understanding of outburst initiation mechanisms and to take appropriate measures to address their threats.Based on the existing two-way sequential coupling ...Numerical modelling is an effective technique to improve the understanding of outburst initiation mechanisms and to take appropriate measures to address their threats.Based on the existing two-way sequential coupling method,two typical types of outbursts,i.e.the gas pocket outburst and the dynamic fracturing outburst,have been successfully simulated using field data from a coalfield in central China.The geological structure commonly observed in the coalfield,known as the‘bedding shear zone’,contributes to the gas pocket outbursts in the region.The model for this type of outburst simulates mininginduced stress and gas pressure distributions during the outburst initiation stage and the subsequent development stage.Both coal ejection and gas release are observed in the model,and the simulation results are consistent with mine site observations,i.e.the amount of ejected coal,outburst cavity profile,and gas release rate change prior to an outburst.The second type of outburst is attributed to gas accumulation and elevated gas pressure due to the gassy floor seam and the heterogeneity in the floor strata,which is explained by the dynamic fracturing theory.While the dynamic coal ejection phenomenon is not captured in the simulation,the abrupt release of retained gas from a floor coal seam is successfully replicated.Both outburst models reveal that abnormal gas emission trends can be used as indicators of an upcoming outburst.The results of this study are expected to provide new insights into the outburst initiation mechanisms and outburst prevention measures.展开更多
A three-dimensional numerical model of sand wave dynamics,incorporating the interaction of currents and waves at various angles,has been developed using the Regional Ocean Modeling System(ROMS).This model accounts for...A three-dimensional numerical model of sand wave dynamics,incorporating the interaction of currents and waves at various angles,has been developed using the Regional Ocean Modeling System(ROMS).This model accounts for both bedload and suspended load sediment transport under combined waves and current conditions.The investigation examines the influence of several key parameters,including the rotation angle of sand waves relative to the main current,tidal current velocity amplitude,residual current,water depth,wave height,wave period,and wave direction,on sand wave evolution.The growth rate and migration rate of sand waves decrease as their rotation angle increases.For rotation angles smaller than 15°,sand wave evolution can be effectively simulated by a vertical 2D model with an error within 10%.The numerical results demonstrate that variations in tidal current velocity amplitude or residual current affect both vertical growth and horizontal migration of sand waves.As tidal current velocity amplitude and residual current increase,the growth rate initially rises to a maximum before decreasing.The migration rate shows a consistent increase with increasing tidal current amplitude and residual current.Under combined waves and current,both growth and migration rates decrease as water depth increases.With increasing wave height and period,the growth rate and migration rate initially rise to maximum values before declining,while showing a consistent increase with wave height and period.The change rate of sand waves reaches its maximum when wave propagation aligns parallel to tidal currents,and reaches its minimum when wave propagation is perpendicular to the currents.This phenomenon can be explained by the fluctuation of total bed shear stress relative to the angle of interaction between waves and current.展开更多
Open channel confluences,where two streams or rivers converge,play a crucial role in hydraulic engineering and river dynamics.These confluences are characterized by complex hydrodynamics influenced by the discharge ra...Open channel confluences,where two streams or rivers converge,play a crucial role in hydraulic engineering and river dynamics.These confluences are characterized by complex hydrodynamics influenced by the discharge ratios of merging water bodies.This study investigated the mixing structure at open channel confluences using three-dimensional numerical modeling.A comprehensive three-dimensional numerical model was developed and validated against a dataset obtained from controlled laboratory experiments.This dataset incorporated three-dimensional time-averaged velocity measurements.The skew-induced and stress-induced equation systems were adopted as the core governing equations,providing a framework for simulating various scenarios.A total of ten different cases were analyzed.The results highlighted the effect of discharge ratios on turbulence,lateral and vertical vorticities,and the distribution of mixing,which intensified with higher magnitudes of discharge ratios.The mixing structure,driven by velocity gradients and vorticity,revealed the significant role of lateral and vertical vorticities in determining hydrodynamic behaviors and mixing distributions at confluences.Specifically,the momentum ratio of incoming flows governed the spatial evolution of mixing processes.This study revealed that the distribution of mixing served as a key indicator for identifying the formation of mid-channel scours.High normalized velocities induced toward the left bank led to the superelevation of the water surface,enhancing the potential for bed material and the formation of significant scour holes beneath the elevated water surface.This novel approach provides a deeper understanding of the mixing patterns at confluences,particularly in scenarios with equilibrated discharge ratios but in different magnitudes.展开更多
Semi-crystalline polymer laser powder bed fusion(L-PBF)has recently attracted increasing interest due to its potential for fabricating complex geometry.However,a more comprehensive understanding of the underlying phys...Semi-crystalline polymer laser powder bed fusion(L-PBF)has recently attracted increasing interest due to its potential for fabricating complex geometry.However,a more comprehensive understanding of the underlying physics during L-PBF is required to better control the properties of the final part.This work proposed a multi-layer numerical model to study the temperature and phase evolution during the polyamide-12(PA12)L-PBF process.The Descend and Parallel Chord methods were introduced to improve the convergence of the non-linear thermal solver.The level-set-based mesh adaptation strategy,governed by multi-physical fields,was applied to alleviate the calculation and accurately track the phase evolution.The processing simulation on the dog-bone model revealed that preheating temperature significantly influences the crystallization behavior.Finally,the multi-layer simulation demonstrated that such a developed numerical model can be used to study the phase transformation during powder layer updating and the cyclic laser sintering phenomena.Moreover,the numerical study suggested that crystallization occurs slowly during the L-PBF process.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.52375340,51975263,52405366).
文摘A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
基金funding received by a grant from the Natural Sciences and Engineering Research Council of Canada(NSERC)(Grant No.CRDPJ 469057e14).
文摘We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.
基金support from the National Key R&D plan(Grant No.2022YFC3004303)the National Natural Science Foundation of China(Grant No.42107161)+3 种基金the State Key Laboratory of Hydroscience and Hydraulic Engineering(Grant No.2021-KY-04)the Open Research Fund Program of State Key Laboratory of Hydroscience and Engineering(sklhse-2023-C-01)the Open Research Fund Program of Key Laboratory of the Hydrosphere of the Ministry of Water Resources(mklhs-2023-04)the China Three Gorges Corporation(XLD/2117).
文摘Rock fragmentation plays a critical role in rock avalanches,yet conventional approaches such as classical granular flow models or the bonded particle model have limitations in accurately characterizing the progressive disintegration and kinematics of multi-deformable rock blocks during rockslides.The present study proposes a discrete-continuous numerical model,based on a cohesive zone model,to explicitly incorporate the progressive fragmentation and intricate interparticle interactions inherent in rockslides.Breakable rock granular assemblies are released along an inclined plane and flow onto a horizontal plane.The numerical scenarios are established to incorporate variations in slope angle,initial height,friction coefficient,and particle number.The evolutions of fragmentation,kinematic,runout and depositional characteristics are quantitatively analyzed and compared with experimental and field data.A positive linear relationship between the equivalent friction coefficient and the apparent friction coefficient is identified.In general,the granular mass predominantly exhibits characteristics of a dense granular flow,with the Savage number exhibiting a decreasing trend as the volume of mass increases.The process of particle breakage gradually occurs in a bottom-up manner,leading to a significant increase in the angular velocities of the rock blocks with increasing depth.The simulation results reproduce the field observations of inverse grading and source stratigraphy preservation in the deposit.We propose a disintegration index that incorporates factors such as drop height,rock mass volume,and rock strength.Our findings demonstrate a consistent linear relationship between this index and the fragmentation degree in all tested scenarios.
基金supported by project XJZ2023050044,A2309002 and XJZ2023070052.
文摘In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.
基金Engineering and Physical Sciences Research Council (EPSRC) is also acknowledged for funding this work under Grant Number EP/N009207/1.
文摘Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.
基金the National Natural Science Foundation of China(No.52074349)the Graduate Research Innovation Project of Hunan Province,China(No.CX20230194)。
文摘This study introduces a coupled electromagnetic–thermal–mechanical model to reveal the mechanisms of microcracking and mineral melting of polymineralic rocks under microwave radiation.Experimental tests validate the rationality of the proposed model.Embedding microscopic mineral sections into the granite model for simulation shows that uneven temperature gradients create distinct molten,porous,and nonmolten zones on the fracture surface.Moreover,the varying thermal expansion coefficients and Young's moduli among the minerals induce significant thermal stress at the mineral boundaries.Quartz and biotite with higher thermal expansion coefficients are subjected to compression,whereas plagioclase with smaller coefficients experiences tensile stress.In the molten zone,quartz undergoes transgranular cracking due to theα–βphase transition.The local high temperatures also induce melting phase transitions in biotite and feldspar.This numerical study provides new insights into the distribution of thermal stress and mineral phase changes in rocks under microwave irradiation.
基金supported by the Research Grants Council of Hong Kong Special Administrative Region Government of China(Grant Nos.15210322 and R5037-18)the financial support(Grant No.86902-00000240)from Shenzhen University.
文摘The vacuum-assisted prefabricated horizontal drain offers a promising method for strengthening soil slurry,allowing simultaneous filling and vacuum-dewatering via staged construction.However,there is limited research on the unique characteristics of staged filling.This study aims to investigate the vacuum consolidation process of staged-filled soil slurry through laboratory model tests and numerical simulations,also assessing the impact of anionic polyacrylamide.Comparative analyses are conducted between vacuum consolidation with and without anionic polyacrylamide,as well as self-weight consolidation without anionic polyacrylamide.Results reveal contour lines of excess pore pressure,water content,and soil strength forming an ellipse around the prefabricated horizontal drain board.During the consolidation process,a higher degree of consolidation,lower water content,and higher soil strength were observed closer to the prefabricated horizontal drain board.After treatment,the uppermost filling layer exhibits an average water content that was approximately 40%higher than the lower filling layer,and its average strength was about 60%lower.This discrepancy is primarily due to the absence of sealing on the top surface and the relatively short vacuum consolidation time caused by staged filling.The introduction of anionic polyacrylamide-induced flocculation significantly improves the initial consolidation rate but minimally affects the dewatering capacity of vacuum preloading.Using flocculant can enhance both the staged filling rate and soil strength(by 1e2 times).Additionally,employing a staggered arrangement between different prefabricated horizontal drain layers is advisable to prevent top-down penetration in areas with low soil strength.
基金The National Natural Science Foundation of China under contract No.42076214.
文摘Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.
基金Natural Science Foundation of China,11925204,Jizeng Wangthe Fundamental Research Funds for the Central Universities,lzujbky-2024-jdzx02,Zhiwen Gao。
文摘Twisted and coiled polymer actuator(TCPA)is a type of artificial muscle that can be driven by heating due to its structure.A key issue with TCPA performance is the low driven frequency due to slow heat transfer in heating and cooling cycles,especially during cooling.We developed a numerical model of coating heating and nitrogen gas cooling that can effectively improve the driven forces and frequencies of the TCPA.Results indicate that natural cooling and electric fan cooling modes used in many experiments cannot restore the TCPA to its initial configuration when driven frequencies are high.Nitrogen gas cooling,at high driven frequencies,can fully restore the TCPA to its initial configuration,which is crucial for maintaining artificial muscle flexibility.In addition,as driven frequency increases,the corresponding driven force decreases.Systematic parametric studies were carried out to provide inspirations for optimizing TCPA design.The integrative computational study presented here provides a fundamental mechanistic understanding of the driven response in TCPA and sheds light on the rational design of TCPA through changing cooling modes.
基金Science and Technology Research Project of Guangdong Meteorological Service(GRMC2021M19,GRMC2022Q16,GRMC2023M29)。
文摘The Pearl River Estuary(PRE)is one of China’s busiest shipping hubs and fishery production centers,as well as a region with abundant island tourism and wind energy resources,which calls for accurate short-term wind forecasts.First,this study evaluated three operational numerical models,i.e.,ECMWF-EC,NCEP-GFS,and CMA-GD,for their ability to predict short-term wind speed over the PRE against in-situ observations during 2018-2021.Overall,ECMWF-EC out-performs other models with an average RMSE of 2.24 m s^(-1)and R of 0.57,but the NCEP-GFS performs better in the case of strong winds.Then,various bias correction and multi-model ensemble(MME)methods are used to perform the deterministic post-processing using a local and lead-specific scheme.Two-factor model output statistics(MOS2)is the optimal bias correction method for reducing(increasing)the overall RMSE(R)to 1.62(0.70)m s^(-1),demonstrating the benefits of considering both initial and lead-specific information.Intercomparison of MME results reveals that Multiple linear regression(MLR)presents superior skills,followed by random forest(RF),but it is slightly inferior to MOS2,particularly for the first few forecasting hours.Furthermore,the incorporation of additional features in MLR reduces the overall RMSE to 1.53 m s^(-1)and increases R to 0.74.Similarly,RF presents comparable results,and both outperform MOS2 in terms of correcting their deficiencies at the first few lead hours and limiting the error growth rate.Despite the satisfactory skill of deterministic post-processing techniques,they are unable to achieve a balanced performance between mean and extreme statistics.This highlights the necessity for further development of probabilistic forecasts.
基金supported by the Project of Stable Support for Youth Team in Basic Research Field,CAS(grant No.YSBR-018)the National Natural Science Foundation of China(grant Nos.42188101,42130204)+4 种基金the B-type Strategic Priority Program of CAS(grant no.XDB41000000)the National Natural Science Foundation of China(NSFC)Distinguished Overseas Young Talents Program,Innovation Program for Quantum Science and Technology(2021ZD0300301)the Open Research Project of Large Research Infrastructures of CAS-“Study on the interaction between low/mid-latitude atmosphere and ionosphere based on the Chinese Meridian Project”.The project was supported also by the National Key Laboratory of Deep Space Exploration(Grant No.NKLDSE2023A002)the Open Fund of Anhui Provincial Key Laboratory of Intelligent Underground Detection(Grant No.APKLIUD23KF01)the China National Space Administration(CNSA)pre-research Project on Civil Aerospace Technologies No.D010305,D010301.
文摘Sporadic E(Es)layers in the ionosphere are characterized by intense plasma irregularities in the E region at altitudes of 90-130 km.Because they can significantly influence radio communications and navigation systems,accurate forecasting of Es layers is crucial for ensuring the precision and dependability of navigation satellite systems.In this study,we present Es predictions made by an empirical model and by a deep learning model,and analyze their differences comprehensively by comparing the model predictions to satellite RO measurements and ground-based ionosonde observations.The deep learning model exhibited significantly better performance,as indicated by its high coefficient of correlation(r=0.87)with RO observations and predictions,than did the empirical model(r=0.53).This study highlights the importance of integrating artificial intelligence technology into ionosphere modelling generally,and into predicting Es layer occurrences and characteristics,in particular.
文摘In this work, we present numerical modelling of coupled heat and mass transfer within porous materials. Our study focuses on cinder block bricks generally used in building construction. The material is assumed to be placed in air. Moisture content and temperature have been chosen as the main transfer drivers and the equations governing these transfer drivers are based on the Luikov model. These equations are solved by an implicit finite difference scheme. A Fortran code associated with the Thomas algorithm was used to solve the equations. The results show that heat and mass transfer depend on the temperature of the air in contact with the material. As this air temperature rises, the temperature within the material increases, and more rapidly at the material surface. Also, thermal conductivity plays a very important role in the thermal conduction of building materials and influences heat and mass transfer in these materials. Materials with higher thermal conductivity diffuse more heat.
基金funded by the applied research and technology project titled“Research on Key Technologies for Optimization Design of Deep Coalbed Methane Development”conducted by PetroChina Company Limited(Grant No.2023ZZ1804).
文摘Deep coalbed methane(CBM)has become one of the most significant potential sources of natural gas in China.However,the exploration and development of deep CBM in China is still in an initial stage,and its accumulation-forming characteristics require further study.Therefore,taking the No.8 deep coal seam in the central-eastern region of Ordos Basin as an example,this study investigated the geologic characteristics of CBM accumulations to establish a numerical model.The evolution of the burial and accumulation of CBM in the area was reconstructed.The modeling results suggest that the No.8 coal seam experienced continuous subsidence from the Late Cretaceous to the Triassic,alternating subsidence and uplift during the Jurassic,rapid burial throughout the Early Cretaceous,and continuous uplift since the Late Cretaceous.The coal reached its maximum maturity at the end of the Early Cretaceous.Furthermore,CBM generation in the region was divided into four stages of thermal eventsdbiogenic and early thermogenic gas,cracking of light oil into gas,cracking of the remaining kerogen into gas,and hydrocarbon generation ceasingdwhich accelerated coal maturity and generation.The adsorption capacity presented an overall declining trend prior to the end of the Cretaceous,followed by a rapid increase since the Late Cretaceous.As for adsorption mass evolution,the CBM successively un-derwent unsaturated minor adsorption,unsaturated rapid-rising adsorption,saturated decreasing adsorption,and saturated rising adsorption.The in-situ gas mass was found to be controlled by a combination of generation,adsorption,and expulsion of hydrocarbons,with its present-day value being 9-29×10^(4)t/km^(2)and the corresponding gas volume per ton of coal being 12-28 m^(3)/t.Moreover,free gas evolution initially showed an increasing trend,followed by a decline,ultimately accounting for 11%-28%of the total gas content.
文摘Atmospheric models are physical equations based on the ideal gas law. Applied to the atmosphere, this law yields equations for water, vapor (gas), ice, air, humidity, dryness, fire, and heat, thus defining the model of key atmospheric parameters. The distribution of these parameters across the entire planet Earth is the origin of the formation of the climatic cycle, which is a normal climatic variation. To do this, the Earth is divided into eight (8) parts according to the number of key parameters to be defined in a physical representation of the model. Following this distribution, numerical models calculate the constants for the formation of water, vapor, ice, dryness, thermal energy (fire), heat, air, and humidity. These models vary in complexity depending on the indirect trigonometric direction and simplicity in the sum of neighboring models. Note that the constants obtained from the equations yield 275.156˚K (2.006˚C) for water, 273.1596˚K (0.00963˚C) for vapor, 273.1633˚K (0.0133˚C) for ice, 0.00365 in/s for atmospheric dryness, 1.996 in<sup>2</sup>/s for humidity, 2.993 in<sup>2</sup>/s for air, 1 J for thermal energy of fire, and 0.9963 J for heat. In summary, this study aims to define the main parameters and natural phenomena contributing to the modification of planetary climate. .
文摘Ice load on underwater vehicles breaking through ice covers from underneath is a significant concern for researchers in polar exploration,and the research on this problem is still in its early stages.Both mechanical experimental measurement and numerical simulation pose research challenges.This study focuses on the ice load of a cylinder structure breaking upward through the ice sheet form underneath in the Small Ice Model Basin of China Ship Scientific Research Center(CSSRC SIMB).A high-speed camera system was employed to observe the ice sheet failure during the tests,in which,with the loading position as center,local radial cracks and circumferential cracks were generated.A load sensor was used to measure the overall ice load during this process.Meanwhile,a numerical model was developed using LS-DYNA for validation and comparison.With this model,numerical simulation was conducted under various ice thicknesses and upgoing speeds to analyze the instantaneous curves of ice load.The calculation results were statistically analyzed under different working conditions to determine the influence of the factors on the ice load of the cylinder.The study explores the measurement method about ice load of objects vertically breaking through model ice sheet and is expected to provide some fundamental insights into the safety design of underwater structures operating in ice waters.
基金the National Natural Science Foundation of China(52304105)National Natural Science Foundation of China-National major scientific research instrument development project(52227901)Jiangsu Province International Collaboration Program-Key national industrial technology research and development cooperation projects(BZ2023050).
文摘Numerical modelling is an effective technique to improve the understanding of outburst initiation mechanisms and to take appropriate measures to address their threats.Based on the existing two-way sequential coupling method,two typical types of outbursts,i.e.the gas pocket outburst and the dynamic fracturing outburst,have been successfully simulated using field data from a coalfield in central China.The geological structure commonly observed in the coalfield,known as the‘bedding shear zone’,contributes to the gas pocket outbursts in the region.The model for this type of outburst simulates mininginduced stress and gas pressure distributions during the outburst initiation stage and the subsequent development stage.Both coal ejection and gas release are observed in the model,and the simulation results are consistent with mine site observations,i.e.the amount of ejected coal,outburst cavity profile,and gas release rate change prior to an outburst.The second type of outburst is attributed to gas accumulation and elevated gas pressure due to the gassy floor seam and the heterogeneity in the floor strata,which is explained by the dynamic fracturing theory.While the dynamic coal ejection phenomenon is not captured in the simulation,the abrupt release of retained gas from a floor coal seam is successfully replicated.Both outburst models reveal that abnormal gas emission trends can be used as indicators of an upcoming outburst.The results of this study are expected to provide new insights into the outburst initiation mechanisms and outburst prevention measures.
基金the National Natural Science Foundation of China(Grant Nos.52371289 and 51979192).
文摘A three-dimensional numerical model of sand wave dynamics,incorporating the interaction of currents and waves at various angles,has been developed using the Regional Ocean Modeling System(ROMS).This model accounts for both bedload and suspended load sediment transport under combined waves and current conditions.The investigation examines the influence of several key parameters,including the rotation angle of sand waves relative to the main current,tidal current velocity amplitude,residual current,water depth,wave height,wave period,and wave direction,on sand wave evolution.The growth rate and migration rate of sand waves decrease as their rotation angle increases.For rotation angles smaller than 15°,sand wave evolution can be effectively simulated by a vertical 2D model with an error within 10%.The numerical results demonstrate that variations in tidal current velocity amplitude or residual current affect both vertical growth and horizontal migration of sand waves.As tidal current velocity amplitude and residual current increase,the growth rate initially rises to a maximum before decreasing.The migration rate shows a consistent increase with increasing tidal current amplitude and residual current.Under combined waves and current,both growth and migration rates decrease as water depth increases.With increasing wave height and period,the growth rate and migration rate initially rise to maximum values before declining,while showing a consistent increase with wave height and period.The change rate of sand waves reaches its maximum when wave propagation aligns parallel to tidal currents,and reaches its minimum when wave propagation is perpendicular to the currents.This phenomenon can be explained by the fluctuation of total bed shear stress relative to the angle of interaction between waves and current.
文摘Open channel confluences,where two streams or rivers converge,play a crucial role in hydraulic engineering and river dynamics.These confluences are characterized by complex hydrodynamics influenced by the discharge ratios of merging water bodies.This study investigated the mixing structure at open channel confluences using three-dimensional numerical modeling.A comprehensive three-dimensional numerical model was developed and validated against a dataset obtained from controlled laboratory experiments.This dataset incorporated three-dimensional time-averaged velocity measurements.The skew-induced and stress-induced equation systems were adopted as the core governing equations,providing a framework for simulating various scenarios.A total of ten different cases were analyzed.The results highlighted the effect of discharge ratios on turbulence,lateral and vertical vorticities,and the distribution of mixing,which intensified with higher magnitudes of discharge ratios.The mixing structure,driven by velocity gradients and vorticity,revealed the significant role of lateral and vertical vorticities in determining hydrodynamic behaviors and mixing distributions at confluences.Specifically,the momentum ratio of incoming flows governed the spatial evolution of mixing processes.This study revealed that the distribution of mixing served as a key indicator for identifying the formation of mid-channel scours.High normalized velocities induced toward the left bank led to the superelevation of the water surface,enhancing the potential for bed material and the formation of significant scour holes beneath the elevated water surface.This novel approach provides a deeper understanding of the mixing patterns at confluences,particularly in scenarios with equilibrated discharge ratios but in different magnitudes.
文摘Semi-crystalline polymer laser powder bed fusion(L-PBF)has recently attracted increasing interest due to its potential for fabricating complex geometry.However,a more comprehensive understanding of the underlying physics during L-PBF is required to better control the properties of the final part.This work proposed a multi-layer numerical model to study the temperature and phase evolution during the polyamide-12(PA12)L-PBF process.The Descend and Parallel Chord methods were introduced to improve the convergence of the non-linear thermal solver.The level-set-based mesh adaptation strategy,governed by multi-physical fields,was applied to alleviate the calculation and accurately track the phase evolution.The processing simulation on the dog-bone model revealed that preheating temperature significantly influences the crystallization behavior.Finally,the multi-layer simulation demonstrated that such a developed numerical model can be used to study the phase transformation during powder layer updating and the cyclic laser sintering phenomena.Moreover,the numerical study suggested that crystallization occurs slowly during the L-PBF process.
