Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 811 records · Page 45

Reliability Sensitivity Analysis and Design Optimization of Composite Structures Based on Response Surface Methodology

This report discusses the development and application of two alternative strategies in the form of global and sequential local response surface (RS) techniques for the solution of reliability-based optimization (RBO) problems. The problem of a thin-walled composite circular cylinder under axial buckling instability is used as a demonstrative example. In this case, the global technique uses a single second-order RS model to estimate the axial buckling load over the entire feasible design space (FDS) whereas the local technique uses multiple first-order RS models with each applied to a small subregion of FDS. Alternative methods for the calculation of unknown coefficients in each RS model are explored prior to the solution of the optimization problem. The example RBO problem is formulated as a function of 23 uncorrelated random variables that include material properties, thickness and orientation angle of each ply, cylinder diameter and length, as well as the applied load. The mean values of the 8 ply thicknesses are treated as independent design variables. While the coefficients of variation of all random variables are held fixed, the standard deviations of ply thicknesses can vary during the optimization process as a result of changes in the design variables. The structural reliability analysis is based on the first-order reliability method with reliability index treated as the design constraint. In addition to the probabilistic sensitivity analysis of reliability index, the results of the RBO problem are presented for different combinations of cylinder length and diameter and laminate ply patterns. The two strategies are found to produce similar results in terms of accuracy with the sequential local RS technique having a considerably better computational efficiency.

Rais-Rohani, Masoud↗

Daily estimates of the earth's pole position with the Global Positioning System

Daily estimates of the earth's pole position have been obtained with measurements from a worldwide network of GPS receivers, obtained during the three week GIG '91 experiment in January-February 1991. For this short-term study, the GPS based polar motion series agrees with the other space based geodetic techniques (Very Long Baseline Interferometry and Satellite Laser Ranging) to about 0.4 mas rms, after the removal of mean biases of order 1-3 mas. The small error in day-to-day variability is not sensitive to the fiducial strategy used, nor are fiducial sites even necessary for monitoring high frequency pole position variability. The small biases indicate that the applied reference frames of the three geodetic techniques are nearly aligned, that the GPS fiducial errors are small, and that systematic errors in GPS are also small (of order 5 ppb). A well determined reference frame is necessary for monitoring the long-term stability of polar motion and for separating it from other long-term signals such as tectonic motion and internal systematic errors.

Lindqwister, Ulf J.↗

Uncertainty and Sensitivity Analysis of TREAT Transient Test #15

The U.S. Department of Energy (DOE) is preparing for the resumption of transient testing using the Transient Reactor Test Facility (TREAT) located at the Idaho National Laboratory (INL). In order to support the analysis of TREAT operation and experiment planning, the MAMMOTH reactor physics team has recently begun development of uncertainty and sensitivity analysis for the TREAT core simulated via MAMMOTH. MAMMOTH is a high-fidelity three-dimensional transient multi-physics code developed at INL under the DOE-NE Advanced Modeling and Simulation program (NEAMS) [1]. This research is part of the development presenting some initial uncertainty and sensitivity analysis of TREAT Transient Test #15 [2]. In this summary, we employ the Sobol sampling strategy implemented in RAVEN [3] to efficiently perform the sensitivity and uncertainty analysis. The idea is using Sobol functional decomposition with generalized polynomial chaos (gPC) expansions to estimate the responses of interest [4]. The gPC coefficients are then computed via the spectral projection method with sparse grid approach. In addition, this strategy includes evaluation of the second or higher order sensitivity indices, allowing the study of interactions between several input variables. The principal focus of this analysis is to quantify the uncertainties in the total power and maximum fuel temperatures with respect to uncertainties in the initial transient conditions, such as rated power, fuel temperature, graphite fuel density and control rod boron concentration. In addition, both the first order and second order time-dependent sensitivity indices are computed and analyzed. The results show that the interactions between initial control rod boron concentration and the rated power will contribute at most 1% to the total uncertainty of the power peak.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainty and Sensitivity Analysis of TREAT Transient Test no. 15

