Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient structure”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Substrate oxygen sponge effect: A parameter for epitaxial manganite thin film growth

The emergent phenomena in complex oxide thin films are strongly tied to the oxygen content, which is often engineered by the oxygen partial pressure during growth. However, such oxygen control by the growth pressure is challenging to synthesize for some oxide films, which requires a subtle control of the oxygen content. A parameter of controlling the oxygen content independent of the growth pressure is desired. Here, we propose a method of controlling the oxygen content of films by engineering the substrate before the growth. The oxide substrate serves as an oxygen sponge, which provides a tunable oxygen environment ranging from oxygen-rich to oxygen-poor for the film growth, depending on the pre-substrate annealing (PSA) conditions. Using manganite as a model system, we demonstrate that this simple PSA method leads to remarkable changes in the structure and physical properties of the as-grown films. Finally, this substrate oxygen sponge effect, driven by the large oxygen concentration gradient at high temperatures, can be applied to explore not only emergent interfacial phenomena but also the growth of a variety of functional oxide thin films and nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Parallelizable adjoint stencil computations using transposed forward-mode algorithmic differentiation

Algorithmic differentiation (AD) is a tool for generating discrete adjoint solvers, which efficiently compute gradients of functions with many inputs, for example for use in gradient-based optimization. AD is often applied to large computations such as stencil operators, which are an important part of most structured-mesh PDE solvers. Stencil computations are often parallelized, for example by using OpenMP, and optimized by using techniques such as cache-blocking and tiling to fully utilize multicore CPUs and many-core accelerators and GPUs. Differentiating these codes with conventional reverse-mode AD results in adjoint codes that cannot be expressed as stencil operations and may not be easily parallelizable. They thus leave most of the compute power of modern architectures unused. We present a method that combines forward-mode AD and loop transformation to generate adjoint solvers that use the same memory access pattern as the original computation that they are derived from and can benefit from the same optimization techniques. The effectiveness of this method is demonstrated by generating a scalable adjoint CFD solver for multicore CPUs and Xeon Phi accelerators.

algorithmic differentiation↗
DOE OSTI2021

Nonlinear dynamics and transport driven by energetic particle instabilities using a gyro-Landau closure model *

Abstract Energetic particle (EP) destabilized Alfvén eigenmode (AE) instabilities are simulated for a DIII-D experimental case with a pulsed neutral beam using a gyro-Landau moments model which introduces EP phase-mixing effects through closure relations. This provides a computationally efficient reduced model which is applied here in the nonlinear regime over timescales that would be difficult to address with more complete models. The long timescale nonlinear evolution and related collective transport losses are examined including the effects of zonal flow/current generation, nonlinear energy cascades, and EP profile flattening. The model predicts frequencies and mode structures that are consistent with experimental observations. These calculations address issues that have not been considered in previous modelling: the EP critical gradient profile evolution in the presence of zonal flows/currents, and the dynamical nature of the saturated state. A strong level of intermittency is present in the predicted instability-driven transport; this is connected to the zonal flow growth and decay cycles and nonlinear energy transfers. Simulation of intermittent AE-enhanced EP transport will be an important issue for the protection of plasma facing components in the next generation of fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Studying the Recrystallization of Cu(InGa)Se 2 Semiconductor Thin Films by Silver Bromide In-situ Treatment

Cu(In,Ga)Se 2 samples were fabricated using a 3-stage thermal co-evaporation process on molybdenum back contact at low temperature. The process of recrystallization was carried out in between the 2 nd and 3 rd stages by flashing 25 mg of AgBr for 2 minutes. A change in morphological structure was observed as small grains transformed into large grains, as confirmed by XRD and SEM measurements. The decrease of the Ga gradient, seen in the SIMS depth profile, suggests Ga interdiffusion due to AgBr treatment. Altogether, the AgBr treatment contributes to a general improvement in device performance as compared to the as-deposited devices.

14 SOLAR ENERGY↗
DOE OSTI2019

Hawaii Play Fairway: Preliminary Core Box Photos, Lanai Island, Hawaii

