Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability”

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

The latent variable proximal point algorithm for variational problems with inequality constraints

The latent variable proximal point (LVPP) algorithm is a framework for solving infinite-dimensional variational problems with pointwise inequality constraints. The algorithm is a saddle point reformulation of the Bregman proximal point algorithm. At the continuous level, the two formulations are equivalent, but the saddle point formulation is more amenable to discretization because it introduces a structure-preserving transformation between a latent function space and the feasible set. Working in this latent space is much more convenient for enforcing inequality constraints than the feasible set, as discretizations can employ general linear combinations of suitable basis functions, and nonlinear solvers can involve general additive updates. LVPP yields numerical methods with observed mesh-independence for obstacle problems, contact, fracture, plasticity, and others besides; in many cases, for the first time. The framework also extends to more complex constraints, providing means to enforce convexity in the Monge–Ampère equation and handling quasi-variational inequalities, where the underlying constraint depends implicitly on the unknown solution. Here, in this paper, we describe the LVPP algorithm in a general form and apply it to ten problems from across mathematics.

Inequality constraints↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗
DOE OSTI2023

Economic evaluation of variable renewable energy participation in U.S. ancillary services markets

Variable renewable energy (VRE) is not yet meaningfully participating in U.S. ancillary services (AS) markets. VRE participation in AS markets could provide a new source of revenue for VRE resource owners to offset declining energy and capacity values and a new tool for power system operators to address emerging system constraints. This paper uses a price-taker dispatch model and historical prices to estimate the economic value of standalone and hybrid (battery-paired) VRE participation in AS markets, from the resource owner and electricity system perspectives, in each of the seven U.S. independent system operator and regional transmission organization (ISO/RTO) markets where ancillary service prices are set. Across ISO/RTO markets, average (2015–2019) simulated incremental revenues from power regulation market participation were 0.0–2.9 USD/MWh (+0–15% of revenue without participation) for standalone VRE owners and 1–33 USD/MWh (+1–69%) for hybrid VRE owners. However, ISO/RTO reserve markets are relatively thin and have the potential to become saturated by energy storage projects that are currently in ISO/RTO interconnection queues. In most markets, standalone and hybrid VRE could provide power regulation reserves during periods with high power regulation prices, suggesting that VRE participation in AS markets could have high system value. The analysis highlights the relevance of separate upward and downward power regulation products and indicates that ISOs/RTOs might consider initially focusing on enabling hybrid VRE provision of AS.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

A regional comparison of sub-daily flow variability in regulated and unregulated rivers in the United States

Regulating rivers for hydropower or other purposes can dramatically alter river flow patterns, including creating substantial changes in flow over short, minutes-to-hours-long timespans known as sub-daily flow variability (SDFV). The impacts of flexible hydropower production on flow and aquatic organisms are increasingly documented in research. However, the degree to which flow alteration relates to different hydropower operational modes in distinct geographical regions and seasons is not well understood. This study offers a methodology for regional- and species-appropriate evaluations of potential impacts of flow on fish based on sub-daily flow characteristics of hydropower operational modes. We analyzed 15-min discharge data between 2018 and 2021 from 69 USGS stream gages to compare SDFV in hydropeaking, run-of-river, and unregulated systems in the US Southeast and Pacific Northwest. Regulated systems exhibited significant SDFV downstream from hydropower facilities relative to unregulated systems, but specific impacts differed between regions. Regulated systems in the Southeast were characterized by high flow coefficients of variation and ratios (hydropeaking only) and extended durations of daily upramping flow phases. Regulated systems in the Pacific Northwest were characterized by many short flow phases per day and large portions of the day spent upramping. Pacific Northwest unregulated systems displayed the strongest seasonal flow patterns while Southeastern hydropeaking systems displayed the greatest SDFV. Given that SDFV impacts multiple dimensions of fish ecology, region-specific sub-daily flow signatures have important implications for understanding and mitigating potential community-, species-, and age-specific effects on fish in different parts of the country.

Fish↗
DOE OSTI2021

Decadal Variability of Southeast US Rainfall in an Eddying Global Coupled Model

Abstract Ocean variability is a dominant source of remote rainfall predictability, but in many cases the physical mechanisms driving this predictability are not fully understood. This study examines how ocean mesoscales (i.e., the Gulf Stream SST front) affect decadal Southeast US (SEUS) rainfall, arguing that the local imprint of large‐scale teleconnections is sensitive to resolved mesoscale features. Based on global coupled model experiments with eddying and eddy‐parameterizing ocean, we find that a resolved Gulf Stream improves localized rainfall and remote circulation response in the SEUS. The eddying model generally improves the air‐sea interactions in the Gulf Stream and the North Atlantic Subtropical High that modulate SEUS rainfall over decadal timescales. The eddy‐parameterizing simulation fails to capture the sharp SST gradient associated with the Gulf Stream and overestimates the role of tropical Pacific SST anomalies in the SEUS rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Spring Dust in Colorado Plateau: Transport Pathways and Interannual Variability Derived From Two Decades of MERRA‐2 Reanalysis

