Determination method of mesh size for numerical simulation of blast load in near-ground detonation

In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.

Full-Scale Numerical Simulation of the Local Scour Under Combined Current and Wave Conditions Based on Field Data

The monopile is the most common foundation to support offshore wind turbines.In the marine environment,local scour due to combined currents and waves is a significant issue that must be considered in the design of wind turbine foundations.In this paper,a full-scale numerical model was developed and validated based on field data from Rudong,China.The scour development around monopiles was investigated,and the effects of waves and the Reynolds number Re were analyzed.Several formulas for predicting the scour depth in the literature have been evaluated.It is found that waves can accelerate scour development even if the KC number is small(0.78<KC<1.57).The formula obtained from small-scale model tests may be unsafe or wasteful when it is applied in practical design due to the scale effect.A new equation for predicting the scour depth based on the average pile Reynolds number(Rea)is proposed and validated with field data.The equilibrium scour depth predicted using the proposed equation is evaluated and compared with those from nine equations in the literature.It is demonstrated that the values predicted from the proposed equation and from the S/M(Sheppard/Melville)equation are closer to the field data.

Three-Dimensional Direct Numerical Simulation of Flow past A Near-Wall Circular Cylinder:Combined Effects of Gap Ratio and Boundary Layer Thickness on Flow Profiles and Pressure Distribution

Three-dimensional direct numerical simulations of the wake flow downstream of a near-wall circular cylinder at different gap ratios and boundary layer thicknesses are carried out by using the iterative immersed boundary method.The non-dimensional gap between the cylinder and the wall,G/D=0.2,0.6 and 1.0,the non-dimensional boundary layer thickness,δ/D=0.0,0.7 and 1.6,the Reynolds number,Re=350,and the aspect ratio of the cylinder,L/D=25are adopted.High-resolution visualizations of the complex vortex structures at differentδ/D and G/D are presented.The transition of the streamwise vortex mode,the combined effects ofδ/D and G/D on the flow statistics,the pressure and shear stress distribution and the hydrodynamic forces are analyzed.Results show that with decreasing G/D and increasingδ/D,the gap flow and its vortex-shedding are significantly weakened,together with an elongated wake and an enlarged low-velocity area near the wall,leading to the wake mode transition from the two-sided to one-sided vortex-shedding.Different relative positions of the cylinder regarding the boundary layer alter the flow features of the shear layers.With an increase inδ/D,the front stagnation point shifts to the upper surface,and the distance between the flow divergence point and the maximum pressure position increases.The mean drag coefficient and r.m.s.values of drag and lift coefficients decrease with a decrease in G/D and an increase inδ/D,while the mean lift coefficient increases with decreasing G/D but decreases with increasingδ/D.

Numerical simulation on size effect of copper with nano-scale twins

To understand the tensile deformation of electro-deposited Cu with nano-scale twins, a numerical study was carried out based on a conventional theory of mechanism-based strain gradient plasticity （CMSG）. The concept of twin lamella strengthening zone was used in terms of the cohesive interface model to simulate grain-boundary sliding and separation. The model included a number of material parameters, such as grain size, elastic modulus, plastic strain hardening exponent, initial yield stress, as well as twin lamellar distribution, which may contribute to size effects of twin layers in Cu polycrystalline. The results provide information to understand the mechanical behaviors of Cu with nano-scale growth twins.

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

Calculation of effective temperature for pavement rutting using numerical simulation methods

In order to predict the long-term rutting of asphalt pavement, the effective temperature for pavement rutting is calculated using the numerical simulation method. The transient temperature field of asphalt pavement was simulated based on actual meteorological data of Nanjing. 24-hour rutting development under a transient temperature field was calculated in each month. The rutting depth accumulated under the static temperature field was also estimated and the relationship between constant temperature parameters was analyzed. Then the effective temperature for pavement rutting was determined based on the rutting equivalence principle. The results show that the monthly effective temperature is above 40 t in July and August, while in June and September it ranges from 30 to 40 Rutting development can be ignored when the monthly effective temperature is less than 30 t. The yearly effective temperature for rutting in Nanjing is around 38. 5 t. The long-term rutting prediction model based on the effective temperature can reflect the influences of meteorological factors and traffic time distribution.

Mesh Size Effect in Numerical Simulation of Blast Wave Propagation and Interaction with Structures

