Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “APPROXIMATIONS”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Network-Cognizant Time-Coupled Aggregate Flexibility of Distribution Systems Under Uncertainties

Increasing integration of distributed energy resources (DERs) within distribution feeders provides unprecedented flexibility at the distribution-transmission interconnection. To exploit this flexibility and to use the capacity potential of aggregate DERs, feasible substation power injection trajectories need to be efficiently characterized. This paper provides an ellipsoidal inner approximation of the set of feasible power injection trajectories at the substation such that for any point in the set, there exists a feasible disaggregation strategy of DERs for any load uncertainty realization. The problem is formulated as one of finding the robust maximum volume ellipsoid inside the flexibility region under uncertainty. Though the problem is NP-hard even in the deterministic case, this paper derives novel approximations of the resulting adaptive robust optimization problem based on optimal second-stage policies. The proposed approach yields less conservative flexibility characterization than existing flexibility region approximation formulations. The efficacy of the proposed method is demonstrated on a realistic distribution feeder.

adaptive robust optimization↗
DOE OSTI2025

Thermal Performance of a Conduction-Cooled CCT Dipole ReBCO Magnet: Several Cycles of Cool-Down and Thermal Gradient Measurements

Here, this paper presents experimental results from conduction-cooled thermal testing of a ReBCO canted cosine theta (CCT) magnet (C2), originally designed and fabricated at LBNL using CORC cables. While the performance of the coil under liquid helium and nitrogen environments has been previously established, this study explores its behavior under conduction cooling using a large test cryostat at The Ohio State University. The magnet, measuring 613 mm in length and weighing ∼75 kg, consists of four helical layers wound with ReBCO-based CORC wire and was thermally anchored to a copper cold ring supported by a G-10 strongback. Cooling was provided by two Sumitomo RDK-415D cryocoolers, offering a combined 3 W at 4.2 K and 150 W at 77 K. Multiple thermal cycles were performed, with cooldown durations of up to 45 hours. Final base temperatures of approximately 10.8 K (at the coil edge) and 12.0 K (at the coil center) were achieved, with an axial temperature difference of approximately 1.2 K. The warm-up period extended over approximately 24.6 hours. Voltage measurements from all four layers were recorded during cooldown and warmup. The system demonstrated stable cooldown performance, repeatable gradients, and good thermal anchoring. These results support the feasibility of conduction cooling in large-scale HTS magnets, aligning with broader goals for “green” cryogen-free accelerator technologies and paving the way for more sustainable, scalable, and energy-efficient high-field magnet systems in next-generation particle accelerators.

Canted cosine theta magnets↗
DOE OSTI2020

Novel AC Distribution Factor for Efficient Outage Analysis

We propose novel line outage distribution factors under AC power flow (AC-LODF) for providing very efficient approximate post-outage power flow solution. The AC-LODF is derived by using holomorphic embedding (HE) approach, and the approximate post-outage power flow solution is obtained with Pad approximation (PA). In this work, tests on IEEE 14-bus system, and Polish 2383-bus system verify that the proposed AC-LODF is effective on most outages, and can significantly speed up outage analysis.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Probabilistic Voltage Sensitivity Analysis to Quantify Impact of High PV Penetration on Unbalanced Distribution System

From an operational and planning perspective, it is important to quantify the impact of increasing penetration of photovoltaics on the distribution system. Most existing impact assessment studies are scenario-based where derived results are scenario specific and not generalizable. Moreover, stochasticity in the temporal behavior of spatially distributed PVs requires a large number of scenarios that increase with the size of the network and the level of penetration. Therefore, we propose a new computationally efficient analytical framework of voltage sensitivity analysis that allows for stochastic analysis of voltage change due to random changes in PV generation. We first derive an analytical approximation for voltage change at any node of the network due to change in power at other nodes in an unbalanced distribution network. The quality of this approximation is reinforced via bounds on the approximation error. Then, we derive the probability distribution of voltage change at a certain node due to random changes in power injections/consumptions at multiple locations of the network. The accuracy of the proposed PVSA is illustrated using a modified version of the IEEE 37 bus test system. As a result, the proposed PVSA can serve as a powerful tool for proactive monitoring/control and ease the computational burden associated with perturbation based cybersecurity mechanisms.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Probing Quantum Telecloning on Superconducting Quantum Processors

