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.

SourceResearch connectionFrom yearThrough year

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

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Effects of expansion ratio on the calculation of parallel-walled backward-facing step flows - Comparison of four models of turbulence

This paper presents a numerical study of the effects of expansion ratio on two-dimensional separating and reattaching flows in plane backward-facing step geometries with parallel walls. Closure of the Reynolds equations was achieved by four different turbulence models: k-epsilon, 'modified' k-epsilon, algebraic stress, and 'modified' algebraic stress models. The k-epsilon model relates the Reynolds stresses to the mean rate of strain through the definition of an isotropic turbulent viscosity. The more advanced algebraic stress model calculates the stresses from implicit algebraic relationships containing the stresses themselves, the mean rate of strain, and the turbulent kinetic energy and its dissipation rate. 'Modified' versions of the models employ a new dissipation rate equation whose production term was made more sensitive to streamwise curvature effects. A new nonequilibrium wall function treatment proposed by Chieng and Launder (1980) was also incorporated into each model.

Sindir, M. M.↗

NASA NTRS1986

Viscoelastic deformation near active plate boundaries

Model deformations near the active plate boundaries of Western North America using space-based geodetic measurements as constraints are discussed. The first six months of this project were spent gaining familarity with space-based measurements, accessing the Crustal Dynamics Data Information Computer, and building time independent deformation models. The initial goal was to see how well the simplest elastic models can reproduce very long base interferometry (VLBI) baseline data. From the Crustal Dynamics Data Information Service, a total of 18 VLBI baselines are available which have been surveyed on four or more occasions. These data were fed into weighted and unweighted inversions to obtain baseline closure rates. Four of the better quality lines are illustrated. The deformation model assumes that the observed baseline rates result from a combination of rigid plate tectonic motions plus a component resulting from elastic strain build up due to a failure of the plate boundary to slip at the full plate tectonic rate. The elastic deformation resulting from the locked plate boundary is meant to portray interseismic strain accumulation. During and shortly after a large interplate earthquake, these strains are largely released, and points near the fault which were previously retarded suddenly catch up to the positions predicted by rigid plate models. Researchers judge the quality of fit by the sum squares of weighted residuals, termed total variance. The observed baseline closures have a total variance of 99 (cm/y)squared. When the RM2 velocities are assumed to model the data, the total variance increases to 154 (cm/y)squared.

Ward, S. N.↗

NASA NTRS2021

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗

NASA NTRS2013

Characterization of Swirl-Venturi Lean Direct Injection Designs for Aviation Gas-Turbine Combustion

Injector geometry, physical mixing, chemical processes, and engine cycle conditions together govern performance, operability and emission characteristics of aviation gas-turbine combustion systems. The present investigation explores swirl-venturi lean direct injection combustor fundamentals, characterizing the influence of key geometric injector parameters on reacting flow physics and emission production trends. In this computational study, a design space exploration was performed using a parameterized swirl-venturi lean direct injector model. From the parametric geometry, 20 three-element lean direct injection combustor sectors were produced and simulated using steady-state, Reynolds-averaged Navier-Stokes reacting computations. Species concentrations were solved directly using a reduced 18-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear κ-ε model. Results demonstrate sensitivities of the geometric perturbations on axially averaged flow field responses. Output variables include axial velocity, turbulent kinetic energy, static temperature, fuel patternation and minor species mass fractions. Significant trends have been reduced to surrogate model approximations, intended to guide future injector design trade studies and advance aviation gas-turbine combustion research.

Combustion Chamber↗

NASA NTRS1984

An elastic failure model of indentation damage

A mechanistically consistent model for indentation damage based on elastic failure at tensile or shear overloads, is proposed. The model accommodates arbitrary crack orientation, stress relaxation, reduction and recovery of stiffness due to crack opening and closure, and interfacial friction due to backward sliding of closed cracks. This elastic failure model was implemented by an axisymmetric finite element program which was used to simulate progressive damage in a silicon nitride plate indented by a tungsten carbide sphere. The predicted damage patterns and the permanent impression matched those observed experimentally. The validation of this elastic failure model shows that the plastic deformation postulated by others is not necessary to replicate the indentation damage of brittle structural ceramics.

Liaw, B. M.↗

NASA NTRS1997

Dynamical System Analysis of Reynolds Stress Closure Equations

