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.

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

The numerical simulation of thesource development of the Tangshan earthquake

According to the precursory data of the 1976 Tangshan earthquake( M =7.8), and combining with the research results of rock fracture, seismic source, tectonics, seismic precursors, dynamic condition and crust structure, the seismogenic process of Tangshan earthquake has been numerically simulated by using the model of water saturated porous medium with a rhombic hard inclusion. First, a system of seismogenic dynamics equations of solid liquid two phase medium has been set up under considering rock inelastic dilation, strain softening and pore water permeation. Then, the finite difference method to solve the equation system in incremental form has been presented. Finally, the seismogenic process of Tangshan earthquake and the temporal spatial characteristics of its precursors have been numerically simulated when the hard inclusion includes a soft fault belt. The numerically simulating result shows that the development of Tangshan earthquake mainly undergoes following processes: elastic accumulation, early inelastic cubic dilation, accelerative softening (failure and creep) of fault first time, inelastic dilation twice in Tianjin—Ninghe area, accelerative softening of fault twice and instability. Corresponding to the fault accelerative softening, elastic recovery and inelastic dilation twice in some parts of seismogenic medium. By comparing with the data of precursory observation, ones can find the preceding simulating result explains the varied characters of practical data better. Meanwhile, the numerical simulation further displays the temporal spatial change law and complex patterns of the field of Tangshan earthquake development as well. It makes a base study for the physical mechanism of change of the seismic precursory field.

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.

Numerical Modeling of a Spar Platform Tethered by a Mooring Cable

Virtual simulation is an economical and efficient method in mechanical system design. Numerical modeling of a spar platform, tethered by a mooring cable with a spherical joint is developed for the dynamic simulation of the floating structure in ocean. The geometry modeling of the spar is created using finite element methods. The submerged part of the spar bears the buoyancy, hydrodynamic drag force, and effect of the added mass and Froude-Krylov force. Strip theory is used to sum up the forces acting on the elements. The geometry modeling of the cable is established based on the lumped-mass-and-spring modeling through which the cable is divided into 10 elements. A new element-fixed local frame is used, which is created by the element orientation vector and relative velocity of the fluid, to express the loads acting on the cable. The bottom of the cable is fixed on the seabed by spring forces, while the top of the cable is connected to the bottom of the spar platform by a modified spherical joint. This system suffers the propagating wave and current in the X-direction and the linear wave theory is applied for setting of the propagating wave. Based on the numerical modeling, the displacement-load relationships are analyzed, and the simulation results of the numerical modeling are compared with those by the commercial simulation code, Proteus DS. The comparison indicates that the numerical modeling of the spar platform tethered by a mooring cable is well developed, which provides an instruction for the optimization of a floating structure tethered by a mooring cable system.

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.

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.

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.

气动优化设计和计算流体力学在无人机设计中的应用

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.

A new traffic model with a lane-changing viscosity term

In this paper, a new continuum traffic flow model is proposed, with a lane-changing source term in the continuity equation and a lane-changing viscosity term in the acceleration equation. Based on previous literature, the source term addresses the impact of speed difference and density difference between adjacent lanes, which provides better precision for free lane-changing simulation; the viscosity term turns lane-changing behavior to a ＂force＂ that may influence speed distribution. Using a flux-splitting scheme for the model discretization, two cases are investigated numerically. The case under a homogeneous initial condition shows that the numerical results by our model agree well with the analytical ones; the case with a small initial disturbance shows that our model can simulate the evolution of perturbation, including propagation,dissipation, cluster effect and stop-and-go phenomenon.

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.

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.

Simulations of vertical jet penetration using a filtered two-fluid model in a gas-solid fluidized bed

The influence of a vertical jet located at the distributor in a cylindrical fluidized bed on the flow behavior of gas and particles was predicted using a filtered two-fluid model proposed by Sundaresan and coworkers. The distributions of volume fraction and the velocity of particles along the lateral direction were investigated for different jet velocities by analyzing the simulated results. The vertical jet penetration lengths at the different gas jet velocities have been obtained and compared with predictions derived from empirical correlations; the predicted air jet penetration length is discussed. Agreement between the numerical simulations and experimental results has been achieved.

Transient simulation of regression rate on thrust regulation process in hybrid rocket motor

The main goal of this paper is to study the characteristics of regression rate of solid grain during thrust regulation process. For this purpose, an unsteady numerical model of regression rate is established. Gas–solid coupling is considered between the solid grain surface and combustion gas.Dynamic mesh is used to simulate the regression process of the solid fuel surface. Based on this model, numerical simulations on a H2O2/HTPB（hydroxyl-terminated polybutadiene） hybrid motor have been performed in the flow control process. The simulation results show that under the step change of the oxidizer mass flow rate condition, the regression rate cannot reach a stable value instantly because the flow field requires a short time period to adjust. The regression rate increases with the linear gain of oxidizer mass flow rate, and has a higher slope than the relative inlet function of oxidizer flow rate. A shorter regulation time can cause a higher regression rate during regulation process. The results also show that transient calculation can better simulate the instantaneous regression rate in the operation process.