SEARCH · Search NASA

Results for “Variational formulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A constitutive model for elastomers tailored by ionic bonds and entanglements

Over the past decade or two, the concept has emerged of using multiple types of weak interactions simultaneously to enhance the mechanical properties of elastomers. These weak interactions include physical entanglements, hydrogen bonds, metal-coordination bonds, dynamic covalent bonds, and ionic bonds. The combination of entanglements and ionic bonding has been minimally explored and is particularly exciting because of the broad application space for polyelectrolytes. In this work, a constitutive model framework is developed to describe the response of elastomers with both ionic bonds and entanglements. Here we formulate a micromechanical model that couples together chain stretching, ionic bond slipping, and entanglement evolution. The ionic bonds provide toughness by enabling plastic deformation in comparison to covalently crosslinked material and add strength compared to a linear polymer. Evolution of the entanglement density is taken as a key mechanism that can govern stiffness, toughness, and self-recovery in elastomers. The model is used to match bulk polyelectrolytes with different fractions of ionic components under a variety of loading histories. The variations in material parameters are then used to help understand the relative importance of different governing mechanisms in the bulk polymers. We show that the theoretical framework can explain our experimental uniaxial tensile experimental results for polyelectrolytes. This model can help to design better material with high stiffness and toughness. We expect that our model can be extended to explain the mechanical behavior of other polyelectrolytes and other soft materials with a wide range of dynamic bonds.

42 ENGINEERING↗

DOE OSTI2022

Online eco-routing for electric vehicles using combinatorial multi-armed bandit with estimated covariance

Identifying energy-efficient routes in real-time has significant implications for the energy-optimal operations of electric vehicles (EVs). Here, this study proposes a novel model for EV online eco-routing problem, which obtains the minimal expected energy consumption paths (MECPs) for multiple origin-destination (OD) pairs simultaneously. Specifically, we formulate the routing problem as a bandit problem and solve it with online algorithms. We extend the algorithms by implementing a path elimination mechanism to reduce the candidate path set and introducing the variance and covariance of the energy consumption to reduce the uncertainties. The numerical results show that the proposed algorithms can efficiently obtain near-optimal MECPs, and the solution is significantly better than the widely used shortest trip time path algorithm (STTP) and shortest trip distance path algorithm (SDP). The variation considering link energy covariance and path elimination generates paths that save 4.1% of energy compared to the SDP and 5.4% to the STTP.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

The photon content of the proton in the CT18 global analysis

Recently, two photon PDF sets based on implementations of the LUX ansatz into the CT18 global analysis were released. In CT18lux, the photon PDF is calculated directly using the LUX master formula for all scales, \mu μ . In an alternative realization, CT18qed, the photon PDF is initialized at the starting scale, \mu_0 μ 0 , using the LUX formulation and evolved to higher scales \mu(>\mu_0) μ ( > μ 0 ) with a combined QED+QCD kernel at \mathcal{O}(\alpha),~\mathcal{O}(\alpha\alpha_s) 𝒪 ( α ) , 𝒪 ( α α s ) and \mathcal{O}(\alpha^2) 𝒪 ( α 2 ) . In the small- x x region, the photon PDF uncertainty is mainly induced by the quark and gluon PDFs, through the perturbative DIS structure functions. In comparison, the large- x x photon uncertainty comes from various low-energy, nonperturbative contributions, including variations of the inelastic structure functions in the resonance and continuum regions, higher-twist and target-mass corrections, and elastic electromagnetic form factors of the proton. We take the production of doubly-charged Higgs pairs, (H^{++}H^{--}) ( H + + H − − ) , as an example of scenarios beyond the Standard Model to illustrate the phenomenological implications of these photon PDFs at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

NFD Methodology Report