Numerical method is popular in analysing the blast wave propagation and interaction with structures.However,because of the extremely short duration of blast wave and energy trans-mission between different grids,the numerical results are sensitive to the finite element mesh size.Previous numerical simulations show that a mesh size acceptable to one blast scenario might not be proper for another case,even though the difference between the two scenarios is very small,indicating a simple numerical mesh size convergence test might not be enough to guarantee accu-rate numerical results.Therefore,both coarse mesh and fine mesh were used in different blast scenarios to investigate the mesh size effect on numerical results of blast wave propagation and interaction with structures.Based on the numerical results and their comparison with field test re-sults and the design charts in TM5-1300,a numerical modification method was proposed to correct the influence of the mesh size on the simulated results.It can be easily used to improve the accu-racy of the numerical results of blast wave propagation and blast loads on structures.

Numerical simulation on borehole breakout and borehole size effect using discrete element method

Estimation of horizontal stress magnitudes from borehole breakouts has been an attractive topic in the petroleum and mining industries,although there are critical research gaps that remain unfilled.In this paper,numerical simulation is conducted on Gosford sandstone to investigate the borehole breakout and its associated borehole size effect,including temperature influence.The discrete element method(DEM)model shows that the borehole breakout angular span is constant after the initial formation,whereas its depth propagates along the minimum horizontal stress direction.This indicates that the breakout angular span is a reliable parameter for horizontal stress estimation.The borehole size effect simulations illustrated the importance of borehole size on breakout geometries in which smaller borehole size leads to higher breakout initiation stress as well as the stress re-distribution from borehole wall outwards through micro-cracking.This implies that the stress may be averaged over a distance around the borehole and breakout initiation occurs at the borehole wall rather than some distance into the rock.In addition,the numerical simulation incorporated the thermal effect which is widely encountered in deep geothermal wells.Based on the results,the higher temperature led to lower breakout initiation stress with same borehole size,and more proportion of shear cracks was generated under higher temperature.This indicates that the temperature might contribute to the micro-fracturing mode and hence influences the horizontal stress estimation results from borehole breakout geometries.Numerical simulation showed that breakout shape and dimensions changed considerably under high stress and high temperature conditions,suggesting that the temperature may need to be considered for breakout stress analysis in deep locations.

AN ANALYSIS OF SUBGRID-RESOLVED SCALE INTERACTIONS WITH USE OF RESULTS FROM DIRECT NUMERICAL SIMULATIONS

Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.

Sheltering effect of punched steel plate sand fences for controlling blown sand hazards along the Golmud-Korla Railway:Field observation and numerical simulation studies

Sand fences made of punched steel plate(PSP)have recently been applied to control wind-blown sand in desertified and Gobi areas due to their strong wind resistance and convenient in situ construction.However,few studies have assessed the protective effect of PSP sand fences,especially through field observations.This study analyzes the effects of double-row PSP sand fences on wind and sand resistance using field observations and a computational fluid dynamics(CFD)numerical simulation.The results of field observations showed that the average windproof efficiencies of the first-row and second-row sand fences were 79.8%and 70.8%,respectively.Moreover,the average windproof efficiencies of the numerical simulation behind the first-row and second-row sand fences were 89.8%and 81.1%,respectively.The sand-resistance efficiency of the double-row PSP sand fences was 65.4%.Sand deposition occurred close to the first-row sand fence;however,there was relatively little sand on the leeward side of the second-row sand fence.The length of sand accumulation near PSP sand fences obtained by numerical simulation was basically consistent with that through field observations,indicating that field observations combined with numerical simulation can provide insight into the complex wind-blown sand field over PSP sand fences.This study indicates that the protection efficiency of the double-row PSP sand fences is sufficient for effective control of sand hazards associated with extremely strong wind in the Gobi areas.The output of this work is expected to improve the future application of PSP sand fences.

Numerical Simulation of the Marangoni Effect with Interphase Mass Transfer Between Two Planar Liquid Layers