In this paper, we establish the causality between the model coefficients in the standard pressure-strain correlation model and the predicted equilibrium states for homogeneous turbulence. We accomplish this by performing a comprehensive fixed point analysis of the modeled Reynolds stress and dissipation rate equations. The results from this analysis will be very useful for developing improved pressure-strain correlation models to yield observed equilibrium behavior.

Girimaji, Sharath S.↗

NASA NTRS1991

Modeling of compressible turbulent shear flows

Despite all the recent developments in computer technologies and numerical algorithms, full numerical simulations of turbulent flows are feasible only at moderate Reynolds numbers and for flows with relatively simple geometries. The main goal of this research is to develop new second order moment closures for compressible turbulence. It has been shown that the models based on the extension of those developed originally for incompressible flows fail to adequately predict turbulent flows at high Mach numbers. In this attempt, the compressibility effects are explicitly considered. A successful development of these models that directly takes into account the compressibility effects may have a range of technological implications in the design of supersonic and hypersonic vehicles.

Liou, William W.↗

NASA NTRS2021

Discrete anisotropic radiative transfer modelling of solar-induced chlorophyll fluorescence: Structural impacts in geometrically explicit vegetation canopies

Solar-induced fluorescence (SIF) is a subtle but informative optical signal of vegetation photosynthesis. Remotely sensed SIF integrates environmental, physiological and structural changes that alter photosynthesis at leaf, plant and canopy scales. Radiative transfer models are ideally suited to investigate the complex sources of variability in the SIF signal to guide the interpretation of SIF retrievals from airborne and space-borne platforms. Here, we coupled the Fluspect-Cx model of leaf optical properties and chlorophyll-a fluorescence with the Discrete Anisotropic Radiative Transfer (DART) model to upscale SIF from individual leaves to three-dimensional (3D) structurally explicit canopies. For one-dimensional homogeneous (turbid-like) canopies, DART-SIF was nearly identical to SIF simulated in two existing models, SCOPE and mSCOPE (RMSE <0.221 W./sq.m.μm.sr). DART simulations in geometrically explicit 3D canopies offered four important insights regarding the influence of vegetation structure on the multi-angular top-of-canopy SIF signal. First, changes in the 3D canopy architecture of maize crops, represented by leaf density (leaf area index), and plant clumping (canopy closure) had a larger impact on SIF than the modelled photosynthetic efficiency distinction between sun-adapted and shade-adapted foliage. Second, clumping of leaves at the crop and stand levels was identified as one of the key driving factors of multi-angular anisotropy of red and far-red SIF (686 and 740 nm) for both maize and eucalyptus canopies. Third, non-photosynthetic woody material had a significant impact on top-of-canopy SIF in modelled 3D forest stands. Wood shadowing decreased the photosynthetically active radiation absorbed by green leaves, and consequently the SIF emissions, by 10% in sparse and 17% in dense eucalyptus stands. The wood obstruction (blocking) effect, quantified as a relative difference of SIF escape probabilities from canopies with and without wood in the nadir viewing direction, decreased far-red SIF by 4–6% but it had a smaller and sometimes positive influence (by less than 2%) on red SIF. Fourth, DART 3D radiative budget profiles revealed that the majority of the SIF signal from a dense eucalyptus stand originated from the top 25% of the simulated canopy. Interestingly, the introduction of bark-covered woody elements did not alter the simulated balance and omnidirectional escape factor of red SIF in this upper canopy part but did raise significantly both of them in case of far-red SIF. These results demonstrate the importance of 3D radiative transfer and radiative budget simulations for investigating SIF interactions in structurally complex plant canopies and for a better understanding of spatiotemporal and multi-angular remote sensing SIF observations.

Zbyněk Malenovský↗

NASA NTRS2017

A Numerical Study on the Edgewise Compression Strength of Sandwich Structures with Facesheet-Core Disbonds

Damage tolerant design approaches require determination of critical damage modes and flaw sizes in order to establish nondestructive evaluation detection requirements. A finite element model is developed to assess the effect of circular facesheet-core disbonds on the strength of sandwich specimens subjected to edgewise compressive loads for the purpose of predicting the critical flaw size for a variety of design parameters. Postbuckling analyses are conducted in which an initial imperfection is seeded using results from a linear buckling analysis. Both the virtual crack closure technique (VCCT) and cohesive elements are considered for modeling disbond growth. Predictions from analyses using the VCCT and analyses using cohesive elements are in good correlation. A series of parametric analyses are conducted to investigate the effect of core thickness and material, facesheet layup, facesheet-core interface properties, and curvature on the criticality of facesheet-core disbonds of various sizes. The results from these analyses provide a basis for determining the critical flaw size for facesheet-core disbonds subjected to edgewise compression loads and, therefore, nondestructive evaluation flaw detection requirements for this configuration.