The NetFlow Dynamics (NFD) model was developed for estimating the availability of a commodity supplied by a national- or regional-scale infrastructure following unexpected disruption of one or more of its components. The large scope of the disruptions of interest produce changes in availability lasting days to weeks. Consequently, the model does not resolve daily variations in system state and does not include the long-term processes that cause infrastructures to evolve as assets are added and removed according to owners ’planning decisions. NFD simulates fluid flow, including petroleum and other incompressible fluids, as well as natural gas and other compressible fluids, through pipeline networks characterized by limits on transmission capacity and storage. It was designed to enable efficient exploration of possible transmission system responses to large-scale disruptions lasting for days or longer. The model formulation reflects constraints on transmission and storage capacity imposed by the physical system assets. Those capacity limits are input parameters and are not derived from more basic system properties such as pipeline diameters and compressor power. A system’s response to a large disruption is controlled by operational decisions as well as damage to physical assets. The NFD model formulation allows users to efficiently consider alternative scenarios about the way remaining capacity might be used so that the analysis result appropriately reflects uncertainties about operator response.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Droplet bioprinting of acellular and cell-laden structures at high-resolutions

Advances in digital light projection(DLP) based (bio) printers have made printing of intricate structures at high resolution possible using a wide range of photosensitive bioinks. A typical setup of a DLP bioprinter includes a vat or reservoir filled with liquid bioink, which presents challenges in terms of cost associated with bioink synthesis, high waste, and gravity-induced cell settling, contaminations, or variation in bioink viscosity during the printing process. Here, we report a vat-free, low-volume, waste-free droplet bioprinting method capable of rapidly printing 3D soft structures at high resolution using model bioinks and model cells. A multiphase many-body dissipative particle dynamics model was developed to simulate the dynamic process of droplet-based DLP printing and elucidate the roles of surface wettability and bioink viscosity. Process variables such as light intensity, photo-initiator concentration, and bioink formulations were optimized to print 3D soft structures (∼0.4–3 kPa) with a typical layer thickness of 50 µm, an XY resolution of 38 ± 1.5 μm and Z resolution of 237 ± 5.4 µm. To demonstrate its versatility, droplet bioprinting was used to print a range of acellular 3D structures such as a lattice cube, a Mayan pyramid, a heart-shaped structure, and a microfluidic chip with endothelialized channels. Droplet bioprinting, performed using model C3H/10T1/2 cells, exhibited high viability (90%) and cell spreading. Additionally, microfluidic devices with internal channel networks lined with endothelial cells showed robust monolayer formation while osteoblast-laden constructs showed mineral deposition upon osteogenic induction. Overall, droplet bioprinting could be a low-cost, no-waste, easy-to-use, method to make customized bioprinted constructs for a range of biomedical applications.

DLP↗

DOE OSTI2021

Multisublattice cluster expansion study of short-range ordering in iron-substituted strontium titanate

Owing to the challenges in obtaining realistic atomic configurations in large chemical phase spaces, it is not straightforward to describe structure–property relations in materials exhibiting configurational disorder. One example is iron-substituted strontium titanate (SrTi 1–x Fe x O 3–d , STF), a promising perovskite-derivative cathode material in solid oxide fuel cells that exhibits full solid solubility 0 ≤ x ≤ 1 and a tendency to exhibit short-range order. Here we demonstrate a multisublattice cluster expansion (CE) framework and apply it to STF across the full composition range. The CE approach is distinct from more traditional CE formulations in that clusters are defined explicitly by the chemical species distributed among multiple sublattices, rather than via cluster functions of occupation variables with decoration. The modified CE approach makes it easy to distinguish meaningful chemical interactions that are harder to extract from conventional CE, since for the latter chemical identity in a cluster is expressed as a product of site occupations. The least absolute shrinkage and selection operator (LASSO) is implemented as a regression analysis tool to select key clusters and avoid overfitting. We demonstrate this formulation on STF, and show that it can accurately predict configurational energies in comparison to conventional CE. From the key clusters, we identify that short-range ordering between substitutional Fe and oxygen vacancies (V O ) results in the formation of Fe–VO strings. In addition, we consider the stability of STF through CE-based Monte Carlo (MC) simulations and confirm the presence of superstructures that were previously observed in transmission electron microscopy. In this work, analysis of atomic configurations from MC samples reveals variations in the oxidation state of Fe atoms, which can be explained by the ordering tendency of Fe and V O . The cluster description and selection formalism described here may be applied to other disordered multisublattice systems for accurate and efficient material modeling.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Galerkin Neural Networks: A Framework for Approximating Variational Equations with Error Control

