Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Forward sensitivity method”

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 37 records · Page 2

On the dual advantage of placing observations through forward sensitivity analysis

The four-dimensional variational data assimilation methodology for assimilating noisy observations into a deterministic model has been the workhorse in forecasting centres for over three decades. While this method provides a computationally efficient framework for dynamic data assimilation, it is largely silent on the important question concerning the minimum number and placement of observations. To answer this question, we demonstrate the dual advantage of placing the observations where the square of the sensitivity of the model solution with respect to the unknown control variables, called forward sensitivities, attains its maximum. By following this approach, we can force the observability Gramian to be of full rank, which in turn guarantees efficient recovery of the optimal values of the control variables, which is the first of the two advantages of this strategy. We further show that the proposed strategy of placing observations has another inherent optimality: the square of the sensitivity of the optimal estimates of the control with respect to the observations (used to obtain these estimates) attains its minimum value, a second advantage that is a direct consequence of the above strategy for placing observations. Our analytical framework and numerical experiments on linear and nonlinear systems confirm the effectiveness of our proposed strategy.

97 MATHEMATICS AND COMPUTING↗

PyApprox: Enabling efficient model analysis

PyApprox is a Python-based one-stop-shop for probabilistic analysis of scientific numerical models. Easy to use and extendable tools are provided for constructing surrogates, sensitivity analysis, Bayesian inference, experimental design, and forward uncertainty quantification. The algorithms implemented represent the most popular methods for model analysis developed over the past two decades, including recent advances in multi-fidelity approaches that use multiple model discretizations and/or simplified physics to significantly reduce the computational cost of various types of analyses. Simple interfaces are provided for the most commonly-used algorithms to limit a user’s need to tune the various hyper-parameters of each algorithm. However, more advanced work flows that require customization of hyper-parameters is also supported. An extensive set of Benchmarks from the literature is also provided to facilitate the easy comparison of different algorithms for a wide range of model analyses. This paper introduces PyApprox and its various features, and presents results demonstrating the utility of PyApprox on a benchmark problem modeling the advection of a tracer in ground water.

97 MATHEMATICS AND COMPUTING↗

Neutron diffraction discriminates between models for the nanoarchitecture of graphene sheets in glassy carbon

Glassy carbon is a chemically inert, biocompatible, disordered material with the graphene sheet as its basic building block. Structural characterisation techniques have so far been unable to provide definitive distinction between the many proposed models for its structure. Computer based simulation methods have made a step forward by predicting structure from interaction potentials, but the results are sensitive to the choice of the potential. Here we use the white neutron beam of the Spallation Neutron Source at Oak Ridge National Laboratory coupled with the SNAP time-of-flight diffractometer to sample the reciprocal space of glassy carbon and calculate accurate radial distribution functions. From the radial distribution information, we determine graphene sheet dimensions, registration, curvature, and assess defect content to rank the proposed structures in terms of their agreement with experimental data. In this work, we find that the recent models derived from empirical potentials give the best agreement.

36 MATERIALS SCIENCE↗

Forward Modeling the Orbits of Companions to Pulsating Stars from Their Light Travel Time Variations

Mutual gravitation between a pulsating star and an orbital companion leads to a time-dependent variation in path length for starlight traveling to Earth. These variations can be used for coherently pulsating stars, such as the δ Scuti variables, to constrain the masses and orbits of their companions. Observing these variations for δ Scuti stars has previously relied on subdividing the light curve and measuring the average pulsation phase in equally sized subdivisions, which leads to undersampling near periapsis. We introduce a new approach that simultaneously forward models each sample in the light curve and show that this method improves upon current sensitivity limits—especially in the case of highly eccentric and short-period binaries. We find that this approach is sensitive enough to observe Jupiter mass planets around δ Scuti stars under ideal conditions, and use gravity-mode pulsations in the subdwarf B star KIC 7668647 to detect its companion without radial velocity data. We further provide robust detection limits as a function of the signal-to-noise ratio of the pulsation mode and determine that the minimum detectable light travel time amplitude for a typical Kepler δ Scuti is around 2 s. This new method significantly enhances the application of light travel time variations to detecting short-period binaries with pulsating components, and pulsating A-type exoplanet host stars, especially as a tool for eliminating false positives.

