Numerical Study of Fluid Dynamics and Heat Transfer Induced by Plasma Discharges

A numerical investigation is conducted to explore the evolution of a plasma discharge and its interaction with the fluid flow based on a self-consistent fluid model which couples the discharge dynamics with the fluid dynamics.The effects of the applied voltage on the distribution of velocity and temperature in initially static air are parainetrically studied.Furthermore,the spatial structure of plasma discharge and the resulting force contours in streamwise and normal directions are discussed in detail.The result shows that the plasma actuator produces a net force that should always be directed away from the exposed electrode,which results in an ionic wind pushing particles into a jet downstream of the actuator.When the energy added by the plasma is taken into account,the ambient air temperature is increased slightly around the electrode,but the velocity is almost not affected.Therefore it is unlikely that the induced flow is buoyancy driven.For the operating voltages considered in this paper,the maximum induced velocity is found to follow a power law,i.e.,it is proportional to the applied voltage to the 3.5 power.This promises an efficient application in the flow control with plasma actuators.

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.

Some numerical experiments of the dynamics of the wintertime circulation in the Northern Huanghai Sea and the Bohai Sea

-Some numerical experiments have been developed to explain dynamically the mechanism of the formation of the wintertime circulation in the Northern Huanghai Sea and the Bohai Sea. From the results of the numerical experiments with Coriolis parameter having positive and negative values respectively, it is found that the current pattern characterized by flowing into the Bohai Sea in the northern part of the Bohai Straits and flowing out of the Bohai Sea in the southern part of the Bohai Straits is controlled basically by the bottom topography with 'the deepness on the northern side and shallowness on the southern side'in the Bohai Straits, and is not related to Coriolis parameter having positive and negative values. It also shows that the Huanghai Warm Current regarded as the inllow condition passing through the in-ports at the water boundary of the numerical experiments only influences the region near the water boundary and the eastern region of the Northern Huanghai Sea. In addition, the mechanism of the formation of several mesoscale eddies and the influence of the bottom topography on the wintertime circulation in the basin have been discussed.

Modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies

We present an efficient, robust computational method for modeling the Newtonian dynamics for rotation curve analysis of thin-disk galaxies. With appropriate mathematical treatments, the apparent numerical difficulties associated with singularities in computing elliptic integrals are completely removed. Using a boundary element discretization procedure, the governing equations are transformed into a linear algebra matrix equation that can be solved by straightforward Gauss elimination in one step without further iterations. The numerical code implemented according to our algorithm can accurately determine the surface mass density distribution in a disk galaxy from a measured rotation curve （or vice versa）. For a disk galaxy with a typical flat rotation curve, our modeling results show that the surface mass density monotonically decreases from the galactic center toward the periphery, according to Newtonian dynamics. In a large portion of the galaxy, the surface mass density follows an approximately exponential law of decay with respect to the galactic radial coordinate. Yet the radial scale length for the surface mass density seems to be generally larger than that of the measured brightness distribution, suggesting an increasing mass-tolight ratio with the radial distance in a disk galaxy. In a nondimensionalized form, our mathematical system contains a dimensionless parameter which we call the ＂galactic rotation number＂ that represents the gross ratio of centrifugal force and gravitational force. The value of this galactic rotation number is determined as part of the numerial solution. Through a systematic computational analysis, we have illustrated that the galactic rotation number remains within 4-10% of 1.70 for a wide variety of rotation curves. This implies that the total mass in a disk galaxy is proportional to V02 Rg, with V0 denoting the characteristic rotation velocity （such as the ＂flat＂ value in a typical ro- tation curve） and Rg the radius of the galactic disk. The predicted total galactic mass of the Milky Way is in good agreement with the star-count data.

The State Key Laboratory of Numerical Modelling for Atmospheric Science and Geophysical Fluid Dynamics

The State Key Laboratory of Numerical Modelling for Atmospheric Science and Geophysical Fluid Dynamics (LASG) was established in 1985. It was considered an Excellent Opening Laboratory in the national assessment organized by the former State Science and Technology Committee (now the Ministry of Science and Technology) in 1988. Because of this, LASG was upgraded to a State Key Laboratory in 1989. It won the honor of "Excellent

SOME DEVELOPMENTS IN VISCOUS FLUID DYNAMICS AND ITS NUMERICAL METHODS

Viscous fluid flows contain abundant "physical phenomena and the viscous fluid dynamics is of wide applications in the fields of natural and engineering sciences. After the basic equations of viscousfluiddynamics (i.e., the Navier-Stokes equations) came out, one of the most important contributions to the discipline was the boundary layer (BL) theory and the BL equations presented by Prandtl

