SEARCH · Engineering Papers

Results for “Finite Volume”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Base-Bleed Effect on X-33 Aerospike Plume Induced Base-Heating Environment During Power-Pack Out

A computational heat transfer methodology was developed to study the dual-engine linear aerospike plume induced base-heating environment during one power-pack out, in ascent flight. One power-pack out results in reduction of power levels for both engines. That, in turn, reduces the amount of base-bleed and changes the distribution of base-bleed on the two pillows. Hence, the concern of increased base-heating during power-pack out. The thermo-flowfield of the entire vehicle was computed. The computational methodology for the convective heating is based on a three-dimensional, finite-volume, viscous, chemically reacting, and pressure-based computational fluid dynamics formulation. The computational methodology for the radiative heating is based on a three-dimensional, finite-volume, and spectral-line-based weighted-sum-of-gray-gases absorption computational radiation heat transfer formulation. A separate radiation model was used for diagnostic purposes. The computational methodology was systematically benchmarked. In this study, near-base radiative heat fluxes were computed and they compared well with those measured from an installed linear aerospike engine tests. The base-heating environment of 18 trajectory points selected from three power-pack out ascent scenarios was computed and is presented here. The power-pack out condition has the most impact on convective base-heating when it happens early in flight. The some of its impact comes from the asymmetric and reduced base-bleed.

Wang, Tee-See↗

NASA NTRS1986

Embedded mesh solution of the 2-D Euler equations - Evaluation of interface formulations

Solution of the steady 2-D Euler equations using mesh embedding, or local grid refinement, with a cell-centered finite volume scheme is investigated. Embedded regions which are topologically similar to the global grid are considered. An isoenergetic model for the governing equations is used in Jameson's finite volume multistage scheme with modifications to the boundary conditions and smoothing. A detailed study of the embedding interface flux and smoothing formulations is conducted. Taylor expansion analysis reveals that local second order spatial accuracy is not possible if a conservative interface flux formulation is used. The analysis also gives constraints for local first order accuracy. An energy stability analysis indicates that downwind weighting of interface fluxes causes local instabilities. Analysis shows that conservative interface smoothing formulations must have a locally convective component, but that correct interface formulations allow globally dissipative smoothing. Embedded mesh solutions obtained with this scheme are presented for a transonic airfoil. They show that if embedding interfaces are close to the shocks, then small modifications in the interface location can have large effects on converge and solution accuracy.

Allmaras, S. R.↗

NASA NTRS1990

Computation and control of asymmetric vortex flow around circular cones using Navier-Stokes equations

The unsteady, compressible, thin-layer and full Navier-Stokes equations are used to numerically simulate steady and unsteady asymmetric, supersonic, locally conical flows around a 5-deg semiapex angle circular cone. The main computational scheme is the implicit, upwind, flux-difference splitting, finite-volume scheme. Comparison of asymmetric flow solutions using the thin-layer and full Navier-Stokes equations is presented and discussed. The implicit, upwind, flux-vector splitting, finite-volume scheme has also been used to solve for the unsteady asymmetric flow with vortex shedding. The unsteady-flow solution using the flux-vector splitting scheme perfectly agrees with the previously obtained solution using the flux-difference splitting scheme. Passive control of asymmetric flows has been demonstrated and studied using sharp- and round-edged, thick and thin strakes.

Kandil, Osama A.↗

NASA NTRS2012

Experimenting with the GMAO 4D Data Assimilation