Quantum information cannot be perfectly cloned, but approximate copies of quantum information can be generated. Quantum telecloning combines approximate quantum cloning, more typically referred to as quantum cloning, and quantum teleportation. Quantum telecloning allows approximate copies of quantum information to be constructed by separate parties, using the classical results of a Bell measurement made on a prepared quantum telecloning state. Quantum telecloning can be implemented as a circuit on quantum computers using a classical coprocessor to compute classical feedforward instructions using if statements based on the results of a midcircuit Bell measurement in real time. We present universal symmetric optimal 1 → M telecloning circuits and experimentally demonstrate these quantum telecloning circuits for M = 2 up to M = 10, natively executed with real-time classical control systems on IBM Quantum superconducting processors, known as dynamic circuits. We perform the cloning procedure on many different message states across the Bloch sphere, on seven IBM Quantum processors, optionally using the error suppression technique X–X sequence digital dynamical decoupling. Two circuit optimizations are utilized: one that removes ancilla qubits for M = 2,3 , and one that reduces the total number of gates in the circuit but still uses ancilla qubits. Parallel single-qubit tomography with maximum likelihood estimation density matrix reconstruction is used in order to compute the mixed-state density matrices of the clone qubits, and clone quality is measured using quantum fidelity. These results present one of the largest and most comprehensive noisy intermediate-scale quantum computer experimental analyses on (single qubit) quantum telecloning to date. The clone fidelity sharply decreases to 0.5 for M > 5 , but for M = 2 , we are able to achieve a mean clone fidelity of up to 0.79 using dynamical decoupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Advanced Multilevel Monte Carlo Methods

This article reviews the application of some advanced Monte Carlo techniques in the context of multilevel Monte Carlo (MLMC). MLMC is a strategy employed to compute expectations, which can be biassed in some sense, for instance, by using the discretization of an associated probability law. The MLMC approach works with a hierarchy of biassed approximations, which become progressively more accurate and more expensive. Additionally, using a telescoping representation of the most accurate approximation, the method is able to reduce the computational cost for a given level of error versus i.i.d. sampling from this latter approximation. All of these ideas originated for cases where exact sampling from couples in the hierarchy is possible. This article considers the case where such exact sampling is not currently possible. We consider some Markov chain Monte Carlo and sequential Monte Carlo methods, which have been introduced in the literature, and we describe different strategies that facilitate the application of MLMC within these methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sorghum Association Panel whole‐genome sequencing establishes cornerstone resource for dissecting genomic diversity

SUMMARY Association mapping panels represent foundational resources for understanding the genetic basis of phenotypic diversity and serve to advance plant breeding by exploring genetic variation across diverse accessions. We report the whole‐genome sequencing (WGS) of 400 sorghum ( Sorghum bicolor (L.) Moench) accessions from the Sorghum Association Panel (SAP) at an average coverage of 38× (25–72×), enabling the development of a high‐density genomic marker set of 43 983 694 variants including single‐nucleotide polymorphisms (approximately 38 million), insertions/deletions (indels) (approximately 5 million), and copy number variants (CNVs) (approximately 170 000). We observe slightly more deletions among indels and a much higher prevalence of deletions among CNVs compared to insertions. This new marker set enabled the identification of several novel putative genomic associations for plant height and tannin content, which were not identified when using previous lower‐density marker sets. WGS identified and scored variants in 5‐kb bins where available genotyping‐by‐sequencing (GBS) data captured no variants, with half of all bins in the genome falling into this category. The predictive ability of genomic best unbiased linear predictor (GBLUP) models was increased by an average of 30% by using WGS markers rather than GBS markers. We identified 18 selection peaks across subpopulations that formed due to evolutionary divergence during domestication, and we found six F st peaks resulting from comparisons between converted lines and breeding lines within the SAP that were distinct from the peaks associated with historic selection. This population has served and continues to serve as a significant public resource for sorghum research and demonstrates the value of improving upon existing genomic resources.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Information theory and machine learning illuminate large‐scale metabolomic responses of Brachypodium distachyon to environmental change

