Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational approach”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory

Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.

Richardson, G. A.↗
DOE OSTI2024

Volatility Basis Set Distributions and Viscosity of Organic Aerosol Mixtures: Insights from Chemical Characterization Using Temperature-Programmed Desorption–Direct Analysis in Real-Time High-Resolution Mass Spectrometry

Quantitative assessment of gas-particle partitioning of individual components within complex atmospheric organic aerosols (OA) mixtures is critical for predicting and comprehending the formation and evolution of OA particles in the atmosphere. This investigation leverages previously documented data obtained through a temperature programmed desorption - direct analysis in real time – high resolution mass spectrometry (TPD-DART-HRMS) platform. This methodology facilitates the bottom-up construction of volatility basis set (VBS) distributions for constituents found in three biogenic secondary organic aerosol (SOA) mixtures produced through the ozonolysis of a-pinene, limonene, and ocimene. The apparent enthalpies (ΔH*, kJ mol -1 ) and saturated vapor mass concentrations (C T *, µg∙m -3 ) of individual SOA components, determined as a function of temperature (T, K), facilitated an assessment of changes in VBS distributions and gas-particle partitioning with respect to T and atmospheric total organic mass loadings (tOM, µg∙m -3 ). Further, the VBS distributions reveal distinct differences in volatilities among monomers, dimers, and trimers, enabling their categorization into separate volatility bins. At the ambient temperature of T = 298 K, only monomers efficiently partition between gas and particle phases across a broad range of atmospherically relevant total organic mass loadings (tOM) values of 1–100 µg∙m -3 . Partitioning of dimers and trimers becomes notable only at T > 360 K and T > 420 K, respectively. The viscosity of SOA mixtures is assessed using a bottom-up calculation approach, incorporating the input of elemental formulas, ΔH*, C T *, and particle-phase mass fractions of the SOA components. Through this approach, we are able to accurately estimate the variations in SOA viscosity that result from the evaporation of its components. These variations are, in turn, influenced by atmospherically relevant changes in tOM and T. Comparison of the calculated SOA viscosity and diffusivity values with literature reported experimental results shows close agreement, thereby validating the employed calculation approach. These findings underscore the significant potential for TPD-DART-HRMS measurements in enabling the untargeted analysis of organic molecules within OA mixtures. This approach facilitates quantitative assessment of their gas-particle partitioning and allows for the estimation of their viscosity and condensed-phase diffusion, thereby contributing valuable insights to atmospheric models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Exploring Quantum State Preparation Using Tensor Networks and Sparse Wavefunction Simulations

The variational quantum eigenvalue solver is a powerful hybrid quantum-classical approach that has been suggested as a candidate method to run on near-term quantum hardware for computing ground state electronic energies of molecular systems. However, even for small molecules, the number of variational parameters and qubits required to minimize the electronic energy is beyond the reach of current quantum computers except for small basis sets. We explore a new paradigm for state preparation where we test how much of the optimization can be approximately prepared with classical computers to reduce the number of optimization steps performed using a quantum device. By adapting a recent algorithm for the factorized form of the UCC ansatz, we can study molecular electronic structure problems with up to 64 qubits. In addition, we also test a related approach of using tensor networks to optimize quantum circuits in order to benchmark various lattice models. We present results using these approaches and discuss strategies for incorporating these ideas into variational algorithms involving near-term quantum computers. Our results help demonstrate the strength of the UCC ansatz and address pressing questions about optimal initial parameterizations and circuit construction.

quantum computing↗
NASA NTRS1991

Efficient real gas upwinded Navier-Stokes computations of high speed flows

An efficient method to account for the chemically frozen thermodynamic and transport properties of air in three-dimensional Navier-Stokes calculations has been demonstrated. This approach uses an explicitly specified equation of state (EOS) so that the fluid pressure, temperature, and transport properties are directly related to the flow variables. The method is efficient since no subiterations are required to deduce the pressure and temperature from the flux variables and is modular by allowing different equations of state to be easily supplied to the code. The flexibility of the EOS approach is shown by its implementation into a high-order total variation diminishing upwinding scheme as well as a standard central-differencing scheme. The EOS approach is then demonstrated by computing the hypersonic flow through the corner region of two mutually perpendicular flat plates using both the upwind and central-differencing schemes.

Coirier, William J.↗
NASA NTRS1975

