Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Puzzling Variation of Gamma Rays from the Sun over the Solar Cycle Revealed with Fermi-LAT

The steady-state gamma-ray emission from the Sun is thought to consist of two emission components due to interactions with Galactic cosmic rays: (1) a hadronic disk component, and (2) a leptonic extended component peaking at the solar edge and extending into the heliosphere. The flux of these components is expected to vary with the 11 yr solar cycle, being highest during solar minimum and lowest during solar maximum, as it varies with the cosmic-ray flux. No study has yet analyzed the flux variation of each component over solar cycles. In this work, we measure the temporal variations of the flux of each component over 15 yr of Fermi Large Area Telescope observations and compare them with the sunspot number and Galactic cosmic-ray flux from AMS-02 near Earth. We find that the flux variation of the disk anticorrelates with the sunspot number and correlates with cosmic-ray protons, as expected, confirming its emission mechanism. In contrast, the extended component exhibits a more complex variation: despite an initial anticorrelation with the sunspot number, we find neither anticorrelation with the sunspot number nor correlation with cosmic-ray electrons over the full 15 yr period. This most likely suggests that cosmic-ray transport and modulation in the inner heliosphere are unexpectedly complex and may differ for electrons and protons or, alternatively, that there is an additional, unknown component of gamma rays or cosmic rays. These findings impact space weather research and emphasize the need for close monitoring of Cycle 25 and the ongoing polarity reversal.

