Engineering Papers⌕ Search

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.

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

At least 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Numerical Modeling of Pressurization of a Propellant Tank

An unsteady finite volume procedure has been developed to predict the history o pressure, temperature and mass flow rate of the pressurant and propellant during the expulsion of the propellant from a tan. The time dependent mass, momentum and energy conservation equations are solved at the ullage space. The model accounts for the change in the ullage volume due to expulsion of the propellant. It also accounts for the heat transfer from the tank wall and propellant to the ullage gas. The procedure was incorporated in the Generalized Fluid System Simulation Program (GFSSP). The results of several test cases were then compared with a published correlation of pressurant requirements for a given displacement of propellant. The agreement between the predictions and the correlation was found to be satisfactory.

Majumdar, Alok↗
NASA NTRS1999

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗
NASA NTRS2000

Three-Dimensional Incompressible Navier-Stokes Flow Computations about Complete Configurations Using a Multiblock Unstructured Grid Approach

A multiblock unstructured grid approach is presented for solving three-dimensional incompressible inviscid and viscous turbulent flows about complete configurations. The artificial compressibility form of the governing equations is solved by a node-based, finite volume implicit scheme which uses a backward Euler time discretization. Point Gauss-Seidel relaxations are used to solve the linear system of equations at each time step. This work employs a multiblock strategy to the solution procedure, which greatly improves the efficiency of the algorithm by significantly reducing the memory requirements by a factor of 5 over the single-grid algorithm while maintaining a similar convergence behavior. The numerical accuracy of solutions is assessed by comparing with the experimental data for a submarine with stem appendages and a high-lift configuration.

Sheng, Chunhua↗
NASA NTRS1999

Adaptive Mesh Refinement for Microelectronic Device Design

Finite element and finite volume methods are used in a variety of design simulations when it is necessary to compute fields throughout regions that contain varying materials or geometry. Convergence of the simulation can be assessed by uniformly increasing the mesh density until an observable quantity stabilizes. Depending on the electrical size of the problem, uniform refinement of the mesh may be computationally infeasible due to memory limitations. Similarly, depending on the geometric complexity of the object being modeled, uniform refinement can be inefficient since regions that do not need refinement add to the computational expense. In either case, convergence to the correct (measured) solution is not guaranteed. Adaptive mesh refinement methods attempt to selectively refine the region of the mesh that is estimated to contain proportionally higher solution errors. The refinement may be obtained by decreasing the element size (h-refinement), by increasing the order of the element (p-refinement) or by a combination of the two (h-p refinement). A successful adaptive strategy refines the mesh to produce an accurate solution measured against the correct fields without undue computational expense. This is accomplished by the use of a) reliable a posteriori error estimates, b) hierarchal elements, and c) automatic adaptive mesh generation. Adaptive methods are also useful when problems with multi-scale field variations are encountered. These occur in active electronic devices that have thin doped layers and also when mixed physics is used in the calculation. The mesh needs to be fine at and near the thin layer to capture rapid field or charge variations, but can coarsen away from these layers where field variations smoothen and charge densities are uniform. This poster will present an adaptive mesh refinement package that runs on parallel computers and is applied to specific microelectronic device simulations. Passive sensors that operate in the infrared portion of the spectrum as well as active device simulations that model charge transport and Maxwell's equations will be presented.

Cwik, Tom↗
NASA NTRS2000

The NASA Tetrahedral Unstructured Software System (TETRUSS)

The NASA Tetrahedral Unstructured Software System (TetrUSS) was developed during the 1990's to provide a rapid aerodynamic analysis and design capability to applied aerodynamicists. The system is comprised of loosely integrated, user-friendly software that enables the application of advanced Euler and Navier-Stokes tetrahedral finite volume technology to complex aerodynamic problems. TetrUSS has matured well because of the generous feedback from many willing users representing a broad cross-section of background and skill levels. This paper presents an overview of the current capabilities of the TetrUSS system along with some representative results from selected applications.

Frink, Neal T.↗
NASA NTRS2000

Model Study of the Sources of Upper Tropospheric Ozone Over The Atlantic, with Emphasis on Interhemispheric Asymmetries

