Search NASA⌕ Search

SEARCH · Search NASA

Results for “sensitivity gradient”

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 199 records · Page 11

Multispectral processing based on groups of resolution elements

Several nine-point rules are defined and compared with previously studied rules. One of the rules performed well in boundary areas, but with reduced efficiency in field interiors; another combined best performance on field interiors with good sensitivity to boundary detail. The basic threshold gradient and some modifications were investigated as a means of boundary point detection. The hypothesis testing methods of closed-boundary formation were also tested and evaluated. An analysis of the boundary detection problem was initiated, employing statistical signal detection and parameter estimation techniques to analyze various formulations of the problem. These formulations permit the atmospheric and sensor system effects on the data to be thoroughly analyzed. Various boundary features and necessary assumptions can also be investigated in this manner.

Richardson, W.↗

Simulation studies of wide and medium field of view earth radiation data analysis

A parameter estimation technique is presented to estimate the radiative flux distribution over the earth from radiometer measurements at satellite altitude. The technique analyzes measurements from a wide field of view (WFOV), horizon to horizon, nadir pointing sensor with a mathematical technique to derive the radiative flux estimates at the top of the atmosphere for resolution elements smaller than the sensor field of view. A computer simulation of the data analysis technique is presented for both earth-emitted and reflected radiation. Zonal resolutions are considered as well as the global integration of plane flux. An estimate of the equator-to-pole gradient is obtained from the zonal estimates. Sensitivity studies of the derived flux distribution to directional model errors are also presented. In addition to the WFOV results, medium field of view results are presented.

Green, R. N.↗

Approach for Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

This paper presents an implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for a quasi 1-D Euler CFD (computational fluid dynamics) code. Given uncertainties in statistically independent, random, normally distributed input variables, a first- and second-order statistical moment matching procedure is performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, the moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Aerodynamic Design Optimization on Unstructured Meshes Using the Navier-Stokes Equations

A discrete adjoint method is developed and demonstrated for aerodynamic design optimization on unstructured grids. The governing equations are the three-dimensional Reynolds-averaged Navier-Stokes equations coupled with a one-equation turbulence model. A discussion of the numerical implementation of the flow and adjoint equations is presented. Both compressible and incompressible solvers are differentiated and the accuracy of the sensitivity derivatives is verified by comparing with gradients obtained using finite differences. Several simplifying approximations to the complete linearization of the residual are also presented, and the resulting accuracy of the derivatives is examined. Demonstration optimizations for both compressible and incompressible flows are given.

Nielsen, Eric J.↗

Approach for Input Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

An implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for quasi 3-D Euler CFD code is presented. Given uncertainties in statistically independent, random, normally distributed input variables, first- and second-order statistical moment procedures are performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, these moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗

Effects of Rangeland Evapotranspiration on Groundwater Recharge, Discharge, and Flow at the Tuba City, Arizona, Disposal Site - 17159