Herein, we present a new approach to using neural networks to approximate the solutions of variational equations, based on the adaptive construction of a sequence of finite-dimensional sub-spaces whose basis functions are realizations of a sequence of neural networks. Here, the finite-dimensional subspaces are then used to define a standard Galerkin approximation of the variational equation. This approach enjoys a number of advantages, including: the sequential nature of the algorithm offers a systematic approach to enhancing the accuracy of a given approximation; the sequential enhancements provide a useful indicator for the error that can be used as a criterion for terminating the sequential updates; the basic approach is largely oblivious to the nature of the partial differential equation under consideration; and, some basic theoretical results are presented regarding the convergence (or otherwise) of the method which are used to formulate basic guidelines for applying the method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Investigation of 3D printed lightweight hybrid composites via theoretical modeling and machine learning

Hybrid composites combine two or more different fillers to achieve multifunctional or advanced material properties, such as lightweight and enhanced mechanical properties. The properties of the composites significantly depend on their microstructures, which can be tailored via advanced 3D printing processes. Understanding the process-structure-property relationships is critical to enable the design and engineering of novel hybrid composites for applications in aerospace, automotive, and protective coatings. Here, for this work, we develop 3D printable and lightweight hybrid composites and leverage the conventional design of experiments, a theoretical hybrid model, and an image-driven machine learning (ML) method to investigate their mechanical behaviors. The hybrid composites are formulated with elastomer matrix, microfillers, and thin-shell particles, enabling a significant degree of design freedom of microstructures with densities and mechanical properties varying up to 70% and 91%, respectively. Our statistical analysis indicates that the 3D printing path direction and the microfibers fraction are dominating process parameters with contribution percentages of 45.3% and 57.7% on the specific stiffness and strength, respectively. A hybrid mechanics model is developed based on a simple Weibull distribution function and classical single-filler models to effectively capture the variations in mechanical properties, however, it overestimates the values due to its statistical constraints and idealization of experimental uncertainty. The image-driven ML model leverages the microscale images directly without losing the structural details, shows more accurate predictions with experimental data, and has 48.6% lower root mean square error than the theoretical model.

3D printing↗

DOE OSTI2021

U.S. Solar Photovoltaic System and Energy Storage Cost Benchmark: Q1 2021 [Slides]

Our benchmarking method includes bottom-up accounting for all necessary system and projectdevelopment costs incurred when installing residential, commercial, and utility-scale systems, and it models the Q1 2021 costs for such systems, excluding any previous supply agreements or contracts. In general, we attempt to model the typical installation techniques and business operations from an installed-cost perspective, and our benchmarks are national averages. The residential PV-only benchmark and the commercial rooftop PV-only benchmark average costs by inverter type (string inverters, string inverters with direct current [DC] optimizers, and microinverters), weighted by inverter market share. The residential PV-only benchmark is further averaged across small installer and national integrator business models, weighted by market share. All benchmarks include variations—accounting for the differences in size, equipment, and operational use (particularly for storage)—that are currently available in the marketplace. All benchmarks assume nonunionized construction labor; residential and commercial PV systems predominantly use nonunionized labor, and the type of labor required for utility-scale PV systems depends heavily on the development process. All benchmarks assume the use of monofacial monocrystalline silicon PV modules. Benchmarks using cadmium telluride or bifacial modules could result in significantly different results. The data in this annual benchmark report inform the formulation of and track progress toward the U.S. Department of Energy (DOE) Solar Energy Technologies Office’s Government Performance and Reporting Act cost targets.