Photos of core samples from Lanai Island. During the third phase of the Hawaii Play Fairway project, further exploration involved drilling a groundwater well in Lanai's Palawai Basin and performing more geophysical surveys. The project deepened an existing water well on Lanai. Drilling occurred 24/7 the entire month of June 2019 over which time Lanai Well 10 was deepened from 427 m to 1057 m, with continuous core collected. The roughly linear temperature gradient was an average of 42 degC/km, and a maximum bottom hole temperature, 66 degC. This gradient is more than twice the background for Hawaii and within a range of gradients measured in this depth range for some exploration wells within KERZ. The Hawaii Play Fairway project seeks to explore the geologic structures that exist in the caldera region of Hawaiian volcanoes; how those structures influence groundwater storage and flow; and how the magmatic heat from Hawaiian shield volcanoes cools over time. The Hawaii Groundwater and Geothermal Resources Center (Hawaii Institute of Geophysics and Planetology, University of Hawaii at Manoa) executed the Hawaii Play Fairway project. For more information, go to HGGRC's website that is linked in the resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Model Reference Adaptive Control (MRAC) for Additive Manufacturing (AM)

Model Reference Adaptive Control (MRAC) is based on the fundamental concept that the process under investigation is to be controlled to follow or “track” a reference system (model) characterized by a state/input/output model employing an adaptive optimization algorithm to adjust the controller parameters in real-time. The generic structure of the MRAC is shown in Fig. 1 consisting of the following primary components: Reference model, Process (system) model, controller and the adaption algorithm. The basic structure of the controller is specified by a linear construct with the corresponding real-time adaption algorithms given by a gradient-type (so-called MIT rule) or based on stability theory (Lyapunov, hyperstability). This approach to adaptive control is termed “direct”, since the controller (parameters) are adjusted based on the component models/algorithm in contrast to the “indirect” approach that adjusts the process model parameters applying real-time system identification techniques.

42 ENGINEERING↗
DOE OSTI2022

BioRT-Flux-PIHM v1.0: a biogeochemical reactive transport model at the watershed scale

Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

Numerical spatial marching techniques in duct acoustics

Direct calculation of the internal structure of a ducted noise source from farfield pressure measurements is regarded as an initial value problem, where the pressure and pressure gradient (farfield impedance) are assumed to be known along a line in the farfield. If pressure and impedance are known at the boundary of the farfield, the pressure can be uniquely determined in the vicinity of the inlet and inside the inlet ducting. A marching procedure is developed which, with this information obtained from measurements, enables a description of a ducted noise source. The technique uses a finite difference representation of the homogeneous Helmholtz equation.

Baumeister, K. J.↗
NASA NTRS1980

A three-dimensional model of corotating streams in the solar wind. II - Hydrodynamic streams

Theoretical aspects of corotating solar wind dynamics on a global scale are examined by means of numerical simulations executed with a nonlinear, inviscid, adiabatic, single fluid, three-dimensional hydrodynamic formulation. The factors most strongly affecting stream evolution, namely, the intrinsic correlations among density, temperature, and velocity existing near the source, the amplitude, longitudinal breadth and latitudinal breadth of the stream, and the heliographic latitude of the centroid of the stream, are considered. Calculations indicate that even in the presence of substantial meridional gradients, the latitudinal transport of mass, energy, and momentum by streams appears to have little impact on the overall evolution of the structure.

Pizzo, V. J.↗
NASA NTRS1989

Radiative forcing of the Venus mesosphere. II - Thermal fluxes, cooling rates, and radiative equilibrium temperatures

A radiative heat-transfer model is presently used to ascertain the way in which radiative forcing contributes to the up to 20 K higher temperature of the Venus polar regions, by comparison with the tropics, in the 60-100 km mesospheric levels. Model global-mean radiative equilibrium temperatures for 55-100 km are compared with observations to show how each opacity source contributes to the thermal structure. The results obtained from latitude-dependent radiative equilibrium experiments indicate that meridional variations in radiative forcing obliterate observed mesospheric temperature gradients and yield polar temperatures up to 40 K cooler than the tropics.

Crisp, David↗
NASA NTRS1990

Multi-disciplinary optimization of aeroservoelastic systems

