A solution for the vertical variation of stress,rather than velocity,in a three-dimensional circulation model

A simple technique is presented that allows a numerical solution to be sought for the vertical variation of shear stress as a substitute for the vertical variation of velocity in a three-dimensional hydrodynamic model. In its most general form the direct stress solution (DSS) method depends only upon the validity of an eddy viscosity relation between the shear stress and the vertical gradient of velocity. The rationale for preferring a numerical solution for shear stress to one for velocity is that shear stress tends to vary more slowly over the vertical than velocity, particularly near boundaries. Consequently, a numerical solution can be obtained much more efficiently for shear stress than for velocity. When needed, the velocity profile can be recovered from the stress profile by solving a one-dimensional integral equation over the vertical. For most practical problems this equation can be solved in closed form. Comparisons are presented between the DSS technique, the standard velocity solution technique and analytical solutions for wind-driven circulation in an unstratified, closed, rectangular channel governed by the linear equations of motion. In no case was the computational effort required by the velocity solution competitive with the DSS when a physically realistic boundary layer was included. The DSS technique should be particularly beneficial in numerical models of relatively shallow water bodies in which the bottom and surface boundary layers occupy a significant portion of the water column.     

Extension of the local domain-free discretization method to large eddy simulation of turbulent flows

In this work, an immersed boundary method, called the local domain-free discretization (DFD) method, is extended to large eddy simulation (LES) of turbulent flows. The discrete form of partial differential equations at an interior node may involve some nodes outside the solution domain. The flow variables at these exterior dependent nodes are evaluated via linear extrapolation along the direction normal to the wall. To alleviate the requirement of mesh resolution in the near-wall region, a wall model based on the turbulence boundary layer equations is introduced. The wall shear stress yielded by the wall model and the no-penetration condition are enforced at the immersed boundary to evaluate the velocity components at an exterior dependent node. For turbulence closure, a dynamic subgrid scale (SGS) model is adopted and the Lagrangian averaging procedure is used to compute the model coefficient. The SGS eddy viscosity at an exterior dependent node is set to be equal to that at the outer layer. To maintain the mass conservation near the immersed boundary, a mass source/sink term is added into the continuity equation. Numerical experiments on relatively coarse meshes with stationary or moving solid boundaries have been conducted to verify the ability of the present LES-DFD method. The predicted results agree well with the published experimental or numerical data.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, “turbulent spots” appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed near sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces. The project supported by the National Natural Science Foundation of China (Grant No.19672070)     

Unsteady compressible boundary layer flow over a circular cone near a plane of symmetry

An analysis has been performed to study the unsteady laminar compressible boundary layer governing the hypersonic flow over a circular cone at an angle of attack near a plane of symmetry with either inflow or outflow in the presence of suction. The flow is assumed to be steady at time t=0 and at t>0 it becomes unsteady due to the time-dependent free stream velocity which varies arbitrarily with time. The nonlinear coupled parabolic partial differential equations under boundary layer approximations have been solved by using an implicit finite-difference method. It is found that suction plays an important role in stabilising the fluid motion and in obtaining unique solution of the problem. The effect of the cross flow parameter is found to be more pronounced on the cross flow surface shear stress than on the streamwise surface shear stress and surface heat transfer. Beyond a certain value of the cross flow parameter overshoot in the cross flow velocity occurs and the magnitude of this overshoot increases with the cross flow parameter. The time variation of the streamwise surface shear stress is more significant than that of the cross flow surface shear stress and surface heat transfer. The suction and the total enthalpy at the wall exert strong influence on the streamwise and cross flow surface shear stresses and the surface heat transfer except that the effect of suction on the cross flow surface shear stress is small.     

Lie-group method of solution for steady two-dimensional boundary-layer stagnation-point flow towards a heated stretching sheet placed in a porous medium

