Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure model”

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 631 records · Page 35

On the prediction of turbulent secondary flows

The prediction of turbulent secondary flows, with Reynolds stress models, in circular pipes and non-circular ducts is reviewed. Turbulence-driven secondary flows in straight non-circular ducts are considered along with turbulent secondary flows in pipes and ducts that arise from curvature or a system rotation. The physical mechanisms that generate these different kinds of secondary flows are outlined and the level of turbulence closure required to properly compute each type is discussed in detail. Illustrative computations of a variety of different secondary flows obtained from two-equation turbulence models and second-order closures are provided to amplify these points.

Speziale, C. G.↗

Shipboard Sunphotometer Measurements of Aerosol Optical Depth Spectra and Columnar Water Vapor During ACE-2 and Comparison with Selected Land, Ship, Aircraft, and Satellite Measurements

Analyses of aerosol optical depth (AOD) and colurnmn water vapor (CWV) measurements acquired with NASA Ames Research Center's 6-channel Airborne Tracking Sunphotometer (AATS-6) operated aboard the R/V Professor Vodyanitskiy during the 2nd Aerosol Characterization Experiment (ACE-2) are discussed. Data are compared with various in situ and remote measurements for selected cases. The focus is on 10 July, when the Pelican airplane flew within 70 km of the ship near the time of a NOAA-14/AVHRR satellite overpass and AOD measurements with the 14-channel Ames Airborne Tracking Sunphotometer (AATS-14) above the marine boundary layer (MBL) permitted calculation of AOD within the MBL from the AATS-6 measurements. A detailed column closure test is performed for MBL AOD on 10 July by comparing the AATS-6 MBL AODs with corresponding values calculated by combining shipboard particle size distribution measurements with models of hygroscopic growth and radiosonde humidity profiles (plus assumptions on the vertical profile of the dry particle size distribution and composition). Large differences (30-80% in the mid-visible) between measured and reconstructed AODs are obtained, in large part because of the high sensitivity of the closure methodology to hygroscopic growth models, which vary considerably and have not been validated over the necessary range of particle size/composition distributions. The wavelength dependence of AATS-6 AODs is compared with the corresponding dependence of aerosol extinction calculated from shipboard measurements of aerosol size distribution and of total scattering mearured by a shipboard integrating nephelometer for several days. Results are highly variable, illustrating further the great difficulty of deriving column values from point measurements. AATS-6 CWV values are shown to agree well with corresponding values derived from radiosonde measurements during 8 soundings on 7 days and also with values calculated from measurements taken on 10 July with the AATS-14 and the University of Washington Passive Humidigraph aboard the Pelican.

Livingston, John M.↗

Shipboard Sunphotometer Measurements of Aerosol Optical Depth Spectra and Columnar Water Vapor During ACE-2, and Comparison with Selected Land, Ship, Aircraft, and Satellite Measurements

Analyses of aerosol optical depth (AOD) and columnar water vapor (CWV) measurements acquired with NASA Ames Research Center's six-channel Airborne Tracking Sunphotometer (AATS-6) operated aboard the R/V (research vehicle) Professor Vodyanitskiy during the second Aerosol Characterization Experiment (ACE-2) are discussed. Data are compared with various in situ and remote measurements for selected cases. The focus is on 10 July, when the Pelican airplane flew within 70 km of the ship near the time of a NOAA (National Oceanographic and Atmospheric Administration)-14/AVHRR (Advanced Very High Resolution Radiometer) satellite overpass and AOD measurements with the 14-channel Ames Airborne Tracking Sunphotometer (AATS-14) above the marine boundary layer (MBL) permitted calculation of AOD within the MBL from the AATS-6 measurements. A detailed column closure test is performed for MBL AOD on 10 July by comparing the AATS-6 MBL AODs with corresponding values calculated by combining shipboard particle size distribution measurements with models of hygroscopic growth and radiosonde humidity profiles (plus assumptions on the vertical profile of the dry particle size distribution and composition). Large differences (30-80% in the mid-visible) between measured and reconstructed AODs are obtained, in large part because of the high sensitivity of the closure methodology to hygroscopic growth models, which vary considerably and have not been validated over the necessary range of particle size/composition distributions. The wavelength dependence of AATS-6 AODs is compared with the corresponding dependence of aerosol extinction calculated from shipboard measurements of aerosol size distribution and of total scattering measured by a shipboard integrating nephelometer for several days. Results are highly variable, illustrating further the great difficulty of deriving column values from point measurements. AATS-6 CWV values are shown to agree well with corresponding values derived from radiosonde measurements during eight soundings on seven days and also with values calculated from measurements taken on 10 July with the AATS-14 and the University of Washington Passive Humidigraph aboard the Pelican.

