Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dimensionality Reduction”

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 217 records · Page 12

Goal-Oriented Optimal Approximations of Bayesian Linear Inverse Problems

We propose optimal dimensionality reduction techniques for the solution of goal-oriented linear-Gaussian inverse problems, where the quantity of interest (QoI) is a function of the inversion parameters. These approximations are suitable for large-scale applications. In particular, we study the approximation of the posterior covariance of the QoI as a low-rank negative update of its prior covariance and prove optimality of this update with respect to the natural geodesic distance on the manifold of symmetric positive definite matrices. Assuming exact knowledge of the posterior mean of the QoI, the optimality results extend to optimality in distribution with respect to the Kullback--Leibler divergence and the Hellinger distance between the associated distributions. We also propose the approximation of the posterior mean of the QoI as a low-rank linear function of the data and prove optimality of this approximation with respect to a weighted Bayes risk. Both of these optimal approximations avoid the explicit computation of the full posterior distribution of the parameters and instead focus on directions that are well informed by the data and relevant to the QoI. These directions stem from a balance among all the components of the goal-oriented inverse problem: prior information, forward model, measurement noise, and ultimate goals. We illustrate the theory using a high-dimensional inverse problem in heat transfer.

97 MATHEMATICS AND COMPUTING↗

Modeling and simulation challenges pursued by the Consortium for Advanced Simulation of Light Water Reactors (CASL)

The Consortium for the Advanced Simulation of Light Water Reactors (CASL), the first Energy Innovation Hub of the Department of Energy, was established in 2010 with the goal of providing modeling and simulation (M&S) capabilities that support and accelerate the improvement of nuclear energy's economic competitiveness and the reduction of spent nuclear fuel volume per unit energy, and all while assuring nuclear safety. To accomplish this requires advances in M&S capabilities in radiation transport, thermal-hydraulics, fuel performance and corrosion chemistry. To focus CASL's R&D, industry challenge problems have been defined, which equate with long standing issues of the nuclear power industry that M&S can assist in addressing. To date CASL has developed a multi-physics “core simulator” based upon pin-resolved radiation transport and subchannel (within fuel assembly) thermal-hydraulics, capitalizing on the capabilities of high performance computing. CASL's fuel performance M&S capability can also be optionally integrated into the core simulator, yielding a coupled multi-physics capability with untapped predictive potential. Material models have been developed to enhance predictive capabilities of fuel clad creep and growth, along with deeper understanding of zirconium alloy clad oxidation and hydrogen pickup. Understanding of corrosion chemistry (e.g., CRUD formation) has evolved at all scales: micro, meso and macro. CFD R&D has focused on improvement in closure models for subcooled boiling and bubbly flow, and the formulation of robust numerical solution algorithms. For multiphysics integration, several iterative acceleration methods have been assessed, illuminating areas where further research is needed. Finally, uncertainty quantification and data assimilation techniques, based upon sampling approaches, have been made more feasible for practicing nuclear engineers via R&D on dimensional reduction and biased sampling. Industry adoption of CASL's evolving M&S capabilities, which is in progress, will assist in addressing long-standing and future operational and safety challenges of the nuclear industry. - Highlights: • Complexity of physics based modeling of light water reactor cores being addressed. • Capability developed to help address problems that have challenged the nuclear power industry. • Simulation capabilities that take advantage of high performance computing developed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Disentangling Ferroelectric Wall Dynamics and Identification of Pinning Mechanisms via Deep Learning

Field-induced domain-wall dynamics in ferroelectric materials underpins multiple applications ranging from actuators to information technology devices and necessitates a quantitative description of the associated mechanisms including giant electromechanical couplings, controlled nonlinearities, or low coercive voltages. While the advances in dynamic piezoresponse force microscopy measurements over the last two decades have rendered visualization of polarization dynamics relatively straightforward, the associated insights into the local mechanisms have been elusive. This work explores the domain dynamics in model polycrystalline materials using a workflow combining deep-learning-based segmentation of the domain structures with nonlinear dimensionality reduction using multilayer rotationally invariant autoencoders (rVAE). The former allows unambiguous identification and classification of the ferroelectric and ferroelastic domain walls. Here, the rVAE discovers the latent representations of the domain wall geometries and their dynamics, consequently providing insight into the intrinsic mechanisms of polarization switching, that can further be compared to simple physical models. The rVAE disentangles the factors affecting the pinning efficiency of ferroelectric walls, offering insights into the correlation of ferroelastic wall distribution and ferroelectric wall pinning.