Spring dust storms in the U.S. Southwest significantly impact environmental and human systems, yet their climatological patterns and driving mechanisms remain poorly understood. Using two decades of MERRA-2 reanalysis data and self-organizing map clustering, we identified four distinct dust transport pathways from surrounding and remote deserts for the top 10% of spring dust events impacting the Colorado Plateau: (a) intense southwesterly near-surface transport where eastward upward transport is blocked by topography, (b) substantial remote transport plus weaker surface emissions, (c) moderate southwesterly transport with stronger upper-level transport, and (d) rare northwesterly pattern linked to the Great Basin sources, which accounts for 38%, 14%, 37%, and 11% of total events, respectively. The frequency and intensity of these events strongly correlate with large-scale climate variability, with increased dust events during La Niña. Our results provide a robust framework for improving seasonal dust forecasts and understanding future dust dynamics under climate change.

aerosol transport↗
DOE OSTI2022

XMM-Newton and TESS observations of the highly variable polar V496 UMa

We studied the temporal and spectral behavior of V496 UMa from the optical to the X-ray regimes. Methods. We used archival XMM-Newton and TESS observations obtained in 2017 and 2019, respectively, to perform a spectral and timing analysis of the highly variable magnetic CV. Results. The light curves of the TESS and XMM-Newton satellites reveal a double-humped pattern modulated with the periodicity of 91.058467 ± 0.00001 min. V496 UMa displays a two-pole accretion geometry in the high accretion state. X-ray spectra from these regions are composed of thermal plasma radiation and soft blackbody components with almost identical temperatures and a total accretion rate of $\dot{M}$ = 1.4(8) × 10 -11 M ⊙ yr -1 . The X-ray centers of the humps show longitudinal shifts of -18° and 4°, and shifts around photometric phase zero of -172° and -186°, for the main hump and second hump, respectively. The long-term ZTF light curves reveal high and low accretions states. Low-state ZTF and SDSS photometric data are consistent with an 0.8 M ⊙ white dwarf at 10 000 K and a main-sequence donor star with a spectral type of M5.0 at a Gaia determined distance of 758 pc. Conclusions. V496 UMa is a very bright polar in X-rays when it is in the high state. Due to its unusual geometric structure, mass accretion onto the second accretion pole is interrupted occasionally. This discontinuous behavior does not follow a certain pattern in time and has been observed so far only in the high state. The X-ray light curves display clear evidence of an accretion stream at the photometric phase of φ = 0.81, which does not show up in optical light curves. An accurate period was derived using the combined TESS and XMM-Newton data, which differs by 3.8 σ from published results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. Here, the non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH 3 demonstrate the accuracy of the non-hierarchical CDVR approach.

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

Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation

Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Dynamic global vegetation models underestimate net CO 2 flux mean and inter-annual variability in dryland ecosystems

Despite their sparse vegetation, dryland regions exert a huge influence over global biogeochemical cycles because they cover more than 40% of the world surface (Schimel 2010 Science 327 418–9). It is thought that drylands dominate the inter-annual variability (IAV) and long-term trend in the global carbon (C) cycle (Poulter et al 2014 Nature 509 600–3, Ahlstrom et al 2015 Science 348 895–9, Zhang et al 2018 Glob. Change Biol. 24 3954–68). Projections of the global land C sink therefore rely on accurate representation of dryland C cycle processes; however, the dynamic global vegetation models (DGVMs) used in future projections have rarely been evaluated against dryland C flux data. Here, we carried out an evaluation of 14 DGVMs (TRENDY v7) against net ecosystem exchange (NEE) data from 12 dryland flux sites in the southwestern US encompassing a range of ecosystem types (forests, shrub- and grasslands). We find that all the models underestimate both mean annual C uptake/release as well as the magnitude of NEE IAV, suggesting that improvements in representing dryland regions may improve global C cycle projections. Across all models, the sensitivity and timing of ecosystem C uptake to plant available moisture was at fault. Spring biases in gross primary production (GPP) dominate the underestimate of mean annual NEE, whereas models' lack of GPP response to water availability in both spring and summer monsoon are responsible for inability to capture NEE IAV. Errors in GPP moisture sensitivity at high elevation forested sites were more prominent during the spring, while errors at the low elevation shrub and grass-dominated sites were more important during the monsoon. We propose a range of hypotheses for why model GPP does not respond sufficiently to changing water availability that can serve as a guide for future dryland DGVM developments. Our analysis suggests that improvements in modeling C cycle processes across more than a quarter of the Earth's land surface could be achieved by addressing the moisture sensitivity of dryland C uptake.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Continuous variable port-based teleportation