Livingston, John M.↗

Scale resolving simulations of the NASA Juncture Flow Model using the LAVA solver

Two distinct simulation methodologies: Delayed Detached Eddy Simulation (DDES) and stress based Wall-modelled Large Eddy Simulation (WMLES) are evaluated using structured overset curvilinear grids for the NASA juncture flow model. While both methodologies are shown to mitigate the primary shortcomings of steady state Reynolds Averaged Navier Stokes(RANS) simulations, several unresolved aspects are identified. Strong sensitivity to RANS-typegrid refinement is observed in the DDES with a substantial deterioration of the solution quality with increasing spatial resolution associated with deficiencies in the shielding function. Mean profiles for attached boundary layers on the fuselage show spurious inflections suggesting modelled stress depletion on finer grids. Lower numerical dissipation in terms of spatial discretization and time step size is seen to improve the solution quality, and the advantage of using RC and QCR2000 corrections in the underlying RANS closure is demonstrated for DDES on RANS-type meshes. Equilibrium Wall-Modelled LES used grids that resolved the tripping dots over the fuselage nose and the wing leading edge, consistent with the experiment setup. These simulations resulted in a cost-competitive approach compared to DDES on RANS-type grids. Although the agreement between WMLES predictions of first and second order single point statistics with experimental measurements is promising, a fundamental shortcoming is noted in terms of an overshoot in streamwise momentum in corner regions of the wing-fuselage juncture. This excess momentum subsequently delays the onset of separation,thereby resulting in an underprediction in length of the separation bubble. Some quantitative sensitivity to numerical discretization is observed; lowering of numerical dissipation shows better agreement with the experiment.

T^3↗

ECLS for the NASA CDG Space Station model

This paper describes the approaches to Environmental Control and Life Support (ECLS) loop closure that have been studied for a permanently manned Space Station, by the NASA Concept Development Group (CDG) and industry. The paper explains the rationale behind the level of the loop closure chosen for ECLS expendables and why this level was chosen. The paper also explains potential synergistic interactions, between ECLS subsystems and other vehicle systems, which can lessen vehicle expendable penalties and help simplify required ECLS hardware, interfaces and control. Of particular significance are the potential interfacing of the ECLS subsystem which concentrates the metabolic CO2 removed from the atmosphere, with propulsion system resistojets and the use of cryogenic boiloff to satisfy crew metabolic needs. The paper shows the trade-offs that must be performed in order to make a selection from available options.

Cushman, R. J.↗

Numerical solution of three-dimensional free turbulent shear flows

Numerical results are presented for three-dimensional compressible turbulent jet and wake flows. An alternating direction implicit numerical procedure is used to solve the finite-difference form of the parabolic elliptic Navier-Stokes equations. A coordinate transformation maps the boundaries at infinity into a finite computational domain in order to properly specify infinity boundary conditions as well as contain the downstream growth of the viscous flow field in a fixed computational grid. Turbulence closure is achieved through an algebraic mixing length eddy viscosity model. Numerical results for supersonic flow are presented for an axisymmetric jet, an elliptical jet, an elliptical wake, and two interacting rectangular jets. Experimental data were not available for comparison with the numerical results. However, the results compare well with empirical results for free shear flows.

Oh, Y. H.↗

A prediction method for velocity and temperature profiles in a two-dimensional nominally steady turbulent boundary layer

