Numerical Simulation and Dynamical Analysis for Low Salinity Water Lens in the Expansion Area of the Changjiang Diluted Water

The low salinity water lenses（LSWLes） in the expansion area of the Changjiang diluted water（CDW） exist in a certain period of time in some years. The impact of realistic river runoff, ocean currents and weather conditions need to be taken into account in the dynamical analysis of LSWL, which is in need of research. In this paper, the POM-σ-z model is used to set up the numerical model for the expansion of the CDW. Then LSWL in summer 1977 is simulated, and its dynamic mechanism driven by wind, tide, river runoff and the Taiwan Warm Current is also analyzed. The simulated results indicate that the isolated LSWL detaches itself from the CDW near the river mouth, and then moves towards the northeast region outside the Changjiang Estuary. Its maintaining period is from July 26 to August 11. Its formation and development is mainly driven by two factors. One is the strong southeasterly wind lasting for ten days. The other is the vertical tidal mixing during the transition from neap tide to spring tide.

Dynamic Interference Photoelectron Spectra in Double Ionization:Numerical Simulation of 1D Helium

We report our numerical simulation on the dynamic interference photoelectron spectra for a one-dimensional （1D） He model exposed to intense ultrashort extreme ultraviolet （XUV） laser pulses. The results demonstrate an unambiguous interference feature in the photoelectron spectra, and the interference is unveiled to originate from the dynamic Stark effect. The interference photoelectron spectra are prompted for intense sub-femtosecond XUV laser pulses in double ionization. The stationary phase picture is corroborated qualitatively in the two-electron system. The ability of probing the dynamic Stark effect by the photoelectron spectra in a pragmatic experiment of single-photon double ionization of He may shed light on further investigation on multi-electron atoms and molecules.

Study on dynamic numerical simulation of overburden movement mining under railways

Focused upon the dynamic changing regularity and tendency of overburden movement and bed separation developing during the strip mining and stoping under rail- ways were simulated by ANSYS model.The results show that movement and rheology developing of key strata is one of the main reasons to cause surface subsidence and to further induce deformation of surface railways.It provids a basis for taking efficient meas- ures to control the deformation of surface railways.

Application of Aerodynamic Optimization Design and Dynamic Numerical Simulation in UAV Design

In the past two decades,the world’s unmanned aerial vehicle(UAV)industry has developed rapidly.Various kinds of UAVs have been used in military and civilian fields.Based on the characteristics of UAVs and the development of aerodynamics,this article analyzes the development of aerodynamic optimization design and dynamic numerical simulation technology,then lists engineering applications.Both aerodynamic optimization design and dynamic numerical simulation have greatly shortened the UAV design period and reduced the research and design cost.These two methods gradually replace traditional methods such as wind tunnel test.

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 study on static and dynamic fracture evolution around rock cavities

In this paper,a numerical code,RFPA2D（rock failure process analysis）,was used to simulate the initiation and propagation of fractures around a pre-existing single cavity and multiple cavities in brittle rocks.Both static and dynamic loads were applied to the rock specimens to investigate the mechanism of fracture evolution around the cavities for different lateral pressure coefficients.In addition,characteristics of acoustic emission（AE） associated with fracture evolution were simulated.Finally,the evolution and interaction of fractures between multiple cavities were investigated with consideration of stress redistribution and transference in compressive and tensile stress fields.The numerically simulated results reproduced primary tensile,remote,and shear crack fractures,which are in agreement with the experimental results.Moreover,numerical results suggested that both compressive and tensile waves could influence the propagation of tensile cracks;in particular,the reflected tensile wave accelerated the propagation of tensile cracks.

Numerical simulation and experimental verification of bubble size distribution in an air dense medium fluidized bed

