Explosion resistance performance of reinforced concrete box girder coated with polyurea:Model test and numerical simulation

To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.

High-Precision Flow Numerical Simulation and Productivity Evaluation of Shale Oil Considering Stress Sensitivity

Continental shale oil reservoirs,characterized by numerous bedding planes and micro-nano scale pores,feature significantly higher stress sensitivity compared to other types of reservoirs.However,research on suitable stress sensitivity characterization models is still limited.In this study,three commonly used stress sensitivity models for shale oil reservoirs were considered,and experiments on representative core samples were conducted.By fitting and comparing the data,the“exponential model”was identified as a characterization model that accurately represents stress sensitivity in continental shale oil reservoirs.To validate the accuracy of the model,a two-phase seepage mathematical model for shale oil reservoirs coupled with the exponential model was introduced.The model was discretely solved using the finite volume method,and its accuracy was verified through the commercial simulator CMG.The study evaluated the productivity of a typical horizontal well under different engineering,geological,and fracture conditions.The results indicate that considering stress sensitivity leads to a 13.57%reduction in production for the same matrix permeability.Additionally,as the fracture half-length and the number of fractures increase,and the bottomhole flowing pressure decreases,the reservoir stress sensitivity becomes higher.

Image processing based three-dimensional model reconstruction for cross-platform numerical simulation

Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.

The numerical simulation for a 3D two-phase anisotropic medium based on BISQ model

Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.

Numerical simulation of residual stress field in green power metallurgy compacts by modified Drucker-Prager Cap model

Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.

Experimental research and numerical simulation of the multi-field performance of cemented paste backfill:Review and future perspectives

Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-field performance interaction.At present,research on the multi-field performance of CPB mainly includes indoor similar simulation experiments,in-situ multi-field performance monitoring experiments,multi-field performance coupling model construction of CPB,and numerical simulation of the multi-field performance of CPB.Because it is hard to study the in-situ multi-field performance of CPB in the real stope,most current research on in-situ multi-field performance adopts the numerical simulation method.By simulating the conditions of CPB in the real stope(e.g.,maintenance environment,stope geometry,drainage conditions,and barricade and backfilling rates),the multi-field performance of CPB is further studied.This paper summarizes the mathematical models employed in the numerical simulation and lists the engineering application cases of numerical simulation in the in-situ multi-field performance of CPB.Finally,it proposes that the multi-field performance of CPB needs to strengthen the theoretical study of multi-field performance,form the strength design criterion based on the multi-field performance of CPB,perform a full-range numerical simulation of the multi-field performance of CPB,develop a pre-warning technology for the CPB safety of CPB,develop automatic and wireless sensors for the multi-field performance monitoring of CPB,and realize the application and popularization of CPB monitoring technology.

Numerical simulation of the flow within and over an intersection model with Reynolds-averaged Navier-Stokes method

In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group （RNG） model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.

2-D Composite Model for Numerical Simulations of Nonlinear Waves

A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.

Numerical Simulation of the Whole Three-Dimensional Flow in a Stirred Tank with Anisotropic Algebraic Stress Model

In accordance to the anisotropic feature of turbulent flow, ananisotropic algebraic stress model is adopted to predict theturbulent flow field and turbulent characteristics generated by aRushton disc turbine with the improved inner-outer iterativeprocedure. The predicted turbulent flow is compared with experimentaldata and the simulation by the standard k-ε turbulence model. Theanisotropic algebraic stress model is found to give better predictionthan the standard k-ε turbulence model. The predicted turbulent flowfield is in accordance to experimental data and the trend of theturbulence intensity can be effectively reflected in the simulation.

Optimization of flow control devices for a T-type five-strand billet caster tundish: water modeling and numerical simulation

The optimization of flow control devices(FCDs) for a T-type five-strand billet caster tundish was carried out by water modeling and numerical simulation. In water modeling experiments, flow characteristics of the bare tundish and tundish conf igurations with designed U-type baff les and a round turbulence inhibitor were analyzed using residence time distribution(RTD) curves. Mathematical models for liquid steel in the real plant tundish were established using the fluid dynamics software package Fluent. The flow field, the temperature field, and the RTD curves of liquid steel in the proposed tundish conf igurations were obtained. The results of numerical simulation and water modeling were validated with each other by the predicted and experimental RTD curves. The results of flow field and temperature field were used to ref lect the actual state of a real plant tundish and to choose the optimal FCD. Finaly, from the whole performance of the multi-strand tundish, the optimal scheme was determined by combining the results of water modeling and numerical simulation. With the optimal tundish equipped with U-type baffle with def lector holes and round turbulence inhibitor, not only was the flow characteristic of each strand improved, but also the difference of flow characteristics between multiple strands was smaller.

Numerical Simulations of Intense Meiyu Rainfall in 1991 over the Changjiang and Huaihe River Valleys by a Regional Climate Model with p-σ Incorporated Coordinate System

Based on the primitive equation model with p- σ incorporated coordinate system originally developed by Qian et al., a one-way nested fine mesh limited area model is developed. This model is nested with ECMWF T42 data to simulate the extra-intensive rainfall event occurring in the Changjiang and Huaihe River valleys in summer of 1991. The results show that the model has certain capacity to fairly reproduce the regional distribution and the movement of the main rainfall belts. Therefore it can be used as a regional climate model to simulate and predict the short-range regional climate changes.

Numerical Simulation of Viscous Flow Through Spherical Particle Assemblage with the Modified Cell Model

