Protective performance of shear stiffening gel-modified foam against ballistic impact:Experimental and numerical study

As one of the most widely used personal protective equipment(PPE),body armors play an important role in protecting the human body from the high-velocity impact of bullets or projectiles.The body torso and critical organs of the wear may suffer severe behind-armor blunt trauma(BABT)even though the impactor is stopped by the body armor.A type of novel composite material through incorporating shear stiffening gel(STG)into ethylene-vinyl acetate(EVA)foam is developed and used as buffer layers to reduce BABT.In this paper,the protective performance of body armors composed of fabric bulletproof layers and a buffer layer made of foam material is investigated both experimentally and numerically.The effectiveness of STG-modified EVA in damage relief is verified by ballistic tests.In parallel with the experimental study,numerical simulations are conducted by LS-DYNA®to investigate the dynamic response of each component and capture the key mechanical parameters,which are hardly obtained from field tests.To fully describe the material behavior under the transient impact,the selected constitutive models take the failure and strain rate effect into consideration.A good agreement between the experimental observations and numerical results is achieved to prove the validity of the modelling method.The tests and simulations show that the impact-induced deformation on the human body is significantly reduced by using STG-modified EVA as the buffering material.The improvement of protective performance is attributed to better dynamic properties and more outstanding energy absorption capability of the composite foam.

Numerical Analysis of Perforation during Hydraulic Fracture Initiation Based on Continuous-Discontinuous Element Method

Perforation is a pivotal technique employed to establish main flow channels within the reservoir formation at the outset of hydraulic fracturing operations.Optimizing perforation designs is critical for augmenting the efficacy of hydraulic fracturing and boosting oil or gas production.In this study,we employ a hybrid finite-discrete element method,known as the continuous–discontinuous element method(CDEM),to simulate the initiation of post-perforation hydraulic fractures and to derive enhanced design parameters.The model incorporates the four most prevalent perforation geometries,as delineated in an engineering technical report.Real-world perforations deviate from the ideal cylindrical shape,exhibiting variable cross-sectional profiles that typically manifest as an initial constriction followed by an expansion,a feature consistent across all four perforation types.Our simulations take into account variations in perforation hole geometries,cross-sectional diameters,and perforation lengths.The findings show that perforations generated by the 39g DP3 HMX perforating bullet yield the lowest breakdown pressure,which inversely correlates with increases in sectional diameter and perforation length.Moreover,this study reveals the relationship between breakdown pressure and fracture degree,providing valuable insights for engineers and designers to refine perforation strategies.

Geometric Constraints for Global Regularity of 3D Shear Thickening Fluids

We consider the relation between the direction of the vorticity and the global regularity of 3D shear thickening fluids.It is showed that a weak solution to the non-Newtonian incompressible fluid in the whole space is strong if the direction of the vorticity is (11-5p)/2-Holder continuous with respect to the space variables when 2<P<11/5.

Motion stability of high-speed maglev systems in consideration of aerodynamic effects: a study of a single magnetic suspension system

In this study, the intrinsic mechanism of aerodynamic effects on the motion stability of a high-speed maglev system was investigated. The concept of a critical speed for maglev vehicles considering the aerodynamic effect is proposed. The study was carried out based on a single magnetic suspension system, which is convenient for proposing relevant concepts and obtaining explicit expressions. This study shows that the motion stability of the suspension system is closely related to the vehicle speed when aerodynamic effects are considered. With increases of the vehicle speed, the stability behavior of the system changes. At a certain vehicle speed,the stability of the system reaches a critical state, followed by instability. The speed corresponding to the critical state is the critical speed. Analysis reveals that when the system reaches the critical state, it takes two forms, with two critical speeds, and thus two expressions for the critical speed are obtained. The conditions of the existence of the critical speed were determined, and the effects of the control parameters and the lift coefficient on the critical speed were analyzed by numerical analysis. The results show that the first critical speed appears when the aerodynamic force is upward,and the second critical speed appears when the aerodynamic force is downward. Moreover, both critical speeds decrease with the increase of the lift coefficient.

Effect of perforation on flow past a conic cylinder at Re=100:vortex-shedding pattern and force history

The flow past a circular-section cylinder with a conic shroud perforated with four holes at the peak was simulated numerically at Re = 100, considering two factors, viz.the angle of attack and the diameter of the holes. The effects of the perforated conic shroud on the vortex shedding pattern in the near wake was mainly investigated, as well as the time history of the drag and lift forces. In the investigated parameter space, three flow regimes were generally identified, corresponding to weak, moderate, and strong disturbance effects.In regime I, the wake can mainly be described by alternately shedding Kármán or Kármán-like vortices. In regime II, the spanwise vortices are obviously disturbed along the span due to the appearance of additional vorticity components and their interactions with the spanwise vortices, but still shed in synchronization along the spanwise direction. In regime III,the typical Kármán vortices partially or totally disappear,and some new vortex shedding patterns appear, such as-type, obliquely shedding, and crossed spanwise vortices with opposite sign. Corresponding to these complex vortex shedding patterns in the near wake, the fluid forces no longer oscillate regularly at a single vortex shedding frequency, but rather with a lower modulation frequency and multiple amplitudes. An overview of these flow regimes is presented.

