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 1,333 records · Page 74

All sources NASA NTRS DOE OSTI

NASA NTRS2007

URANS Application with CFL3D

This case was run using CFL3D, a multi-zone Reynolds-averaged Navier-Stokes code developed at NASA Langley [1]. It solves the thin-layer form of the Navier-Stokes equations in each of the (selected) coordinate directions. It can use 1-to-1, patched, or overset grids, and employs local time step scaling, grid sequencing, and multigrid to accelerate convergence to steady state. In time-accurate mode, CFL3D has the option to employ dual-time stepping with subiterations and multigrid, and it achieves second order temporal accuracy. CFL3D is a finite volume method. It uses third-order upwind-biased spatial differencing on the convective and pressure terms, and second-order differencing on the viscous terms; it is globally second-order spatially accurate. The flux difference-splitting (FDS) method of Roe is employed to obtain fluxes at the cell faces. It is advanced in time with an implicit three-factor approximate factorization method.

Rumsey, C. L.↗

NASA NTRS2007

Simulation of a Periodic Jet in a Crossflow with a RANS Solver Using an Unstructured Grid

A second-order unstructured-grid code, developed and used primarily for steady aerodynamic simulations, is applied to the synthetic jet in a cross flow. The code, FUN3D, is a vertex-centered finite-volume method originally developed by Anderson[1, 2], and is currently supported by members of the Fast Adaptive Aerospace Tools team at NASA Langley. Used primarily for design[3] and analysis[4] of steady aerodynamic configurations, FUN3D incorporates a discrete adjoint capability, and supports parallel computations using MPI. A detailed description of the FUN3D code can be found in the references given above. The code is under continuous development and contains a variety of flux splitting algorithms for the inviscid terms, two methods for computing gradients, several turbulence models, and several solution methodologies; all in varying states of development. Only the most robust and reliable components, based on experiences with steady aerodynamic simulations, were employed in this work. As applied in this work, FUN3D solves the Reynolds averaged Navier-Stokes equations using the one equation turbulence model of Spalart and Allmaras[5]. The spatial discretization is formed on unstructured meshes using a vertex-centered approach. The inviscid terms are evaluated by a flux-difference splitting formulation using least-squares reconstruction and Roe-type approximate Riemann fluxes. Green-Gauss gradient evaluations are used for viscous and turbulence modeling terms. The discrete spatial operator is combined with a backward time operator which is then solved iteratively using point or line Gauss-Seidel and local time stepping in a pseudo time. For steady flows, the physical time step is set to infinity and the pseudo time step is ramped up with the iteration count. A second-order backward in time operator is used for time accurate flows with 20 to 50 steps in the pseudo time applied at each physical time step. For this effort, FUN3D was modified to support spatially varying boundary and initial conditions, and unsteady boundary conditions. Also, a specialized in/out flow boundary condition was implemented to model the action of the diaphragm. This boundary condition is described below in more detail. The grids were generated using the internally developed codes GridEX[6] for meshing the surfaces and inviscid regions of the domain, and for CAD access; and MesherX[7] for meshing the viscous regions. Grid spacing in on the surfaces and in the inviscid regions are indirectly controlled by specifying sources. The viscous layers are generated using an advancing layer technique. MeshersX allows the user to control the spatial variation of the first step off the surface, growth rates, and the termination criterion by providing small problem dependent subroutines.

Atkins, H. L.↗

NASA NTRS2007

RANS and URANS Application with CFL3D

Rumsey, C. L.↗