Measurements of ozone and aerosols from ships transecting the Atlantic Ocean reveal considerable interhemispheric asymmetries in the vertical structure of the ozone distribution. There is much more upper tropospheric ozone to the South of the Equator than in the northern hemisphere. These profiles suggest that the air in the southern hemispheric tropics and subtropics has a source in the stratosphere, while that in the northern hemisphere originates in the boundary layer. A model experiment designed to investigate this problem will be discussed. They are based on a new general circulation model, based on a Finite-Volume dynamical core and the physical parametrizations from the Community Climate Model at the National Center for Atmospheric Research (NCAR). The model simulates a climatological state in reasonable agreement with observations, comparable to other contemporary climate models, but has a much more accurate transport capability than most other present-day models. The experiment involves using a suite of tracers, designed to highlight air from different source regions. Some of these (water vapor and ozone) represent real atmospheric quantities, while others are markers which indicate whether the air originates in the stratosphere, the boundary layer, or has traversed regions of convection (for example). The experiments were run for periods of up to one month, initialized from snapshots from a long climate integration of the model. A sample of results will be presented, chosen to illustrate the different sources of air in the upper troposphere and the physical mechanisms by which these air parcels arrived there.

Pawson, Steven↗
NASA NTRS2000

Assimilation of (A) TOVS data at the NASA Goddard Data Assimilation Office

At the NASA Goddard Data Assimilation Office (DAO), a 1D variational radiance assimilation system has been developed. This system, called DAOTOVS (Tiros operational vertical sounder (TOVS)), uses (A)TOVS level 1b radiances. It has been implemented within the DAO's semi-operational system as well as within the next generation data assimilation system that uses a finite-volume dynamical core. We will show results from (A)TOVS assimilation, including stratospheric analyses and validation. We will also describe our systematic error correction scheme which is based on collocated radiosondes.

Joiner, Joanna↗
NASA NTRS2000

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗
NASA NTRS2000

Using TRMM Field Campaign Data for Assessing GEOS Forecast and Assimilation Products

The wealth of in-situ measurements gathered during Tropical Rain Measuring Mission (TRMM) field campaigns over a wide range of tropical conditions constitute an important data source for evaluating the quality of global model forecasts and assimilated datasets. In this study we use selected observations of cloud microphysics and atmospheric sounding from TEFLUN-1998, SCEMEX-1998, and TRMMLBA-1999 to examine the assimilation and forecast fields produced by the operational GEOS-3 (Goddard Earth Observing System - version 3) global data assimilation system (DAS) and a new finite-volume DAS under development at the Data Assimilation Office. Additionally, TRMM field campaign measurements are used to verify the impact of assimilating rainfall and moisture data derived from TRMM Microwave Imager and Special Sensor Microwave/Imager instruments on the GEOS analysis. We will also explore issues concerning the 'error of representativeness' in using in-situ observations of quantities with large spatial and temporal variability such as precipitation for validating gridded global data products.

Hou, Arthur Y.↗
NASA NTRS2000

A Technique of Treating Negative Weights in WENO Schemes

High order accurate weighted essentially non-oscillatory (WENO) schemes have recently been developed for finite difference and finite volume methods both in structural and in unstructured meshes. A key idea in WENO scheme is a linear combination of lower order fluxes or reconstructions to obtain a high order approximation. The combination coefficients, also called linear weights, are determined by local geometry of the mesh and order of accuracy and may become negative. WENO procedures cannot be applied directly to obtain a stable scheme if negative linear weights are present. Previous strategy for handling this difficulty is by either regrouping of stencils or reducing the order of accuracy to get rid of the negative linear weights. In this paper we present a simple and effective technique for handling negative linear weights without a need to get rid of them.

Shi, Jing↗
NASA NTRS2001

Regional Climate Studies with Variable-Resolution Stretched-Grid GCMs

A variable resolution GCM using a global stretched grid with fine resolution over the area(s) of interest, is a viable new approach to regional and subregional climate studies and applications. It is an alternative to the widely used nested grid approach introduced a decade ago as a pioneering step in regional climate modeling. The first version of the SG-GCM based on the GEOS (Goddard Earth Observing System) GCM using a finite-difference approximation, has been developed and thoroughly tested during the last few years. Successful simulations have been performed with the SG-GCM for the anomalous regional climate events of the U.S. 1988 summer drought and 1993 summer flood. They have shown the practical feasibility of the SG-approach for regional climate modeling. The GEOS SG-DAS (Data Assimilation System) incorporating the SG-GCM has also been developed and tested. The assimilated regional fields and diagnostics are used for validating the SG-GCM regional simulations. Two new SG-GCMs are being developed. The first is the SG-version of the new NASANCAR FV-GCM (with the finite-volume (FV) dynamics), and the second is the SG-version of the new GCM with spectral-element dynamics. Both GCMs use the WAR CCM4 physics. Using these advanced numerics will provide increased computational efficiency for the new the SG-GCMs, and will allow us to employ more flexible stretching strategies beneficial for the efficient regional down-scaling. The major current developments are focused on: simulating the 1997-1999 (and beyond) ENSO cycle and related monsoonal circulations, with enhanced regional resolution; studying intraseasonal and interannual regional climate variability for the extended multiyear (AMIP-type) SG-GCM simulations; and studying the impact of ensemble integrations.