Acharyya, A. [Univ. of Southern Denmark, Odense (D↗
DOE OSTI2021

Strong temporal variation in treefall and branchfall rates in a tropical forest is related to extreme rainfall: results from 5 years of monthly drone data for a 50 ha plot

Abstract. A mechanistic understanding of how tropical-tree mortality responds to climate variation is urgently needed to predict how tropical-forest carbon pools will respond to anthropogenic global change, which is altering the frequency and intensity of storms, droughts, and other climate extremes in tropical forests. We used 5 years of approximately monthly drone-acquired RGB (red–green–blue) imagery for 50 ha of mature tropical forest on Barro Colorado Island, Panama, to quantify spatial structure; temporal variation; and climate correlates of canopy disturbances, i.e., sudden and major drops in canopy height due to treefalls, branchfalls, or the collapse of standing dead trees. Canopy disturbance rates varied strongly over time and were higher in the wet season, even though wind speeds were lower in the wet season. The strongest correlate of monthly variation in canopy disturbance rates was the frequency of extreme rainfall events. The size distribution of canopy disturbances was best fit by a Weibull function and was close to a power function for sizes above 25 m2. Treefalls accounted for 74 % of the total area and 52 % of the total number of canopy disturbances in treefalls and branchfalls combined. We hypothesize that extremely high rainfall is a good predictor because it is an indicator of storms having high wind speeds, as well as saturated soils that increase uprooting risk. These results demonstrate the utility of repeat drone-acquired data for quantifying forest canopy disturbance rates at fine temporal and spatial resolutions over large areas, thereby enabling robust tests of how temporal variation in disturbance relates to climate drivers. Further insights could be gained by integrating these canopy observations with high-frequency measurements of wind speed and soil moisture in mechanistic models to better evaluate proximate drivers and with focal tree observations to quantify the links to tree mortality and woody turnover.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

An Introduction to Variational Calculus

Variational Calculus is an advanced topic in mathematics. This document exists to compile the basics of Variational Calculus in a readable and digestible manner for application by engineers. This document is a living document and will be updated as time moves forward. It is the hope that by reading this document, the basics of Variational Calculus can be readily applied to any system/function of interest. This does not serve as a replacement for textbooks or other learned sources but should act as a companion piece to better serve the reader in times of confusion and distress.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗
DOE OSTI2024

Investigation of inter-subject variation in ultrafine particle deposition across human nasal airways: A study involving children, adults, and the elderly

Occupational and environmental exposure to toxic nanoparticles, driven by the rapid expansion of nanotechnology, raises significant respiratory health concern. Numerous studies have explored airflow and particle dynamics in adult nasal airways, but understanding the impact of age-related anatomical changes in children and the elderly remains limited. This study systematically investigates age-related anatomical variations and associated influence on nasal airflow dynamics and ultrafine particle deposition characteristics. Using Computational Fluid-Particle Dynamics (CFPD) method, simulation was conducted under diverse inhalation conditions spanning a wide age range, including: two children (5 years old), two young adults (in their twenties), and two elderly (over 77 years old). Our results reveal distinctive variations across age groups in anatomical dimensions, which affect distribution of wall shear stress where the elderly and children display unique patterns distinct from the young adults. While total deposition efficiency differs significantly between children and adults, filtration efficiency in the subregion with most deposition, main respiratory, remains consistent. However, inter-subject differences are observed in the vestibular and olfactory regions,emphasizing nuanced impact of age-related anatomical variations. Overall and subregional empirical equations for deposition efficiency were developed by incorporating the combined diffusion parameter, Sc a Δ b , corroborating the use of geometrical characteristic parameters for each specific subject in predicting nasal deposition efficiency across age groups. Our findings contribute to predictive nanoparticle exposure analysis in nasal airways across different age groups, thereby enhancing respiratory healthcare for individuals across the life span.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Variational quantum state eigensolver

Extracting eigenvalues and eigenvectors of exponentially large matrices will be an important application of near-term quantum computers. The variational quantum eigensolver (VQE) treats the case when the matrix is a Hamiltonian. Here, we address the case when the matrix is a density matrix ρ. We introduce the variational quantum state eigensolver (VQSE), which is analogous to VQE in that it variationally learns the largest eigenvalues of ρ as well as a gate sequence V that prepares the corresponding eigenvectors. VQSE exploits the connection between diagonalization and majorization to define a cost function C=Tr(ρ~H) where H is a non-degenerate Hamiltonian. Due to Schur-concavity, C is minimized when ρ~=VρV† is diagonal in the eigenbasis of H. VQSE only requires a single copy of ρ (only n qubits) per iteration of the VQSE algorithm, making it amenable for near-term implementation. We heuristically demonstrate two applications of VQSE: (1) Principal component analysis, and (2) Error mitigation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Numerical integration of stochastic contact Hamiltonian systems via stochastic Herglotz variational principle

Within this work, we establish a stochastic contact variational integrator and its discrete version via stochastic Herglotz variational principle for stochastic contact Hamiltonian systems. A general structure-preserving stochastic contact method is provided to seek the stochastic contact variational integrators. Numerical experiments are performed to verify the validity of this approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Characterizing the gene–environment interaction underlying natural morphological variation in Neurospora crassa conidiophores using high-throughput phenomics and transcriptomics

Abstract Neurospora crassa propagates through dissemination of conidia, which develop through specialized structures called conidiophores. Recent work has identified striking variation in conidiophore morphology, using a wild population collection from Louisiana, United States of America to classify 3 distinct phenotypes: Wild-Type, Wrap, and Bulky. Little is known about the impact of these phenotypes on sporulation or germination later in the N. crassa life cycle, or about the genetic variation that underlies them. In this study, we show that conidiophore morphology likely affects colonization capacity of wild N. crassa isolates through both sporulation distance and germination on different carbon sources. We generated and crossed homokaryotic strains belonging to each phenotypic group to more robustly fit a model for and estimate heritability of the complex trait, conidiophore architecture. Our fitted model suggests at least 3 genes and 2 epistatic interactions contribute to conidiophore phenotype, which has an estimated heritability of 0.47. To uncover genes contributing to these phenotypes, we performed RNA-sequencing on mycelia and conidiophores of strains representing each of the 3 phenotypes. Our results show that the Bulky strain had a distinct transcriptional profile from that of Wild-Type and Wrap, exhibiting differential expression patterns in clock-controlled genes (ccgs), the conidiation-specific gene con-6, and genes implicated in metabolism and communication. Combined, these results present novel ecological impacts of and differential gene expression underlying natural conidiophore morphological variation, a complex trait that has not yet been thoroughly explored.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Variational quantum and neural quantum states algorithms for the linear complementarity problem

Variational quantum algorithms (VQAs) are promising hybrid quantum-classical methods designed to leverage the computational advantages of quantum computing while mitigating the limitations of current noisy intermediate-scale quantum (NISQ) hardware. Although VQAs have been demonstrated as proofs of concept, their practical utility in solving real-world problems—and whether quantum-inspired classical algorithms can match their performance—remains an open question. We present a novel application of the variational quantum linear solver (VQLS) and its classical neural quantum states-based counterpart, the variational neural linear solver (VNLS), as key components within a minimum map Newton solver for a complementarity-based rigid-body contact model. We demonstrate using the VNLS that our solver accurately simulates the dynamics of rigid spherical bodies during collision events. These results suggest that quantum and quantum-inspired linear algebra algorithms can serve as viable alternatives to standard linear algebra solvers for modelling certain physical systems.

neural quantum states↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Variational quantum algorithm for estimating the quantum Fisher information

The quantum Fisher information (QFI) quantifies the ultimate precision of estimating a parameter from a quantum state and can be regarded as a reliability measure of a quantum system as a quantum sensor. However, estimation of the QFI for a mixed state is in general a computationally demanding task. In this paper we present a variational quantum algorithm called variational quantum Fisher information estimation (VQFIE) to address this task. By estimating lower and upper bounds on the QFI, based on bounding the fidelity, VQFIE outputs a range in which the actual QFI lies. This result can then be used to variationally prepare the state that maximizes the QFI, for the application of quantum sensing. In contrast to previous approaches, VQFIE does not require knowledge of the explicit form of the sensor dynamics. We simulate the algorithm for a magnetometry setup and demonstrate the tightening of our bounds as the state purity increases. For this example, we compare our bounds with literature bounds and show that our bounds are tighter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Ensemble approximate control variate estimators: Applications to multi-fidelity importance sampling.

The recent growth in multifidelity uncertainty quantification has given rise to a large set of variance reduction techniques that leverage information from model ensembles to provide variance reduction for estimates of the statistics of a high-fidelity model. In this paper we provide two contributions: (1) we utilize an ensemble estimator to account for uncertainties in the optimal weights of approximate control variate (ACV) approaches and derive lower bounds on the number of samples required to guarantee variance reduction; and (2) we extend an existing multifidelity importance sampling (MFIS) scheme to leverage control variates. Our approach directly addresses a limitation of many multifidelity sampling strategies that require the usage of pilot samples to estimate covariances. As such we make significant progress towards both increasing the practicality of approximate control variates—for instance, by accounting for the effect of pilot samples—and using multifidelity approaches more effectively for estimating low-probability events. The numerical results indicate our hybrid MFIS-ACV estimator achieves up to 50% improvement in variance reduction over the existing state-of-the-art MFIS estimator, which had already shown an outstanding convergence rate compared to the Monte Carlo method, on several problems of computational mechanics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Covariance-Free Bifidelity Control Variates Importance Sampling for Rare Event Reliability Analysis

Multifidelity modeling has been steadily gaining attention as a tool to address the problem of exorbitant model evaluation costs that makes the estimation of failure probabilities a significant computational challenge for complex real-world problems, particularly when failure is a rare event. To implement multifidelity modeling, estimators that efficiently combine information from multiple models/sources are necessary. In past works, the variance reduction techniques of control variates (CV) and importance sampling (IS) have been leveraged for this task. In this paper, we present the CVIS framework—a creative take on a coupled CV and IS estimator for bifidelity reliability analysis. The framework addresses some of the practical challenges of the CV method by using an estimator for the control variate mean and sidestepping the need to estimate the covariance between the original estimator and the control variate through a clever choice for the tuning constant. Furthermore, the task of selecting an efficient IS distribution is also considered, with a view towards maximally leveraging the bifidelity structure and maintaining expressivity. Additionally, a diagnostic is provided that indicates both the efficiency of the algorithm as well as the relative predictive quality of the models utilized. Finally, the behavior and performance of the framework is explored through analytical and numerical examples.

Markov chain Monte Carlo↗
DOE OSTI2023

Soil fertility and drought interact to determine large variations in wood production for a hyperdominant Amazonian tree species

The productivity of the Amazon Rainforest is related to climate and soil fertility. However, the degrees to which these interactions influence multiannual to decadal variations in tree diameter growth are still poorly explored. To fill this gap, we used radiocarbon measurements to evaluate the variation in tree growth rates over the past decades in an important hyperdominant species, Eschweilera coriacea (Lecythidaceae), from six sites in the Brazilian Amazon that span a range of soil properties and climate. Using linear mixed-effects models, we show that temporal variations in mean annual diameter increment evaluated over a specific time period reflect interactions between soil fertility and the drought index (SPEI-Standardized Precipitation and Evapotranspiration Index). Our results indicate that the growth response of trees to drought is strongly dependent on soil conditions, a facet of forest productivity that is still underexplored, and which has great potential for improving predictions of future tropical tree growth in the face of projected climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Little vertical and circumferential variations in stem xylem water Δ 2 H and Δ 18 O in three tree species

Vertical and circumferential patterns of xylem water isotope signatures in large trees are rarely characterized but may influence interpretation of water source investigations and soil/xylem water isotope offsets. Furthermore, to examine vertical and circumferential variation in xylem water δ 2 H and δ 18 O, we collected xylem tissue at 5 heights (0, 1, 4, 7 and 10 m) from three replicate trees of three species with contrasting xylem anatomy, two angiosperms (Liquidambar styraciflua and Quercus nigra) and one gymnosperm (Pinus taeda). Concurrently, we also determined δ 2 H and δ 18 O of groundwater, soil water, and recent precipitation. On a different day, we collected circumferential samples at 1 m from the same trees to test for sectorality effects. Water from stem and soil samples were cryogenically extracted and analyzed for δ 2 H and δ 18 O. Mean xylem water δ 2 H and δ 18 O were significantly different between species, both vertically and circumferentially. We did not find significant (p = 0.05) systematic variation in δ 2 H or δ 18 O with height. We found no significant evidence for sectorality effects on δ 2 H, δ 18 O. Variances of vertical synoptic xylem water δ 2 H or δ 18 O were similar between species, 7.2–10.4‰ for δ 2 H and 0.58–0.81‰ for δ 18 O. Circumferential variances were also similar between species, 4.0–6.0 for δ 2 H and 0.37–0.44 for δ 18 O. A mixing model showed that sweetgum, water oak and loblolly pine, were drawing most of their water from deep soil from 45 to 190 cm (84.7, 68.4 and 53.2%, respectively) however, soil water-excess values indicate δ 2 H fractionation effects on these estimates. Dual isotope mixing model evaluation with single and multiple sample configurations showed that source water estimates were not affected by within tree variability of xylem water signatures. Xylem water δ 2 H and δ 18 O variability with height or circumference, was 3.2 and 2.7 times less than between tree variation because of transient temporal and spatial processes and is, therefore, not likely to affect interpretations of water sourcing in these three species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A Type II Hamiltonian Variational Principle and Adjoint Systems for Lie Groups

We present a novel Type II variational principle on the cotangent bundle of a Lie group which enforces Type II boundary conditions, i.e., fixed initial position and final momentum. In general, such Type II variational principles are only globally defined on vector spaces or locally defined on general manifolds; however, by left translation, we are able to define this variational principle globally on cotangent bundles of Lie groups. Type II boundary conditions are particularly important for adjoint sensitivity analysis, which is our motivating application. As such, we additionally discuss adjoint systems on Lie groups, their properties, and how they can be used to solve optimization problems subject to dynamics on Lie groups.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

What explains the year-to-year variation in growing season timing of boreal black spruce forests?

Amplified climate warming in high latitudes is expected to affect growing season timing of the vast boreal biome. It is unclear whether the presence of permafrost (perennially frozen ground) might have an influence on changes in growing season timing. Here, this study examined how different environmental variables explained, either directly or indirectly, the variation in growing season timing of boreal forest stands with and without permafrost. We expected that environmental variables explaining the variation in growing season timing differed or had different explanatory power depending on permafrost presence or absence. The growing season was delineated from daily gross primary productivity (GPP) time series derived from 40 site-year data of net ecosystem carbon dioxide exchange measured with eddy covariance techniques over five black spruce (Picea mariana [Mill.])-dominated boreal forest stands in North America. In permafrost-free forest stands, a combination of start in canopy ‘green-up’ in spring and the timing of air and soil temperature increasing above freezing explained the start-of-season (SOS GPP ). Results from commonality analysis and structural equation modeling suggest that canopy ‘green-up’ and air temperature directly affected SOS GPP in permafrost-free forest stands. In addition, soil temperature acted as mediator for an indirect effect of air temperature on SOS GPP . In contrast, none of the environmental variables, or their combination, explained the variation in SOS GPP in forest stands with permafrost. The explanatory power of environmental variables was more consistent regarding the end-of-season (EOS GPP ). In both, forest stands with and without permafrost, EOS GPP was directly explained by mean soil water content in the fall and the first day of continuous snowpack formation. A better understanding how environmental variables control SOS GPP and EOS GPP in forest stands with and without permafrost will help to refine parameterizations of the boreal biome in Earth system models.

54 ENVIRONMENTAL SCIENCES↗