On the prediction of critical flutter conditions from subcritical response data and some related wind-tunnel experience

Methods of interpreting response measurements which could be amenable to flight flutter testing procedures were studied analytically and in the wind tunnel. One suggested scheme, which requires evaluation, is an iterative technique in which derivatives obtained from subcritical response data are used to indicate the approach to flutter. A simplification of this procedure was considered by examining the manner in which a single characteristic of the subcritical response behaves in relation to variations of the density or dynamic pressure in the approach to flutter. The use of this single parameter scheme was examined for random excitation as well as for sinusoidal forcing. The feasibility of the method is illustrated by several examples and the relative merits of random and sinusoidal excitation are discussed.

Houbolt, J. C.↗
NASA NTRS1984

A variational method for the calculation of dynamic polarizabilities and two-photon transition moments - The dressed molecule approach

The solution of the time-dependent Schroedinger equation of the molecule + radiation field system is analyzed. A quantized radiation field is used. The relationship between the oscillatory wave function and the dressed state wave function is established, the oscillatory wave function being the solution if the radiation field is turned on at t = 0 and the dressed state being the stationary solution satisfying the boundary condition only if the field is present at t = negative infinity. In general, the oscillatory wave function can be expressed using a complete set of dressed states. However, in the presence of a nonresonant radiation field, the system is well represented by a single dressed state. Molecular properties such as dynamic polarizabilities and two-photon transition moments can be deduced from the dressed state wave function instead of the oscillatory wave function as in previous methods. Because of its stationary character, the dressed state is more amenable to approximate calculations. A CI method is developed for this purpose. The CI matrix is simple in form and can be readily constructed using existing computer codes. The present method can also be adapted to calculate other optical properties.

Huo, W. M.↗
NASA NTRS2015

A Complementary Note to 'A Lag-1 Smoother Approach to System-Error Estimation': The Intrinsic Limitations of Residual Diagnostics

Recently, this author studied an approach to the estimation of system error based on combining observation residuals derived from a sequential filter and fixed lag-1 smoother. While extending the methodology to a variational formulation, experimenting with simple models and making sure consistency was found between the sequential and variational formulations, the limitations of the residual-based approach came clearly to the surface. This note uses the sequential assimilation application to simple nonlinear dynamics to highlight the issue. Only when some of the underlying error statistics are assumed known is it possible to estimate the unknown component. In general, when considerable uncertainties exist in the underlying statistics as a whole, attempts to obtain separate estimates of the various error covariances are bound to lead to misrepresentation of errors. The conclusions are particularly relevant to present-day attempts to estimate observation-error correlations from observation residual statistics. A brief illustration of the issue is also provided by comparing estimates of error correlations derived from a quasi-operational assimilation system and a corresponding Observing System Simulation Experiments framework.

Fixed-lag smoothing↗
NASA NTRS1994

(abstract) Optimal Low Thrust Trajectories Using Differential Inclusion Concepts

Low thrust propulsion systems typically have their greatest benefit for high energy missions or missions with large post-launch maneuver requirements. Missions which have been examined include main belt asteroid rendezvous, comet rendezvous, outer planet and Mercury orbiters, Pluto flyby, and solar probe missions. Low thrust mission design software used to determine these trajectories is based on two distinct formulations of the optimal control problem: the indirect and direct methods. The traditional approach (indirect) is to use the calculus of variations to obtain first order necessary conditions on the states and costates. In contrast, direct methods are conceptually different in that no explicit integration takes place. A direct method based on differential inclusion concepts has been developed and successfully used to compute low thrust trajectories. This new approach removes explicit control dependence from the problem thereby reducing the dimension of the parameter space for the nonlinear programming problem. Also when compared to other direct methods, fewer nonlinear constraints are required to represent the dynamics of the problem.

low↗
NASA NTRSDate not supplied

Variational Coupled Loads Analysis using the Hybrid Parametric Variation Method

