The mesoscale numerical simulation of the flow field of the Hainan Island and the Leizhou Peninsula

The flow field over Hainan Island and the Leizhou Peninsula in summer and winter is discussed with three-dimensional mesoscale model developed in the University of Virginia and using the representative meteorological data of January and July.Simulation results indicate that the local weather characteristics over the Hainan Island are distinctly influenced by theWuzhi Mountain terrain. The cloudy or rainfall weather over the northeast of the Wuzhi Mountain occurs easily, under proper large-scale conditions of flow, temperature and humidity. while west wind prevails. The overcast or rainfall weather is often induced by strong convection in the afternoon over west of the Hainan Island under easterly prevailing wind.

Numerical simulation of blood flow and interstitial fluid pressure in solid tumor microcirculation based on tumor-induced angiogenesis

A coupled intravascular-transvascular-interstitial fluid flow model is developed to study the distributions of blood flow and interstitial fluid pressure in solid tumor microcirculation based on a tumor-induced microvascular network. This is generated from a 2D nine-point discrete mathematical model of tumor angiogenesis and contains two parent vessels. Blood flow through the microvascular network and interstitial fluid flow in tumor tissues are performed by the extended Poiseuille＇s law and Darcy＇s law, respectively, transvascular flow is described by Starling＇s law; effects of the vascular permeability and the interstitial hydraulic conductivity are also considered. The simulation results predict the heterogeneous blood supply, interstitial hypertension and low convection on the inside of the tumor, which are consistent with physiological observed facts. These results may provide beneficial information for anti-angiogenesis treatment of tumor and further clinical research.

NUMERICAL SIMULATION ANALYSIS OF EXTERNAL FLOW FIELD OF WAGON-SHAPED CAR AT THE MOMENT OF PASSING

In the course of studying on aerodynamic change and its effect on steering stability and controllability of an automobile in passing, because of multi interaction streams, it is difficult to use traditional methods, such as wind tunnel test and road test. If the passing process of an automobile is divided into many time segments, so as to avoid the use of moving mesh which takes large calculation resource and CPU processing time in calculating, the segments are simulated with computational fluid dynamics （CFD） method, then the approximate computational results about external flow field will be obtained. On the basis of the idea, the change of external flow field of wagon-shaped car at the moment of passing is simulated through solving three-dimensional, steady and uncompressible N-S equations with finite volume method. Numerical simulation analysis of side force coefficient, stream lines, body surface pressure distribution of wagon-shaped car are presented and a preliminary discussion of aerodynamic characteristics of correlative situations is obtained. Finally, the C3 -x/l curve of side force coefficient（C3） of car following relative distance （x/l） between cars is obtained. By comparison, the curve is coincident well with the experimental data, which shows creditability of numerical simulation methods presented.

Wind and Snow Impacts on Urban Space Planning:A Multiphase Numerical Simulation

The wind and snow environment outside the planned space plays a key role in the comfort and safety of the human habitat in severe cold regions. Traditional studies of the external environment of human settlements, however, frequently overlook the combined impacts of wind and snow environments. Furthermore, in urban meteorological studies, it is impossible to accurately assess the wind and snow environment in specific areas or locations. In this study, a refined Computational Fluid Dynamics(CFD) multiphase flow numerical method was used to simulate a planning space's wind and snow environment. The study classified the Snowstorm Weather Grade(SWG) by incorporating the Snowstorm Weather Index(SWI) to generate calculation results of the wind environment and snow environment. In particular, 150 measurement points in the planning space were chosen for analysis and evaluation of their wind and snow environments. The results demonstrated that the SWI index can effectively correlate to the wind and snow environment calculation results. In addition, the graph of SWI showed that 55% of the measurement points had a moderate wind and snow grade SWI, which exceeds the average grade for the entire region. The practical application shows that the wind and snow environment assessment indexes and technical methods developed in this paper can be successfully applied to wind and snow environment studies in other cold cities.

Uncertainty Quantification of Numerical Simulation of Flows around a Cylinder Using Non-intrusive Polynomial Chaos

The uncertainty quantification of flows around a cylinder is studied by the non-intrusive polynomial chaos method. Based on the validation with benchmark results, discussions are mainly focused on the statistic properties of the peak lift and drag coefficients and base pressure drop over the cylinder with the uncertainties of viscosity coefficient and inflow boundary velocity. As for the numerical results of flows around a cylinder, influence of the inflow boundary velocity uncertainty is larger than that of viscosity. The results indeed demonstrate that a five-order degree of polynomial chaos expansion is enough to represent the solution of flow in this study.