The Global Modeling and Assimilation Office (GMAO) has been working to promote its prototype four-dimensional variational (4DVAR) system to a version that can be exercised at operationally desirable configurations. Beyond a general circulation model (GeM) and an analysis system, traditional 4DV AR requires availability of tangent linear (TL) and adjoint (AD) models of the corresponding GeM. The GMAO prototype 4DVAR uses the finite-volume-based GEOS GeM and the Grid-point Statistical Interpolation (GSI) system for the first two, and TL and AD models derived ITom an early version of the finite-volume hydrodynamics that is scientifically equivalent to the present GEOS nonlinear GeM but computationally rather outdated. Specifically, the TL and AD models hydrodynamics uses a simple (I-dimensional) latitudinal MPI domain decomposition, which has consequent low scalability and prevents the prototype 4DV AR ITom being used in realistic applications. In the near future, GMAO will be upgrading its operational GEOS GCM (and assimilation system) to use a cubed-sphere-based hydrodynamics. This versions of the dynamics scales to thousands of processes and has led to a decision to re-derive the TL and AD models for this more modern dynamics, thus taking advantage of a two-dimensional MPI decomposition and improved scalability properties. With the aid of the Transformation of Algorithms in FORTRAN (l'AF) automatic adjoint generation tool and some hand-coding, a version of the cubed-sphere-based TL and AD models, with a simplified vertical diffusion scheme, is now available, enabling multiple configurations of standard implementations of 4DV AR in GEOS. Concurrent to this development, collaboration with the National Centers for Environmental Prediction (NCEP) and the Earth System Research Laboratory (ESRL) has allowed GMAO to implement a hybrid-ensemble capability within the GEOS data assimilation system. Both 3Dand 4D-ensemble capabilities are presently available thus allowing GMAO to now evaluate the performance and benefit of various ensemble and variational assimilation strategies. This presentation will cover the most recent developments taking place at GMAO and show results from various comparisons from traditional techniques to more recent ensemble-based ones.

Todling, R.↗

NASA NTRS2003

Retrospective Data Assimilation for GEOS-4

The NASA/Goddard Earth Observing System 3 (GEOS-3) has recently been replaced operationally by GEOS-4. One of the major components of the new system is the finite-volume general circulation model (GCM) of Lin and Rood. Updating the GEOS-3 retrospective data assimilation system of Zhu et al. to GEOS-4 requires two main steps: (i) developing the adjoint of the finite-volume GCM; and (ii) developing the adjoint of the procedures interfacing the GCM with the Physical-space Statistical Analysis System. These required adjoints are being developed automatically through the use of the Transformation of Algorithms in Fortran compiler of Giering and Kaminski. In this presentation we discuss the current status of development and plans for the near future as well as show preliminary results.

Todling, Ricardo↗

NASA NTRS1995

Thermocapillary convection in a rectangular cavity with a deformable interface

A finite-volume method and a boundary element technique are used to compute two-dimensional (2D) thermocapillary convection in a rectangular cavity. The free surface is deformable and the deviations from a flat interface, h, are assumed small in the finite-volume calculations as appropriate for small Capillary (Ca). Two asymptotic approaches are employed; the first is an expansion valid as Ca tends towards 0 and the second assumes h tends towards 0, retaining Ca explicitly as a parameter. On the other hand, the boundary elements approach can be used with O(1) surface deformations. These three different formulations are used to calculate thermocapillary motions in fluids with small Prandtl numbers of about 0.01, Ca less than 0.05, aspect ratios (width/height) 1, 2 and 4, and various values of the Maragoni number (Ma). The same solutions are calculated with these different approaches and are found in good agreement for values of Ca up to 0.05. All solutions calculated are steady, which is both in agreement and disagreement with recently published results by different authors employing different numerical techniques.

Mundrane, M.↗

NASA NTRS1994

Probability density function approach for compressible turbulent reacting flows

The objective of the present work is to extend the probability density function (PDF) tubulence model to compressible reacting flows. The proability density function of the species mass fractions and enthalpy are obtained by solving a PDF evolution equation using a Monte Carlo scheme. The PDF solution procedure is coupled with a compression finite-volume flow solver which provides the velocity and pressure fields. A modeled PDF equation for compressible flows, capable of treating flows with shock waves and suitable to the present coupling scheme, is proposed and tested. Convergence of the combined finite-volume Monte Carlo solution procedure is discussed. Two super sonic diffusion flames are studied using the proposed PDF model and the results are compared with experimental data; marked improvements over solutions without PDF are observed.

Hsu, A. T.↗

NASA NTRS2005

A Robust Absorbing Boundary Condition for Compressible Flows

An absorbing non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented with theoretical proof. This paper is a continuation and improvement of a previous paper by the author. The absorbing NRBC technique is based on a first principle of non reflecting, which contains the essential physics that a plane wave solution of the Euler equations remains intact across the boundary. The technique is theoretically shown to work for a large class of finite volume approaches. When combined with the hyperbolic conservation laws, the NRBC is simple, robust and truly multi-dimensional; no additional implementation is needed except the prescribed physical boundary conditions. Several numerical examples in multi-dimensional spaces using two different finite volume schemes are illustrated to demonstrate its robustness in practical computations. Limitations and remedies of the technique are also discussed.

Loh, Ching Y.↗

NASA NTRS1992

Gigaflop (billion floating point operations per second) performance for computational electromagnetics

Accurate and rapid evaluation of radar signature for alternative aircraft/store configurations would be of substantial benefit in the evolution of integrated designs that meet radar cross-section (RCS) requirements across the threat spectrum. Finite-volume time domain methods offer the possibility of modeling the whole aircraft, including penetrable regions and stores, at longer wavelengths on today's gigaflop supercomputers and at typical airborne radar wavelengths on the teraflop computers of tomorrow. A structured-grid finite-volume time domain computational fluid dynamics (CFD)-based RCS code has been developed at the Rockwell Science Center, and this code incorporates modeling techniques for general radar absorbing materials and structures. Using this work as a base, the goal of the CFD-based CEM effort is to define, implement and evaluate various code development issues suitable for rapid prototype signature prediction.

Shankar, V.↗

NASA NTRS1989

Grid generation and inviscid flow computation about aircraft geometries

Grid generation and Euler flow about fighter aircraft are described. A fighter aircraft geometry is specified by an area ruled fuselage with an internal duct, cranked delta wing or strake/wing combinations, canard and/or horizontal tail surfaces, and vertical tail surfaces. The initial step before grid generation and flow computation is the determination of a suitable grid topology. The external grid topology that has been applied is called a dual-block topology which is a patched C (exp 1) continuous multiple-block system where inner blocks cover the highly-swept part of a cranked wing or strake, rearward inner-part of the wing, and tail components. Outer-blocks cover the remainder of the fuselage, outer-part of the wing, canards and extend to the far field boundaries. The grid generation is based on transfinite interpolation with Lagrangian blending functions. This procedure has been applied to the Langley experimental fighter configuration and a modified F-18 configuration. Supersonic flow between Mach 1.3 and 2.5 and angles of attack between 0 degrees and 10 degrees have been computed with associated Euler solvers based on the finite-volume approach. When coupling geometric details such as boundary layer diverter regions, duct regions with inlets and outlets, or slots with the general external grid, imposing C (exp 1) continuity can be extremely tedious. The approach taken here is to patch blocks together at common interfaces where there is no grid continuity, but enforce conservation in the finite-volume solution. The key to this technique is how to obtain the information required for a conservative interface. The Ramshaw technique which automates the computation of proportional areas of two overlapping grids on a planar surface and is suitable for coding was used. Researchers generated internal duct grids for the Langley experimental fighter configuration independent of the external grid topology, with a conservative interface at the inlet and outlet.

Smith, Robert E.↗

NASA NTRS1991

Shock capturing schemes for multidimensional flow

Progress made in the development of a genuinely multidimensional finite volume algorithm for problems in inviscid gas dynamics is presented. The approach entails: (1) the reconstruction of flowfield data using a planar wave pattern in which the strengths and orientations of the component waves are derived independently of the mesh geometry; and (2) the development of a flux formula which provides a numerical approximation to the flux at a finite volume cell face during the passage of waves which are in general oblique to the face. Several algorithms are outlined, and the most recent developments are included. The results of several numerical cases are also included.

Parpia, Ijaz H.↗

NASA NTRS1990

Flow field analysis for a class of waverider configurations

A package of computer codes for analysis of flow fields for waverider configurations is described. The package consists of a surface/volume grid generator and a finite-volume flow solver. The grid generator defines body geometries and computational grids by an algebraic homotopy procedure. The algebraic procedure is versatile in its application and can readily generate configurations in the class of blended wing-body geometries. This code has the ability to produce a wide variety of geometries in the given class with varying geometrical attributes. The flow solver employs a finite-volume formation and solves the explicit, Runge-Kutta integration technique. The method or flow simulation incorporates several techniques for acceleration of the convergence of the interaction process and an entropy corrected enthalpy damping procedure for efficient computation of high Mach number flows.

Moitra, Anutosh↗

NASA NTRS2003

Numerical Modeling of Fluid Transient in Cryogenic Fluid Network of Rocket Propulsion System

Fluid transients, also known as water hammer, can have a significant impact on the design and operation of both spacecraft and launch vehicles propulsion systems. These transients often occur at system activation and shut down. For ground safety reasons, many spacecrafts are launched with the propellant lines dry. These lines are often evacuated by the time the spacecraft reaches orbit. When the propellant isolation valve opens during propulsion system activation, propellant rushes into lines creating a pressure surge. During propellant system shutdown, a pressure surge is created due to sudden closure of a valve. During both activation and shutdown, pressure surges must be predicted accurately to ensure structural integrity of the propulsion system fluid network. The method of characteristics is the most widely used method of calculating fluid transients in pipeline [ 1,2]. The method of characteristics, however, has limited applications in calculating flow distribution in complex flow circuits with phase change, heat transfer and rotational effects. A robust cryogenic propulsion system analyzer must have the capability to handle phase change, heat transfer, chemical reaction, rotational effects and fluid transients in conjunction with subsystem flow model for pumps, valves and various pipe fittings. In recent years, such a task has been undertaken at Marshall Space Flight Center with the development of the Generalized Fluid System Simulation Program (GFSSP), which is based on finite volume method in fluid network [3]. GFSSP has been extensively verified and validated by comparing its predictions with test data and other numerical methods for various applications such as internal flow of turbo-pump [4], propellant tank pressurization [5,6], chilldown of cryogenic transfer line [7] and squeeze film damper rotordynamics [8]. The purpose of the present paper is to investigate the applicability of the finite volume method to predict fluid transient in cryogenic flow circuits.

Majumdar, Alok↗

NASA NTRS1991

Free stream capturing in fluid conservation law for moving coordinates in three dimensions

The free-stream capturing technique for both the finite-volume (FV) and finite-difference (FD) framework is summarized. For an arbitrary motion of the grid, the FV analysis shows that volumes swept by all six surfaces of the cell have to be computed correctly. This means that the free-stream capturing time-metric terms should be calculated not only from a surface vector of a cell at a single time level, but also from a volume swept by the cell surface in space and time. The FV free-stream capturing formulation is applicable to the FD formulation by proper translation from an FV cell to an FD mesh.

Obayashi, Shigeru↗

NASA NTRS2008

The Global Modeling and Assimilation Office (GMAO) 4d-Var and its Adjoint-based Tools

The fifth generation of the Goddard Earth Observing System (GEOS-5) Data Assimilation System (DAS) is a 3d-var system that uses the Grid-point Statistical Interpolation (GSI) system developed in collaboration with NCEP, and a general circulation model developed at Goddard, that includes the finite-volume hydrodynamics of GEOS-4 wrapped in the Earth System Modeling Framework and physical packages tuned to provide a reliable hydrological cycle for the integration of the Modern Era Retrospective-analysis for Research and Applications (MERRA). This MERRA system is essentially complete and the next generation GEOS is under intense development. A prototype next generation system is now complete and has been producing preliminary results. This prototype system replaces the GSI-based Incremental Analysis Update procedure with a GSI-based 4d-var which uses the adjoint of the finite-volume hydrodynamics of GEOS-4 together with a vertical diffusing scheme for simplified physics. As part of this development we have kept the GEOS-5 IAU procedure as an option and have added the capability to experiment with a First Guess at the Appropriate Time (FGAT) procedure, thus allowing for at least three modes of running the data assimilation experiments. The prototype system is a large extension of GEOS-5 as it also includes various adjoint-based tools, namely, a forecast sensitivity tool, a singular vector tool, and an observation impact tool, that combines the model sensitivity tool with a GSI-based adjoint tool. These features bring the global data assimilation effort at Goddard up to date with technologies used in data assimilation systems at major meteorological centers elsewhere. Various aspects of the next generation GEOS will be discussed during the presentation at the Workshop, and preliminary results will illustrate the discussion.

Todling, Ricardo↗

NASA NTRSDate not supplied

Validation of the Hypersolve CFD Solver for Entry Descent and Landing Applications

The functional equivalence of the HyperSolve unstructured edge-based, finite-volume computational fluid dynamics code to the Langley Aerothermodynamic Upwind Relaxation Algorithm multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing community. A suite of cases using a range of thermochemical gas models on relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flowfield quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. The functional equivalence of the HyperSolve unstructured edge-based finite-volume computational fluid dynamics (CFD) code to the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) multiblock structured grid code is documented for applications of interest to the Entry, Descent, and Landing (EDL) community. A suite of cases using a range of thermochemical gas models on EDL-relevant vehicle configurations were analyzed with both codes and the results compared. A tolerance of ±4% difference in surface pressure and surface heat flux from a benchmark LAURA solution was used as the criterion for functional equivalence, and comparisons of flow field quantities are also included to verify that the thermochemical nonequilibrium capabilities in HyperSolve match those of the LAURA code. In general, HyperSolve predictions for surface pressure and surface heat flux are in close agreement with those predicted by LAURA.

hypersolve↗

NASA NTRSDate not supplied

Validation of the HyperSolve CFD Solver for Entry Descent and Landing Applications

hypersolve↗

NASA NTRS1992

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