The boundary-layer equations for two-dimensional steady flow of an incompressible, viscous fluid near a stagnation point at a heated stretching sheet placed in a porous medium are considered. We apply Lie-group method for determining symmetry reductions of partial differential equations. Lie-group method starts out with a general infinitesimal group of transformations under which the given partial differential equations are invariant. The determining equations are a set of linear differential equations, the solution of which gives the transformation function or the infinitesimals of the dependent and independent variables. After the group has been determined, a solution to the given partial differential equations may be found from the invariant surface condition such that its solution leads to similarity variables that reduce the number of independent variables of the system. The effect of the velocity parameter λ, which is the ratio of the external free stream velocity to the stretching surface velocity, permeability parameter of the porous medium k ₁, and Prandtl number Pr on the horizontal and transverse velocities, temperature profiles, surface heat flux and the wall shear stress, has been studied.     

Explicit formulations for evaluation of velocity gradients using boundary‐domain integral equations in 2D and 3D viscous flows

In this paper, explicit boundary‐domain integral equations for evaluating velocity gradients are derived from the basic velocity integral equations. A free term is produced in the new strongly singular integral equation, which is not included in recent formulations using the complex variable differentiation method (CVDM) to compute velocity gradients (Int. J. Numer. Meth. Fluids 2004; 45 :463–484; Int. J. Numer. Meth. Fluids 2005; 47 :19–43). The strongly singular domain integrals involved in the new integral equations are accurately evaluated using the radial integration method (RIM). Considerable computational time for evaluating integrals of velocity gradients can be saved by using present formulation than using CVDM. The formulation derived in this paper together with those presented in reference (Int. J. Numer. Meth. Fluids 2004; 45 :463–484) for 2D and in (Int. J. Numer. Meth. Fluids 2005; 47 :19–43) for 3D problems constitutes a complete boundary‐domain integral equation system for solving full Navier–Stokes equations using primitive variables. Three numerical examples for steady incompressible viscous flow are given to validate the derived formulations. Copyright © 2007 John Wiley & Sons, Ltd.     

Transient Dynamic Response of Elastoplastic Thick Plates Resting on Winkler-Type Foundation

A direct domain/boundary element method (D/BEM) for the dynamicanalysis of elastoplastic thick plates resting on Winkler-typefoundation is presented. The formulation presented includes the effectsof shear deformation and rotatory inertia following Reissner–Mindlin'sdeformation theory assumption. The method employs the elastostaticfundamental solution of the problem resulting in both boundary anddomain integrals due to inertia, plasticity and foundation pressureterms. By discretizing the integral equations and integrating theresulting matrix equation of motion by an explicit time integrationscheme the history of the plate dynamic response can be obtained.Applications are shown for impacted elastoplastic Reissner–Mindlinplates with smooth boundaries on several different foundations.     

Gauge finite element method for incompressible flows

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J‐G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher‐order) finite elements. This method can achieve high‐order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Efficient solution techniques for implicit finite element schemes with flux limiters

The algebraic flux correction (AFC) paradigm is equipped with efficient solution strategies for implicit time‐stepping schemes. It is shown that Newton‐like techniques can be applied to the nonlinear systems of equations resulting from the application of high‐resolution flux limiting schemes. To this end, the Jacobian matrix is approximated by means of first‐ or second‐order finite differences. The edge‐based formulation of AFC schemes can be exploited to devise an efficient assembly procedure for the Jacobian. Each matrix entry is constructed from a differential and an average contribution edge by edge. The perturbation of solution values affects the nodal correction factors at neighbouring vertices so that the stencil for each individual node needs to be extended. Two alternative strategies for constructing the corresponding sparsity pattern of the resulting Jacobian are proposed. For nonlinear governing equations, the contribution to the Newton matrix which is associated with the discrete transport operator is approximated by means of divided differences and assembled edge by edge. Numerical examples for both linear and nonlinear benchmark problems are presented to illustrate the superiority of Newton methods as compared to the standard defect correction approach. Copyright © 2007 John Wiley & Sons, Ltd.     

Solution of the incompressible Navier–Stokes equations on domains with one or several open boundaries1

The aim of this paper is to develop a methodology for solving the incompressible Navier–Stokes equations in the presence of one or several open boundaries. A new set of open boundary conditions is first proposed. This has been developed in the context of the velocity–vorticity formulation, but it is also emphasized how it can be formally extended to the equations in primitive variables. The case of a domain involving several independent open boundaries is considered next. An influence matrix technique is applied such that the inlet mass flux is split onto the several outlets in order to enforce the prescribed mean pressure at each outlet. Both approaches are validated by numerical test cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Three-dimensional incompressible Navier–Stokes equations on non-orthogonal staggered grids using the velocity–vorticity formulation