36 MATERIALS SCIENCE↗

Celestial Topology, Symmetry Theories, and Evidence for a NonSUSY D3‐Brane CFT

Symmetry Theories (SymThs) provide a flexible framework for analyzing the global categorical symmetries of a D -dimensional QFT D in terms of a (D + 1)-dimensional bulk system SymTh D+1 . In QFTs realized via local string backgrounds, these SymThs naturally arise from dimensional reduction of the linking boundary geometry. To track possible time dependent effects we introduce a celestial generalization of the standard “boundary at infinity” of a SymTh. As an application of these considerations we revisit large N quiver gauge theories realized by spacetime filling D3-branes probing a non-supersymmetric orbifold $\mathbb{R}$ 6 /Γ. Comparing the imprint of symmetry breaking on the celestial geometry at small and large ‘t Hooft coupling we find evidence for an intermediate symmetry preserving conformal fixed point.

conformal field theory↗

Bayesian Calibration of Stochastic Agent Based Model via Random Forest

Agent-based models (ABM) provide an excellent framework for modeling outbreaks and interventions in epidemiology by explicitly accounting for diverse individual interactions and environments. However, these models are usually stochastic and highly parametrized, requiring precise calibration for predictive performance. When considering realistic numbers of agents and properly accounting for stochasticity, this high-dimensional calibration can be computationally prohibitive. This paper presents a random forest-based surrogate modeling technique to accelerate the evaluation of ABMs and demonstrates its use to calibrate an epidemiological ABM named CityCOVID via Markov chain Monte Carlo (MCMC). The technique is first outlined in the context of CityCOVID's quantities of interest, namely hospitalizations and deaths, by exploring dimensionality reduction via temporal decomposition with principal component analysis (PCA) and via sensitivity analysis. The calibration problem is then presented, and samples are generated to best match COVID-19 hospitalization and death numbers in Chicago from March to June in 2020. Further, these results are compared with previous approximate Bayesian calibration (IMABC) results, and their predictive performance is analyzed, showing improved performance with a reduction in computation.

60 APPLIED LIFE SCIENCES↗

Statistical learning for predicting density–matrix-based electron dynamics

In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.

59 BASIC BIOLOGICAL SCIENCES↗

Black hole spectroscopy and AdS 2 holography

We compute the spectrum of extremal nonBPS black holes in four dimensions by studying supergravity on their AdS 2 × S2 near horizon geometry. We find that the spectrum exhibits significant simplifications even though supersymmetry is completely broken. We interpret our results in the framework of nAdS 2 /nCFT 1 correspondence and by comparing with dimensional reduction from AdS 3 /CFT 2 duality. As an additional test we compute quantum corrections to extremal black hole entropy on the nonBPS branch and recover results previously determined using very different methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Light-front representation of chiral dynamics with Δ isobar and large-N c relations

Transverse densities describe the spatial distribution of electromagnetic current in the nucleon at fixed light-front time. At peripheral distances b = O(M π –1 ) the densities are governed by chiral dynamics and can be calculated model-independently using chiral effective field theory (EFT). Recent work has shown that the EFT results can be represented in first-quantized form, as overlap integrals of chiral light-front wave functions describing the transition of the nucleon to soft-pion-nucleon intermediate states, resulting in a quantum-mechanical picture of the peripheral transverse densities. We now extend this representation to include intermediate states with Δ isobars and implement relations based on the large-N c limit of QCD. We derive the wave function overlap formulas for the Δ contributions to the peripheral transverse densities by way of a three-dimensional reduction of relativistic chiral EFT expressions. Our procedure effectively maintains rotational invariance and avoids the ambiguities with higher-spin particles in the light-front time-ordered approach. We study the interplay of πN and πΔ intermediate states in the quantum-mechanical picture of the densities in a transversely polarized nucleon. We show that the correct N c -scaling of the charge and magnetization densities emerges as the result of the particular combination of currents generated by intermediate states with degenerate N and Δ. The off-shell behavior of the chiral EFT is summarized in contact terms and can be studied easily. As a result, the methods developed here can be applied to other peripheral densities and to moments of the nucleon's generalized parton distributions.

1/N expansion↗

The sub-leading scattering waveform from amplitudes