The US Department of Energy Office of Legacy Management is modeling groundwater flow and contaminant transport at a former uranium mill site near Tuba City, Arizona. A goal is to project groundwater travel times and flow volumes between the former mill site and Moenkopi Wash, a down-gradient stream. Aquifer recharge and discharge are sensitive model parameters; however, assigning representative rates is inexact and involves approximations. We applied a remote sensing algorithm to determine spatially variable evapotranspiration (ET), precipitation recharge, and groundwater discharge rates for use in a large scale groundwater flow model for the site. ET is the combination of evaporation and plant transpiration from soil and groundwater. We estimated landscape-scale ET over a 13-year period (2000 to 2012) for distinct plant communities within a 3531-hectare groundwater model domain (GMD), and then evaluated effects of ET on groundwater recharge and discharge within the GMD. Our empirical algorithm was derived from ground ET measurements, multispectral satellite imagery, and temperature data. Groundwater recharge or discharge rates were calculated for each plant community (or ET Zone) as the difference between precipitation (PPT) and ET rates. Recharge occurred in plant communities where PPT exceeded ET; discharge occurred where ET exceeded PPT. Estimates of groundwater recharge rates for upland plant communities that survive on meteoric water ranged from 0 to 88 mm yr{sup -1}. Estimates of groundwater discharge for phreatophyte plant communities ranged from 23 mm yr{sup -1} to 150 mm yr{sup -1}. Phreatophytes are plants that survive by extracting groundwater. Zero recharge occurred in plant communities where ET rates equaled annual average precipitation. An increase in groundwater recharge (PPT > ET) was associated with past land disturbances and heavy livestock grazing in upland areas. Groundwater discharge (ET > PPT) was highest in riparian phreatophyte communities, but lower than optimal in upland phreatophyte communities impacted by heavy grazing. Protecting phreatophyte communities from grazing at Tuba City could potentially increase ET from 153 mm yr{sup -1} to at least 500 mm yr{sup -1}. The validity of using the ET algorithm to estimate recharge and discharge was evaluated by comparing the net volumetric outflow rate for the model domain with base flow measurements collected from Moenkopi Wash. Modeled outflow approximated measured flow gains in the Wash, increasing our confidence in the ET algorithm and in the model. Results suggest that rangeland management practices that reduce groundwater recharge rates and increase groundwater discharge rates should be evaluated as part of an overall groundwater remediation strategy. (authors)

54 ENVIRONMENTAL SCIENCES↗

Scale selection and feedback loops for patterns in drift wave-zonal flow turbulence

The scale selection and feedback loops for the formation and sustainment of a mesoscopic staircase profile structure are investigated for drift wave-zonal flow turbulence. A mean field model derived from the Hasegawa-Wakatani system and including the evolution of mean density, mean vorticity and perturbed potential enstrophy (PE), is used. Here, it is found that a quasi-periodic zonal staircase forms from self-sharpening of modulation. The principle feedback loop is through the nonlinear dependence of mixing length on electron density gradient, which enters by way of the potential vorticity (PV) gradient. Counterintuitively, E×B shearing is not effective. Moreover, the number of steps in the staircase is sensitive to both the drive (production rate of PE and initial density gradient) and damping (flow viscosity and collisional diffusivity) factors. The minimal step scale is selected by competition between the initial density gradient and diffusive dissipation. Finite turbulence spreading is necessary to form the staircase, but moderate enhancement of turbulence spreading tends to wash out the pattern. The staircase retains a memory of its initial state. Both the mean E×B shear and zonal shear affect the staircase evolution. A strong mean shear quenches the pattern by suppressing the drift wave turbulence. Lastly, the implications of these findings are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved Cloud Condensation Nucleus Spectrometer

An improved thermal-gradient cloud condensation nucleus spectrometer (CCNS) has been designed to provide several enhancements over prior thermal- gradient counters, including fast response and high-sensitivity detection covering a wide range of supersaturations. CCNSs are used in laboratory research on the relationships among aerosols, supersaturation of air, and the formation of clouds. The operational characteristics of prior counters are such that it takes long times to determine aerosol critical supersaturations. Hence, there is a need for a CCNS capable of rapid scanning through a wide range of supersaturations. The present improved CCNS satisfies this need. The improved thermal-gradient CCNS (see Figure 1) incorporates the following notable features: a) The main chamber is bounded on the top and bottom by parallel thick copper plates, which are joined by a thermally conductive vertical wall on one side and a thermally nonconductive wall on the opposite side. b) To establish a temperature gradient needed to establish a supersaturation gradient, water at two different regulated temperatures is pumped through tubes along the edges of the copper plates at the thermally-nonconductive-wall side. Figure 2 presents an example of temperature and supersaturation gradients for one combination of regulated temperatures at the thermally-nonconductive-wall edges of the copper plates. c) To enable measurement of the temperature gradient, ten thermocouples are cemented to the external surfaces of the copper plates (five on the top plate and five on the bottom plate), spaced at equal intervals along the width axis of the main chamber near the outlet end. d) Pieces of filter paper or cotton felt are cemented onto the interior surfaces of the copper plates and, prior to each experimental run, are saturated with water to establish a supersaturation field inside the main chamber. e) A flow of monodisperse aerosol and a dilution flow of humid air are introduced into the main chamber at the inlet end. The inlet assembly is designed to offer improved (relative to prior such assemblies) laminar-flow performance within the main chamber. Dry aerosols are subjected to activation and growth in the supersaturation field. f) After aerosol activation, at the outlet end of the main chamber, a polished stainless-steel probe is used to sample droplets into a laser particle counter. The probe features an improved design for efficient sampling. The counter has six channels with size bins in the range of 0.5- to 5.0-micron diameter. g) To enable efficient sampling, the probe is scanned along the width axis of the main chamber (thereby effecting scanning along the temperature gradient and thereby, further, effecting scanning along the supersaturation gradient) by means of a computer-controlled translation stage.