Efficient analytical and computational tools for simultaneous optimal design of the structural and control components of aeroservoelastic systems are presented. The optimization objective is to achieve aircraft performance requirements and sufficient flutter and control stability margins with a minimal weight penalty and without violating the design constraints. Analytical sensitivity derivatives facilitate an efficient optimization process which allows a relatively large number of design variables. Standard finite element and unsteady aerodynamic routines are used to construct a modal data base. Minimum State aerodynamic approximations and dynamic residualization methods are used to construct a high accuracy, low order aeroservoelastic model. Sensitivity derivatives of flutter dynamic pressure, control stability margins and control effectiveness with respect to structural and control design variables are presented. The performance requirements are utilized by equality constraints which affect the sensitivity derivatives. A gradient-based optimization algorithm is used to minimize an overall cost function. A realistic numerical example of a composite wing with four controls is used to demonstrate the modeling technique, the optimization process, and their accuracy and efficiency.

Karpel, Mordechay↗
NASA NTRS1992

Multidisciplinary optimization of aeroservoelastic systems using reduced-size models

Efficient analytical and computational tools for simultaneous optimal design of the structural and control components of aeroservoelastic systems are presented. The optimization objective is to achieve aircraft performance requirements and sufficient flutter and control stability margins with a minimal weight penalty and without violating the design constraints. Analytical sensitivity derivatives facilitate an efficient optimization process which allows a relatively large number of design variables. Standard finite element and unsteady aerodynamic routines are used to construct a modal data base. Minimum State aerodynamic approximations and dynamic residualization methods are used to construct a high accuracy, low order aeroservoelastic model. Sensitivity derivatives of flutter dynamic pressure, control stability margins and control effectiveness with respect to structural and control design variables are presented. The performance requirements are utilized by equality constraints which affect the sensitivity derivatives. A gradient-based optimization algorithm is used to minimize an overall cost function. A realistic numerical example of a composite wing with four controls is used to demonstrate the modeling technique, the optimization process, and their accuracy and efficiency.

Karpel, Mordechay↗
NASA NTRS1992

Second order closure modeling of turbulent buoyant wall plumes

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Zhu, Gang↗
NASA NTRS1991

Stochastic inflation lattice simulations - Ultra-large scale structure of the universe

Non-Gaussian fluctuations for structure formation may arise in inflation from the nonlinear interaction of long wavelength gravitational and scalar fields. Long wavelength fields have spatial gradients, a (exp -1), small compared to the Hubble radius, and they are described in terms of classical random fields that are fed by short wavelength quantum noise. Lattice Langevin calculations are given for a toy model with a scalar field interacting with an exponential potential where one can obtain exact analytic solutions of the Fokker-Planck equation. For single scalar field models that are consistent with current microwave background fluctuations, the fluctuations are Gaussian. However, for scales much larger than our observable Universe, one expects large metric fluctuations that are non-Gaussian. This example illuminates non-Gaussian models involving multiple scalar fields which are consistent with current microwave background limits.

Salopek, D. S.↗
NASA NTRS1992

Oregon transect: Comparison of leaf-level reflectance with canopy-level and modelled reflectance

The Oregon Transect Ecosystem Research (OTTER) project involves the collection of a variety of remotely-sensed and in situ measurements for characterization of forest biophysical and biochemical parameters. The project includes nine study plots located along an environmental gradient in west-central Oregon, extending from the Pacific coast inland approximately 300km. These plots represent a broad range in ecosystem structure and function. Within the OTTER project, the sensitivity of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) signal to absorption by foliar biochemicals is being examined. AVIRIS data were acquired over all plots in conjunction with the four OTTER Multi-sensor Aircraft Campaigns spanning the growing season. Foilage samples were gathered during each campaign for biochemical determination (at Ames Research Center), to estimate stand-level constituency at each plot. Directional-hemispheric leaf reflectance throughout the 400-2400nm region was measured in the laboratory as an aid to interpreting concurrent AVIRIS data. Obtaining leaf spectra in this manner reduces or eliminates the confounding influences of atmosphere, canopy architecture, and reflectance by woody components, understory, and exposed soils which are present in airborne observations. These laboratory spectra were compared to simulated spectra derived by inverting the PROSPECT leaf-level canopy reflectance derived from AVIRIS data by use of the LOWTRAN-7 atmospheric radiative-transfer model.

