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,279 records · Page 71

All sources NASA NTRS DOE OSTI

NASA NTRS2003

TOMS and SBUV Data: Comparison to 3D Chemical-Transport Model Results

We have updated our merged ozone data (MOD) set using the TOMS data from the new version 8 algorithm. We then analyzed these data for contributions from solar cycle, volcanoes, QBO, and halogens using a standard statistical time series model. We have recently completed a hindcast run of our 3D chemical-transport model for the same years. This model uses off-line winds from the finite-volume GCM, a full stratospheric photochemistry package, and time-varying forcing due to halogens, solar uv, and volcanic aerosols. We will report on a parallel analysis of these model results using the same statistical time series technique as used for the MOD data.

Stolarski, Richard S.↗

NASA NTRS2004

Polar Processes in a 50-year Simulation of Stratospheric Chemistry and Transport

The unique chemical, dynamical, and microphysical processes that occur in the winter polar lower stratosphere are expected to interact strongly with changing climate and trace gas abundances. Significant changes in ozone have been observed and prediction of future ozone and climate interactions depends on modeling these processes successfully. We have conducted an off-line model simulation of the stratosphere for trace gas conditions representative of 1975-2025 using meteorology from the NASA finite-volume general circulation model. The objective of this simulation is to examine the sensitivity of stratospheric ozone and chemical change to varying meteorology and trace gas inputs. This presentation will examine the dependence of ozone and related processes in polar regions on the climatological and trace gas changes in the model. The model past performance is base-lined against available observations, and a future ozone recovery scenario is forecast. Overall the model ozone simulation is quite realistic, but initial analysis of the detailed evolution of some observable processes suggests systematic shortcomings in our description of the polar chemical rates and/or mechanisms. Model sensitivities, strengths, and weaknesses will be discussed with implications for uncertainty and confidence in coupled climate chemistry predictions.

Kawa, S.R.↗

NASA NTRS2004

Ozone Measurements and a 3D Chemical Transport Model

We have used our three-dimensional chemical transport model (CTM) to calculate the expected reponse of stratospheric composition over the past 30 years to forcing by chlorine and bromine compounds, solar ultraviolet, and volcanic aerosols. The CTM uses off-line winds and temperatures fiom a 50-year run of the finite volume general circulation model (FVGCM). We compare the total column ozone and the ozone profile fiom the CTM output to a variety of data sources. These include a merged total ozone data set from TOMS and SBUV using the new version 8 algorithm. Total ozone fiom the CTM are compared to ground-station measurements of total ozone at specific locations. Ozone profiles are compared to satellite meausrements fiom SBUV, SAGE, and HALOE. Profiles are also compared to ozonesondes over several locations. The results of the comparisons are quantified by using a time-series statistical analysis to determine trends, solar cycle, and volcanic reponse in both the model and in the data. Initial results indicate that the model responds to forcings in a way that is similar to the observed atmospheric response. The model does seem to be more sensitive to the chlorine and bromine perturbation ihan is the data. Further details and comparisons wiii be discussed.

Stolarski, Richard S.↗

NASA NTRS1999

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

NASA NTRS2002

Heat Transfer in High-Temperature Fibrous Insulation

The combined radiation/conduction heat transfer in high-porosity, high-temperature fibrous insulations was investigated experimentally and numerically. The effective thermal conductivity of fibrous insulation samples was measured over the temperature range of 300-1300 K and environmental pressure range of 1.33 x 10(exp -5)-101.32 kPa. The fibrous insulation samples tested had nominal densities of 24, 48, and 72 kilograms per cubic meter and thicknesses of 13.3, 26.6 and 39.9 millimeters. Seven samples were tested such that the applied heat flux vector was aligned with local gravity vector to eliminate natural convection as a mode of heat transfer. Two samples were tested with reverse orientation to investigate natural convection effects. It was determined that for the fibrous insulation densities and thicknesses investigated no heat transfer takes place through natural convection. A finite volume numerical model was developed to solve the governing combined radiation and conduction heat transfer equations. Various methods of modeling the gas/solid conduction interaction in fibrous insulations were investigated. The radiation heat transfer was modeled using the modified two-flux approximation assuming anisotropic scattering and gray medium. A genetic-algorithm based parameter estimation technique was utilized with this model to determine the relevant radiative properties of the fibrous insulation over the temperature range of 300-1300 K. The parameter estimation was performed by least square minimization of the difference between measured and predicted values of effective thermal conductivity at a density of 24 kilograms per cubic meters and at nominal pressures of 1.33 x 10(exp -4) and 99.98 kPa. The numerical model was validated by comparison with steady-state effective thermal conductivity measurements at other densities and pressures. The numerical model was also validated by comparison with a transient thermal test simulating reentry aerodynamic heating conditions.

