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.

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 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.

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.

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 SIMULATIONS OF TEMPERATURE FIELD IN DIRECT METAL LASER SINTERING PROCESS

A mathematical model is developed for simulating the heat transferring behavior in a direct metal laser sintering process. The model considers the thermal phenomena involved in the process, including conduction, radiation, and convection. A formula for the calculation of the heat conductivity of a sintering system containing solid phase, liquid phase, and gas phase is given. Due to the continuous movement of the laser beam, a local coordinate system centered on the laser beam is used to simplify the analytical calculation. Assuming that it is approximately a Gaussian laser beam, the heat conduction equation is resolved based on the assumption of the thermal insulating boundary conditions and the fixed thermal physical parameters. The FORTRAN language is employed to compile the program to simulate the temperature field in the direct copper powder sintering process. It shows a good agreement with the preliminary experimental results.[KH3/4D]

Numerical Simulations of Polymer Melt Conveying in Co-Rotating Twin Screw Extruder

The 3 D non isothermal flow of non Newtonian viscous polymer melt in a co rotating twin screw extruder is modeled. The distributions of the velocity, temperature, pressure and the viscous dissipation in the flow domain are presented by using a fluid dynamics analysis package (Polyflow). The numerical results show that the temperatures are high in the intermeshing region and on the screw surface, the maximum pressure and the minimum pressure occur in the intermeshing region, and the flow rate is almost proportional to the screw speed.

Some Analyses and Numerical Simulations of Meiyu in East Asia in 1983

In this paper the differences between Meiyu and Baiu front in 1983 have firstly been analysed, the trajectories of air on and to the north side of Meiyu and Baiu fronts during the Meiyu period have then been traced, and the forecasting and simulating of 4 sets of Meiyu onset of the year have finally been run utilizing the global model at UK Me-leorological Office. The results show: 1) Meiyu fronts are different from Baiu ones in temperature, humidity and stratification fields in lower atmosphere; and the possibly reasons for it are explained. 2) The Bay of Bengal is the main moisture source for Meiyu front, the South China Sea and the Pacific, for Baiu ones; and some existed arguments on it are also discussed. 3) The onset of Meiyu and its rainfall and rain belts are sensitive to the Tibetan Plateau, and the water vapour conditions over the Bay of Bengal and the South China Sea, but not sensitive to the SST over the equatorial area or to the East of Japan.

Numerical simulations of kinetic formation mechanism of Tangjiashan landslide

Tangjiashan landslide is a typical high-speed landslide hosted on consequent bedding rock. The landslide was induced by Wenchuan earthquake at a medium-steep hill slope. The occurrence of Tangjiashan landslide was basically controlled by the tectonic structure, topography, stratum lithology, slope structure, seismic waves, and strike of river. Among various factors, the seismic loading with great intensity and long duration was dominant. The landslide initiation exhibited the local amplification effect of seismic waves at the rear of the slope, the dislocation effect on the fault, and the shear failure differentiating effect on the regions between the soft and the hard layers. Based on field investigations and with the employment of the distinct element numerical simulation program UDEC （universal distinct element code）, the whole kinetic sliding process of Tan iashan landslide was represented and the formation mechanism of the consequent rock landslide under seismic loading was studied. The results are helpful for understanding seismic dynamic responses of consequent bedding rock slopes, where the slope stability could be governed by earthquakes.

STRESS-STRAIN HYSTERESIS OF Cu-Zn-Al POLYCRYSTALLINE SMAs─Ⅱ: EXPERIMENTS AND NUMERICAL SIMULATIONS

In this paper, the stress-strain bys torests at various temperatures was investigated experimentally for Cu-Zn-Al polycrystalline shape memory alloys (SMAs). Numerical simulations of the pseudoelastic hysteresis were performed based on a model, which had been proposed by the authors. As observed in the experiments, the shapes of the outer loops of the hysteresis varied strongly from specimen to specimen, which have the same chemical composition but different heat-treatment. Rather complicated inner hysteretic curves were obtained at testing temperatures higher than Af. The numerical simulations of the stress-strain hysteresis and the inner hysteretic curves agreed quite well with the experimental results.

EFFECTS OF PACIFIC SSTA ON SUMMER PRECIPITATION OVER EASTERN CHINA, PART II: NUMERICAL SIMULATIONS