NASA NTRS2006

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing during Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on numerical simulations with the NASA Global Modeling and Assimilation Office finite-volume general circulation model (fvGCM) with Microphyics of Clouds with the Relaxed Arakawa Schubert Scheme (McRAS), using aerosol forcing functions derived from the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The authors find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excite a planetary-scale teleconnection pattern in sea level pressure, temperature, and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, that is, South Asia, East Asia, and northern and western Africa. Significant atmospheric heating is found in regions with large loading of dust (over northern Africa and the Middle East) and black carbon (over Southeast Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east west dipole anomaly with strong cooling over the Caspian Sea and warming over central and northeastern Asia, where aerosol concentrations are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection pattern driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations. The surface temperature signature associated with the aerosol-induced teleconnection bears striking resemblance to the spatial pattern of observed long-term trend in surface temperature over Eurasia. Additionally, the boreal spring wave train pattern is similar to that reported by Fukutomi et al. associated with the boreal summer precipitation seesaw between eastern and western Siberia. The results of this study raise the possibility that global aerosol forcing during boreal spring may play an important role in spawning atmospheric teleconnections that affect regional and global climates.

Kim, Maeng-Ki↗

NASA NTRS2006

Asian Summer Monsoon Anomalies Induced by Aerosol Direct Forcing: The Role of the Tibetan Plateau

In this paper we present results of a numerical study using the NASA finite-volume GCM to elucidate a plausible mechanism for aerosol impact on the Asian summer monsoon involving interaction with physical processes over the Tibetan Plateau (TP). During the premonsoon season of March April, dusts from the deserts of western China, Afghanistan/Pakistan, and the Middle East are transported into and stacked up against the northern and southern slopes of the TP. The absorption of solar radiation by dust heats up the elevated surface air over the slopes. On the southern slopes, the atmospheric heating is reinforced by black carbon from local emission. The heated air rises via dry convection, creating a positive temperature anomaly in the mid-to-upper troposphere over the TP relative to the region to the south. In May through early June in a manner akin to an elevated heat pump , the rising hot air forced by the increasing heating in the upper troposphere, draws in warm and moist air over the Indian subcontinent, setting the stage for the onset of the South Asia summer monsoon. Our results suggest that increased dust loading coupled with black carbon emission from local sources in northern India during late spring may lead to an advance of the rainy periods and subsequently an intensification of the Indian summer monsoon. The enhanced rainfall over India is associated with the development of an aerosol-induced large-scale sea level pressure anomaly pattern, which causes the East Asia (Mei-yu) rain belt to shift northwestward, suppressing rainfall over East Asia and the adjacent oceanic regions.

Lau, K. M.↗

NASA NTRS2007

Modeling of Compressible Flow with Friction and Heat Transfer Using the Generalized Fluid System Simulation Program (GFSSP)

The present paper describes the verification and validation of a quasi one-dimensional pressure based finite volume algorithm, implemented in Generalized Fluid System Simulation Program (GFSSP), for predicting compressible flow with friction, heat transfer and area change. The numerical predictions were compared with two classical solutions of compressible flow, i.e. Fanno and Rayleigh flow. Fanno flow provides an analytical solution of compressible flow in a long slender pipe where incoming subsonic flow can be choked due to friction. On the other hand, Raleigh flow provides analytical solution of frictionless compressible flow with heat transfer where incoming subsonic flow can be choked at the outlet boundary with heat addition to the control volume. Nonuniform grid distribution improves the accuracy of numerical prediction. A benchmark numerical solution of compressible flow in a converging-diverging nozzle with friction and heat transfer has been developed to verify GFSSP's numerical predictions. The numerical predictions compare favorably in all cases.

Bandyopadhyay, Alak↗

NASA NTRS2006

Liquid Annular Seal CFD Analysis for Rotordynamic Force Prediction

A commercially available code is utilized to analyze a plain and grooved liquid annular seal. These type seals are commonly used in modern turbopumps and have a pronounced effect on the rotordynamic behavior of these systems. Accurate prediction of both leakage and dynamic reaction forces is vital to ensure good performance and sound mechanical operation. The code SCISEAL developed by CFDRC is a generic 3-D, finite volume based CFD code solving the 3-D Reynolds averaged Navier Stokes equations. The code allows body-fitted, multi-blocked structured grids, turbulence modeling, rotating coordinate frames, as well as integration of dynamic pressure and shear forces on the rotating journal. The code may be used with the commercially available pre-and post-processing codes from CFDRC as well.

Moore, Jeff↗

NASA NTRS2007

Robustness of Modeling of Out-of-Service Gas Mechanical Face Seal

Gas lubricated mechanical face seal are ubiquitous in many high performance applications such as compressors and gas turbines. The literature contains various analyses of seals having orderly face patterns (radial taper, waves, spiral grooves, etc.). These are useful for design purposes and for performance predictions. However, seals returning from service (or from testing) inevitably contain wear tracks and warped faces that depart from the aforementioned orderly patterns. Questions then arise as to the heat generated at the interface, leakage rates, axial displacement and tilts, minimum film thickness, contact forces, etc. This work describes an analysis of seals that may inherit any (i.e., random) face pattern. A comprehensive computer code is developed, based upon the Newton- Raphson method, which solves for the equilibrium of the axial force and tilting moments that are generated by asperity contact and fluid film effects. A contact mechanics model is incorporated along with a finite volume method that solves the compressible Reynolds equation. Results are presented for a production seal that has sustained a testing cycle.

Green, Itzhak↗

NASA NTRS2008

Towards Verification of Unstructured-Grid Solvers

New methodology for verification of computational methods using unstructured grids is presented. The discretization order properties are studied in computational windows, easily constructed within a collection of grids or a single grid. The windows can be adjusted to isolate the interior discretization, the boundary discretization, or singularities. A major component of the methodology is the downscaling test, introduced previously for studying the convergence rates of truncation and discretization errors of finite-volume discretization schemes on general unstructured grids. Demonstrations of the method are shown, including a comparative accuracy assessment of commonly-used schemes on general mixed grids and the identification of local accuracy deterioration at intersections of tangency and inflow/outflow boundaries. Recommendations for the use of the methodology in large-scale computational simulations are given.

Thomas, James L.↗

NASA NTRS2008

Direct-View Multi-Point Two-Component Interferometric Rayleigh Scattering Velocimeter

This paper describes an instantaneous velocity measurement system based on the Doppler shift of elastically scattered laser light from gas molecules (Rayleigh scattering) relative to an incident laser. The system uses a pulsed laser as the light source, direct-viewing optics to collect the scattered light, an interferometer to analyze spectrally the scattered light mixed with the incident laser light, and a CCD camera to capture the resulting interferogram. The system is capable of simultaneous, spatially (approximately 0.2 mm(exp 3)) and temporally (approximately 40 ns) resolved, multiple point measurements of two orthogonal components of flow velocity in the presence of background scattered light, acoustic noise and vibrations, and flow particulates. Measurements in a large-scale axi-symmetric Mach 1.6 H2-air combustion-heated jet running at a flow sensible enthalpy specific to Mach 5.5 hypersonic flight are performed to demonstrate the technique. The measurements are compared with CFD calculations using a finite-volume discretization of the Favre-averaged Navier-Stokes equations (VULCAN code).

Bivolaru, Daniel↗

NASA NTRS2008

Implementation of a Blowing Boundary Condition in the LAURA Code

Preliminary steps toward modeling a coupled ablation problem using a finite-volume Navier-Stokes code (LAURA) are presented in this paper. Implementation of a surface boundary condition with mass transfer (blowing) is described followed by verification and validation through comparisons with analytic results and experimental data. Application of the code to a carbon-nosetip ablation problem is demonstrated and the results are compared with previously published data. It is concluded that the code and coupled procedure are suitable to support further ablation analyses and studies.

Thompson, Richard a.↗

NASA NTRS2008

Near-Field Noise Computation for a Subsonic Coannular Jet

A high-Reynolds-number, subsonic coannular jet is simulated, using a three-dimensional finite-volume LES method, with emphasis on the near field noise. The nozzle geometry used is the NASA Glenn 3BB baseline model. The numerical results are generally in good agreement with existing experimental findings.

Loh, Ching Y.↗

NASA NTRS2007

Heat Transfer in High Temperature Multilayer Insulation

High temperature multilayer insulations have been investigated as an effective component of thermal-protection systems for atmospheric re-entry of reusable launch vehicles. Heat transfer in multilayer insulations consisting of thin, gold-coated, ceramic reflective foils and Saffil(TradeMark) fibrous insulation spacers was studied both numerically and experimentally. A finite volume numerical thermal model using combined conduction (gaseous and solid) and radiation in porous media was developed. A two-flux model with anisotropic scattering was used for radiation heat transfer in the fibrous insulation spacers between the reflective foils. The thermal model was validated by comparison with effective thermal conductivity measurements in an apparatus based on ASTM standard C201. Measurements were performed at environmental pressures in the range from 1x10(exp -4) to 760 torr over the temperature range from 300 to 1300 K. Four multilayer samples with nominal densities of 48 kg/cu m were tested. The first sample was 13.3 mm thick and had four evenly spaced reflective foils. The other three samples were 26.6 mm thick and utilized either one, two, or four reflective foils, located near the hot boundary with nominal foil spacing of 1.7 mm. The validated thermal model was then used to study relevant design parameters, such as reflective foil spacing and location in the stack-up and coating of one or both sides of foils.

Daryabeigi, Kamran↗

NASA NTRS2008

A Goddard Multi-Scale Modeling System with Unified Physics

A multi-scale modeling system with unified physics has been developed at NASA Goddard Space Flight Center (GSFC). The system consists of an MMF, the coupled NASA Goddard finite-volume GCM (fvGCM) and Goddard Cumulus Ensemble model (GCE, a CRM); the state-of-the-art Weather Research and Forecasting model (WRF) and the stand alone GCE. These models can share the same microphysical schemes, radiation (including explicitly calculated cloud optical properties), and surface models that have been developed, improved and tested for different environments. The following is presented in this report: (1) a brief review of the GCE model and its applications on the impact of aerosols on deep precipitation processes, (2) the Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), and preliminary results (the comparison with traditional GCMs), and (3) a discussion on the Goddard WRF version (its developments and applications).