Daryabeigi, Kamran↗

NASA NTRS2004

Transient Two-Dimensional Analysis of Side Load in Liquid Rocket Engine Nozzles

Two-dimensional planar and axisymmetric numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to develop a computational methodology to identify nozzle side load physics using simplified two-dimensional geometries, in order to come up with a computational strategy to eventually predict the three-dimensional side loads. The computational methodology is based on a multidimensional, finite-volume, viscous, chemically reacting, unstructured-grid, and pressure-based computational fluid dynamics formulation, and a transient inlet condition based on an engine system modeling. The side load physics captured in the low aspect-ratio, two-dimensional planar nozzle include the Coanda effect, afterburning wave, and the associated lip free-shock oscillation. Results of parametric studies indicate that equivalence ratio, combustion and ramp rate affect the side load physics. The side load physics inferred in the high aspect-ratio, axisymmetric nozzle study include the afterburning wave; transition from free-shock to restricted-shock separation, reverting back to free-shock separation, and transforming to restricted-shock separation again; and lip restricted-shock oscillation. The Mach disk loci and wall pressure history studies reconfirm that combustion and the associated thermodynamic properties affect the formation and duration of the asymmetric flow.

Wang, Ten-See↗

NASA NTRS2004

Multidimensional Unstructured-Grid Liquid Rocket Engine Nozzle Performance and Heat Transfer Analysis

The objectives of this study are to conduct a unified computational analysis for computing the design parameters such as axial thrust, convective and radiative wall heat fluxes for liquid rocket engine nozzles, so as to develop a computational strategy for computing those design parameters through parametric investigations. The computational methodology is based on a multidimensional, finite-volume, turbulent, chemically reacting, radiating, unstructured-grid, and pressure-based formulation, with grid refinement capabilities. Systematic parametric studies on effects of wail boundary conditions, combustion chemistry, radiation coupling, computational cell shape, and grid refinement were performed and assessed. Comparisons of the computed axial thrust performance, flow features, and wail heat fluxes with those of available test data and design calculations are presented.

Wang, Ten-See↗

NASA NTRS1999

A Pseudo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

White, J. A.↗

NASA NTRS1999

X-33 Aerodynamic and Aeroheating Computations for Wind Tunnel and Flight Conditions

This report provides an overview of hypersonic Computational Fluid Dynamics research conducted at the NASA Langley Research Center to support the Phase II development of the X-33 vehicle. The X-33, which is being developed by Lockheed-Martin in partnership with NASA, is an experimental Single-Stage-to-Orbit demonstrator that is intended to validate critical technologies for a full-scale Reusable Launch Vehicle. As part of the development of the X-33, CFD codes have been used to predict the aerodynamic and aeroheating characteristics of the vehicle. Laminar and turbulent predictions were generated for the X 33 vehicle using two finite- volume, Navier-Stokes solvers. Inviscid solutions were also generated with an Euler code. Computations were performed for Mach numbers of 4.0 to 10.0 at angles-of-attack from 10 deg to 48 deg with body flap deflections of 0, 10 and 20 deg. Comparisons between predictions and wind tunnel aerodynamic and aeroheating data are presented in this paper. Aeroheating and aerodynamic predictions for flight conditions are also presented.

Hollis, Brian R.↗

NASA NTRS2003

Mathematical model of diffusion-limited evolution of multiple gas bubbles in tissue

