A CFD Based Investigation of the Unsteady Hydrodynamic Coefficients of 3-D Fins in Viscous Flow

The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly unsteady, generating vortices and requiring detailed analysis of fluid-structure interactions.An understanding of the complexities of such flows is of interest to engineers developing vehicles capable of high dynamic performance in their propulsion and maneuvering.In the present study, a CFD based RANS simulation of a 3-D fin body moving in a viscous fluid was developed.It investigated hydrodynamic performance by evaluating the hydrodynamic coefficients (lift, drag and moment) at two different oscillating frequencies.A parametric analysis of the factors that affect the hydrodynamic performance of the fin body was done, along with a comparison of results from experiments.The results of the simulation were found in close agreement with experimental results and this validated the simulation as an effective tool for evaluation of the unsteady hydrodynamic coefficients of 3-D fins.This work can be further be used for analysis of the stability and maneuverability of fin actuated underwater vehicles.

Conical Sonic-Supersonic Solutions for the 3-D Steady Full Euler Equations

This paper concerns the sonic-supersonic structures of the transonic crossflow generated by the steady supersonic flow past an infinite cone of arbitrary cross section.Under the conical assumption,the three-dimensional(3-D)steady Euler equations can be projected onto the unit sphere and the state of fluid can be characterized by the polar and azimuthal angles.Given a segment smooth curve as a conical-sonic line in the polar-azimuthal angle plane,we construct a classical conical-supersonic solution near the curve under some reasonable assumptions.To overcome the difficulty caused by the parabolic degeneracy,we apply the characteristic decomposition technique to transform the Euler equations into a new degenerate hyperbolic system in a partial hodograph plane.The singular terms are isolated from the highly nonlinear complicated system and then can be handled successfully.We establish a smooth local solution to the new system in a suitable weighted metric space and then express the solution in terms of the original variables.

Simulation of subcritical flow pattern in 180° uniform and convergent open-channel bends using SSIIM 3-D model

In meandering rivers, the flow pattern is highly complex, with specific characteristics at bends that are not observed along straight paths. A numerical model can be effectively used to predict such flow fields. Since river bends are not uniform-some are divergent and others convergent-in this study, after the SSIIM 3-D model was calibrated using the result of measurements along a uniform 180° bend with a width of 0.6 m, a similar but convergent 180v bend, 0.6 m to 0.45 m wide, was simulated using the SSI1M 3-D numerical model. Flow characteristics of the convergent 180° bend, including lengthwise and vertical velocity profiles, primary and secondary flows, lengthwise and widtbwise slopes of the water surface, and the helical flow strength, were compared with those of the uniform 180° bend. The verification results of the model show that the numerical model can effectively simulate the flow field in the uniform bend. In addition, this research indicates that, in a convergent channel, the maximum velocity path at a plane near the water surface crosses the channel＇s centerline at about a 30° to 40° cross-section, while in the uniform bend, this occurs at about the 50° cross-section. The varying range of the water surface elevation is wider in the convergent channel than in the uniform one, and the strength of the helical flow is generally greater in the uniform channel than in the convergent one. Also, unlike the uniform bend, the convergent bend exhibits no rotational cell against the main direction of secondary flow rotation at the 135° cross-section.

EFFECT OF SUBMERGED ENTRY NOZZLE (SEN) PARAMETERS AND SHAPE ON 3-D FLUID FLOW IN MOULD FOR BEAM BLANK CONTINUOUS CASTING

According to turbulent theory and characteristics of beam blank continuous casting, 3-D model to represent the flow of beam blank mould is established. The predicted results indicate that the exit obliquity of up 15°(+15°) should be adopted, which will benefit the floatation of non-metallic inclusion and purification of the molten steel. When the nozzle angle is 120°, the flow pattern is reasonable. Proper nozzle depth can be 200mm. Turbulent kinetic of meniscus can be reduced by adopting the square nozzle and suitable area of side outlet when casting speed increases. The results are consistent with those of water model experiment, so the model is exact and reasonable. The model can provide important information for design of SEN and defining of immersion depth.

Predication of 3-D Viscous Flowfield of a Centrifugal Impeller

A three-dimensional viscous code has been developed to solve Reynolds-averaged Navier-Stokes equations. The governing equations in finite volume form are solved by two-step Runge-Kutta scheme with implicit residual smoothing. The eddy viscous is obtained using the Baldwin-Lomax model. A prediction of the 3-D turbulent flow and the performance in the “all-over controlled vortex distribution” centrifugal impeller with a vaneless diffuser has been made for the compressor at design and off-design condition. The predicted effi-ciency is a little higher than the experiment data. These results suggest that the present calculation code is able to determine the flow development in the impeller and also the turbulence model in the centrifugal im-peller should be improved.

A Higher-Efficient Non-Hydrostatic Finite Volume Model for Strong Three-Dimensional Free Surface Flows and Sediment Transport

In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional （2-D） grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.

3-D Finite-Element Numerical Simulation of Centrifugal Induction Electrosalg Casting Solidification Process

A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.

Analysis of Combined Radiation and Forced Convection Heat Transfer in 3D Laminar Flow over an Inclined Forward Facing Step

