Engineering topics

Leonard, A.

Publications and source records attributed to Leonard, A..

At least 19 records

NASA NTRS1989

Improved vortex methods for three-dimensional flows

Robust numerical methods are developed for three-dimensional incompressible vortical flows, using Lagrangian vortex elements. A successful scheme must be able to handle regions of intense vortex stretching and vortex reconnection with reasonable accuracy (without diverging). Here, consideration is given to vortex particles, also commonly called vortons or vortex sticks. The following issues are discussed: (1) use of delta-function elements and weak solutions of the vorticity equation; (2) use of smoothed elements and the choice of the smoothing function; (3) representation of viscous effects and the redistribution of element strength; and (4) conservation laws (are they satisfied?). The various proposed schemes have been tested on flows involving a strong interaction between two vortex rings.

Winckelmans, G.↗

NASA NTRS1988

Weak solutions of the three-dimensional vorticity equation with vortex singularities

The extension of the concept of vortex singularities, developed by Saffman and Meiron (1986) for the case of two-dimensional point vortices in an incompressible vortical flow, to the three-dimensional case of vortex sticks (vortons) is investigated analytically. The derivation of the governing equations is explained, and it is demonstrated that the formulation obtained conserves total vorticity and is a weak solution of the vorticity equation, making it an appropriate means for representing three-dimensional vortical flows with limited numbers of vortex singularities.

Winckelmans, G.↗

NASA NTRS1988

Robust vortex methods for three-dimensional incompressible flows

Vortex methods for the numerical simulation of incompressible three-dimensional flows at high Reynolds number are discussed. Particular emphasis is placed on schemes that prevent the excessive generation of mesh elements, include viscous effects, and are relatively inexpensive even as the number of elements becomes large. Several simulations are presented involving the interaction of vortex rings.

Chua, K.↗

NASA NTRS1986

Vortex simulation of forced mixing layers

Two-dimensional, spatially growing, turbulent mixing layers are simulated numerically by a vortex method and the results are compared with those determined experimentally. The effects of artificial forcing on flow development are also studied. Many of the flow features which have been observed experimentally are reproduced, and good quantitative agreements between experiments and computations are obtained.

Inoue, O.↗

NASA NTRS1986

Internal structure of a vortex breakdown

An axisymmetric vortex breakdown was well simulated by the vortex filament method. The agreement with the experiment was qualitatively good. In particular, the structure in the interior of the vortex breakdown was ensured to a great degree by the present simulation. The second breakdown, or spiral type, which occurs downstream of the first axisymmetric breakdown, was simulated more similarly to the experiment than before. It shows a kink of the vortex filaments and strong three-dimensionality. Furthermore, a relatively low velocity region was observed near the second breakdown. It was also found that it takes some time for this physical phenomenon to attain its final stage. The comparison with the experiment is getting better as time goes on. In this paper, emphasis is placed on the comparison of the simulated results with the experiment. The present results help to make clear the mechanism of a vortex breakdown.

Nakamura, Y.↗

NASA NTRS1986

Evolution of a curved vortex filament into a vortex ring

The deformation of a hairpin-shaped vortex filament under self-induction and in the presence of shear is studied numerically using the Biot-Savart law. It is shown that the tip region of an elongated hairpin vortex evolves into a vortex ring and that the presence of mean shear impedes the process. Evolution of a finite-thickness vortex sheet under self-induction is also investigated using the Navier-Stokes equations. The layer evolves into a hairpin vortex which in turn produces a vortex ring of high Reynolds stress content. These results indicate a mechanism for the generation of ring vortices in turbulent shear flows, and a link between the experimental and numerical observation of hairpin vortices and the observation of ring vortices in the outer regions of turbulent boundary layers.

Moin, P.↗

NASA NTRS1985

Evolution of a curved vortex filament into a vortex ring

Moin, P.↗

NASA NTRS1985

Vortex breakdown simulation

A vortex breakdown was simulated by the vortex filament method, and detailed figures are presented based on the results. Deformations of the vortex filaments showed clear and large swelling at a particular axial station which implied the presence of a recirculation bubble at that station. The tendency for two breakdowns to occur experimentally was confirmed by the simulation, and the jet flow inside the bubble was well simulated. The particle paths spiralled with expansion, and the streamlines took spiral forms at the breakdown with expansion.

Nakamura, Y.↗

NASA NTRS1985

Computing three-dimensional incompressible flows with vortex elements

The techniques, capabilities and applicability of numerical models of three-dimensional, unsteady vortical flows with high Re are assessed. Vorticity is calculated only in appropriate regions and the velocity field is derived from the boundary conditions. Vorticity is assumed to take the shape of tubes with uniform core structures in the case of turbulence. The efforts being made to simplify equations for dense collections of vortex filaments in order to make them tractable to computer simulations are described. The effectiveness of vorticity arrow representations for accurately describing vorticity fields near surfaces is discussed, along with Lagrangian vortex elements, which may be of use in modelling the rotational part of flows around bluff bodies, nonuniform density flows and chemically reacting flows.

Leonard, A.↗

NASA NTRS1985

Direct numerical simulation of equilibrium turbulent boundary layers

This paper describes the simulation of turbulent boundary layers by direct numerical solution of the three-dimensional, time-dependent Navier-Stokes equations, using a spectral method. The flow is incompressible, with Re sub delta = 1000. The equations are written in the self-similar coordinate system and periodic streamwise and spanwise boundary conditions are imposed. A family of nine 'equilibrium' boundary layers, from the strongly accelerated 'sink' flow to Stratford's separating boundary layer is treated. Good general agreement with experiments is observed. The effects of pressure gradients on the structures and statistics, both in the wall and wake regions are discussed.

