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.

Investigations on High-Speed Flash Boiling Atomization of Fuel Based on Numerical Simulations

Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.

Numerical simulations for radon migration and exhalation behavior during measuring radon exhalation rate with closed-loop method

Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.

Numerical Simulations of the Flow Field around a Cylindrical Lightning Rod

As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.

Numerical simulations of earthquake rupture induced by pressure perturbation

The subsurface fluid injection can cause pressure increase within faults,leading to earthquake occurrences.However,the factors controlling earthquake rupture due to pressure perturbation remain poorly understood.To resolve this problem,we simulate the physical processes of earthquake nucleation and rupture on strike-slip faults perturbated by pressure migration based on the slip-weakening law.Multiple kinds of factors,including background stress,fluid injection rates,the area of the pressurized region,fault geometry,and fault friction coefficients,are considered in our simulations.Our simulation results reveal that the ratio of shear stress to normal stress rather than their absolute values controls the rupture behavior.With the large stress ratios,high injection rates,and large pressurized areas,earthquakes are prone to propagate as runaway ruptures.Additionally,faults with large aspect ratios of length to width are also favorable for causing runaway ruptures.In contrast,the factors of fault strike,dip angles and friction coefficients have minor influence on rupture behavior.

Optimal Design of High-Speed Partial Flow Pumps using Orthogonal Tests and Numerical Simulations

To investigate the influence of structural parameters on the performances and internal flow characteristics of partial flow pumps at a low specific speed of 10000 rpm,special attention was paid to the first and second stage impeller guide vanes.Moreover,the impeller blade outlet width,impeller inlet diameter,blade inclination angle,and number of blades were considered for orthogonal tests.Accordingly,nine groups of design solutions were formed,and then used as a basis for the execution of numerical simulations(CFD)aimed at obtaining the efficiency values and heads for each design solution group.The influence of impeller geometric parameters on the efficiency and head was explored,and the“weight”of each factor was obtained via a range analysis.Optimal structural parameters were finally chosen on the basis of the numerical simulation results,and the performances of the optimized model were verified accordingly(yet by means of CFD).Evidence is provided that the increase in the efficiency and head of the optimized model was 12.11%and 23.5 m,respectively,compared with those of the original model.

Vertical impedance functions of pile groups under low-to-high loading amplitudes:numerical simulations and experimental validation

Piles in a group experience additional displacements in soil due to pile-to-pile interactions apart from those resulting from the external loading.The effect of these interactions determined assuming soil as an elastic and/or viscoelastic material on pile head impedance functions of the pile group is studied by relating the group stiffness to the static stiffness of a single pile.However,the prevailing elastic solutions may misestimate the resulting pile group response due to the lack of consideration for either soil(material)and/or soil-pile interface nonlinearities.It is well established that soil behaves nonlinearly under moderate-to-high loading amplitudes,and besides,the soil-pile interface nonlinearity can exist even at small loading amplitudes.This study addresses the effects of these nonlinearities on the vertical impedance functions of a 3×3-pile group using numerical methods by direct analyses and superposition using pile-to-pile interaction factors.The numerical results are validated using scaled model tests under 1 g conditions.The results highlight the overestimation of pile-to-pile interactions in the pile group when assuming elastic soil conditions.The cases either by direct analyses or superposition approach involving soil and soil-pile interface nonlinearities agree well with the experimental pile group responses under close-to-elastic and nonlinear conditions.

Image Representations of Numerical Simulations for Training Neural Networks

A large amount of data can partly assure good fitting quality for the trained neural networks.When the quantity of experimental or on-site monitoring data is commonly insufficient and the quality is difficult to control in engineering practice,numerical simulations can provide a large amount of controlled high quality data.Once the neural networks are trained by such data,they can be used for predicting the properties/responses of the engineering objects instantly,saving the further computing efforts of simulation tools.Correspondingly,a strategy for efficiently transferring the input and output data used and obtained in numerical simulations to neural networks is desirable for engineers and programmers.In this work,we proposed a simple image representation strategy of numerical simulations,where the input and output data are all represented by images.The temporal and spatial information is kept and the data are greatly compressed.In addition,the results are readable for not only computers but also human resources.Some examples are given,indicating the effectiveness of the proposed strategy.

