A Taylor-Galerkin Method for Modeling Unsteady Fluid Flow Under its Own Weight

The free fluid-surface of incompressible creeping flows is analyzed using a finite element method. A pseudo-concentration (PC) function is introduced to determine the position of the free surface. The Taylor-Galerkin finite element method (TGFEM) is applied to solve the equation of the PC function. Nine-node quadratic interpolation is used for both PC and velocity. The unsteady flows of fluids moving of their own weight are analyzed using the proposed method.     

Stability of plane-parallel vibrational flow in a two-layer system

The stability of the interface separating two immiscible incompressible fluids of different densities and viscosities is considered in the case of fluids filling a cavity which performs horizontal harmonic oscillations. There exists a simple basic state which corresponds to the unperturbed interface and plane-parallel unsteady counter flows; the properties of this state are examined. A linear stability problem for the interface is formulated and solved for both (a) inviscid and (b) viscous fluids. A transformation is found which reduces the linear stability problem under the inviscid approximation to the Mathieu equation. The parametric resonant regions of instability associated with the intensification of capillary-gravity waves at the interface are examined and the results are compared to those found in the viscous case in a fully numerical investigation.     

Functional and generalized separable solutions to unsteady Navier–Stokes equations

The paper studies unsteady Navier–Stokes equations with two space variables. It shows that the non-linear fourth-order equation for the stream function with three independent variables admits functional separable solutions described by a system of three partial differential equations with two independent variables. The system is found to have a number of exact solutions, which generate new classes of exact solutions to the Navier–Stokes equations. All these solutions involve two or more arbitrary functions of a single argument as well as a few free parameters. Many of the solutions are expressed in terms of elementary functions, provided that the arbitrary functions are also elementary; such solutions, having relatively simple form and presenting significant arbitrariness, can be especially useful for solving certain model problems and testing numerical and approximate analytical hydrodynamic methods. The paper uses the obtained results to describe some model unsteady flows of viscous incompressible fluids, including flows through a strip with permeable walls, flows through a strip with extrusion at the boundaries, flows onto a shrinking plane, and others. Some blow-up modes, which correspond to singular solutions, are discussed.     

Analytical solutions for equations of unsteady flow of non-newtonian fluids in tube

This paper presents analytical solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse the behaviour of axial velocity and shear stress for unsteady flow of non-Newtonian visco-elastic fluids in tube, and to provide a theoretical base for the projection of pipe-line engineering.     

Transient flows of a second grade fluid

Exact analytical solutions for a class of unsteady unidirectional flows of an incompressible second-order fluid are constructed. The flows are generated impulsively from rest by motion of a plate or two plates or by sudden application of a pressure gradient. Expressions for velocity, flux and skin friction are obtained for both large and small times. It is found that large and small times solutions are dependent on the coefficient of viscoelasticity. The solutions corresponding to Newtonian fluids can be easily obtained from those for fluids of second order by letting the viscoelastic parameter to be zero.     

Four-field Galerkin/least-squares formulation for viscoelastic fluids

A new Galerkin/Least-Squares (GLS) stabilized finite element method is presented for computing viscoelastic flows of complex fluids described by the conformation tensor; it extends the well-established GLS method for computing flows of incompressible Newtonian fluids. GLS methods are attractive for large-scale computations because they yield linear systems that can be solved easily with iterative solvers (e.g., the Generalized Minimum Residual method) and because they allow simple combinations of interpolation functions that can be conveniently and efficiently implemented on modern distributed-memory cache-based clusters.Like other state-of-the-art methods for computing viscoelastic flows (e.g., DEVSS-TG/SUPG), the new GLS method introduces a separate variable to represent the velocity gradient; with the aid of this variable, the conservation equations of mass, momentum, conformation, and the definition of velocity gradient are converted into a set of first-order partial differential equations in four unknown fields—pressure, velocity, conformation, and velocity gradient. The unknown fields are represented by low-order (continuous piecewise linear or bilinear) finite element basis functions.The method is applied to the Oldroyd-B constitutive equation and is tested in two benchmark problems—flow in a planar channel and flow past a cylinder in a channel. Results show that (1) the mesh-convergence rate of GLS is comparable to the DEVSS-TG/SUPG method; (2) the LS stabilization permits using equal-order basis functions for all fields; (3) GLS handles effectively the advective terms in the evolution equation of the conformation tensor; and (4) GLS yields accurate results at lower computational costs than DEVSS-type methods.     

Applications of URANS on predicting unsteady turbulent separated flows

Accurate prediction of unsteady separated turbulent flows remains one of the toughest tasks and a practi cal challenge for turbulence modeling. In this paper, a 2D flow past a circular cylinder at Reynolds number 3,900 is numerically investigated by using the technique of unsteady RANS （URANS）. Some typical linear and nonlinear eddy viscosity turbulence models （LEVM and NLEVM） and a quadratic explicit algebraic stress model （EASM） are evaluated. Numerical results have shown that a high-performance cubic NLEVM, such as CLS, are superior to the others in simulating turbulent separated flows with unsteady vortex shedding.     