14 SOLAR ENERGY↗

DOE OSTI2023

Exploring Parameter Redundancy in the Unitary Coupled-Cluster Ansätze for Hybrid Variational Quantum Computing

One of the commonly used chemical-inspired approaches in variational quantum computing is the unitary coupled-cluster (UCC) ansatze. Despite being a systematic way of approaching the exact limit, the number of parameters in the standard UCC ansatze exhibits unfavorable scaling with respect to the system size, hindering its practical use on near-term quantum devices. Efforts have been taken to propose some variants of UCC ansatze with better scaling. In this paper we explore the parameter redundancy in the preparation of unitary coupled-cluster singles and doubles (UCCSD) ansatze employing spin-adapted formulation, small amplitude filtration, and entropy-based orbital selection approaches. Numerical results of using our approach on some small molecules have exhibited a significant cost reduction in the number of parameters to be optimized and in the time to convergence compared with conventional UCCSD-VQE simulations. Further, we also discuss the potential application of some machine learning techniques in further exploring the parameter redundancy, providing a possible direction for future studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Accurate Ultrasonic Thickness Measurement for Arbitrary Time-Variant Thermal Profile

Ultrasonic thickness measurement of mechanical structures is one of the most popular and commonly used nondestructive methods for various kinds of process control and corrosion monitoring. With ultrasonic propagation speed being temperature-dependent, the thickness measurement can be performed reliably only when the thermal profile is completely known. Most conventional techniques assume the temperature of the test structure is uniform and at room temperature across its thickness. Such assumptions may lead to large errors in the thickness measurement, especially when there are significant temperature variations across the thickness. State-of-the-art techniques use external temperature measurements or implement iterative methods to compensate for the unknown thermal profiles. However, such techniques produce unsatisfactory results when the heat distribution is complex or varies rapidly with time. In this work, we propose a two-sensors technique, using both compressive and shear excitations, with a non-iterative rapid data processing method for accurate thickness measurement under arbitrary time-variant thermal profile. The independent behavior of shear and compressive waves is used to formulate a real-time thickness estimation technique. The developed technique is experimentally validated on a steel plate with fixed acoustic sensors. Test results show that the error in thickness estimation can be reduced by up to 98% compared to conventional thickness gauging methods.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Nanoscale Photoexcited Carrier Dynamics in Perovskites

The optoelectronic properties of lead halide perovskite thin films can be tuned through compositional variations and strain, but the associated nanocrystalline structure makes it difficult to untangle the link between composition, processing conditions, and ultimately material properties and degradation. Here, we study the effect of processing conditions and degradation on the local photoconductivity dynamics in [(CsPbI 3 ) 0.05 (FAPbI 3 ) 0.85 (MAPbBr 3 ) 0.15 ] and (FA 0.7 Cs 0.3 PbI 3 ) perovskite thin films using temporally and spectrally resolved microwave near-field microscopy with a temporal resolution as high as 5 ns and a spatial resolution better than 50 nm. For the latter FACs formulation, we find a clear effect of the process annealing temperature on film morphology, stability, and spatial photoconductivity distribution. After exposure of samples to ambient conditions and illumination, we find spectral evidence of halide segregation-induced degradation below the instrument resolution limit for the mixed halide formulation, while we find a clear spatially inhomogeneous increase in the carrier lifetime for the FACs formulation annealed at 180 degrees C.

14 SOLAR ENERGY↗

DOE OSTI2022

CAPG: comprehensive allopolyploid genotyper

