Search NASA⌕ Search

SEARCH · Search NASA

Results for “Navier Stokes solver”

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.

At least 487 records · Page 27

Experimental and analytical study of close-coupled ventral nozzles for ASTOVL aircraft

Flow in a generic ventral nozzle system was studied experimentally and analytically with a block version of the PARC3D computational fluid dynamics program (a full Navier-Stokes equation solver) in order to evaluate the program's ability to predict system performance and internal flow patterns. For the experimental work a one-third-size model tailpipe with a single large rectangular ventral nozzle mounted normal to the tailpipe axis was tested with unheated air at steady-state pressure ratios up to 4.0. The end of the tailpipe was closed to simulate a blocked exhaust nozzle. Measurements showed about 5 1/2 percent flow-turning loss, reasonable nozzle performance coefficients, and a significant aftward axial component of thrust due to flow turning loss, reasonable nozzle performance coefficients, and a significant aftward axial component of thrust due to flow turning more than 90 deg. Flow behavior into and through the ventral duct is discussed and illustrated with paint streak flow visualization photographs. For the analytical work the same ventral system configuration was modeled with two computational grids to evaluate the effect of grid density. Both grids gave good results. The finer-grid solution produced more detailed flow patterns and predicted performance parameters, such as thrust and discharge coefficient, within 1 percent of the measured values. PARC3D flow visualization images are shown for comparison with the paint streak photographs. Modeling and computational issues encountered in the analytical work are discussed.

Mcardle, Jack G.↗

CFD analysis for high speed inlets

The increased national interest in high speed flight has increased research for high speed propulsion components. The highly 3-D flows present in supersonic/hypersonic inlets are currently being studied at NASA-Lewis both experimentally and computationally using a family of steady Parabolized Navier-Stokes (PNS) and Navier-Stokes (NS) solvers and unsteady NS solvers. Some of the results of these efforts are presented with an emphasis on the comparison of the computational and experimental results. The flow in high speed inlets typically involves the interaction of compression shock waves and boundary layers on the internal surfaces. The fundamentals of these interactions have been studied experimentally for many years, while more recently, computations have been used to study these complex 3-D flow fields. Attempts to control the flow through boundary layer bleed are being investigated computationally prior to wind tunnel experiments. The ultimate goal is the higher performing inlets required for high speed flight.

Benson, Tom↗

Numerical investigation of an internal layer in turbulent flow over a curved hill

A numerical investigation of incompressible and compressible turbulent flows over strongly curved surfaces is presented. The turbulent flow equations are solved by a pressure based Navier-Stokes equations solver. In the method, the conservation of mass equation is replaced by a pressure correction equation applicable for both compressible and incompressible flows. The turbulence is described by a multiple time scale turbulence model supplemented with a near-wall turbulence model. The numerical results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. The development of the internal layer is attributed to the enormous mean flow strain rate caused by the streamline curvature. The external boundary layer flow responds rather slowly to the streamline curvature. Thus, the turbulence field of the forward corner of the curved hill is characterized by two turbulence fields interacting with each other. The turbulence intensity of the internal layer is much stronger than that of the external boundary layer, so that the development of a new boundary layer in the downstream region of the curved hill depends mostly on the internal layer. These numerical results are in good agreement with the measured data, and show that the turbulence model can resolve the turbulence field subjected to the strong streamline curvature.

Kim, S.-W.↗

Flow studies in close-coupled ventral nozzles for STOVL aircraft

Flow in a generic ventral nozzle system was studied experimentally and analytically with the PARC3D computational fluid dynamics program (a full Navier-Stokes equations solver) in order to evaluate the program's ability to predict system performance and internal flow patterns. A generic model of a tailpipe with a rectangular ventral nozzle, about one-third of full size, was tested with unheated air at steady-state pressure ratios up to 4.0. Measurements showed about 5.5 percent flow-turning loss and reasonable nozzle performance coefficients. The flow turned more than the designed 90 deg, causing an aftward axial component in the total thrust. Flow behavior into and through the ventral duct is discussed and illustrated with paint streak flow visualization photographs. PARC3D graphic images are shown for comparison with the experiment photographs. The program successfully predicted internal flow patterns; it also computed thrust and discharge coefficients within 1 percent of measured values.

Mcardle, Jack G.↗