Bergan, Andrew C.↗

NASA NTRS2012

NASA's Carbon Monitoring System Flux-Pilot Project: A Multi-Component Analysis System for Carbon-Cycle Research and Monitoring

The importance of greenhouse gas increases for climate motivates NASA s observing strategy for CO2 from space, including the forthcoming Orbiting Carbon Observatory (OCO-2) mission. Carbon cycle monitoring, including attribution of atmospheric concentrations to regional emissions and uptake, requires a robust modeling and analysis infrastructure to optimally extract information from the observations. NASA's Carbon-Monitoring System Flux-Pilot Project (FPP) is a prototype for such analysis, combining a set of unique tools to facilitate analysis of atmospheric CO2 along with fluxes between the atmosphere and the terrestrial biosphere or ocean. NASA's analysis system is unique, in that it combines information and expertise from the land, oceanic, and atmospheric branches of the carbon cycle and includes some estimates of uncertainty. Numerous existing space-based missions provide information of relevance to the carbon cycle. This study describes the components of the FPP framework, assessing the realism of computed fluxes, thus providing the basis for research and monitoring applications. Fluxes are computed using data-constrained terrestrial biosphere models and physical ocean models, driven by atmospheric observations and assimilating ocean-color information. Use of two estimates provides a measure of uncertainty in the fluxes. Along with inventories of other emissions, these data-derived fluxes are used in transport models to assess their consistency with atmospheric CO2 observations. Closure is achieved by using a four-dimensional data assimilation (inverse) approach that adjusts the terrestrial biosphere fluxes to make them consistent with the atmospheric CO2 observations. Results will be shown, illustrating the year-to-year variations in land biospheric and oceanic fluxes computed in the FPP. The signals of these surface-flux variations on atmospheric CO2 will be isolated using forward modeling tools, which also incorporate estimates of transport error. The results will be discussed in the context of interannual variability of observed atmospheric CO2 distributions.

Pawson, S.↗

NASA NTRS2015

Revisiting Turbulence Model Validation for High-Mach Number Axisymmetric Compression Corner Flows

Two axisymmetric shock-wave/boundary-layer interaction (SWBLI) cases are used to benchmark one- and two-equation Reynolds-averaged Navier-Stokes (RANS) turbulence models. This validation exercise was executed in the philosophy of the NASA Turbulence Modeling Resource and the AIAA Turbulence Model Benchmarking Working Group. Both SWBLI cases are from the experiments of Kussoy and Horstman for axisymmetric compression corner geometries with SWBLI inducing flares of 20 and 30 degrees, respectively. The freestream Mach number was approximately 7. The RANS closures examined are the Spalart-Allmaras one-equation model and the Menter family of kappa − omega two equation models including the Baseline and Shear Stress Transport formulations. The Wind-US and CFL3D RANS solvers are employed to simulate the SWBLI cases. Comparisons of RANS solutions to experimental data are made for a boundary layer survey plane just upstream of the SWBLI region. In the SWBLI region, comparisons of surface pressure and heat transfer are made. The effects of inflow modeling strategy, grid resolution, grid orthogonality, turbulent Prandtl number, and code-to-code variations are also addressed.

Shock Wave↗

NASA NTRSDate not supplied

Exploring Radiation, Water and Energy Closure in Key Deep Tropical Regions Over the Last 3 Decades

The closure of global energy and water budgets at various time scales depends heavily upon the on the processes of water and energy exchange in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance in key tropical regions is observed of water budget terms between observations and atmospheric reanalysis systems, implying that processes are not simulated well. For instance, periods of imbalance in the water budget may be at least partly attributed to precipitation estimates due to the transition between cloud and precipitation systems in which ice cloud affects microwave radar estimates (Lietmann-Niimi et al. 2024, submitted). A recent paper by Hsaio et al (2023, pre-published version) found that the transition process is characterized by wind shear and longwave cloud radiative feedback signatures as well. This is just one example how these processes also affect the radiative fluxes and thus the energy budget at the TOA. In this presentation, we assess the net radiative budget and components in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES EBAF/SYN1deg Ed4.2, ISSCP FH, etc.) and describing the temporal variability of energy balance terms across the deep tropics for the period from 1988 to near present in six key oceanic regions. These six key tropical oceanic regions are West Indian Ocean, the south Indian Ocean east, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean. In these regions we assess the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear to assess the uncertainties of energy closure in these regions. These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. These differences also affect closure terms and inferences on systematic uncertainties within the observational estimates. We conclude with recommendations for continuing work