Time-domain coupled loads analysis (CLA)is used to determine the response of a launch vehicle and payload system to transient forces, such as liftoff, engine ignitions and shutdowns, jettison events, and atmospheric flight loads, such as buffet. CLA, using Hurty/Craig-Bampton (HCB)component models, is the accepted method for the establishment of design-level loads for launch systems. However, uncertainty in the component models flows into uncertainty in predicted system results. Uncertainty in the structural responses during launch is a significant concern because small variations in launch vehicle and payload mode shapes and their interactions can result in significant variations in system loads. Uncertainty quantification (UQ)is used to determine statistical bounds on prediction accuracy based on model uncertainty. In this paper uncertainty is treated at the HCB component-model level. In an effort to account for model uncertainties and statistically bound their effect on CLA predictions, this work combines CLA with UQ in a process termed variational coupled loads analysis (VCLA). The modeling of uncertainty using a parametric approach, in which input parameters are represented by random variables, is common, but its major drawback is the resulting uncertainty is limited to the form of the nominal model. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. Model-form uncertainty can be represented using a nonparametric approach based on random matrix theory (RMT). In this work, UQ is performed using the hybrid parametric variation (HPV)method, which combines parametric with nonparametric uncertainty at the HCB component model level. The HPV method requires the selection of dispersion values for the HCB fixed-interface (FI)eigenvalues, and the HCB mass and stiffness matrices. The dispersions are based upon component test-analysis modal correlation results. During VCLA, random component models are assembled into an ensemble of random systems using a Monte Carlo (MC)approach. CLA is applied to each of the ensemble members to produce an ensemble of system-level responses for statistical analysis. The proposed methodology is demonstrated through its application to a buffet loads analysis of NASA’s Space Launch System (SLS)during the transonic regime fifty seconds after liftoff. Core stage (CS)section shears and moments are recovered, and statistics are computed.

Uncertainty Quantification↗
DOE OSTI2025

Frontal Slice Approaches for Tensor Linear Systems

Inspired by the row and column action methods for solving large-scale linear systems, in this work, we explore the use of frontal slices for solving tensor linear systems. In particular, this paper presents a novel approach for using frontal slices of a tensor $\mathcal{A}$ to solve tensor linear systems $\mathcal{A} ∗\mathcal{X} = \mathcal{B}$ where ∗ denotes the $t$-product. In addition, we consider variations of this method, including cyclic, block, and randomized approaches, each designed to optimize performance in different operational contexts. Our primary contribution lies in the development and convergence analysis of these methods. Experimental results on synthetically generated and real-world data, including applications such as image and video deblurring, demonstrate the efficacy of our proposed approaches and validate our theoretical findings.

Luo, Hengrui↗
DOE OSTI2025

Meta Biome: a multiscale model integrating agent-based and metabolic networks to reveal spatial regulation in gut mucosal microbial communities

ABSTRACT Mucosal microbial communities (MMCs) are complex ecosystems near the mucosal layers of the gut essential for maintaining health and modulating disease states. Despite advances in high-throughput omics technologies, current methodologies struggle to capture the dynamic metabolic interactions and spatiotemporal variations within MMCs. In this work, we presentMetaBiome, a multiscale model integrating agent-based modeling (ABM), finite volume methods, and constraint-based models to explore the metabolic interactions within these communities. Integrating ABM allows for the detailed representation of individual microbial agents each governed by rules that dictate cell growth, division, and interactions with their surroundings. Through a layered approach—encompassing microenvironmental conditions, agent information, and metabolic pathways—we simulated different communities to showcase the potential of the model. Using ourin-silicoplatform, we explored the dynamics and spatiotemporal patterns of MMCs in the proximal small intestine and the cecum, simulating the physiological conditions of the two gut regions. Our findings revealed how specific microbes adapt their metabolic processes based on substrate availability and local environmental conditions, shedding light on spatial metabolite regulation and informing targeted therapies for localized gut diseases.MetaBiome provides a detailed representation of microbial agents and their interactions, surpassing the limitations of traditional grid-based systems. This work marks a significant advancement in microbial ecology, as it offers new insights into predicting and analyzing microbial communities. IMPORTANCE Our study presents a novel multiscale model that combines agent-based modeling, finite volume methods, and genome-scale metabolic models to simulate the complex dynamics of mucosal microbial communities in the gut. This integrated approach allows us to capture spatial and temporal variations in microbial interactions and metabolism that are difficult to study experimentally. Key findings from our model include the following: (i) prediction of metabolic cross-feeding and spatial organization in multi-species communities, (ii) insights into how oxygen gradients and nutrient availability shape community composition in different gut regions, and (iii) identification of spatiallyregulated metabolic pathways and enzymes inE. coli. We believe this work represents a significant advance in computational modeling of microbial communities and provides new insights into the spatial regulation of gut microbiome metabolism. The multiscale modeling approach we have developed could be broadly applicable for studying other complex microbial ecosystems.