This paper is concerned with the numerical resolution of the incompressible Navier–Stokes equations in the velocity–vorticity form on non-orthogonal structured grids. The discretization is performed in such a way, that the discrete operators mimic the properties of the continuous ones. This allows the discrete equivalence between the primitive and velocity–vorticity formulations to be proved. This last formulation can thus be seen as a particular technique for solving the primitive equations. The difficulty associated with non-simply connected computational domains and with the implementation of the boundary conditions are discussed. One of the main drawback of the velocity–vorticity formulation, relative to the additional computational work required for solving the additional unknowns, is alleviated. Two- and three-dimensional numerical test cases validate the proposed method. © 1998 John Wiley & Sons, Ltd.     

Boundary element solution of unsteady magnetohydrodynamic duct flow with differential quadrature time integration scheme

A numerical scheme which is a combination of the dual reciprocity boundary element method (DRBEM) and the differential quadrature method (DQM), is proposed for the solution of unsteady magnetohydrodynamic (MHD) flow problem in a rectangular duct with insulating walls. The coupled MHD equations in velocity and induced magnetic field are transformed first into the decoupled time‐dependent convection–diffusion‐type equations. These equations are solved by using DRBEM which treats the time and the space derivatives as nonhomogeneity and then by using DQM for the resulting system of initial value problems. The resulting linear system of equations is overdetermined due to the imposition of both boundary and initial conditions. Employing the least square method to this system the solution is obtained directly at any time level without the need of step‐by‐step computation with respect to time. Computations have been carried out for moderate values of Hartmann number (M?50) at transient and the steady‐state levels. As M increases boundary layers are formed for both the velocity and the induced magnetic field and the velocity becomes uniform at the centre of the duct. Also, the higher the value of M is the smaller the value of time for reaching steady‐state solution. Copyright © 2005 John Wiley & Sons, Ltd.     

A complete boundary integral formulation for steady compressible inviscid flows governed by non-linear equations

A complete boundary integral formulation for steady compressible inviscid flows governed by non-linear equations is established by using the specific mass flux as a dependent variable. Thus, the dimensionality of the problem to be solved is reduced by one and the computational mesh to be generated is needed only on the boundary of the domain. It is shown that the boundary integral formulation developed in this paper is equivalent of the results of distributions of the fundamental solutions of the Laplacian operator equation with a different order along the boundaries of the domain. Hence, we have succeeded in establishing the fundamental-solution method for compressible inviscid flows governed by non-linear equations.     

Unsteady mixed convection on the stagnation-point flow adjacent to a vertical plate with a magnetic field

An analysis is performed to study the unsteady combined forced and free convection flow (mixed convection flow) of a viscous incompressible electrically conducting fluid in the vicinity of an axisymmetric stagnation point adjacent to a heated vertical surface. The unsteadiness in the flow and temperature fields is due to the free stream velocity, which varies arbitrarily with time. Both constant wall temperature and constant heat flux conditions are considered in this analysis. By using suitable transformations, the Navier–Stokes and energy equations with four independent variables (x, y, z, t) are reduced to a system of partial differential equations with two independent variables (, ). These transformations also uncouple the momentum and energy equations resulting in a primary axisymmetric flow, in an energy equation dependent on the primary flow and in a buoyancy-induced secondary flow dependent on both primary flow and energy. The resulting system of partial differential equations has been solved numerically by using both implicit finite-difference scheme and differential-difference method. An interesting result is that for a decelerating free stream velocity, flow reversal occurs in the primary flow after certain instant of time and the magnetic field delays or prevents the flow reversal. The surface heat transfer and the surface shear stress in the primary flow increase with the magnetic field, but the surface shear stress in the buoyancy-induced secondary flow decreases. Further the heat transfer increases with the Prandtl number, but the surface shear stress in the secondary flow decreases.     

A least-squares method for solving the mixed form of the groundwater flow equations