CFD modeling of turbulent flows around the SSME main injector assembly using porosity formulation

Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle Main Engine (SSME) and liquid oxygen (LOX) flow distributions through the LOX posts have great effect on the combustion phenomenon inside the main combustion chamber. An advanced computational fluid dynamics (CFD) analysis will help to provide more accurate and efficient characterization of this type of flow field. In order to design a CFD model to be an effective engineering analysis tool with good computational turn-around time and still maintain good accuracy in describing the flow features, the concept of porosity is employed to describe the effects of blockage and drag force due to the presence of the LOX posts in the turbulent flow field around the main injector assembly of the SSME. A validated non-isotropic porosity model is developed and incorporated into an existing Navier-Stokes flow solver (FDNS). Volume and surface porosity parameters, which are based on the configurations of local LOX post clustering, are introduced into the governing equations.

Cheng, Gary C.↗

CFD analysis of pump consortium impeller

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design impact in a productive manner. The main goal of this study is to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A Navier-Stokes flow solver, FDNS, embedded with the extended k-epsilon turbulence model and with appropriate moving interface boundary conditions, is developed to analyze turbulent flows in the turbomachinery devices. The FDNS code was benchmarked with its numerical predictions of the pump consortium inducer, and provides satisfactory results. In the present study, a CFD analysis of the pump consortium impeller will be conducted with the application of the FDNS code. The pump consortium impeller, with partial blades, is the new design concept of the advanced rocket engine.

Cheng, Gary C.↗

Time accurate Navier-Stokes computations of solid rocket motor internal field joint cavity flows

The low-speed flow of combustion gases in the Space Shuttle Solid Rocket Motor (SRM) following ignition are investigated with a time-accurate CFD technique. The flow of the combustion gases is examined in the forward joint section of the SRM at 80 s after ignition with attention given to the amplification of existing 'organ-pipe' modes of 15, 30, and 45 Hz resulting from frequency lock-in with vortex-shedding modes. The CFD code employed is a multizone compressible Navier-Stokes flow solver using a finite-volume, upwind-biased, total-variation-diminishing scheme. The vortex-generated oscillations generated by the inhibitor and the cavity at the forward field joint are described in terms of amplitude and frequency and related to the flow dynamics in the SRM. Oscillations are shown to occur spontaneously when the trailing edge of the cavity at the joint is sharp and can be amplified if the unburned portion of the inhibitor protrudes into the flow.

Farr, Rebecca A.↗

Experimental and analytical study of close-coupled ventral nozzles for ASTOVL aircraft

Flow in a generic ventral nozzle system was studied experimentally and analytically with a block version of the PARC3D computational fluid dynamics program (a full Navier-Stokes equation solver) in order to evaluate the program's ability to predict system performance and internal flow patterns. Measurements showed about 5.5-percent flow-turning loss, reasonable nozzle performance coefficients, and a significant aftward axial component of thrust due to flow turning more than 90 deg. Flow behavior into and through the ventral duct is discussed and illustrated with paint streak flow visualization photographs. PARC3D flow visualization images are shown for comparison with the paint streak photographs. Modeling and computational issues encountered in the analytical work are discussed.

Mcardle, Jack G.↗

Analysis of a pneumatic forebody flow control concept about a full aircraft geometry

A full aircraft geometry is used to computationally analyze the effectiveness of a pneumatic forebody flow control concept. An overset grid technique is employed to model the aircraft and slot geometry. Steady-state solutions for both isolated forebody and full aircraft configurations are carried out using a thin-layer Navier-Stokes flow solver. A solution obtained using the full aircraft geometry and a flight sideslip condition investigates the effect of sideslip on the leading edge extention vortex burst point. A no-sideslip blowing solution using the isolated forebody at full-scale wind tunnel test conditions is compared with experimental data to determine the accuracy of the numerical method. A solution employing the full geometry and slot blowing at flight conditions is obtained.

Gee, Ken↗

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗

CFD analysis on control of secondary losses in STME LOX turbines with endwall fences

The rotor blade in the newly designed LOX turbine for the future Space Transportation Main Engine (STME) has a severe flow turning angle, nearly 160 degrees. The estimated secondary loss in the rotor alone accounts for nearly 50 percent of the total loss over the entire stage. To reduce such a loss, one of the potential methods is to use fences attached on the turbine endwall (hub). As a prelude to examining the effects of endwall fence with actual STME turbine configuration, the present study focuses on similar issues with a different, but more generic, geometry - a rectangular duct with a 160-degree bend. The duct cross-section has a 2-to-1 aspect ratio and the radii of curvature for the inner and outer wall are 0.25 and 1.25 times the duct width, respectively. The present emphasis lies in examining the effects of various fence-length extending along the streamwise direction. The flowfield is numerically simulated using the FDNS code developed earlier by Wang and Chen. The FDNS code is a pressure based, finite-difference, Navier-Stokes equations solver.