Influence of original microstructure on the transformation behavior and mechanical properties of ultra-high-strength TRIP-aided steel

The transformation behavior and tensile properties of an ultra-high-strength transformation-induced plasticity （TRIP） steel （0.2C-2.0Si-I.SMn） were investigated by different heat treatments for automobile applications. The results show that F-TRIP steel, a tradi- tional TRIP steel containing as-cold-rolled ferfite and pearlite as the original microstructure, consists of equiaxed grains of intercritical ferrite surrounded by discrete particles of M/RA and B. In contrast, M-TRIP steel, a modified TRiP-aided steel with martensite as the original mi- crostlucture, containing full martensite as the original microstructure is comprised of lath-shaped grains of ferrite separated by lath-shaped martensite/retained austenite and bainite. Most of the austenite in F-TRIP steel is granular, while the austenite in M-TRIP steel is lath-shaped. The volume fraction of the retained austenite as well as its carbon content is lower in F-TRIP steel than in M-TRIP steel, and austenite grains in M-TRIP steel are much finer than those in F-TRIP steel. Therefore, M-TRIP steel was concluded to have a higher austenite stability, re- sulting in a lower transformation rate and consequently contributing to a higher elongation compared to F-TRIP steel. Work hardening be- havior is also discussed for both types of steel.

Effect of perforation on flow past a conic cylinder at Re =100: wavy vortex and sign laws

In order to find the intrinsic physical mechanism of the original Kármán vortex wavily distorted across the span due to the introduction of three-dimensional (3-D) geometric disturbances,a flow past a peak-perforated conic shroud is numerically simulated at a Reynolds number of 100.Based on previous work by Meiburg and Lasheras (1988),the streamwise and vertical interactions with spanwise vortices are introduced and analyzed.Then vortex-shedding patterns in the near wake for different flow regimes are reinspected and illustrated from the view of these two interactions.Generally,in regime Ⅰ,spanwise vortices are a little distorted due to the weak interaction.Then in regime Ⅱ,spanwise vortices,even though curved obviously,are still shed synchronously with moderate streamwise and vertical interactions.But in regime Ⅲ,violently wavy spanwise vortices in some vortex-shedding patterns,typically an Ω-type vortex,are mainly attributed to the strong vertical interactions,while other cases,such as multiple vortex-shedding patterns in sub-regime Ⅲ-D,are resulted from complex streamwise and vertical interactions.A special phenomenon,spacial distribution of streamwise and vertical components of vorticity with specific signs in the near wake,is analyzed based on two models of streamwise and vertical vortices in explaining physical reasons of top and bottom shear layers wavily varied across the span.Then these two models and above two interactions are unified.Finally two sign laws are summarized:the first sign law for streamwise and vertical components of vorticity is positive in the upper shear layer,but negative in the lower shear layer,while the second sign law for three vorticity components is always negative in the wake.

Effect of Peak Perforation on Flow Past A Conic Cylinder at Re=100:Drag,Lift and Strouhal Number

A flow past a circular-section cylinder with a perforated conic shroud, in which the perforation is located at the peak of the conic disturbance as the shroud installed on the cylinder and uniformly distributed with several circular holes, is numerically simulated at a Reynolds number of 100. Two factors in the perforation are taken into account, i.e. the attack angle relative to the direction of incoming flow and diameter of holes. The effect of such perforation on the drag, lift and vortex-shedding frequency is mainly investigated. Results have shown that variation of the attack angle has a little effect, especially on the drag and vortex-shedding frequency, except in certain cases due to the varied vortex-shedding patterns in the near wake. The increasing hole diameter still exhibits a little effect on the drag and frequency of vortex shedding, but really reduces the lift, in particular at larger wavelength, such as the lift reduction reaching almost 66% 68% after introducing the perforation.

Seepage-stress coupled modeling for rainfall induced loess landslide

Although rainfall is rare on the Loess Plateau of western China, landslides occur frequently there in rainy season. Surveys report that landslide hazards always follow heavy rains. In this study, a seepage-stress coupling model for rainfall induced landslide is used to examine an actual disastrous event in Yulin by the end of July, 2017. The effects of rainfall duration, rainfall intensity and soil weakening on slope stability are studied in detail. The results illustrate that the safety factor drops sharply at first and then is gradually declining to below 1.05 during additional two days of heavy rain. With soil strength softening considered, the slope would be more unstable, in which the weakening in soil cohesion is found to be a more sensitive factor.

