Numerical investigation of friction-heating-pressurization and its control parameters in the shear band of high-speed landslides

High-speed sliding often leads to catastrophic landslides,many of which,in the initial sliding phase before disintegration,experience a friction-induced thermal pressurization effect in the bottom shear band,accelerating the movement of the overlying sliding mass.To quantitatively investigate this complex multiphysical phenomenon,we established a set of equations that describe the variations in temperature and excess pore pressure within the shear band,as well as the conservation of momentum equation for the overlying sliding mass.With a simplified landslide model,we investigated the variations of temperature and excess pore pressure within the shear band and their impacts on the velocity of the overlying sliding mass.On this basis,we studied the impact of seven key parameters on the maximum temperature and excess pore pressure in the shear band,as well as the impact on the velocity of the overlying sliding mass.The simulation results of the standard model show that the temperature and excess pore pressure in the shear band are significantly higher than those in the adjacent areas,and reach the maximum values in the center.Within a few seconds after the start,the maximum excess pore pressure in the shear zone is close to the initial stress,and the shear strength loss rate exceeds 90%.The thermal pressurization mechanism significantly increases the velocity of the overlying sliding mass.The results of parameter sensitivity analysis show that the thermal expansion coefficient has the most significant impact on the temperature and excess pore pressure in the shear band,and the sliding surface dip angle has the most significant impact on the velocity of the overlying sliding mass.The results of this study are of great significance for clarifying the mechanism of thermal pressurization-induced high-speed sliding.

Numerical investigation of hydrodynamic tractor-retarder assembly under traction work condition

To obtain the performance characteristics of hydrodynamic tractor-retarder assembly under traction work condition,a numerical simulation model of flow channel was established and tetrahedron unstructured grids were adopted in the meshing stage.The racing rotating speed of the brake wheels was calculated by computational fluid dynamics（CFD） calculation and interpolation,and then accurate boundary condition was applied to the CFD simulation to study the pressure and velocity distribution of internal flow field in hydrodynamic tractor-retarder assembly.Finally,the original characteristics were calculated by CFD post-processing analysis.Comparison of experimental data and flow field analysis results showed that the calculation tolerance of the torque ratio K and the efficiency η was less than 5%,and the calculation tolerance of the pump torque coefficient λ was less than 10%.

Numerical Investigation of the Performance of an Axial-Flow Pump with Tandem Blades

The performance characteristics of an axial-flow pump with tandem blades are studied based on the numerical computations. The arrangement of the pump impellers is established through the analysis of velocity triangles. With the commercial computational fluid dynamics (CFD) software NUMECA, the turbulent flow in the tandem axial-flow pump is simulated in various flow conditions. The detail flow structure in the leading edge region of the rear impeller is described, and the influence of the deflection angle of the rear blade on the head performance is studied. According to the simulation, the performance comparison is made between the tandem axial-flow pump and the conventional two-stage axial-flow pump with a uniform impeller size. Results of the study indicate that the tandem axial-flow pump can work in a wider range with high efficiency.

Numerical investigations on effects of bluff body in flat plate micro thermo photovoltaic combustor with sudden expansion

In order to reveal combustion characteristics of H_2/air mixture in a micro-combustor with and without bluff body, the effects of inlet velocities, equivalence ratios and bluff body's blockage ratios on the temperature field, pressure of the combustor wall, combustion efficiency and blow-off limit were investigated. The numerical results indicate that the sudden expansion plate micro combustor with bluff body could enhance the turbulent disturbance of the mixed gas in the combustion chamber and the combustion condition is improved. Moreover, a low-speed and high temperature recirculation region was formed between the sudden expansion step and the bluff body so that the high and uniform wall temperature(>1000 K) could be gotten. As a result, it could strengthen the mixing process, prolong the residence time of gas, control the flame position effectively and widen the operation range by the synergistic effect of the bluff body and steps. When the blockage ratio ranged from 0.3 to 0.6, it could be found that the bluff body could play a stabilizing effect and expand combustion blow burning limit, and combustion efficiency firstly was increased with the inlet velocity and equivalence ratio, and then was decreased.

Numerical Investigation on the Propagation Mechanism of Steady Cellular Detonations in Curved Channels

The propagation mechanism of steady cellular detonations in curved channels is investigated numerically with a detailed chemical reaction mechanism, The numerical results demonstrate that as the radius of the curvature decreases, detonation fails near the inner wall due to the strong expansion effect. As the radius of the curvature increases, the detonation front near the inner wall can sustain an underdriven detonation. In the case where deto- nation fails, a transverse detonation downstream forms and re-initiates the quenched detonation as it propagates toward the inner wall. Two kinds of propagation modes exist as the detonation is propagating in the curved channel. One is that the detonation fails first, and then a following transverse detonation initiates the quenched detonation and this process repeats itself. The other one is that without detonation failure and re-initiation, a steady detonation exists which consists of an underdriven detonation front near the inner wall subject to the diffraction and an overdriven detonation near the outer wall subject to the compression.