In the current study, a numerical investigation of three-dimensional combined convection-radiation heat transfer over an inclined forward facing step (FFS) in a horizontal rectangular duct is presented. The fluid is treated as a gray, absorbing, emitting and scattering medium. To simulate the incline surface of FFS, the blocked-off method is employed in this study. The set of governing equations for gas flow are solved numerically using the CFD technique to obtain the temperature and velocity fields. Since the gas is considered as a radiating medium, all of the convection, conduction and radiation heat transfer mechanisms are presented in the energy equation. For computation of radiative term in energy equation, the radiative transfer equation (RTE) is solved numerically by the discrete ordinates method (DOM) to find the divergence of radiative heat flux distribution inside the radiating medium. The effects of optical thickness, radiation-conduction parameter and albedo coefficient on heat transfer behavior of the system are carried out.

Formation of the Tibetan Plateau during the India-Eurasia Convergence:Insight from 3-D Multi-Terrane Thermomechanical Modeling

Various models exist to explain the formation of the Tibetan Plateau,including“tectonic escape”,“pure shear thickening”,“convective removal of the lithospheric mantle”,and“lower crustal flow”model.The first two models are primarily constructed on pure mechanical models but are unable to reasonably explain the tension and shear phenomena inside the plateau.The latter two are rheological dynamic models based on deep geophysical observations.However,the spatial range of the lower crustal flow and its role in the plateau formation/uplift remain controversial.Five multi-terrane viscoplastic thermomechanical models were constructed to simulate the uplift and lithospheric structure change of the Tibetan Plateau during the post-collision stage(since 35 Ma)under the convergence of the Indian Plate.Results show that the plateau's formation begins with crustal thickening,blocked by strong terranes at the northern plateau,and expanded laterally to the east.The lithosphere thickens gradually and experiences delamination at its base,elevating temperature within the crust and forming partial melting layers in the central plateau.As convergence persists on the southern side,the northern plateau's lithosphere bends downward and undergoes delamination,further heating the crust and promoting the northward and eastward flow of partial melting layers,leading to secondary uplift around the plateau.

Numerical simulation and analysis of solid-liquid two-phase threedimensional unsteady flow in centrifugal slurry pump

Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.

Simulation of Polymer Melt Flow Fields in Intermeshing Co-Rotating Three-Screw Extruders

Three-dimension isothermal flows of polymer melt in the kneading blocks of triangularly-arranged and parallelly-arranged intermeshing co-rotating three-screw extruders are simulated using the finite element package POLYFLOW. Based on the velocity fields calculated, the particle trajectories in both machines are visualized using particle tracking technique. The numerical results indicate that the flow patterns in three-screw extruders are similar to those in twin-screw extruders. The triangularly-arranged three-screw extruder has the largest pumping capacity and also the highest extrusion stability in terms of flowrate fluctuation with screw rotation. The instantaneous mixing and cumulative residence time distribution （RTD） characteristics are also analyzed and compared with traditional intermeshing co-rotating twin-screw extruders. It is shown that the start section of the cumulative RTD curve for the triangularly-arranged machine has a small shoulder, which is attributed to the faster flow in the central region of this type of extruder.

Investigation of semi-solid aluminum alloy filling-plastic flowing in thixoforging

The semi-solid filling-plastic flowing integrated forging process of semi-solid 6061 Al alloy was simulated by commercial finite element software DEFORM-3D.Temperature,fluid and stress-strain fields were considered in numerical simulation.The simulation results show that the plastic deformation of billet of the ends is higher than that of billet in the straight cylinder.The value of plastic deformation varies with loading mode and plastic deformation fields at the stage of increasing pressure to constant value.When the thixoforging experiments were performed at 590 ℃,15 mm/s of punch velocity and 46 MPa of pressure side urn,it gets the filling wholly and dense internal organization of semi-solid thixoforging parts is gotten.Finite element analysis results are compatible with experimental ones.

A New Efficient Finite Volume Modeling of Small Amplitude Free Surface Flows with Unstructured Grid

A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.

Corner Flow Control in High Through-Flow Axial Commercial Fan/Booster Using Blade 3-D Optimization

This study is aimed at using blade 3-D optimization to control corner flows in the high through-flow fan/booster of a high bypass ratio commercial turbofan engine. Two kinds of blade 3-D optimization, end-bending and bow, are focused on. On account of the respective operation mode and environment, the approach to 3-D aerodynamic modeling of rotor blades is different from stator vanes. Based on the understanding of the mechanism of the corner flow and the consideration of intensity problem for rotors, this paper uses a variety of blade 3-D optimization approaches, such as loading distribution optimization, perturbation of departure angles and stacking-axis manipulation, which are suitable for rotors and stators respectively. The obtained 3-D blades and vanes can improve the corner flow features by end-bending and bow effects. The results of this study show that flows in corners of the fan/booster, such as the fan hub region, the tip and hub of the vanes of the booster, are very complex and dominated by 3-D effects. The secondary flows there are found to have a strong detrimental effect on the compressor performance. The effects of both end-bending and bow can improve the flow separation in corners, but the specific ways they work and application scope are somewhat different. Redesigning the blades via blade 3-D optimization to control the corner flow has effectively reduced the loss generation and improved the stall margin by a large amount.

