SEARCH · Search NASA

Results for “Variational formulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Influence of extreme temperature conditions on CO 2 direct air capture using amino-acid solutions

Geological features play a pivotal role in determining the feasibility of deploying CO₂ direct air capture (DAC) technologies, primarily because they influence the availability of cost-effective energy sources, such as natural gas and geothermal energy, and also due to the potential for CO₂ sequestration. Many regions face challenges due to variable weather conditions including seasonal temperature fluctuations, high or low humidity, and sub-ambient temperatures. These extremes can reduce DAC performance or even lead to catastrophic events. Aqueous solvents considered for DAC systems are particularly vulnerable to seasonal variations in colder climates, where the solvent may underperform or freeze. It is therefore essential to investigate the CO₂ capture efficiency of aqueous solvents across a broad range of environmental temperatures, spanning sub-zero to hot conditions (>30 °C). In this study, DAC operation is examined using a high-flux solvent–air crossflow contactor under two major weather scenarios: (i) cold conditions below 0 °C and (ii) hot conditions above 30 °C. A parametric study is conducted to investigate the contactor performance regarding CO₂ removal efficiency, uptake capacity, and reaction kinetics versus temperature when the air velocity through the contactor exceeds 1 m/s. The efficacy of the contactor is systematically investigated using various anti-freeze amino-acid solvent formulations. A mass-transfer mechanistic model is developed to assess the process performance over a wide temperature range and propose scalable design guidelines. Machine learning is also employed to identify key parameters affecting the CO₂ capture efficiency. It is shown that air velocity and temperature are the primary factors influencing CO₂ uptake. Based on performance data obtained under subfreezing temperatures, a technoeconomic analysis is conducted to evaluate the feasibility of using aqueous solvents in seasonal cold regions. In conclusion, the findings of this study provide valuable insights into siting considerations for deploying solvent-based DAC, thereby contributing to the advancement of sustainable carbon removal solutions.

Air–liquid contactor↗

DOE OSTI2025

Scaling Exponents of Turbulent Static Pressure Structure Function in the Inertial Subrange

The measured variations in the turbulent static pressure structure function D pp (r) with scale $r$ in the roughness sublayer above a subarctic forest are empirically shown to exhibit exponents that are smaller than $r$ 4/3 predicted for the inertial subrange (ISR). Three hypotheses are offered to explain these deviations. The first is based on conventional intermittency correction to the averaged turbulent kinetic energy dissipation rate, the second is based on shearing introducing deviations from locally isotropic state that must be sensed by both velocity and pressure structure functions, and the third is based on large and inertial scale pressure interactions that persist at values of $r$ within the resolvable ISR. The third hypothesis is shown to yield superior results, which allows a new formulation for D pp (r) to be derived that accommodates such finite interactions.

structure function↗

DOE OSTI2026

Multi-strain analysis of Pseudomonas putida reveals the metabolic and genetic diversity of the species

Pseudomonas putida is a gram-negative bacterial species increasingly utilized in biotechnology due to its robust growth, ability to degrade aromatic compounds, solvent tolerance, and genetic tractability. In this study, we report a comprehensive multi-strain analysis of 164 P. putida strains based on the reconstruction of a pan-putida metabolic network and the formulation of strain-specific genome-scale metabolic models (GEMs). We performed whole-genome sequencing and hybrid assembly for 40 strains, contributing a ~8% increase to the available genomic data for P. putida . Furthermore, high-throughput phenotypic profiling using the Biolog phenotype microarray system for 24 strains on 190 unique carbon sources, along with 15 aromatic compounds not present on Biolog plates, yielded 4,920 unique strain-phenotype measurements. These data were leveraged to curate GEMs for 24 representative strains, including a refined model for strain KT2440, which comprised 1,480 genes and 2,191 metabolites, achieving a prediction accuracy of 91.2% in carbon utilization. Systematic comparison of genomes and GEMs revealed both conserved core pathways and significant allelic and functional divergence across strains, highlighting strain-specific variation in aromatic degradation. While pathways for protocatechuate and phenylacetate degradation were widely conserved, metabolic capabilities for compounds such as ferulate, phenol, and cresols varied markedly, suggesting adaptation to distinct ecological niches. Alleleome analysis of enzymes, such as PcaI and PcaJ, revealed distinct, functionally similar clades, indicating possible convergent evolution or horizontal gene transfer. These results provide computable resources and informative models for selecting P. putida strains with desired traits for biomanufacturing and bioremediation and offer insights into the evolution and phylogeny of the P. putida species.

aromatics utilization↗

DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

A fast multi-fidelity method with uncertainty quantification for complex data correlations: Application to vortex-induced vibrations of marine risers

Here we develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.

42 ENGINEERING↗

DOE OSTI2023

BEYONDPLANCK VIII. Efficient sidelobe convolution and corrections through spin harmonics

We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BEYONDPLANCK, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3–10, depending on the frequency bandlimits, l max and m max . The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI are in good agreement with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Surviving Under Pressure: The Role of Solvent, Crystal Size, and Morphology During Pelletization of Metal–Organic Frameworks

