Search NASA⌕ Search

SEARCH · Search NASA

Results for “rational approximation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Expanded capabilities of the CarMa code in modeling resistive wall mode dynamics with 3-D conductors

In this work, an improved version of the CarMa code is presented, called CarMa-D, for the analysis of resistive wall modes (RWMs) in fusion devices, simultaneously considering the effects of volumetric three-dimensional conducting structures, plasma dynamics, toroidal rotation or drift-kinetic damping. Unlike static CarMa, CarMa-D does not rely on the simplifying assumptions such as neglecting the plasma mass, toroidal rotation and kinetic damping. The new coupling strategy is based on matrix-based Padé rational functions approximation of plasma a response. The arising mathematical model is formally equivalent to the original CarMa model, but with a higher number of degrees of freedom to model the dynamics of the plasma. CarMa-D overcomes the main limitations of the original CarMa, in particular: (i) the massless assumption for the plasma is removed, allowing modeling of global modes growing both on ideal kink time scales and in the typical RWM growth rate regime, with a suitable treatment of the model; (ii) the effects of toroidal plasma flow and drift kinetic damping can be included into the new model, providing a powerful tool to study macroscopic phenomena where both plasma dynamics and 3-D conducting structures play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Biexcitons are bound in CsPbBr 3 perovskite nanocrystals

Here, we study the energetics of quasiparticle excitations in CsPbBr 3 perovskite nanocrystals using path integral molecular dynamics simulations. Employing detailed molecular models, we elucidate the interplay of anharmonic lattice degrees of freedom, dielectric confinement, and electronic correlation on exciton and biexciton binding energies of a range of nanocrystal sizes. We find generally good agreement with some experimental observations of binding energies and additionally explain the observed size-dependent Stokes shift. The explicit model calculations are compared with simplified approximations to rationalize the lattice contributions to binding. We find that polaron formation significantly reduces exciton binding energies, whereas these effects are negligible for biexciton interactions. While experimentally the binding energy of biexcitons is uncertain, based on our study, we conclude that biexcitons are bound in CsPbBr 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Computational Math Problems for a Clean Energy Future

Cutting edge computational mathematics are ubiquitous in renewable energy research. Problems in resilient and reliable electric grid operations, infrastructure planning, wind farm yaw control, and more demand sophisticated and scalable computational tools that enable the transition of renewable energy technologies from proof of concept to deployment into our energy system. The mission of the Computational Science Center at NREL is to lead the lab's efforts to solve energy challenges using high-performance computing (HPC), computational science, applied mathematics, scientific data management, visualization, and informatics. In this poster, we provide a short overview of three areas of computational mathematics research at NREL: wind power scenario generation for stochastic grid operations and infrastructure planning, improved rational function approximations for electromagnetic transients codes, and wind farm yaw control using a combination of the Alternating Direction Method of Multipliers (ADMM) and reinforcement learning (RL). Increasing penetrations of renewable energy into power grids motivate the investigation of new approaches to characterizing uncertainty for five-minute economic dispatch problems. Similarly, as the penetration of distributed energy resources on power grids increases, it becomes important to revisit our methods of modelling transient phenomena, i.e. electromagnetic transients programs. Finally, the combination of ADMM and RL for wind farm yaw control presented here can potentially increase the efficiency of the deployed distributed controllers by orders of magnitude.

ADMM↗
DOE OSTI2022

Ingrained: An Automated Framework for Fusing Atomic-Scale Image Simulations into Experiments

To fully leverage the power of image simulation to corroborate and explain patterns and structures in atomic resolution microscopy, an initial correspondence between the simulation and experimental image must be established at the outset of further high accuracy simulations or calculations. Furthermore, if simulation is to be used in context of highly automated processes or high-throughput optimization, the process of finding this correspondence itself must be automated. In this work, "ingrained," an open-source automation framework which solves for this correspondence and fuses atomic resolution image simulations into the experimental images to which they correspond, is introduced. Here, the overall "ingrained" workflow, focusing on its application to interface structure approximations, and the development of an experimentally rationalized forward model for scanning tunneling microscopy simulation are described.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Strong-Weak Duality via Jordan-Wigner Transformation: Using Fermionic Methods for Strongly Correlated su(2) Spin Systems

The Jordan-Wigner transformation establishes a duality between $su(2)$ and fermionic algebras. We present qualitative arguments and numerical evidence that when mapping spins to fermions, the transformation makes strong correlation weaker, as demonstrated by the Hartree-Fock approximation to the transformed Hamiltonian. This result can be rationalized in terms of rank reduction of spin shift terms when transformed to fermions. Conversely, the mapping of fermions to qubits makes strong correlation stronger, complicating its solution when one uses qubit-based correlators. The presence of string operators poses challenges to the implementation of quantum chemistry methods on classical computers, but these can be dealt with using established techniques of low computational cost. Here, our proof of principle results for XXZ and J1-J2 Heisenberg (in 1D and 2D) indicate that the JW transformed fermionic Hamiltonian has reduced complexity in key regions of their phase diagrams, and provides a better starting point for addressing challenging spin problems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Water Uptake by Silica Nanopores: Impacts of Surface Hydrophilicity and Pore Size

Water uptake at 298 K by two mesoporous silicas with different pore sizes was studied using volumetric vapor sorption. Here, through variation of sample pretreatment temperature and time, the number of surface hydroxyl groups was varied, leading to marked changes in the water sorption behavior. With increasing surface hydroxyl density, the adsorption of water to the pore surfaces increases, leading to formation of a thicker and denser water sorption layer. The BET model was used to measure surface hydroxyl density from low pressure parts of water adsorption isotherms. Adsorption of water molecules scales with the number density of surface hydroxyl groups, and monolayer coverage is reached at reduced pressure of ca. 0.3. The onset of adsorption pore condensation of water shifts to lower reduced pressure with increasing surface hydroxyl density, indicating growing thickness of adsorption layers. However, the water pore condensation step of the desorption branch shifts to smaller reduced pressure with increasing surface hydroxylation, reducing the adsorption-desorption hysteresis width. An adsorbed phase model is utilized to calculate the distribution of pore water molecules between adsorbed and pore condensation phases, and characterize the density and thickness of the water sorption phase. In the smaller pores, the adsorbed phase forms a sparse monolayer, while approximately a bilayer is formed in the wider pores. The previously reported silica nanopore underfilling by water is confirmed and rationalized by a significantly reduced water density in the adsorbed phase with respect to bulk, while the pore core is filled with water of approximately bulk density.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Isogeometric large-eddy simulations of turbulent particle-laden flows

In recent years, isogeometric analysis (IGA) has attracted significant attention from the computational mechanics community due to its ability to integrate design and analysis. Besides, IGA is also a higher-order discretization technique for solving partial differential equations, showing high approximation capability per degree of freedom. In this paper, we extend the application realm of IGA to particle-laden flows based on Eulerian–Eulerian description that couples Navier–Stokes equations with a density transport equation through a Boussinesq approximation. The coupled systems are solved by using quadratic non-uniform rational B-spline (NURBS) functions and a recently developed residual-based variational multiscale (VMS) formulation, which introduces coupling between the fine velocity scales and density equation residuals. We deploy the proposed approach to perform large-eddy simulations (LES) of dilute particle-laden flows over a flat surface at Reynolds number = 10,000. We compare the simulation results against direct numerical simulation (DNS) results from the literature. We find that combining VMS and IGA, the proposed approach enables accurate prediction of a wide range of flow/particle statistics with a relatively lower mesh resolution.

Mathematics↗
DOE OSTI2020

IGA-MPM: The Isogeometric Material Point Method

In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.

42 ENGINEERING↗
DOE OSTI2026

Surface Reconstruction in Hydrated Amphiphilic Block Copolymer Thin Films Probed by Fluid Cell Atomic Force Microscopy

In many thin film materials, nuanced interplays of interfacial energies control the surface morphology and rearrangement. This work evaluates polymer−solvent interactions and solvent-driven surface reconstructions via ex situ and in situ fluid cell Atomic Force Microscopy (fc-AFM) analysis of amphiphilic block copolymer (BCP) thin films upon exposure to deionized (DI) water. We examine the differences in surface morphology, whole-film swelling, and force response in thin films of polystyrene-block-poly(ethylene oxide) (PS-b-PEO) and polystyrene- block-poly[(allyl glycidyl ether)-co-(ethylene oxide)] (PS-b- P[AGE-co-EO]) processed into standing-up cylinder morphologies perpendicular to a silicon substrate (⊥C). Using Amplitude Modulation AFM (AM-AFM) and Amplitude-Phase Distance (APD) force spectroscopy, this work probes the mechanoresponsive nature of the dynamic surface layers of these films, unveiling surface layer stratification and surface chain rearrangement via minimal tip−sample stimulation. To help rationalize the observed reconfigurations, the energetic driving forces were estimated using the harmonic mean approximations of interfacial energies. Given the nonionizable nature of the minority P(AGE-co-EO) block and the energetic driving forces for chain mobility, this work shows how the elimination of unfavorable PS−water interfaces drives chain rearrangement and coverage of the PS surface by chains of the hydrophilic block. This work highlights considerations for increasing the heterogeneity and complexity of BCP thin films via random blocks and how those changes to local interfacial energies may drive larger scale film morphology reconstructions, with broader implications for tuning interface hydrophilicity.