The mixed form of the areal groundwater flow equations is solved with a least-squares finite element procedure (LESFEM). Hydraulic head and x- and y-directed fluxes are state variables. Physical parameters and state variables are approximated using a bilinear basis. Grid refinements and irregular domain boundaries are implemented on rectangular meshes. Residuals are constructed at collocation points for conservation of mass and Darcy's law. Boundary condition residuals are constructed at discrete points along the boundary. The residuals are weighted, squared and summed. A set of algebraic equations is formed by taking the derivatives of the weighted sum of the squares of the residuals with respect to each unknown parameter in the approximation for the state variable and setting them to zero. Proper choice of a potential scaling parameter and residual weights is essential for the effective application of the algorithm. Test problem results demonstrate that the method is effective for both transient and steady state cases. The LESFEM algorithm generates a C°-continuous velocity field. The continuous velocity field and the rectangular mesh simplify the implementation of algorithms that require tracking. In addition, rectangular meshes simplify mesh and boundary generation.     

The interplay between boundary conditions and flow geometries in shear banding: Hysteresis, band configurations, and surface transitions

We study shear banding flows in models of wormlike micelles or polymer solutions, and explore the effects of different boundary conditions for the viscoelastic stress. These are needed because the equations of motion are inherently non-local and include “diffusive” or square-gradient terms. Using the diffusive Johnson–Segalman model and a variant of the Rolie-Poly model for entangled micelles or polymer solutions, we study the interplay between different boundary conditions and the intrinsic stress gradient imposed by the flow geometry. We consider prescribed gradient (Neumann) or value (Dirichlet) of the viscoelastic stress tensor at the boundary, as well as mixed boundary conditions in which an anchoring strength competes with the gradient contribution to the stress dynamics. We find that hysteresis during shear rate sweeps is suppressed if the boundary conditions favor the state that is induced by the sweep. For example, if the boundaries favor the high shear rate phase then hysteresis is suppressed at the low shear rate edges of the stress plateau. If the boundaries favor the low shear rate state, then the high shear rate band can lie in the center of the flow cell, leading to a three-band configuration. Sufficiently strong stress gradients due to curved flow geometries, such as that of cylindrical Couette flow, can convert this to a two-band state by forcing the high shear rate phase against the wall of higher stress, and can suppress the hysteresis loop observed during a shear rate sweep.     

COVOLUME-DUAL VARIABLE METHOD FOR THERMALLY EXPANDABLE FLOW ON UNSTRUCTURED TRIANGULAR GRIDS

Finite difference like discretizations are developed for the time dependent Navier-Stokes equations and the thermal energy equation. The flow is assumed to be thermally expandable, that is, the density varies only with temperature. A new pointwise first order upwind scheme for convection is presented which is of nonnegative type. Also presented are new approaches to reconstructing the velocity vector field from the covolume primitive variables. The resulting difference equations reproduce linear flow fields.     

A NUMERICAL STUDY OF THE FLOW OVER ELLIPSOIDAL OBJECTS INSIDE A CYLINDRICAL TUBE

A numerical scheme is developed to obtain the flow field around one, two and five ellipsoidal objects inside a cylindrical tube. The scheme uses the Galerkin finite element technique and the primitive variable(u–v–p) formulation. The two-dimensional incompressible Navier–Stokes equations are solved numerically by using the direct mixed interpolation method. A Picard iteration scheme is used for the solution of the resulting system of non-linear algebraic equations. The computer code is verified by checking with known analytical solutions for the flow past a sphere. Results for the shear stress distributions along the ellipsoids, forces and drag coefficients are obtained for different geometric ratios and Reynolds numbers. Some of the intermediate computational results on the velocity fields developed are also reported.     

Numerical study on a vertical plate with variable viscosity and thermal conductivity

Free convection over an isothermal vertical plate immersed in a fluid with variable viscosity and thermal conductivity is studied in this paper. We consider the two-dimensional, laminar and unsteady boundary layer equations. Using the appropriate variables, the basic governing equations are transformed to non-dimensional governing equations. These equations are then solved numerically using a very efficient implicit finite difference scheme known as Crank–Nicolson scheme. The fluid considered in this study is of viscous incompressible fluid of temperature dependent viscosity and thermal conductivity. The effect of varying viscosity and thermal conductivity on velocity, temperature, shear stress and heat transfer rate are discussed. The velocity and temperature profiles are compared with previously published works and are found to be in good agreement.