As metal–organic frameworks (MOFs) gain traction for applications, such as hydrogen storage, it is essential to form the as-synthesized powder materials into shaped bodies with high packing densities to maximize their volumetric performance. Mechanical compaction, which involves compressing the materials at high pressure, has been reported to yield high monolith density but often results in a significant loss in accessible porosity. Herein, we sought to systematically control (1) crystal size, (2) solvation, and (3) compacting pressure in the pelletization process to achieve high packing density without compromising the porosity that makes MOFs functional. It was determined that solvation is the most critical factor among the three factors examined. Solvation that exceeds the pore volume prevents the framework from collapsing, allowing for porosity to be maintained through pelletization. Higher pelletization pressure results in higher packing density, with extensive loss of porosity being observed at a higher pressure if the solvation is below the pore volume. Lastly, we observed that the morphology and size of the MOF particles result in variation in the highest achievable packing efficiency, but these numbers (75%) are still greater than many existing techniques used to form MOFs. Furthermore, we concluded that the application of pressure through pelletization is a suitable and widely applicable technique for forming high-density MOF-monoliths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems

Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

Calculus of variations↗

DOE OSTI2024

Quantum Davidson algorithm for excited states

Abstract Excited state properties play a pivotal role in various chemical and physical phenomena, such as charge separation and light emission. However, the primary focus of most existing quantum algorithms has been the ground state, as seen in quantum phase estimation and the variational quantum eigensolver (VQE). Although VQE-type methods have been extended to explore excited states, these methods grapple with optimization challenges. In contrast, the quantum Krylov subspace (QKS) method has been introduced to address both ground and excited states, positioning itself as a cost-effective alternative to quantum phase estimation. However, conventional QKS methodologies depend on a pre-generated subspace through real or imaginary-time evolutions. This subspace is inherently expansive and can be plagued with issues like slow convergence or numerical instabilities, often leading to relatively deep circuits. Our research presents an economic QKS algorithm, which we term the quantum Davidson (QDavidson) algorithm. This innovation hinges on the iterative expansion of the Krylov subspace and the incorporation of a pre-conditioner within the Davidson framework. By using the residues of eigenstates to expand the Krylov subspace, we manage to formulate a compact subspace that aligns closely with the exact solutions. This iterative subspace expansion paves the way for a more rapid convergence in comparison to other QKS techniques, such as the quantum Lanczos. Using quantum simulators, we employ the novel QDavidson algorithm to delve into the excited state properties of various systems, spanning from the Heisenberg spin model to real molecules. Compared to the existing QKS methods, the QDavidson algorithm not only converges swiftly but also demands a significantly shallower circuit. This efficiency establishes the QDavidson method as a pragmatic tool for elucidating both ground and excited state properties on quantum computing platforms.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗

DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗

DOE OSTI2020

Probing Conformational Evolution and Associated Dynamics of Mg(N(SO 2 CF 3 ) 2 ) 2 ·Dimethoxyethane Adduct Using Solid-State 19 F and 1 H NMR

Bis(trifluoromethanesulfonimide) or TFSI is widely used as a counter anion in electrolyte design due to its structural flexibility and chemical stability. Here we studied the conformational variations and associated dynamics of TFSI in adduct of Mg(TFSI) 2 with dimethoxyethane (DME), a solvate crystalline material using solid-state 1 H and 19 F NMR. TFSI molecular motion in this solvate structure falls within the timescale of the 19 F NMR experiment, yielding spectroscopic signatures for unique TFSI conformers under the coordination environment of Mg 2+ cation. Within the temperature range of -5 to 82 °C, we observe nine distinct TFSI sites in both crystalline and disordered regions using 19 F NMR, reflecting complexity of structural and dynamics of TFS1 anions within solvate structure. The four distinguishable sites in the disordered region for the two CF 3 groups of the same TFSI molecule are identified using chemical shift analysis. The exchange rate constants from site to site are calculated through variable temperature 19 F NMR and two-dimensional (2D) exchange spectroscopy (EXSY) experiments, along with respective activation enthalpies using Eyring's formulation. The flip rate of CF 3 around the S-C bond is estimated as ~ 15 s -1 at 8 °C with ΔH ≠ ~ 22 kJ/mol, but the rotation of the entire TFSI is 4.8 s -1 at 8 °C with a significantly greater ΔH ≠ = 98 ± 10 kJ/mol. Furthermore, the slow conversion of trans to cis conformers at a lower temperature (T ≤ 1 °C) in the crystalline region is monitored, with a conversion rate of ~ 2 x 10 -5 s -1 at -5 °C. Density functional theory (DFT)-based calculations were performed to support further the assignment of experimental chemical shifts, and the activation energy E a = 21.1 kJ/mol obtained for the cis to trans transition is consistent with experimental values. The combined set of 19 F and 1 H under both one-dimensional (1D) and 2D NMR methods demonstrated here can be further used for examining electrode-electrolyte interfaces to probe the motions of various constituents that can enable detailed studies of interfacial processes and dynamics. Ultimately, such studies will aid in the design and discovery of interfacial constructs in which directed defect chemistry, chemical moiety distribution, and nanostructure are employed to drive efficient charge transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Topologically derived dislocation theory for twist and stretch moiré superlattices in bilayer graphene