The cell model developed since 1950s is a useful tool forexploring the behavior of particle assemblages, but it demandsfurther careful development of the outer boundary conditions so thatinteraction in a particle swarm is better represented. In this paper,the cell model and its development were reviewed, and themodifications of outer cell boundary conditions were suggested. Atthe cell outer boundary, the restriction of uniform liquid flow wasremoved in our simulation conducted in the reference frame fixed withthe particle.

Numerical simulation for macrosegregation in direct-chill casting of 2024 aluminum alloy with an extended continuum mixture model

An extended continuum mixture model for macrosegregation is applied to predicting Cu and Mg segregation in large-size ingot of 2024 aluminum alloy during direct chill casting（DC）. A microsegregation model using the approximate phase diagram data was coupled with macroscopic transport equations for macrosegregation profiles. Then, the impacts of transport mechanisms on the formation of macrosegregation were discussed. It is found that copper and magnesium have a similar segregation configuration from the billet center to surface. Negative segregation is observed in the centerline and subsurface, whereas positive segregation is obtained in the surface and somewhat underestimated positive segregation in the middle radius. Further, the discrepancy between the predicted and experimental results was discussed in detail. The results show that the magnesium to some extent alleviates the copper segregation in ternary alloy, compared with that in binary alloy. The predicted results show good agreement with measured experimental data obtained from literatures.

Numerical Simulation of Shrinkage Cavity Formation in Spheroidal Graphite Iron Castings Based on Micro Modeling

A micro-modeling method (MM) for the quantitative prediction of the shrinkage cavity formation in SGiron castings is proposed. The mathematical models describing the volume changes during the solidification ofspheroidal graphite cast iron are established based on the models of solidification kinetics. The shrinkage cavityformation of T-shaped SG iron castings is calculated with MM method. The calculated results are compared with theexperimental results. It is shown that the predicted size, shape and distribution of shrinkage cavity by MM methodare in good agreement with the measured results.

Multiple car-following model of traffic flow and numerical simulation

On the basis of the full velocity difference （FVD） model, an improved multiple car-following （MCF） model is proposed by taking into account multiple information inputs from preceding vehicles. The linear stability condition of the model is obtained by using the linear stability theory. Through nonlinear analysis, a modified Korteweg-de Vries equation is constructed and solved. The traffic jam can thus be described by the klnk-antikink soliton solution for the mKdV equation. The improvement of this new model over the previous ones lies in the fact that it not only theoretically retains many strong points of the previous ones, but also performs more realistically than others in the dynamical evolution of congestion. Furthermore, numerical simulation of traffic dynamics shows that the proposed model can avoid the disadvantage of negative velocity that occurs at small sensitivity coefficients λ in the FVD model by adjusting the information on the multiple leading vehicles. No collision occurs and no unrealistic deceleration appears in the improved model.

NUMERICAL SIMULATION OF METHANE-AIR TURBULENT JET FLAME USING A NEW SECOND-ORDER MOMENT MODEL

A new second-order moment model for turbulent combustion is applied in the simulation of methane-air turbulent jet flame. The predicted results are compared with the experimental results and with those predicted using the well-known EBU-Arrhenius model and the original second-order moment model. The comparison shows the advantage of the new model that it requires almost the same computational storage and time as that of the original second-order moment model, but its modeling results are in better agreement with experiments than those using other models. Hence, the new second-order moment model is promising in modeling turbulent combustion with NOx formation with finite reaction rate for engineering application.

A COUPLED MODEL FOR MUILTIPHASE FLUID FLOW AND SEDIMENTATION DEFORMATION IN OIL RESERVOIR AND ITS NUMERICAL SIMULATION

A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.

NUMERICAL SIMULATION OF WET STEAM CONDENSING FLOW WITH AN EULERIAN/EULERIAN MODEL

Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which the wet steam is regarded as mixture comprisingtwo coupled systems: the vapor phase and the liquid phase. These two systems are both described byconservation equations. High resolution TVD scheme is employed to capture condensing phenomena inwet steam flow. This model has been validated by numerical simulations of condensing flows in 1D and2D nozzles. Compared with experimental data, a good agreement is observed. This Eulerian/Eulerianmodel can be extended to 3D calculation of condensing flow.

Mathematical model for the thermal process of controlled cooling of wires and its numerical simulation

The mathematical model for the thermal process of billets rolling has beenestablished, including transporting in air and temperature-holding cover, descaling withhigh-pressure water, and the process of rolling and cooling in water box. The calculated data by themodel have been compared with the measured data and the results show that the model is right andcreditable. Based on the model, the main thermal characters of rolling line have been simulated andthe influence of all the parameters on the temperature of rolling has been analyzed.

Estimating the Operation Status of Steam Cracking Furnace Using Numerical Simulation with Combustion Models

An accurate and complete geometric model was constructed to simulate the combustion, flow and temperature environment in the radiant section of the steam cracking furnace. Simulation of flow and radiation status has utilized the standard k-ε model and P1 model. The finite-rate/eddy-dissipation (finite-rate/ED) combustion model and non-premixed combustion model were both used to simulate accurately the combustion and the operation status of the steam cracking furnace. Three different surfaces of the steam cracking furnace were obtained from the simulation, namely:the flue gas temperature field of the entrance surface in long flame burners, the central surface location of tubes, and the crossover section surface. Detailed information on the flue gas temperature and the mass concentration fraction of these different surfaces in the steam cracking furnace can also be obtained by the simulation. This paper analyzed and compared the simulation results with the two combustion models, estimated the operation status of the steam cracking furnace, and reported that the finite-rate/ED model is appropriate to simulate the steam cracking furnace by comparing key simulation data with actual test data. This work has also provided a theoretical basis for simulating and operating the steam cracking furnace.