Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble 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 109 records · Page 6

Disturbance Dynamics in Transitional and Turbulent Flows

The dynamics of an ensemble of linear disturbances with a known probability distribution associated with the initial mode amplitudes are studied in boundary-layer flows through an analysis of the transport equations for the mean disturbance kinetic energy and disturbance energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis. Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field in proximity to the wall, the laminar regime fluctuations studied here are influenced across the boundary layer by the solid boundary. In addition to the low Reynolds number, early stage transition regime, the dynamics of these disturbance fields can be related in some respects to the near-wall dynamics of the fully turbulent regime.

Grosch, C. E.↗

A New Tropical Cyclone Dynamic Initialization Technique Using High Temporal and Spatial Density Atmospheric Motion Vectors and Airborne Field Campaign Data

Background: Initialization of tropical cyclones in numerical weather prediction (NWP) systems is a great challenge: Mass-wind field balance; Secondary circulation and heating; Asymmetries. There can be large adjustments in structure and intensity in the first 24 hours if the initial vortex is not in balance: Spurious gravity waves; Spin-up (model and physics). Existing mesoscale NWP model TC (Tropical Cyclone) initialization strategies: Bogus vortex, cold start from global analyses; 3DVAR or 4DVAR, possibly with synthetic observations; EnKF (Ensemble Kalman Filter); Dynamic initialization. Dynamic initialization allows vortex to have improved balance and physics spin-up at the initial time (e.g., Hendricks et al. 2013, 2011; Nguyen and Chen 2011; Fiorino and Warner 1981; Hoke and Anthes 1976). Himawari-8 geostationary satellite has capability of continuous imagery (10-minutes) over the full disk: New GOES-R satellites will have same capability. This will allow for unprecedented observations of tropical cyclones. However, current data assimila1on systems are not capable of ingesting such high temporal observations (Atmospheric Mo1on Vectors - AMVs). Hourly AMVs are produced, and thinned to 100-kilometer spacing in the horizontal. An entirely new data assimilation concept is required to utilize these observations.

remote sensing↗

Interactive Vegetation Phenology, Soil Moisture, and Monthly Temperature Forecasts

The time scales that characterize the variations of vegetation phenology are generally much longer than those that characterize atmospheric processes. The explicit modeling of phenological processes in an atmospheric forecast system thus has the potential to provide skill to subseasonal or seasonal forecasts. We examine this possibility here using a forecast system fitted with a dynamic vegetation phenology model. We perform three experiments, each consisting of 128 independent warm-season monthly forecasts: 1) an experiment in which both soil moisture states and carbon states (e.g., those determining leaf area index) are initialized realistically, 2) an experiment in which the carbon states are prescribed to climatology throughout the forecasts, and 3) an experiment in which both the carbon and soil moisture states are prescribed to climatology throughout the forecasts. Evaluating the monthly forecasts of air temperature in each ensemble against observations, as well as quantifying the inherent predictability of temperature within each ensemble, shows that dynamic phenology can indeed contribute positively to subseasonal forecasts, though only to a small extent, with an impact dwarfed by that of soil moisture.

Seasonal Forecasting↗

Electrochemically scrambled nanocrystals are catalytically active for CO 2 -to-multicarbons

Significance The electrocatalytic conversion of CO 2 to value-added products, especially valuable multicarbon products, is a pathway toward sustainable formation of chemicals and fuels typically derived from fossil fuels, while mitigating CO 2 emissions. Fundamental understanding and development of more efficient catalysts for this reaction require deep investigation into structures with high intrinsic activity, which are limited at present. This work comprehensively investigates a dynamic copper nanoparticle ensemble catalyst that significantly improves intrinsic activity of copper for multicarbon formation. Through concerted ex situ and in situ characterization techniques, it illustrates an electrochemically induced fusion of copper nanoparticles that result in a catalytically active disordered structure, motivating closer study of disordered metal nanostructures for C–C coupling electrocatalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling the Mechanisms and Implications of a Stronger Mid-Pliocene Atlantic Meridional Overturning Circulation (AMOC) in PlioMIP2