79 ASTRONOMY AND ASTROPHYSICS↗

PyApprox: A software package for sensitivity analysis, Bayesian inference, optimal experimental design, and multi-fidelity uncertainty quantification and surrogate modeling

PyApprox is a Python-based one-stop-shop for probabilistic analysis of numerical models such as those used in the earth, environmental and engineering sciences. Easy to use and extendable tools are provided for constructing surrogates, sensitivity analysis, Bayesian inference, experimental design, and forward uncertainty quantification. The algorithms implemented represent a wide range of methods for model analysis developed over the past two decades, including recent advances in multi-fidelity approaches that use multiple model discretizations and/or simplified physics to significantly reduce the computational cost of various types of analyses. An extensive set of Benchmarks from the literature is also provided to facilitate the easy comparison of new or existing algorithms for a wide range of model analyses. Here, this paper introduces PyApprox and its various features, and presents results demonstrating the utility of PyApprox on a benchmark problem modeling the advection of a tracer in groundwater.

54 ENVIRONMENTAL SCIENCES↗

Three-dimensional trajectory analyses of two drop sizing instruments - PMS OAP and PMS FSSP

Flow-induced distortions of water drop fluxes and speeds as seen by the PMS optical array probe (OAP) and the PMS forward scattering spectrometer probe (FSSP) are estimated via three-dimensional flow and trajectory calculation methods. The sensitivities of the instruments to water drop diameter, angle of attack, and free stream air speed are determined. The instruments are first placed in isolation and then mounted under the wing of a Twin Otter airplane. For the wing-mounted OAP at 4-deg angle of attack, partial flow stagnation under the uptilted wing causes a significant decrease in both the flux and speed for small water drops. For the wing-mounted FSSP, sensitivity is found to both angle of attack and free stream air speed.

Norment, H. G.↗

Computational Response Theory for Dynamics

Quantifying the sensitivity - how a quantity of interest (QoI) varies with respect to a parameter – and response – the representation of a QoI as a function of a parameter - of a computer model of a parametric dynamical system is an important and challenging problem. Traditional methods fail in this context since sensitive dependence on initial conditions implies that the sensitivity and response of a QoI may be ill-conditioned or not well-defined. If a chaotic model has an ergodic attractor, then ergodic averages of QoIs are well-defined quantities and their sensitivity can be used to characterize model sensitivity. The response theorem gives sufficient conditions such that the local forward sensitivity – the derivative with respect to a given parameter - of an ergodic average of a QoI is well-defined. We describe a method based on ergodic and response theory for computing the sensitivity and response of a given QoI with respect to a given parameter in a chaotic model with an ergodic and hyperbolic attractor. This method does not require computation of ensembles of the model with perturbed parameter values. The method is demonstrated and some of the computations are validated on the Lorenz 63 and Lorenz 96 models.

97 MATHEMATICS AND COMPUTING↗

An easily implemented static condensation method for structural sensitivity analysis

A black-box approach to static condensation for sensitivity analysis is presented with illustrative examples of a cube and a car structure. The sensitivity of the structural response with respect to joint stiffness parameter is calculated using the direct method, forward-difference, and central-difference schemes. The efficiency of the various methods for identifying joint stiffness parameters from measured static deflections of these structures is compared. The results indicate that the use of static condensation can reduce computation times significantly and the black-box approach is only slightly less efficient than the standard implementation of static condensation. The ease of implementation of the black-box approach recommends it for use with general-purpose finite element codes that do not have a built-in facility for static condensation.

Gangadharan, S. N.↗

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗

On Study of Air/Space-borne Dual-Wavelength Radar for Estimates of Rain Profiles