Based on an observational analysis, seven numerical experiments are designed to study the impacts of Pacific SSTA on summer precipitation over eastern China and relevant physical mechanism by NCAR CCM3. The numerical simulation results show that preceding winter SSTA in the Kuroshio region leads to summer precipitation anomaly over the Yangtze River valleys by modifying atmospheric general circulation over eastern Asia and middle-high latitude. West Pacific subtropical high is notably affected by preceding spring SSTA over the middle and east of Equator Pacific; SSTA of the central region of middle latitude in the corresponding period causes the summer rainfall anomaly over eastern China so as to trigger the atmospheric Eurasia-Pacific teleconnection pattern.

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy,i.e.,they have a non-homogeneous mass distribu-tion along their symmetry axis.We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate.The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravi-tational torque stability number and by the particle aspect ratio.We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential ori-entation and tumbling rate of the particles.We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.

MULTI-GRANULE IMPACT MODEL OF NOPD AND ITS NUMERICAL SIMULATIONS

Aiming at the impaction among granules of non obstructive particle damping(NOPD), the vibration absorption model for vertical impact of granules is established by adopting Hertz contact theory. The numerical simulation of the granules movement process is proceeded, and the vibration response of a free free uniform beam is obtained for the case when all granules act on it. Through this method, the effect on vibration absorption of impaction is investigated. The simulational data show that multi gra nule vertical impaction is not sensitive to the movement clearance. The vibration absorption is also very well when the clearance changes within a large range. Therefore, the phenomenon that the vibration magnitude may increase if the clearance in a single impact body is improperly selected will not happen. The effect of vibration suppression in the range of middle and high frequencies(2 500~6 000 Hz) is better than that in the range of low frequency(<2 500 Hz). It indicates that the effect on vibration absorption of multi granule can well restrain the vibration of middle and high frequencies.

Parallel numerical simulations for quantized vortices in Bose-Einstein condensates

We employ the parallel computing technology to study numerically the three-dimensional structure of quantized vortices of Bose-Einstein condensates, For anisotropic cases, the bending process of vortices is described in detail by the decrease of Gross-Pitaevskii energy. A completely straight vortex and the steady and symmetrical multiple-vortex configurations are obtained. We analyse the effect of initial conditions and angular velocity on the number and shape of vortices.

Numerical Simulations of One-Directional Fractional Pharmacokinetics Model

In this paper,we present a three-compartment of pharmacokinetics model with irreversible rate constants.The compartment consists of arterial blood,tissues and venous blood.Fick’s principle and the law of mass action were used to develop the model based on the diffusion process.The model is modified into a fractional pharmacokinetics model with the sense of Caputo derivative.The existence and uniqueness of the model are investigated and the positivity of the model is established.The behaviour of the model is investigated by implementing numerical algorithms for the numerical solution of the system of fractional differential equations.MATLAB software is used to plot the graphs for illustrating the variation of drug concentration concerning time.Therefore,the numerical simulations of the model are presented for different values ofαwhich verified the theoretical analysis.Besides,we also observed the pattern of the simulations at the three-compartment of the model by using different values of initial conditions.

Identification of Forcing Mechanisms of Convective Initiation over Mountains through High-Resolution Numerical Simulations

Convection and its ensuing severe weather, such as heavy rainfall, hail, tornado, and high wind, have significant im- pacts on our society and economy （e.g., Cao et al., 2004; Fritsch and Carbone, 2004; Verbout et al., 2006; Ashley and Black, 2008; Cao, 2008; Cao and Ma, 2009; Zhang et al., 2014）. Due to its localized and transient nature, the initiation of convection or convective initiation remains one of the least understood aspects of convection in the scientific communi- ties, and it is a significant challenge to accurately predict the exact timing and location of convective initiation （e.g., Cai et al., 2006; Wilson and Roberts, 2006; Xue and Martin, 2006; Cao and Zhang, 2016）.

NUMERICAL SIMULATIONS OF β-GYRES IN TROPICAL CYCLONES

The circulation of β-gyres in tropical cyclones is studied using numerical simulations. As shown in the result, there is clear circulation of β-gyres in the deviation flow field of the middle layer of the model,i.e. there is cyclone current west of the vortex center but anticyclone current east of it. The theory analysis shows that the circulation of β-gyres is formed by the advection of geostrophic vorticity.

Numerical simulations of full-wave fi elds and analysis of channel wave characteristics in 3-D coal mine roadway models

Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.