The U.S. Department of Energy (DOE) is preparing for the resumption of transient testing using the Transient Reactor Test Facility (TREAT) located at the Idaho National Laboratory (INL). In order to support the analysis of TREAT operation and experiment planning, the MAMMOTH reactor physics team has recently begun development of uncertainty and sensitivity analysis for the TREAT core simulated via MAMMOTH. MAMMOTH is a high-fidelity three-dimensional transient multi-physics code developed at INL under the DOE-NE Advanced Modeling and Simulation program (NEAMS). This research is part of the development presenting some initial uncertainty and sensitivity analysis of TREAT Transient Test no. 15. In this summary, we employ the Sobol sampling strategy implemented in RAVEN to efficiently perform the sensitivity and uncertainty analysis. The idea is using Sobol functional decomposition with generalized polynomial chaos (gPC) expansions to estimate the responses of interest. The gPC coefficients are then computed via the sparse grid spectral projection method. In addition, this strategy includes evaluation of the second or higher order sensitivity indices, allowing the study of interactions between several input variables. The principal focus of this analysis is to quantify the uncertainties in the total power and maximum fuel temperatures with respect to uncertainties in the initial transient conditions, such as rated power, fuel temperature, graphite fuel density and control rod boron concentration. In addition, both the first order and second order time-dependent sensitivity indices are computed and analyzed. The results show that the interactions between initial control rod boron concentration and the rated power will contribute at most 1% to the total uncertainty of the power peak.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Low-symmetry polymorph of GaP upends bonding paradigms of metallic high-pressure III–V compounds

Here, the pressure-induced polymorphism of binary octet compounds has long been considered a settled problem although the possible atomic disordering of some phases remains a puzzling observation. Taking GaP as a case study, we conclude, through x-ray microdiffraction and first-principles calculations, that its high-pressure metallic phase (previously reported as being disordered) adopts in fact an ordered base-centered monoclinic structure previously unknown in this class of compounds. The formation of layered patterns with variable degrees of interlayer dimerization, as observed in GaP, marks a paradigm shift of our understanding of ordering in octet high-pressure phases which calls for a more extensive re-examination. A rich polymorphism with fine tuning of chemical and physical properties can be envisioned.

36 MATERIALS SCIENCE↗

Finding simplicity: unsupervised discovery of features, patterns, and order parameters via shift-invariant variational autoencoders *

Abstract Recent advances in scanning tunneling and transmission electron microscopies (STM and STEM) have allowed routine generation of large volumes of imaging data containing information on the structure and functionality of materials. The experimental data sets contain signatures of long-range phenomena such as physical order parameter fields, polarization, and strain gradients in STEM, or standing electronic waves and carrier-mediated exchange interactions in STM, all superimposed onto scanning system distortions and gradual changes of contrast due to drift and/or mis-tilt effects. Correspondingly, while the human eye can readily identify certain patterns in the images such as lattice periodicities, repeating structural elements, or microstructures, their automatic extraction and classification are highly non-trivial and universal pathways to accomplish such analyses are absent. We pose that the most distinctive elements of the patterns observed in STM and (S)TEM images are similarity and (almost-) periodicity, behaviors stemming directly from the parsimony of elementary atomic structures, superimposed on the gradual changes reflective of order parameter distributions. However, the discovery of these elements via global Fourier methods is non-trivial due to variability and lack of ideal discrete translation symmetry. To address this problem, we explore the shift-invariant variational autoencoders (shift-VAEs) that allow disentangling characteristic repeating features in the images, their variations, and shifts that inevitably occur when randomly sampling the image space. Shift-VAEs balance the uncertainty in the position of the object of interest with the uncertainty in shape reconstruction. This approach is illustrated for model 1D data, and further extended to synthetic and experimental STM and STEM 2D data. We further introduce an approach for training shift-VAEs that allows finding the latent variables that comport to known physical behavior. In this specific case, the condition is that the latent variable maps should be smooth on the length scale of the atomic lattice (as expected for physical order parameters), but other conditions can be imposed. The opportunities and limitations of the shift VAE analysis for pattern discovery are elucidated.

97 MATHEMATICS AND COMPUTING↗

Examination of Mean Stress Calculation Approaches in Rock Mechanics

