Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random fields”

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 991 records · Page 55

Global field reconstruction from sparse sensors with Veronoi tessellation-assisted deep learning

Achieving accurate and robust global situational awareness of a complex time-evolving field from a limited number of sensors has been a longstanding challenge. This reconstruction problem is especially difficult when sensors are sparsely positioned in a seemingly random or unorganized manner, which is often encountered in a range of scientific and engineering problems. Moreover, these sensors can be in motion and can become online or offline over time. The key leverage in addressing this scientific issue is the wealth of data accumulated from the sensors. As a solution to this problem, we propose a data-driven spatial field recovery technique founded on a structured grid-based deep-learning approach for arbitrary positioned sensors of any numbers. It should be noted that the naïve use of machine learning becomes prohibitively expensive for global field reconstruction and is furthermore not adaptable to an arbitrary number of sensors. In the present work, we consider the use of Voronoi tessellation to obtain a structured-grid representation from sensor locations enabling the computationally tractable use of convolutional neural networks. One of the central features of the present method is its compatibility with deep-learning based super-resolution reconstruction techniques for structured sensor data that are established for image processing. The proposed reconstruction technique is demonstrated for unsteady wake flow, geophysical data, and three-dimensional turbulence. The current framework is able to handle an arbitrary number of moving sensors, and thereby overcomes a major limitation with existing reconstruction methods. The presented technique opens a new pathway towards the practical use of neural networks for real-time global field estimation.

Fukami, Kai↗

Probability Distribution Estimated From the Minimum, Maximum, and Most Likely Values: Applied to Turbine Inlet Temperature Uncertainty

Modern engineering design practices are tending more toward the treatment of design parameters as random variables as opposed to fixed, or deterministic, values. The probabilistic design approach attempts to account for the uncertainty in design parameters by representing them as a distribution of values rather than as a single value. The motivations for this effort include preventing excessive overdesign as well as assessing and assuring reliability, both of which are important for aerospace applications. However, the determination of the probability distribution is a fundamental problem in reliability analysis. A random variable is often defined by the parameters of the theoretical distribution function that gives the best fit to experimental data. In many cases the distribution must be assumed from very limited information or data. Often the types of information that are available or reasonably estimated are the minimum, maximum, and most likely values of the design parameter. For these situations the beta distribution model is very convenient because the parameters that define the distribution can be easily determined from these three pieces of information. Widely used in the field of operations research, the beta model is very flexible and is also useful for estimating the mean and standard deviation of a random variable given only the aforementioned three values. However, an assumption is required to determine the four parameters of the beta distribution from only these three pieces of information (some of the more common distributions, like the normal, lognormal, gamma, and Weibull distributions, have two or three parameters). The conventional method assumes that the standard deviation is a certain fraction of the range. The beta parameters are then determined by solving a set of equations simultaneously. A new method developed in-house at the NASA Glenn Research Center assumes a value for one of the beta shape parameters based on an analogy with the normal distribution (ref.1). This new approach allows for a very simple and direct algebraic solution without restricting the standard deviation. The beta parameters obtained by the new method are comparable to the conventional method (and identical when the distribution is symmetrical). However, the proposed method generally produces a less peaked distribution with a slightly larger standard deviation (up to 7 percent) than the conventional method in cases where the distribution is asymmetric or skewed. The beta distribution model has now been implemented into the Fast Probability Integration (FPI) module used in the NESSUS computer code for probabilistic analyses of structures (ref. 2).

Holland, Frederic A., Jr.↗

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗

On the wave excitation in the turbulent meteor trace

The solution to the problem of excitation of longitudinal and transverse electromagnetic waves in randomly inhomogeneous media is reduced to the derivation of a complex effective dielectric constant (EDC) tensor which nonlocally connects together the average macroscopic electromagnetic fields in continuous media, it is implied that the field values which are rapidly fluctuating on a microscopic scale in space and time become smoothed out in a specified way due to the inhomogeneous mixing of diffusion. Proceedings from the derived general expression for the EDC tensor new modes of longitudinal and transverse electromagnetic wave generation due to fluctuation in the parameters of the medium was predicted. In this connection, it is of interest to investigate the peculiarities of electromagnetic longitudinal and transverse wave propagation in such randomly inhomogeneous media where, apart from the charged particle concentration change, the random spatial and temporal changes of natural frequency of closely located oscillators take place.

Jandieri, G. V.↗

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗

Magnetic damping of rotation

Based on Wilson's (1977) article on the magnetic effects on space vehicles and other celestial bodies, the magnetic damping of rotation is considered. The inadequacy of the interstellar magnetic field in overcoming solar wind shielding and thus influencing the rotation of bodies is described. The ionospheric shielding of the interstellar field is discussed along with the permeability and magnetic damping by the solar or stellar wind. Star formation and angular momentum is discussed and attention is given to the magnetic damping of unshielded small bodies. Calculations of the rate for damping through random particle impact are made. Theories concerning the rotation of asteroids and the origin of meteorites are reviewed. The shielding process of ionospheric plasmas is outlined and the damping effect of the geomagnetic field on the rotation of artificial satellites is evaluated.

