Experimental and Three-Dimensional Numerical Simulation of Phenomena Induced by Submerged Oblique Jet Scouring

Scouring experiments were conducted using a three-dimensional laser scanning technology for angles of the jet spanning the interval from 0°to 30°,and the characteristics of the scour hole in equilibrium conditions were investigated accordingly.The results indicate that the optimal scouring effects occur when the jet angle is in the ranges between 15°and 20°.Moreover,the dimensionless profiles of the scour hole exhibit a high degree of similarity at different jet angles.Numerical simulations conducted using the Flow-3D software to investigate the bed shear stress along the jet impingement surface have shown that this stress is influenced by both the resultant force and the jet impingement surface area.It reaches its maximum value when the jet is vertical,decreases rapidly as the jet starts to tilt,then increases slightly,and decreases again significantly when the angle exceeds 20°.

Bulking factor of the strata overlying the gob and a three-dimensional numerical simulation of the air leakage flow field

The present study examines the results of the researches related to the gob bulking factor carried out at home and abroad.A mathematical function of a three-dimensional gob bulking factor is described based on a three-dimensional gob model.The method of taking value for interstice and permeability ratios is also proposed.The law of air leakage of fully mechanized top coal is researched in this study.The results show that the speed of air flow near the upper and lower crossheadings is higher than that in the central section of the gob at the same distance from the working face.When the amount of air at the working face exceeds a critical amount,the width of the spontaneous combustion zone in the upper and lower crossheadings is also larger than that in the central section.In this situation,the key is preventing the coal left in the upper and lower crossheadings from self-igniting.Reducing the amount of air at the working face can decrease the width of the spontaneous combustion zone,especially the width near the upper and lower crossheadings.This also moves the spontaneous combustion zone in the direction of the working face.It can prevent the coal in the gob from self-igniting by making the coal left in the crossheadings to be inert and by effectively controlling the amount of air at the working face.

Three-dimensional numerical simulation of crown spike due to coupling effect between bubbles and free surface

The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.

Three-dimensional numerical simulation of M_2 internal tides in the Luzon Strait

A three-dimensional isopycnic-coordinate internal tidal model is employed to investigate the generation, propagation, vertical structure and energy conversion ofM2 internal tides in the Luzon Strait （LS） with mooring observations. Simulated results, especially the tidal current amplitudes, agree well with observations, demonstrating the reasonability and accuracy of the model. Results indicate that M2 internal tides mainly propagate into three directions horizontally, i.e., eastward towards the western Pacific Ocean, westward towards the Dongsha Island and southwestward towards the South China Sea Basin. In the horizontal direction, tidal current amplitudes decrease as distance increases away from the LS; in the vertical direction, they show an obvious decreasing tendency with depth. Between the double ridges of the LS, a clockwise gyre of M2 baroclinic energy flux appears, which is caused by reflections of M2 internal tides at supercritical topographies, and resonance of M2 internal tides happens along 19.5° and 21.5°N due to the heights and separation distance of the double ridges. The total energy conversion in the LS is about 14.20 GW.

Three-Dimensional Numerical Simulations of the Effects of a Cold Water Surface on the Evolution and Propagation of Thunderstorms

The influences of large areas of semi-unbounded cold water surface on the evolution, propagation and precipitation production of thunderstorms are simulated by using a fully elastic three-dimensional numerical hailstorm model. Real sounding profiles for temperature, humidity and wind are employed. The model has successfully simulated the significant modification of the propagation path of thunderstorms near the cold water area. The path change can be either' along-bank' or' toward-bank', depending on the position of the storm system relative to convergence zone of the water-land circulation. The simulations also show that thunderstorms developing or propagating within the convergence zone of local circulation will be intensified and produce much heavier hail, whereas those over cold water surface or the air that has been cooled by the water will be strongly inhibited.The influence of the cold water surface on thunderstorm characters is largely dependent upon the direction and intensity of the low-level wind.

Three-dimensional numerical simulation of internal tides that radiated from the Luzon Strait into the Western Pacific

Recent satellite altimeter observations have indicated that internal tides （ITs） from the Luzon Strait （LS） propagate more than 2 500 km into the Western Pacific （WP）. This study used a high-resolution three-dimensional numerical model to reproduce and examine the ITs radiation process. The propagation of diurnal and semidiurnal ITs showed different patterns and variations. Diurnal ITs with lower frequency were affected more by the earth＇s rotation and they were bent more toward the equator than semidiurnal ITs. ITs phase speeds are functions of latitude and diurnal ITs showed greater distinctions of phase speeds during propagation. For M2 ITs, the wavelength remained nearly unchanged but the beam width increased significantly during propagation away from the LS. For diurnal ITs （K1 and O0, the wavelength decreased noticeably with latitude, while the beam width varied little during propagation because of blocking by land. Baroclinic energy was also examined as a complement to satellite results reported by Zhao （2014）. The magnitude of the generated baroclinic energy flux reduced remarkably within 300 km from the generation site but it then decayed slowly when propagating into abyssal sea. Baroclinic energy of diurnal ITs was found to dissipate at a slower rate than semidiurnal ITs along the main propagation path in the WP.