Genotyping by sequencing is a powerful tool for investigating genetic variation in plants, but many economically important plants are allopolyploids, where homoeologous similarity obscures the subgenomic origin of reads and confounds allelic and homoeologous SNPs. Recent polyploid genotyping methods use allelic frequencies, rate of heterozygosity, parental cross or other information to resolve read assignment, but good subgenomic references offer the most direct information. The typical strategy aligns reads to the joint reference, performs diploid genotyping within each subgenome, and filters the results, but persistent read misassignment results in an excess of false heterozygous calls. We introduce the Comprehensive Allopolyploid Genotyper (CAPG), which formulates an explicit likelihood to weight read alignments against both subgenomic references and genotype individual allopolyploids from whole-genome resequencing data. We demonstrate CAPG in allotetraploids, where it performs better than Genome Analysis Toolkit’s HaplotypeCaller applied to reads aligned to the combined subgenomic references.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Experimental test of model predictive control in a variable air volume system

Model predictive control (MPC) has been widely studied as a promising approach for improving energy efficiency and operational flexibility in buildings, yet its real-world performance for commercial variable air volume (VAV) systems remains insufficiently characterized. In particular, the impacts of model mismatch on control robustness, real-time computational burden, and device-level operation are rarely evaluated using long-term field data. Here, this study presents a comprehensive experimental evaluation of MPC applied to a full-scale VAV system in Oak Ridge National Laboratory’s Flexible Research Platform-2 building with constant cooling/heating temperature setpoints and no occupancy. The study offers three key advantages over existing work: (1) it uses a representative building in a full-scale experimental test, capturing realistic system dynamics and complexity; (2) it evaluates a relatively sophisticated MPC formulation using two different optimization solvers (Gurobi and PSO), fully accounting for computational complexity and methodological diversity; and (3) it systematically assesses potential negative impacts on various building devices, benchmark against a well-established baseline, ASHRAE Guideline 36 (G36). To isolate zone- and air-handling-unit–level supervisory control effects, the supply fan was operated with a fixed static pressure setpoint under all strategies, and the trim-and-response static pressure reset in G36 was not enabled. Results show that MPC maintained thermal comfort while improving energy efficiency. Abrupt solar radiation variations degraded performance. Computation times ranged from ∼1 s (Gurobi) to ∼ 70 s (PSO). Compared with G36, MPC achieves 33% energy savings and reduces median reheat coil output by approximately a factor of 5–10 for a representative cooling day under matched weather conditions. However, it increases the maximum discomfort deviation from 0.5 to 1°C and results in a 32% increase in staging frequency. In addition, PSO-based MPC introduced damper oscillations, also affecting actuator longevity.

ASHRAE guideline 36↗

DOE OSTI2022

VQE method: a short survey and recent developments

Abstract The variational quantum eigensolver (VQE) is a method that uses a hybrid quantum-classical computational approach to find eigenvalues of a Hamiltonian. VQE has been proposed as an alternative to fully quantum algorithms such as quantum phase estimation (QPE) because fully quantum algorithms require quantum hardware that will not be accessible in the near future. VQE has been successfully applied to solve the electronic Schrödinger equation for a variety of small molecules. However, the scalability of this method is limited by two factors: the complexity of the quantum circuits and the complexity of the classical optimization problem. Both of these factors are affected by the choice of the variational ansatz used to represent the trial wave function. Hence, the construction of an efficient ansatz is an active area of research. Put another way, modern quantum computers are not capable of executing deep quantum circuits produced by using currently available ansatzes for problems that map onto more than several qubits. In this review, we present recent developments in the field of designing efficient ansatzes that fall into two categories—chemistry–inspired and hardware–efficient—that produce quantum circuits that are easier to run on modern hardware. We discuss the shortfalls of ansatzes originally formulated for VQE simulations, how they are addressed in more sophisticated methods, and the potential ways for further improvements.

Fedorov, Dmitry A. (ORCID:0000000316598580)↗

DOE OSTI2025

Study of Quadrupole Fringe Fields in the Interaction Region of the Hadron Storage Ring of the Electron Ion Collider