Bubble size distribution is the basic apparent performance and obvious characteristics in the air dense medium fluidized bed (ADMFB). The approaches of numerical simulation and experimental verification were combined to conduct the further research on the bubble generation and movement behavior. The results show that ADMFB could display favorable expanded characteristics after steady fluidization. With different particle size distributions of magnetite powder as medium solids, we selected an appropriate prediction model for the mean bubble diameter in ADMFB. The comparison results indicate that the mean bubble diameters along the bed heights are 35 mm < D b < 66 mm and 40 mm < D b < 69 mm with the magnetite powder of 0.3 mm+0.15mm and 0.15mm+0.074mm, respectively. The prediction model provides good agreements with the experimental and simulation data. Based on the optimal operating gas velocity distribution, the mixture of magnetite powder and <1mm fine coal as medium solids were utilized to carry out the separation experiment on 6-50mm raw coal. The results show that an optimal separation density d P of 1.73g/cm 3 with a probable error E of 0.07g/cm 3 and a recovery efficiency of 99.97% is achieved, which indicates good separation performance by applying ADMFB.

NUMERICAL SIMULATION STUDY OF THE INNER-CORE STRUCTURES AND THE MECHANISM FOR INSHORE STRENGTHENING OF SOUTH CHINA SEA TYPHOON VONGFONG (0214) DURING LANDFALL

An explicit simulation with a fine mesh at intervals of 6 km is used to explore the inner-core structures of Vongfong (0214). The dynamic mechanism for the inshore strengthening of Vongfong is examined. It is found as follows. (1) The radius of maximum wind of the axisymmetric structures of the typhoon decreased with height during its mature stage. When Vongfong was inshore, the strongest low-layer inflow located in front of it and the outflow was to the rear of it, which was just reversed from the Atlantic hurricanes and other Pacific typhoons. (2) The dynamic and thermodynamic fields were highly asymmetric in structure. Convection was stronger in the northwest quadrant of the typhoon than in the southeast; the strongest convective cloud bands were consistent with the maximum wind region. During its strengthening stage, it was cold west of and warm east of the eye in the lower layer but warm in the west and cold in the east of the mid-upper layer. During its mature stage, a warm-core structure was evident in the lower and middle-upper layer. (3) The interactions between a mid-latitude cold low in the middle-upper troposphere and the typhoon were responsible for the latter to strengthen inshore. Firstly, the outer circulation of the cold low entered the typhoon from the middle troposphere when an outer cold airflow from the cold low flowed into the northwest quadrant of the typhoon so that geopotentially instable energy increased and convection developed. Secondly, the downdraft in the cold low was just the corresponding branch of the secondary circulation of the typhoon system; when the cold low weakened while moving south, the typhoon strengthened inshore. Due to the CISK mechanism, these two phenomena might be realized.

Numerical Modeling of a Spherical Buoy Moored by a Cable in Three Dimensions

Floating facilities have been studied based on the static analysis of mooring cables over the past decades. To analyze the floating system of a spherical buoy moored by a cable with a higher accuracy than before, the dynamics of the cables are considered in the construction of the numerical modeling. The cable modeling is established based on a new element frame through which the hydrodynamic loads are expressed efficiently. The accuracy of the cable modeling is verified with an experiment that is conducted by a catenary chain moving in a water tank. In addition, the modeling of a spherical buoy is established with respect to a spherical coordinate in three dimensions, which can suffers the gravity, the variable buoyancy and Froude-Krylov loads. Finally, the numerical modeling for the system of a spherical buoy moored by a cable is established, and a virtual simulation is proceeded with the X- and Y-directional linear waves and the X-directional current. The comparison with the commercial simulation code Proteus DS indicates that the system is accurately analyzed by the numerical modeling. The tensions within the cable, the motions of the system, and the relationship between the motions and waves are illustrated according to the defined sea state. The dynamics of the cables should be considered in analyzing the floating system of a spherical buoy moored by a cable.

Towards full predictions of temperature dynamics in McNary Dam forebay using OpenFOAM

Hydroelectric facilities impact water temperature; low velocities in a reservoir increase residence time and enhance heat exchange in surface layers. In this study, an unsteady three-dimensional model was developed to predict the temperatm＇e dynamics in the McNary Dam forebay. The model is based on the open-source code OpenFOAM. RANS equations with the Boussinesq approximation were used to solve the flow field. A： realizable k-ε model that accounts for the production of wind turbulence was developed. Solar radiation and convective heat transfer at the free surface were included. The result of the model was compared with the field data collected on August 18, 2004. Changes in diurnal stratification were adequately predicted by the model. Observed vertical and lateral temperature distributions were accurately captured. Results indicate that the model can be used as a numerical tool to assess structural and operational alternatives to reduce the forebay temperature.

