Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained dynamics”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Exploring the Longwave Radiative Effects of Dust Aerosols

Dust aerosols not only affect air quality and visibility where they pose a significant health and safety risk, but they can also play a role in modulating the energy balance of the Earth-atmosphere system by directly interacting with local radiative fields. Consequently, dust aerosols can impact regional climate patterns such as changes in precipitation and the evolution of the hydrological cycle. Assessing the direct effect of dust aerosols at the solar wavelengths is fairly straightforward due in part to the relatively large signal-to-noise ratio in broadband irradiance measurements. The longwave (LW) impacts, on the other hand, are rather difficult to ascertain since the measured dust signal level (10 Wm-2) is on the same order as the instrumental uncertainties. Moreover, compared to the shortwave (SW), limited experimental data on the LW optical properties of dust makes it a difficult challenge for constraining the LW impacts. Owing to the strong absorption features found in many terrestrial minerals (e.g., silicates and clays), the LW effects, although much smaller in magnitude compared to the SW, can still have a sizeable impact on the energetics of the Earth-atmosphere system, which can potentially trigger changes in the heat and moisture surface budgets, and dynamics of the atmosphere. The current endeavor is an integral part of an on-going research study to perform detailed assessments of dust direct aerosol radiative effects (DARE) using comprehensive global datasets from NASA Goddards mobile ground-based facility (cf. http://smartlabs.gsfc.nasa.gov/) during previous field experiments near key dust source regions. Here we examine and compare the results from two of these studies: the 2006 NASA African Monsoon Multidisciplinary Activities and the 2008 Asian Monsoon Years. The former study focused on transported Saharan dust at Sal Island (16.73N, 22.93W), Cape Verde along the west coast of Africa while the latter focused on Asian dust at Zhangye China (39.082N, 100.276E), a semi-arid region between the Taklimakan and Gobi deserts. NASA Goddards Giovanni system is used to help map out the spatial distribution of retrieved aerosol optical depths across the latter desert regions. 1-D radiative transfer model constrained by local measurements, including spectral photometry/interferometry and lidar for characterizing the spatiotemporal variability in dust properties and atmospheric conditions, is employed to evaluate the local instantaneous LW DARE of dust both at the surface and at the top of the atmosphere along with heating rate profiles for cloud-free atmospheres. The efficiency in LW DARE and its significance relative to the diurnally averaged SW effects are explored and compared in both studies. Found to be non-negligible, LW DARE is an important component in the study of regional climate variation with important implications for more detailed global assessments.

Hansell, Richard A., Jr.↗
DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Observations and Modeling of the Onset of Fast Reconnection in the Solar Transition Region

Magnetic reconnection is a fundamental plasma process that plays a critical role not only in energy release in the solar atmosphere, but also in fusion, astrophysics, and other space plasma environments. One of the challenges in explaining solar observations in which reconnection is thought to play a critical role is to account for the transition of the dynamics from a slow quasi-continuous phase to a fast and impulsive energetic burst of much shorter duration. Despite the theoretical progress in identifying mechanisms that might lead to rapid onset, a lack of observations of this transition has left models poorly constrained. High-resolution spectroscopic observations from NASA's Interface Region Imaging Spectrograph now reveal tell-tale signatures of the abrupt transition of reconnection from a slow phase to a fast, impulsive phase during UV bursts or explosive events in the Sun's atmosphere. Our observations are consistent with numerical simulations of the plasmoid instability, and provide evidence for the onset of fast reconnection mediated by plasmoids and new opportunities for remote-sensing diagnostics of reconnection mechanisms on the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Active suppression of an 'apparent shock induced instability'

A control law was designed, using constrained optimization techniques, to suppress an apparent shock induced instability of a sweptback, aeroelastic wing with supercritical airfoil sections. The controller design was based on an approximate linear plant representation obtained using forced response data from a previous entry in the Langley Transonic Dynamics tunnel. During a second tunnel entry, it was found that there was not an instability in the uncontrolled case but there was a region of very low damping (high dynamic response) near a Mach number of 0.92. Controller performance was obtained during the test in near real-time and revealed that the controller attenuated the open-loop response and provided a small but significant amount of damping over a Mach number range from M = 0.70 to M = 0.92.

Adams, William M., Jr.↗
NASA NTRS2006

Distributed Constrained Optimization with Semicoordinate Transformations

Recent work has shown how information theory extends conventional full-rationality game theory to allow bounded rational agents. The associated mathematical framework can be used to solve constrained optimization problems. This is done by translating the problem into an iterated game, where each agent controls a different variable of the problem, so that the joint probability distribution across the agents moves gives an expected value of the objective function. The dynamics of the agents is designed to minimize a Lagrangian function of that joint distribution. Here we illustrate how the updating of the Lagrange parameters in the Lagrangian is a form of automated annealing, which focuses the joint distribution more and more tightly about the joint moves that optimize the objective function. We then investigate the use of "semicoordinate" variable transformations. These separate the joint state of the agents from the variables of the optimization problem, with the two connected by an onto mapping. We present experiments illustrating the ability of such transformations to facilitate optimization. We focus on the special kind of transformation in which the statistically independent states of the agents induces a mixture distribution over the optimization variables. Computer experiment illustrate this for &sat constraint satisfaction problems and for unconstrained minimization of NK functions.

Macready, William↗
DOE OSTI2023

Efficient Super-Resolution of Near-Surface Climate Modeling Using the Fourier Neural Operator

Downscaling methods are critical in efficiently generating high-resolution atmospheric data. However, state-of-the-art statistical or dynamical downscaling techniques either suffer from the high computational cost of running a physical model or require high-resolution data to develop a downscaling tool. Here, we demonstrate a recently proposed zero-shot super-resolution method, the Fourier neural operator (FNO), to efficiently perform downscaling without the need for high-resolution data. Because the FNO learns dynamics in Fourier space, FNO is a resolution-invariant emulator; it can be trained at a coarse resolution and produces emulation at any high resolution. We applied FNO to downscale a 4-km resolution Weather Research and Forecasting (WRF) Model simulation of near-surface heat-related variables over the Great Lakes region. The FNO is driven by the atmospheric forcings, and topographic features used in the WRF model at the same resolution. We incorporated a physics-constrained loss in FNO by using the Clausius-Clapeyron relation to better constrain the relations among the emulated states. Trained on merely 600 WRF snapshots at 4-km resolution, the FNO shows comparable performance with a widely used convolutional network, U-Net, achieving averaged modified Kling-Gupta Efficiency of 0.88 and 0.94 on the test dataset for temperature and pressure, respectively. We then employed the FNO to produce 1-km emulations to reproduce the fine climate features. Further, by taking the WRF simulation as ground truth, we show consistent performances at the two resolutions, suggesting the reliability of FNO in producing high-resolution dynamics. Our study demonstrates the potential of using FNO for zero-shot super-resolution in generating first-order estimation on atmospheric modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Warming induces divergent stomatal dynamics in co-occurring boreal trees

Abstract Climate warming will alter photosynthesis and respiration not only via direct temperature effects on leaf biochemistry but also by increasing atmospheric dryness, thereby reducing stomatal conductance and suppressing photosynthesis. Our knowledge on how climate warming affects these processes is mainly derived from seedlings grown under highly controlled conditions. However, little is known regarding temperature responses of trees growing under field settings. We exposed mature tamarack and black spruce trees growing in a peatland ecosystem to whole-ecosystem warming of up to +9°C above ambient air temperatures in an ongoing long-term experiment (SPRUCE: Spruce and Peatland Responses Under Changing Environments). Here, we report the responses of leaf gas exchange after the first two years of warming. We show that the two species exhibit divergent stomatal responses to warming and vapor pressure deficit. Warming of up to 9°C increased leaf N in both spruce and tamarack. However, higher leaf N in the warmer plots translate into higher photosynthesis in tamarack but not in spruce, with photosynthesis being more constrained by stomatal limitations in spruce than in tamarack under warm conditions. Surprisingly, dark respiration did not acclimate to warming in spruce, and thermal acclimation of respiration was only seen in tamarack once changes in leaf N were considered. Our results highlight how warming can lead to differing stomatal responses to warming in co-occurring species, with consequent effects on both vegetation carbon and water dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics

Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. In conclusion, we present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Obstruction and Interference in Low-Energy Models for Twisted Bilayer Graphene

The electronic bands of twisted bilayer graphene (TBLG) with a large-period moiré superlattice fracture to form narrow Bloch minibands that are spectrally isolated by forbidden energy gaps from remote dispersive bands. When these gaps are sufficiently large, one can study a band-projected Hamiltonian that correctly represents the dynamics within the minibands. This inevitably introduces nontrivial geometrical constraints that arise from the assumed form of the projection. Here we show that this choice has a profound consequence in a low-energy experimentally observable signature that therefore can be used to tightly constrain the analytic form of the appropriate low-energy theory. Additionally, we find that this can be accomplished by a careful analysis of the electron density produced by backscattering of Bloch waves from an impurity potential localized on the moiré superlattice scale. We provide numerical estimates of the effect that can guide experimental work to clearly discriminate between competing models for the low-energy band structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1993

Scheduling with partial orders and a causal model

In an ongoing project at Honeywell SRC, we are constructing a prototype scheduling system for a NASA domain using the 'Time Map Manager' (TMM). The TMM representations are flexible enough to permit the representation of precedence constraints, metric constraints between activities, and constraints relative to a variety of references (e.g., Mission Elapsed Time vs. Mission Day). The TMM also supports a simple form of causal reasoning (projection), dynamic database updates, and monitoring specified database properties as changes occur over time. The greatest apparent advantage to using the TMM is the flexibility added to the scheduling process: schedules are constructed by a process of 'iterative refinement,' in which scheduling decisions correspond to constraining an activity either with respect to another activity or with respect to one time line. The schedule becomes more detailed as activities and constraints are added. Undoing a scheduling decision means removing a constraint, not removing an activity from a specified place on the time line. For example, we can move an activity around on the time line by deleting constraints and adding new ones.

Boddy, Mark↗
NASA NTRS2003

Jupiter Atmospheric Science with JIMO: Linking Science Objectives and Measurement Goals

Although the primary focus of the Jupiter Icy Moons Orbiter (JIMO) mission will be the characterization and study of Jupiter's icy moons, there will be opportunities throughout the mission for unprecendented observations of Jupiter. With an adaptable suite of payload instruments, the atmospheric data collected by JIMO can help to answer fundamental questions about the largest planet in our solar system that remain after (or were generated by) previous spacecraft reconnaissance (e.g. Voyager, Galileo, and Cassini). Near-IR (0.7-4 micron) spectral imaging will most likely be used to identify mineralogies and ices on the Jovian satellites by virtue of their spectral signatures. This same capability is very well tailored for studies of Jovian atmospheric dynamics and structure. Near-IR methane absorption bands allow 2-D mapping of the horizontal wind field at size scales to tens of kms, as well as the height dependence of this field above the ammonia cloud deck (700 to a few mbar), constraining current models of atmospheric vertical structure. Likewise, atmospheric ice aerosols with unique spectroscopic signatures (ammonia ice near 1.5, 2.0, and 2.8 microns and water ice between 3.0 - 3.5 microns) can be detected and mapped using spectral difference imaging or spectrally inclusive principal-component methods. Spectral imaging of the Jovian aurora via (3)H(+) emission lines between 3 - 4 microns can be used to spatially map the interplay between the satellites) Jupiter's magnetosphere, and Jupiter's atmosphere. Each of these measurements addresses one or more fundamental questions related to the energy balance in Jupiter's atmosphere. All of these tunable imaging objectives can be achieved using acousto-optic tunable filters (AOTF's), which have been used for years in ground-based observing instruments and which have been proposed for numerous planetary missions. The application of this technology to the science objectives of both the icy satellites and Jovian atmospheric components of the JIMO mission will be discussed.

Chanover, N. J.↗
NASA NTRS2020

Capturing an Evolving Nebular Environment: A Petrographic and Geochemical Study of a Type A, B & C CAI

Calcium, Aluminum-rich Inclusions (CAIs) were the first formed solids in our Solar System, with mineral assemblages reflecting the first phases predicted to condense out of a hot nebular gas of Solar composition. Geochemical, textural and crystallographic information in CAIs can be used to constrain the temperature, pressure, and composition (e.g., oxygen fugacity) of the gaseous reservoir(s) from which they formed, as well as any secondary (nebular and parent body) processes they underwent. Coordinated geochemical and textural analyses provide information on nebular conditions (i.e., astrophysical environments and dynamics of nebular gas reservoirs) in which these CAIs formed. In order to better understand the evolution of nebular reservoirs at the time of CAI formation, we analyzed a Type A, B and C CAI using Electron Probe Micro-Analyzer (EPMA) and Electron BackScatter Diffraction (EBSD) at NASA Johnson Space Center (JSC).

Comins, M. B.↗
NASA NTRSDate not supplied

Temporal Changes in Astronauts’ Muscle and Cardiorespiratory Physiology Before, During and After Spaceflight

Background: NASA’s planned space exploration missions will require astronauts to safely perform extravehicular activity (EVA) and to safely egress vehicles in a variety of landing scenarios. Prolonged exposure to spaceflight can diminish tolerance for physical activity, decrease cardiovascular and sensorimotor function, and cause loss of bone mineral density, and reduced muscle mass and strength. Although exercise can mitigate these spaceflight-induced physiological decrements, little is known regarding the time-course of changes in muscle and aerobic performance during spaceflight. Furthermore, these exercise countermeasures are not fully protective. For example, maximal aerobic capacity (VO2pk), lower body muscle cross-sectional area, and strength decrease by about 10% to 15% after short- (~ 14 days) and long-duration (~ 6 months) missions on the International Space Station (ISS). Future space missions that are longer in duration and further from Earth will employ exploration vehicles that will have exercise hardware with less robust and more constrained exercise capabilities than of those available on the ISS. Thus, countermeasures will need to be optimized to protect crew health and performance on exploration-class missions that will last up to 3 years. This will require a more detailed understanding of the dynamic effects of spaceflight on human health and performance, and the ability of exercise to protect against this deconditioning, and the interaction of exercise with interrelated factors like nutrition, sleep, and environmental conditions. Methods: We will use standardized research and medical testing protocols previously validated in 1-G and 0-G to quantify the time course and the inter-individual variability of changes in physical performance, including cardiorespiratory fitness, and muscle strength, and endurance, before, during and after spaceflight missions lasting 2 months, 6 months, and 1 year. Additionally, we will use an extrapolation model to predict changes associated with 2–3-year exploration missions. Additionally, we will monitor in-flight exercise, nutrition, and sleep using in-flight assessment tools. Significance: Our testing protocols will provide valuable information for determining time course of change and the interindividual variability of spaceflight-induced deconditioning of aerobic capacity and muscle strength and endurance over the course of spaceflight missions up to and beyond 1 year. This information will be vital to assess whether humans can be physically ready for deep space exploration, such as on a mission to Mars, using current technology, or if additional mitigation strategies are necessary.

countermeasures↗
DOE OSTI2020

PyCDFT: A Python package for constrained density functional theory

In this paper, we present PyCDFT, a Python package to compute diabatic states using constrained density functional theory (CDFT). PyCDFT provides an object-oriented, customizable implementation of CDFT, and allows for both single-point self-consistent-field calculations and geometry optimizations. PyCDFT is designed to interface with existing density functional theory (DFT) codes to perform CDFT calculations where constraint potentials are added to the Kohn–Sham Hamiltonian. Here, we demonstrate the use of PyCDFT by performing calculations with a massively parallel first-principles molecular dynamics code, Qbox, and we benchmark its accuracy by computing the electronic coupling between diabatic states for a set of organic molecules. We show that PyCDFT yields results in agreement with existing implementations and is a robust and flexible package for performing CDFT calculations. The program is available at https://dx.doi.org/10.5281/zenodo.3821097.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Nonintrusive projection-based reduced order modeling using stable learned differential operators

Nonintrusive projection-based reduced order models (ROMs) are essential for dynamics prediction in multi-query applications where underlying governing equations are known but the access to the source of the underlying full order model (FOM) is unavailable; that is, FOM is a glass-box. This article proposes a learn-then-project approach for nonintrusive model reduction. In the first step of this approach, high-dimensional stable sparse learned differential operators (S-LDOs) are determined using the generated data. In the second step, the ordinary differential equations, comprising these S-LDOs, are used with suitable dimensionality reduction and low-dimensional subspace projection methods to provide equations for the evolution of reduced states. This approach allows easy integration into the existing intrusive ROM framework to enable nonintrusive model reduction while allowing the use of Petrov–Galerkin projections. The applicability of the proposed approach is demonstrated for Galerkin and LSPG projection-based ROMs through four numerical experiments: 1-D scalar advection, 1-D Burgers, 2-D scalar advection and 1-D scalar advection–diffusion–reaction equations. In conclusion, the results indicate that the proposed nonintrusive ROM strategy provides accurate and stable dynamics prediction.

42 ENGINEERING↗
DOE OSTI2024

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Leaf Trait Plasticity Alters Competitive Ability and Functioning of Simulated Tropical Trees in Response to Elevated Carbon Dioxide

The response of tropical ecosystems to elevated carbon dioxide (CO 2 ) remains a critical uncertainty in projections of future climate. Here, we investigate how leaf trait plasticity in response to elevated CO 2 alters projections of tropical forest competitive dynamics and functioning. We use vegetation demographic model simulations to quantify how plasticity in leaf mass per area and leaf carbon to nitrogen ratio alter the responses of carbon uptake, evapotranspiration, and competitive ability to a doubling of CO 2 in a tropical forest. Observationally constrained leaf trait plasticity in response to CO 2 fertilization reduces the degree to which tropical tree carbon uptake is affected by a doubling of CO 2 (up to -14.7% as compared to a case with no plasticity; 95% confidence interval [CI95%] -14.4 to -15.0). It also diminishes evapotranspiration (up to -7.0%, CI95% -6.4 to -7.7), and lowers competitive ability in comparison to a tree with no plasticity. Consideration of leaf trait plasticity to elevated CO 2 lowers tropical ecosystem carbon uptake and evapotranspirative cooling in the absence of changes in plant-type abundance. However, “plastic” responses to high CO 2 which maintain higher levels of plant productivity, many of which fall outside of the observed range of response, are potentially more competitively advantageous, thus, including changes in plant type abundance may mitigate these decreases in ecosystem functioning. Models that explicitly represent competition between plants with alternative leaf trait plasticity in response to elevated CO 2 are needed to capture these influences on tropical forest functioning and large-scale climate.

54 ENVIRONMENTAL SCIENCES↗