A matrix‐free preconditioned Newton/GMRES method for unsteady Navier–Stokes solutions

The unsteady compressible Reynolds‐averaged Navier–Stokes equations are discretized using the Osher approximate Riemann solver with fully implicit time stepping. The resulting non‐linear system at each time step is solved iteratively using a Newton/GMRES method. In the solution process, the Jacobian matrix–vector products are replaced by directional derivatives so that the evaluation and storage of the Jacobian matrix is removed from the procedure. An effective matrix‐free preconditioner is proposed to fully avoid matrix storage. Convergence rates, computational costs and computer memory requirements of the present method are compared with those of a matrix Newton/GMRES method, a four stage Runge–Kutta explicit method, and an approximate factorization sub‐iteration method. Effects of convergence tolerances for the GMRES linear solver on the convergence and the efficiency of the Newton iteration for the non‐linear system at each time step are analysed for both matrix‐free and matrix methods. Differences in the performance of the matrix‐free method for laminar and turbulent flows are highlighted and analysed. Unsteady turbulent Navier–Stokes solutions of pitching and combined translation–pitching aerofoil oscillations are presented for unsteady shock‐induced separation problems associated with the rotor blade flows of forward flying helicopters. Copyright © 2000 John Wiley & Sons, Ltd.     

Reduced order state‐space models from the pulse responses of a linearized CFD scheme

This paper describes a method for obtaining a time continuous reduced order model (ROM) from a system of time continuous linear differential equations. These equations are first put into a time discrete form using a finite difference approximation. The unit sample responses of the discrete system are calculated for each system input and these provide the Markov parameters of the system. An eigenvalue realization algorithm (ERA) is used to construct a discrete ROM. This ROM is then used to obtain a continuous ROM of the original continuous system. The focus of this paper is on the application of this method to the calculation of unsteady flows using the linearized Euler equations on moving meshes for aerofoils undergoing heave or linearized pitch motions. Applying a standard cell‐centre spatial discretization and taking account of mesh movement a continuous system of differential equations is obtained which are continuous in time. These are put into discrete time form using an implicit finite difference approximation. Results are presented demonstrating the efficiency of the system reduction method for this system. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulations of non-stationary flows of non-Newtonian fluids between concentric and eccentric cylinders by stream-tube method and domain decomposition

This study involves a theoretical formulation of the stream-tube method in non-stationary flows. Initially, this approach allowed flow computations by determining an unknown transformation between the physical domain and a mapped domain where the streamlines are rectilinear and parallel. To take into account vortex zones, we define local transformations of subregions of the physical domain that are mapped into rectangular domains where the transformed streamlines are still parallel and straight. The local functions must be determined numerically from the governing equations and boundary conditions put together with compatibility equations. The method enables to compute streamlines and flow data at every time, using distinguishing properties, as verification of mass conservation and definition of rectangular meshes allowing to adopt finite-difference schemes. The numerical simulations concern different non-Newtonian fluids under various geometrical and kinematic specifications related to flows between concentric and eccentric cylinders, leading to comparisons with literature data. The results also highlight the influence of the rheological properties on the flow characteristics in unsteady conditions.     

A general formula for the drag on a sphere placed in a creeping unsteady micropolar fluid flow

In the present work, we investigate the creeping unsteady motion of an infinite micropolar fluid flow past a fixed sphere. The technique of Laplace transform is used. The drag formula is obtained in the physical domain analytically by using the complex inversion formula of the Laplace transform. The well known formula of Basset for the drag on a sphere placed in an unsteady viscous fluid flow and that of Ramkissoon and Majumdar for steady motion in the case of micropolar fluids are recovered as special cases. The obtained formula is employed to calculate the drag force for some micropolar fluid flows. Numerical results are obtained and represented graphically.     

Steady-state solutions for some simple flows of generalized Burgers fluids

The steady-state solutions for three types of unsteady oscillating flows of generalized Burgers fluids are determined by means of the Fourier sine transforms. These solutions are also presented in equivalent forms in terms of elementary functions exp, sine, cosine, hyperbolic sine and hyperbolic cosine. The similar solutions for Burgers, Oldroyd-B, Maxwell, Second grade and Navier-Stokes fluids can be also obtained as limiting cases of our solutions.     

Analysis on creeping channel flows of compressible fluids subject to wall slip