In this study, a framework is discussed to apply air/space-borne dual-wavelength radar for the estimation of characteristic parameters of hydrometeors. The focus of our study is on the Global Precipitation Measurements (GPM) precipitation radar, a dual-wavelength radar that operates at Ku (13.8 GHz) and Ka (35 GHz) bands. As the droplet size distributions (DSD) of rain are expressed as the Gamma function, a procedure is described to derive the median volume diameter (D(sub 0)) and particle number concentration (N(sub T)) of rain. The correspondences of an important quantity of dual-wavelength radar, defined as deferential frequency ratio (DFR), to the D(sub 0) in the melting region are given as a function of the distance from the 0 C isotherm. A self-consistent iterative algorithm that shows a promising to account for rain attenuation of radar and infer the DSD without use of surface reference technique (SRT) is examined by applying it to the apparent radar reflectivity profiles simulated from the DSD model and then comparing the estimates with the model (true) results. For light to moderate rain the self-consistent rain profiling approach converges to unique and correct solutions only if the same shape factors of Gamma functions are used both to generate and retrieve the rain profiles, but does not converges to the true solutions if the DSD form is not chosen correctly. To further examine the dual-wavelength techniques, the self-consistent algorithm, along with forward and backward rain profiling algorithms, is then applied to the measurements taken from the 2nd generation Precipitation Radar (PR-2) built by Jet Propulsion Laboratory. It is found that rain profiles estimated from the forward and backward approaches are not sensitive to shape factor of DSD Gamma distribution, but the self-consistent method is.

Liao, Liang↗

Sensitivity of rotor blade vibration characteristics to torsional oscillations.

A theoretical investigation of dynamic response characteristics of helicopter rotor blades in forward flight was carried out with special emphasis on the torsional degrees-of-freedom. The finite element method was applied in the formulation of the coupled equations of motion for flapwise bending and torsion for blades with non-collinear elastic, mass and aerodynamic axes. The sensitivity of blade vibration characteristics with respect to structural, geometric and aerodynamic properties as well as flight conditions was evaluated. Numerical results for sample blades were plotted to show the variation of the coupling between bending and torsional components of the response.

Bratanow, T.↗

Stability sensitivity analysis for the aeroelastic optimization of a helicopter rotor

A sensitivity study of blade stability in forward flight for a hingeless rotor with respect to structural design variables is carried out using a direct analytical method. Structural design variables include nonstructural mass distribution (spanwise and chordwise), chordwise offset of center of gravity, and blade bending stiffnesses (flap, lag and torsion). The formulation for blade steady response is based on a finite element method in space and time. The vehicle trim and blade steady response are calculated iteratively as one coupled solution using a modified Newton method. Eigenvalues corresponding to different blade modes are calculated using Floquet transition matrix theory. The formulation for derivatives of the eigenvalues with respect to design variables is implemented using a direct analytical approach (chain rule differentiation), and constitutes an integral part of the regular stability analysis. The stability sensitivity derivatives were obtained at a fraction of computation time compared to the frequently adopted finite difference method. A parametric study showed that nonstructural mass and chordwise cg offset of outboarad elements, and lag bending stiffness of inboard elements, have powerful influence on blade stability.

Lim, Joon W.↗

PMA-PhyloChip DNA Microarray to Elucidate Viable Microbial Community Structure

Since the Viking missions in the mid-1970s, traditional culture-based methods have been used for microbial enumeration by various NASA programs. Viable microbes are of particular concern for spacecraft cleanliness, for forward contamination of extraterrestrial bodies (proliferation of microbes), and for crew health/safety (viable pathogenic microbes). However, a "true" estimation of viable microbial population and differentiation from their dead cells using the most sensitive molecular methods is a challenge, because of the stability of DNA from dead cells. The goal of this research is to evaluate a rapid and sensitive microbial detection concept that will selectively estimate viable microbes. Nucleic acid amplification approaches such as the polymerase chain reaction (PCR) have shown promise for reducing time to detection for a wide range of applications. The proposed method is based on the use of a fluorescent DNA intercalating agent, propidium monoazide (PMA), which can only penetrate the membrane of dead cells. The PMA-quenched reaction mixtures can be screened, where only the DNA from live cells will be available for subsequent PCR reaction and microarray detection, and be identified as part of the viable microbial community. An additional advantage of the proposed rapid method is that it will detect viable microbes and differentiate from dead cells in only a few hours, as opposed to less comprehensive culture-based assays, which take days to complete. This novel combination approach is called the PMA-Microarray method. DNA intercalating agents such as PMA have previously been used to selectively distinguish between viable and dead bacterial cells. Once in the cell, the dye intercalates with the DNA and, upon photolysis under visible light, produces stable DNA adducts. DNA cross-linked in this way is unavailable for PCR. Environmental samples suspected of containing a mixture of live and dead microbial cells/spores will be treated with PMA, and then incubated in the dark. Thereafter, the sample is exposed to visible light for five minutes, so that the DNA from dead cells will be cross-linked. Following this PMA treatment step, the sample is concentrated by centrifugation and washed (to remove excessive PMA) before DNA is extracted. The 16S rRNA gene fragments will be amplified by PCR to screen the total microbial community using PhyloChip DNA microarray analysis. This approach will detect only the viable microbial community since the PMA intercalated DNA from dead cells would be unavailable for PCR amplification. The total detection time including PCR reaction for low biomass samples will be a few hours. Numerous markets may use this technology. The food industry uses spore detection to validate new alternative food processing technologies, sterility, and quality. Pharmaceutical and medical equipment companies also detect spores as a marker for sterility. This system can be used for validating sterilization processes, water treatment systems, and in various public health and homeland security applications.