SUMMARY Plant responses to environmental change are mediated via changes in cellular metabolomes. However, <5% of signals obtained from liquid chromatography tandem mass spectrometry (LC‐MS/MS) can be identified, limiting our understanding of how metabolomes change under biotic/abiotic stress. To address this challenge, we performed untargeted LC‐MS/MS of leaves, roots, and other organs of Brachypodium distachyon (Poaceae) under 17 organ–condition combinations, including copper deficiency, heat stress, low phosphate, and arbuscular mycorrhizal symbiosis. We found that both leaf and root metabolomes were significantly affected by the growth medium. Leaf metabolomes were more diverse than root metabolomes, but the latter were more specialized and more responsive to environmental change. We found that 1 week of copper deficiency shielded the root, but not the leaf metabolome, from perturbation due to heat stress. Machine learning (ML)‐based analysis annotated approximately 81% of the fragmented peaks versus approximately 6% using spectral matches alone. We performed one of the most extensive validations of ML‐based peak annotations in plants using thousands of authentic standards, and analyzed approximately 37% of the annotated peaks based on these assessments. Analyzing responsiveness of each predicted metabolite class to environmental change revealed significant perturbations of glycerophospholipids, sphingolipids, and flavonoids. Co‐accumulation analysis further identified condition‐specific biomarkers. To make these results accessible, we developed a visualization platform on the Bio‐Analytic Resource for Plant Biology website ( https://bar.utoronto.ca/efp_brachypodium_metabolites/cgi‐bin/efpWeb.cgi ), where perturbed metabolite classes can be readily visualized. Overall, our study illustrates how emerging chemoinformatic methods can be applied to reveal novel insights into the dynamic plant metabolome and stress adaptation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Mesoscale Effects of Composition and Calendering in Lithium-Ion Battery Composite Electrodes

Macrohomogeneous battery models are widely used to predict battery performance, necessarily relying on effective electrode properties, such as specific surface area, tortuosity, and electrical conductivity. While these properties are typically estimated using ideal effective medium theories, in practice they exhibit highly non-ideal behaviors arising from their complex mesostructures. In this paper, we computationally reconstruct electrodes from X-ray computed tomography of 16 nickel–manganese–cobalt-oxide electrodes, manufactured using various material recipes and calendering pressures. Due to imaging limitations, a synthetic conductive binder domain (CBD) consisting of binder and conductive carbon is added to the reconstructions using a binder bridge algorithm. Reconstructed particle surface areas are significantly smaller than standard approximations predicted, as the majority of the particle surface area is covered by CBD, affecting electrochemical reaction availability. Finite element effective property simulations are performed on 320 large electrode subdomains to analyze trends and heterogeneity across the electrodes. Significant anisotropy of up to 27% in tortuosity and 47% in effective conductivity is observed. Electrical conductivity increases up to 7.5× with particle lithiation. We compare the results to traditional Bruggeman approximations and offer improved alternatives for use in cell-scale modeling, with Bruggeman exponents ranging from 1.62 to 1.72 rather than the theoretical value of 1.5. We also conclude that the CBD phase alone, rather than the entire solid phase, should be used to estimate effective electronic conductivity. This study provides insight into mesoscale transport phenomena and results in improved effective property approximations founded on realistic, image-based morphologies.

25 ENERGY STORAGE↗
DOE OSTI2022

Designing Cladded Components for High Temperature Nuclear Service: Part I—Analysis Methods

Use of corrosion-resistant cladding can greatly extend the design life of structural components in many advanced reactor systems. However, there are currently no ASME design rules for cladded components to guard against creep-fatigue failure and ratcheting strain accumulation in elevated temperature nuclear service. This paper, presented in two parts, addresses this gap by proposing a design strategy for cladded components that do not require long-term testing of clad materials. The proposed approach relies on approximate design analysis methods for two types of clad materials—a soft clad that creeps much faster than and has lower yield stress than the class A base material and a hard clad that creeps much slower than and has higher yield stress than the class A base material. Part I discusses the approximate analysis strategies for the clad materials—treat a soft clad as perfectly compliant and a hard clad as linear elastic—and Part II develops a complete set of design rules for each of the two types of cladded components. In conclusion, finite element analyses of representative high temperature reactor components show that the proposed design analysis methods can bound the design quantities in soft cladded components and approximate the design quantities in hard cladded components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Designing Cladded Components for High Temperature Nuclear Service: Part II—Design Rules

The challenge of using existing ASME Section III, Division 5, class A metallic materials for the construction of structural components of advanced reactors with corrosive coolants could be mitigated by allowing designers to use cladding to protect the base material from corrosion. However, the existing Section III, Division 5 rules provide no guidance on the evaluation of strain accumulation and creep-fatigue damage in cladded components. The availability of design rules for cladded components that do not require long-term clad material testing could promote the application of the cladding approach to accelerate the deployment schedule of these advanced reactor systems. To avoid long-term properties for the clad materials Part I of this work proposes two approximate design analysis methods for two types of clad materials—soft clads that creep much faster and have lower yield stress than the class A base material, and hard clads that creep much slower and have higher yield stress than the class A base material. The proposed analysis methods approximate the response of a soft clad by treating it as perfectly compliant and of a hard clad by treating it as linear elastic. Based on these approximate design analysis strategies this Part II develops a complete set of design rules for class A components cladded with either soft or hard clad materials. In conclusion, Part II discusses the reasoning behind the proposed design rules and uses example finite element analyses of representative reactor components to illustrate the use of these design methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Reliability of Copper Inverse Opal Surfaces for Extreme-Heat-Flux Micro-Coolers in Low-Global-Warming-Potential Refrigerant R-1233zd Pool Boiling Experiments

This presentation provides a brief snapshot of the InterPACK paper InterPACK2023-113781. The paper explores copper inverse opal (CIO) surface reliability in pool boiling experiments in water and a new, low-global-warming-potential (GWP = 1) hydrofluoroolefin (HFO) refrigerant R-1233zd. The CIO-based structure is intended to develop enhanced two-phase heat transfer surfaces for extreme-heat-flux (approximately 1 kW/cm2) micro-coolers. In this study, a limited number of pool boiling experiments were performed using water and HFO-1233zd fluid, and the reliability of the CIO-based surfaces was evaluated. Critical heat flux (CHF) values in HFO-1233zd at 40 Degrees Celsius to 45 Degrees Celsius saturation temperatures and the corresponding saturation pressures were also measured. The CHF values with the refrigerant are significantly lower compared to those with water, but the refrigerant allows for a wider usable temperature range in the end application of the micro-coolers and is not limited to data centers with controlled ambient conditions. Reliability experiments with CIO surface samples - involving pool boiling with water on the CIO surfaces for approximately 48 hours and with HFO-1233zd for 144 hours - showed no structural degradation of the enhanced surface or any significant performance drop in heat transfer coefficients. The CIO surface samples in water were oxidized, most likely due to the presence of air in water and in the experimental vessel.

critical heat flux↗
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 OSTI2023

Generalized analytical and numerical modeling of optical second harmonic generation in anisotropic crystals and complex heterostructures using #SHAARP package

Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. The closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, here we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP) for single interface (si) and in multilayers (ml) for homogeneous crystals. The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials. SHAARP.si and SHAARP.ml are available through GitHub https://github.com/Rui-Zu/SHAARP and https://github.com/bzw133/SHAARP.ml, respectively.