Creeping channel flows of compressible fluids subject to wall slip are widely encountered in industries. This paper analyzes such flows driven by pressure in planar as well as circular channels. The analysis elucidates unsteady flows of Newtonian fluids subject to the Navier slip condition, followed by steady flows of viscoplastic fluids, in particular, Herschel–Bulkley fluids and their simplifications including power law and Newtonian fluids, that slip at wall with a constant coefficient or a coefficient inversely proportional to pressure. Under the lubrication assumption, analytical solutions are derived, validated, and discussed over a wide range of parameters. Analysis based on the derived solutions indicates that unsteadiness alters cross-section velocity profiles. It is demonstrated that compressibility of the fluids gives rise to a concave pressure distribution in the longitudinal direction, whereas wall slip with a slip coefficient that is inversely proportional to pressure leads to a convex pressure distribution. Energy dissipation resulting from slippage can be a significant portion in the total dissipation of such a flow. A distinctive feature of the flow is that, in case of the pressure-dependent slip coefficient, the slip velocity increases rapidly in the flow direction and the flow can evolve into a pure plug flow at the exit.     

Unsteady flow of an elastico-viscous fluid in tubes of uniform cross-section

In this paper, the uniqueness of solution for internal bounded unsteady flows of a shortmemory fluid is first established. Closed-form solutions are then obtained for the equations characterizing flows of such fluids in circular and rectangular tubes of uniform cross-section under an arbitrary pressure gradient. Special cases including the oscillatory flow between two parallel plates are discussed.     

An incompressible SPH method for simulation of unsteady viscoelastic free-surface flows

In this paper, an incompressible smoothed particle hydrodynamics (SPH) method is presented to solve unsteady free-surface flows. Both Newtonian and viscoelastic fluids are considered. In the case of viscoelastic fluids, both the Maxwell and Oldroyd-B models are investigated. The proposed SPH method uses a Poisson pressure equation to satisfy the incompressibility constraints. The solution algorithm is an explicit predictor-corrector scheme and employs an adaptive smoothing length based on density variations. To alleviate the numerical difficulties encountered when fluid is highly stretched, an artificial stress term is incorporated into the momentum equation which reduces the risk of unrealistic fractures in the material. Two challenging test cases, the impacting drop and the jet buckling problems, are solved to demonstrate the capability of the proposed scheme in handling viscoelastic flows with complex free surfaces. The jet buckling test case was solved for a wide range of Weissenberg numbers. It was shown that in all cases the method is stable and fairly accurate and agrees well with the available data.     

UNSTEADY CONFINED VISCOUS FLOWS WITH OSCILLATING WALLS AND MULTIPLE SEPARATION REGIONS OVER A DOWNSTREAM-FACING STEP

This paper presents computational solutions for unsteady viscous flows in channels with a downstream-facing step, followed by an oscillating floor. These solutions of the unsteady Navier–Stokes equations are obtained with a time-integration method using artificial compressibility in a fixed computational domain, which is obtained via a time-dependent coordinate transformation from the fluid domain with moving boundaries. The computational method is first validated for steady viscous flows past a downstream-facing step by comparison with previous numerical solutions and experimental results. This method is then used to obtain solutions for unsteady viscous flows with multiple separation regions over a downstream-facing step with oscillating walls, for which there are no previously known solutions. Thus, the present results may be used as benchmark solutions for the unsteady viscous flows with multiple separation regions between fixed and oscillating walls.     

Reduced-order modeling of three-dimensional unsteady partial cavity flows

An efficient reduced-order modeling to analyze three-dimensional unsteady partial cavity flows is proposed. The proposed approach is based on the boundary element method along with the potential flow assumption. To this end, a novel non-iterative method based on the flow eigenmodes of three-dimensional partial cavity flows is applied. Eigenanalysis and reduced-order modeling for unsteady flows over a three-dimensional hydrofoil with various sections are performed. The results obtained from the present analysis are compared with those reported in the literature to verify the strength of the proposed approach. In order to examine the performance of the introduced algorithm for unsteady cavitating flows, various simulations for several reduced frequencies, hydrofoil geometries and different cavitation numbers are also investigated. Comparison between the obtained results using the novel and conventional methods indicates that the present algorithm works very well with sufficient accuracy. Moreover, it is shown that the proposed method is computationally more efficient than the conventional ones for unsteady sheet cavitation analysis on three-dimensional hydrofoils.     

PEGASE: A NAVIER–STOKES SOLVER FOR DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS

A hybrid conservative finite difference/finite element scheme is proposed for the solution of the unsteady incompressible Navier–Stokes equations. Using velocity–pressure variables on a non-staggeredgrid system, the solution is obtained with a projection method basedon the resolution of a pressure Poisson equation. The new proposed scheme is derived from the finite element spatial discretization using the Galerkin method with piecewise bilinear polynomial basis functions defined on quadrilateral elements. It is applied to the pressure gradient term and to the non-linear convection term as in the so-called group finite element method. It ensures strong coupling between spatial directions, inhibiting the development of oscillations during long-term computations, as demonstrated by the validation studies. Two- and three-dimensional unsteady separated flows with open boundaries have been simulated with the proposed method using Cartesian uniform mesh grids. Several examples of calculations on the backward-facing step configuration are reported and the results obtained are compared with those given by other methods. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids 24: 833–861, 1997.     

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids.