Search NASA⌕ Search

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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Validation of a Generalized Formulation for Load-Sharing Behavior in Epicyclic Gears for Wind Turbines: Preprint

In an ideal epicyclic gear set every parallel gear path transmits the same amount of torque. However, it is well known that certain manufacturing variations result in unequal load sharing between the parallel gear paths. Previous works have developed and validated a general closed-form analytical model of this phenomenon that describes the load sharing characteristics of epicyclic gear sets from three to six planets at any torque level. More recently, this analytical model has been reformulated to include the effects of gravity, carrier bearing clearance, and external applied moments, all of which are relevant to most gearboxes and their mounting configuration in horizontal-axis wind turbines. In this paper, the reformulated model is compared to load measurements collected from two similar wind turbine gearboxes with three-planet epicyclic gear sets. The resulting load sharing values are also compared to the mesh load factor requirements in the AGMA 6006 and IEC 61400-4 wind turbine gearbox design standards. Load sharing factors as high as 1.16 at extreme rotor moments and 1.08 with no moment were measured for the gearbox with cylindrical roller bearings, but the load sharing factor remained below 1.10 for the gearbox with tapered roller bearings The results show that in the wind turbine application, load sharing is not equal - even for three-planet systems with a floating central member because of the effects of gravity, rotor moments, and the resulting relative motions among the epicyclic gear components within the carrier bearing clearance.

closed-form analytical model↗
DOE OSTI2020

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Semicoherent symmetric quantum processes: Theory and applications

Discovering pragmatic and efficient approaches to construct ε-approximations of quantum operators such as real (imaginary) time-evolution propagators in terms of the basic quantum operations (gates) is challenging. Prior ε-approximations are invaluable, in that they enable the compilation of classical and quantum algorithm modeling of, e.g., dynamical and thermodynamic quantum properties. In parallel, symmetries are powerful tools concisely describing the fundamental laws of nature; the symmetric underpinnings of physical laws have consistently provided profound insights and substantially increased predictive power. In this work, we consider the interplay between the ε-approximate processes and the exact symmetries in a semicoherent context—where measurements occur at each logical clock cycle. Here we draw inspiration from Pascual Jordan's groundbreaking formulation of nonassociative, but commutative, symmetric algebraic form. Our symmetrized formalism is then applied in various domains such as quantum random walks, real-time evolutions, variational algorithm ansatzes, and efficient entanglement verification. Our work paves the way for a deeper understanding and greater appreciation of how symmetries can be used to control quantum dynamics in settings where coherence is a limited resource.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Optimizing Power Generation of a Bottom-Raised Oscillating Surge Wave Energy Converter Using a Theoretical Model: Preprint

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗
DOE OSTI2021

Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗
DOE OSTI2026

Variations of Alloying Site Density in Pd 1 Cu Single‐Atom Alloy Catalysts Lead to Shifted Product Yields in Electrochemical CO Reduction

Single-atom alloy (SAA) catalysis research often reports that a SAA catalyst, in the general formulation of a single-atom metal M1 alloyed on the surface of the host metal M2, facilitates a probe reaction. However, for catalytic reactions that present decoupled rate- and selectivity-limiting steps, the alloying site density may significantly manipulate these independent steps, but it has rarely been explicitly examined for any SAA systems. Herein, using the electrocatalytic CO reduction as a probe reaction, we report that the nominal Pd 1 Cu cube SAA catalysts exhibit distinctive high reactivity toward ethylene or ethanol, respectively, depending on whether the Pd atoms are in dilute or crowded forms. Although the presence of single-atom Pd embedded on Cu uniformly promotes CHO* formation and C─C coupling, the dilute-Pd 1 Cu favors ethylene formation by enabling low-barrier C─O cleavage from a flat CH 2 CH 2 OH* intermediate, whereas the crowded-Pd 1 Cu promotes ethanol formation by stabilizing an upright hydrogenation transition state of the same intermediate. Furthermore, we present evidence that the catalytic chemistry of crowded Pd 1 species differs from that of the Pd 2 -dimer; the latter, albeit unstable, steers reaction selectivity to acetate instead. In conclusion, these results uncovered the underappreciated importance of controlling SAA catalytic chemistry from the perspective of single-atom site densities.

Jin, Zehua [Clemson University, SC (United States)↗
DOE OSTI2022

Effect of lot microstructure variations on detonation performance of the triaminotrinitrobenzene (TATB)-Based insensitive high explosive PBX 9502

PBX 9502 is an important insensitive high explosive due to its combination of safety properties and detonation performance. It is a polymer-bonded formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the high explosive crystal, bound with Kel F-800 (FK-800), a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Two different types are used, one known as virgin that uses only pristine manufactured TATB, and a second known as recycled that has 50 wt.% of its TATB reclaimed from machining scraps of previously pressed virgin PBX 9502. Recycled lots have a higher percentage of fine particles compared to virgin lots, due to the fracturing and damage sustained by TATB crystals during pressing. We examine the influence that this TATB microstructure difference between virgin and recycled PBX 9502 has on detonation performance properties. New rate-stick geometry diameter effect, detonation front shape and cylinder expansion test data are obtained for two previously uncharacterized virgin lots and one recycled lot of PBX 9502. This is combined with previously published data for one virgin and one recycled lot for evaluation purposes. Detonation shock dynamics model calibrations are conducted on each of the five lots to provide an assessment of the detonation timing characteristics of virgin versus recycled PBX 9502 lots. For two of the virgin lots, detonations propagate slower in the rate-stick geometry than those in the two recycled lots. However, the other virgin lot tested has propagation rates comparable to that of the recycled lots for larger diameter rate-sticks. The latter result, though, is shown to be geometry dependent and depends on the range of detonation curvatures accessed in different geometries. New copper-confined, cylinder expansion tests are conducted on each of the three virgin and two recycled lots to obtain detonation product Jones-Wilkins-Lee equations of state, enabling an assessment of the metal push capabilities for each of the five lots. We find that the metal push capabilities, characterized by the evolution of the heat of detonation with volume, are similar between the virgin and recycled lots. Thus, changes in microstructure between different PBX 9502 lots seemingly affect the rate of reaction, but not the overall energy content.

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

Exploring the impact of nucleotide length on lipid nanoparticle structure and properties

Lipid nanoparticles (LNPs) are versatile carriers for nucleic acid (NA) therapeutics, including ASOs, siRNA, mRNA, and poly-IC. While lipid composition is known to influence LNP properties, the impact of NA length on morphology and internal structure is less understood, particularly during the stages of carrier–cargo assembly. Here, we examine NA length and lipid composition immediately after mixing using high-throughput SAXS, dynamic light scattering, and cryogenic electron microscopy. All LNPs form ordered NA/lipid compartments, with longer NAs promoting inverse hexagonal (H II ) phases and larger intercompartment distances. In contrast, short NAs, especially in formulations with SM102 ionizable lipid, favor lamellar phases. SAXS peak deconvolution quantifies ordered versus disordered phases via a Robustness of Ordered Phase factor, which correlates with particle size and encapsulation efficiency. Formulations with MC3 ionizable and DOPE helper lipids exhibit the most stable H II -phase packing, highlighting the role of helper-lipid curvature in compartment stabilization. Variations in NA compartmentalization indicate differences in payload capacity, offering a framework for rational LNP design across diverse nucleic acid cargos.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Plug-and-Play Methods for Integrating Physical and Learned Models in Computational Imaging: Theory, algorithms, and applications

Plug-and-play (PnP) priors constitute one of the most widely used frameworks for solving computational imaging problems through the integration of physical models and learned models. PnP leverages high-fidelity physical sensor models and powerful machine learning methods for prior modeling of data to provide state-of-the-art reconstruction algorithms. PnP algorithms alternate between minimizing a data fidelity term to promote data consistency and imposing a learned regularizer in the form of an image denoiser. Recent highly successful applications of PnP algorithms include biomicroscopy, computerized tomography (CT), magnetic resonance imaging (MRI), and joint ptychotomography. This article presents a unified and principled review of PnP by tracing its roots, describing its major variations, summarizing main results, and discussing applications in computational imaging. Additionally, we also point the way toward further developments by discussing recent results on equilibrium equations that formulate the problem associated with PnP algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Compositional Variation Tests on DuraMelter 100 with LAW Sub-Envelope C1 Feed (LAWC22) Glass in Support of the LAW Pilot Melter (Final Report)

The Vitreous State Laboratory at The Catholic University of America (VSL) is developing and testing glass formulations for RPP-WTP waste envelopes to provide data to meet the RPP-WTP contract requirements and to support system design activities. That work is based upon small-scale batch melts (“crucible melts”) using waste envelope simulants. Selected formulations have also been tested in small-scale continuously-fed joule-heated melters (DM10 and DM100 systems) and, ultimately, in the LAW Pilot Melter. Such melter tests provide information on key process factors such as feed processing behavior, dynamic effects during processing, sulfate incorporation, processing rates, off-gas amounts and compositions, foaming control, etc., that cannot be reliably obtained from crucible melts. This sequential scale-up approach in the vitrification testing program ensures that maximum benefit is obtained from the more costly pilot-scale tests and that the most effective use is made of that resource. The principal objective of the work described in this report was to collect the necessary small-scale melter test data with a LAW Sub-Envelope C1 waste simulant in order to support the next series of tests with the same melter feed on the LAW Pilot Melter. This work was conducted under a corresponding Test Specification and Test Plan.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Herbaceous Feedstock 2018 State of Technology Report

The U.S. Department of Energy (DOE) promotes the production of advanced liquid transportation fuels from lignocellulosic biomass by funding fundamental and applied research that advances the state of technology (SOT). As part of its involvement with this mission, Idaho National Laboratory (INL) completes an annual SOT report for biomass feedstock logistics. This report summarizes supply system impacts of Bioenergy Technologies Office (BETO)-funded research and development efforts at INL and elsewhere (such as the High-Tonnage Feedstock Logistics projects (Webb et al. 2013a, Webb et al. 2013b, Webb et al. 2013c, Webb and Sokhansanj 2014, Sokhansanj et al. 2014) that lead to improvements in feedstock supply systems. These include improvements to and observed performance of innovative harvest and collection methods, storage technologies, transportation and handling approaches, and advanced preprocessing technologies. Biomass quality and variability, and the interface between feedstock quality and conversion performance are key drivers in addition to delivered feedstock cost. In this report, we estimate the benefits of R&D improvements to individual supply system unit operations, and present the status of feedstock logistics technology development for converting biomass into biofuels. These analyses are supported by experimental data where possible, and help to align the SOT relative to the cost goals defined in the Multi-Year Program Plan. The 2018 Herbaceous SOT aligned feedstock logistic design with current biorefinery’s design capacity utilized by biochemical conversion platform. Currently biochemical conversion platform utilizes a 725,000 dry ton/year biorefiney design for the techno economic analysis. Hence, feedstock delivered cost in the 2018 Herbaceous SOT is calculated based on biorefinery’s 725,000 dry ton design capacity instead of 800, 000 dry ton capacity utilized in the 2017 Herbaceous SOT. Biomass availabilities in this SOT were updated to year 2018 data from the 2016 Billion-Ton Report (BT16) (DOE 2016a), with the exception of switchgrass, for which the 2018 Herbaceous SOT utilized the 2019 switchgrass availability data from BT16. The BT16 report (DOE 2016a) does not project switchgrass availability in 2018; the soonest switchgrass is available in the BT16 report is 2019. Therefore, availability of switchgrass for this analysis was that projected for 2019. The 2018 Herbaceous SOT incorporates same technologies utilized in the 2017 Herbaceous SOT. However, a sensitivity analysis is performed to understand the impact of variation of process parameters on those technologies on feedstock logistic cost. New R&D data that shows the variations of process parameters affecting process performance is incorporated in the 2018 SOT to measure the variations in delivered feedstock cost. The 2018 Herbaceous SOT has also provided projected delivered feedstock of 2022 design case based on near term technical target under BETO funded R&D project. Finally, updated biorefinery size of 725,000 dry ton/year was incorporated within least-cost formulation model to select optimal siting and depot scales during optimization of the least cost blend. This modification to the optimization algorithm allows the trade-off between the cost of increased supply radius and the savings from selecting biomass from higher producing counties to be assessed. Such optimization has also showed the economic benefit of decentralized depots in comparison to centralized preprocessing co-located with the biorefinery by decoupling the biorefinery and feedstock locations. The 2018 Herbaceous SOT report documents the current modeled cost of a herbaceous feedstock supply system (from harvest to the pretreatment reactor throat, including grower payment) for hydrocarbon fuel production via biochemical conversion, based on equipment and processes now available or potentially available in the near term. The modeled cost also considers both the required quality and the availability of the biomass resources. The 2018 Herbaceous SOT predicts a modeled delivered feedstock cost of $83.67/dry ton (2016$); this is a $0.23/dry ton (2016$) decrease from the 2017 Herbaceous SOT. The modification of biorefinery’s designed capacity and increased projected biomass availability in the same supply shed contributed to this modeled cost reduction. The least-cost formulation model to optimally site and scale local distributed preprocessing depots also contributed to the cost reduction by considering county-level grower payment and distance from the biorefinery as variables in the optimization algorithm. Sensitivity analysis on various process parameters that affect delivered feedstock cost in the 2018 Herbaceous SOT shows that the delivered cost could varies from $80.45-$88.83/dry ton. The top factors that causes such variations are: effective baling rate, bale density, hammer mill throughput, interest rate and storage dry matter loss.

09 BIOMASS FUELS↗
DOE OSTI2024

LAW Glass Formulation to Support Melter Runs with Simulants

The major objective of this work was to prepare and characterize glass formulations in order to develop and select formulations that are compliant with contractual and processing requirements to support melter tests with each of the LAW Sub-Envelopes. Glasses that cover composition ranges of interest were prepared and characterized to ensure that selected glass compositions can tolerate expected composition variations resulting from process errors without adversely affecting processing characteristics or product quality. Glass compositions that are expected to result during transitions from one waste stream to another were also prepared and characterized. Other objectives of the work included preparation and characterization of glasses with respect to the properties of interest, optimization of sulfate loading in the glasses, evaluation of ability to achieve waste loading limits, testing to demonstrate compatibility of glass melts with melter materials of construction, and identification of glass formulation issues with respect to contract specifications and processing requirements.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Quarkyonic matter equation of state in beta-equilibrium

Quark matter may appear due to a hadronic-quark transition in the core of a hybrid star. Quarkyonic matter is an approach in which both quarks and nucleons appear as quasiparticles in a crossover transition, and provides an explicit realization of early ideas concerning quark matter (e.g., the MIT bag model). This description has recently been employed by McLerran and Reddy to model chargeless (pure neutron) matter with an approach that has the virtue that the speed of sound rises quickly at a neutron-quark transition so as to satisfy observational constraints on the neutron star maximum mass (≳ 2 $M_⊙$) and the radius of a 1.4 $M_⊙$ star ($R_{1.4}$ ≲ 13.5 km). Traditional models involving first-order transitions result in softer pressure-energy density relations that have difficulty satisfying these constraints except with very narrow choices of parameters. Here, we propose a variation of quarkyonic matter involving protons and leptons whose energy can be explicitly minimized to achieve both chemical and beta equilibrium, which cannot be done in the chargeless formulation. Quarkyonic stellar models are able to satisfy observed mass and radius constraints with a wide range of model parameters, avoiding the obligatory fine-tuning of conventional hybrid star models, including requiring the transition density to be very close to the nuclear saturation density. Our formulation fits experimental and theoretical properties of the nuclear symmetry energy and pure neutron matter, and contains as few as three free parameters. This makes it an ideal tool for the study of high-density matter that is an efficient alternative to piecewise polytrope or spectral decomposition methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A constitutive structural parameter c b for the work hardening behavior of additively manufactured Ti-6Al-4V

The mechanical behavior of Ti-6Al-4V produced by additive manufacturing processes is assessed as based on a model derived from the Kocks–Mecking relationship. A constitutive parameter c b is derived from a linear Kocks–Mecking relationship for the microstructure that is characteristic of the work hardening behavior. The formulation for c b is determined by considering the plastic strain between the strengths at the proportional limit and the plastic instability. In this way, the model accommodates the variation in work hardening behavior observed when evaluating material as produced and tested along different orientations. The modeling approach is presented and evaluated for the case of Ti-6Al-4V additively manufactured materials as tested under quasi-static uniaxial tension. It is found that different test specimen orientations, along with postbuild heat treatments, produce a change in the microstructure and plasticity behavior which can be accounted for in the corresponding change of the c b values.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Partitioned exponential methods for coupled multiphysics systems

Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