Search NASA⌕ Search

SEARCH · Search NASA

Results for “complex ensemble”

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 307 records · Page 17

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗

Screening of Therapeutic Agents for COVID-19 Using Machine Learning and Ensemble Docking Studies

The current pandemic demands a search for therapeutic agents against the novel coronavirus SARS-CoV-2. In this work, we present an efficient computational strategy that combines machine learning (ML)-based models and high-fidelity ensemble docking studies to enable rapid screening of possible therapeutic ligands. Targeting the binding affinity of molecules for either the isolated SARS-CoV-2 S-protein at its host receptor region or the S-protein:human ACE2 interface complex, we screen ligands from drug and biomolecule data sets that can potentially limit and/or disrupt the host-virus interactions. Top scoring one hundred eighty-seven ligands (with 75 approved by the Food and Drug Administration) are further validated by all atom docking studies. Important molecular descriptors ( 2 χ n , topological surface area, and ring count) and promising chemical fragments (oxolane, hydroxy, and imidazole) are identified to guide future experiments. Overall, this work expands our knowledge of small-molecule treatment against COVID-19 and provides a general screening pathway (combining quick ML models with expensive high-fidelity simulations) for targeting several chemical/biochemical problems.

60 APPLIED LIFE SCIENCES↗

Revealing the evolution of order in materials microstructures using multi-modal computer vision

The development of high-performance materials for microelectronics, energy storage, and extreme environments depends on our ability to describe and direct property-defining microstructural order. Our present understanding is typically derived from laborious manual analysis of imaging and spectroscopy data, which is difficult to scale, challenging to reproduce, and lacks the ability to reveal latent associations needed for mechanistic models. Here, we demonstrate a multi-modal machine learning (ML) approach to describe order from electron microscopy analysis of the complex oxide La 1−x Sr x FeO 3 . We construct a hybrid pipeline based on fully and semi-supervised classification, allowing us to evaluate both the characteristics of each data modality and the value each modality adds to the ensemble. We observe distinct differences in the performance of uni- and multi-modal models, from which we draw general lessons in describing crystal order using computer vision.

36 MATERIALS SCIENCE↗

Impact of the Earth’s Density Profile on Atmospheric Neutrino Oscillations

As atmospheric neutrinos traverse the Earth, the matter potential influences their oscillation probabilities in intricate ways, enhancing and suppressing the conversion of one neutrino’s flavor to another along the propagation. Comprehending their behavior within Earth’s complex density profile, primarily described by the Preliminary Reference Earth Model (PREM) [1], is essential for accurately describing oscillations and determining expected atmospheric neutrino event rates in the DUNE Far Detector. In this study, we address this challenge by considering an ensemble of Earth models constrained by astronomical measurements of the planet’s mass and moment of inertia. We evaluate how variations in densities and layer boundaries can affect oscillation probabilities and event rates, and what implications this could have on various physics analyses.

Ismerio Oliveira, Marcelo [Rio de Janeiro, Pont. U↗

Modeling human target acquisition in ground-to-air weapon systems

The problems associated with formulating and validating mathematical models for describing and predicting human target acquisition response are considered. In particular, the extension of the human observer model to include the acquisition phase as well as the tracking segment is presented. Relationship of the Observer model structure to the more complex Standard Optimal Control model formulation and to the simpler Transfer Function/Noise representation is discussed. Problems pertinent to structural identifiability and the form of the parameterization are elucidated. A systematic approach toward the identification of the observer acquisition model parameters from ensemble tracking error data is presented.

Phatak, A. V.↗

Microphysics in Goddard Multi-Scale Modeling Systems

Advances in computing power allow atmospheric prediction and general circulation models to be run at progressively finer scales of resolution, using increasingly more sophisticated physical parameterizations. The representation of cloud microphysical processes is one of key components of these models. In addition, over the past decade both research and operational numerical weather prediction models have started using more complex microphysical schemes that were originally developed for high-resolution cloud-resolving models (CRMs). In the paper, we described different microphysics schemes that are used in Goddard Multi-scale Modeling System. There are three major models, Goddard Cumulus Ensemble (GCE), NASA Unified Weather Research Forecast (NU-WRF) and Multi-scale Modeling Framework (MMF) model, in this modeling system. The microphysics schemes are Goddard three class ice (3ICE) and four class (4ICE) scheme, Morrison two moments (2M) 3ICE, Colorado State University Regional Atmospheric Modeling System (RAMS) 2M five class ice (5ICE) and spectral bin microphysics schemes. The performance of these schemes are examined and compared with radar and satellite observation. In addition, the inter-comparison with different microphysics schemes are conducted. Current and future observations needed for microphysics schemes evaluation as well as major characteristics of current microphysics are discussed.

Tao, W.-K.↗

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale of post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of the impacts of fire on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify the impacts of wildfire on runoff generation. We consider the disturbances in the vegetation and soil as the two main factors contributing to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We consider the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration for two case studies. These case studies demonstrate the general applicability of hydrophobicity formulation to serve as a guideline for exploring the range of hydrologic responses post-fire.

Fire Hydrology↗

Toward Resilient Heterogeneous Computing Workflow through Kokkos-DataSpaces Integration

With the growing number of applications designed for heterogeneous HPC devices, application programmers and users are finding it challenging to compose scalable workflows as ensembles of these applications, that are portable, performant and resilient. The Kokkos C++ library has been designed to simplify this cumbersome procedure by providing an intra-application uniform programming model and portable performance. However, assembling multiple Kokkos-enabled applications into a complex workflow is still a challenge. Although Kokkos enables a uniform programming model, the inter-application data exchange still remains a challenge from both performance and software development cost perspectives. In order to address this issue, we propose a Kokkos-DataSpaces Integration, with the goal of providing a virtual shared-space abstraction that can be accessed concurrently by all applications in an Kokkos workflow, thus extending Kokkos to support inter-application data exchange.

97 MATHEMATICS AND COMPUTING↗

Design and Performance of Kokkos Staging Space toward Scalable Resilient Application Couplings

With the growing number of applications designed for heterogeneous HPC devices, application programmers and users are finding it challenging to compose scalable workflows as ensembles of these applications, that are portable, performant and resilient. The Kokkos C++ library has been designed to simplify this cumbersome procedure by providing an intra-application uniform programming model and portable performance. However, assembling multiple Kokkos-enabled applications into a complex workflow is still a challenge. Although Kokkos enables a uniform programming model, the inter-application data exchange still remains a challenge from both performance and software development cost perspectives. In order to address this issue, we propose Kokkos data staging memory space, an extension of Kokkos' data abstraction (memory space) for heterogeneous computing systems. This new abstraction allows to express data on a virtual shared-space for multiple Kokkos applications, thus extending Kokkos to support inter-application data exchange to build an efficient application workflow. Additionally, we study the effectiveness of asynchronous data layout conversions for applications requiring different memory access patterns for the shared data. Our preliminary evaluation with a synthetic benchmark indicate the effectiveness of this conversion adapted to three different scenarios representing access frequency and use patterns of the shared data.

97 MATHEMATICS AND COMPUTING↗

Robust Strategy for Rocket Engine Health Monitoring

Monitoring the health of rocket engine systems is essentially a two-phase process. The acquisition phase involves sensing physical conditions at selected locations, converting physical inputs to electrical signals, conditioning the signals as appropriate to establish scale or filter interference, and recording results in a form that is easy to interpret. The inference phase involves analysis of results from the acquisition phase, comparison of analysis results to established health measures, and assessment of health indications. A variety of analytical tools may be employed in the inference phase of health monitoring. These tools can be separated into three broad categories: statistical, rule based, and model based. Statistical methods can provide excellent comparative measures of engine operating health. They require well-characterized data from an ensemble of "typical" engines, or "golden" data from a specific test assumed to define the operating norm in order to establish reliable comparative measures. Statistical methods are generally suitable for real-time health monitoring because they do not deal with the physical complexities of engine operation. The utility of statistical methods in rocket engine health monitoring is hindered by practical limits on the quantity and quality of available data. This is due to the difficulty and high cost of data acquisition, the limited number of available test engines, and the problem of simulating flight conditions in ground test facilities. In addition, statistical methods incur a penalty for disregarding flow complexity and are therefore limited in their ability to define performance shift causality. Rule based methods infer the health state of the engine system based on comparison of individual measurements or combinations of measurements with defined health norms or rules. This does not mean that rule based methods are necessarily simple. Although binary yes-no health assessment can sometimes be established by relatively simple rules, the causality assignment needed for refined health monitoring often requires an exceptionally complex rule base involving complicated logical maps. Structuring the rule system to be clear and unambiguous can be difficult, and the expert input required to maintain a large logic network and associated rule base can be prohibitive.

Santi, L. Michael↗

Calibration and Uncertainty Estimation Using the Ensemble Kalman Filter with a Large Subsurface Flow and Transport Model - 20321

At routinely monitored groundwater contamination sites, periodically measured environmental conditions such as groundwater levels and contaminant concentrations are used to inform and confirm a conceptual site model (CSM) and guide the development and calibration of a numerical groundwater flow and transport model. The calibration of groundwater flow and transport models after each measurement (sampling) event can illuminate deficiencies in a CSM, identify areas where additional monitoring is warranted, and predict the behavior of the system to guide decision making. However, manual and automated (e.g. PEST) model calibration tools can be time-consuming and computationally expensive to implement after each sampling event. Perhaps as a result, such calibration tools generally utilize all available monitoring data simultaneously rather than sequentially assimilating monitoring data one sampling event at a time as the results from sampling become available. A more real-time data assimilation approach may reduce parameter uncertainty, quantify the value of additional monitoring data, and produce a usable model more quickly and with less effort. To mitigate the potential time-consuming aspects of manual and widely applied automated calibration techniques, a data assimilation algorithm called the ensemble Kalman filter (EnKF) was evaluated as a relatively efficient method of model calibration and uncertainty assessment via the sequential integration of monitoring data into a model. The EnKF was able to successfully and efficiently assimilate monitoring and modeling data to calibrate a complex flow and transport model at a real-world site with significant subsurface heterogeneity, uncertainty, and 12 years of monitoring data (over 4,000 individual measurements of groundwater levels and over 2,500 measurements of contaminant concentrations). Starting with an uncalibrated model data from annual sampling events were sequentially assimilated, and the resultant predication errors and estimated parameter uncertainties were tracked. After all monitoring data were assimilated, both flow and transport residuals at the end of the EnKF process were comparable to those produced via a concurrent PEST calibration effort but required fewer model simulations. Both uncertainty and prediction errors decreased over time. In a real-time application, the adequacy of the model could be assessed after each sampling event. The benefits of such a real-time approach to utilizing monitoring data include reduced costs (in the form of model updates or site characterization efforts), early flagging of possible errors in the CSM, and a reduced risk of overfitting and corresponding increased confidence in model predictions. This tool may be particularly useful compared to other calibration techniques (e.g. manual, PEST) when model runtimes are long, calibration parameters are many, or parameter uncertainty is large. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Determining Infrared Optical Constants for Dolomite using Single Angle Reflectance and Spectroscopic Ellipsometry

The optical constants, namely the real (n) and imaginary (k) parts of the complex refractive index, are of particular interest to generate the infrared (IR) spectra of liquid and solid materials in different morphologies. To obtain n/k, however, most materials are typically not found in the monolithic forms necessary to easily measure n/k, and thus require the use of other methods such as pressing powders to form planar, specularly-reflective pellets. In this work dolomite crystals are measured using fixed-angle IR reflectance spectroscopy in both 1) a monolithic form with the crystal fixed in epoxy and polished, but data are also recorded from 2) pressed-pellet forms of powders of the same mineral using the same spectroscopic protocols. First results for the two methods are compared. It was found that the measured reflectance can vary by as much as a factor of two between the dolomite crystals and the pressed powder forms. For the two-sample preparation approaches the preliminary spectra are compared and the implications and limitations of each method for determining the optical constants of given materials are discussed. Comparison to literature data suggest that polarization effects likely account for the differing amplitude results for reflectance (and hence the k -vectors): Dolomite is known to be biaxial with significantly differing ordinary and extraordinary rays; the pressed pellet method measures an ensemble of microcrystals in randomly oriented positions, whereas the single crystal maintains just one orientation relative to the optical axis.

ellipsometry, dolomite, complex refractive index, ↗

Reactive Sequencing for Autonomous Navigation Evolving from Phoenix Entry, Descent, and Landing

Virtual Machine Language (VML) is an award-winning advanced procedural sequencing language in use on NASA deep-space missions since 1997, and was used for the successful entry, descent, and landing (EDL) of the Phoenix spacecraft onto the surface of Mars. Phoenix EDL utilized a state-oriented operations architecture which executed within the constraints of the existing VML 2.0 flight capability, compatible with the linear "land or die" nature of the mission. The intricacies of Phoenix EDL included the planned discarding of portions of the vehicle, the complex communications management for relay through on-orbit assets, the presence of temporally indeterminate physical events, and the need to rapidly catch up four days of sequencing should a reboot of the spacecraft flight computer occur shortly before atmospheric entry. These formidable operational challenges led to new techniques for packaging and coordinating reusable sequences called blocks using one-way synchronization via VML sequencing global variable events. The coordinated blocks acted as an ensemble to land the spacecraft, while individually managing various elements in as simple a fashion as possible. This paper outlines prototype VML 2.1 flight capabilities that have evolved from the one-way synchronization techniques in order to implement even more ambitious autonomous mission capabilities. Target missions for these new capabilities include autonomous touch-and-go sampling of cometary and asteroidal bodies, lunar landing of robotic missions, and ultimately landing of crewed lunar vehicles. Close proximity guidance, navigation, and control operations, on-orbit rendezvous, and descent and landing events featured in these missions require elaborate abort capability, manifesting highly non-linear scenarios that are so complex as to overtax traditional sequencing, or even the sort of one-way coordinated sequencing used during EDL. Foreseeing advanced command and control needs for small body and lunar landing guidance, navigation and control scenarios, work began three years ago on substantial upgrades to VML that are now being exercised in scenarios for lunar landing and comet/asteroid rendezvous. The advanced state-based approach includes coordinated state transition machines with distributed decision-making logic. These state machines are not merely sequences - they are reactive logic constructs capable of autonomous decision making within a well-defined domain. Combined with the JPL's AutoNav software used on Deep Space 1 and Deep Impact, the system allows spacecraft to autonomously navigate to an unmapped surface, soft-contact, and either land or ascend. The state machine architecture enabled by VML 2.1 has successfully performed sampling missions and lunar descent missions in a simulated environment, and is progressing toward flight capability. The authors are also investigating using the VML 2.1 flight director architecture to perform autonomous activities like rendezvous with a passive hypothetical Mars sample return capsule. The approach being pursued is similar to the touch-and-go sampling state machines, with the added complications associated with the search for, physical capture of, and securing of a separate spacecraft. Complications include optically finding and tracking the Orbiting Sample Capsule (OSC), keeping the OSC illuminated, making orbital adjustments, and physically capturing the OSC. Other applications could include autonomous science collection and fault compensation.