Numerical simulations of flow and sediment transport within the Ning-Meng reach of the Yellow River,northern China

Effective management of a river reach requires a sound understanding of flow and sediment transport generated by varying natural and artificial runoff conditions. Flow and sediment transport within the Ning-Meng reach of the Yellow River（NMRYR）, northern China are controlled by a complex set of factors/processes, mainly including four sets of factors：（1） aeolian sediments from deserts bordering the main stream;（2） inflow of water and sediment from numerous tributaries;（3） impoundment of water by reservoir/hydro-junction; and（4） complex diversion and return of irrigation water. In this study, the 1-D flow ＆ sediment transport model developed by the Yellow River Institute of Hydraulic Research was used to simulate the flow and sediment transport within the NMRYR from 2001 to 2012. All four sets of factors that primarily control the flow and sediment transport mentioned above were considered in this model. Compared to the measured data collected from the hydrological stations along the NMRYR, the simulated flow and sediment transport values were generally acceptable, with relative mean deviation between measured and simulated values of 〈15%. However, simulated sediment concentration and siltation values within two sub-reaches（i.e., Qingtongxia Reservoir to Bayan Gol Hydrological Station and Bayan Gol Hydrological Station to Toudaoguai Hydrological Station） for some periods exhibited relatively large errors（the relative mean deviations between measured and simulated values of 18% and 25%, respectively）. These errors are presumably related to the inability to accurately determine the quantity of aeolian sediment influx to the river reach and the inflow of water from the ten ephemeral tributaries. This study may provide some valuable insights into the numerical simulations of flow and sediment transport in large watersheds and also provide a useful model for the effective management of the NMRYR.

A note on the numerical simulations of flow past a wavy square-section cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Kairmain vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

Application of Stochastic Fracture Network with Numerical Fluid Flow Simulations to Groundwater Flow Modeling in Fractured Rocks

The continuum approach in fluid flow modeling is generally applied to porous geological media, but has limited applicability to fractured rocks. With the presence of a discrete fracture network relatively sparsely distributed in the matrix, it may be difficult or erroneous to use a porous medium fluid flow model with continuum assumptions to describe the fluid flow in fractured rocks at small or even large field scales. A discrete fracture fluid flow approach incorporating a stochastic fracture network with numerical fluid flow simulations could have the capability of capturing fluid flow behaviors such as inhomogeneity and anisotropy while reflecting the changes of hydraulic features at different scales. Moreover, this approach can be implemented to estimate the size of the representative elementary volume (REV) in order to find out the scales at which a porous medium flow model could be applied, and then to determine the hydraulic conductivity tensor for fractured rocks. The following topics are focused on in this study: (a) conceptual discrete fracture fluid flow modeling incorporating a stochastic fracture network with numerical flow simulations; (b) estimation of REV and hydraulic conductivity tensor for fractured rocks utilizing a stochastic fracture network with numerical fluid flow simulations; (c) investigation of the effect of fracture orientation and density on the hydraulic conductivity and REV by implementing a stochastic fracture network with numerical fluid flow simulations, and (d) fluid flow conceptual models accounting for major and minor fractures in the 2 D or 3 D flow fields incorporating a stochastic fracture network with numerical fluid flow simulations.

NUMERICAL SIMULATION OF THE FLOW PAST WING-BODY USING EULER EQUATIONS

AbstractThis paper gives H-O type grid generation and applies finite-volume method forwing-body flow.A 3-D Euler code for wing-body has been developed. We improve in thiscode the distribution of the number of grid points over wing and body surfaces,and keepthe training edge of the wing being one of the grid lines even when the edge has sweepbackor sweep forward.The code developed can provide not only the spanwise pressure distribu-tion as usual ou cross-flow planes,but also the chordwise Pressure distributions.The com-putation for NASA TN D-712 wing-body model with the present code shows that thecomputed pressure distributions are in very good agreement with the experiment.

A numerical simulation of unsteady flow in small diameter helical grafts

The application of small diameter arterial grafts is limited due to the fact of relatively poor long-time patency which is caused by thrombosis formation in the short term and intimal hyperplasia(IH) in the medium and long term.Thrombosis,obstructing the flow of blood

NUMERICAL SIMULATION OF UNSTEADY FLOW AROUND A CIRCULAR CYLINDER AT HIGH REYNOLDS NUMBERS