The mean stress, as a fundamental statistical property of a group of stress data, is essential for stress variability characterisation. However, currently in rock mechanics, the mean stress is customarily and erroneously calculated by separately averaging the principal stress magnitudes and orientations. In order to draw the attention of our community to the appropriate approach for stress variability characterisation, here we compare the customary scalar/vector mean with that obtained by the mathematically rigorous tensorial approach—the Euclidean mean. Calculation of mean stress using both a small group of actual in situ stress measurement results and a large group of simulated stress data (obtained using the combined finite– discrete element method, FEMDEM) demonstrates that the two approaches yield different results. Further investigation of these results shows that the scalar/vector approach may yield non-unique and non-orthogonal mean principal stresses, and these may deviate significantly from those of the Euclidean mean. Our calculations and comparisons reveal that the scalar/ vector approach is deficient because it processes the principal stress magnitudes and orientations separately as independent quantities and ignores the connection between them. Conversely, the tensorial approach appropriately averages the tensors that simultaneously carry not only the information of stress magnitudes and orientations, but also the inherent relations between them. Therefore, arbitrarily using scalar/vector mean stress of in situ stress measurements as input in further rock engineering analyses may yield significantly erroneous results. Finally, we advise that stress data should be statistically processed in a tensorial manner using tensors referred to a common Cartesian coordinate system.

58 GEOSCIENCES↗

Regulating Multiple Variables To Understand the Nucleation and Growth and Transformation of PbS Nanocrystal Superlattices

Nanocrystals (NCs) can self-assemble into ordered superlattices with collective properties, but the ability for controlling NC assembly remains poorly understandable toward achievement of desired superlattice. This paper regulates several key variables of PbS NC assembly (e.g., NC concentration and solubility, solvent type, evaporation rate, seed mediation and thermal treatment), and thoroughly exploits the nucleation and growth as well as subsequent superlattice transformation of NC assembles and underneath mechanisms. PbS NCs in toluene self-assemble into a single face-centered-cubic (fcc) and body-centered-cubic (bcc) superlattice, respectively, at concentrations ≤17.5 and ≥70 mg/mL, but an intermediate concentration between them causes the coexistence of the two superlattices. Differently, NCs in hexane or chloroform self-assemble into only a single bcc superlattice. Distinct controls of NC assembly in solvent with variable concentrations confirm the NC concentration/solubility mediated nucleation and growth of superlattice, in which an evaporation-induced local gradient of NC concentration causes simultaneous nucleation of the two superlattices. The observation for the dense packing planes of NCs in fast growing fcc rather than bcc reveals the difference of entropic driving forces responsible for the two distinct superlattices. Decelerating the solvent evaporation does not amend the superlattice symmetry, but improves the superlattice crystallinity. In addition to shrinking the superlattice volume, thermal treatment also transforms the bcc to an fcc superlattice at 175 °C. Through a seed-meditated growth, the concentration-dependent superlattice does not change lattice symmetry over the course of continuous growth, whereas the newly nucleated secondary small nuclei through a concentration change have relatively higher surface energy and quickly dissolve in solution, providing additional NC sources for the ripening of the primarily nucleated larger and stable seeds. Finally, the observations under multiple controls of assembly parameters not only provide insights into the nucleation and growth as well as transformation of various superlattice polymorphs but also lay foundation for controlled fabrication of desired superlattice with tailored property.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advancing river corridor science beyond disciplinary boundaries with an inductive approach to catalyse hypothesis generation

Abstract A unified conceptual framework for river corridors requires synthesis of diverse site‐, method‐ and discipline‐specific findings. The river research community has developed a substantial body of observations and process‐specific interpretations, but we are still lacking a comprehensive model to distill this knowledge into fundamental transferable concepts. We confront the challenge of how a discipline classically organized around the deductive model of systematically collecting of site‐, scale‐, and mechanism‐specific observations begins the process of synthesis. Machine learning is particularly well‐suited to inductive generation of hypotheses. In this study, we prototype an inductive approach to holistic synthesis of river corridor observations, using support vector machine regression to identify potential couplings or feedbacks that would not necessarily arise from classical approaches. This approach generated 672 relationships linking a suite of 157 variables each measured at 62 locations in a fifth order river network. Eighty four percent of these relationships have not been previously investigated, and representing potential (hypothetical) process connections. We document relationships consistent with current understanding including hydrologic exchange processes, microbial ecology, and the River Continuum Concept, supporting that the approach can identify meaningful relationships in the data. Moreover, we highlight examples of two novel research questions that stem from interpretation of inductively‐generated relationships. This study demonstrates the implementation of machine learning to sieve complex data sets and identify a small set of candidate relationships that warrant further study, including data types not commonly measured together. This structured approach complements traditional modes of inquiry, which are often limited by disciplinary perspectives and favour the careful pursuit of parsimony. Finally, we emphasize that this approach should be viewed as a complement to, rather than in place of, more traditional, deductive approaches to scientific discovery.