Venkateswaran, Kasthuri J.↗

Retrieving the Vertical Structure of the Effective Aerosol Complex Index of Refraction from a Combination of Aerosol in Situ and Remote Sensing Measurements During TARFOX

The largest uncertainty in estimates of the effects of atmospheric aerosols on climate stems from uncertainties in the determination of their microphysical properties, including the aerosol complex index of refraction, which in turn determines their optical properties. A novel technique is used to estimate the aerosol complex index of refraction in distinct vertical layers from a combination of aerosol in situ size distribution and remote sensing measurements during the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX). In particular, aerosol backscatter measurements using the NASA Langley LASE (Lidar Atmospheric Sensing Experiment) instrument and in situ aerosol size distribution data are utilized to derive vertical profiles of the "effective" aerosol complex index of refraction at 815 nm (i.e., the refractive index that would provide the same backscatter signal in a forward calculation on the basis of the measured in situ particle size distributions for homogeneous, spherical aerosols). A sensitivity study shows that this method yields small errors in the retrieved aerosol refractive indices, provided the errors in the lidar-derived aerosol backscatter are less than 30% and random in nature. Absolute errors in the estimated aerosol refractive indices are generally less than 0.04 for the real part and can be as much as 0.042 for the imaginary part in the case of a 30% error in the lidar-derived aerosol backscatter. The measurements of aerosol optical depth from the NASA Ames Airborne Tracking Sunphotometer (AATS-6) are successfully incorporated into the new technique and help constrain the retrieved aerosol refractive indices. An application of the technique to two TARFOX case studies yields the occurrence of vertical layers of distinct aerosol refractive indices. Values of the estimated complex aerosol refractive index range from 1.33 to 1.45 for the real part and 0.001 to 0.008 for the imaginary part. The methodology devised in this study provides, for the first time, a complete set of vertically resolved aerosol size distribution and refractive index data. yielding the vertical distribution of aerosol optical properties required for the determination of aerosol-induced radiative flux changes.

Redemann, J.↗

Retrieving the Vertical Structure of the Effective Aerosol Complex Index of Refraction from a Combination of Aerosol in Situ and Remote Sensing Measurements During TARFOX

The largest uncertainty in estimates of the effects of atmospheric aerosols on climate stems from uncertainties in the determination of their microphysical properties, including the aerosol complex index of refraction, which in turn determines their optical properties. A novel technique is used to estimate the aerosol complex index of refraction in distinct vertical layers from a combination of aerosol in situ size distribution and remote sensing measurements during the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX). In particular, aerosol backscatter measurements using the NASA Langley LASE (Lidar Atmospheric Sensing Experiment) instrument and in situ aerosol size distribution data are utilized to derive vertical profiles of the 'effective' aerosol complex index of refraction at 815 nm (i.e., the refractive index that would provide the same backscatter signal in a forward calculation on the basis of the measured in situ particle size distributions for homogeneous, spherical aerosols). A sensitivity study shows that this method yields small errors in the retrieved aerosol refractive indices, provided the errors in the lidar derived aerosol backscatter are less than 30% and random in nature. Absolute errors in the estimated aerosol refractive indices are generally less than 0.04 for the real part and can be as much as 0.042 for the imaginary part in the case of a 30% error in the lidar-derived aerosol backscatter. The measurements of aerosol optical depth from the NASA Ames Airborne Tracking Sunphotometer (AATS-6) are successfully incorporated into the new technique and help constrain the retrieved aerosol refractive indices. An application of the technique to two TARFOX case studies yields the occurrence of vertical layers of distinct aerosol refractive indices. Values of the estimated complex aerosol refractive index range from 1.33 to 1.45 for the real part and 0.001 to 0.008 for the imaginary part. The methodology devised in this study provides, for the first time a complete set of vertically resolved aerosol size distribution and refractive index data, yielding the vertical distribution of aerosol optical properties required for the determination of aersol-induced radiative flux changes