Fox-Rabinovitz, Michael↗
NASA NTRS2000

High Order Discontinuous Gelerkin Methods for Convection Dominated Problems with Application to Aeroacoustics

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. On the analysis side, we have studied the efficient and stable discontinuous Galerkin framework for small second derivative terms, for example in Navier-Stokes equations, and also for related equations such as the Hamilton-Jacobi equations. This is a truly local discontinuous formulation where derivatives are considered as new variables. On the applied side, we have implemented and tested the efficiency of different approaches numerically. Related issues in high order ENO and WENO finite difference methods and spectral methods have also been investigated. Jointly with Hu, we have presented a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the RungeKutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method. Jointly with Hu, we have constructed third and fourth order WENO schemes on two dimensional unstructured meshes (triangles) in the finite volume formulation. The third order schemes are based on a combination of linear polynomials with nonlinear weights, and the fourth order schemes are based on combination of quadratic polynomials with nonlinear weights. We have addressed several difficult issues associated with high order WENO schemes on unstructured mesh, including the choice of linear and nonlinear weights, what to do with negative weights, etc. Numerical examples are shown to demonstrate the accuracies and robustness of the methods for shock calculations. Jointly with P. Montarnal, we have used a recently developed energy relaxation theory by Coquel and Perthame and high order weighted essentially non-oscillatory (WENO) schemes to simulate the Euler equations of real gas. The main idea is an energy decomposition under the form epsilon = epsilon(sub 1) + epsilon(sub 2), where epsilon(sub 1) is associated with a simpler pressure law (gamma)-law in this paper) and the nonlinear deviation epsilon(sub 2) is convected with the flow. A relaxation process is performed for each time step to ensure that the original pressure law is satisfied. The necessary characteristic decomposition for the high order WENO schemes is performed on the characteristic fields based on the epsilon(sub l) gamma-law. The algorithm only calls for the original pressure law once per grid point per time step, without the need to compute its derivatives or any Riemann solvers. Both one and two dimensional numerical examples are shown to illustrate the effectiveness of this approach.

Shu, Chi-Wang↗
NASA NTRS2001

Application of Aeroelastic Solvers Based on Navier Stokes Equations

The propulsion element of the NASA Advanced Subsonic Technology (AST) initiative is directed towards increasing the overall efficiency of current aircraft engines. This effort requires an increase in the efficiency of various components, such as fans, compressors, turbines etc. Improvement in engine efficiency can be accomplished through the use of lighter materials, larger diameter fans and/or higher-pressure ratio compressors. However, each of these has the potential to result in aeroelastic problems such as flutter or forced response. To address the aeroelastic problems, the Structural Dynamics Branch of NASA Glenn has been involved in the development of numerical capabilities for analyzing the aeroelastic stability characteristics and forced response of wide chord fans, multi-stage compressors and turbines. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading is available. To obtain the steady and unsteady aerodynamic forces for the complex flows around the engine components, for the flow regimes encountered by the rotor, an advanced compressible Navier-Stokes solver is required. A finite volume based Navier-Stokes solver has been developed at Mississippi State University (MSU) for solving the flow field around multistage rotors. The focus of the current research effort, under NASA Cooperative Agreement NCC3- 596 was on developing an aeroelastic analysis code (entitled TURBO-AE) based on the Navier-Stokes solver developed by MSU. The TURBO-AE code has been developed for flutter analysis of turbomachine components and delivered to NASA and its industry partners. The code has been verified. validated and is being applied by NASA Glenn and by aircraft engine manufacturers to analyze the aeroelastic stability characteristics of modem fans, compressors and turbines.

Keith, Theo G., Jr.↗
NASA NTRS1994

Hyperbolic Prismatic Grid Generation and Solution of Euler Equations on Prismatic Grids