Mechanism on simulation and experiment of pre-crack seam formation in stope roof

The pre-crack blast technology has been used to control the induction caving area in the roof. The key is to form the pre-crack seam and predict the effect of the seam. The H-J-C blast model was built in the roof. Based on the theories of dynamic strength and failure criterion of dynamic rock, the rock dynamic damage and the evolution of pre-crack seam were simulated by the tensile damage and shear failure of the model. According to the actual situation of No. 92 ore body test stope at Tongkeng Mine, the formation process of the pre-crack blast seam was simulated by Ansys/Ls-dyna software, the pre-crack seam was inspected by a system of digital panoramic borehole camera. The pre-crack seam was inspected by the system of digital panoramic borehole in the roof. The results of the numerical simulation and inspection show that in the line of centers of pre-hole, the minimum of the tensile stress reaches 20 MPa, which is much larger than 13.7 MPa of the dynamic tensile strength of rock. The minimum particle vibration velocity reaches 50 cm/s, which is greater than 30-40 cm/s of the allowable vibration velocity. It is demonstrated that the rock is destroyed near the center line and the pre-crack is successfully formed by the large diameters and large distances pre-crack holes in the roof.

Mixing processes and patterns of fluids in alkane-CO_(2)-water systems under high temperature and high pressure——Microscopic visual physical thermal simulations and molecular dynamics simulations

Various types of geofluids exist in deep and ultra-deep layers in petroliferous basins.The geofluids are much more active under high-temperature and high-pressure(HTHP)conditions,but their properties are unclear.We simulated the mixing of different fluids in CH_(4)/C_(3)H_(8)/C_(6)H_(14)/C_(8)H_(18)-water systems and C_(6)H_(14)/C_(8)H_(18)-CO_(2)-H_(2)O systems at temperatures of 25℃ to 425℃ and pressures of 5 MPa to 105 MPa,using an in-situ micron quartz capillary tube thermal simulation system and molecular dynamics numerical simulation software.The mixing processes,patterns,and mechanisms of various fluids were analyzed at microscale under increasing temperature and pressure conditions.The results show that the miscibility of fluids in the different alkane-H_(2)O and alkane-CO_(2)-H_(2)O systems is not instantaneous,but the miscibility degree between different fluid phases increases as the temperature and pressure rise during the experiments.The physical thermal experiments(PTEs)show that the mixing process can be divided into three stages:initial miscibility,segmented dynamic miscibility,and complete miscibility.The molecular dynamics numerical simulations(MDNSs)indicate that the mixing process of fluids in the alkane-H_(2)O and alkane CO_(2)-H_(2)O systems can be divided into seven and eight stages,respectively.The carbon number affects the miscibility of alkanes and water,and the temperature and pressure required to reach the same miscibility stage with water increase with the carbon number(C_(3)H_(8),C_(6)H_(14),CH_(4),C_(8)H_(18)).CO_(2) has a critical bridge role in the miscibility of alkanes and water,and its presence significantly reduces the temperatures required to reach the initial,dynamic,and complete miscibility of alkanes and water.The results are of great significance for analyzing and understanding the miscibility of geofluids in deep and ultra-deep HTHP systems.

Numerical simulation of heat transfer process in cement grate cooler based on dynamic mesh technique

A grate cooler is key equipment in quenching clinker and recovering heat in cement production. A two-dimensional numerical model based on a 5000 t/d cement plant is proposed to for a study on the gas-solid coupled heat transfer process between the cooling air and clinker in a grate cooler. In this study, we use the Fluent dynamic mesh technique and porous media model through which the transient temperature distribution with the clinker motion process and steady temperature and pressure distribution are obtained. We validate the numerical model with the operating data of the cooling air outlet temperature. Then, we discuss the amount of mid-temperature air outlet and average diameter of clinker particles, which affect the heat effective utilization and cooling air pressure drop in clinker layer. We found that after adding one more mid-temperature air outlet, the average temperature of the air flowing into the heat recovery boiler increases by 29.04℃ and the ratio of heat effective utilization increases by 5.3%. This means heat recovery is more effective on adding one more mid-temperature air outlet. Further, the smaller the clinker particles, the more is the pressure drop in clinker layer; thus more power consumption is needed by the cooling fan.