Fringe fields in quadrupole magnets are usually neglected in studies of beam dynamics at accelerators. However, the extreme optical parameters present in the final focus of a collider such as the Electron–Ion Collider (EIC) may give rise to effects that should not be overlooked. The calculation of quadrupole fringe fields presented in this study follows the procedure outlined in Ref. [1], specialized to the case of a straight reference orbit (i.e., with no dipole field component). A right-handed Cartesian coordinate system is employed, with the $z$-axis aligned with the quadrupole axis and $x$ and $y$ denoting the horizontal and vertical transverse coordinates, respectively. The magnetic quadrupole field gradient, $G = \partial B_y / \partial x$, transitions from its peak value inside the quadrupole—where it is nearly independent of the longitudinal coordinate $z$—to zero at some distance beyond the magnet edge. Consequently, $G$ is treated as a function of $z$. The region over which this variation occurs is defined as the quadrupole fringe field region. The study begins with the development of a description of the magnetic field in the fringe region using a power-series expansion in the transverse coordinates $x$and $y$, consistent with the longitudinally varying gradient. A model for the $z$-dependence of the gradient is then proposed and adjusted to reproduce magnetic field data obtained from three-dimensional field calculations. To evaluate the impact of the fringe fields on beam dynamics, the corresponding vector potential is derived and incorporated into a Hamiltonian formulation of particle motion. The significance of the fringe fields is quantified by calculating the amplitude-dependent tune shift from the Hamiltonian. Using linear beam optics parameters of the Hadron Storage Ring (HSR) of the EIC, the tune shift due to the fringe fields of all quadrupole magnets in the IR-6 interaction region is evaluated. Finally, the resulting tune shifts are compared with those arising from other nonlinear field components present in the HSR.

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

Effects of Media Nutrient Variation on Microalgae Productivity and Economics During Semi-Continuous Cultivation

The development of large-scale microalgae growth for biofuel production is currently limited by the cost of biomass production. However, new approaches to infrastructure and cultivation practices are bringing the field closer to realization. Macronutrients in the cultivation media contribute significant costs, especially since their concentrations have not been optimized for specific strains and conditions. Environmental photobioreactors (ePBRs) were used to simulate cultivation under outdoor conditions, during which the nitrogen and phosphorus levels in the media were varied. The growth of two potential biofuel production strains, Picochlorum celeri and Tetraselmis striata, with varying nutrient inputs during summer and winter scripts, respectively, was studied. This study demonstrated that nitrogen and phosphorus in f/2 media could be reduced by more than 60% from the standard formulation, while maintaining growth rates in a semi-continuous harvesting approach. Experiments comparing the standard and reduced nutrient input concentrations were also conducted for both species in 820 L outdoor raceway ponds, in Mesa, AZ. P. celeri grown in these ponds in October had a growth rate of 10.6 ± 0.7 g/m 2 /day and 10.6 ± 0.3 g/m 2 /day for the standard and low-nutrient P. celeri ponds, respectively. T. striata grown in April–May had a growth rate of 16.6 ± 1.4 g/m 2 /day for the standard nutrient input ponds and 17.4 ± 1.1 g/m 2 /day for the low-nutrient input ponds, and in October 14.5 ± 0.6 g/m 2 /day for standard nutrient ponds and 14.4 ± 0.6 g/m 2 /day for low-nutrient ponds. These outdoor data therefore confirmed the indoor ePBR data. Techno-economic analysis shows that, if high growth rates can be attained at lower nutrient concentrations, a reduction of at least 60% in nutrient costs can be achieved. Such results highlight the importance of managing macronutrient media inputs, as these have a considerable contribution to biomass production costs in large-scale facilities. The analysis also points to the importance of maintaining high spent medium recycling rates in an industrial deployment, so as to minimize the losses of nitrogen and phosphorus compounds.

09 BIOMASS FUELS↗