Leu, Ming-Taun↗

A molecular dynamics simulation study of thermal conductivity anisotropy in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX)

Molecular dynamics (MD) simulations were used to predict the thermal conductivity of β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) along directions normal to the (011), (110), and (010) crystal planes. These directions were selected based on the measured morphological importance of the corresponding crystal surfaces. A reverse non-equilibrium MD approach was used wherein a constant heat flux is imposed along a prescribed direction and the resulting steady-state temperature gradient determined. The coefficient of thermal conductivity λ is the quotient of heat flux and temperature gradient (i.e. Fourier's law). Finite-size effects and sensitivity to imposed heat flux were investigated. The results reveal a modest dependence of the conductivity on crystal orientation, significant finite-size effects, and low sensitivity to imposed flux so long as the Fourier's law analysis is limited to the spatial interval in the simulation cell for which the temperature gradient is constant. Infinite-length thermal conductivities were estimated by linear regression of λ -1 (L) versus reciprocal cell length L -1 for each direction. The predicted values are systematically larger, but within a factor of two, than most published experimental determinations, the latter of which were obtained for pressed-powder or composite samples rather than oriented single crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vibration isolation technology: Sensitivity of selected classes of experiments to residual accelerations

The solution was sought of a 2-D axisymmetric moving boundary problem for the sensitivity of isothermal and nonisothermal liquid columns and the sensitivity of thermo-capillary flows to buoyancy driven convection caused by residual accelerations. The sensitivity of a variety of space experiments to residual accelerations are examined. In all the cases discussed, the sensitivity is related to the dynamic response of a fluid. In some cases the sensitivity can be defined by the magnitude of the response of the velocity field. This response may involve motion of the fluid associated with internal density gradients, or the motion of a free liquid surface. For fluids with internal density gradients, the type of acceleration to which the experiment is sensitive will depend on whether buoyancy driven convection must be small in comparison to other types of fluid motion (such as thermocapillary flow), or fluid motion must be suppressed or eliminated (such as in diffusion studies, or directional solidification experiments). The effect of the velocity on the composition and temperature field must be considered, particularly in the vicinity of the melt crystal interface. As far as the response to transient disturbances is concerned the sensitivity is determined by both the magnitude and frequency the acceleration and the characteristic momentum and solute diffusion times.

Alexander, J. Iwan D.↗

Process modelling for Space Station experiments

Examined here is the sensitivity of a variety of space experiments to residual accelerations. In all the cases discussed the sensitivity is related to the dynamic response of a fluid. In some cases the sensitivity can be defined by the magnitude of the response of the velocity field. This response may involve motion of the fluid associated with internal density gradients, or the motion of a free liquid surface. For fluids with internal density gradients, the type of acceleration to which the experiment is sensitive will depend on whether buoyancy driven convection must be small in comparison to other types of fluid motion, or fluid motion must be suppressed or eliminated. In the latter case, the experiments are sensitive to steady and low frequency accelerations. For experiments such as the directional solidification of melts with two or more components, determination of the velocity response alone is insufficient to assess the sensitivity. The effect of the velocity on the composition and temperature field must be considered, particularly in the vicinity of the melt-crystal interface. As far as the response to transient disturbances is concerned, the sensitivity is determined by both the magnitude and frequency of the acceleration and the characteristic momentum and solute diffusion times. The microgravity environment, a numerical analysis of low gravity tolerance of the Bridgman-Stockbarger technique, and modeling crystal growth by physical vapor transport in closed ampoules are discussed.

Alexander, J. Iwan D.↗

Crack tip microplasticity mediated by microstructure gradients

Abstract Traditional fracture theories infer damage at cracks (local field) by surveying loading conditions away from cracks (far field). This approach has been successful in predicting ductile fracture, but it normally assumes isotropic and homogeneous materials. However, myriads of manufacturing procedures induce heterogeneous microstructural gradients that can affect the accuracy of traditional fracture models. This work presents a microstructure‐sensitive finite element approach to explore the shielding effects of grain size and crystallographic orientation gradients on crack tip microplasticity and blunting. A dislocation density‐based crystal plasticity model conveys texture evolution, grain size effects, and directional hardening by computing the constraint from dislocation structures. The results demonstrate that the microstructure can act as a buffer between the local and far fields that affects the crack tip microplasticity variability. For nominal opening loading, grain size and texture affect the local ductility and induce a non‐negligible multiaxial plastic deformation. Furthermore, driving forces based on measuring displacements away from the crack tip are less affected by the microstructure, which suggests that traditional experimental methods smear out important crack tip variability.

Castelluccio, Gustavo M.↗

Optimal control concepts in design sensitivity analysis

A close link is established between open loop optimal control theory and optimal design by noting certain similarities in the gradient calculations. The resulting benefits include a unified approach, together with physical insights in design sensitivity analysis, and an efficient approach for simultaneous optimal control and design. Both matrix displacement and matrix force methods are considered, and results are presented for dynamic systems, structures, and elasticity problems.

Belegundu, Ashok D.↗

Late quaternary time series of Arabian Sea productivity: Global and regional signals

Modern annual floral and faunal production in the northwest Arabian Sea derives primarily from upwelling induced by strong southwest winds during June, July, and August. Indian Ocean summer monsoon winds are, in turn, driven by differential heating between the Asian continent and the Indian ocean to the south. This differential heating produces a strong pressure gradient resulting in southwest monsoon winds and both coastal and divergent upwelling off the Arabian Peninsula. Over geologic time scales (10(exp 4) to 10(exp 6) years), monsoon wind strength is sensitive to changes in boundary conditions which influence this pressure gradient. Important boundary conditions include the seasonal distribution of solar radiation, global ice volume, Indian Ocean sea surface temperature, and the elevation and albedo of the Asian continent. To the extent that these factors influence monsoon wind strength, they also influence upwelling and productivity. In addition, however, productivity associated with upwelling can be decoupled from the strength of the summer monsoon winds via ocean mechanisms which serve to inhibit or enhance the nutrient supply in the intermediate waters of the Indian Ocean, the source for upwelled waters in the Arabian Sea. To differentiate productivity associated with wind-induced upwelling from that associated with other components of the system such as nutrient sequestering in glacial-age deep waters, we employ a strategy which monitors independent components of the oceanic and atmospheric subsystems. Using sediment records from the Owen Ridge, northwest Arabian Sea, we monitor the strength of upwelling and productivity using two independent indicators, percent G. bulloides and opal accumulation. We monitor the strength of southwest monsoon winds by measuring the grain-size of lithogenic dust particles blown into the Arabian Sea from the surrounding deserts of the Somali and Arabian Peninsulas. Our current hypothesis is that the variability associated with the 41 kyr power in the G. bulloides and opal accumulation records derive from nutrient availability in the intermediate waters which are upwelled via monsoon winds. This hypothesis is testable by comparison with Cd records of intermediate and deep waters of the Atlantic and Indian Ocean.

Clemens, Steven C.↗

Propagation of input parameter uncertainties in transport models

Here, the many sources of uncertainty in validation studies of plasma turbulence in magnetically confined fusion devices are well-known. In this paper, we investigate how to efficiently transform uncertainties in experimentally derived transport model inputs into model prediction uncertainties, using the quasilinear trapped-gyro-Landau-fluid (TGLF) turbulent transport model. We use the rapidly converging and computationally inexpensive non-intrusive probabilistic collocation method (PCM) to propagate input parameter uncertainty probability distribution functions (PDFs) through TGLF, yielding PDFs of predicted transport fluxes. We observe in many cases that the flux PDFs exhibit significant non-normal features such as strong skewness, even when the input distributions were normal. To illustrate the utility of the PCM approach, we apply this methodology to transport predictions for a DIII-D ITER baseline plasma in which the mix of neutral beam injection (NBI) and electron cyclotron heating (ECH) was varied. The model predictions show clear changes in the parametric dependencies and sensitivities of the turbulence between the two heating mixes. Specifically, when only NBI heating was used, the transport fluxes responded significantly only to the ion temperature gradient scale length. However, when both NBI and ECH were applied, the electron transport channels demonstrate a strong sensitivity to the electron temperature and density gradients not observed in the NBI-only case. Additional context for the PCM approach is provided by comparing its predictions with those obtained via a local flux-matching approach. A new set of validation metrics based on the Wasserstein distance is proposed for PDF-based comparisons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Environmental controls on simulated deep moist convection initiation occurring during RELAMPAGO-CACTI

This study synthesizes the results of 13 high-resolution simulations of deep convective updrafts forming over idealized terrain using environments observed during the RELAMPAGO and CACTI field projects. Using composite soundings from multiple observed cases, and variations upon them, we explore the sensitivity of updraft properties (e.g., size, buoyancy, and vertical pressure gradient forces) to influences of environmental relative humidity, wind shear, and mesoscale orographic forcing that support or suppress deep convection initiation (CI). Emphasis is placed on differentiating physical processes affecting the development of updrafts (e.g., entrainment-driven dilution of updrafts) in environments typifying observed successful and null (i.e., no CI despite affirmative operational forecasts) CI events. Thermally-induced mesoscale orographic lift favors the production of deep updrafts originating from ~1–2-km-wide boundary layer thermals. Simulations without terrain forcing required much larger ( ~5-km-wide) thermals to yield precipitating convection. CI outcome was quite sensitive to environmental relative humidity; updrafts with increased buoyancy, depth, and intensity thrived in otherwise inhospitable environments by simply increasing the free tropospheric relative humidity. This implicates the entrainment of free-tropospheric air into updrafts as a prominent governor of CI, consistent with previous studies. Sensitivity of CI to the environmental wind is manifested by: 1) low-level flow affecting the strength and depth of mesoscale convergence along the terrain, and 2) clouds encountering updraft-suppressing pressure gradient forces while interacting with vertical wind shear in the free-troposphere. Among the ensemble of thermals occurring in each simulation, the widest deep updrafts in each simulation were the most sensitive to environmental influences.

54 ENVIRONMENTAL SCIENCES↗

Highly Rectifying Fluidic Diodes Based on Asymmetric Layer-by-Layer Nanofilms on Nanochannel Membranes

Nanochannel-based fluidic diodes display ion selectivity and ion-current rectification (ICR), which may prove important in energy harvesting and biosensors. This paper reports asymmetric functionalization of the outer surface of a flexible nanochannel polymer membrane to create fluidic diodes that give ICR. Layer-by-layer (LbL) adsorption with cross-linking of only the underlying part of the polyelectrolyte nanofilm leads to a porosity step across the film. The combination of a high effective surface-charge density and the porosity step in the film leads to a remarkable maximum ICR factor of ~200 with a pH gradient across the film. Incorporation of pH-sensitive polyelectrolyte components enables the ICR factor to increase an order of magnitude on going from pH 8 to pH 3. Furthermore, the coated membrane shows excellent anion-selectivity. Thus, LbL adsorption with partial cross-linking provides a simple method for creating coated nanochannel membranes that serve as pH-responsive ionic diodes for potential chemical/biosensors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