54 ENVIRONMENTAL SCIENCES↗

Watersheds and stream networks viewed longitudinally: Example insights from novel spatial portrayals of watershed characteristics

Abstract Longitudinal depictions of watershed structure and characteristics, including topography, stream networks, wetlands, ground water levels, and land use, can provide watershed knowledge and understanding unavailable from standard plan view maps. Three case studies provide examples of knowledge gained by applying longitudinal views of stream networks, watershed hydrologic behavior, and land use distributions. Longitudinal views of mountain stream networks show extreme variability in the slope‐area relationships of low Strahler order streams, large discontinuities in drainage area (large parts of drainage area space are absent in networks), and large variations in network curvature. Longitudinal views of a groundwater‐dominated headwater watershed increase the inference available from limited groundwater observations and clearly reveal how groundwater connections affect the permanence of surface water features and the distribution of vadose zone storage in the landscape. Plotting land uses longitudinally illuminates and allows a quantitative analysis of how land uses are distributed relative to topographic position. Viewing watersheds and stream networks longitudinally can provide new insights into watershed forms and processes and motivate new questions and research.

60 APPLIED LIFE SCIENCES↗

Why is it still too warm or cold in my house? Examining the relationships between energy efficient capital and household energy insecurity

Here, this paper examines the relationships between energy efficient (EE) capital technology and household energy insecurity in the United States. The theoretical model of these relationships employs household production theory to capture the demand for and production of household energy services, and a stochastic production frontier approach to describe how having access to and the usage intensity of EE capital technology could help alleviate inefficiency in the production of household energy services. A working hypothesis formulated from the theoretical model posits that having EE capital technology in the home will reduce the level of household energy insecurity experienced. The extent of energy insecurity experienced is inferred from an energy insecurity index value assigned to each household, generated via the application of a dichotomous Rasch model to questions contained in the 2015 Residential Energy Consumption Survey. Noting the potential simultaneous relationship that exists between a household having access to and the usage intensity of EE capital technology and the experience of being energy insecure, an instrumental variables approach was employed to estimate a series of ordered logit models. Results suggest access to EE capital technology in the form of Energy Star® appliances, Energy Star® windows, or a SMART thermostat does not reduce the probability of experiencing a greater level of energy insecurity. Nor does the usage intensity of EE capital. Thus, policy instruments designed to alleviate household energy insecurity may need to go beyond simply helping households obtain EE capital technology.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Do We Need to Parameterize Mesoscale Convective Organization to Mitigate the MJO-Mean State Trade-Off?

Modifications in cumulus parameterizations that improve Madden-Julian oscillation (MJO) simulation tend to degrade the mean state—known as the MJO-mean state trade-off. The impacts of parameterizing mesoscale convective organization on the relationship between the MJO and the mean state simulation fidelity are examined. A series of experiments are made with a general circulation model that parameterizes the mesoscale convective organization and simulates well both the MJO and the mean state. In the control simulation, a prognostic nondimensional variable (Ω) represents the degree of convective organization. In order to examine the effect of the parameterized convective organization, fixed Ω values are imposed in a series of constrained experiments. The fixed-Ω simulations show a negative relationship between the MJO and the mean state simulation fidelity. The control simulation is deviated from the negative relationship, suggesting the parameterized mesoscale convective organization helps general circulation models to mitigate the MJO-mean state trade-off.

58 GEOSCIENCES↗

Extended MHD modeling of tearing-driven magnetic relaxation