asteroid↗

Experimental realization of an extended Fermi-Hubbard model using a 2D lattice of dopant-based quantum dots

The Hubbard model is an essential tool for understanding many-body physics in condensed matter systems. Artificial lattices of dopants in silicon are a promising method for the analog quantum simulation of extended Fermi-Hubbard Hamiltonians in the strong interaction regime. However, complex atom-based device fabrication requirements have meant emulating a tunable two-dimensional Fermi-Hubbard Hamiltonian in silicon has not been achieved. Here, we fabricate 3 × 3 arrays of single/few-dopant quantum dots with finite disorder and demonstrate tuning of the electron ensemble using gates and probe the many-body states using quantum transport measurements. By controlling the lattice constants, we tune the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from metallic to Mott insulating behavior. We simulate thermally activated hopping and Hubbard band formation using increased temperatures. As atomically precise fabrication continues to improve, these results enable a new class of engineered artificial lattices to simulate interactive fermionic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Programmable interactions and emergent geometry in an array of atom clouds

Interactions govern the flow of information and the formation of correlations between constituents of many-body quantum systems, dictating phases of matter found in nature and forms of entanglement generated in the laboratory. Typical interactions decay with distance and thus produce a network of connectivity governed by geometry—such as the crystalline structure of a material or the trapping sites of atoms in a quantum simulator. However, many envisioned applications in quantum simulation and computation require more complex coupling graphs including non-local interactions, which feature in models of information scrambling in black hole, and mappings of hard optimization problems onto frustrated classical magnets. Here we describe the realization of programmable non-local interactions in an array of atomic ensembles within an optical cavity, in which photons carry information between atomic spins. By programming the distance dependence of the interactions, we access effective geometries for which the dimensionality, topology and metric are entirely distinct from the physical geometry of the array. As examples, we engineer an antiferromagnetic triangular ladder, a Möbius strip with sign-changing interactions and a treelike geometry inspired by concepts of quantum gravity. The tree graph constitutes a toy model of holographic duality, in which the quantum system lies on the boundary of a higher-dimensional geometry that emerges from measured correlations. Overall, our work provides broader prospects for simulating frustrated magnets and topological phases, investigating quantum optimization paradigms and engineering entangled resource states for sensing and computation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Direct numerical simulation of reacting flows