Numerical Investigation on the Flow and Temperature Fields in an Inductively Coupled Plasma Reactor

This paper gives a numerical study on the flow and temperature fields in an induced plasma reactor, which worked in 0.5 ATM with air as a working gas. We employed a two-dimensional mode of an inductively coupled plasma to calculate the temperature and flow field of the reactor as well as the generator. The algorithm is based on the solutions of the two-dimensional continuity, momentum, and energy equations in term of vorticity, stream function and enthalpy. An upwind finite-difference scheme was adopted to solve those equations with appropriate boundary conditions. The computed results show that there is a flat region with little parameter change in the reactor, that the diameter of the region is not much larger than that of the generator and that a deep change of parameter exists in the outer side of the region.

NUMERICAL INVESTIGATION OF THE MECHANISM OF ASYMMETRIC VORTEX FLOWS

A robust iterative method suitable for the numerical simulation of high angle-of-attack vortex flows is established based upon the multiple line-vortex model(MLVM).With symmetric or asymmetric positions of sep- aration lines given,the first converged solution at an angle of attack as high as 60 degree is obtained by means of the present method.Numerical experiments for a tangent-ogive forebody indicate the viscous onset mechanism of asymmetric vortex flows over a body of revolution at high angles of attack and zero sideslip.

NUMERICAL INVESTIGATION OF THE LIMIT LOADS FOR PRESSURE VESSELS WITH PART-THROUGH SLOTS

A numerical investigation of the limit loads is carried out for pressure vessels with part-through slots using a general computational method for the limit analysis of 3-D structures. The limit pressures are given for a comprehensive range of geometric parameters. Some of the calculated results are compared with the results of 3-D elastic-plastic finite element analysis and existing numerical solutions. The effects of various shapes and sizes of part-through slots on the load carrying capacity of cylindrical shells are investigated and evaluated. Two kinds of typical failure modes corresponding to different dimensions of slots are studied. Based on the numerical results, a geometric parameter G which combines the slot dimensions and the cylinder geometry is presented. It reasonably reflects the overall effect of slots on the limit loads of cylinders. An empirical formula for estimating the limit pressures of cylindrical shells with part-through slots is obtained.

Hydrodynamics numerical investigation of hoistable masts for underwater vehicles

Using the unsteady incompressible Navier-Stokes equation as the governing equation, the large eddy simulation (LES) model is implemented to investigate the shedding of vortices, the flow pattern of turbulence, the unsteady pressure fluctuation and the time history of the lift coefficient and drag coefficient of hoistable masts with various mast shapes and various arrangements in this paper. Combining the FFT, combined time-frequency transform and wavelet power spectrum analysis, the characteristics of unsteady pressure can be obtained in both time and frequency domain. It shows that the main frequency of pressure fluctuation is near the frequency of vortex shedding in time domain using the FFT method. It can be inferred from the combined time-frequency transform that the unsteady pressure fluctuation has obviously the peak value and the second peak value in time domain. It could indicate that the fluctuation power varies from the fluctuation frequency through the power spectrum analysis. By the data analysis, it shows that the vortex shedding is the dominant cause of the periodically pressure fluctuation. And the interaction pattern of wake and interplay between wake and the walls of masts under different arrangements are also discussed in this paper.

Numerical Investigations of Laminar Air Flow and Heat Transfer Characteristics in a Square Channel Inserted with Discrete X-V Baffles (XVB)

Thermal performance enhancement in a square channel heat exchanger(HX)using a passive technique is presented.Vortex turbulator insertion in a square channel HX as a passive technique is selected for thermal improvement.The vortex turbulator of interest is discrete X-V baffles(XVB).The discrete XVBs are inserted in the square channel with the main aim of generating vortex flow.The vortex flow generated can support the enhanced convective heat transfer coefficient and also enhance HX performance.Effects of baffle configuration(type A and B),baffle size(w/H=0.05,0.10,0.15 and 0.20),baffle distance(e/H=1,1.5 and 2)and flow direction(±x air flow paths)on fluid flow and thermal topologies are numerically investigated by using a commercial code.As shown by the numerical results,the predicted flow configuration with the discrete XVB insertions,which include impinging and vortex streams,is found through the HX channel.The perturbing thermal boundary layer and greater air blending are also found through the HX channel inserted with the discrete XVB.These mechanisms promote and augment the convection heat transfer coefficient,heat transfer rate and rise thermal potentiality.The maximum Nusselt number of the channel with the baffles inserted is 11.01 times upper than that of the smooth channel,while the greatest thermal performance factor(TPF)is observed to be around 3.45.