Radiation budget↗

NASA NTRS1988

Theoretical study of turbulent channel flow: Bulk properties, pressure fluctuations, and propagation of electromagnetic waves

In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr--Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O--S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend stongly on the structure of the turbulence spectrun at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.

Canuto, V. M.↗

NASA NTRS1990

Theoretical study of turbulent channel flow - Bulk properties, pressure fluctuations, and propagation of electromagnetic waves

In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr-Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O-S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend strongly on the structure of the turbulence spectrum at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.

Canuto, V. M.↗

NASA NTRS1990

On near-wall turbulent flow modelling

The behavior of various terms in the Reynolds-stress transport equations is evaluated in view results of near-wall turbulence characterizations which show that all components of the velocity-pressure gradient correlation vanish at the wall. A methodology for deriving an asymptotically correct model of the velocity-pressure gradient correlation is proposed, and a model capable of approaching the high-Reynolds number model for pressure distribution far from the wall is derived. A parametric study of the model constants introduced by the near-wall closure reveals that one constant in the dissipation-rate equation is Reynolds number-dependent; with this modification, excellent agreement is obtained with measured and simulated near-wall turbulence statistics.

Lai, Y. G.↗

NASA NTRS1992

Long-term orbital evolution of short-period comets found in Project Cosmo-DICE

Orbital evolutions of about 160 short-period (SP) comets are numerically integrated for 4400 years in the framework of a realistic dynamical model. By the round-trip error in closure test, a reliable time space of the integrated orbits is estimated for each comet. Majority of the SP comets with their Tisserand's constant(J) between 2.8 and 3.1 are found to evolve within the past 1000-2000 years from the orbits whose perihelia are near the Jovian orbit to the orbits with perihelia of 1-2 AU. This evolution is much more rapid than that expected from Monte Carlo simulations based on symmetric distribution of planetary perturbations, thus suggesting that asymmetry of perturbation distribution play an important role in cometary evolution. Several comets are shown to evolve from the near-Saturn orbits and then to be handed over under the control of Jupiter. We also find that a few comets were captured from long-period orbits (a = 75-125 AU) via only a few close encounters with Jupiter. It is confirmed that the captured SP comets of low-inclination with 2.7 less than J less than 3.1 show more or less strong chaotic behavior. On the other hand, comets with longer orbital period and/or of high inclination reveal slow or quasi-periodic orbital evolution.

Nakamura, Tsuko↗

NASA NTRS1994

Rotorcraft control system design for uncertain vehicle dynamics using quantitative feedback theory

Quantitative Feedback Theory describes a frequency-domain technique for the design of multi-input, multi-output control systems which must meet time or frequency domain performance criteria when specified uncertainty exists in the linear description of the vehicle dynamics. This theory is applied to the design of the longitudinal flight control system for a linear model of the BO-105C rotorcraft. Uncertainty in the vehicle model is due to the variation in the vehicle dynamics over a range of airspeeds from 0-100 kts. For purposes of exposition, the vehicle description contains no rotor or actuator dynamics. The design example indicates the manner in which significant uncertainty exists in the vehicle model. The advantage of using a sequential loop closure technique to reduce the cost of feedback is demonstrated by example.

Hess, R. A.↗

NASA NTRS2000

Long-Term Global Morphology of Gravity Wave Activity Using UARS Data

Gravity waves in satellite data from CRISTA and MLS are studied in depth this quarter. Results this quarter are somewhat limited due to the PI'S heavy involvement throughout this reporting period in on-site forecasting of mountain wave-induced turbulence for the NASA's ER-2 research aircraft at Kiruna, Sweden during the SAGE Ill Ozone Loss and Validation Experiment (SOLVE). Results reported concentrate on further mesoscale modeling studies of mountain waves over the southern Andes, evident in CRISTA and MLS data. Two-dimensional mesoscale model simulations are extended through generalization of model equations to include both rotation and a first-order turbulence closure scheme. Results of three experiments are analyzed in depth and submitted for publication. We also commence simulations with a three-dimensional mesoscale model (MM5) and present preliminary results for the CRISTA 1 period near southern South America. Combination of ground-based temperature data at 87 km from two sites with global HRDl data was continued this quarter, showing stationary planetary wave structures. This work was also submitted for publication.

Eckermann, Stephen D.↗