Microbiology↗
NASA NTRS2002

Computational Relativistic Astrophysics Using the Flow Field-Dependent Variation Theory

We present our method for solving general relativistic nonideal hydrodynamics. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks which may lead to the study of gamma-ray bursts. Nonideal flows are present where radiation, magnetic forces, viscosities, and turbulence play an important role. Our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flow field-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for computational relativistic astrophysics (CRA) are demonstrated.

Richardson, G. A.↗
DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗
DOE OSTI2026

Flow synthesis of iron oxide nanoparticles: using multiple precursor additions to improve size control

Controlling the size of iron oxide nanoparticles while maintaining uniformity in flow-based synthesis systems has been a great challenge in nanoparticle synthesis. Using an extended LaMer mechanism, we improve both the size and shape uniformity as compared to a conventional flow synthesis. The key to this approach is the injection of additional Fe precursor during the flow reaction, providing extra precursor during the growth process leading to a size focusing step. This approach has also allowed for the systematic variation of the nanoparticle size produced in the flow reaction. This work represent a new synthetic concept for improving the size uniformity of iron oxide nanoparticles in flow-based synthesis and has potential for application across a wide range of nanoparticle systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamics of argon metastables in Ar–CH 4 radio frequency capacitively-coupled plasma: real-time monitoring with neural network-augmented broadband optical emission spectroscopy

In moderate-pressure radio frequency (RF) capacitively coupled plasmas generated in argon–methane mixtures, the density of argon metastable atoms (Ar 1 s 5 ) exhibits a non-monotonic dependence on methane (CH 4 ) concentration. Laser-induced fluorescence (LIF) was used to measure and compare local Ar 1 s 5 densities in Ar and Ar–CH 4 plasmas at 2.6 Torr and RF powers of 17–117 W. The addition of 1% CH 4 increases the metastable density, and 2% CH 4 triggers a strong depletion by an order of magnitude, compared to 1% CH 4 case. This non-monotonic behavior demonstrates the sensitivity of metastable populations to small gas admixtures, which is critical for processes where metastables drive precursor dissociation. For real-time monitoring of metastable population, broadband optical emission spectroscopy (OES) is augmented with a feedforward neural network (NN) to predict Ar s1 5 densities from spectral features. When trained on LIF data, the NN replicates the absolute densities and the dynamic trends of Ar 1 s 5 density variation. The NN-augmented broadband OES approach can be used as a simple and cost-effective tool for tracking Ar metastables in Ar-rich plasmas, facilitating industrial-scale optimization.

Yatom, Shurik [Princeton Plasma Physics Laboratory↗
DOE OSTI2025

Analysis of Polarized Dust Emission Using Data from the First Flight of SPIDER

Using data from the first flight of Spider and from the Planck High Frequency Instrument, we probe the properties of polarized emission from interstellar dust in the Spider observing region. Component-separation algorithms operating in both the spatial and harmonic domains are applied to probe their consistency and to quantify modeling errors associated with their assumptions. Analyses of diffuse Galactic dust emission spanning the full Spider region demonstrate (i) a spectral energy distribution that is broadly consistent with a modified-blackbody (MBB) model with a spectral index of β d = 1.45 ± 0.05 (1.47 ± 0.06) for E (B)-mode polarization, slightly lower than that reported by Planck for the full sky; (ii) an angular power spectrum broadly consistent with a power law; and (iii) no significant detection of line-of-sight polarization decorrelation. Tests of several modeling uncertainties find only a modest impact (~10% in σ r ) on Spider's sensitivity to the cosmological tensor-to-scalar ratio. The size of the Spider region further allows for a statistically meaningful analysis of the variation in foreground properties within it. Assuming a fixed dust temperature T d = 19.6 K, an analysis of two independent subregions of that field results in inferred values of β d = 1.52 ± 0.06 and β d = 1.09 ± 0.09, which are inconsistent at the 3.9σ level. Furthermore, a joint analysis of Spider and Planck 217 and 353 GHz data within one subregion is inconsistent with a simple MBB at more than 3σ, assuming a common morphology of polarized dust emission over the full range of frequencies. This evidence of variation may inform the component-separation approaches of future cosmic microwave background polarization experiments.

79 ASTRONOMY AND ASTROPHYSICS↗