This paper describes a recently developed boundary-layer prediction method for a variable property compressible flow, in which heat transfer takes place primarily by forced convection and for which the mainstream Mach number is small. The leading order terms, in asymptotic expansions for large Reynolds numbers, are obtained for the mean velocity and temperature distribution in both the inner and outer layer of the turbulent boundary layer. Closure in the inner layer is achieved using an analytical model for the mean profiles which is based on the observed coherent structure of the time-dependent inner layer flow. For the outer layer, simple eddy viscosity and conductivity models are developed without recourse to the Reynolds analogy. In the prediction method a numerical solution of the outer layer equations is matched to the analytical inner layer profiles as the computation procedes downstream. Calculations are presented for a range of adverse and favorable pressure gradient flows and the predicted results compare well with existing data.

Weigand, G. G.↗

Current QSO statistics - Implications for the intergalactic medium

The results of numerous QSO surveys have been compiled and analyzed to form a single spatially averaged QSO ionizing function that is independent of evolution mode (number density or luminosity). An intergalactic medium (IGM) that satisfies the observational constraints may now have a density with respect to closure of approximately 0.1, which is substantially less than hitherto modeled. The aggregate of QSO data also indicate that if evolution is due to number density, then, when split into luminosity classes, an exponential in look-back time fits the data better than a power law, and the evolution rates increase roughly with absolute luminosity if analyzed by the exponential.

Sherman, R. D.↗

Three-dimensional flow over a conical afterbody containing a centered propulsive jet - A numerical simulation

The supersonic flow field over a body of revolution incident to the free stream is simulated numerically on a large, array processor (the CDC Cyber 205). The configuration is composed of a cone-cylinder forebody followed by a conical afterbody from which emanates a centered, supersonic propulsive jet. The free-stream Mach number is 2, the jet-exit Mach number is 2.5, and the jet-to-free-stream static pressure ratio is 3. Both the external flow and the exhaust are ideal air at a common total temperature. The thin-layer approximation to the time-dependent, compressible, Reynolds-averaged Navier-Stokes equations are solved using an implicit finite-difference algorithm. The data base, of 5 million words, is structured in a 'pencil' format so that efficient use of the array processor can be realized. The computer code is completely vectorized to take advantage of the data structure. Turbulence closure is achieved using an empirical algebraic eddy-viscosity model. The configuration and flow conditions correspond to published experimental tests and the computed solutions are consistent with the experimental data.

Deiwert, G. S.↗

Numerical Analysis Of Interlaminar-Fracture Toughness

Finite-element analysis applied in conjunction with strain-energy and micromechanical concepts. Computational procedure involves local, local-crack-closure, and/or the "unique" local-crack-closure method developed at NASA Lewis Research Center, for mathematical modeling of ENF and MMF. Methods based on three-dimensional finite-element analysis in conjunction with concept of strain-energy-release rate and with micromechanics of composite materials. Assists in interpretation of ENF and MMF fracture tests performed to obtain fracture-toughness parameters, by enabling evaluation of states of stress likely to induce interlaminar fractures.

Chamis, C. C.↗

Navier-Stokes calculation of solid-propellant rocket motor internal flowfields

A comprehensive numerical analysis has been carried out to study the detailed physical and chemical processes involved in the combustion of homogeneous propellant in a rocket motor. The formulation is based on the time-dependent full Navier-Stokes equations, with special attention devoted to the chemical reactions in both gas and condensed phases. The turbulence closure is achieved using both the Baldwin-Lomax algebraic model and a modified k-epsilon two-equation scheme with a low Reynolds number and near-wall treatment. The effects of variable thermodynamic and transport properties are also included. The system of governing equations are solved using a multi-stage Runge-Kutta shceme with the source terms treated implicitly. Preliminary results clearly demonstrate the presence of various combustion regimes in the vicinity of propellant surface. The effects of propellant combustion on the motor internal flowfields are investigated in detail.

Hsieh, Kwang-Chung↗

Fatigue crack propagation in aerospace aluminum alloys