Opik, E. J.↗

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sensitivity Enhancement and Random Telegraph Noise in Magnetic Tunnel Junctions with Compensated Anisotropy

We demonstrate the manipulation of perpendicular magnetic anisotropy (PMA) in magnetic tunnel junctions (MTJs) by varying the layer stacks and temperature. PMA is tuned to compensate the shape anisotropy, giving a nonhysteretic magnetic response, a noteworthy sensitivity enhancement, and a field detectability of 1.8 nT/$\sqrt{Hz}$ at 100 kHz. Such a method is further exemplified in multiple MTJs, providing a solution to obtain desired sensitivities and operating temperatures. Additionally, the electronic noise of this MTJ is revealed as a random telegraph noise (RTN) due to a generation-recombination process. In conclusion, the observed voltage-dependent RTN could potentially be applied to true random number generators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Colorado Front Range Disasters: Understanding the Impact of Forest Management on the Cameron Peak and CalWood Fire

Along the Colorado Front Range, forest management has gained significant attention due to uncharacteristically large fires that burned late in 2020. The Cameron Peak Fire (largest in Colorado recorded history) and the CalWood Fire collectively burned an estimated 219,019 acres from August through December of 2020. Project partners at the Coalition for the Poudre River Watershed, Colorado State Forest Service, Ben Delatour Scout Ranch, The Nature Conservancy, Colorado Forest Restoration Institute, and Colorado State University were interested in understanding the effectiveness of previous forest treatments in reducing burn severity within the Cameron Peak Fire and the CalWood Fire. We first collated a forest treatment dataset from pre-existing datasets by reclassifying over 29,000 treatments, which occurred across the Northern Colorado Front Range between 1970-2020. Secondly, we mapped three burn severity indices using Landsat 8 OLI and Sentinel-2 MSI Earth observations and compared them to soil burn severity field data. Thirdly, a total of 35 topographic, disturbance, forest structure, and treatment predictor variables were generated across the fires. Finally, we assessed relationships between these predictor variables and burn severity using the random forest algorithm. Model results indicate that the primary drivers of burn severity were elevation and distance to treatment edge for the Cameron Peak Fire and fire area and forest canopy cover for the CalWood Fire. Further analysis of these variables paired with field data is necessary to understand the relationship between burn severity and treatments to guide future restoration efforts, improve forest resiliency, and mitigate fire risks.

Neal Swayze↗

Thermodynamics of order and randomness in dopant distributions inferred from atomically resolved imaging

Abstract Exploration of structure-property relationships as a function of dopant concentration is commonly based on mean field theories for solid solutions. However, such theories that work well for semiconductors tend to fail in materials with strong correlations, either in electronic behavior or chemical segregation. In these cases, the details of atomic arrangements are generally not explored and analyzed. The knowledge of the generative physics and chemistry of the material can obviate this problem, since defect configuration libraries as stochastic representation of atomic level structures can be generated, or parameters of mesoscopic thermodynamic models can be derived. To obtain such information for improved predictions, we use data from atomically resolved microscopic images that visualize complex structural correlations within the system and translate them into statistical mechanical models of structure formation. Given the significant uncertainties about the microscopic aspects of the material’s processing history along with the limited number of available images, we combine model optimization techniques with the principles of statistical hypothesis testing. We demonstrate the approach on data from a series of atomically-resolved scanning transmission electron microscopy images of Mo x Re 1- x S 2 at varying ratios of Mo/Re stoichiometries, for which we propose an effective interaction model that is then used to generate atomic configurations and make testable predictions at a range of concentrations and formation temperatures.

25 ENERGY STORAGE↗

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗

Laser particles with omnidirectional emission for cell tracking

The ability to track individual cells in space over time is crucial to analyzing heterogeneous cell populations. Recently, microlaser particles have emerged as unique optical probes for massively multiplexed single-cell tagging. However, the microlaser far-field emission is inherently direction-dependent, which causes strong intensity fluctuations when the orientation of the particle varies randomly inside cells. Here, we demonstrate a general solution based on the incorporation of nanoscale light scatterers into microlasers. Two schemes are developed by introducing either boundary defects or a scattering layer into microdisk lasers. The resulting laser output is omnidirectional, with the minimum-to-maximum ratio of the angle-dependent intensity improving from 0.007 (-24 dB) to > 0.23 (-6 dB). After transfer into live cells in vitro, the omnidirectional laser particles within moving cells could be tracked continuously with high signal-to-noise ratios for 2 h, while conventional microlasers exhibited frequent signal loss causing tracking failure.

59 BASIC BIOLOGICAL SCIENCES↗

Autonomous fabrication of tailored defect structures in 2D materials using machine learning-enabled scanning transmission electron microscopy