A hyperbolic grid generation method is used to generate prismatic grids and an approach using prismatic grids to solve the Euler equations is presented. The theory of the stability and feasibility of the hyperbolic grid generation method is presented. The hyperbolic grid generation method of Steger et al for structured grids is applied to a three dimensional triangularized surface definition to generate a grid that is unstructured on each successive layer. The grid, however, retains structure in the body-normal direction and has a computational cell shaped like a triangular prism. In order to take advantage of the structure in the normal direction, a finite-volume scheme that treats the unknowns along the normal direction implicitly is introduced and the flow over a sphere is simulated.

Pandya, S. A.↗
NASA NTRS2001

Linearized Unsteady Aerodynamic Analysis of the Acoustic Response to Wake/Blade-Row Interaction

The three-dimensional, linearized Euler analysis, LINFLUX, is being developed to provide a comprehensive and efficient unsteady aerodynamic scheme for predicting the aeroacoustic and aeroelastic responses of axial-flow turbomachinery blading. LINFLUX couples a near-field, implicit, wave-split, finite-volume solution to far-field acoustic eigensolutions, to predict the aerodynamic responses of a blade row to prescribed structural and aerodynamic excitations. It is applied herein to predict the acoustic responses of a fan exit guide vane (FEGV) to rotor wake excitations. The intent is to demonstrate and assess the LINFLUX analysis via application to realistic wake/blade-row interactions. Numerical results are given for the unsteady pressure responses of the FEGV, including the modal pressure responses at inlet and exit. In addition, predictions for the modal and total acoustic power levels at the FEGV exit are compared with measurements. The present results indicate that the LINFLUX analysis should be useful in the aeroacoustic design process, and for understanding the three-dimensional flow physics relevant to blade-row noise generation and propagation.

Verdon, Joseph M.↗
NASA NTRS2000

Numerical Solutions to the Third CAA Workshop Benchmark Problems

This paper presents numerical solutions to the problems of propagation of sound waves through a transonic nozzle, shock-sound interactions, and automobile noise involving feedback of the third NASA Computational Aeroacoustics (CAA) Workshop on benchmark problems. The numerical algorithm is based on a dual time scheme for temporal discretization and a third-order finite volume scheme for spatial discretization. The aims of this study are to apply a dual time stepping scheme to treat aeroacoustic problems of sound propagation and to validate our CAA solver with the benchmark problems for developing a numerical tool for noise analysis and control.

Loh, Roy H.↗
NASA NTRS2001

Heat Transfer in Adhesively Bonded Honeycomb Core Panels

The Swann and Pittman semi-empirical relationship has been used as a standard in aerospace industry to predict the effective thermal conductivity of honeycomb core panels. Recent measurements of the effective thermal conductivity of an adhesively bonded titanium honeycomb core panel using three different techniques, two steady-state and one transient radiant step heating method, at four laboratories varied significantly from each other and from the Swann and Pittman predictions. Average differences between the measurements and the predictions varied between 17 and 61% in the temperature range of 300 to 500 K. In order to determine the correct values of the effective thermal conductivity and determine which set of the measurements or predictions were most accurate, the combined radiation and conduction heat transfer in the honeycomb core panel was modeled using a finite volume numerical formulation. The transient radiant step heating measurements provided the best agreement with the numerical results. It was found that a modification of the Swann and Pittman semi-empirical relationship which incorporated the facesheets and adhesive layers in the thermal model provided satisfactory results. Finally, a parametric study was conducted to investigate the influence of adhesive thickness and thermal conductivity on the overall heat transfer through the panel.

Daryabeigi, Kamran↗
NASA NTRS2001

Thermal Analysis and Design of Multi-layer Insulation for Re-entry Aerodynamic Heating

The combined radiation/conduction heat transfer in high-temperature multi-layer insulations was modeled using a finite volume numerical model. The numerical model was validated by comparison with steady-state effective thermal conductivity measurements, and by transient thermal tests simulating re-entry aerodynamic heating conditions. A design of experiments technique was used to investigate optimum design of multi-layer insulations for re-entry aerodynamic heating. It was found that use of 2 mm foil spacing and locating the foils near the hot boundary with the top foil 2 mm away from the hot boundary resulted in the most effective insulation design. A 76.2 mm thick multi-layer insulation using 1, 4, or 16 foils resulted in 2.9, 7.2, or 22.2 percent mass per unit area savings compared to a fibrous insulation sample at the same thickness, respectively.

Daryabeigi, Kamran↗