The Marangoni effect induced by mass transfer at the interface between two immiscible liquids displays important influence on laboratory and industrial operation of solvent extraction. A systematic numerical study of the two-dimensional Marangoni effect in a two liquid layer system was conducted. The linear relationship of the inter- facial tension versus the solute concentration was incorporated into a mathematical model accounting for liquid flow and mass transfer in both phases. The typical cases analyzed by Sternling ＆ Scriven （AIChE J., 1959） using the linear instability theory were simulated bv the finite difference method and good agreement between the theory and the numerical simulation was observed. The simulation suggests that the Marangoni convection needs certain time to develop sufficiently in strength and scale to enhance the interphase mass transfer, the Marangoni effect is dynamic and transient, and remains at some stabilized level as long as the mass transfer driving force is kept con- stant. When certain level of shear is imposed at the interface as in most cases of practical significance, the Maran- goni effect is suppressed slightly but progressively as the shear is increased gradually. The present two-dimensional simulation of the Marangoni effect provides some insight into the underlying mechanism and also the basis for further theoretical study of the three-dimensional Marangoni effect in the real world and in chemical engineering applications.

Numerical simulation of mean interference effects for the interfering buildings with different heights

The mean wind-induced interference effects between two high-rise buildings,in which the interfering buildings have different heights,were numerical simulated in the terrain roughness of B and D types by the Reynolds stress equation model(RSM)of fluent.The results are in good agreement with those of the wind tunnel test.Influences of the relative arrangement of two buildings,the height of the interfering buildings and the terrain roughness upon the mean interference effects were analyzed,and the space distributions of IFCPs on the principal building under tandem arrangement were studied.The results indicate that the lower interfering buildings can always bring larger interference factors comparing to the higher ones under tandem arrangement except that the height is larger than 1.25h,and the heights' influence on the mean interference effects will increase as the reduced spacing of two buildings.The influence of heights will be little under stagger arrangement.

Dynamic characteristics of charging effects on the dielectric constant due to E-beam irradiation:a numerical simulation

A series of synthetic variations of material intrinsic properties always come with charging phenomena due to electron beam irradiation.The effects of charging on the dielectric constant will influence the charging dynamic in return.In this paper,we propose a numerical simulation for investigating the dynamic characteristics of charging effects on the dielectric constant due to electron beam irradiation.The scattering process between electrons and atoms is calculated considering elastic and inelastic collisions via the Rutherford model and the fast secondary electron model,respectively.Internal charge drift due to E-field,density gradient caused diffusion,charges trap by material defect,free electron and hole neutralization,and variation in the internal dielectric constant are considered when simulating the transport process.The dynamics of electron and hole distributions and charging states are demonstrated during E-beam irradiation.As a function of material nonlinear susceptibility and primary energy,the dynamics of charging states and dielectric constants are then presented in the charging process.It is found that the variation in the internal dielectric constant is more with respect to the depth and irradiation time.Material with a larger nonlinear susceptibility corresponds a faster charging enhancement.In addition,the effective dielectric constant and the surface potential have a linear relationship in the charging balance.Nevertheless,with shrinking charging affect range,the situation with a higher energy primary electron comes with less dielectric constant variation.The proposed numerical simulation mode of the charging process and the results presented in this study offer a comprehensive insight into the complicated charging phenomena in electron irradiation related fields.

NUMERICAL SIMULATION STUDY OF THE EFFECTS OF TERRAINAROUND TAIWAN ON TYPHOON DOT (1990)

With the use of analyzed TCM-90 data,Typhoon Dot (1990)and complex evolutions are successfullyreproduced in numerical sboulation as it travels over the island of Taiwan. The simulation includes theformation of secondary highs in both circulation and geopotential fields, trajectory jumps during mer-gence with main centers, anomalous northward movement of a low center west of Taiwan, evolutionalprocess of a low-level jet over the Taiwan Straits and significant deviations of circulation center fromgeopotentia1 center between upper and lower levels. By examining ev0lutions of 3-h interval simulatedresults, detailed processes of changes in structure and track before and after Dot’s passage of Taiwanis given, whose evolutional images are otherwise impossible with conventionai observations at intervalsof 6 h. A number of control experboents are conducted in the end of the work for understanding causesand mechanisms behind various properties.

Numerical Simulation of Microphysics in Meso-β-Scale Convective Cloud System Associated with a Mesoscale Convective Complex

