The numerical solution of the transient two‐phase flow in rigid pipelinesHadj‐Taieb, Ezzeddine; Lili, Taieb
doi: 10.1002/(SICI)1097-0363(19990315)29:5<501::AID-FLD716>3.0.CO;2-8pmid: N/A
Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.
Convergence acceleration of segregated algorithms using dynamic tuning additive correction multigrid strategyZdravistch, Franz; Fletcher, Clive A.J.; Behnia, Masud
doi: 10.1002/(SICI)1097-0363(19990315)29:5<515::AID-FLD798>3.0.CO;2-Ipmid: N/A
A convergence acceleration method based on an additive correction multigrid–SIMPLEC (ACM‐S) algorithm with dynamic tuning of the relaxation factors is presented. In the ACM‐S method, the coarse grid velocity correction components obtained from the mass conservation (velocity potential) correction equation are included into the fine grid momentum equations before the coarse grid momentum correction equations are formed using the additive correction methodology. Therefore, the coupling between the momentum and mass conservation equations is obtained on the coarse grid, while maintaining the segregated structure of the single grid algorithm. This allows the use of the same solver (smoother) on the coarse grid. For turbulent flows with heat transfer, additional scalar equations are solved outside of the momentum–mass conservation equations loop. The convergence of the additional scalar equations is accelerated using a dynamic tuning of the relaxation factors. Both a relative error (RE) scheme and a local Reynolds/Peclet (ER/P) relaxation scheme methods are used. These methodologies are tested for laminar isothermal flows and turbulent flows with heat transfer over geometrically complex two‐ and three‐dimensional configurations. Savings up to 57% in CPU time are obtained for complex geometric domains representative of practical engineering problems. Copyright © 1999 John Wiley & Sons, Ltd.
Numerical simulation of a flow around an impulsively started radially deforming circular cylinderHanchi, Samir; Askovic, Radomir; Ta Phuoc, Loc
doi: 10.1002/(SICI)1097-0363(19990315)29:5<555::AID-FLD800>3.0.CO;2-7pmid: N/A
The development of a two‐dimensional viscous incompressible flow generated by a deformable circular cylinder impulsively started into rectilinear motion is studied numerically for the Reynolds numbers equal to 550 and 3000. The vorticity transport equation is solved by a second‐order finite difference method in both directions of the domains. The Poisson equation for the streamfunction is solved by a Fourier–Galerkin method in the direction of the flow that is assumed to remain symmetrical and a second‐order finite difference for the radial direction. The advance in time is achieved by a second‐order Adams–Bashforth scheme. The computed results are compared qualitatively with experimental and numerical results done before in the particular non‐deformable case. The comparison is found to be satisfactory. The influence of the deformation of the cylinder on the flow structure and the drag coefficient is then analyzed. Copyright © 1999 John Wiley & Sons, Ltd.
Computation of unsteady transonic flow over a fighter wing using a zonal Navier–Stokes/full‐potential methodde Faria Mello, Olympio Achilles; Sankar, Lakshmi N.
doi: 10.1002/(SICI)1097-0363(19990315)29:5<575::AID-FLD801>3.0.CO;2-2pmid: N/A
An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the inner zone, the Navier–Stokes equations are solved using a diagonal form of an alternating‐direction implicit (ADI) approximate factorisation procedure. In the outer zone, the unsteady full‐potential equation (FPE) is solved. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic‐based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are about 60% of the full Navier–Stokes CPU times for unsteady flows in non‐vector processing machines. Applications of the method are presented for a F‐5 wing in steady and unsteady transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to the full Navier–Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier–Stokes equations can be retained with significant savings in computational time. Copyright © 1999 John Wiley & Sons, Ltd.
A discontinuous Galerkin method for the Navier–Stokes equationsLomtev, Igor; Karniadakis, George Em
doi: 10.1002/(SICI)1097-0363(19990315)29:5<587::AID-FLD805>3.0.CO;2-Kpmid: N/A
The foundations of a new discontinuous Galerkin method for simulating compressible viscous flows with shocks on standard unstructured grids are presented in this paper. The new method is based on a discontinuous Galerkin formulation both for the advective and the diffusive contributions. High‐order accuracy is achieved by using a recently developed hierarchical spectral basis. This basis is formed by combining Jacobi polynomials of high‐order weights written in a new co‐ordinate system. It retains a tensor‐product property, and provides accurate numerical quadrature. The formulation is conservative, and monotonicity is enforced by appropriately lowering the basis order and performing h‐refinement around discontinuities. Convergence results are shown for analytical two‐ and three‐dimensional solutions of diffusion and Navier–Stokes equations that demonstrate exponential convergence of the new method, even for highly distorted elements. Flow simulations for subsonic, transonic and supersonic flows are also presented that demonstrate discretization flexibility using hp‐type refinement. Unlike other high‐order methods, the new method uses standard finite volume grids consisting of arbitrary triangulizations and tetrahedrizations. Copyright © 1999 John Wiley & Sons, Ltd.
Discrete non‐local absorbing boundary condition for exterior problems governed by Helmholtz equationBonet, Ruperto P.; Nigro, Norberto; Storti, Mario A.; Idelsohn, Sergio R.
doi: 10.1002/(SICI)1097-0363(19990315)29:5<605::AID-FLD808>3.0.CO;2-Mpmid: N/A
The finite element method is employed to approximate the solutions of the Helmholtz equation for water wave radiation and scattering in an unbounded domain. A discrete, non‐local and non‐reflecting boundary condition is specified at an artificial external boundary by the DNL method, yielding an equivalent problem that is solved in a bounded domain. This procedure formulates a boundary value problem in a bounded region by imposing a relation in the discrete medium between the nodal values at the two last layers. For plane geometry, this relation can be found by straightforward eigenvalue decomposition. For circular geometry, the plane condition is applied at the external layer and this condition is condensed through a structured annular region, resulting in a condition at an inner radius. Exterior problems with a bounded internal physical obstacle are considered. It is well‐known that these kind of problems are well‐posed, and have a unique solution. Numerical studies based on standard Galerkin methodology examine the dependence of the DNL condition with respect to the circular annular region width. The DNL condition is compared with local boundary conditions of several orders. Numerical examples confirm the important improvement in accuracy obtained by the DNL method over standard conditions. Copyright © 1999 John Wiley & Sons, Ltd.