Models of gas bubble dynamics employed in probabilistic analyses of decompression sickness incidence in man must be theoretically consistent and simple, if they are to yield useful results without requiring excessive computations. They are generally formulated in terms of ordinary differential equations that describe diffusion-limited gas exchange between a gas bubble and the extravascular tissue surrounding it. In our previous model (Ann. Biomed. Eng. 30: 232-246, 2002), we showed that with appropriate representation of sink pressures to account for gas loss or gain due to heterogeneous blood perfusion in the unstirred diffusion region around the bubble, diffusion-limited bubble growth in a tissue of finite volume can be simulated without postulating a boundary layer across which gas flux is discontinuous. However, interactions between two or more bubbles caused by competition for available gas cannot be considered in this model, because the diffusion region has a fixed volume with zero gas flux at its outer boundary. The present work extends the previous model to accommodate interactions among multiple bubbles by allowing the diffusion region volume of each bubble to vary during bubble evolution. For given decompression and tissue volume, bubble growth is sustained only if the bubble number density is below a certain maximum.

Evaluation Studies↗

NASA NTRS1997

Hypersonic Flow Computations on Unstructured Meshes

A method for computing inviscid hypersonic flow over complex configurations using unstructured meshes is presented. The unstructured grid solver uses an edge{based finite{volume formulation. Fluxes are computed using a flux vector splitting scheme that is capable of representing constant enthalpy solutions. Second{order accuracy in smooth flow regions is obtained by linearly reconstructing the solution, and stability near discontinuities is maintained by locally forcing the scheme to reduce to first-order accuracy. The implementation of the algorithm to parallel computers is described. Computations using the proposed method are presented for a sphere-cone configuration at Mach numbers of 5.25 and 10.6, and a complex hypersonic re-entry vehicle at Mach numbers of 4.5 and 9.8. Results are compared to experimental data and computations made with established structured grid methods. The use of the solver as a screening tool for rapid aerodynamic assessment of proposed vehicles is described.

Bibb, K. L.↗

NASA NTRS2004

Testing the PV-Theta Mapping Technique in a 3-D CTM Model Simulation

Mapping lower stratospheric ozone into potential vorticity (PV)- potential temperature (Theta) coordinates is a common technique employed to analyze sparse data sets. Ozone transformed into a flow-following dynamical coordinate system is insensitive to meteorological variations. Therefore data from a wide range of times/locations can be compared, so long as the measurements were made in the same airmass (as defined by PV). Moreover, once a relationship between ozone and PV/Theta is established, a full 3D ozone field can be estimated from this relationship and the 3D analyzed PV field. However, ozone data mapped in this fashion can be hampered by noisy PV fields, or "mis-matches" in the resolution and/or exact location of the ozone and PV measurements. In this study, we investigate the PV-ozone relationship using output from a recent 50-year run of the Goddard 3D chemical transport model (CTM). Model constituents are transported using off-line dynamics from the finite volume general circulation model (FVGCM). By using the internally consistent model PV and ozone fields, we minimize noise due to mis-matching and resolution issues. We calculate correlations between model ozone and PV throughout the stratosphere, and test the sensitivity of the technique to initial data resolution. To do this we degrade the model data to that of various satellite instruments, then compare the mapped fields derived from the sub-sampled data to the full resolution model data. With these studies we can determine appropriate limits for the PV-theta mapping technique in latitude, altitude, and as a function of original data resolution.

Frith, Stacey M.↗

NASA NTRS1997

Multi-Body Analysis of a Tiltrotor Configuration

The paper describes the aeroelastic analysis of a tiltrotor configuration. The 1/5 scale wind tunnel semispan model of the V-22 tiltrotor aircraft is considered. The analysis is performed by means of a multi-body code, based on an original formulation. The differential equilibrium problem is stated in terms of first order differential equations. The equilibrium equations of every rigid body are written, together with the definitions of the momenta. The bodies are connected by kinematic constraints, applied in form of Lagrangian multipliers. Deformable components are mainly modelled by means of beam elements, based on an original finite volume formulation. Multi-disciplinar problems can be solved by adding user-defined differential equations. In the presented analysis the equations related to the control of the swash-plate of the model are considered. Advantages of a multi-body aeroelastic code over existing comprehensive rotorcraft codes include the exact modelling of the kinematics of the hub, the detailed modelling of the flexibility of critical hub components, and the possibility to simulate steady flight conditions as well as wind-up and maneuvers. The simulations described in the paper include: 1) the analysis of the aeroelastic stability, with particular regard to the proprotor/pylon instability that is peculiar to tiltrotors, 2) the determination of the dynamic behavior of the system and of the loads due to typical maneuvers, with particular regard to the conversion from helicopter to airplane mode, and 3) the stress evaluation in critical components, such as the pitch links and the conversion downstop spring.