The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period in which the atmospheric CO2 concentration was approximately equal to the concentration we measure today (ca. 400 ppm). Sea surface temperature (SST) proxies indicate above-average warming over the North Atlantic in the mid-Pliocene with respect to the pre-industrial period, which may be linked to an intensified Atlantic Meridional Overturning Circulation (AMOC). Earlier results from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) show that the ensemble simulates a stronger AMOC in the mid-Pliocene than in the pre-industrial. However, no consistent relationship between the stronger mid-Pliocene AMOC and either the Atlantic northward ocean heat transport (OHT) or average North Atlantic SSTs has been found. In this study, we look further into the drivers and consequences of a stronger AMOC in mid-Pliocene compared to pre-industrial simulations in PlioMIP2. We find that all model simulations with a closed Bering Strait and Canadian Archipelago show reduced freshwater transport from the Arctic Ocean into the North Atlantic. This contributes to an increase in salinity in the subpolar North Atlantic and Labrador Sea that can be linked to the stronger AMOC in the mid-Pliocene. To investigate the dynamics behind the ensemble's variable response of the total Atlantic OHT to the stronger AMOC, we separate the Atlantic OHT into two components associated with either the overturning circulation or the wind-driven gyre circulation. While the ensemble mean of the overturning component is increased significantly in magnitude in the mid-Pliocene, it is partly compensated by a reduction in the gyre component in the northern subtropical gyre region. This indicates that the lack of relationship between the total OHT and AMOC is due to changes in OHT by the subtropical gyre. The overturning and gyre components should therefore be considered separately to gain a more complete understanding of the OHT response to a stronger mid-Pliocene AMOC. In addition, we show that the AMOC exerts a stronger influence on North Atlantic SSTs in the mid-Pliocene than in the pre-industrial, providing a possible explanation for the improved agreement of the PlioMIP2 ensemble mean SSTs with reconstructions in the North Atlantic.

AMOC↗

CryoFold: Determining protein structures and data-guided ensembles from cryo-EM density maps

Cryoelectron microscopy requires molecular modeling for refinement of structures. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins. Here, we introduce CryoFold, a pipeline of molecular dynamics simulations that determines ensembles of protein structures by integrating density data of varying sparsity at 3–5 Å resolution with sequence information and coarse-grained topological knowledge of the protein folds. We present six examples, folding proteins between 72 and 2,000 residues, including large membrane and multi-domain systems, and results from two Electron Microscopy Data Bank (EMDB) competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python graphical user interface (GUI).

59 BASIC BIOLOGICAL SCIENCES↗

Improved Predictability of the Indian Ocean Dipole Using a Stochastic Dynamical Model Compared to the North American Multimodel Ensemble Forecast

This study assesses the predictive skill of eight North American Multimodel Ensemble (NMME) models in predicting the Indian Ocean dipole (IOD). We find that the forecasted ensemble-mean IOD–El Niño–Southern Oscillation (ENSO) relationship deteriorates away from the observed relationship with increasing lead time, which might be one reason that limits the IOD predictive skill in coupled models. We are able to improve the IOD predictive skill using a recently developed stochastic dynamical model (SDM) forced by forecasted ENSO conditions. The results are consistent with the previous result that operational IOD predictability beyond persistence at lead times beyond one season is mostly controlled by ENSO predictability and the signal-to-noise ratio of the Indo-Pacific climate system. The multimodel ensemble (MME) investigated here is found to be of superior skill compared to each individual model at most lead times. Importantly, the skill of the SDM IOD predictions forced with forecasted ENSO conditions were either similar or better than those of the MME IOD forecasts. Moreover, the SDM forced with observed ENSO conditions exhibits significantly higher IOD prediction skill than the MME at longer lead times, suggesting the large potential skill increase that could be achieved by improving operational ENSO forecasts. We find that both cold and warm biases of the predicted Niño-3.4 index may cause false alarms of negative and positive IOD events, respectively, in NMME models. Many false alarms for IOD forecasts at lead times longer than one season in the original forecasts disappear or are significantly reduced in the SDM forced by forecasted ENSO conditions.

54 ENVIRONMENTAL SCIENCES↗

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scenario Storyline Discovery for Planning in Multi‐Actor Human‐Natural Systems Confronting Change