Spalart, P. R.↗

NASA NTRS1984

Simulating a Three-Dimensional Flow in Pipes

Numerical method simulates in unsteady, incompressible flow in axisymmetric pipe in three dimensions. Method relies on velocity-vector expansion technique in which each vector function in expansion set is divergence-free and satisfies boundary conditions for viscous flow.

Leonard, A.↗

NASA NTRS1983

A spectral numerical method for the Navier-Stokes equations with applications to Taylor-Couette flow

A new spectral method for solving the incompressible Navier-Stokes equations in a plane channel and between concentric cylinders is presented. The method uses spectral expansions which inherently satisfy the boundary conditions and the continuity equation and yield banded matrices which are efficiently solved at each time step. In addition, the number of dependent variables is reduced, resulting in a reduction in computer memory requirements. Several test problems have been computed for the channel flow and for flow between concentric cylinders, including Taylor-Couette flow with axisymmetric Taylor vortices and wavy vortices. In all cases, agreement with available experimental and theoretical results is very good.

Moser, R. D.↗

NASA NTRS1983

Numerical simulation of vortex breakdown by the vortex-filament method

The vortex filament method was applied to the simulation of vortex breakdown. The principal vortex region was represented by multiple filaments, and an axial velocity component was induced by a spiral winding of the filaments. First, an accuracy check was performed for a cylindrical swirling flow with simple analytical expressions for the axial and theta velocities. The result suggests that the flow field is simulated to any accuracy by increasing the number of filaments. Second, an axisymmetric type vortex breakdown was simulated, with experimental data serving as upstream conditions. The calculated axial and theta velocity contours show the breakdown of the vortex, including a rapid change in the vortex core, followed axially by a recovery zone and then a second breakdown. When three dimensional initial data are used the second breakdown appears to be of the spiral type in correspondence with experimental observations. The present method is easily used to simulate other types of vortex breakdown or other vortex flows with axial velocity.

Nakamura, Y.↗

NASA NTRS1983

Numerical simulation of vortex breakdown by the vortex-filament method

The vortex-filament method was applied to the simulation of vortex breakdown. The principal vortex region was represented by multiple filaments, and an axial velocity component was induced by a spiral winding of the filaments. An accuracy check was performed for a cylindrical swirling flow with simple analytical expressions for the axial and theta velocities. The result suggests that the flow field can be simulated to any accuracy by increasing the number of filaments. An axisymmetric-type vortex breakdown was simulated, with experimental data serving as upstream conditions. The calculated axial- and theta-velocity contours show the breakdown of the vortex, including a rapid change in the vortex core, followed axially by a recovery zone and then a second breakdown. When three dimensional initial data are used the second breakdown appears to be of the spiral type in correspondence with experimental observations. The present method can easily be used to simulate other types of vortex breakdown or other vortex flows with axial velocity.

Nakamura, Y.↗

NASA NTRS1983

Numerical Simulation of Turbulent Fluid Flows

Numerical simulation of turbulent flows is discussed. Computational requirements for the direct simulaton of turbulence, simulation of arbitrary homogeneous flows, an expansion technique for wall bounded flows with application to pipe flow, and possibilities of flow representations or modeling techniques that allow the simulation of high Reynolds number flows with a relatively small number of dependent variables are included.

Leonard, A.↗

NASA NTRS1983

Numerical simulation of separated flows

A new numerical method, based on the Vortex Method, for the simulation of two-dimensional separated flows, was developed and tested on a wide range of gases. The fluid is incompressible and the Reynolds number is high. A rigorous analytical basis for the representation of the Navier-Stokes equation in terms of the vorticity is used. An equation for the control of circulation around each body is included. An inviscid outer flow (computed by the Vortex Method) was coupled with a viscous boundary layer flow (computed by an Eulerian method). This version of the Vortex Method treats bodies of arbitrary shape, and accurately computes the pressure and shear stress at the solid boundary. These two quantities reflect the structure of the boundary layer. Several versions of the method are presented and applied to various problems, most of which have massive separation. Comparison of its results with other results, generally experimental, demonstrates the reliability and the general accuracy of the new method, with little dependence on empirical parameters. Many of the complex features of the flow past a circular cylinder, over a wide range of Reynolds numbers, are correctly reproduced.

Spalart, P. R.↗

NASA NTRS1982

A new numerical method for the simulation of three dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction.

Leonard, A.↗

NASA NTRS1982

Vortex simulation of an inviscid shear layer

The accuracy of the vortex-blob method was tested by simulating a free-shear-layer instability, Kirchhoff's elliptical vortex, and a circular vortex. The main numerical parameters in the vortex-blob method are the density of the vortices, and the distribution of vorticity within each vortex core. The growth rate of a periodic unstable mode of the shear layer was calculated numerically and compared with the exact result. The error is only a few percent for about 10 rows of vortex blobs. The error is reduced by decreasing the spacing between vortices and, correspondingly, the core size. In the simulation of the motion of the elliptical vortex, the rotation of the boundary, without change of shape, and the circular particle paths of the vortical fluid were well simulated. For the circular vortex, optimum sets of parameters were obtained by comparing them with the exact velocity. The results are consistent with convergence theories of the vortex-blob method. In particular, second-order convergence is observed with a Gaussian core from velocity calculation.

Nakamura, Y.↗