State Key Lab of Numerical Modeling for Atmospheric Science and Geophysical Fluid Dynamics

The State Key Laboratory of Numerical Modeling for Atmospheric Science and Geophysical Fluid Dynamics (LASG) was set up as early as in 1985. It has been appraised as an Excellent State Key Lab in each of the three national assessments organized by State Economic Planning Commission, National Natural Science Foundation (NSFC) and the Ministry of Science and Technology (MOST) in 1991, 1996 and 2000, respectively.

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

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.

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 modeling and dynamic analysis of the 2017 Xinmo landslide in Maoxian County, China

A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eighty-three people were killed or went missing and more than 103 houses were destroyed. In this paper, the geological conditions of the landslide are analyzed via field investigation and high-resolution imagery. The dynamic process and runout characteristics of the landslide are numerically analyzed using a depth-integrated continuum method and Mac Cormack-TVD finite difference algorithm.Computational results show that the evaluated area of the danger zone matchs well with the results of field investigation. It is worth noting that soil sprayed by the high-speed blast needs to be taken into account for such kind of large high-locality landslide. The maximum velocity is about 55 m/s, which is consistent with most cases. In addition, the potential danger zone of an unstable block is evaluated. The potential risk area evaluated by the efficient depthintegrated continuum method could play a significant role in disaster prevention and secondary hazard avoidance during rescue operations.

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 of Dynamic Response of an Existing Subway Station Subjected to Internal Blast Loading

In order to design and retrofit a subway station to resist an internal blast, the distribution of blast loading and its effects on structures should be investigated firstly. In this paper, the behavior of a typical subway station subjected to different internal blast Ioadings was analyzed. It briefly introduced the geometric characteristics and material constitutive model of an existing two-layer and three-span frame subway station. Then three cases of different explosive charges were consid- ered to analyze the dynamic responses of the structure. Finally, the maximum principal stress, dis- placement and velocity of the columns in the three cases were obtained and discussed. It con- cluded that the responses of the columns are sensitive to the charge of explosive and the distance from the detonation. It＇s also found that the stairs between the two layers have significant effects on the distribution of the maximum principal stress of the columns in the upper layer. The explicit dynamic nonlinear finite element software ANSYS/LS-DYNA was used in this study.

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

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.

Dynamic responses of deep underground explosions based on improved Grigorian model

It is important to investigate the dynamic behaviors of deep rocks near explosion cavity to reveal the mechanisms of deformations and fractures. Some improvements are carried out for Grigorian model with focuses on the dilation effects and the relaxation effects of deep rocks, and the high pressure equations of states with Mie-Grüneisen form are also established. Numerical calculations of free field parameters for deep underground explosions are carried out based on the user subroutines which are compiled by means of the secondary development functions of LS-DYNA9703 D software. The histories of radial stress, radial velocity and radial displacement of rock particles are obtained, and the calculation results are compared with those of U.S. Hardhat nuclear test. It is indicated that the dynamic responses of free field for deep underground explosions are well simulated based on improved Grigorian model, and the calculation results are in good agreement with the data of U.S. Hardhat nuclear test. The peak values of particle velocities are consistent with those of test, but the waveform widths and the rising times are obviously greater than those without dilation effects. The attenuation rates of particle velocities are greater than the calculation results with classic plastic model, and they are consistent with the results of Hardhat nuclear test. The attenuation behaviors and the rising times of stress waves are well shown by introducing dilation effects and relaxation effects into the calculation model. Therefore, the defects of Grigorian model are avoided. It is also indicated that the initial stress has obvious influences on the waveforms of radial stress and the radial displacements of rock particles.

On numerical earthquake prediction

Can earthquakes be predicted？ How should people overcome the difficulties encountered in the study of earthquake prediction？ This issue can take inspiration from the experiences of weather forecast. Although weather forecasting took a period of about half a century to advance from empirical to numerical forecast, it has achieved significant success. A consensus has been reached among the Chinese seismological community that earth- quake prediction must also develop from empirical fore- casting to physical prediction. However, it is seldom mentioned that physical prediction is characterized by quantitatively numerical predictions based on physical laws. This article discusses five key components for numerical earthquake prediction and their current status. We conclude that numerical earthquake prediction should now be put on the planning agenda and its roadmap designed, seismic stations should be deployed and observations made according to the needs of numerical prediction, and theoretical research should be carried out.

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.