We compute the next-to-leading order term in the scattering waveform of uncharged black holes in classical general relativity and of half-BPS black holes in $\mathcal{N}$ = 8 supergravity. We propose criteria, generalizing explicit calculations at next-to-leading or der, for determining the terms in amplitudes that contribute to local observables. For general relativity, we construct the relevant classical integrand through generalized unitary in two distinct ways, (1) in a heavy-particle effective theory and (2) in general relativity minimally-coupled to scalar fields. With a suitable prescription for the matter propagator in the former, we find agreement between the two methods, thus demonstrating the absence of interference of quantum and classically-singular contributions. The classical $\mathcal{N}$ = 8 integrand for massive scalar fields is constructed through dimensional reduction of the known five-point one-loop integrand. Our calculation exhibits novel features compared to conservative calculations and inclusive observables, such as the appearance of master integrals with intersecting matter lines and the appearance of a classical infrared divergence whose absence from classical observables requires a suitable definition of the retarded time.

79 ASTRONOMY AND ASTROPHYSICS↗

Islands in non-minimal dilaton gravity: exploring effective theories for black hole evaporation

We start from (3 + 1)-dimensional Einstein gravity with minimally coupled massless scalar matter, through spherical dimensional reduction, the matter theory is non-minimally coupled with the dilaton in (1 + 1)-dimensions. Despite its simplicity, constructing a self-consistent one-loop effective theory for this model remains a challenge, partially due to a Weyl-invariant ambiguity in the effective action. With a universal splitting property for the one-loop action, the ambiguity can be identified with the state-dependent part of the covariant quantum stress tensor. By introducing on-shell equivalent auxiliary fields to construct minimal candidates of Weyl-invariant terms, we derive a one-parameter family of one-loop actions with unique, regular, and physical stress tensors corresponding to the Boulware, Hartle-Hawking and Unruh states. We further study the back-reacted geometry and the corresponding quantum extremal islands that were inaccessible without a consistent one-loop theory. Along the way, we elaborate on the implications of our construction for the non-minimal dilaton gravity model.

2D Gravity↗

Scaling similarities and quasinormal modes of D0 black hole solutions

We study the gravity solution dual to the D0 brane quantum mechanics, or BFSS matrix model, in the ’t Hooft limit. The classical physics described by this gravity solution is invariant under a scaling transformation, which changes the action with a specific critical exponent, sometimes called the hyperscaling violating exponent. We present an argument for this critical exponent from the matrix model side, which leads to an explanation for the peculiar temperature dependence of the entropy in this theory, S ∝ T 9/5 . We also present a similar argument for all other Dp-brane geometries. We then compute the black hole quasinormal modes. This involves perturbing the finite temperature geometry. These perturbations can be easily obtained by a mathematical trick where we view the solution as the dimensional reduction of an Ad S2+9/5 × S 8 geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

M5-brane and D-brane scattering amplitudes

We present tree-level $\mathcal{n}$-particle on-shell scattering amplitudes of various brane theories with 16 conserved supercharges. These include the world-volume theory of a probe D3-brane or D5-brane in 10D Minkowski spacetime as well as a probe M5-brane in 11D Minkowski spacetime, which describes self interactions of an abelian tensor supermultiplet with 6D (2, 0) supersymmetry. Twistor-string-like formulas are proposed for tree-level scattering amplitudes of all multiplicities for each of these theories. The R symmetry of the D3-brane theory is shown to be SU(4) × U(1), and the U(1) factor implies that its amplitudes are helicity conserving. Each of 6D theories (D5-brane and M5-brane) reduces to the D3-brane theory by dimensional reduction. Conclusively, as special cases of the general M5-brane amplitudes, we present compact formulas for examples involving only the self-dual B field with $\mathcal{n}$ = 4, 6, 8.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sharpening the Distance Conjecture in diverse dimensions

