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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

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.↗
NASA NTRS1981

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.↗
NASA NTRS1983

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.↗
NASA NTRS1988

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.↗
NASA NTRS1988

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↗
NASA NTRS1990

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.↗
NASA NTRS1991

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↗
NASA NTRS1993

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↗
NASA NTRS2000

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.↗
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2006

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.↗
NASA NTRS2016

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.↗
DOE OSTI2023

Utah FORGE 5-2615: Thermo-poromechanical Response of Fractured Rock - Workshop Presentation

This is a presentation on the Experimental Determination and Modeling-Informed Analysis of Thermo-poromechanical Response of Fractured Rock for Application to Utah FORGE project by the University of Oklahoma, presented by Dr. Ahmad Ghassemi, McCasland Chair Prof. The project objective is to improve understanding and control of coupled thermo-poromechanical (or thermo-hydro-mechanical- TPM) processes in reservoir development, and to study any role in interpretations of the fracture closure. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
NASA NTRSDate not supplied

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↗
DOE OSTI2021

Aerosol Direct Radiative Effects at the ARM SGP and TWP Sites: Clear Skies

The clear-sky aerosol direct radiative effect (DRE) was estimated at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) and Tropical Western Pacific (TWP) sites. The NASA Langley Fu-Liou radiation model was used with observed inputs including aerosol vertical extinction profile from the Raman lidar; spectral aerosol optical depth (AOD), single-scattering albedo and asymmetry factor from Aerosol Robotic Network; temperature and water vapor profiles from radiosondes; and surface shortwave (SW) spectral albedo from radiometers. A radiative closure experiment was conducted for clear-sky conditions. The mean differences of modeled and observed surface downwelling SW total fluxes were 1 W m-2 at SGP and 2 W m -2 at TWP, which are within observational uncertainty. At SGP, the estimated annual mean clear-sky aerosol DRE is -3.00 W m -2 at the top of atmosphere (TOA) and -6.85 W m -2 at the surface. The strongest aerosol DRE of -4.81 (-10.77) W m -2 at the TOA (surface) are in the summer when AODs are largest. The weakest aerosol DRE of -.28 (-2.77) W m -2 at the TOA (surface) are in November–January when AODs and single-scattering albedos are lowest. At TWP, the annual mean clear-sky DRE is -2.82 W m-2 at the TOA and -10.34 W m -2 at the surface. The strongest aerosol DRE of -5.95 (-22.20) W m -2 at the TOA (surface) are in November (October) due to the biomass burning season’s peak. The weakest aerosol DRE of -0.96 (-4.16) W m -2 at the TOA (surface) are in March (April) when AODs are smallest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

Interaction of Delaminations and Matrix Cracks in a CFRP Plate, Part II: Simulation Using an Enriched Shell Finite Element Model

Numerical simulations are presented of a recently developed test which creates multiple delaminations in a CFRP laminate specimen that grow and interact via transverse matrix cracks. A novel shell element enriched with the Floating Node Method, and a damage algorithm based on the Virtual Crack Closure Technique, were used to successfully simulate the tests. Additionally, a 3D high mesh fidelity model based on cohesive zones and continuum damage mechanics was used to simulate the tests and act as a representative of other similar state-of-the-art high mesh fidelity modeling techniques to compare to the enriched shell element. The enriched shell and high mesh fidelity models had similar levels of accuracy and generally matched the experimental data. With runtimes of 36 minutes for the shell model and 55 hours for the high mesh fidelity model, the shell model is 92 times faster than the high- fidelity simulation.

Mark W McElroy↗
NASA NTRS1987

The role of turbulent convection in the primitive solar nebula. I - Theory. II - Results

A model of convective turbulence which takes radiative dissipation, rotation, and convective motion anisotropy into account, on the basis of a closure for the nonlinear interactions that employs the growth rates of hydrodynamic instabilities, is used to obtain a theoretical framework for modeling the primordial solar nebula. It is assumed that convection is the sole source of turbulence causing the solar nebula to evolve. Vertical structure equations in the thin disk approximation are developed and a detailed comparison with the previous solar nebula convective models of such workers as Lin et al. (1981, 1982) is undertaken. The present values for the turbulent efficiency are much lower and more sensitive to opacity and surface density, resulting in low turbulent speeds, a more massive disk, a lower accretion rate 'best value', and a longer characteristic dispersal time for the disk. It is concluded that convection may not be the dominant source of turbulence needed to evolve young solar/stellar nebulae.

Cabot, W.↗
DOE OSTI2023

Memory-based parameterization with differentiable solver: Application to Lorenz ’96

Physical parameterizations (or closures) are used as representations of unresolved subgrid processes within weather and global climate models or coarse-scale turbulent models, whose resolutions are too coarse to resolve small-scale processes. These parameterizations are typically grounded on physically based, yet empirical, representations of the underlying small-scale processes. Machine learning-based parameterizations have recently been proposed as an alternative solution and have shown great promise to reduce uncertainties associated with the parameterization of small-scale processes. Yet, those approaches still show some important mismatches that are often attributed to the stochasticity of the considered process. This stochasticity can be due to coarse temporal resolution, unresolved variables, or simply to the inherent chaotic nature of the process. To address these issues, we propose a new type of parameterization (closure), which is built using memory-based neural networks, to account for the non-instantaneous response of the closure and to enhance its stability and prediction accuracy. We apply the proposed memory-based parameterization, with differentiable solver, to the Lorenz ’96 model in the presence of a coarse temporal resolution and show its capacity to predict skillful forecasts over a long time horizon of the resolved variables compared to instantaneous parameterizations. This approach paves the way for the use of memory-based parameterizations for closure problems.

Mathematics↗