Scenarios have emerged as valuable tools in managing complex human-natural systems, but the traditional approach of limiting focus on a small number of predetermined scenarios can inadvertently miss consequential dynamics, extremes, and diverse stakeholder impacts. Exploratory modeling approaches have been developed to address these issues by exploring a wide range of possible futures and identifying those that yield consequential vulnerabilities. However, vulnerabilities are typically identified based on aggregate robustness measures that do not take full advantage of the richness of the underlying dynamics in the large ensembles of model simulations and can make it hard to identify key dynamics and/or storylines that can guide planning or further analyses. This study introduces the FRamework for Narrative Storylines and Impact Classification (FRNSIC; pronounced “forensic”): a scenario discovery framework that addresses these challenges by organizing and investigating consequential scenarios using hierarchical classification of diverse outcomes across actors, sectors, and scales, while also aiding in the selection of scenario storylines, based on system dynamics that drive consequential outcomes. We present an application of this framework to the Upper Colorado River Basin, focusing on decadal droughts and their water scarcity implications for the basin's diverse users and its obligations to downstream states through Lake Powell. We show how FRNSIC can explore alternative sets of impact metrics and drought dynamics and use them to identify drought scenario storylines, that can be used to inform future adaptation planning.

54 ENVIRONMENTAL SCIENCES↗

Mapping the global distribution of C 4 vegetation using observations and optimality theory

Plants with the C 4 photosynthesis pathway typically respond to climate change differently from more common C 3 -type plants, due to their distinct anatomical and biochemical characteristics. These different responses are expected to drive changes in global C 4 and C 3 vegetation distributions. However, current C 4 vegetation distribution models may not predict this response as they do not capture multiple interacting factors and often lack observational constraints. Here, we used global observations of plant photosynthetic pathways, satellite remote sensing, and photosynthetic optimality theory to produce an observation-constrained global map of C 4 vegetation. We find that global C 4 vegetation coverage decreased from 17.7% to 17.1% of the land surface during 2001 to 2019. This was the net result of a reduction in C 4 natural grass cover due to elevated CO 2 favoring C 3 -type photosynthesis, and an increase in C 4 crop cover, mainly from corn (maize) expansion. Using an emergent constraint approach, we estimated that C 4 vegetation contributed 19.5% of global photosynthetic carbon assimilation, a value within the range of previous estimates (18–23%) but higher than the ensemble mean of dynamic global vegetation models (14 ± 13%; mean ± one standard deviation). Our study sheds insight on the critical and underappreciated role of C 4 plants in the contemporary global carbon cycle.

54 ENVIRONMENTAL SCIENCES↗

Variable domain mutational analysis to probe the molecular mechanisms of high viscosity of an IgG 1 antibody

Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG 1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~200 anti-GCGR IgG 1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2–10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.

59 BASIC BIOLOGICAL SCIENCES↗

Quantum turbulence, superfluidity, non-Markovian dynamics, and wave function thermalization

While quantum turbulence has been addressed both experimentally (predominantly for superfluid He 4 and He 3 ) and theoretically, the dynamics of various ensembles of quantized vortices has been followed in time only until the vortices have decayed into phonons. How this “thermalization” is achieved is still an unaddressed and thus an unelucidated question. The unitary Fermi gas (UFG) is a unique quantum system, which has no classical counterpart and is of relevance to neutron stars, cold atoms, condensed-matter and nuclear many-body systems. The non-Markovian evolution of an isolated UFG is put in evidence and its entire nonequilibrium evolution can be studied theoretically within a unified theoretical framework. The initial lattice of quantum vortices and antivortices evolves through a couple of vortex tangles and excitation of Kelvin waves, where vortices cross and reconnect, until very slowly thermalization sets in. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantification of Sub-Pixel Dynamics in High-Speed Neutron Imaging

The high penetration depth of neutrons through many metals and other common materials makes neutron imaging an attractive method for non-destructively probing the internal structure and dynamics of objects or systems that may not be accessible by conventional means, such as X-ray or optical imaging. While neutron imaging has been demonstrated to achieve a spatial resolution below 10 μm and temporal resolution below 10 μs, the relatively low flux of neutron sources and the limitations of existing neutron detectors have, until now, dictated that these cannot be achieved simultaneously, which substantially restricts the applicability of neutron imaging to many fields of research that could otherwise benefit from its unique capabilities. In this work, we present an attenuation modeling approach to the quantification of sub-pixel dynamics in cyclic ensemble neutron image sequences of an automotive gasoline direct injector at a 5 μs time scale with a spatial noise floor in the order of 5 μm.

47 OTHER INSTRUMENTATION↗

Analysis of instrumentation error effects on the identification accuracy of aircraft parameters

