Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

Recent progresses in the development of tubular segmented-in-series solid oxide fuel cells:Experimental and numerical study

Solid oxide fuel cells(SOFCs)have attracted a great deal of interest because they have the highest efficiency without using any noble metal as catalysts among all the fuel cell technologies.However,traditional SOFCs suffer from having a higher volume,current leakage,complex connections,and difficulty in gas sealing.To solve these problems,Rolls-Royce has fabricated a simple design by stacking cells in series on an insulating porous support,resulting in the tubular segmented-in-series solid oxide fuel cells(SIS-SOFCs),which achieved higher output voltage.This work systematically reviews recent advances in the structures,preparation methods,perform-ances,and stability of tubular SIS-SOFCs in experimental and numerical studies.Finally,the challenges and future development of tubular SIS-SOFCs are also discussed.The findings of this work can help guide the direction and inspire innovation of future development in this field.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

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.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Numerical study on local scour characteristics around submarine pipelines in the Yellow River Delta silty sandy soil under waves and currents

Due to their high reliability and cost-efficiency,submarine pipelines are widely used in offshore oil and gas resource engineering.Due to the interaction of waves,currents,seabed,and pipeline structures,the soil around submarine pipelines is prone to local scour,severely affecting their operational safety.With the Yellow River Delta as the research area and based on the renormalized group(RNG)k-εturbulence model and Stokes fifth-order wave theory,this study solves the Navier-Stokes(N-S)equation using the finite difference method.The volume of fluid(VOF)method is used to describe the fluid-free surface,and a threedimensional numerical model of currents and waves-submarine pipeline-silty sandy seabed is established.The rationality of the numerical model is verified using a self-built waveflow flume.On this basis,in this study,the local scour development and characteristics of submarine pipelines in the Yellow River Delta silty sandy seabed in the prototype environment are explored and the influence of the presence of pipelines on hydrodynamic features such as surrounding flow field,shear stress,and turbulence intensity is analyzed.The results indicate that(1)local scour around submarine pipelines can be divided into three stages:rapid scour,slow scour,and stable scour.The maximum scour depth occurs directly below the pipeline,and the shape of the scour pits is asymmetric.(2)As the water depth decreases and the pipeline suspension height increases,the scour becomes more intense.(3)When currents go through a pipeline,a clear stagnation point is formed in front of the pipeline,and the flow velocity is positively correlated with the depth of scour.This study can provide a valuable reference for the protection of submarine pipelines in this area.

Numerical simulation of microwave-induced cracking and melting of granite based on mineral microscopic models

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.

Numerical analysis of geosynthetic-reinforced embankment performance under moving loads

The performance of geosynthetic-reinforced embankments under traffic moving loads is always a hotspot in the geotechnical engineering field.A three-dimensional(3D)model of a geosynthetic-reinforced embankment without drainage consolidation was established using the finite element software ABAQUS.In this model,the traffic loads were simulated by two moving loads of rectangular pattern,and their amplitude,range,and moving speed were realized by a Fortran subroutine.The embankment fill was simulated by an equivalent linear viscoelastic model,which can reflect its viscoelasticity.The geogrid was simulated by the truss element,and the geocell was simulated by the membrane element.Infinite elements were utilized to weaken the boundary effect caused by the model geometry at the boundaries.Validation of the established numerical model was conducted by comparing the predicted deformations in the cross-section of the geosynthetic-reinforced embankment with those from the existing literature.On this basis,the dynamic stress and strain distribution in the pavement structure layer of the geosynthetic-reinforced embankment under a moving load was also analyzed.Finally,a parametric study was conducted to examine the influences of the different types of reinforcement,overload,and the moving load velocity on the geosynthetic-reinforced embankment.

Numerical simulation on the multiphase flow and reoxidation of the molten steel in a two-strand tundish during ladle change

A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.

Multi-physical fields distribution in billet during helical electromagnetic stirring:A numerical simulation research

Electromagnetic stirring is one of the widely applied techniques to modify the quality of casting billets.Different from conventional rotate stirring,the helical stirring is more professional in assisting multi-dimensional flow of molten metal and eliminating solidification defects.In this study,the single-winding helical stirring(SWHS)was introduced,offering advantages such as smaller volume and lower electromagnetic shielding compared to traditional helical stirring methods.Following a comprehensive numerical simulation,the stirring parameters of SWHS were adjusted to yoke inclination angle of 43°and frequency of 12 Hz.The higher electromagnetic force and flow velocity in drawing direction,as well as the lower temperature gradient induced by the SWHS,are positive factors for homogeneous solidification of billet.The experimental results on Al-8%Si alloy and 0.4%C-1.1%Mn steel demonstrate that compared to rotate stirring,the SWHS process can induce better billet quality and is more effective in accelerating the equiaxed expansion and reducing element segregation.The SWHS process can enhance the equiaxed ratio of the billet by 58.3%and reduce segregation degree of carbon element by 10.97%.Consequently,SWHS holds great promise as a potential approach for improving the quality of continuous casting billets.