This paper reviews fracture mechanics based, damage tolerant characterizations and predictions of fatigue crack growth in aerospace aluminum alloys. The results of laboratory experimentation and modeling are summarized in the areas of: (1) fatigue crack closure, (2) the wide range crack growth rate response of conventional aluminum alloys, (3) the fatigue behavior of advanced monolithic aluminum alloys and metal matrix composites, (4) the short crack problem, (5) environmental fatigue, and (6) variable amplitude loading. Remaining uncertainties and necessary research are identified. This work provides a foundation for the development of fatigue resistant alloys and composites, next generation life prediction codes for new structural designs and extreme environments, and to counter the problem of aging components.

Gangloff, R. P.↗

Magnetosphere-ionosphere coupling during plasmoid evolution - First results

The influence of magnetosphere-ionosphere coupling on the dynamic evolution of the magnetotail is investigated by a three-dimensional resistive MHD code that includes the effects of the closure of field-aligned currents in a simple resistive model ionosphere. Particular emphasis is on the role of this coupling during substorm evolution and the modification of the latter by the convection driven by the ionospheric electric fields. For comparison, results are presented from a simulation which uses an infinitely conducting ionosphere but is otherwise identical. Comparison of the two simulations shows that the major impact of magnetosphere-ionospheric communication is an acceleration of magnetotail evolution. Otherwise, phenomena in the two models are qualitatively similar. It is concluded that ionospheric effects do not significantly affect substorm associated magnetotail dynamics.

Hesse, Michael↗

Adaptation of the Advanced Spray Combustion Code to Cavitating Flow Problems

A very important consideration in turbopump design is the prediction and prevention of cavitation. Thus far conventional CFD codes have not been generally applicable to the treatment of cavitating flows. Taking advantage of its two-phase capability, the Advanced Spray Combustion Code is being modified to handle flows with transient as well as steady-state cavitation bubbles. The volume-of-fluid approach incorporated into the code is extended and augmented with a liquid phase energy equation and a simple evaporation model. The strategy adopted also successfully deals with the cavity closure issue. Simple test cases will be presented and remaining technical challenges will be discussed.

Liang, Pak-Yan↗

Aerosol-Induced Radiative Flux Changes Off the United States Mid-Atlantic Coast: Comparison of Values Calculated from Sunphotometer and In Situ Data with Those Measured by Airborne Pyranometer

The Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX) measured a variety of aerosol radiative effects (including flux changes) while simultaneously measuring the chemical, physical, and optical properties of the responsible aerosol particles. Here we use TARFOX-determined aerosol and surface properties to compute shortwave radiative flux changes for a variety of aerosol situations, with midvisible optical depths ranging from 0.06 to 0.55. We calculate flux changes by several techniques with varying degrees of sophistication, in part to investigate the sensitivity of results to computational approach. We then compare computed flux changes to those determined from aircraft measurements. Calculations using several approaches yield downward and upward flux changes that agree with measurements. The agreement demonstrates closure (i.e. consistency) among the TARFOX-derived aerosol properties, modeling techniques, and radiative flux measurements. Agreement between calculated and measured downward flux changes is best when the aerosols are modeled as moderately absorbing (midvisible single-scattering albedos between about 0.89 and 0.93), in accord with independent measurements of the TARPOX aerosol. The calculated values for instantaneous daytime upwelling flux changes are in the range +14 to +48 W/sq m for midvisible optical depths between 0.2 and 0.55. These values are about 30 to 100 times the global-average direct forcing expected for the global-average sulfate aerosol optical depth of 0.04. The reasons for the larger flux changes in TARFOX include the relatively large optical depths and the focus on cloud-free, daytime conditions over the dark ocean surface. These are the conditions that produce major aerosol radiative forcing events and contribute to any global-average climate effect.

Russell, P. B.↗

A Damage Model for the Simulation of Delamination in Advanced Composites under Variable-Mode Loading

A thermodynamically consistent damage model is proposed for the simulation of progressive delamination in composite materials under variable-mode ratio. The model is formulated in the context of Damage Mechanics. A novel constitutive equation is developed to model the initiation and propagation of delamination. A delamination initiation criterion is proposed to assure that the formulation can account for changes in the loading mode in a thermodynamically consistent way. The formulation accounts for crack closure effects to avoid interfacial penetration of two adjacent layers after complete decohesion. The model is implemented in a finite element formulation, and the numerical predictions are compared with experimental results obtained in both composite test specimens and structural components.