Three-dimensional numerical simulation of geothermal fi eld in space-wavenumber domain

Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.

Three-dimensional Numerical Simulation of Combustion Field in the Combustion Chamber

In order to study the effect of rotation on the combustion in the underwater vehicle,a two-phase turbulent combustion process is described with Reynolds stress turbulence model,eddy-dissipation turbulent combustion model,P-1 radiation model and particle tracking model of liquid. The flow in the rotating combustion chamber is simulated at two different working speeds,0?r/min and 1?000?r/min by Fluent software. The temperature,gas velocity,static pressure of wall and fuel concentration are computed and compared. The results show that the combustion in rotating combustor is faster and more effective.

Three-Dimensional Finite Element Numerical Simulation and Physical Experiment for Magnetism-Stress Detecting in Oil Casing

The casing damage has been a big problem in oilfield production. The current detection methods mostly are used after casing damage, which is not very effective. With the rapid development of China's offshore oil industry, the number of offshore oil wells is becoming larger and larger. Because the cost of offshore oil well is very high, the casing damage will cause huge economic losses. What's more, it can also bring serious pollution to marine environment. So the effective methods of detecting casing damage are required badly. The accumulation of stress is the main reason for the casing damage. Magnetic anisotropy technique based on counter magnetostriction effect can detect the stress of casing in real time and help us to find out the hidden dangers in time. It is essential for us to prevent the casing damage from occurring. However, such technique is still in the development stage. Previous studies mostly got the relationship between stress and magnetic signals by physical experiment, and the study of physical mechanism in relative magnetic permeability connecting the stress and magnetic signals is rarely reported. The present paper uses the ANSYS to do the three-dimensional finite element numerical simulation to study how the relative magnetic permeability works for the oil casing model. We find that the quantitative relationship between the stress' s variation and magnetic induction intensity's variation is: Δδ =K* ΔB, K = 8.04×109, which is proved correct by physical experiment.

A three-dimensional Eulerian method for the numerical simulation of high-velocity impact problems

In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.

Numerical Simulation of Water Exchange Characteristics of the Jiaozhou Bay Based on A Three-Dimensional Lagrangian Model

Based on theory of three-dimensional hydrodynamics, an Euler-Lagrangian particle model is established to study the transport and water exchange capability in the Jiaozhou Bay. The three-dimensional hydrodynamic model, driven by tide and wind, is used to study the effects of wetting and drying of estuarine intertidal flats by the dry-wet grid technology based on the Estuarine, Coastal and Ocean Model （ECOM）. The particle model includes the advection and the diffusion processes, of which the advection process is simulated with a certain method, and the diffusion process is simulated with the random walk method. The effect of the intertidal zone, the turbulent diffusion and the timescales of the water exchange are also discussed. The results show that a moving boundary model can simulate the transport process of the particle in the intertidal zone, where the particles are transported for a longer distance than that of the stationary result. Simulations with and without the turbulent random walk show that the effect of turbulent diffusion is very effective at spreading particles throughout the estuary and speeding up the particle movement. The spatial distribution of residence time is given to quantify the water exchange capability that has very important ramifications to water quality. The effect of wind on the water exchange is also examined and the southeasterly wind in summer tends to block the water exchange near the northeast coast, while the northerly wind in winter speeds up the transport process. These results indicate that the Lagrangian particle model is applicable and has a large potential to help understanding the water exchange capability in estuaries, which can also be useful to simulate the transport process of contaminant.

Numerical Simulation of Nonlinear Three-Dimensional Waves in Water of Arbitrary Varying Topography

The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.

The Effect of Three-Dimensional Variational Data Assimilation of QuikSCAT Data on the Numerical Simulation of Typhoon Track and Intensity

In this paper, the three-dimensional variational data assimilation scheme (3DVAR) in the mesoscale model version 5 (MM5) of the US Pennsylvania State University/National Center for Atmospheric Research is used to study the effect of assimilating the sea-wind data from QuikSCAT on the prediction of typhoon track and intensity. The case of Typhoon Dujuan (2003) is first tested and the results show appreciable improvements. Twelve other cases in 2003 are then evaluated. The assimilation of the QuikSCAT data produces significant impacts on the structure of Dujuan in terms of the horizontal and vertical winds, sea-level pressure and temperature at the initial time. With the assimilation, the 24-h (48-h) track prediction of 11 (10) out of the 12 typhoons is improved. The 24-h (48-h) prediction of typhoon intensity is also improved in 10 (9) of the 12 cases. These experiments therefore demonstrate that assimilation of the QuikSCAT sea-wind data can increase the accuracy of typhoon track and intensity predictions through modification of the initial fields associated with the typhoon.