Chyu, Mingking K.↗

Radiation/convection coupling in rocket motor and plume analysis

A method for describing radiation/convection coupling to a flow field analysis was developed for rocket motors and plumes. The three commonly used propellant systems (H2/O2, RP-1/O2, and solid propellants) radiate primarily as: molecular emitters, non-scattering small particles (soot), and scattering larger particles (Al2O3), respectively. For the required solution, the divergence of the radiation heat flux was included in the energy equation, and the local, volume averaged intensity was determined by a solution to the radiative transfer equation. A rigorous solution to this problem is intractable, therefore, solution methods which use the ordinary and improved differential approximation were developed. This radiation model was being incorporated into the FDNS code, a Navier-Stokes flowfield solver for multiphase, turbulent combusting flows.

Saladino, A. J.↗

Juncture flow improvement for wing/pylon configurations by using CFD methodology

Transonic flow field around a fighter wing/pylon configuration was simulated by using an implicit upwinding Navier-Stokes flow solver (F3D) and overset grid technology (Chimera). Flow separation and local shocks near the wing/pylon junction were observed in flight and predicted by numerical calculations. A new pylon/fairing shape was proposed to improve the flow quality. Based on numerical results, the size of separation area is significantly reduced and the onset of separation is delayed farther downstream. A smoother pressure gradient is also obtained near the junction area. This paper demonstrates that computational fluid dynamics (CFD) methodology can be used as a practical tool for aircraft design.

Gea, Lie-Mine↗

Numerical vorticity capturing for vortex-solid body interaction problems

A new numerical method based on adding a term to the Euler/Navier-Stokes equations has been demonstrated it can effectively treat vortex-dominated flows using low-order numerical schemes and coarse grids. The numerical diffusion introduced by convection schemes can be eliminated by modifying the velocity field as a result of solving the modified flow governing equations. The modification of the velocity field, which conserves the total vorticity, essentially convects the vorticity toward its local extreme to offset the numerical diffusion. This method preserves the vortex structure even when vortices travel on coarse grid region. The method is implemented with two distinctive codes: a Navier-Stokes flow solver based on vorticity-velocity formulation and an Euler code based on velocity-pressure formulation. Problems presented in this paper include airfoil dynamic stall, vortex airfoil interaction, and vortex-fuselage interactions.

Wang, Clin M.↗

Numerical study of 3-D inducer and impeller for pump model development

Current design of high-performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study is to develop a robust and effective CFD pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, two key components of the turbopump, the inducer and impeller, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne.

Cheng, G. C.↗

Space marching calculations about hypersonic configurations using a solution-adaptive mesh algorithm

A solution-adaptive marching algorithm is developed and applied to a three-dimensional parabolized Navier-Stokes equation solver. The resulting algorithm obtains accurate solutions by using a spatial-marching/adaptive grid procedure. The adaptation step redistributes grid points line by line in both crossflow directions, with grid point motion controlled by forces analogous to tensional and torsional spring forces with the tensional force proportional to the error measure or weighting functions. The solution-adaptive marching procedure is applied to the hypersonic flow about two generic aircraft configurations. The first of these is an all-body-type geometry with elliptical cross sections and is studied at angles of attack of 0.5, and 15 deg. The second geometry is a generic blended-wing-body design. Results are presented that demonstrate the improvements in flowfield resolution obtainable with the solution-adaptive marching procedure over conventional fixed grid techniques. Computed pitot pressure profiles obtained using the solution-adaptive algorithm show improved agreement with experimental data compared to predictions obtained using a fixed grid.

Harvey, Albert D.↗

Thin-layer and full Navier-Stokes calculations for turbulent supersonic flow over a cone at an angle of attack

The proper use of a computational fluid dynamics code requires a good understanding of the particular code being applied. In this report the application of CFL3D, a thin-layer Navier-Stokes code, is compared with the results obtained from PARC3D, a full Navier-Stokes code. In order to gain an understanding of the use of this code, a simple problem was chosen in which several key features of the code could be exercised. The problem chosen is a cone in supersonic flow at an angle of attack. The issues of grid resolution, grid blocking, and multigridding with CFL3D are explored. The use of multigridding resulted in a significant reduction in the computational time required to solve the problem. Solutions obtained are compared with the results using the full Navier-Stokes equations solver PARC3D. The results obtained with the CFL3D code compared well with the PARC3D solutions.

Smith, Crawford F.↗

Time-accurate Navier-Stokes calculations with multigrid acceleration

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.

Melson, N. Duane↗