A Lagrangian-Eulerian hybrid scheme to solve unsteady N-S equation in two-dimensional incompressible fluid at high Reynolds numbers is presented in this paper. A random walk is imposed to simulate the viscous diffusion, the vortex-in-cell method is used to obtain the convection velocity, and nascent vortices are created on a cylinder to satisfy the zero-slip condition. The impulsively started flow around a circular cylinder and the separation induced by a pair of incident vortices symmetrically approaching a circular cylinder have been successfully simulated by the hybrid scheme. The impulsively started flow from rest has been computed at Reynolds numbers 3000 and 9500. Comparisons are made with those results of finite-difference method, vortex method and flow visualization. Agreement is good. The particular attention has been paid to the evolutions of flow pattern. A topological analysis has been proposed in the region of the near wake. The bulge, isolated secondary vortex, a pair of secondary vortices, ' forewake phenomenon and other patterns are simulated numerically. The separation induced by a pair of incident vortices approaching a circular cylinder has been investigated by using the same scheme. The rebounding phenomenon of the incident vortex is observed and is attributed to the effect of the secondary vortex. In particular, we have found that a tertiary vortex can be formed near the surface; this phenomenon has been verified by flow visualization reported recently.

THREE-DIMENSIONAL COUPLED IMPELLER-VOLUTE SIMULATION OF FLOW IN CENTRIFUGAL PUMP AND PERFORMANCE PREDICTION

A three-dimensional turbulent flow through an entire centrifugal pump is simulated using k-ε turbulence model modified by rotation and curvature, SIMPLEC method and body-fitted coordinate. The velocity and pressure fields are obtained for the pump under various working conditions, which is used to predict the head and hydraulic efficiency of the pump, and the results correspond well with the measured values. The calculation results indicate that the pressure is higher on the pressure side than that on the suction side of the blade; The relative velocity on the suction side gradually decreases from the impeller inlet to the outlet, while increases on the pressure side, it finally results in the lower relative velocity on the suction side and the higher one on the pressure side at the impeller outlet; The impeller flow field is asymmetric, i.e. the velocity and pressure fields arc totally different among all channels in the impeller; In the volute, the static pressure gradually increases with the flow route, and a large pressure gratitude occurs in the tongue; Secondary flow exists in the rear part of the spiral.

DEM simulation of particle flow on a single deck banana screen

A mathematical study of particle flow on a banana screen deck using the discrete element method (DEM) was presented in this paper. The motion characteristics and penetrating mechanisms of particles on the screen deck were studied. Effects of geometric parameters of screen deck on banana screening process were also investigated. The results show that when the values of inclination of discharge and increment of screen deck inclination are 10° and 5° respectively, the banana screening process get a good screening performance in the simulation. The relationship between screen deck length and screening efficiency was further confirmed. The conclusion that the screening efficiency will not significantly increase when the deck length L≥430 mm (L/B ≥ 3.5) was obtained, which can provide theoretical basis for the optimization of banana screen.

Numerical Investigation of Flow Motion and Performance of A Horizontal Axis Tidal Turbine Subjected to A Steady Current

Horizontal axis tidal turbines have attracted more and more attentions nowadays, because of their convenience and low expense in construction and high efficiency in extracting tidal energy. The present study numerically investigates the flow motion and performance of a horizontal axis tidal turbine with a supporting vertical cylinder under steady current. In the numerical model, the continuous equation and incompressible Reynolds-averaged Navier-Stokes equations are solved, and the volume of fluid method is employed to track free surface motion. The RNG k-ε model is adopted to calculate turbulence transport while the fractional area/volume obstacle representation method is used to describe turbine characteristics and movement. The effects of installation elevation of tidal turbine and inlet velocity on the water elevation, and current velocity, rotating speed and resultant force on turbine are discussed. Based on the comparison of the numerical results, a better understanding of flow structure around horizontal axis tidal turbine and turbine performance is achieved.

NUMERICAL STUDY OF FLOW IN CONICAL DIFFUSER WITH VORTEX GENERATOR JETS

To develop vortex generator jet （VGJ） method for flow control, the turbulence flow in a 14° conical diffuser with and without vortex generator jets are simulated by solving Navier-Stokes equations with k-ε turbulence model. The diffuser performance, based on different velocity ratio （ratio of the jet speed to the mainstream velocity）, is investigated and compared with the experimental study. On the basis of the flow characteristics using computation fluid dynamics （CFD） method observed in the conical diffuser and the downstream development of the longitudinal vortices, attempt is made to correlate the pressure recovery coefficient with the behavior of vortices produced by vortex generator jets.