Numerical investigation on the aerodynamic drag reduction based on bottom deflectors and streamlined bogies of a high-speed train

The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In this study,computational fluid dynamics(CFD)based on three-dimensional steady incompressible Reynolds-average Naiver-Stokes(RANS)equations and Realizable k-ε turbulence model were utilized for numerical simulations.Inspired by the concept of streamlined design and the idea of bottom flow field control,this study iteratively designed the bogies in a streamlined shape and combined them with the bottom deflectors to investigate the joint drag reduction mechanism.Three models,i.e.,single-bogie model,simplified train model,and eight-car high-speed train model,were created and their aerodynamic characteristics were analyzed.The results show that the single-bogie model with streamlined design shows a noticeable drag reduction,whose power bogie and trailer bogie experience 13.92%and 7.63%drag reduction,respectively.The range of positive pressure area on the bogie is reduced.The aerodynamic drag can be further reduced to 15.01%by installing both the streamlined bogie and the deflector on the simplified train model.When the streamlined bogies and deflectors are used on the eight-car model together,the total drag reduction rate reaches 2.90%.Therefore,the proposed aerodynamic kit for the high-speed train bottom is capable to improve the flow structure around the bogie regions,reduce the bottom flow velocity,and narrow the scope of the train’s influence on the surrounding environment,achieving the appreciable reduction of aerodynamic drag.This paper can provide a new idea for the drag reduction of high-speed trains.

In-situ measurement via the flow-through method and numerical simulations for radon exhalation during measurements of the radon exhalation rate

Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.

Numerical investigation of friction-heating-pressurization and its control parameters in the shear band of high-speed landslides

High-speed sliding often leads to catastrophic landslides,many of which,in the initial sliding phase before disintegration,experience a friction-induced thermal pressurization effect in the bottom shear band,accelerating the movement of the overlying sliding mass.To quantitatively investigate this complex multiphysical phenomenon,we established a set of equations that describe the variations in temperature and excess pore pressure within the shear band,as well as the conservation of momentum equation for the overlying sliding mass.With a simplified landslide model,we investigated the variations of temperature and excess pore pressure within the shear band and their impacts on the velocity of the overlying sliding mass.On this basis,we studied the impact of seven key parameters on the maximum temperature and excess pore pressure in the shear band,as well as the impact on the velocity of the overlying sliding mass.The simulation results of the standard model show that the temperature and excess pore pressure in the shear band are significantly higher than those in the adjacent areas,and reach the maximum values in the center.Within a few seconds after the start,the maximum excess pore pressure in the shear zone is close to the initial stress,and the shear strength loss rate exceeds 90%.The thermal pressurization mechanism significantly increases the velocity of the overlying sliding mass.The results of parameter sensitivity analysis show that the thermal expansion coefficient has the most significant impact on the temperature and excess pore pressure in the shear band,and the sliding surface dip angle has the most significant impact on the velocity of the overlying sliding mass.The results of this study are of great significance for clarifying the mechanism of thermal pressurization-induced high-speed sliding.

Numerical Simulation and Entropy Production Analysis of Centrifugal Pump with Various Viscosity

The fluid’s viscosity significantly affects the performance of a centrifugal pump.The entropy production method and leakage are employed to analyze the performance changes under various viscosities by numerical simulation and validated by experiments.The results showed that increasing viscosity reduces both the pump head and efficiency.In addition,the optimal operating point shifts to the left.Leakage is influenced by vortex distribution in the front chamber and boundary layer thickness in wear-ring clearance,leading to an initial increase and subsequent decrease in leakage with increasing viscosity.The total entropy production Spro,Total inside the pump rises with increasing viscosity.The different mechanisms dominate under varying conditions:Turbulent dissipation dominates at low viscosity.Under high-viscosity conditions,energy loss is primarily caused by direct dissipation Spro,D and wall entropy production Spro,W.This study provides a deeper and more objective understanding of the energy characteristics of centrifugal pumps handling fluids of various viscosity,potentially aiding in optimizing pump design and improving energy conversion efficiency.

Experimental and numerical study of buckling-restrained braces configured with out-of-plane eccentricity under cyclic loading