Copolymers↗
DOE OSTI2025

Precision Reconstruction of Rational Conformal Field Theory from Exact Fixed-Point Tensor Network

The novel concept of entanglement renormalization and its corresponding tensor network renormalization technique have been highly successful in developing a controlled real-space renormalization group (RG) scheme. Numerically approximate fixed-point (FP) tensors are widely used to extract the conformal data of the underlying conformal field theory (CFT) describing critical phenomena. In this paper, we present an explicit analytical construction of the FP tensor for 2D rational CFT. We define it as a correlation function between the “boundary-changing operators” (BCO) on triangles. Our construction fully captures all the real-space RG conditions. We also provide concrete examples, such as Ising, Yang-Lee, and tricritical Ising models, to compute the scaling dimensions explicitly based on the corresponding FP tensor. The BCO descendants turn out to be an optimal basis such that truncation in bond dimensions naturally produces comparable accuracies with the leading existing FP algorithms. Interestingly, our construction of FP tensors is closely related to a strange correlator, where the holographic picture naturally emerges. Our results also open a new door toward understanding CFT in higher dimensions. Published by the American Physical Society 2025

Cheng, Gong (ORCID:0009000891587404)↗
DOE OSTI2023

Mutations in the Second Alternative Oxidase Gene: A New Approach to Group Aspergillus niger Strains

Alternative oxidase is a terminal oxidase in the branched mitochondrial electron transport chain of most fungi including Aspergillus niger (subgenus Circumdati, section Nigri). A second, paralogous aox gene (aoxB) is extant in some A. niger isolates but also present in two divergent species of the subgenus Nidulantes—A. calidoustus and A. implicatus—as well as in Penicillium swiecickii. Black aspergilli are cosmopolitan opportunistic fungi that can cause diverse mycoses and acute aspergillosis in immunocompromised individuals. Amongst the approximately 75 genome-sequenced A. niger strains, aoxB features considerable sequence variation. Five mutations were identified that rationally affect transcription or function or terminally modify the gene product. One mutant allele that occurs in CBS 513.88 and A. niger neotype strain CBS 554.65 involves a chromosomal deletion that removes exon 1 and intron 1 from aoxB. Another aoxB allele results from retrotransposon integration. Three other alleles result from point mutations: a missense mutation of the start codon, a frameshift, and a nonsense mutation. A. niger strain ATCC 1015 has a full-length aoxB gene. The A. niger sensu stricto complex can thus be subdivided into six taxa according to extant aoxB allele, which may facilitate rapid and accurate identification of individual species

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Parsimonious Inference Information-Theoretic Foundations for a Complete Theory of Machine Learning (CIS-LDRD Project 218313 Final Technical Report)

This work examines how we may cast machine learning within a complete Bayesian framework to quantify and suppress explanatory complexity from first principles. Our investigation into both the philosophy and mathematics of rational belief leads us to emphasize the critical role of Bayesian inference in learning well-justified predictions within a rigorous and complete extended logic. The Bayesian framework allows us to coherently account for evidence in the learned plausibility of potential explanations. As an extended logic, the Bayesian paradigm regards probability as a notion of degrees of truth. In order to satisfy critical properties of probability as a coherent measure, as well as maintain consistency with binary propositional logic, we arrive at Bayes' Theorem as the only justifiable mechanism to update our beliefs to account for empiracle evidence. Yet, in the machine learning paradigm, where explanations are unconstrained algorithmic abstractions, we arrive at a critical challenge: Bayesian inference requires prior belief. Conventional approaches fail to yield a consistent framework in which we could compare prior plausibility among the infinities of potential choices in learning architectures. The difficulty of articulating well-justified prior belief over abstract models is the provinence of memorization in traditional machine learning training practices. This becomes exceptionally problematic in the context of limited datasets, when we wish to learn justifiable predictions from only a small amount of data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Accurate Kohn-Sham auxiliary system from the ground-state density of solids

The Kohn-Sham (KS) system is an auxiliary system whose effective potential is unknown in most cases. It is in principle determined by the ground-state density and it has been found numerically for some low-dimensional systems by inverting the KS equations starting from a given accurate density. For solids, only approximate results are available. In this work, we determine accurate exchange-correlation (xc) potentials for Si and NaCl using the ground-state densities obtained from auxiliary field quantum Monte Carlo calculations. We show that these xc potentials can be rationalized as an ensemble of a few local functions of the density, whose form depends on the specific environment and can be well characterized by the gradient of the density and the local kinetic energy density. Further, the KS band structure can be obtained with high accuracy. The true KS band gap turns out to be larger than the prediction of the local density approximation, but significantly smaller than the measurable photoemission gap, which confirms previous estimates. Finally, our findings show that the conjecture that very different xc potentials can lead to very similar densities and other KS observables is true also in solids, which questions the meaning of details of the potentials and, at the same time, confirms the stability of the KS system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb–exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.