Abstract Port-based teleportation (PBT) is a generalisation of the standard teleportation protocol which does not require unitary operations by the receiver. This comes at the price of requiring N > 1 entangled pairs, while N = 1 for the standard teleportation protocol. The lack of correction unitaries allows PBT to be used as a fundamental theoretical tool to simulate arbitrary channels with a general resource, with applications to study fundamental limits of quantum communication, cryptography and sensing, and to define general programmable quantum computers. Here we introduce a general formulation of port-based teleportation in continuous variable systems and study in detail the N = 2 case. In particular, we interpret the resulting channel as an energy truncation and analyse the kinds of channels that can be naturally simulated after this restriction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Assessing decadal variability of subseasonal forecasts of opportunity using explainable AI

Abstract Identifying predictable states of the climate system allows for enhanced prediction skill on the generally low-skill subseasonal timescale via forecasts with higher confidence and accuracy, known as forecasts of opportunity. This study takes a neural network approach to explore decadal variability of subseasonal predictability, particularly during forecasts of opportunity. Specifically, this work quantifies subseasonal prediction skill provided by the tropics within the Community Earth System Model Version 2 (CESM2) Large Ensemble and assesses how this skill evolves on decadal timescales. Utilizing the networks’ confidence and explainable artificial intelligence, physically meaningful sources of predictability associated with periods of enhanced skill are identified. Using these networks, we find that tropically-driven subseasonal predictability varies on decadal timescales during forecasts of opportunity. Further, we investigate the drivers of the low frequency modulation of the tropical-extratropical teleconnection and discuss the implications. Analysis is extended to ECMWF Reanalysis v5 data, revealing that the relationships learned within the CESM2-Large Ensemble holds in modern reanalysis data. These results indicate that the neural networks are capable of identifying predictable decadal states of the climate system within CESM2 that are useful for making confident, accurate subseasonal precipitation predictions in the real world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Variable stars in Local Group galaxies – V. The fast and early evolution of the low-mass Eridanus II dSph galaxy

ABSTRACT We present a detailed study of the variable star population of Eridanus II (Eri II), an ultra-faint dwarf galaxy that lies close to the Milky Way virial radius. We analyse multi-epoch g, r, i ground-based data from Goodman and the Dark Energy Camera, plus F475W, F606W, F814W space data from the Advanced Camera for Surveys. We report the detection of 67 RR Lyrae (RRL) stars and 2 Anomalous Cepheids, most of them new discoveries. With the RRL stars, we measure the distance modulus of Eri II, μ0 = 22.84 ± 0.05 mag (D⊙ = 370 ± 9 kpc), and derive a metallicity spread of 0.3 dex (0.2 dex intrinsic). The colour distribution of the horizontal branch (HB) and the period distribution of the RRL stars can be nicely reproduced by a combination of two stellar models of [Fe/H] = (−2.62, −2.14). The overall low metallicity is consistent with the red giant branch bump location, 0.65 mag brighter than the HB. These results are in agreement with previous spectroscopic studies. The more metal-rich RRL and the RRab stars have greater central concentration than the more metal-poor RRL and the RRc stars that are mainly located outside ∼1 rh. This is similar to what is found in larger dwarf galaxies such as Sculptor, and in agreement with an outside-in galaxy formation scenario. This is remarkable in such a faint dwarf galaxy with an apparently single and extremely short (<1 Gyr) star formation burst. Finally, we have derived new and independent structural parameters for Eri II and its star cluster using our new data that are in very good agreement with previous estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Wavelength conversion for single-photon polarization qubits through continuous-variable quantum teleportation

A quantum internet connects remote quantum processors that need to interact and exchange quantum signals over a long distance through photonic channels. However, these quantum nodes operate at the wavelength ranges unsuitable for long-distance transmission. Therefore, quantum wavelength conversion to telecom bands is crucial for long-distance quantum networks based on optical fiber. Here, we propose wavelength conversion devices for single-photon polarization qubits using continuous-variable quantum teleportation that can efficiently convert qubits between near-infrared (780–795 nm suitable for interacting with atomic quantum nodes) and telecom wavelength (1300–1500 nm suitable for long-distance transmission). The teleportation uses entangled photon fields (i.e., nondegenerate two-mode squeezed state) that can be generated by four-wave mixing in a rubidium atomic gas using a diamond configuration of atomic transitions. The entangled fields can be emitted in two orthogonal polarizations with locked relative phase, making them especially suitable for interfacing with single- photon polarization qubits. Furthermore, our work may pave the way for the realization of long-distance quantum networks.

74 ATOMIC AND MOLECULAR PHYSICS↗