An analytical investigation is presented of the effect of unmodeled measurement system errors on the accuracy of aircraft stability and control derivatives identified from flight test data. Such error sources include biases, scale factor errors, instrument position errors, misalignments, and instrument dynamics. Two techniques (ensemble analysis and simulated data analysis) are formulated to determine the quantitative variations to the identified parameters resulting from the unmodeled instrumentation errors. The parameter accuracy that would result from flight tests of the F-4C aircraft with typical quality instrumentation is determined using these techniques. It is shown that unmodeled instrument errors can greatly increase the uncertainty in the value of the identified parameters. General recommendations are made of procedures to be followed to insure that the measurement system associated with identifying stability and control derivatives from flight test provides sufficient accuracy.

Sorensen, J. A.↗

Mathematical and Computational Challenges in Population Biology and Ecosystems Science

Mathematical and computational approaches provide powerful tools in the study of problems in population biology and ecosystems science. The subject has a rich history intertwined with the development of statistics and dynamical systems theory, but recent analytical advances, coupled with the enhanced potential of high-speed computation, have opened up new vistas and presented new challenges. Key challenges involve ways to deal with the collective dynamics of heterogeneous ensembles of individuals, and to scale from small spatial regions to large ones. The central issues-understanding how detail at one scale makes its signature felt at other scales, and how to relate phenomena across scales-cut across scientific disciplines and go to the heart of algorithmic development of approaches to high-speed computation. Examples are given from ecology, genetics, epidemiology, and immunology.

Levin, Simon A.↗

Electronic Transport through Self Assembled Thiol Molecules: Effect of Monolayer Order, Dynamics and Temperature

We present the charge transport and tunneling conductance of self assembled organic thiol molecules and discuss the influence of order and dynamics in the monolayer on the transport behavior and the effect of temperature. Conjugated thiol molecular wires and organometals such as terpyridine metal complexes provide a new platform for molecular electronic devices and we study their self assembly on Au(111) substrates by the scanning tunneling microscope. Determining the organization of the molecule and the ability to control the nature of its interface with the substrate is important for reliable performance of the molecular electronic devices. By concurrent scanning tunneling microscopy and spectroscopy studies on SAMs formed from oligo (phenelyne ethynelyne) monolayers with and without molecular order, we show that packing and order determine the response of a self assembled monolayer (SAM) to competing interactions. Molecular resolution STM imaging in vacuum shows that the OPES adopt an imcommensurate SAM structure on Au(111) with a rectangular unit cell. Tunneling spectroscopic measurements were performed on the SAM as a function of junction resistance. STS results show that the I-Vs are non linear and asymmetric due to the inherent asymmetry in the molecular structure, with larger currents at negative sample biases. The asymmetry increases with increasing junction resistance due to the asymmetry in the coupling to the leads. This is brought out clearly in the differential conductance, which also shows a gap at the Fermi level. We also studied the effect of order and dynamics in the monolayer on the charge transport and found that competing forces between the electric field, intermolecular interactions, tip-molecule physisorption and substrate-molecule chemisorption impact the transport measurements and its reliability and that the presence of molecular order is very important for reproducible transport measurements. Thus while developing new electronic platforms based on molecules, it is important to have a good control of the molecule-substrate interface, for the devices to perform reliably. While such a control would minimize fluctuations and dynamics in the ensemble, the real challenge is to develop device architectures that are tolerant to fluctuations, since they cannot be totally eliminated in these low dimensional soft systems. Results of temperature dependent STS measurements will also be discussed.

Dholakia, Geetha↗

Information-Guided Dynamic Nuclear Polarization

A spin ensemble in thermal equilibrium continuously undergoes random fluctuations analogous to those observed in Brownian motion, a process known as spin noise. Here, we investigate the dynamics of a system comprising a spin-1 paramagnetic center and a hyperfine-coupled spin-1/2 nucleus in the vicinity of a level crossing. We theoretically show that nuclear spins polarize efficiently under the combined action of thermal fluctuations and a closed-loop feedback protocol. The latter articulates periodic observations of the electronic magnetization and radio-frequency pulses connecting hybrid states with opposite nuclear spin alignment. Since nuclear polarization emerges from electronic spin fluctuations, not spin order, this microwave-free technique generically benefits from warmer, not colder, operation temperatures. Furthermore, because the spin dynamics at play near a level crossing is rather insensitive to the absolute value of the magnetic field, our work promises opportunities for high-field dynamic nuclear polarization, which is difficult to attain through present methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