The Distance Conjecture holds that any infinite-distance limit in the scalar field moduli space of a consistent theory of quantum gravity must be accompanied by a tower of light particles whose masses scale exponentially with proper field distance ||$\phi$|| as m ~ exp(–λ||$\phi$||), where λ is order-one in Planck units. While the evidence for this conjecture is formidable, there is at present no consensus on which values of λ are allowed. In this paper, we propose a sharp lower bound for the lightest tower in a given infinite-distance limit in d dimensions: λ ≥ $1/\sqrt{d-2}$. In support of this proposal, we show that (1) it is exactly preserved under dimensional reduction, (2) it is saturated in many examples of string/M-theory compactifications, including maximal supergravity in d = 4 – 10 dimensions, and (3) it is saturated in many examples of minimal supergravity in d = 4 – 10 dimensions, assuming appropriate versions of the Weak Gravity Conjecture. We argue that towers with λ < $1/\sqrt{d-2}$ discussed previously in the literature are always accompanied by even lighter towers with λ ≥ $1/\sqrt{d-2}$, thereby satisfying our proposed bound. We discuss connections with and implications for the Emergent String Conjecture, the Scalar Weak Gravity Conjecture, the Repulsive Force Conjecture, large-field inflation, and scalar field potentials in quantum gravity. In particular, we argue that if our proposed bound applies beyond massless moduli spaces to scalar fields with potentials, then accelerated cosmological expansion cannot occur in asymptotic regimes of scalar field space in quantum gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗

Equivariant Verlinde Formula from Fivebranes and Vortices

We study complex Chern–Simons theory on a Seifert manifold M 3 by embedding it into string theory. We show that complex Chern–Simons theory on M 3 is equivalent to a topologically twisted supersymmetric theory and its partition function can be naturally regularized by turning on a mass parameter. We find that the dimensional reduction of this theory to 2d gives the low energy dynamics of vortices in fourdimensional gauge theory, the fact apparently overlooked in the vortex literature. We also generalize the relations between (1) the Verlinde algebra, (2) quantum cohomology of the Grassmannian, (3) Chern–Simons theory on Σ × S 1 and (4) index of a spin c Dirac operator on the moduli space of flat connections to a new set of relations between (1) the “equivariant Verlinde algebra” for a complex group, (2) the equivariant quantum K-theory of the vortex moduli space, (3) complex Chern–Simons theory on Σ × S 1 and (4) the equivariant index of a spinc Dirac operator on the moduli space of Higgs bundles.

97 MATHEMATICS AND COMPUTING↗

Toward memory-efficient melt pool monitoring: a classification framework using event-based imaging and sparse sensing technique

Vision sensors like CMOS and CCD cameras are often used for in-process monitoring of melt pools in laser-based additive and welding processes, but they require transferring large amounts of data and computational processing resources. Event-based neuromorphic imagery, on the other hand, detects only the change in pixel intensity, thus potentially reducing the data amount and latency. With an event imager, this study develops a framework for melt pool condition classification, including image construction, time scale selection, optimal pixel selection, and sparse classification, to achieve a highly memory-efficient scheme. These are based on sparse sensing techniques with singular value decomposition (SVD) and QR pivoting, the two fundamental matrix transformations for linear dimensionality reduction. The framework is then validated by classifying a controlled experiment by exciting various mode shapes of liquid gallium pools of varying depths (3, 6, and 8 mm). At 200 pixels, the classifier can reach overall accuracy of 75%, while at 2000 pixels (0.013% of the total possible pixels), the accuracy is nearly 90% (89.86%). At the same number of pixels, random selection can only achieve 46% and 67%, respectively. The memory savings of the sparsely sampled event data compared to a conventional imager is about 500 times. In addition to performance, implementation and limitations of the framework are also discussed.

42 ENGINEERING↗

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗

Scaling Limits and Generic Bounds for Exploration Processes

We consider exploration algorithms of the random sequential adsorption type both for homogeneous random graphs and random geometric graphs based on spatial Poisson processes. At each step, a vertex of the graph becomes active and its neighboring nodes become blocked. Given an initial number of vertices N growing to infinity, we study statistical properties of the proportion of explored (active or blocked) nodes in time using scaling limits. We obtain exact limits for homogeneous graphs and prove an explicit central limit theorem for the final proportion of active nodes, known as the jamming constant, through a diffusion approximation for the exploration process which can be described as a unidimensional process. We then focus on bounding the trajectories of such exploration processes on random geometric graphs, i.e., random sequential adsorption. As opposed to exploration processes on homogeneous random graphs, these do not allow for such a dimensional reduction. Instead we derive a fundamental relationship between the number of explored nodes and the discovered volume in the spatial process, and we obtain generic bounds for the fluid limit and jamming constant: bounds that are independent of the dimension of space and the detailed shape of the volume associated to the discovered node. Lastly, using coupling techinques, we give trajectorial interpretations of the generic bounds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