Turon, A.↗

Evaluating Model Parameterizations of Submicron Aerosol Scattering and Absorption with in situ Data from ARCTAS 2008

Accurate modeling of the scattering and absorption of ultraviolet and visible radiation by aerosols is essential for accurate simulations of atmospheric chemistry and climate. Closure studies using in situ measurements of aerosol scattering and absorption can be used to evaluate and improve models of aerosol optical properties without interference from model errors in aerosol emissions, transport, chemistry, or deposition rates. Here we evaluate the ability of four externally mixed, fixed size distribution parameterizations used in global models to simulate submicron aerosol scattering and absorption at three wavelengths using in situ data gathered during the 2008 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) campaign. The four models are the NASA Global Modeling Initiative (GMI) Combo model, GEOS-Chem v9- 02, the baseline configuration of a version of GEOS-Chem with online radiative transfer calculations (called GC-RT), and the Optical Properties of Aerosol and Clouds (OPAC v3.1) package. We also use the ARCTAS data to perform the first evaluation of the ability of the Aerosol Simulation Program (ASP v2.1) to simulate submicron aerosol scattering and absorption when in situ data on the aerosol size distribution are used, and examine the impact of different mixing rules for black carbon (BC) on the results. We find that the GMI model tends to overestimate submicron scattering and absorption at shorter wavelengths by 10-23 percent, and that GMI has smaller absolute mean biases for submicron absorption than OPAC v3.1, GEOS-Chem v9-02, or GC-RT. However, the changes to the density and refractive index of BC in GCRT improve the simulation of submicron aerosol absorption at all wavelengths relative to GEOS-Chem v9-02. Adding a variable size distribution, as in ASP v2.1, improves model performance for scattering but not for absorption, likely due to the assumption in ASP v2.1 that BC is present at a constant mass fraction throughout the aerosol size distribution. Using a core-shell mixing rule in ASP overestimates aerosol absorption, especially for the fresh biomass burning aerosol measured in ARCTAS-B, suggesting the need for modeling the time-varying mixing states of aerosols in future versions of ASP.

Alvarado, Matthew J.↗

Kinetic Model Development of the Carbon Vapor Deposition Reactor to Predict Performance Versus Design

To enable complete closure of a Sabatier-based atmosphere revitalization system where oxygen recovery is predicated by carbon dioxide reduction to methane and water, there must be a means of recovering the hydrogen from the byproduct methane. Honeywell’s carbon vapor deposition (CVD) reactor is a potential technology that provides a thermochemical means by which methane decomposes at high temperatures (> 1,000 °C) into hydrogen, solid carbon, and a small percentage of hydrocarbons. This advanced methane pyrolysis reactor utilizes a carbon fiber substrate to capture and immobilize the solid pyrolytic carbon rather than allowing continuation toward undesirable soot. The CVD reactor is a highly dynamic system whose characteristic geometry changes while operating. The reactor model was developed to provide predictive performance assessment capabilities at different operating conditions, geometries, and configurations to aid in technology development. Modeling of the CVD reactor considered the multitude of simultaneous and interdependent physicochemical processes that are entailed, including: (1) the homogeneous gas-phase pyrolysis reactions, (2) the heterogeneous deposition reactions, and (3) the transport of gaseous species in the reactor domains. A 1-D CVD reactor model was developed based on a general packed-bed, axial-dispersion plug flow reactor model with a solid substrate balance to account for carbon deposition. The model includes 59 species and 243 reactions constituting the methane pyrolysis reaction network, carbon deposition, and interconversions of oxygenate species. The reactor model was compared to brassboard test data, where it was shown to be able to capture the pertinent physicochemical phenomena as evidenced by its agreement with short-duration run data (i.e., looking at different temperature and flow rate conditions), as well as trend agreement with the brassboard reactor’s long-duration test data. Recommendations are provided on further model development steps as well as its use cases.

Air Revitalization↗