NUMERICAL INVESTIGATION OF AXIAL FLOW COMPRESSOR CASING TREATMENT

Using an inviscid model with inlet total pressure gradient and a J. D. Denton scheme, this paper for the first time numerically solves the 3D flow field of compressor casing treatment, and also explores some boundary singularities and numerical stability. Agreement is attained in qualitative explanations of some, casing treatment test results and its mechanism.

NUMERICAL INVESTIGATION OF THREE-DIMENSIONAL VISCOUS INCOMPRESSIBLE FLOWS IN DIVERGENT CURVED CHANNELS AND TURBULENT MODEL STUDY

In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.

Numerical Investigation on Flow and Heat Transfer Performance of Supercritical Carbon Dioxide Based on Variable Turbulent Prandtlnumber Model

Flow and heat transfer characteristic of supercritical carbon dioxide(SCO_(2))are numerically investigated in the horizontal and vertical tubes.TWL turbulent Prandtl number model could well describe the behavior of SCO_(2) affected by the buoyancy.Under the cooling condition,the heat transfer performance of SCO_(2) along the upward direction is best and that along the downward direction is worst when bulk fluid temperatures are below the pseudocritical temperature.Reducing the ratio of heat flux to mass flux could decrease the difference of convective heat transfer coefficient in three flow directions.Under the heating condition,heat transfer deterioration only occurs in vertical upward and horizontal flow directions.Heat transfer deterioration of SCO_(2) could be delayed by increasing the mass flux and the deterioration degree is weakened in the second half of tube along the vertical upward flow direction.Compared with the straight tube,the corrugated tube shows better comprehensive thermal performance.

Numerical investigation of the concentric annulus flow around a cylindrical body with contrasted effecting factors

This paper studies the wall-bounded flow around a cylindrical at a high Reynolds numbers body in a determined computational domain, with simulations of the 3-D, turbulent concentric annulus flow in a straight pipe. Numerical results show that a reversing zone, appearing as a tongue zone with nested velocities higher than the surrounding area, exists behind the cylindrical body. The annulus space is a region of high velocity and low pressure. The zero velocity, of combined the X- velocity and the Y- velocity, exists in the cross sections and no vortex shedding is formed behind the attaching cylinders. Among all investigated effecting factors, the diameters of the attaching and the main cylinders affect the wake feature behind the cylindrical body while the main cylinder length does not affect the distribution tendency of the flow field. The diameters of the main cylinder and the pipe affect the pressure values and the distribution tendencies on the main cylinder surface. Obviously, the increase of the pipe diameter reduces the drag coefficient of the cylindrical body and the increase of the diameter of the main cylinder increases the drag coefficient greatly. The numerical investigation of the concentric annulus flow provides foundations for further improvements of the intricate flow studies.

Numerical Investigation of a High-subsonic Axial-flow Compressor Rotor with Non-axisymmetric Hub Endwall

The major source of loss in modem compressors is the secondary loss. Non-axisymmetric endwall profile contouring is now a well established design methodology in axial flow turbines. However, flow development in axial compressors is differ from turbines, the effects of non-axisymmetric endwall to axial compressors requires flow analysis in detail. This paper presents both experimental and numerical data to deal with the application of a non-axisymmetric hub endwall in a high-subsonic axial-flow compressor. The aims of the experiment here were to make sure the numerically obtained flow fields is the physical mechanism responsible for the improvement in efficiency, due to the non-axisymmetric hub endwall. The computational results were first compared with avail- able measured data of axisymmetric hub endwall. The results agreed well with the experimental data for estima- tion of the global performance. The coupled flow of the compressor rotor with non-axisymmetric hub endwall was simulated by a state-of-the-art multi-block flow solver. The non-axisymmetric hub endwall was designed for a subsonic compressor rotor with the help of sine and cosine functions. This type of non-axisymmetric hub end- wall was found to have a significant improvement in efficiency of 0.45% approximately and a slightly increase for the total pressure ratio. The fundamental mechanisms of non-axisymmetric hub endwall and their effects on the subsonic axial-flow compressor endwall flow field were analyzed in detail. It is concluded that the non-axisymmetric endwall profiling, though not optimum, can mitigate the secondary flow in the vicinity of the hub endwall, resulting in the improvement of aerodynamic performance of the compressor rotor.

Using air curtains to reduce short-range infection risk in consulting ward: A numerical investigation

