Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical manifold method”

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.

82 records · Page 5

A long short-term memory embedding for hybrid uplifted reduced order models

In this paper, we introduce an uplifted reduced order modeling (UROM) approach through the integration of standard projection based methods with long short-term memory (LSTM) embedding. Our approach has three modeling layers or components. In the first layer, we utilize an intrusive projection approach to model dynamics represented by the largest modes. The second layer consists of an LSTM model to account for residuals beyond this truncation. This closure layer refers to the process of including the residual effect of the discarded modes into the dynamics of the largest scales. However, the feasibility of generating a low rank approximation tails off for higher Kolmogorov n -width systems due to the underlying nonlinear processes. The third uplifting layer, called super-resolution, addresses this limited representation issue by expanding the span into a larger number of modes utilizing the versatility of LSTM. Therefore, our model integrates a physics-based projection model with a memory embedded LSTM closure and an LSTM based super-resolution model. In several applications, we exploit the use of Grassmann manifold to construct UROM for unseen conditions. We performed numerical experiments by using the Burgers and Navier-Stokes equations with quadratic nonlinearity. Finally, our results show robustness of the proposed approach in building reduced order models for parameterized systems and confirm the improved trade-off between accuracy and efficiency.

42 ENGINEERING↗

Fluid & Thermal Analysis of a Manifold Microchannel Heat Sink

A three-dimensional (3-D) numerical model has been developed to study the fluid flow and heat transfer through a manifold-microchannel (MMC) heat sink, a microfluidic cooling system designed as a heat dissipation solution for compact, high-powered electronic systems. An MMC consists of a 3-D manifold structure that distributes working fluid through alternating inlet and outlet manifold channels that guide the fluid to and from the microchannels, shortening the flow length and thereby, reducing the pressure drop across the system. Wide bandgap semiconductors such as silicon carbide (SiC) are increasingly utilized due to their ability to enable smaller electronic devices and components to operate more efficiently at higher frequencies. These semiconductors can greatly benefit from the MMC, which can extract all the heat produced by compact, high-powered electronics, enabling them to achieve optimal performance. The MMC also utilizes an embedded cooling approach which allows for fabrication directly within SiC-based electronics, resulting in significant mass and weight savings and lower thermal resistance than conventional heat sink designs. This paper outlines the methods used to conduct fluid and thermal analysis of different MMC configurations for flow parameters such as flow rate and inlet flow conditions. This study involves both steady state and transient thermal analysis as a gradually increasing heat flux is applied to the system over time, which is representative of the heat dissipation from the electronics over a given period. Results from the analysis show that the MMC is capable of effectively dissipating heat flux of ~100 W/cm2 while maintaining considerably low pressure drop (< 1 kPa) and thermal resistance (< 1 K/W) for certain single-phase working fluids.

manifold microchannel↗

Fluid & Thermal Analysis of a Manifold Microchannel Heat Sink

A three-dimensional (3-D) numerical model has been developed to study the fluid flow and heat transfer through a manifold-microchannel (MMC) heat sink, a microfluidic cooling system designed as a heat dissipation solution for compact, high-powered electronic systems. An MMC consists of a 3-D manifold structure that distributes working fluid through alternating inlet and outlet manifold channels that guide the fluid to and from the microchannels, shortening the flow length and thereby, reducing the pressure drop across the system. Wide bandgap semiconductors such as silicon carbide (SiC) are increasingly utilized due to their ability to enable smaller electronic devices and components to operate more efficiently at higher frequencies. These semiconductors can greatly benefit from the MMC, which can extract all the heat produced by compact, high-powered electronics, enabling them to achieve optimal performance. The MMC also utilizes an embedded cooling approach which allows for fabrication directly within SiC-based electronics, resulting in significant mass and weight savings and lower thermal resistance than conventional heat sink designs. This paper outlines the methods used to conduct fluid and thermal analysis of different MMC configurations for flow parameters such as flow rate and inlet flow conditions. This study involves both steady state and transient thermal analysis as a gradually increasing heat flux is applied to the system over time, which is representative of the heat dissipation from the electronics over a given period. Results from the analysis show that the MMC is capable of effectively dissipating heat flux of ~100 W/cm2 while maintaining considerably low pressure drop (< 1 kPa) and thermal resistance (< 1 K/W) for certain single-phase working fluids.

manifold microchannel↗

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗

Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