In this work, we develop a continuum dislocation description of twist and stretch moiré superlattices in two-dimensional material bilayers. The continuum formulation is based on the topological constraints introduced by the periodic dislocation network associated with the moiré structure. The approach is based on solving analytically for the structural distortion and displacement fields that satisfy the topological constraints and that minimize the total energy. The total energy is described by both the strain energy of each individual distorted layer and a Peierls-Nabarro-like interfacial contribution arising from stacking disregistry. The dislocation core emerges naturally within the formalism as a result of the competition between the two contributions. The approach presented here captures the structure and energetics of twist and stretch moiré superlattices of dislocations with arbitrary direction and character, without assuming an analytical solution a priori and while accounting naturally for dislocation-dislocation image interactions. In comparisons to atomistic simulations using classical potentials, the maximum structure deviation is 6%, while the maximum line energy deviation is 0.019 eV/Å. Several applications of our model are shown, including predicting the variation of structure with twist angle and describing dislocation line tension and junction energies.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗

DOE OSTI2021

Dissolution of periodic arrays of grains: Upscaling of pore-scale simulations with fast reactions

Upscaling methods are frequently used to derive transport equations at the macroscopic scale from more fundamental equations formulated at the pore scale. These methods typically give a suitable structure for the macroscopic equations and can also provide explicit expressions for the constitutive parameters, such as permeability and dispersion coefficients. Introducing chemical reactions complicates upscaling in at least two important ways. First, the interplay between chemical reactions and transport processes introduces a new length scale, which can be much smaller than the convective or dispersive length scales. A small reactive length scale breaks one of the key assumptions in upscaling; that there is a significant separation in length between the pore-scale and macro-scale processes. The second complication is that if reactions take place at mineral surfaces (dissolution or precipitation) then the pore space itself is evolving in time. In this paper we suggest ways in which these difficulties can be approached, based on analysis of pore-scale simulation data. First, we noticed that the concentration field in successive unit cells has an almost identical spatial variation, with a single scaling factor for each unit cell that is proportional to the incoming reactant flux. Using pore-scale simulations to determine the mass transfer coefficient in a few unit cells, we can calculate the concentration field in the whole domain, even when dissolution is entirely transport limited. Second, we have noticed a time-dependent mapping of the grain shapes from different unit cells. From these observations, we can deduce constitutive relations where the only time-varying parameter is the porosity. Here, we show that a model based on these ideas can quantitatively account for the pore-scale simulation data.

58 GEOSCIENCES↗

DOE OSTI2020

Quantum maximin surfaces

We formulate a quantum generalization of maximin surfaces and show that a quantum maximin surface is identical to the minimal quantum extremal surface, introduced in the EW prescription. We discuss various subtleties and complications associated to a maximinimization of the bulk von Neumann entropy due to corners and unboundedness and present arguments that nonetheless a maximinimization of the UV-finite generalized entropy should be well-defined. We give the first general proof that the EW prescription satisfies entanglement wedge nesting and the strong subadditivity inequality. In addition, we apply the quantum maximin technology to prove that recently proposed generalizations of the EW prescription to nonholographic subsystems (including the so-called “quantum extremal islands”) also satisfy entanglement wedge nesting and strong subadditivity. Our results hold in the regime where backreaction of bulk quantum fields can be treated perturbatively in G N h, but we emphasize that they are valid even when gradients of the bulk entropy are of the same order as variations in the area, a regime recently investigated in new models of black hole evaporation in AdS/CFT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Antioxidant Activity, Microbiological Quality, and Acceptability of Spontaneously Fermented Shrimp Sausage (Litopenaeus vannamei)

Fermented shrimp sausages are prepared spontaneously with Litopenaeus vannamei as raw material. Shrimp is one of the marine sources with rich chitosan as bioactive compounds, antioxidants, vitamin E, and probiotic lactic acid bacteria (LAB) produced by fermentation processes. This study aimed to analyze the variations in antioxidant activity, vitamin E content, total LAB, total pathogenic bacteria, pH, and acceptability of shrimp sausage produced at different fermentation times. A completely randomized experimental design study was performed using four levels of exposure time, including control (0 days), 1, 2, and 3 days. The treatment was conducted spontaneously with 1.2% salt concentration, a drying temperature of 50°C for 3 h, and fermentation at 35°C. The evaluated parameters included the antioxidant activity measured with the radical 2,2-azino-bis (3-ethylbenzothiazoline-6-sulphonic acid) method, vitamin E levels by high-performance liquid chromatography, total LAB with total plate count, and Escherichia coli bacteria by the most probable number method. Salmonella sp. and Staphylococcus aureus were estimated by the identification method. A pH meter was used to assess acidity, and hedonic organoleptic testing was performed for taste, aroma, color, and texture. The results show significant differences in antioxidant activity, vitamin E, and shrimp sausage acceptability at varied fermentation times p ≤ 0.05 . However, the best formulation was obtained with the first-day fermented shrimp sausage, as shown by its bioactive content and the level of acceptability.

Afifah, Diana N.↗