Materials with tailored quantum properties can be engineered from atomic-scale assembly techniques, but existing methods often lack the agility and accuracy to precisely and intelligently control the manufacturing process. Here, we demonstrate a fully autonomous approach for fabricating atomic-level defects using electron beams in scanning transmission electron microscopy (STEM) that combines advanced machine learning and automated beam control. As a proof of concept, we achieved controlled fabrication of MoS-nanowire (MoS-NW) edge structures by iterative and targeted exposure of MoS 2 monolayer to a focused electron beam to selectively eject sulfur atoms, utilizing high-angle annular dark-field (HAADF) imaging for feedback-controlled monitoring of structural evolution of defects. A machine learning framework combining a random forest model and a convolutional neural network (CNN) was developed to decode the HAADF image and accurately identify atomic positions and species. This atomic-level information was then integrated into an autonomous decision-making platform, which applied predefined fabrication strategies to instruct beam control about atomic sites to be ejected. The selected sites were subsequently exposed to a localized electron beam using an FPGA-controlled scan routine with precise control over beam positioning and duration. While the MoS-NW edge structures produced exhibit promising mechanical and electronic properties, the proposed methods to build the autonomous fabrication framework is material-agnostic and can be extended to other 2D materials for the creation of diverse defect structures and heterostructures beyond Mo S2 .

Engineering↗

Is Terzan 5 the remnant of a building block of the Galactic bulge? Evidence from APOGEE

ABSTRACT It has been proposed that the globular cluster-like system Terzan 5 is the surviving remnant of a primordial building block of the Milky Way bulge, mainly due to the age/metallicity spread and the distribution of its stars in the α–Fe plane. We employ Sloan Digital Sky Survey data from the Apache Point Observatory Galactic Evolution Experiment to test this hypothesis. Adopting a random sampling technique, we contrast the abundances of 10 elements in Terzan 5 stars with those of their bulge field counterparts with comparable atmospheric parameters, finding that they differ at statistically significant levels. Abundances between the two groups differ by more than 1σ in Ca, Mn, C, O, and Al, and more than 2σ in Si and Mg. Terzan 5 stars have lower [α/Fe] and higher [Mn/Fe] than their bulge counterparts. Given those differences, we conclude that Terzan 5 is not the remnant of a major building block of the bulge. We also estimate the stellar mass of the Terzan 5 progenitor based on predictions by the Evolution and Assembly of GaLaxies and their Environments suite of cosmological numerical simulations, concluding that it may have been as low as ∼3 × 108 M⊙ so that it was likely unable to significantly influence the mean chemistry of the bulge/inner disc, which is significantly more massive (∼1010 M⊙). We briefly discuss existing scenarios for the nature of Terzan 5 and propose an observational test that may help elucidate its origin.

79 ASTRONOMY AND ASTROPHYSICS↗

SDSS-IV MaNGA: The MaNGA Dwarf Galaxy Sample Presentation

We present the MaNGA Dwarf galaxy (MaNDala) Value Added Catalog (VAC), from the final release of the Sloan Digital Sky Survey-IV program. MaNDala consists of 136 randomly selected bright dwarf galaxies with M * < 10 9.1 M ⊙ and M g > -18.5, making it the largest integral field spectroscopy homogeneous sample of dwarf galaxies. We release a photometric analysis of the g, r, and z broadband imaging based on the DESI Legacy Imaging Surveys, as well as a spectroscopic analysis based on the Pipe3D SDSS-IV VAC. Our release includes the surface brightness (SB), geometric parameters, and color profiles, Sérsic fits as well as stellar population properties (such as stellar ages, metallicities, and star formation histories), and emission lines' fluxes within the FOV and the effective radii of the galaxies. We find that the majority of the MaNDala galaxies are star-forming late-type galaxies with $\langle$n Sersic,r $\rangle$ ~ 1.6 that are centrals (central/satellite dichotomy). MaNDala covers a large range of SB values (we find 11 candidate ultra-diffuse galaxies and three compact ones), filling the gap between classical dwarfs and low-mass galaxies in the Kormendy Diagram and in the size–mass/luminosity relation, which seems to flatten at 10 8 < M * /M ⊙ < 10 9 with $\langle$R e,r $\rangle$ ~ 2.7 kpc. A large fraction of MaNDala galaxies formed from an early low-metallicity burst of SF, but also from late SF events from more metal-enriched gas: half of the MaNDala galaxies assembled 50% of their mass at $\langle$z$\rangle$ > 2, while the last 20% was at $\langle$z$\rangle$ < 0.3. Finally, a bending of the sSFR-M * relation at M * ~ 10 9 M ⊙ for the main-sequence galaxies seems to be supported by MaNDala.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic stability

The field of stochastic stability is surveyed, with emphasis on the invariance theorems and their potential application to systems with randomly varying coefficients. Some of the basic ideas are reviewed, which underlie the stochastic Liapunov function approach to stochastic stability. The invariance theorems are discussed in detail.

Kushner, H. J.↗

New Instrumentation for Phase Partitioning

Cells and molecules can be purified by partitioning between the two immiscible liquid phases formed by aqueous solutions of poly/ethylene glycol and dextran. Such purification can be more selective, higher yielding, and less destructive to sensitive biological materials than other available techniques. Earth's gravitational field is a hindering factor as it causes sedimentation of particles to be purified and shear-induced particle randomization. The present proposal is directed toward developing new instrumentation for performing phase partitioning both on Earth and in microgravity.

Harris, J. M.↗