This study focuses on variations in the hysteretic behavior of buckling-restrained braces(BRBs)configured with or without out-of-plane eccentricity under cyclic loading.Quasi-static experiments and numerical simulations were carried out on concentrically and eccentrically loaded BRB specimens to investigate the mechanical properties,energy dissipation performance,stress distribution,and high-order deformation pattern.The experimental and numerical results showed that compared to the concentrically loaded BRBs,the stiffness,yield force,cumulated plastic ductility(CPD)coefficient,equivalent viscous damping coefficient and energy dissipation decreased,and the yield displacement and compression strength adjustment factor increased for the eccentrically loaded BRBs.With the existence of the out-of-plane eccentricity,the initial yield position changes from the yield segment to the junction between the yield segment and transition segment under a tensile load,while the initial high-order buckling pattern changes from a first-order C-shape to a secondorder S-shape under a compressive load.

Progressive fragmentation of granular assemblies within rockslides: Insights from discrete-continuous numerical modeling

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.

Numerical study on the cavity dynamics for vertical water entries of twin spheres

In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbulence structure is described using the shear-stress transport k-ω(SST k-ω) model, and the volume of fluid(VOF) method is used to track the complex air-liquid interface. The motion of spheres during water entry is simulated using an independent overset grid. The numerical model is verified by comparing the cavity evolution results from simulations and experiments. Numerical results reveal that the time interval between the twin water entries evidently affects cavity expansion and contraction behaviors in the radial direction. However, this influence is significantly weakened by increasing the lateral distance between the two spheres. In synchronous water entries, pressure is reduced on the midline of two cavities during surface closure, which is directly related to the cavity volume. The evolution of vortexes inside the two cavities is analyzed using a velocity vector field, which is affected by the lateral distance and time interval of water entries.

Resistance of full-scale beams against close-in explosions.Numerical modeling and field tests

This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.

Numerical simulation of melt flow and temperature field during DC casting 2024 aluminium alloy under different casting conditions

Casting speed,casting temperature and secondary cooling water flow rate are the main process parameters affecting the DC casting process.These parameters significantly influence the flow and temperature fields during casting,which are crucial for the quality of the ingot and can determine the success or failure of the casting operation.Numerical simulation,with the advantages of low cost,rapid execution,and visualized results,is an important method to study and optimize the DC casting process.In the present work,a simulation model of DC casting 2024 aluminum alloy was established,and the reliability of the model was verified.Then,the influence of casting parameters on flow field and temperature field was studied in detail by numerical simulation method.Results show that with the increase of casting speed,the melt flow becomes faster,the depths of slurry zone and mushy zone increase,and the variation of slurry zone depth is greater than that of mushy zone.With an increase in casting temperature,the melt flow rate increases,the depth of the slurry zone becomes shallower,and the depth of the mushy zone experiences only minor changes.The simulation results further indicate that the increase of the flow rate of the secondary cooling water slightly reduces the depths of both slurry and mushy zone.

Numerical simulation of flow field deposition and erosion characteristics around bridge-road transition section

Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flow field around roadbeds, overlooking detailed examinations of sand particle erosion and deposition patterns near bridges and roadbeds. This study employs numerical simulations to analyze the influence of varying heights and wind speeds on sand deposition and erosion characteristics at different locations: the bridge-road transition section(side piers), middle piers, and roadbeds. The results show that the side piers, experience greater accumulation than the middle piers. Similarly, the leeward side of the roadbed witnesses more deposition compared to the windward side. Another finding reveals a reduced sand deposition length as the vertical profile, in alignment with the wind direction, moves further from the bridge abutments at the same clearance height. As wind speeds rise, there’s a decline in sand deposition and a marked increase in erosion around the side piers, middle piers and roadbeds. In conclusion, a bridge clearance that’s too low can cause intense sand damage near the side piers, while an extremely high roadbed may lead to extensive surface sand deposition. Hence, railway bridges in areas prone to sandy winds should strike a balance in clearance height. This research provides valuable guidelines for determining the most suitable bridge and roadbed heights in regions affected by wind and sand.

Numerical Study on the Effect of Gap Diffraction on the Hydrodynamic Performance of A Floating Breakwater

Two-dimensional(2D)flume experiments are useful in investigating the performances of floating breakwaters(FBs),including hydrodynamic performances,motion responses,and mooring forces.Designing a reasonable gap between the flume wall and the FBs is a critical step in 2D flume tests.However,research on the effect of the gap on the accuracy of 2D FB experimental results is scarce.To address this issue,a numerical wave tank is developed using CFD to estimate the wave-FB interaction of a moored dual-cylindrical FB,and the results are compared to experimental data from a previously published work.There is good agreement between them,indicating that the numerical model is sufficiently accurate.The numerical model is then applied to explore the effect of gap diffraction on the performance of FBs in2D experiments.It was discovered that the nondimensional gap length L_(Gap)/W_(Pool)should be smaller than 7.5%to ensure that the relative error of the transmission coefficient is smaller than 3%.The influence of the gap is also related to the entering wave properties,such as the wave height and period.