Numerical Simulation of the Whole Three-Dimensional Flow in a Stirred Tank with Anisotropic Algebraic Stress Model

In accordance to the anisotropic feature of turbulent flow, ananisotropic algebraic stress model is adopted to predict theturbulent flow field and turbulent characteristics generated by aRushton disc turbine with the improved inner-outer iterativeprocedure. The predicted turbulent flow is compared with experimentaldata and the simulation by the standard k-ε turbulence model. Theanisotropic algebraic stress model is found to give better predictionthan the standard k-ε turbulence model. The predicted turbulent flowfield is in accordance to experimental data and the trend of theturbulence intensity can be effectively reflected in the simulation.

Numerical Simulation of Three-dimensional Heat and Mass Transfer in Spray Cooling of Converter Gas in a Venturi Scrubber

In order to predict the pressure drop, collection efficiency, velocity, temperature and mole fraction of vapor in an industrial venturi scrubber with water spraying for converter gas cooling, a three-dimensional model of heat and mass transfer with phase change is established. The gas flow and liquid droplets are treated as a continuous phase with a Eulerian approach and as a discrete phase with a Lagrangian approach, respectively. The coupled problem of heat, force, and mass transfers between gas flow and liquid droplets is solved by a commercial computational fluid dynamics（CFD） package, FLUENT. The numerical results show that the water injections have an important influence on the distributions of pressure, velocity, temperature, and mole fraction of vapor, especially for the spraying region in the throat. In the spraying region, the pressure drop is higher and the velocity is lower than in other regions due to the gas-droplet drag, while the temperature is lower because the droplet absorbs large amounts of heat from the high temperature gas and the mole fraction of vapor is higher due to the phase change of the liquid droplet. A number of cases with different Water-to-gas volume flow ratios and baffle openings were simulated. The dependence of pressure drop, velocity, temperature, mole fraction of vapor, and collection efficiency on both the water-to-gas volume flow ratio and baffle opening are analyzed. The good agreements between simulation results and experiment data of pressure drop, temperature, and collection efficiency validate the model. The model should facilitate optimization of the venturi scrubber design in order to give better performance with lower pressure drops and higher collection efficiency.

THREE-DIMENSIONAL NUMERICAL SIMULATION OF AERATED FLOWS DOWNSTREAM SUDDEN FALL AERATOR EXPANSION-IN A TUNNEL

Air entrainment is known to be one of efficient and inexpensive methods to prevent cavitation damages in hydropower projects.The shape of sudden expansion-fall is used as a common device for mitigating cavitation erosions.The complex flow patterns with cavitation are numerically simulated by using the realizable k-εturbulence model and the air-water mixture model.The calculated results are compared well with the experimental results as well as those obtained with the k-εturbulence model with the Volume Of Fluid（VOF）Model.The calculated results agree well with the experimental data for the aeration cavity and wall pressure.Moreover,the air concentration near sidewall is simulated by a mixture model.It is found that the mixture turbulence model is superior to the VOF turbulence model.

Numerical Simulation for Evolutionary History of Three-Dimensional Basin

Numerical simulation or evolutionary history of an oil and gas-bearing basin is to repeat geological and thermodyanomic history of basin evolution on a computer and then to quantitate petroleum generation,accumulation and migration.The mathemat-ical model describing geological and thermodynamic history of the basin evolution ischaracterised by an initial-boundary value problem of a system of nonlinear partial dif-ferential equations. In the present paper, a numerical method for three-dimensionalproblem and the analysis of its stability are established and a numerical result for apractical model is given, which shows that the abnormal pressure and paleo-temperat-ure computed are reasonable and display physical characteristics clearly as well.

Numerical simulation of asphalt mixture based on three-dimensional heterogeneous specimen

In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.

Three-dimensional discrete element numerical simulation of Paleogene salt structures in the western Kuqa foreland thrust belt

Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.

Three-Dimensional Numerical Simulation of Ice Crystal Melting in Jet Engine

A method was presented for calculating the melting rate of ice crystals in different axial positions in a compressor. The volume of fluid model and enthalpy-porosity model were used to simulate the melting process of single ice crystal. The result was validated by experimental data from literature. The change of the liquid fraction of the ice crystals with the axial position was obtained when the temperature distribution inside the compressor was introduced into the melting process. The liquid fraction of the spherical ice crystals was fitted with the ice crystal diameter as a variable, and the general distribution function for computing liquid fraction in any axial position in the compressor was obtained. The function of axial position as a variable for non-spherical ice crystal was also got.