SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A new Navier-Stokes solver is developed by combining the efficiency of the LU-SSOR scheme and the accuracy of the flux-limited dissipation scheme. Application to laminar and turbulent flows and hypersonic flows proves the reliability of the new algorithm.
A rapidly converging relaxation technique for the parabolized Navier-Stokes equations has been devised. The scheme is applicable in both supersonic and subsonic flows, but it is discussed here in the context of supersonic flows. The upstream propagating acoustic influence in the subsonic part of the flow is introduced semi-implicitly through the streamwise momentum equation applied on the body, and through a forward-differencing on the streamwise pressure gradient term in the interior. This procedure yields a new boundary condition on the energy in the total energy equation. The pressure-velocity system in the subsonic layer is coupled, but the positive time-like marching characteristic of the governing equations is still maintained. The relaxation technique is demontrated to work for a three-dimensional flow over a cone-flare in supersonic flight.
The Bardina and Lombard (1985) three-dimensional CSCM Navier-Stokes method is presently extended to the simulation of complex hypersonic reentry vehicle external flows at angle of attack. The robust stability of the method derives from the combination of conservative implicit upwind flux difference splitting with a three-dimensional diagonally-dominant approximate factorization and relaxation scheme and characteristic-based implicit boundary approximations. The method's efficiency derives from an implicit symmetric Gauss-Seidel 'method of planes' relaxation scheme with alternating directional space marching sweeps along the flow coordinate direction.
The solution of the incompressible Navier-Stokes equations by means of a spectral method is used in the present numerical investigation of viscous axisymmetric vortex rings, in the cases of a single ring over the 0.001-1000 Re number range, and a pair of interacting rings. The propagation speeds of vortex rings of varying Re are computed and compared with the results of Saffman's (1970) theory; the present error estimate is found to be smaller than Saffman's, and to decrease with increasing Re at fixed core/ring radius ratio. The 'leapfrogging' of two vortex rings of equal sign and Re=1000 each is observed, indicating severe strain of the vortex core through the first pass and merging during the second pass.
A real gas, upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the two-dimensional/axisymmetric hypersonic flow of equilibrium air around various body shapes. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence, Tannehill and Chaussee. The upwind algorithm is based on Roe's flux-difference splitting scheme which has been modified to account for real gas effects. Simplified curve fits are used to obtain the thermodynamic and transport properties of equilibrium air. The new code has been validated by computing the hypersonic laminar flow of air over a flat plate, a wedge, a ramp, and a cone. The results of these computations are compared with the results from a conventional centrally-differenced, real gas, PNS code and the agreement is excellent, except in the vicinity of shock waves where the present code exhibits far superior shock capturing capabilities.
The paper presents procedures for the solution of the Navier-Stokes equations in the vorticity-stream function form. The difficulties involved are related to the convection term in the vorticity transport equation and to the lack of boundary conditions for voritcity at no-slip surfaces. The use of a streamline-upwind/Petrov-Galerkin finite element formulation for the solution of the vorticity transport equation. In the present scheme, the weighting functions are dependent on both spatial and temporal discretizations. A proper numerical treatment of the boundary conditions leading to an implicit treatment of the vorticity at no-slip surfaces is presented. These procedures have successfully been employed to simulate various flows of engineering interest.
A new Navier-Stokes solver is developed by combining the efficiency of the LU-SSOR scheme and the accuracy of the flux-limited dissipation scheme. Application to laminar and turbulent flows and hypersonic flows proves the reliability of the new algorithm.
A second-order parabolized Navier-Stokes algorithm based on the MacCormack (1969) explicit scheme is used to study three-dimensional chemically reacting flows with finite-rate chemistry. The method can treat the chemical source term implicitly, and it accounts for the multicomponent diffusion and convection of the chemical species. The method is demonstrated with the nonreacting case of a Mach-3 flow over a double-wedge compression corner and the case of streamwise hydrogen injection at sonic velocity in a Mach-2.44 vitiated air stream.
A real gas, upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional hypersonic flow of equilibrium air around various body shapes. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence et al. (1986). The upwind algorithm is based on Roe's (1981) flux-difference splitting scheme which has been modified to account for real gas effects using the nearly exact approach of Vinokur and Liu (1988). Simplified curve fits are employed to obtain the thermodynamic and transport properties of equilibrium air. The new code has been validated by computing the M-infinity = 25 laminar flow of air over cones at various angles of attack. The results of these computations are compared with the results from a conventional centrally-differenced, real gas PNS code and the previous axisymmetric, upwind, real gas code. The agreement is excellent in all cases.
Finite-difference approximations for steady-state compressible Navier-Stokes equations, whose two spatial dimensions are written in generalized curvilinear coordinates and strong conservation-law form, are presently solved by means of Newton's method in order to obtain a lifting-airfoil flow field under subsonic and transonnic conditions. In addition to ascertaining the computational requirements of an initial guess ensuring convergence and the degree of computational efficiency obtainable via the approximate Newton method's freezing of the Jacobian matrices, attention is given to the need for auxiliary methods assessing the temporal stability of steady-state solutions. It is demonstrated that nonunique solutions of the finite-difference equations are obtainable by Newton's method in conjunction with a continuation method.
A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.
Report discusses research on dual-potential formulation of Navier-Stokes equations. Purpose of study to extend capability of dual-potential method to computation of unsteady flows of viscous, compressible fluids.
Full 3D Navier-Stokes computational results are presented for compressible flows within nondiffusing and diffusing S-ducts. The present study provides an understanding of the performance characteristics of typical S-ducts with attached and separated flows and provides a frame of reference for future computational fluid dynamic studies of internal flows with strong secondary flows. The predicted results, which were obtained using both H- and O-grids, are compared with the experimental wall pressures, static and total pressure fields, and velocity vectors. In addition, computed boundary-layer thickness, velocity profiles in wall coordinates and skin friction values are presented. The inviscid contributions to the secondary flows are quantified. The S-duct entrance Mach number was 0.6, and the Reynolds number was 1.76 x 10 exp 6 based on the upstream duct diameter.
A detailed comparison is made between Navier-Stokes and DSMC calculations for flows near the continuum limit to assess the accuracy of the continuum equations in this regime. Meaningful comparisons require the use of similar physical models. This necessitates the inclusion of a separate rotational energy equation and use of slip boundary conditions. Inclusion of slip boundary conditions resulted in improved agreement between surface properties. Moreover, good agreement was obtained for the various temperatures in the nonequilibrium portion of the flow field that does not contain the shock region. Departures are noted in the shock region and in regions where thermal diffusion effects are important.
Axisymmetric subsonic diffuser flows were calculated with the NPARC Navier-Stokes code in order to determine the effects various code features have on the flow solutions. The code features examined in this work were turbulence models and boundary conditions. Four turbulence models available in NPARC were used: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, and the Chien kappa-epsilon and Wilcox kappa-omega two-equation models. The three boundary conditions examined were the free boundary, the mass flux boundary and the subsonic outflow with variable static pressure. In addition to boundary condition type, the geometry downstream of the diffuser was varied to see if upstream influences were present. The NPARC results are compared with experimental data and recommendations are given for using NPARC to compute similar flows.
A multigrid-mask method for solution of incompressible Navier-Stokes equations in primitive variable form has been developed. The main objective is to apply this method in conjunction with the pseudospectral element method solving flow past multiple objects. There are two key steps involved in calculating flow past multiple objects. The first step utilizes only Cartesian grid points. This homogeneous or mask method step permits flow into the interior rectangular elements contained in objects, but with the restriction that the velocity for those Cartesian elements within and on the surface of an object should be small or zero. This step easily produces an approximate flow field on Cartesian grid points covering the entire flow field. The second or heterogeneous step corrects the approximate flow field to account for the actual shape of the objects by solving the flow field based on the local coordinates surrounding each object and adapted to it. The noise occurring in data communication between the global (low frequency) coordinates and the local (high frequency) coordinates is eliminated by the multigrid method when the Schwarz Alternating Procedure (SAP) is implemented. Two dimensional flow past circular and elliptic cylinders will be presented to demonstrate the versatility of the proposed method. An interesting phenomenon is found that when the second elliptic cylinder is placed in the wake of the first elliptic cylinder a traction force results in a negative drag coefficient.
This paper presents an Unstructured Navier-Stokes Analysis of Full TCA (Technology Concept Airplane) Configuration. The topics include: 1) Motivation; 2) Milestone and approach; 3) Overview of the unstructured-grid system; 4) Results on full TCA W/B/N/D/E configuration; 5) Concluding remarks; and 6) Future directions.
A parallel version of an unstructured-grid based Navier-Stokes solver, USM3Dns, previously developed for efficient operation on a variety of parallel computers, has been enhanced to incorporate upgrades made to the serial version. The resultant parallel code has been extensively tested on a variety of problems of aerospace interest and on two sets of parallel computers to understand and document its characteristics. An innovative grid renumbering construct and use of non-blocking communication are shown to produce superlinear computing performance. Preliminary results from parallelization of a recently introduced "porous surface" boundary condition are also presented.