complex systems↗
DOE OSTI2020

Physics-Informed Generative Adversarial Networks for Stochastic Differential Equations

We developed a new class of physics-informed generative adversarial networks (PI-GANs) to solve forward, inverse, and mixed stochastic problems in a unified manner based on a limited number of scattered measurements. Unlike standard GANs relying solely on data for training, here we encoded into the architecture of GANs the governing physical laws in the form of stochastic differential equations (SDEs) using automatic differentiation. In particular, we applied Wasserstein GANs with gradient penalty (WGAN-GP) for its enhanced stability compared to vanilla GANs. We first tested WGAN-GP in approximating Gaussian processes of different correlation lengths based on data realizations collected from simultaneous reads at sparsely placed sensors. We obtained good approximation of the generated stochastic processes to the target ones even if there is a mismatch between the input noise dimensionality and the effective dimensionality of the target stochastic processes. We also studied the overfitting issue for both the discriminator and the generator, and we found that overfitting occurs also in the generator in addition to the discriminator as previously reported. Subsequently, we considered the solution of elliptic SDEs requiring approximations of three stochastic processes, namely the solution, the forcing, and the diffusion coefficient. Here again, we assumed data collected from simultaneous reads at a limited number of sensors for the multiple stochastic processes. Furthermore, three generators were used for the PI-GANs: two of them were feed forward deep neural networks (DNNs), while the other one was the neural network induced by the SDE. For the case where we have one group of data, we employed one feed forward DNN as the discriminator, while for the case of multiple groups of data we employed multiple discriminators in PI-GANs. We solved forward, inverse, and mixed problems without changing the framework of PI-GANs, obtaining both the means and the standard deviations of the stochastic solution and the diffusion coefficient in good agreement with benchmarks. In this work, we have demonstrated the effectiveness of PI-GANs in solving SDEs for about 120 dimensions. In principle, PI-GANs could tackle very high dimensional problems given more sensor data with low-polynomial growth in computational cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Estimating Higher-Order Moments Using Symmetric Tensor Decomposition