Chemistry↗
DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Competing Intermolecular and Molecule–Surface Interactions: Dipole–Dipole-Driven Patterns in Mixed Carborane Self-Assembled Monolayers

Carboranedithiol isomers adsorbing with opposite orientations of their dipoles on surfaces are self-assembled together to form mixed monolayers where both lateral dipole–dipole and lateral thiol–thiolate (S–H···S) interactions provide enhanced stability over single-component monolayers. Here, we demonstrate the first instance of the ability to map individual isomers in a mixed monolayer using the model system carboranedithiols on Au{111}. The addition of methyl groups to one isomer provides both an enhanced dipole moment and extra apparent height for differentiation via scanning tunneling microscopy (STM). Associated computational investigations rationalize favorable interactions of mixed pairs and the associated stability changes that arise from these interactions. Both STM images and Monte Carlo simulations yield similarly structured mixed monolayers, where approximately 10% of the molecules have reversed dipole moment orientations but no direct chemical attachment to the surface, leading to homogeneous monolayers with no apparent phase separation. Deprotonating the thiols by depositing the molecules under basic conditions eliminates the lateral S–H···S interactions while accentuating the dipole–dipole forces. The molecular system investigated is composed of isomeric molecules with opposite orientations of dipoles and identical surface packing, which enables the mapping of individual molecules within the mixed monolayers and enables analyses of the contributions of the relatively weak lateral interactions to the overall stability of the assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Interpreting Tafel behavior of consecutive electrochemical reactions through combined thermodynamic and steady state microkinetic approaches

Assessing the reaction pathway of multi-electron-transfer reactions is an essential yet difficult task for the rational design of electrocatalysts. In this work, we develop a heuristic approach that combines thermodynamic adsorption energetics calculated through density functional theory with microkinetic modeling using the steady state approximation to interpret the potential-dependent Tafel behavior of consecutive electrochemical reactions. In doing so, we introduce a kinetic framework for ab initio calculations that ensures self-consistent adsorption energetics based on kinetically limited adsorbate coverages. The approach is applied to experimental results on CoO x (OH) 2-x single crystal electrocatalyst particles yielding coverage dependent mechanistic information and identification of the rate-limiting step with standard rate constants for the oxygen evolution reaction on the ($11\bar{2}0$) surfaces of the β-Co(OH) 2 , β-CoOOH, and CoO 2 bulk phases. Finally, this generalizable method enables catalyst benchmarking based on determining the active species involved and associated intrinsic reaction rate constants in consecutive multi-electron-transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mechanical behavior and texture evolution of WE43 magnesium-rare earth alloy in Split-Hopkinson Pressure Bar and Taylor Impact Cylinder Testing

Mechanical behavior and texture evolution of Mg rare-earth alloy WE43 is investigated for strain-rates 10 -3 /s to upwards of 10 5 /s for the two material conditions - as-cast (AC) and T6 age hardened, rolled plate (RT6). The high strain-rate behavior is tested using both Taylor cylinder impact tests (TC) and split Hopkinson pressure bar tests (SHB) and bulk textures are obtained using neutron diffraction. Unlike the quasi-static strained material, AC and RT6 SHB retained high hardening rates throughout the test, even up to 30% true strain. Moreover, the high strain-rate data revealed that the RT6 material has a much higher strength than the AC material, but similar hardening rates despite significantly different initial texture. The flow stress near yield increased up to I 0% for RT6 and up to 30% for AC as the strain-rate increased six orders of magnitude from quasi-static rates 10 -3 /s to 10 3 /s. Neither material exhibited significant plastic anisotropy over the broad range of strain rates, despite the fact that the RT6 material had a moderately strong initial texture. In the TC tests, the geometric cross-sectional changes and texture along the cylinder from the cylindrical sample foot to head are measured and from the neutron diffraction texture analysis, upper-bound estimates of twin volume fraction are obtained as well as dislocation density from analyzing diffraction peak broadening. Recorded geometrical changes along several loading directions show that the material has deformed remarkably homogeneously under impact. Analysis of deformed textures indicates that {101¯2} extension deformation twinning occurred in the RT6 condition over the range of strain rates, with an upper bound estimate of 40% twin volume fraction for approximately 0.10-0.25 true strain. The peak texture components after the impact have their c-axes closely aligned with the impact direction. As a result, these observations are presented and rationalized in the paper.

42 ENGINEERING↗