Modified asymptotic Adomian decomposition method for solving Boussinesq equation of groundwater flow

The Adomian decomposition method （ADM） is an approximate analytic method for solving nonlinear equations. Generally, an approximate solution can be ob- tained by using only a few terms. However, in applications, we need to use it flexibly according to the real problem. In this paper, based on the ADM, we give a modified asymptotic Adomian decomposition method and use it to solve the nonlinear Boussinesq equation describing groundwater flows. The example shows effectiveness of the modified asymptotic Adomian decomposition method.

Physical mechanism for origin of streamwise vortices in mode A of a square-section cylinder

In the present paper, physical mechanism responsible for origin of streamwise vortices in mode A appeared in the threedimensional (3-D) wake transition of a square-section cylinder is investigated. Direct numerical simulations at a Reynolds number of 180 firstly show that such streamwise vorticity is not originated from lateral surfaces. Then through the analysis of local flow field in the immediate neighborhood of rear surface, based on the theory of vortex-induced vortex, a new physical mechanism is identified. At first, the vertical vorticity on rear surface is generated by the intrinsic three-dimensional instability with the same instability wavelength of mode A. Then the streamwise vorticity at a specific sign is induced by such vertical vorticity, convected and concentrated in the shear layers. Finally, streamwise vortices are formed and shed with alternatively shedding spanwise vortices in the near wake. Moreover, the effect of induced spanwise vorticity on original two-dimensional (2-D) spanwise vorticity is also presented in detail.

Frame-invariance in finite element formulations of geometrically exact rods

This article is concerned with finite element implementations of the three- dimensional geometrically exact rod. The special attention is paid to identifying the con- dition that ensures the frame invariance of the resulting discrete approximations. From the perspective of symmetry, this requirement is equivalent to the commutativity of the employed interpolation operator I with the action of the special Euclidean group SE（3）, or I is SE（3）-equivariant. This geometric criterion helps to clarify several subtle issues about the interpolation of finite rotation. It leads us to reexamine the finite element for- mulation first proposed by Simo in his work on energy-momentum conserving algorithms. That formulation is often mistakenly regarded as non-objective. However, we show that the obtained approximation is invariant under the superposed rigid body motions, and as a corollary, the objectivity of the continuum model is preserved. The key of this proof comes from the observation that since the numerical quadrature is used to compute the integrals, by storing the rotation field and its derivative at the Gauss points, the equiv- ariant conditions can be relaxed only at these points. Several numerical examples are presented to confirm the theoretical results and demonstrate the performance of this al- gorithm.

A modified criterion for wave-induced momentary liquefaction of sandy seabed

The assessment of the wave-induced soil liquefaction plays a key role in the geotechnical design for offshore foundations. The underlying shortcomings of the existing momentary liquefaction criteria are identified and clarified by mechanism analyses and the recent field observations. A modified criterion for the wave-induced momentary liquefaction of a sandy seabed is given to describe the vertical pore- pressure distributions. An improved approximation of the momentary liquefaction depth is further presented. Parametric study of the effects of the saturation degree of soils indicates that this modification is significant for the evaluation of wave-induced momentary liauefaction.

Experiment and numerical simulation on the characteristics of fluid–structure interactions of non-rigid airships

Fluid-structure interaction is an important issue for non-rigid airships with inflated envelopes. In this study, a wind tunnel test is conducted, and a loosely coupled procedure is correspondingly established for numerical simulation based on computational fluid dynamics and nonlinear finite element analysis methods. The typical results of the numerical simulation and wind tunnel experiment, including the overall lift and deformation, are in good agreement with each other. The results obtained indicate that the effect of fluid-structure interaction is noticeable and should be considered for non-rigid airships. Flow- induced deformation can further intensify the upward lift force and pitching moment, which can lead to a large deformation. Under a wind speed of 15 m/s, the lift force of the non-rigid model is increased to approximatelv 60% compared with that of the rigid model under a high angle of attack.

The effect of gravity on self-similarity of Worthington jet after water entry of a two-dimensional wedge

The effect of gravity on the self-similarity of jet shape at late stage of Worthington jet development is investigated by experiment in the study.In addition,the particle image velocimetry(PIV)method is introduced to analyze the development of flow field.There is a linear scaling regarding the axial velocity of the jet and the scaling coefficient increases with the Froude number.

Numerical Simulation of Solitary Wave Forces on A Vertical Cylinder on A Slope Beach