Ghiringhelli, G. L.↗

NASA NTRS1998

On Some Numerical Dissipation Schemes

Several schemes for introducing an artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar dissipation and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical solutions are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme. The coarse-grid accuracy for the original CUSP scheme is improved by modifying the limiter function used with the scheme, giving comparable accuracy to that obtained with the MATD scheme. The modifications reduce the background dissipation and provide control over the regions where the scheme can become first order.

Swanson, R. C.↗

NASA NTRS2004

Computation of Tone Noises Generated in Viscous Flows

Three benchmark problems from the current and previous CAA workshops involving tone noise generated in viscous flows are investigated using the CE/SE finite volume method. The CE/SE method is first briefly reviewed. Then, the benchmark problems, namely, flow past a single cylinder (CAA Workshop II problem), flow past twin cylinders (from the current CAA Workshop IV, Category 5, Problem 1) and flow past a deep cavity with overhang (CAA Workshop III problem) are investigated. Generally good results are obtained in comparison with the experimental data.

Loh, Ching Y.↗

NASA NTRS2004

A Fast Upwind Solver for the Euler Equations on Three-Dimensional Unstructured Meshes

An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a novel cell reconstruction process which results in computational times per cell comparable to those of structured codes. The approach yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicit limiting. Solutions are advanced in time by a 3-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Solutions are presented for a range of configurations in the transonic speed regime to demonstrate code accuracy, speed, and robustness. The results include an assessment of grid sensitivity and convergence acceleration by mesh sequencing.

Frink, Neal T.↗

NASA NTRS2004

Stratospheric Impact of Varying Sea Surface Temperatures

The Finite-Volume General Circulation Model (FVGCM) has been run in 50 year simulations with the: 1) 1949-1999 Hadley Centre sea surface temperatures (SST), and 2) a fixed annual cycle of SSTs. In this presentation we first show that the 1949-1999 FVGCM simulation produces a very credible stratosphere in comparison to an NCEP/NCAR reanalysis climatology. In particular, the northern hemisphere has numerous major and minor stratospheric warming, while the southern hemisphere has only a few over the 50-year simulation. During the northern hemisphere winter, temperatures are both warmer in the lower stratosphere and the polar vortex is weaker than is found in the mid-winter southern hemisphere. Mean temperature differences in the lower stratosphere are shown to be small (less than 2 K), and planetary wave forcing is found to be very consistent with the climatology. We then will show the differences between our varying SST simulation and the fixed SST simulation in both the dynamics and in two parameterized trace gases (ozone and methane). In general, differences are found to be small, with subtle changes in planetary wave forcing that lead to reduced temperatures in the SH and increased temperatures in the NH.

Newman, Paul A.↗

NASA NTRS2004

Sensitivity of the FVGCM to Changes in Ozone

We have carried out an experiment with the finite volume general circulation model (FVGCM). This experiment consisted of two different imposed changes in the climatological ozone fields assumed in the radiation code. for conditions with no significant ozone hole. This distribution was obtained from a 50-year simulation of the full stratospheric ozone chemistry, with a time-dependent chlorine loading, done with our off-line chemical transport model (CTM). Three years (1978-1980) of this simulation were averaged to form a monthly, zonal-mean ozone distribution that was used in the 20-year integration of the FVGCM for "unperturbed" conditions. The second 20-year GCM integration included a fully-developed ozone hole. This ozone distribution was from three years, 1998-2000, from the same CTM simulation. The goal of this work is to determine the coupled response of the chemistry and dynamics of the stratosphere. These experiments are the first step in understanding the coupled response. An important initial question concerns the significance of the signals: if 20-year integrations turn out to be too short, the runs will be extended.

Stolarski↗