In this paper, we consider the problem of decomposing higher-order moment tensors, i.e., the sum of symmetric outer products of data vectors. Such a decomposition can be used to estimate the means in a Gaussian mixture model and for other applications in machine learning. The dth-order empirical moment tensor of a set of p observations of n variables is a symmetric d-way tensor. Our goal is to nd a low-rank tensor approximation comprising r $\ll$ p symmetric outer products. The challenge is that forming the empirical moment tensor costs O(pn d ) operations and O(n d ) storage, which may be prohibitively expensive; additionally, the algorithm to compute the low-rank approximation costs O(n d ) per iteration. Our contribution is avoiding formation of the moment tensor, computing the low-rank tensor approximation of the moment tensor implicitly using O(pnr) operations per iteration and no extra memory. This advance opens the door to more applications of higher-order moments since they can now be efficiently computed. We present numerical evidence of the computational savings and show an example of estimating the means for higher-order moments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Modifying the Asynchronous Jacobi Method for Data Corruption Resilience

Moving scientific computation from high-performance computing (HPC) and cloud computing (CC) environments to devices on the edge, i.e., physically near instruments of interest, has received tremendous interest in recent years. Such edge computing environments can operate on data in situ, offering enticing benefits over data aggregation to HPC and CC facilities that include avoiding costs of transmission, increased data privacy, and real-time data analysis. Because of the inherent unreliability of edge computing environments, new fault-tolerant approaches must be developed before the benefits of edge computing can be realized. Motivated by algorithm-based fault tolerance, a variant of the asynchronous Jacobi (ASJ) method is developed that achieves resilience to data corruption by rejecting solution approximations from neighbor devices according to a bound derived from convergence theory. Numerical results on a two-dimensional Poisson problem show that the new rejection criterion, along with a novel approximation to the shortest path length on which the criterion depends, restores convergence for the ASJ variant in the presence of certain types data corruption. Numerical results are obtained for when the singular values in the analytic bound are approximated. Additional linear systems are also explored, one with a more dense sparsity pattern and one that includes advection. All results indicate that successful resilience to data corruption depends on whether the bound tightens fast enough to reject corrupted data before the iteration evolution deviates significantly from that predicted by the convergence theory defining the bound. This observation generalizes to future work on algorithm-based fault tolerance for other asynchronous algorithms, including upcoming approaches that leverage Krylov subspaces.

97 MATHEMATICS AND COMPUTING↗