Wave forces acting on a vertical cylinder at different locations on a slope beach in the near-shore region are investigated considering solitary waves as incoming waves.Based on the Reynolds-averaged Navier-Stokes equations and the k-ε turbulence model,wave forces due to the interaction between the solitary wave and cylinder are simulated and analyzed with different incident wave heights and cylinder locations.The numerical results are first compared with previous theoretical and experimental results to validate the model accuracy.Then,the wave forces and characteristics around the cylinder are studied,including the velocity field,wave surface elevation and pressure.The effects of relative wave height,Keulegan-Carpenter(KC)number and cylinder locations on the wave forces are also discussed.The results show that the wave forces exerted on a cylinder exponentially increase with the increasing incident wave height and KC number.Before the wave force peaks,the growth rate of the wave force shows an increasing trend as the cylinder moves onshore.The cylinder location has a notable effect on the wave force on the cylinder in the near-shore region.As the cylinder moves onshore,the wave force on the cylinder initially increases and then decreases.For the cases considered here,the maximum wave force appears when the cylinder is located one cylinder diameter below the still-water shoreline.Furthermore,the fluid velocity peaks when the maximum wave force appears at the same location.

Analytical solutions for sediment concentration in waves based onlinear diffusivity

Two kinds of analytical solutions are derived through Laplace transform for the equation that governswave-induced suspended sediment concentration with linear sediment diffusivity under two kinds ofbottom boundary conditions,namely the reference concentration(Dirichlet)and pickup function(Nu-mann),based on a variable transformation that is worked out to transform the governing equation intoa modified Bessel equation.The ability of the two analytical solutions to describe the profiles of sus-pended sediment concentration is discussed by comparing with different experimental data.And it isdemonstrated that the two analytical solutions can well describe the process of wave-induced suspendedsediment concentration,including the amplitude and phase and vertical profile of sediment concentra-tion.Furthermore,the solution with boundary condition of pickup function provides better results thanthat of reference concentration in terms of the phase-dependent variation of concentration.

Triangular temporal-distribution law for disintegrating internal solitons over a step

Internal solitary waves have been found to disintegrate into a series of solitons over variablebathymetry, with important applications for offshore engineering. Considering realisticbackground stratification in the South China Sea, internal solitary waves propagating over a stepare studied here. By assuming disintegrated solitons propagate independently, a theoreticalmodel, namely a triangular temporal-distribution law based on the Korteweg–de Vries theory, isproposed to describe the fission process of internal solitary waves undergoing disintegration. Aparameter is then introduced to quantify the accuracy of the theoretical model. The resultsindicate that the triangular law predicts the fission process better for a longer travelling distanceand a larger amplitude of internal solitary waves.

Linear analysis of the dynamic response of a riser subject to internal solitary waves

The flow field induced by internal solitary waves(ISWs)is peculiar wherein water motion occurs in the whole water depth,and the strong shear near the pycnocline can be generated due to the opposite flow direction between the upper and lower layers,which is a potential threat to marine risers.In this paper,the flow field of ISWs is obtained with the Korteweg-de Vries(Kd V)equation for a two-layer fluid system.Then,a linear analysis is performed for the dynamic response of a riser with its two ends simply supported under the action of ISWs.The explicit expressions of the deflection and the moment of the riser are deduced based on the modal superposition method.The applicable conditions of the theoretical expressions are discussed.Through comparisons with the finite element simulations for nonlinear dynamic responses,it is proved that the theoretical expressions can roughly reveal the nonlinear dynamic response of risers under ISWs when the approximation for the linear analysis is relaxed to some extent.

Numerical study on wake characteristics of high-speed trains

Intensive turbulence exists in the wakes of high speed trains, and the aerodynamic performance of the trailing car could deteriorate rapidly due to complicated features of the vortices in the wake zone. As a result, the safety and amenity of high speed trains would face a great challenge. This paper considers mainly the mechanism of vortex formation and evolution in the train flow field. A real CRH2 model is studied, with a leading car, a middle car and a trailing car included. Different running speeds and cross wind conditions are considered, and the approaches of un- steady Reynold-averaged Navier-Stokes （URANS） and de- tached eddy simulation （DES） are utilized, respectively. Re- suits reveal that DES has better capability of capturing small eddies compared to URANS. However, for large eddies, the effects of two approaches are almost the same. In conditions without cross winds, two large vortex streets stretch from the train nose and interact strongly with each other in the wake zone. With the reinforcement of the ground, a complicated wake vortex system generates and becomes strengthened as the running speed increases. However, the locations of flow separations on the train surface and the separation mechanism keep unchanged. In conditions with cross winds, three large vortices develop along the leeward side of the train, among which the weakest one has no obvious influence on the wake flow while the other two stretch to the tail of the train and combine with the helical vortices in the train wake. Thus, optimization of the aerodynamic performance of the trailing car should be aiming at reducing the intensity of the wake vortex system.