Solving generalized lattice Boltzmann model for 3-D cavity flows using CUDA-GPU

The generalized lattice Boltzmann equation(GLBE),with the addition of the standard Smagorinsky subgrid-stress(SGS) model,has been proved that it is more suitable for simulating high Reynolds number turbulent flows when compared with the lattice BGK Boltzmann equation(LBGK).However,the computing efficiency of lattice Boltzmann method(LBM) is too low to make it for practical applications,unless using a massive parallel computing clusters facility.In this study,the massive parallel computing power from an inexpensive graphic processor unit(GPU) and a typical personal computer has been developed for improving the computing efficiency,more than 100 times.This developed three-dimensional(3-D) GLBE-SGS model,with the D3Q19 scheme for simplifying collision and streaming courses,has been successfully used to study 3-D rectangular cavity flows with Reynolds number up to 10000.

APPLICATION OF VIDEO ANAGLYPH MAKER FOR 3-D FLOW SIMULATION

In this article,the anaglyph video maker is employed for generating realistic 3-D flows and the software FlowAnimator is developed using that technology.Based on Microsoft Windows Presentation Foundation（WPF）,the real 3-D scene is set up and marker particles are distributed in it randomly in order to create a more natural flow scenario.The trajectory of the particle motion is calculated with Lagrangian description in 3-D space.During the simulation,the viewport can be changed in order to focus on different parts of the model by panning,zooming,rotating and inclination variation etc.Marker particles may appear in different shapes：spheres,tracking-balls,cylinders and ribbons in order to fit different flows.It is the first time that the video anaglyph technology is employed in the 3-D hydrodynamic simulation,which removes the obstacles for 3-D scenes to be rendered on a flat-panel display.

3-D simulation of transient flow patterns in a corridor-shaped air-cushion surge chamber based on computational fluid dynamics

The 3-D characteristics of the water-air flow patterns in a corridor-shaped air-cushion surge chamber during hydraulic transients need to be considered in the shape optimization. To verify the reliability of the water-air two-phase model, namely, the volume of fluid model, the process of charging water into a closed air chamber is successfully simulated. Using the model, the 3-D flow characteristics under the load rejection and acceptance conditions within the air-cushion surge chamber of a specific hydro- power station are studied. The free surface waves, the flow patterns, and the pressure changes during the surge wave process are analyzed in detail. The longitudinal flow of water in the long corridor-shaped surge chamber is similar to the open channel flow with respect to the wave propagation, reflection and superposition characteristics. The lumped parameters of the 3-D numerical simulation agree with the results of a 1-D calculation of hydraulic transients in the whole water conveying system, which validates the 3-D method. The 3-D flow structures obtained can be applied to the shape optimization of the chamber.

EXPERIMENTAL STUDY OF 3-D TURBULENT BEND FLOWS IN OPEN CHANNEL

A generalized hend flow model, treating a 90° single bend and 60° continuous hends, was designed to quantitatively describe 3-D turhulenee mechanism of circulating notfully-developed flow in open channels with hends. The 3-D fluctuating veloeities of turbulent flow were measured and analyzed with a 3 D acoustic-Doppler velocimeter. Formula for 3 D turbulent intensity was derived using the dimension analysis approaeh. Expressions of vertical turbulent intensity distributions were obtained with the multivariant-rcgression theo ry, whieh agree with experiment data. Distrihutions of turbulent intensity and turbulent stress were characterized, and their relationships were concluded. In the bend-turbulent flow core region, longitudinal and lateral turbulent-intensity distri hutions are coincident with linear distribution, hut in nearwall region are coincident with the Gamma distribution. Verotical turbulent intensity distributions are coincident with the Rayleigh distribution. Herein, it is concluded that the bend turbulence is anisotropic.

A NUMERICAL MODEL FOR 3-D FLOW IN OPEN CHANNELS

A 3-D model based on the Reynolds equations with closed k-ε turbulence model is presented in this paper,which can be used to predict surface water flow in open channels.In- stead of the“rigid lid”approximation,the solution of the free surface equation is implemented in the velocity-pressure iterative procedure on the basis of the conventional SIMPLE method.This model was used to compute the flow in rectangular channels with trenches dredged across the bot- tom.The velocity,eddy viscosity coefficient,turbulent shear stress,turbulent kinetic energy and elevation of the free surface over the trenches dredged in the main channel,can be obtained. The computed results are in good agreement with existing experimentaing data.

NUMERICAL COMPUTATION OF 3-D TURBULENT FLOW IN TURBINE CASCADE

Ⅰ. INTRODUCTIONThe flow phenomena that exist in the advanced turbomachinery are extremely complex and propose a challenge to the engineers and scientists to improve the design procedure. The rapid progress of high speed and large capacity computers has encouraged the development of computational fluid dynamics that has had an increasingly important influence on turbomachinery blade design and analysis. Compared with the inviscid solver, the N.S.