THE NUMERICAL SIMULATION OF TWO-DIMENSIONAL DYNAMICAL CLIMATE MODEL WITH MOUNTAIN FORCING

In order to investigate the effects of trace gases on climate variation in the atmosphere, we have devel- oped a primitive equation two-dimensional dynamical climate model with five levels. A series of simula- tion results and discussions are shown in this paper, indicating that the model is useful and can correctly reproduced the main feature of the general atmospheric circulation and its seasonal changes. In addition, we have discussed the role of the Qinghai-Xizang Plateau on the formation process of summer monsoon in South Asia and found that the thermal effect of the Qjnghai-Xizang Plateau may not be the main factor controlling the onset and the variation of the summer monsoon in South Asia.

新一代运载火箭发射燃气动力学数值模拟

By using the mesh resolution control method based on the nozzle scale,a paralleled super numerical simulation and high-quality mesh model of the launch jet dynamics for new-generation launch vehicles was developed.Based upon this,a transient numerical simulation method,combining the pressure and velocity,tightly coupled algorithm and SST turbulence model,was used to complete the unsteady numerical simulation of the launch jet dynamics of the new-generation launch vehicles.The numerical simulation results of the launch jet dynamics,for the new-generation launch vehicles,demonstrated that despite the complex structure of the launch platform,the jet flows of the core stage and booster engines were generally smoothly channeled into the double deflecting trench through the launch platform’s diversion hole at the initial stage of ignition.After the lift off,the jet flows of the core stage and the booster engines began to affect and ablate the grillage-shaped beam and the adjoined surface of the launch platform adjacent to the booster engines.At a higher altitude after lift off,it could be seen for the new-generation launch vehicles the ablation range of high temperature and high-speed jet flows on the launch platform further expanded,which would have a severe ablation effect on the fuel filling tower near the booster engines and even all the support arms.The numerical simulation of launch jet dynamics also established that the jet flows embers at the bottom of the core stage rocket body continued to be affected for an extended period of time due to the large number of nozzles in the new-generation launch vehicles engine and the weak suction effect of the jet flows in the core-stage engines.

A Crashworthy Device Against Ship-OWT Collision and Its Protection Effects on the Tower of Offshore Wind Farms

At present, more and more offshore wind farms have been built anti ntnnerous projects are on the drawing tables. Therefore, the study on the safety of collision between ships and offshore wind turbines （OWT） is of great practical signifieance. The present study takes the advantage of the famous LS-DYNA explicit code to simulate the dynamic proeess of the collision between a typical 3MW offshore wind turbine model with monopile fi）undation and a simplified 2000t-class ship model. In the simulation, the added mass effect of the ship, contact nonlinearity of collision, material nonlinearity of steel and aluminum foam and adaptive mesh tectmique for large structure deformation have been taken into considera- tion. Proposed is a crashworthy device for OWF of new conceptual steel sphere shell-cireular ring aluminum foam pad, and the good pe.rfurmanee of the device under the conditions of ship-OWT front impact and side impact has been verified from the views of theoretical analysis and numerical results. The new crashworthy device can effectively smooth the contact force and reduce the top structure dynamic response, using its own structure plastic deformation to absorb most of the ship collision enerty. As a result, the main structure of the OWF and the inside key electric control equipments can be saved by scarifying the structural plastic deformation of new sphere crashworthy device. What is more, the sphere configuratiun design of the crashworthy device can effectively guide the ship to run away from the main OWT structure and reduce the damage of the ship and OWT to some degree during side impact.

Process Modeling of Ferrofluids Flow for Magnetic Targeting Drug Delivery

Among the proposed techniques for delivering drugs to specific sites within the human body, magnetic targeting drug delivery surpasses due to its non-invasive character and its high targeting efficiency. Although there have been some analyses theoretically for magnetic drug targeting, very few researchers have addressed the hydrodynamic models of magnetic fluids in the blood vessel of human body. This paper presents a mathematical model to describe the hydrodynamics of ferrofluids as drug carriers flowing in a blood vessel under the applied magnetic field. A 3D flow field of magnetic particles in a blood vessel model is numerically simulated in order to further understand clinical application of magnetic targeting drug delivery. Simulation results show that magnetic nanoparticles can be enriched in a target region depending on the applied magnetic field intensity. Magnetic resonance imaging confirms the enrichment of ferrofluids in a desired body tissue of Sprague-Dawley rats. The simulation results coincide with those animal experiments. Results of the analysis provide the important information and can suggest strategies for improving delivery in favor of the clinical application.