Tao, Wei-Kuo↗

NASA NTRS2006

A Multi-scale Modeling System: Developments, Applications and Critical Issues

A multi-scale modeling framework (MMF), which replaces the conventional cloud parameterizations with a cloud-resolving model (CRM) in each grid column of a GCM, constitutes a new and promising approach. The MMF can provide for global coverage and two-way interactions between the CRMs and their parent GCM. The GCM allows global coverage and the CRM allows explicit simulation of cloud processes and their interactions with radiation and surface processes. A new MMF has been developed that is based the Goddard finite volume GCM (fvGCM) and the Goddard Cumulus Ensemble (GCE) model. This Goddard MMF produces many features that are similar to another MMF that was developed at Colorado State University (CSU), such as an improved .surface precipitation pattern, better cloudiness, improved diurnal variability over both oceans and continents, and a stronger, propagating Madden-Julian oscillation (MJO) compared to their parent GCMs using conventional cloud parameterizations. Both MMFs also produce a precipitation bias in the western Pacific during Northern Hemisphere summer. However, there are also notable differences between two MMFs. For example, the CSU MMF simulates less rainfall over land than its parent GCM. This is why the CSU MMF simulated less overall global rainfall than its parent GCM. The Goddard MMF overestimates global rainfall because of its oceanic component. Some critical issues associated with the Goddard MMF are presented in this paper.

