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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic th...Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.展开更多
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.展开更多
With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ...With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.展开更多
Understanding the behaviour of composite marine propellers during operating conditions is a need of the present era since they emerge as a potential replacement for conventional propeller materials such as metals or a...Understanding the behaviour of composite marine propellers during operating conditions is a need of the present era since they emerge as a potential replacement for conventional propeller materials such as metals or alloys.They offer several benefits,such as high specific strength,low corrosion,delayed cavitation,improved dynamic stability,reduced noise levels,and overall energy efficiency.In addition,composite materials undergo passive deformation,termed as“bend-twist effect”,under hydrodynamic loads due to their inherent flexibility and anisotropy.Although performance analysis methods were developed in the past for marine propellers,there is a significant lack of literature on composite propellers.This article discusses the recent advancements in experimental and numerical modelling,state-of-the-art computational technologies,and mutated mathematical models that aid in designing,analysing,and optimising composite marine propellers.In the initial sections,performance evaluation methods and challenges with the existing propeller materials are discussed.Thereafter,the benefits of composite propellers are critically reviewed.Numerical and experimental FSI coupling methods,cavitation performance,the effect of stacking sequence,and acoustic measurements are some critical areas discussed in detail.A two-way FSI-coupled simulation was conducted in a non-cavitating regime for four advanced ratios and compared with the literature results.Finally,the scope for future improvements and conclusions are mentioned.展开更多
Despite the long-standing success of waste tailings recycling at shallower lithosphere,heat influx from the stifling environment of deep strata and ground squeezing due to aggravating deposit rheology have posed subst...Despite the long-standing success of waste tailings recycling at shallower lithosphere,heat influx from the stifling environment of deep strata and ground squeezing due to aggravating deposit rheology have posed substantial challenges to traditional mine backfilling.To explore the minefill behavior in difficult geology settings,a new thermo-poroelastic model is developed in this study for characterizing the porepressure evolution in hydrating backfill subjected to coupled thermal perturbation and mechanical deformation.By scrutinizing the undrained pressure responses to a hierarchy of thermo-mechanical loading configurations,we investigated the relative contribution of heat exchange and wall convergence to the macroscopic minefill behavior.The result suggests that the backfill response could exhibit unique rate and path-dependence upon thermal loading due to the competing processes that govern the pressure evolution.Moreover,it shows that while rapid wall convergence would always induce significant pressure,thermal straining via heat exchange might still modulate substantially the pressure distribution when the temperature catalyzes considerable water expansivity.These findings could facilitate understanding of the complex backfill behavior in difficult geology settings,and thus have significant implications for adaptive tailings management to deep mining development.展开更多
In the realm of all-electric aircraft research,the integration of cathode-open proton exchange membrane fuel cells(PEMFC)with lithiumbatteries as a hybrid power source for small to medium-sized unmanned aerial vehicle...In the realm of all-electric aircraft research,the integration of cathode-open proton exchange membrane fuel cells(PEMFC)with lithiumbatteries as a hybrid power source for small to medium-sized unmanned aerial vehicles(UAVs)has garnered significant attention.The PEMFC,serving as the primary energy supply,markedly extends the UAV’s operational endurance.However,due to payload limitations and spatial constraints in the airframe layout of UAVs,the stack requires customized adaptation.Moreover,the implementation of auxiliary systems to facilitate cold starts of the PEMFC under low-temperature conditions is not feasible.Relying solely on thermal insulation measures also proves inadequate to address the challenges posed by complex low-temperature startup scenarios.To overcomethis,our study leverages the UAV’s lithium battery to heat the cathode inlet airflow,aiding the cathode-open PEMFC cold start process.To validate the feasibility of the proposed air-assisted heating strategy during the conceptual design phase,this study develops a transient,non-isothermal 3Dcathode-open PEMF Cunitmodel incorporating cathode air-assisted heating and gas-ice phase change.The model’s accuracy was verified against experimental cold-start data from a stack composed of identical single cells.This computational framework enables quantitative analysis of temperature fields and ice fraction distributions across domains under varying air-assisted heating powers during cold starts.Building upon this model,the study further investigates the improvement in cold start performance by heating the cathode intake air with varying power levels.The results demonstrate that the fuel cell achieves self-startup at temperatures as low as−13℃ under a constant current density of 100mA/cm^(2) without air-assisted heating.At an ambient temperature of−20℃,a successful start-up can be achieved with a heating power of 0.45 W/cm^(2).The temperature variation overtime during the cold start process can be represented by a sum of two exponential functions.The air-assisted heating scheme proposed in this study has significantly improved the cold start performance of fuel cells in low-temperature environments.Additionally,it provides critical reference data and validation support for component selection and feasibility assessment of hybrid power systems.展开更多
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.展开更多
To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of t...To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.展开更多
This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in r...This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The ...In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The grid can be discretionarily refined as a non-structure triangle or quadrilateral so that piers can be treated as one or several impermeable elements with an area of 20 to 30 km^2 in a model domain over 85 700 km^2. The computational results of tidal levels and horizontal velocities are in good agreement with the field data. Based on the computed results by the model, the layout of an open 105 DWT liquefied natural gas （LNG）terminal in Caofeidian, Bohal Sea is effectively and reasonably optimized. It can be concluded that the model is suitable and reasonable for direct simulation of tidal currents around small structures in projects.展开更多
基金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.
基金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.
基金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.
基金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.
基金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.
基金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.
文摘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.
基金The National Key R&D Program of China under contract No.2023YFE0126300the National Natural Science Foundation of China under contract Nos 42066002 and U20A2099.
文摘Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.
文摘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.
文摘With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.
基金Supporting by the project‘FILE NO.CRG/2022/001718’.
文摘Understanding the behaviour of composite marine propellers during operating conditions is a need of the present era since they emerge as a potential replacement for conventional propeller materials such as metals or alloys.They offer several benefits,such as high specific strength,low corrosion,delayed cavitation,improved dynamic stability,reduced noise levels,and overall energy efficiency.In addition,composite materials undergo passive deformation,termed as“bend-twist effect”,under hydrodynamic loads due to their inherent flexibility and anisotropy.Although performance analysis methods were developed in the past for marine propellers,there is a significant lack of literature on composite propellers.This article discusses the recent advancements in experimental and numerical modelling,state-of-the-art computational technologies,and mutated mathematical models that aid in designing,analysing,and optimising composite marine propellers.In the initial sections,performance evaluation methods and challenges with the existing propeller materials are discussed.Thereafter,the benefits of composite propellers are critically reviewed.Numerical and experimental FSI coupling methods,cavitation performance,the effect of stacking sequence,and acoustic measurements are some critical areas discussed in detail.A two-way FSI-coupled simulation was conducted in a non-cavitating regime for four advanced ratios and compared with the literature results.Finally,the scope for future improvements and conclusions are mentioned.
基金jointly supported by the Natural Science Foundation of Xinjiang Uygur Autonomous Region(Grant No.2022D01E32)Opening Fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(Grant No.SKLGP 2023K012)the National Natural Science Foundation of China(Grant No.U20A20266).
文摘Despite the long-standing success of waste tailings recycling at shallower lithosphere,heat influx from the stifling environment of deep strata and ground squeezing due to aggravating deposit rheology have posed substantial challenges to traditional mine backfilling.To explore the minefill behavior in difficult geology settings,a new thermo-poroelastic model is developed in this study for characterizing the porepressure evolution in hydrating backfill subjected to coupled thermal perturbation and mechanical deformation.By scrutinizing the undrained pressure responses to a hierarchy of thermo-mechanical loading configurations,we investigated the relative contribution of heat exchange and wall convergence to the macroscopic minefill behavior.The result suggests that the backfill response could exhibit unique rate and path-dependence upon thermal loading due to the competing processes that govern the pressure evolution.Moreover,it shows that while rapid wall convergence would always induce significant pressure,thermal straining via heat exchange might still modulate substantially the pressure distribution when the temperature catalyzes considerable water expansivity.These findings could facilitate understanding of the complex backfill behavior in difficult geology settings,and thus have significant implications for adaptive tailings management to deep mining development.
基金funded by Zhejiang Province Spearhead and Leading Goose Research and Development Key Program,grant number 2023C01239.
文摘In the realm of all-electric aircraft research,the integration of cathode-open proton exchange membrane fuel cells(PEMFC)with lithiumbatteries as a hybrid power source for small to medium-sized unmanned aerial vehicles(UAVs)has garnered significant attention.The PEMFC,serving as the primary energy supply,markedly extends the UAV’s operational endurance.However,due to payload limitations and spatial constraints in the airframe layout of UAVs,the stack requires customized adaptation.Moreover,the implementation of auxiliary systems to facilitate cold starts of the PEMFC under low-temperature conditions is not feasible.Relying solely on thermal insulation measures also proves inadequate to address the challenges posed by complex low-temperature startup scenarios.To overcomethis,our study leverages the UAV’s lithium battery to heat the cathode inlet airflow,aiding the cathode-open PEMFC cold start process.To validate the feasibility of the proposed air-assisted heating strategy during the conceptual design phase,this study develops a transient,non-isothermal 3Dcathode-open PEMF Cunitmodel incorporating cathode air-assisted heating and gas-ice phase change.The model’s accuracy was verified against experimental cold-start data from a stack composed of identical single cells.This computational framework enables quantitative analysis of temperature fields and ice fraction distributions across domains under varying air-assisted heating powers during cold starts.Building upon this model,the study further investigates the improvement in cold start performance by heating the cathode intake air with varying power levels.The results demonstrate that the fuel cell achieves self-startup at temperatures as low as−13℃ under a constant current density of 100mA/cm^(2) without air-assisted heating.At an ambient temperature of−20℃,a successful start-up can be achieved with a heating power of 0.45 W/cm^(2).The temperature variation overtime during the cold start process can be represented by a sum of two exponential functions.The air-assisted heating scheme proposed in this study has significantly improved the cold start performance of fuel cells in low-temperature environments.Additionally,it provides critical reference data and validation support for component selection and feasibility assessment of hybrid power systems.
基金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.
基金supported by the National High Technology Research and Development Program (863 Program)(2009AA06Z108)
文摘To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.
文摘This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
文摘In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The grid can be discretionarily refined as a non-structure triangle or quadrilateral so that piers can be treated as one or several impermeable elements with an area of 20 to 30 km^2 in a model domain over 85 700 km^2. The computational results of tidal levels and horizontal velocities are in good agreement with the field data. Based on the computed results by the model, the layout of an open 105 DWT liquefied natural gas （LNG）terminal in Caofeidian, Bohal Sea is effectively and reasonably optimized. It can be concluded that the model is suitable and reasonable for direct simulation of tidal currents around small structures in projects.