Numerical simulation on sand sedimentation and erosion characteristics around HDPE sheet sand barrier under different wind angles

For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In this study,using numerical simulations,we examined the behavior of High Density Polyethylene(HDPE)sheet sand barriers under different wind angles,focusing on flow field distribution,windproof efficiency,and sedimentation erosion dynamics.This study discovered that at a steady wind speed,airflow velocity varies as the angle between the airflow and the HDPE barrier changes.Specifically,a 90°angle results in the widest low-speed airflow area on the barrier’s downwind side.If the airflow is not perpendicular to the barrier,it prompts a lateral airflow movement which decreases as the angle expands.The windproof efficiency correlates directly with this angle but inversely with the wind’s speed.Notably,with a wind angle of 90°,wind speed drops by 81%.The minimum wind speed is found at 5.1H(the sand barrier height)on the barrier’s downwind side.As the angle grows,the barrier’s windproof efficiency improves,extending its protective reach.Sedimentation is most prominent on the barrier’s downwind side,as the wind angle shifts from 30°to 90°,the sand sedimentation area on the barrier’s downwind side enlarges by 14.8H.As the angle grows,sedimentation intensifies,eventually overtakes the forward erosion and enlarges the sedimentation area.

Elastic-viscoplastic behaviors of polymer-blend geocell sheets:Numerical and experimental investigations

Polymer-blend geocell sheets(PBGS)have been developed as substitute materials for manufacturing geocells.Various attempts have been made to test and predict the behaviors of commonly used geogrids,geotextiles,geomembranes,and geocells.However,the elastic-viscoplastic behaviors of novel-developed geocell sheets are still poorly understood.Therefore,this paper investigates the elastic-viscoplastic behaviors of PBGS to gain a comprehensive understanding of their mechanical properties.Furthermore,the tensile load-strain history under various loading conditions is simulated by numerical calculation for widespread utilization.To achieve this goal,monotonic loading tests,short-term creep and stress relaxation tests,and multi-load-path tests(also known as arbitrary loading history tests)are performed using a universal testing machine.The results are simulated using the nonlinear three-component(NLTC)model,which consists of three nonlinear components,i.e.a hypo-elastic component,a nonlinear inviscid component,and a nonlinear viscid component.The experimental and numerical results demonstrate that PBGS exhibit significant elastic-viscoplastic behavior that can be accurately predicted by the NLTC model.Moreover,the tensile strain rates significantly influence the tensile load,with higher strain rates resulting in increased tensile loads and more linear load-strain curves.Also,parametric analysis of the rheological characteristics reveals that the initial tensile strain rates have negligible impact on the results.The rate-sensitivity coefficient of PBGS is approximately 0.163,which falls within the typical range observed in most geosynthetics.

Analysis of debris flow control effect and hazard assessment in Xinqiao Gully,Wenchuan M_(s)8.0 earthquake area based on numerical simulation

Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the effectiveness of the debris flow control project and evaluated the debris flow hazards.Through field investigation and numerical simulation methods,the indicators of flow intensity reduction rate and storage capacity fullness were proposed to quantify the effectiveness of the engineering measures in the debris flow event.The simulation results show that the debris flow control project reduced the flow intensity by41.05%to 64.61%.The storage capacity of the dam decreases gradually from upstream to the mouth of the gully,thus effectively intercepting and controlling the debris flow.By evaluating the debris flow of different recurrence intervals,further measures are recommended for managing debris flow events.

Flow field, sedimentation, and erosion characteristics around folded linear HDPE sheet sand fence: Numerical simulation study

Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.

Magnetopause properties at the dusk magnetospheric flank from global magnetohydrodynamic simulations,the kinetic Vlasov equilibrium,and in situ observations--Potential implications for SMILE