Tao, Wei-Kuo↗

NASA NTRS2008

Adjoint-Based Methodology for Time-Dependent Optimization

This paper presents a discrete adjoint method for a broad class of time-dependent optimization problems. The time-dependent adjoint equations are derived in terms of the discrete residual of an arbitrary finite volume scheme which approximates unsteady conservation law equations. Although only the 2-D unsteady Euler equations are considered in the present analysis, this time-dependent adjoint method is applicable to the 3-D unsteady Reynolds-averaged Navier-Stokes equations with minor modifications. The discrete adjoint operators involving the derivatives of the discrete residual and the cost functional with respect to the flow variables are computed using a complex-variable approach, which provides discrete consistency and drastically reduces the implementation and debugging cycle. The implementation of the time-dependent adjoint method is validated by comparing the sensitivity derivative with that obtained by forward mode differentiation. Our numerical results show that O(10) optimization iterations of the steepest descent method are needed to reduce the objective functional by 3-6 orders of magnitude for test problems considered.

Yamaleev, N. K.↗

NASA NTRS2009

Design Evolution and Verification of the A-3 Chemical Steam Generator

Following is an overview of the Chemical Steam Generator system selected to provide vacuum conditions for a new altitude test facility, the A-3 Test Stand at Stennis Space Center (SSC) in Bay St. Louis, MS. A-3 will serve as NASA s primary facility for altitude testing of the J-2X rocket engine, to be used as the primary propulsion device for the upper stages of the Ares launch vehicles. The Chemical Steam Generators (CSGs) will produce vacuum conditions in the test cell through the production and subsequent supersonic ejection of steam into a diffuser downstream of the J-2X engine nozzle exit. The Chemical Steam Generators chosen have a rich heritage of operation at rocket engine altitude test facilities since the days of the Apollo program and are still in use at NASA White Sands Test Facility (WSTF) in New Mexico. The generators at WSTF have been modified to a degree, but are still very close to the heritage design. The intent for the A-3 implementation is to maintain this heritage design as much as possible, making minimal updates only where necessary to substitute for obsolete parts and to increase reliability. Reliability improvements are especially desired because the proposed system will require 27 generators, which is nine times the largest system installed in the 1960s. Improvements were suggested by the original design firm, Reaction Motors, by NASA SSC and NASA WSTF engineers, and by the A-3 test stand design contractor, Jacobs Technology, Inc. (JTI). This paper describes the range of improvements made to the design to date, starting with the heritage generator and the minor modifications made over time at WSTF, to the modernized configuration which will be used at A-3. The paper will discuss NASA s investment in modifications to SSC s E-2 test facility fire a full-scale Chemical Steam Generator in advance of the larger steam system installation at A-3. Risk mitigation testing will be performed in early 2009 at this test facility to verify that the CSGs operate as expected. The generator which will undergo this testing is of the most recent A-3 configuration, and will be instrumented far in excess of what is normally required for operation. The extra data will allow for easier troubleshooting and more complete knowledge of expected generator performance. In addition, the early testing will give SSC personnel experience in operating the CSG systems, which will expedite the process of installation and activation at A-3. Each Chemical Steam Generator is supported by a complement of valves, instruments, and flow control devices, with the entire assembly called a "module." The generators will be installed in groups of three, historically called "units". A module is so called because of its modular ability to be replaced or serviced without disturbing the other two modules installed on the same unit. A module is pictured in Figure 1, shown with its generator secured by white bands in its shipping (vs. installed) configuration. The heritage system at WSTF is composed of a single unit (three generator modules), pictured in Figure 2 as it was installed in 1965. In contrast, A-3 will have nine units operating in parallel to achieve vacuum conditions appropriate for testing the J-2X engine. Each of the combustors operates in two modes and achieves the so-called "full-steam" mode after all three of its stages ignite. Ignition of the first stage is achieved by exciting a spark plug; the second stage and main stage are lit by the flame front of the previous stage. The main stage burns approximately 97% of the total propellant flow and uses the heat energy to vaporize water into superheated steam. While the main stage remains unlit, the combustor is in so-called "idle" mode. In the WSTF system, this idle mode is not optimized for water usage, and does not need to be, as the water is pumped from a large reservoir. The water supply at A-3 will be contained in tanks with finite volume, so water optimization is preferred for the modnized configuration. Multiple solutions for this issue have been proposed, with the leading concept being a change to the operational definition of "idle mode," with the generator running in a lower heat flux condition.

Kirchner, Casey K.↗

NASA NTRS2008

Towards Verification of Unstructured-Grid Solvers

New methodology for verification of finite-volume computational methods using unstructured grids is presented. The discretization order properties are studied in computational windows, easily constructed within a collection of grids or a single grid. Tests are performed within each window and address a combination of problem-, solution-, and discretization/grid-related features affecting discretization error convergence. The windows can be adjusted to isolate particular elements of the computational scheme, such as the interior discretization, the boundary discretization, or singularities. Studies can use traditional grid-refinement computations within a fixed window or downscaling, a recently-introduced technique in which computations are made within windows contracting toward a focal point of interest. Grids within the windows are constrained to be consistently refined, allowing a meaningful assessment of asymptotic error convergence on unstructured grids. Demonstrations of the method are shown, including a comparative accuracy assessment of commonly-used schemes on general mixed grids and the identification of local accuracy deterioration at boundary intersections. Recommendations to enable attainment of design-order discretization errors for large-scale computational simulations are given.

Thomas, James L.↗