UNIFIED COMPUTATION OF FLOW WITH COMPRESSIBLE AND INCOMPRESSIBLE FLUID BASED ON ROE''''S SCHEME

A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe's type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe's scheme.     

MHD flow and heat transfer of micropolar fluid between two porous disks

A numerical study is carried out for the axisymmetric steady laminar incompressible flow of an electrically conducting micropolar fluid between two infinite parallel porous disks with the constant uniform injection through the surface of the disks. The fluid is subjected to an external transverse magnetic field. The governing nonlinear equations of motion are transformed into a dimensionless form through von Karman’s similarity transformation. An algorithm based on a finite difference scheme is used to solve the reduced coupled ordinary differential equations under associated boundary conditions. The effects of the Reynolds number, the magnetic parameter, the micropolar parameter, and the Prandtl number on the flow velocity and temperature distributions are discussed. The results agree well with those of the previously published work for special cases. The investigation predicts that the heat transfer rate at the surfaces of the disks increases with the increases in the Reynolds number, the magnetic parameter, and the Prandtl number. The shear stresses decrease with the increase in the injection while increase with the increase in the applied magnetic field. The shear stress factor is lower for micropolar fluids than for Newtonian fluids, which may be beneficial in the flow and thermal control in the polymeric processing.     

UNSTEADY／STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.     

Slip effects on shearing flows in a porous medium

This paper deals with the magnetohydrodynamic （MHD） flow of an Oldroyd 8-constant fluid in a porous medium when no-slip condition is no longer valid. Modified Darcy＇s law is used in the flow modelling. The non-linear differential equation with non-linear boundary conditions is solved numerically using finite difference scheme in combination with an iterative technique. Numerical results are obtained for the Couette, Poiseuille and generalized Couette flows. The effects of slip parameters on the velocity profile are discussed.     

Numerical simulation of flow in Hartmann resonance tube and flow in ultrasonic gas atomizer

The gas flow in the Hartmann resonance tube is numerically investigated by the finite volume method based on the Roe solver.The oscillation of the flow is studied with the presence of a needle actuator set along the nozzle axis.Numerical results agree well with the theoretical and experimental results available.Numerical results indicate that the resonance mode of the resonance tube will switch by means of removing or adding the actuator.The gas flow in the ultrasonic gas atomization (USGA) nozzle is also studied by the same numerical methods.Oscillation caused by the Hartmann resonance tube structure,coupled with a secondary resonator,in the USGA nozzle is investigated.Effects of the variation of parameters on the oscillation are studied.The mechanism of the transition of subsonic flow to supersonic flow in the USGA nozzle is also discussed based on numerical results.     

Numerical investigations of asymmetric flow of a micropolar fluid between two porous disks

Numerical solution is presented for the two- dimensional flow of a micropolar fluid between two porous coaxial disks of different permeability for a range of Reynolds number Re （-300≤ Re 〈 0） and permeability parameter A （1.0≤A ≤2.0）. The main flow is superimposed by the injection at the surfaces of the two disks. Von Karman＇s similarity transformations are used to reduce the governing equations of motion to a set of non-linear coupled ordinary differential equations （ODEs） in dimensionless form. An algorithm based on the finite difference method is employed to solve these ODEs and Richardson＇s extrapolation is used to obtain higher order accuracy. The results indicate that the parameters Re and A have a strong influence on the velocity and microrotation profiles, shear stresses at the disks and the position of the viscous/shear layer. The micropolar material constants cl, c2, c3 have profound effect on microrotation as compared to their effect on streamwise and axial velocity profiles. The results of micropolar fluids are compared with the results for Newtonian fluids.     

MHD non-Newtonian micropolar fluid flow and heat transfer in channel with stretching walls

A study is presented for magnetohydrodynamics （MHD） flow and heat transfer characteristics of a viscous incompressible electrically conducting micropolar fluid in a channel with stretching walls. The micropolar model introduced by Eringen is used to describe the working fluid. The transformed self similar ordinary differential equations together with the associated boundary conditions are solved numerically by an algorithm based on quasi-linearization and multilevel discretization. The effects of some physical parameters on the flow and heat transfer are discussed and presented through tables and graphs. The present investigations may be beneficial in the flow and thermal control of polymeric processing.     

A study of the upper limit of solid scatters density for gray Lattice Boltzmann Method

The upper limit of the solid scatters density ns （x）, a key parameter for the simulation of flows in porous media with a gray Lattice Boltzmann Method, is studied by an analytical way for the infiltration Poiseuille flow between two infinite parallel plates. Analyses of three different gray Lattice Boltzmann schemes, separately proposed by Gao and Sharma et al., Dardis and McCloskey, and Thorne and Sukop, indicate that the effective domain of Gao and Sharma＇s scheme is restricted to ns 〈 1/2√3≈0.289, Dardis and McCloskey＇s scheme is restricted to ns 〈 （√57-1）/28≈0.234, and that there is no extra restriction on ns（x） with Thorne and Sukop＇s scheme. These results are obtained for the dimensionless relaxation time τ= 1. The above analytical results are verified by our numerical simulations. The use of a gray LBM is further illustrated by simulating the flow at the interface of a porous medium. Simulation results yield velocity profiles which agree very well with Brinkman＇s prediction.     

Robust design of natural laminar flow supercritical airfoil by multi-objective evolution method

A transonic, high Reynolds number natural laminar flow airfoil is designed and studied. The γ-θ transition model is combined with the shear stress transport(SST)k-w turbulence model to predict the transition region for a laminar-turbulent boundary layer. The non-uniform free-form deformation(NFFD) method based on the non-uniform rational B-spline(NURBS) basis function is introduced to the airfoil parameterization.The non-dominated sorting genetic algorithm-II(NSGA-II) is used as the search algorithm, and the surrogate model based on the Kriging models is introduced to improve the efficiency of the optimization system. The optimization system is set up based on the above technologies, and the robust design about the uncertainty of the Mach number is carried out for NASA0412 airfoil. The optimized airfoil is analyzed and compared with the original airfoil. The results show that natural laminar flow can be achieved on a supercritical airfoil to improve the aerodynamic characteristic of airfoils.     

Formulation of thermo-hydro-mechanical coupling behavior of unsaturated soils based on hybrid mixture theory

Thermo-Hydro-Mechanical （THM） coupling pro- cesses in unsaturated soils are very important in both theoretical researches and engineering applications. A coupled formulation based on hybrid mixture theory is derived to model the THM coupling behavior of unsaturated soils. The free-energy and dissipative functions for different phases are derived from Taylor＇s series expansions. Constitutive relations for THM coupled behaviors of unsaturated soils, which include deformation, entropy change, fluid flow, heat conduction, and dynamic compatibility conditions on the interfaces, are then established. The number of field equations is shown to be equal to the number of unknown variables; thus, a closure of this coupling problem is established. In addition to modifications of the physical conservation equations with coupling effect terms, the constitutive equations, which consider the coupling between elastoplastic deformation of the soil skeleton, fluid flow, and heat transfer, are also derived.     

Mixed convection stagnation-point flow on vertical stretching sheet with external magnetic field

The problem of steady laminar magnetohydrodynamic （MHD） mixed con- vection stagnation-point flow of an incompressible viscous fluid over a vertical stretch- ing sheet is studied. The effect of an externally magnetic field is taken into account. The transformed boundary layer equations are solved numerically by using an implicit finite-difference scheme. Numerical results are obtained for various values of the mixed convection parameter, Hartmann number, and Prandtl number. The effects of an exter- nally magnetic field on the skin friction coefficient, local Nusselt number, velocity, and temperature profiles for both A 〉 1 and A ~ 1, where A is the velocity ratio parameter, are presented graphically and discussed in detail. Both assisting and opposing flows are considered, and it is found that dual solutions exist for the opposing flow.     

Investigation of the turbulent model for pressure fluctuations with spectral theory

The pressure fluctuations in turbulent shear flows are investigated with the theory of spectral analysis.An expression for pressure spectra is analytically derived in terms of velocity spectra.This derivation is based on a formal solution of the Navier-Stokes equation and quasi-normal assumption to express the third and fourth order velocity correlations in terms of double velocity correlation.Then,a turbulent model for the computation of pressure fluctuation intensity with Renolds stress and mean flow velocity gradients is established.The turbulent constants in the model are calculated from the assumptions about the general behaviour of velocity spectra in high Renolds number flows.Comparison with direct simulation of turbulent boundary layer is made.It is found that the turbulent-turbulent,cross correlation,and turbulent-shear source terms for mean square value of pressure fluctuation are about the same magnitude.     

Regimes of vortex shedding from an in-line oscillating circular cylinder in the uniform flow

A finite volume method for the time dependent viscous incompressible flow around an in-line oscillating circular cylinder at Reynolds number of 200, 855 is presented in this paper. The Navier-Stokes equations in a finite volume form are solved with a moving grid system, based on a time dependent coordinate transformation. To investigate the vortex-shedding characteristics behind the circular cylinder and the effects of Reynolds number and other non-dimensional parameters such as reduced amplitude and reduced frequency， several numerical schemes have been tested with different amplitude and frequency close to Sto and a harmonic at each Reynolds number. Present numerical results indicate several types of vortex shedding mode which is known mainly depending on the reduced frequency and also the reduced amplitude, which is called synchronization or lock-on.     

A gas-kinetic theory based multidimensional high-order method for the compressible Navier–Stokes solutions

This paper presents a gas-kinetic theory based multidimensional high-order method for the compressible Naiver–Stokes solutions. In our previous study, a spatially and temporally dependent third-order flux scheme with the use of a third-order gas distribution function is employed.However, the third-order flux scheme is quite complicated and less robust than the second-order scheme. In order to reduce its complexity and improve its robustness, the secondorder flux scheme is adopted instead in this paper, while the temporal order of method is maintained by using a two stage temporal discretization. In addition, its CPU cost is relatively lower than the previous scheme. Several test cases in two and three dimensions, containing high Mach number compressible flows and low speed high Reynolds number laminar flows, are presented to demonstrate the method capacity.     

Fully coupled flow-induced vibration of structures under small deformation with GMRES method

Lagrangian-Eulerian formulations based on a generalized variational principle of fluid-solid coupling dynamics are established to describe flow-induced vibration of a structure under small deformation in an incompressible viscous fluid flow. The spatial discretization of the formulations is based on the multi-linear interpolating functions by using the finite element method for both the fluid and solid structures. The generalized trapezoidal rule is used to obtain apparently non-symmetric linear equations in an incremental form for the variables of the flow and vibration. The nonlinear convective term and time factors are contained in the non-symmetric coefficient matrix of the equations. The generalized minimum residual （GMRES） method is used to solve the incremental equations. A new stable algorithm of GMRES-Hughes-Newmark is developed to deal with the flow-induced vibration with dynamical fluid-structure interaction in complex geometries. Good agreement between the simulations and laboratory measurements of the pressure and blade vibration accelerations in a hydro turbine passage was obtained, indicating that the GiViRES-Hughes-Newmark algorithm presented in this paper is suitable for dealing with the flow-induced vibration of structures under small deformation.     

A SOLVER USING NEWTON‘S METHOD FOR UNSTEADY VISCOUS FLOWS

A solver is developed for time-accurate computations of viscous flows based on the conception of Newton‘s method. A set of pseudo-time derivatives are added into governing equations and the discretized system is solved using GMRES algorithm. Due to some special properties of GMRES algorithm, the solution procedure for unsteady flows could be regarded as a kind of Newton iteration. The physical-time derivatives of governing equations are discretized using two different approaches, I.e., 3-point Euler backward, and Crank-Nicolson formulas, both with 2nd-order accuracy in time but with different truncation errors. The turbulent eddy viscosity is calculated by using a version of Spalart~Allmaras one-equation model modified by authors for turbulent flows. Two cases of unsteady viscous flow are investigated to validate and assess the solver, I.e., low Reynolds number flow around a row of cylinders and transonic bi-circular-arc airfoil flow featuring the vortex shedding and shock buffeting problems, respectively. Meanwhile, comparisons between the two schemes of timederivative discretizations are carefully made. It is illustrated that the developed unsteady flow solver shows a considerable efficiency and the Crank-Nicolson scheme gives better results compared with Euler method.     

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.     

Perturbation solution to unsteady flow in a porous channel with expanding or contracting walls in the presence of a transverse magnetic field

An incompressible flow in a porous channel with expanding or contacting walls in the presence of a transverse magnetic field is considered. Using similarity transformations, the governing equations are reduced to the nonlinear ordinary differential equations. The exact similar solutions for the different cases of the expansion ratio and the Hartmann number are obtained with a singular perturbation method, and the associated behavior is discussed in detail.     

ON GALERKIN DISCRETIZATION OF AXIALLY MOVING NONLINEAR STRINGS

A computational technique is proposed for the Galerkin discretization of axially moving strings with geometric nonlinearity. The Galerkin discretization is based on the eigenfunctions of stationary strings. The discretized equations are simplified by regrouping nonlinear terms to reduce the computation work. The scheme can be easily implemented in the practical programming. Numerical results show the effectiveness of the technique. The results also highlight the feature of Galerkin＇s discretization of gyroscopic continua that the term number in Galerkin＇s discretization should be even. The technique is generalized from elastic strings to viscoelastic strings.     

Effects of temperature-dependent viscosity variation on entropy generation,heat and fluid flow through a porous-saturated duct of rectangular cross-section

Effect of temperature-dependent viscosity on fully developed forced convection in a duct of rectangular cross-section occupied by a fluid-saturated porous medium is investigated analytically. The Darcy flow model is applied and the viscosity-temperature relation is assumed to be an inverse-linear one. The case of uniform heat flux on the walls, i.e. the H boundary condition in the terminology of Kays and Crawford [12], is treated. For the case of a fluid whose viscosity decreases with temperature, it is found that the effect of the variation is to increase the Nusselt number for heated walls. Having found the velocity and the temperature distribution, the second law of thermodynamics is invoked to find the local and average entropy generation rate. Expressions for the entropy generation rate, the Bejan number, the heat transfer irreversibility, and the fluid flow irreversibility are presented in terms of the Brinkman number, the Peclet number, the viscosity variation number, the dimensionless wall heat flux, and the aspect ratio (width to height ratio). These expressions let a parametric study of the problem based on which it is observed that the entropy generated due to flow in a duct of square cross-section is more than those of rectangular counterparts while increasing the aspect ratio decreases the entropy generation rate similar to what previously reported for the clear flow case by Ratts and Raut [14].