Johnson, Lee F.↗
NASA NTRS1995

Solution of the three-dimensional Helmholtz equation with nonlocal boundary conditions

The Helmholtz equation is solved within a three-dimensional rectangular duct with a nonlocal radiation boundary condition at the duct exit plane. This condition accurately models the acoustic admittance at an arbitrarily-located computational boundary plane. A linear system of equations is constructed with second-order central differences for the Helmholtz operator and second-order backward differences for both local admittance conditions and the gradient term in the nonlocal radiation boundary condition. The resulting matrix equation is large, sparse, and non-Hermitian. The size and structure of the matrix makes direct solution techniques impractical; as a result, a nonstationary iterative technique is used for its solution. The theory behind the nonstationary technique is reviewed, and numerical results are presented for radiation from both a point source and a planar acoustic source. The solutions with the nonlocal boundary conditions are invariant to the location of the computational boundary, and the same nonlocal conditions are valid for all solutions. The nonlocal conditions thus provide a means of minimizing the size of three-dimensional computational domains.

Hodge, Steve L.↗
NASA NTRS1998

Ionospheric Instabilities and Shock Capturing Algorithms

We describe the basic physics of the gradient-drift instability and the collisional Rayleigh-Taylor instability, the two plasma instabilities which dominate the production of small-scale structure in the earth's ionosphere. We then focus our attention on the nonlinear evolution of these instabilities for two particular cases: 1) the recursive bifurcation and "freezing" of artificial ion clouds released in the ionosphere; and 2) equatorial spread F. Finally, we discuss the numerical challenges faced in attempting to numerically simulate these instabilities, and the role played by shock-capturing algorithms therein.

Zalesak, Steven T.↗
NASA NTRS2000

Analysis of the Radar Reflectivity of Aircraft Vortex Wakes

Radar has been proposed as a way to track wake vortices to reduce aircraft spacing and tests have revealed radar echoes from aircraft wakes in clear air. The results are always interpreted qualitatively using Tatarski's theory of weak scattering by isotropic atmospheric turbulence. The goal of the present work was to predict the value of the radar cross-section (RCS) using simpler models. This is accomplished in two steps. First, the refractive index is obtained. Since the structure of the aircraft wakes is different from atmospheric turbulence, three simple mechanisms specific to vortex wakes are considered: (1) Radial density gradient in a two-dimensional vortex, (2) three-dimensional fluctuations in the vortex cores, and (3) Adiabatic transport of the atmospheric fluid in a two-dimensional oval surrounding the pair of vortices. The index of refraction is obtained more precisely for the two-dimensional mechanisms than for the three-dimensional ones. In the second step, knowing the index of refraction, a scattering analysis is performed. Tatarski's weak scattering approximation is kept but the usual assumptions of a far-field and a uniform incident wave are dropped. Neither assumption is generally valid for a wake that is coherent across the radar beam. For analytical insight, a simpler approximation that invokes, in addition to weak scattering, the far-field and wide cylindrical beam assumptions, is also developed and compared with the more general analysis. The predicted RCS values for the oval surround the vortices (mechanism C) agree with the experiments of Bilson conducted over a wide range of frequencies. However, the predictions have a cut-off away from normal incidence which is not present in the measurements. Estimates suggest that this is due to turbulence in the baroclinic vorticity generated at the boundary of the oval. The reflectivity of a vortex itself (mechanism A) is comparable to that of the oval (mechanism C) but cuts-off at frequencies lower than those considered in all the experiments to date. The RCS of a vortex happens to peak at the frequency (about 49 MHz) where atmospheric radars (known as ST radars) operate and so the present prediction could be verified in the future. Finally , we suggest that hot engine exhaust could increase RCE by 40 db and reveal vortex circulation, provided its mixing with the surroundings is prevented in the laminarising flow of the vortices.

Shariff, Karim↗