We derived the properties of the terrestrial magnetopause(MP)from two modeling approaches,one global–fluid,the other local–kinetic,and compared the results with data collected in situ by the Magnetospheric Multiscale 2(MMS2)spacecraft.We used global magnetohydrodynamic(MHD)simulations of the Earth’s magnetosphere(publicly available from the NASA-CCMC[National Aeronautics and Space Administration–Community Coordinated Modeling Center])and local Vlasov equilibrium models(based on kinetic models for tangential discontinuities)to extract spatial profiles of the plasma and field variables at the Earth’s MP.The global MHD simulations used initial solar wind conditions extracted from the OMNI database at the time epoch when the MMS2 observes the MP.The kinetic Vlasov model used asymptotic boundary conditions derived from the same in situ MMS measurements upstream or downstream of the MP.The global MHD simulations provide a three-dimensional image of the magnetosphere at the time when the MMS2 crosses the MP.The Vlasov model provides a one-dimensional local view of the MP derived from first principles of kinetic theory.The MMS2 experimental data also serve as a reference for comparing and validating the numerical simulations and modeling.We found that the MP transition layer formed in global MHD simulations was generally localized closer to the Earth(roughly by one Earth radius)from the position of the real MP observed by the MMS.We also found that the global MHD simulations overestimated the thickness of the MP transition by one order of magnitude for three analyzed variables:magnetic field,density,and tangential speed.The MP thickness derived from the local Vlasov equilibrium was consistent with observations for all three of these variables.The overestimation of density in the Vlasov equilibrium was reduced compared with the global MHD solutions.We discuss our results in the context of future SMILE(Solar wind Magnetosphere Ionosphere Link Explorer)campaigns for observing the Earth’s MP.

Quasi-direct numerical simulations of the flow characteristics of a thermal plasma reactor with counterflow jet

Three-dimensional quasi-direct numerical simulations have been performed to investigate a thermal plasma reactor with a counterflow jet. The effects of the momentum flux ratio and distance between the counterflow jet and the thermal plasma jet on the flow characteristics are addressed. The numerical results show that the dimensionless location of the stagnation layer is significantly affected by the momentum flux ratio, but it is not dependent on the distance.Specifically, the stagnation layer is closer to the plasma torch outlet with the increase of the momentum flux ratio. Furthermore, the flow regimes of the stagnation layer and the flow characteristics of the thermal plasma jet are closely related to the momentum flux ratio. The characteristic frequencies associated with the different regimes are identified. The deflecting oscillation flow regimes are found when the momentum flux ratio is low, which provokes axial velocity fluctuations inside the thermal plasma jet. By contrast, for cases with a high momentum flux ratio, flapping flow regimes are distinguished. The thermal plasma jets are very stable and the axial velocity fluctuations mainly exist in the stagnation layer.

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.

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.

Desulfurization characteristics of slaked lime and regulation optimization of circulating fluidized bed flue gas desulfurization process--A combined experimental and numerical simulation study

Circulating fluidized bed flue gas desulfurization(CFB-FGD) process has been widely applied in recent years. However, high cost caused by the use of high-quality slaked lime and difficult operation due to the complex flow field are two issues which have received great attention. Accordingly, a laboratory-scale fluidized bed reactor was constructed to investigate the effects of physical properties and external conditions on desulfurization performance of slaked lime, and the conclusions were tried out in an industrial-scale CFB-FGD tower. After that, a numerical model of the tower was established based on computational particle fluid dynamics(CPFD) and two-film theory. After comparison and validation with actual operation data, the effects of operating parameters on gas-solid distribution and desulfurization characteristics were investigated. The results of experiments and industrial trials showed that the use of slaked lime with a calcium hydroxide content of approximately 80% and particle size greater than 40 μm could significantly reduce the cost of desulfurizer. Simulation results showed that the flow field in the desulfurization tower was skewed under the influence of circulating ash. We obtained optimal operating conditions of 7.5 kg·s^(-1)for the atomized water flow, 70 kg·s^(-1)for circulating ash flow, and 0.56 kg·s^(-1)for slaked lime flow, with desulfurization efficiency reaching 98.19% and the exit flue gas meeting the ultraclean emission and safety requirements. All parameters selected in the simulation were based on engineering examples and had certain application reference significance.

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.

Operation optimization of prefabricated light modular radiant heating system:Thermal resistance analysis and numerical study

The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.

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.