Numerical simulation of meso-β-scale convective cloud systems associated with a PRE-STORM MCC case has been carried out using a 2-D version of the CSU Regional Atmospheric Modeling System (RAMS) nonhydrostatic model with parameterized microphysics. It is found that the predicted meso-r-scale convective phenomena arc basically unsteady under the situation of strong shear at low-levels, while the meso-β-scale convective system is maintained up to 3 hours or more. The meso -β- scale cloud system exhibits characteristics of a multi-celled convective storm in which the meso-r-scale convective cells have lifetime of about 30 min. Pressure perturbation depicts a meso-low after a half hour in the low levels. As the cloud system evolves, the meso-low intensifies and extends to the upshear side and covers the entire domain in the mid-lower levels with the peak values of 5-8 hPa. Temperature perturbation depicts a warm region in the middle levels through the entire simulation period. The meso-r-scale warm cores with peak values of 4-8 ℃ are associated with strong convective cells. The cloud top evaporation causes a stronger cold layer around the cloud top levels.Simulation of microphysics exhibits that graupel is primarily concentrated in the strong convective cells forming the main source of convective rainfall after one hour of simulation time. Aggregates are mainly located in the stratiform region and decaying convective cells which produce the stratiform rainfall. Riming of the ice crystals is the predominant precipitation formation mechanism in the convection region, whereas aggregation of ice crystals is the predominant one in the stratiform region, which is consistent with observations. Sensitivity experiments of ice-phase mierophysical processes show that the microphysical structures of the convective cloud system can be simulated better with the diagnosed aggregation collection efficiencies.

Numerical simulation of gap effect in supersonic flows

The gap effect is a key factor in the design of the heat sealing in super- sonic vehicles subjected to an aerodynamic heat load. Built on S-A turbulence model and Roe discrete format, the aerodynamic environment around a gap on the surface of a supersonic aircraft was simulated by the finite volume method. As the presented results indicate, the gap effect depends not only on the attack angle, but also on the Mach number.

Numerical simulation of displacement characteristics of CO_2 injected in pore-scale porous media

Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous media, affecting displacement results. Direct observation of the flow patterns in the porous media is difficult, and therefore knowledge about the two-phase displacement flow is insufficient. In this paper, a two-dimensional(2D) pore structure was extracted from a sandstone sample, and the flow process that CO_2 displaces resident brine in the extracted pore structure was simulated using the Navier eStokes equation combined with the conservative level set method. The simulation results reveal that the pore throat is a crucial factor for determining CO_2 displacement process in the porous media. The two-phase meniscuses in each pore throat were in a self-adjusting process. In the displacement process,CO_2 preferentially broke through the maximum pore throat. Before breaking through the maximum pore throat, the pressure of CO_2 continually increased, and the curvature and position of two-phase interfaces in the other pore throats adjusted accordingly. Once the maximum pore throat was broken through by the CO_2, the capillary force in the other pore throats released accordingly; subsequently, the interfaces withdrew under the effect of capillary fore, preparing for breaking through the next pore throat.Therefore, the two-phase displacement in CO_2 injection is accompanied by the breaking through and adjusting of the two-phase interfaces.

An Application of the RAMS/FLUENT System on the Multi-Scale Numerical Simulation of the Urban Surface Layer—A Preliminary Study

The Regional Atmospheric Modeling System （RAMS） and the computational fluid dynamics （CFD） codes known as FLUENT are combinatorially applied in a multi-scale numerical simulation of the urban surface layer （USL）. RAMS and FLUENT are combined as a multi-scale numerical modeling system, in which the RAMS simulated data are delivered to the computational model for FLUENT simulation in an offline way. Numerical simulations are performed to present and preliminarily validate the capability of the multi-scale modeling system, and the results show that the modeling system can reasonably provide information on the meteorological elements in an urban area from the urban scale to the city-block scale, especially the details of the turbulent flows within the USL.

Simulation of large-scale numerical substructure in real-time dynamic hybrid testing

A solution scheme is proposed in this paper for an existing RTDHT system to simulate large-scale finite element （FE） numerical substructures. The analysis of the FE numerical substructure is split into response analysis and signal generation tasks, and executed in two different target computers in real-time. One target computer implements the response analysis task, wherein a large time-step is used to solve the FE substructure, and another target computer implements the signal generation task, wherein an interpolation program is used to generate control signals in a small time-step to meet the input demand of the controller. By using this strategy, the scale of the FE numerical substructure simulation may be increased significantly. The proposed scheme is initially verified by two FE numerical substructure models with 98 and 1240 degrees of freedom （DOFs）. Thereafter, RTDHTs of a single frame-foundation structure are implemented where the foundation, considered as the numerical substructure, is simulated by the FE model with 1240 DOFs. Good agreements between the results of the RTDHT and those from the FE analysis in ABAQUS are obtained.