Air curtains is promising in reducing the short-range infection risk in hospitals.To quantitatively evaluate its performance,this paper explores air curtains equipped on normal consulting desk to avoid doctor’s direct exposure to the patient exhaled pollutants.A numerical investigation is conducted to evaluate the effects of supply air velocity and angle on cutting off performance.Simulation results show that the average mass fraction of exhaled pollutants decreases significantly(70%–90%)in the consulting ward,indicating satisfying performance of air curtains.Increasing supply air velocity is demonstrated to be conducive in forming full air curtains,whereas an excessively high supply air velocity may be of adverse effects by entraining exhaled flow.Besides,the supply air angle is also critical due to its coupling with supply air velocity.It is found that larger angle(0°–40°)is better where velocity is less than 3 m/s,otherwise a small angle(20°)is preferable where velocity is larger than 3 m/s.Exhaled flow could be well suppressed at the supply air angle 20°but moves over air curtains at 40°.This study can provide effective and intuitive guidance in applying air curtains in consulting wards.

Numerical investigations of an optical switch based on a silicon stripe waveguide embedded with vanadium dioxide layers

A novel scheme for the design of an ultra-compact and high-performance optical switch is proposed and investigated numerically. Based on a standard silicon(Si) photonic stripe waveguide, a section of hyperbolic metamaterials(HMM) consisting of 20-pair alternating vanadium dioxide (VO_2)∕Si thin layers is inserted to realize the switching of fundamental TE mode propagation. Finite-element-method simulation results show that, with the help of an HMM with a size of 400 nm × 220 nm × 200 nm(width × height × length), the ON/OFF switching for fundamental TE mode propagation in an Si waveguide can be characterized by modulation depth(MD) of5.6 d B and insertion loss(IL) of 1.25 dB. It also allows for a relatively wide operating bandwidth of 215 nm maintaining MD > 5 dB and IL < 1.25 dB. Furthermore, we discuss that the tungsten-doped VO_2 layers could be useful for reducing metal-insulator-transition temperature and thus improving switching performance. In general, our findings may provide some useful ideas for optical switch design and application in an on-chip all-optical communication system with a demanding integration level.

Further numerical investigation on concrete dynamic behaviors with considering stress non-equilibrium in SHPB test based on the waveform features

In this study,with the meso-scale model reliably validated in our previous work(Construction and Building Materials,2018),the waveform features of plain concrete under various loading conditions and especially with considering stress non-equilibrium are reliably reproduced and predicted.Associating with waveform features,the violation indicator of the specimen stress equilibrium in the split Hopkinson pressure bar test is identified for concrete-like damage softening materi-als.The concrete material behaviors for stress non-equilibrium are further analyzed,e.g.the dynamic increase factor(DIF)and damage development,etc.The conception of“damage failure volume”is introduced,and a new method of defining the development of concrete dynamic damage is given in the nimierical study.What’s more,the“compression wave”and“double peak”phenomena observed in the experiment are further interpreted based on the means of numerical simulation.Waveform features how to reflect the concrete material properties is also concluded.The results show that,the disappearance of the“double peak” phenomenon of reflection curve under high strain rate can be regarded as the indicator of the violation of stress equilibrium.After the violation of the stress equilibrium,the relevant DIFs of the concrete specimen will not change significantly.Especially,the concrete specimen will turn into structural response from material response.The conception of“damage failure volume”can well explain the generation of the“double peak”phenomenon of the reflection curve.The “compression wave” phenomenon of reflection curve under lower strain rates is derived from the unloading expansion recovery of the concrete specimen.Furthermore,under the same loading condition,the amplitude of the first peak of the reflection curve can be used as the evaluation standard of the bonding quality between mortar and aggregates.

NUMERICAL INVESTIGATION OF FLOW OVER A WEIR

Water flow over a weir was simulated numerically in this paper. The numerical model involved the Reynolds equation for mean flow filed with the k-ε turbulent model. To track free surface movements the VOF method with geometric reconstruction approach was employed. Three flow patterns including surface jet, surface wave and plunging jet were simulated in this paper. The free surface profile, velocity field and the distribution of shear stress on the bottom at downstream of the weir were obtained. The results of present numerical model, inviscid model and the rigid-lid assumption were compared with experimental data. It is shown that the present numerical model has great advantage to simulate the flow over a weir. The validities of the inviscid model and the rigid-lid assumption were also discussed.

Investigation of the Mechanism of Grout Penetration in Intersected Fractures

To study the penetration mechanism of cement-based slurry in intersected fractures during grouting and the related pressure distribution,we have used two different variants of cement,namely,basic cement slurry and fast-setting cement slurry.The influence of a retarder,time-varying viscosity,fracture width and location of injection hole is also considered.A finite element software is used to implement two and three-dimensional numerical models for grouting of intersected fractures in hydrostatic conditions.Results show that there are significant differences in the diffusion morphology and pressure distribution depending on the considered cement slurry.Retarder can effectively slow down the rising rate of injection pressure and extend the diffusion distance of grout.The influence of the branch fracture is more important when basic cement slurry is considered,indicating that the change of grout pressure is correlated with the slurry viscosity.The faster the viscosity increases,the less evident is the effect.