Characteristics and generation of elastic turbulence in a three-dimensional parallel plate channel using direct numerical simulation

Direct numerical simulations（DNSs） of purely elastic turbulence in rectilinear shear flows in a three-dimensional（3D） parallel plate channel were carried out,by which numerical databases were established.Based on the numerical databases,the present paper analyzed the structural and statistical characteristics of the elastic turbulence including flow patterns,the wall effect on the turbulent kinetic energy spectrum,and the local relationship between the flow motion and the microstructures＇ behavior.Moreover,to address the underlying physical mechanism of elastic turbulence,its generation was presented in terms of the global energy budget.The results showed that the flow structures in elastic turbulence were 3D with spatial scales on the order of the geometrical characteristic length,and vortex tubes were more likely to be embedded in the regions where the polymers were strongly stretched.In addition,the patterns of microstructures＇ elongation behave like a filament.From the results of the turbulent kinetic energy budget,it was found that the continuous energy releasing from the polymers into the main flow was the main source of the generation and maintenance of the elastic turbulent status.

Numerical Study on the Hydrodynamic Characteristics of Biofouled Full-Scale Net Cage

The effect of biofouling on the hydrodynamic characteristics of the net cage is of particular interest as biofouled nettings can significantly reduce flow of well-oxygenated water reaching the stocked fish. For computational efficiency, the porous-media fluid model is proposed to simulate flow through the biofouled plane net and full-scale net cage. The porous coefficients of the porous-media fluid model can be determined from the quadratic-function relationship between the hydrodynamic forces on a plane net and the flow velocity using the least squares method. In this study, drag forces on and flow fields around five plane nets with different levels of biofouling are calculated by use of the proposed model. The numerical results are compared with the experimental data of Swift et al.（2006） and the effectiveness of the numerical model is presented. On that basis, flow through full-scale net cages with the same level of biofouling as the tested plane nets are modeled. The flow fields inside and around biofouled net cages are analyzed and the drag force acting on a net cage is estimated by a control volume analysis method. According to the numerical results, empirical formulas of reduction in flow velocity and load on a net cage are derived as function of drag coefficient of the corresponding biofouled netting.

Uniformity evaluation and optimization of fluid flow characteristics in a seven-strand tundish

The effect of flow control devices（FCDs） on the uniformity of flow characteristics in a seven-strand symmetrical trapezoidal tundish was studied using both an experimental 1：2.5 hydraulic model and a numerical simulation of a 1：1 geometric model.The variation coefficient（CV） was defined to evaluate the flow uniformity of the seven-strand tundish.An optimized FCD configuration was proposed on the basis of the evaluation of experimental results.It is concluded that a turbulence inhibitor（TI） and U-type dam are essential to improve the uniformity of fluid flow in the seven-strand tundish.In addition,the configuration of inclination T-type dams with a height of 200 mm between the second and third strands and with a height of 300 mm between the third and fourth strands can minimize the proportion of dead zone.After optimizing the configuration of FCDs,the variation coefficient reduces below 20%of the mean value,and the average proportion of dead zone is just 14.6%;in addition,the temperature fluctuation between the strands could be controlled within 0.6 K.In summary,the uniformity of flow and temperature in the seven-strand tundish is greatly improved.

Numerical simulation of flow hydrodynamics of struvite pellets in a liquid–solid fluidized bed

Phosphorus recovery in the form of struvite has been aroused in recent decades for its dual advantages in eutrophication control and resource protection.The usage of the struvite products is normally determined by the size which is largely depended on the hydrodynamics.In this study,flow behavior of struvite pellets was simulated by means of Eulerian–Eulerian two-fluid model combining with kinetic theory of granular flow in a liquid–solid fluidized bed reactor（FBR）.A parametric study including the mesh size,time step,discretization strategy,turbulent model and drag model was first developed,followed by the evaluations of crucial operational conditions,particle characteristics and reactor shapes.The results showed that a cold model with the mesh resolution of 16 × 240,default time step of 0.001 sec and first order discretization scheme was accurate enough to describe the fluidization.The struvite holdup profile using Syamlal–O＇Brien drag model was best fitted to the experimental data as compared with other drag models and the empirical Richardson–Zaki equation.Regarding the model evaluation,it showed that liquid velocity and particle size played important roles on both solid holdups and velocities.The reactor diameter only influenced the solid velocity while the static bed height almost took no effect.These results are direct and can be applied to guide the operation and process control of the struvite fluidization.Moreover,the model parameters can also be used as the basic settings in further crystallization simulations.