The objectives of this work are: (1) to extend the technique of direct numerical simulations to turbulent, chemically reacting flows, (2) to test the validity of the method by comparing computational results with laboratory data, and (3) to use the simulations to gain a better understanding of the effects of turbulence on chemical reactions. The effects of both the large scale structure and the smaller scale turbulence on the overall reaction rates are addressed. The relationship between infinite reaction rate and finite reaction rate chemistry is compared with some of the results of calculations with existing theories and laboratory data. The direct numerical simulation method involves the numerical solution of the detailed evolution of the complex turbulent velocity and concentration fields. Using very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear (possibly low pass filtered) equations of motion are solved and no closure assumptions or turbulence models are used. Statistical data are obtained by performing spatial, temporal, and/or ensemble averages over the computed flow fields.

Riley, J. J.↗

Visualization and Quantification of Rotor Tip Vortices in Helicopter Flows

Helicopter aeromechanics encompasses a highly vortical flow field. The vortices generated at each blade tip contain unsteady, complex, three-dimensional structures, which interact with each other, other blades, the fuselage and various components of the helicopter. It is crucial to understand vortex kinematics and their subsequent dynamic evolution. Much research has been devoted to the understanding of helicopter vortex dynamics, including a number of experimental studies.1-6 In May 2010 Particle Image Velocimetry (PIV) measurements of a full-scale UH-60A rotor were acquired in the National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel.1 These measurements were taken at a plane just downstream of the advancing blade in the vicinity of the blade tip—the so-called PIV plane. The resulting PIV data were then processed using an ensemble-average approach to create graphical representations of the vortical wake velocity and vorticity fields, which, in turn, have enhanced the understanding of rotorcraft vortical wake flow field physics and have provided a more detailed validation of vortical wake computer simulations.7 A common approach used to analyze flow field features is to compute and plot color contour maps of various scalar quantities such as pressure, velocity magnitude and vorticity magnitude. For example, the color map of the vorticity magnitude is typically used to determine vortical flow structure. With this approach the vortex core may appear larger or smaller, depending on the contour levels that are selected. Thus, the resulting visualization is sensitive to user-specified contour levels. For vortex core radius measurements, it is more accurate to calculate the vortex core radius using the cross-flow velocity profile across the vortex core. The task of extracting the cross-flow velocity profile can be time consuming with existing tools since the user needs to manually select the core center then specify sampling points along the profile axis. The task becomes even more challenging when the associated grid system uses AMR (Adaptive Mesh Refinement) where the profile axis could span multiple grid blocks. There are a number of existing techniques for profiling of vortex core attributes;8-9 however, these techniques are not fully automatic in that the user still needs to select the vortex core center to compute the cross-flow velocity profile. The present study introduces a new color map scheme that is based on the vortex core radius, which is fully automatic and does not require user intervention. Analysis and visualization of blade tip vortices on the PIV plane using the proposed new color map scheme are described in Section II. The new approach is evaluated using two case studies, which are described in Section III. The paper ends with a summary in Section IV.

Quantification↗

Synergy of Observations and Dynamo Models to Understand and Predict Solar Activity Cycles

The long-standing problem of understanding the evolution of the global magnetic fields that drive solar activity through different temporal scales is becoming more tractable because, in addition to 400 years of sunspot records, we now have almost 4 solar cycles of magnetic field observations. These observations allow us to discern physical connections between dynamo model variables and observations using data assimilation analysis. In particular, the Ensemble Kalman Filter approach takes into account uncertainties in both observations and modeling and allows us to make reliable forecasts of solar cycle activity by using a relatively simple non-linear dynamical model of the solar dynamo. To expand this approach for more complex 2D and 3D dynamo modeling, it is necessary to decompose the observed synoptic magnetograms into poloidal and toroidal field components. In this presentation I will present initial results on magnetogram decomposition and assimilation of magnetogram data into dynamo modeling.

Kitiashvili, Irina N.↗