Redemann, J.↗

Integrated aerodynamic-structural-control wing design

The aerodynamic-structural-control design of a forward-swept composite wing for a high subsonic transport aircraft is considered. The structural analysis is based on a finite-element method. The aerodynamic calculations are based on a vortex-lattice method, and the control calculations are based on an output feedback control. The wing is designed for minimum weight subject to structural, performance/aerodynamic and control constraints. Efficient methods are used to calculate the control-deflection and control-effectiveness sensitivities which appear as second-order derivatives in the control constraint equations. To suppress the aeroelastic divergence of the forward-swept wing, and to reduce the gross weight of the design aircraft, two separate cases are studied: (1) combined application of aeroelastic tailoring and active controls; and (2) aeroelastic tailoring alone. The results of this study indicated that, for this particular example, aeroelastic tailoring is sufficient for suppressing the aeroelastic divergence, and the use of active controls was not necessary.

Rais-Rohani, M.↗

NASA/American Cancer Society High-Resolution Flow Cytometry Project-I

BACKGROUND: The NASA/American Cancer Society (ACS) flow cytometer can simultaneously analyze the electronic nuclear volume (ENV) and DNA content of cells. This study describes the schematics, resolution, reproducibility, and sensitivity of biological standards analyzed on this unit. METHODS: Calibrated beads and biological standards (lymphocytes, trout erythrocytes [TRBC], calf thymocytes, and tumor cells) were analyzed for ENV versus DNA content. Parallel data (forward scatter versus DNA) from a conventional flow cytometer were obtained. RESULTS: ENV linearity studies yielded an R value of 0.999. TRBC had a coefficient of variation (CV) of 1.18 +/- 0.13. DNA indexes as low as 1.02 were detectable. DNA content of lymphocytes from 42 females was 1.9% greater than that for 60 males, with a noninstrumental variability in total DNA content of 0.5%. The ENV/DNA ratio was constant in 15 normal human tissue samples, but differed in the four animal species tested. The ENV/DNA ratio for a hypodiploid breast carcinoma was 2.3 times greater than that for normal breast tissue. CONCLUSIONS: The high-resolution ENV versus DNA analyses are highly reliable, sensitive, and can be used for the detection of near-diploid tumor cells that are difficult to identify with conventional cytometers. ENV/DNA ratio may be a useful parameter for detection of aneuploid populations.

NASA Center JSC↗

Some Advanced Concepts in Discrete Aerodynamic Sensitivity Analysis

An efficient incremental iterative approach for differentiating advanced flow codes is successfully demonstrated on a two-dimensional inviscid model problem. The method employs the reverse-mode capability of the automatic differentiation software tool ADIFOR 3.0 and is proven to yield accurate first-order aerodynamic sensitivity derivatives. A substantial reduction in CPU time and computer memory is demonstrated in comparison with results from a straightforward, black-box reverse-mode applicaiton of ADIFOR 3.0 to the same flow code. An ADIFOR-assisted procedure for accurate second-rder aerodynamic sensitivity derivatives is successfully verified on an inviscid transonic lifting airfoil example problem. The method requires that first-order derivatives are calculated first using both the forward (direct) and reverse (adjoinct) procedures; then, a very efficient noniterative calculation of all second-order derivatives can be accomplished. Accurate second derivatives (i.e., the complete Hesian matrices) of lift, wave drag, and pitching-moment coefficients are calculated with respect to geometric shape, angle of attack, and freestream Mach number.

Taylor, Arthur C., III↗