Discrete relaxation events in reversed-field pinch relevant configurations are investigated numerically with nonlinear extended magnetohydrodynamic modeling, including the Hall term in Ohm’s law and first-order ion finite Larmor radius effects. Our results show variability among relaxation events, where the Hall dynamo effect may help or impede the MHD dynamo effect in relaxing the parallel current density profile. The competitive behavior arises from multi-helicity conditions where the dominant magnetic fluctuation is relatively small. The resulting changes in parallel current density and parallel flow are aligned in the core, consistent with experimental observations. Analysis of simulation results also confirms that force density from fluctuation-induced Reynolds stress arises subsequent to the drive from fluctuation-induced Lorentz force density. Transport of momentum density is found to be dominated by the fluctuation-induced Maxwell stress over most of the cross section with viscous and gyroviscous contributions being large in the edge region. The findings resolve a discrepancy with respect to the relative orientation of current density and flow relaxation, which had not been realized or investigated in Ref. [King et. al. Phys. Plasmas 19, 055905 (2012)], where only the magnitude of flow relaxation is actually consistent with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Substitution Models of Protein Evolution with Selection on Enzymatic Activity

Abstract Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme–substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.

Ferreiro, David↗

Heralded Multiplexed High-Efficiency Cascaded Source of Dual-Rail Entangled Photon Pairs Using Spontaneous Parametric Down-Conversion

Deterministic sources of high-fidelity entangled qubit pairs encoded in the dual-rail photonic basis, i.e., presence of a single photon in one of two orthogonal modes, are a key enabling technology of many applications of quantum information processing, including high-rate, high-fidelity quantum communications over long distances. The most popular and mature sources of such photonic entanglement, e.g., those that leverage spontaneous parametric down-conversion (SPDC) or spontaneous four-wave mixing, generate an entangled (so-called continuous-variable) quantum state that contains contributions from high-order photon terms that lie outside the span of the dual-rail basis, which is detrimental to most applications. One often uses low pump power to mitigate the effects of those high-order terms. However, that reduces the pair generation rate, and the source becomes inherently probabilistic. We investigate a cascaded source that performs a linear-optical entanglement swap between two SPDC sources, to generate a heralded photonic entangled state that has a higher fidelity (to the ideal Bell state) compared to a free-running SPDC source. Furthermore, with the Bell swap providing a heralding trigger, we show how to build a multiplexed source, which despite reasonable switching losses and detector loss and noise, yields a fidelity versus success probability trade-off of a high-efficiency source of high-fidelity dual-rail photonic entanglement. We find, however, that there is a threshold of 1.5 dB of loss per switch, beyond which multiplexing hurts the fidelity versus success probability trade-off.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗

SDE_quark

This program provides a fixed-grid quadrature algorithm to compute integrals in the self-energy of the quark propagator within the Maris-Tandy model. The quark propagator both on t he spacelike real axis and at complex-valued momenta is determined from its Schwinger-Dyson equation (SDE). We first apply an iterative solver to find the quark propagator on the spacelike real axis. The propagator at complex-valued momenta is then computed from its self-energy based on this solution, where demanding integrals are encountered. In order to compute of these integrals, we apply customized variable transformations for the radial integral after subtracting the asymptotics. We subsequently apply an optional compound of quadrature rules for the angular integral.

Jia, Shaoyang↗

Identification of SPR Caverns with Multiple Oil Layers

This analysis shows that when lower density crude oil is injected into the top of an underground salt storage cavern containing more dense crude, separate oil phases can form and coexist indefinitely. This has been observed at the U.S. Strategic Petroleum Reserve in spite of geothermal heating and natural convection, which tend to mix the contents of containers with significant vertical extent subjected to wall and bottom heating. Such persistent layering can create operational challenges for meeting delivery specifications if high-value, low-vapor pressure oil becomes trapped below incoming low-density, high-vapor pressure oil, effectively blocking access to the lower layers until the top layer is removed. Previous conceptual models assumed that the oil injection process mixed incoming oil with resident oil in a storage cavern, forming a single oil phase with relatively homogeneous properties. Here, a review of historical data from the Strategic Petroleum Reserve revealed that several caverns contain multiple oil layers. As a result, oil layering needs to be another variable considered when planning oil movements at SPR in order to optimize low-vapor pressure oil availability to assist in oil delivery blending.

02 PETROLEUM↗