Here, we develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems, including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input–output and the hierarchical equations of motion formalisms into a double hierarchy of density matrix equations. We show analytically that under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an N-photon Fock state input induces no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon with or without the influence of phonons. For short pulses (or, equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δ and an exciton–light coupling constant Γ. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyze the absorption probability in the opposite long pulse (narrow bandwidth) regime. We then derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective vs independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy equations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Broadband cross polarization for ultra-wideline magic-angle spinning NMR

Over the past decade, there has been a sustained interest in using frequency-swept (FS) pulses for the efficient acquisition of wideline and ultra-wideline (UW) NMR powder patterns. Such experiments are typically conducted under static conditions, employing both direct- and indirect-excitation methods (i.e., WCPMG and BRAIN-CP/WCPMG, respectively). Recently, Koppe et al. demonstrated that the WCPMG pulse sequence can be used to efficiently acquire wideline and UW NMR spectra with spinning sideband (SSB) manifolds under magic-angle spinning (MAS) conditions, capitalizing on the increased signal-to-noise ratios (SNR) afforded by MAS. To date, there have been only a few instances of broadband cross-polarization (CP) experiments using FS pulses under MAS conditions and no applications to systems exhibiting wideline and/or ultra-wideline powder patterns, despite the clear advantages these experiments could offer. Herein, we demonstrate that FS pulses selectively applied to a single sideband of the S spin can be used for efficient 1 H-S polarization transfer to S = 1/2 nuclides with large anisotropic chemical shift interactions at slow to moderate MAS rates. The Hartmann–Hahn matching conditions in BRAIN-CP/WCPMG-MAS experiments bear similarity to those of standard CP sequences, yet operate over UW frequency ranges and only require low-amplitude RF pulses on the S channel. Crucial to the success of the BRAIN-CP/WCPMG-MAS experiment is careful calibration of the RF amplitude, transmitter offset, and effective frequency sweep of the FS pulse applied to the S spins at a given MAS rate. Thus, by means of numerical simulations and experimental testing, we provide recommendations for the parameterization and setup of BRAIN-CP/WCPMG-MAS experiments for their most efficient use. Results showcasing the capability of the BRAIN-CP/WCPMG-MAS pulse sequence are presented, including applications to 119 Sn, 195 Pt, and 103 Rh NMR.

Kimball, James J. [Florida State University, Talla↗

Deep learning to generate in silico chemical property libraries and candidate molecules for small molecule identification in complex samples

Comprehensive and unambiguous identification of small molecules in complex samples will revolutionize our understanding of the role of metabolites in biological systems. Existing and emerging technologies have enabled measurement of chemical properties of molecules in complex mixtures and, in concert, are sensitive enough to resolve even stereoisomers. Despite these experimental advances, small molecule identification is inhibited by (i) chemical reference libraries (e.g. mass spectra, collision cross section, and other measurable property libraries) representing <1% of known molecules, limiting the number of possible identifications, and (ii) the lack of a method to generate candidate matches directly from experimental features (i.e. without a library). To this end, we developed a variational autoencoder (VAE) to learn a continuous numerical, or latent, representation of molecular structure to expand reference libraries for small molecule identification. We extended the VAE to include a chemical property decoder, trained as a multitask network, in order to shape the latent representation such that it assembles according to desired chemical properties. The approach is unique in its application to metabolomics and small molecule identification, with its focus on properties that can be obtained from experimental measurements (m/z, CCS) paired with its training paradigm, which involved a cascade of transfer learning iterations. First, molecular representation is learned from a large dataset of structures with m/z labels. Next, in silico property values are used to continue training, as experimental property data is limited. Finally, the network is further refined by being trained with the experimental data. This allows the network to learn as much as possible at each stage, enabling success with progressively smaller datasets without overfitting. Once trained, the network can be used to predict chemical properties directly from structure, as well as generate candidate structures with desired chemical properties. Our approach is orders of magnitude faster than first-principles simulation for CCS property prediction. Additionally, the ability to generate novel molecules along manifolds, defined by chemical property analogues, positions DarkChem as highly useful in a number of application areas, including metabolomics and small molecule identification, drug discovery and design, chemical forensics, and beyond.

Colby, Sean M.↗

Dynamics of a class of vortex rings

The contour dynamics method is extended to vortex rings with vorticity varying linearly from the symmetry axis. An elliptic core model is also developed to explain some of the basic physics. Passage and collisions of two identical rings are studied focusing on core deformation, sound generation and stirring of fluid elements. With respect to core deformation, not only the strain rate but how rapidly it varies is important and accounts for greater susceptibility to vortex tearing than in two dimensions. For slow strain, as a passage interaction is completed and the strain relaxes, the cores return to their original shape while permanent deformations remain for rapidly varying strain. For collisions, if the strain changes slowly the core shapes migrate through a known family of two-dimensional steady vortex pairs up to the limiting member of the family. Thereafter energy conservation does not allow the cores to maintain a constant shape. For rapidly varying strain, core deformation is severe and a head-tail structure in good agreement with experiments is formed. With respect to sound generation, good agreement with the measured acoustic signal for colliding rings is obtained and a feature previously thought to be due to viscous effects is shown to be an effect of inviscid core deformation alone. For passage interactions, a component of high frequency is present. Evidence for the importance of this noise source in jet noise spectra is provided. Finally, processes of fluid engulfment and rejection for an unsteady vortex ring are studied using the stable and unstable manifolds. The unstable manifold shows excellent agreement with flow visualization experiments for leapfrogging rings suggesting that it may be a good tool for numerical flow visualization in other time periodic flows.

Shariff, Karim↗

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