In-situ remediation of deep petroleum-contaminated soil injection

A computational fluid dynamics(CFD)numerical simulation and field experiment were used to investigate optimal operating parameters of high-pressure jet grouting equipment and clarify the boundary law of the injection area in the remediation process.The response surface optimization design results show that the optimal injection pressure is 30 MPa,rotation speed is 23 r/min,commission speed is 30 cm/min,and the optimal injection diameter is 147.3 cm.Based on the CFD numerical simulation,the ratio of the injection core,turbulent zone,and seepage zone is approximately 1∶4∶2.The distribution law of jet core,turbulence zone and seepage zone at different cross-sections under 30 MPa operating conditions is as follows:The jet core radius is approximately 100 mm,the turbulence zone is mainly distributed at 100 to 500 mm,the seepage zone is mainly distributed at 500 to 700 mm,the seepage zone could be completed within 2 h,and the proportion of the three boundary zones in the injection zone is similar to that of the numerical simulation.This study provides theoretical parameters and practical reference for the remediation of deep pollution via in-situ chemical oxidation in the Loess Plateau soil environment.

The effects of caudal fin deformation on the hydrodynamics of thunniform swimming under self-propulsion

To investigate the effects of the caudal fin deformation on the hydrodynamic performance of the self-propelled thunniform swimming,we perform fluid-body interaction simulations for a tuna-like swimmer with thunniform kinematics.The 3-D vortices are visualized to reveal the role of the leading-edge vortex(LEV)in the thrust generation.By comparing the swimming velocity of the swimmer with different caudal fin flexure amplitudes fa,it is shown that the acceleration in the starting stage of the swimmer increases with the increase of fa,but its cruising velocity decreases.The results indicate that the caudal fin deformation is beneficial to the fast start but not to the fast cruising of the swimmer.During the entire swimming process,the undulation amplitudes of the lateral velocity and the yawing angular velocity decrease as fa increases.It is found that the formation of an attached LEV on the caudal fin is responsible for generating the low-pressure region on the surface of the caudal fin,which contributes to the thrust.Furthermore,the caudal fin deformation can delay the LEV shedding from the caudal fin,extending the duration of the low pressure on the caudal fin,which will cause the caudal fin to generate a drag-type force over a time period in one swimming cycle and reduce the cruising speed of the swimmer.

Numerical study on morphological characteristics of rotational natural supercavitation by rotational supercavitating evaporator with optimized blade shape

In view of the supercavitation effect, a novel device named the rotational supercavitating evaporator (RSCE) has been designed for the desalination. In order to improve the blade shape of the rotational cavitator in the RSCE for the performance optimization, the blade shapes of different sizes are designed by utilizing the improved calculation method for the blade shape and the validated empirical formulae based on previous two-dimensional numerical simulations, from which the optimized blade shape with the wedge angle of 45° and the design speed of 5 000 r/min is selected. The estimation method for the desalination performance parameters is developed to validate the feasibility of the utilization of the results obtained by the two-dimensional numerical simulations in the design of the three-dimensional blade shape. Three-dimensional numerical simulations are then conducted for the supercavitating flows around the rotational cavitator with the optimized blade shape at different rotational speeds to obtain the morphological characteristics of the rotational natural supercavitation. The results show that the profile of the supercavity tail is concaved toward the inside of the supercavity due to the re-entrant jet. The empirical formulae for estimating the supercavity size with consideration of the rotation are obtained by fitting the data, with the exponents different from those obtained by the previous two-dimensional numerical simulations. The influences of the rotation on the morphological characteristics are analyzed from the perspectives of the tip and hub vortices and the interaction between the supercavity tail and the blade. Further numerical simulation of the supercavitating flow around the rotational cavitator made up by the blades with exit edge of uniform thickness illustrate that the morphological characteristics are also affected by the blade shape.