Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model reduction”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Anomalous Sound Velocities of Antigorite at High Pressure and Implications for Detecting Serpentinization at Mantle Wedges

Elastic properties of antigorite are important for interpretation of seismic mapping of serpentinization in the mantle wedge above subducting slabs. The compressional (V P ) and shear (V S ) wave velocities in pure antigorite aggregates were measured simultaneously up to 8.4 GPa by ultrasonic interferometry. We found that V P increases monotonically with pressure while V S increases with pressure up to about 3 GPa but undergoes a negative pressure dependence above 4 GPa. Compared to other mantle minerals, antigorite exhibits significantly lower P and S wave velocities as well as a higher V P /V S ratio at upper mantle pressures. We modeled velocity reductions manifesting through the formation of antigorite in mantle peridotite and provide compelling evidence that seismic anomalies with low-velocity and high V P /V S ratios is caused by varying degrees of serpentinization in subduction zones.

58 GEOSCIENCES↗
DOE OSTI2018

Sterically controlled mechanochemistry under hydrostatic pressure

Mechanical stimuli can modify the energy landscape of chemical reactions and enable reaction pathways, offering a synthetic strategy that complements conventional chemistry. These mechanochemical mechanisms have been studied extensively in one-dimensional polymers under tensile stress using ring-opening and reorganization, polymer unzipping and disulfide reduction as model reactions. In these systems, the pulling force stretches chemical bonds, initiating the reaction. Additionally, it has been shown that forces orthogonal to the chemical bonds can alter the rate of bond dissociation. Furthermore, these bond activation mechanisms have not been possible under isotropic, compressive stress (that is, hydrostatic pressure). Here we show that mechanochemistry through isotropic compression is possible by molecularly engineering structures that can translate macroscopic isotropic stress into molecular-level anisotropic strain.

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

Forest residue harvest optimization: spanning the bridge between plant biology and biorefinery performance

Forestry residues have immense potential as alternative feedstocks to petroleum, yet their inherent complexity remains a major challenge to widespread use. Pairing the temporal rhythms of plant biology with biorefinery performance is critical to industrial-scale biorefinery development. Here, we provide the first report of a techno-economic analysis (TEA) and life cycle assessment (LCA) for a model integrated reductive catalytic fractionation (RCF)–molten salt hydrolysis process for forestry residues varying in tree part, species, and phenophase. All forestry residues resulted in net-negative greenhouse gas (GHG) emissions vs. comparable petroleum feedstocks, with GHG emissions potentially reduced >4.0× through composition-based feedstock selection (e.g., harvesting American beech bark in spring vs. summer). Moreover, American beech twigs/branchlets and bark in leafed and emergence phenophases, respectively, had 7.9× lower predicted phenolic minimum selling prices (MSPs) vs. other feedstocks and MSPs within the current global phenolic market range. Hemicellulose content and RCF yield emerged as key parameters impacting GHG emissions and biorefinery revenue, identifying hardwood twigs/branchlets in the leafed phenophase as optimal biofeedstocks. Biorefinery expenses were dominated by purchased equipment, raw materials, and utility costs, highlighting essential areas for future study. Notably, RCF reactor pressures drove 85–90% of equipment costs, but sensitivity analysis revealed that decreasing the pressure 20% could reduce the phenol MSP 4-fold. Structural carbohydrate dynamics were also investigated using a two-step acid hydrolysis method to resolve tissue- and species-level patterns in biomass composition throughout the year to enable harvest optimization based on TEA/LCA findings. Ultimately, elucidating the impact of biofeedstock dynamics on biorefinery performance enables harvest optimization, informed engineering design, and progress towards an expanded bioeconomy.

Shapiro, Alison J. [University of Delaware, Newark↗
DOE OSTI2018

A posteriori snapshot location for POD in optimal control of linear parabolic equations

In this paper we study the approximation of an optimal control problem for linear parabolic PDEs with model order reduction based on Proper Orthogonal Decomposition (POD-MOR). POD-MOR is a Galerkin approach where the basis functions are obtained upon information contained in time snapshots of the parabolic PDE related to given input data. In the present work we show that for POD-MOR in optimal control of parabolic equations it is important to have knowledge about the controlled system at the right time instances. We propose to determine the time instances (snapshot locations) by an a posteriori error control concept. The proposed method is based on a reformulation of the optimality system of the underlying optimal control problem as a second order in time and fourth order in space elliptic system which is approximated by a space-time finite element method. Finally, we present numerical tests to illustrate our approach and show the effectiveness of the method in comparison to existing approaches.

Mathematics↗
DOE OSTI2016

Data-driven parameterization of the generalized Langevin equation

Significance The generalized Langevin equation (GLE) provides a precise description of coarse-grained variable dynamics in reduced dimension models. However, computation of the memory kernel poses a major challenge to the practical use of the GLE. This paper presents a data-driven approach to compute the memory kernel, relying on a hierarchy of parameterized rational approximations in terms of the Laplace transform, which can be expanded to arbitrarily high order as necessary. This parameterization makes it convenient to represent the GLE via an extended stochastic model where the memory term is eliminated by properly introducing auxiliary variables. The present method is well-suited for constructing reduced models for nonequilibrium properties of complex systems such as biomolecules, chemical reaction networks, and climate simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗
DOE OSTI2022

Improving Schrödinger Equation Implementations with Gray Code for Adiabatic Quantum Computers

We reformulate the continuous-space Schrödinger equation in terms of spin Hamiltonians. For the kinetic energy operator, the critical concept facilitating the reduction in model complexity is the idea of position encoding. A binary encoding of position produces a spin-1/2 Heisenberg-like model and yields exponential improvement in space complexity when compared to classical computing. Encoding with a binary reflected Gray code (BRGC), and a Hamming-distance-2 Gray code (H2GC) reduces the model complexity down to the 𝑋⁢𝑍 and transverse Ising model, respectively. For 𝐴 qubits BRGC yields 2 𝐴 positions and is reduced to its 2-local form with O⁡(𝐴) ancillary qubits. H2GC yields 2 𝐴/2+1 positions with O⁡(𝐴 2 ) three-local penalty terms. We also identify the bijective mapping between diagonal unitaries and the Walsh series, producing the mapping of any real potential to a series of 𝑘 -local Ising models through the fast Walsh transform. Finally, in a finite volume, we provide some numerical evidence to support the claim that the total time needed for adiabatic evolution is protected by the infrared cutoff of the system. As a result, initial state preparation from a free-field wave function to an interacting system is expected to exhibit polynomial time complexity with volume and constant scaling with respect to lattice discretization for all encodings. For H2GC, if the evolution starts with the transverse Hamiltonian due to hardware restrictions, then penalties are dynamically introduced such that the low-lying spectrum reproduces the energy levels of the Laplacian. The adiabatic evolution of the penalty Hamiltonian is therefore sensitive to the ultraviolet scale. It is expected to exhibit polynomial time complexity with lattice discretization, or exponential time complexity with respect to the number of qubits given a fixed volume.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Event Detection and Classification Using Machine Learning Applied to PMU Data for the Western US Power System

Smart grid technology enhances our comprehension and reliability of the power grid, leveraging Phasor Measurement Unit (PMU) data—time-synchronized, high-frequency measurements gathered across the US power grid. This paper employs machine learning techniques to effectively analyze the vast PMU data in Wide Area Monitoring Systems (WAMS) for power grid event detection and classification. Analyzing several months of real-world PMU data, the paper focuses on machine learning for fast, precise event detection and classification, corroborated by utility event logs. Practical challenges like feature extraction, dimensionality reduction, and model selection are addressed. A novel feature yielding improved results is discovered, and a supplementary algorithm for detecting small power grid faults is developed. The final algorithm is validated using a month-long real PMU data set, demonstrating its capability in accurately identifying power grid events in near real-time.

machine learning, event detection, PMU↗
DOE OSTI2018

Reduced-Order Aggregate Model for Parallel-Connected Single-Phase Inverters

This paper outlines a reduced-order aggregate dynamical model for parallel-connected single-phase grid-connected inverters. For each inverter, we place no restrictions on the converter topology and merely assume that the ac-side switch-averaged voltage can be controlled via pulsewidth modulation. The ac output of each inverter interfaces through an LCL filter to the grid. The closed-loop system contains a phase locked loop for grid synchronization, and real- and reactive-power control are realized with inner and outer PI current- and power-control loops. We derive a necessary and sufficient set of parametric relationships to ensure that a reduced-order aggregated state-space model for an arbitrary number of such paralleled inverters has the same model order and structure as any single inverter. We also present reduced-order models for the settings where the real- and reactive-power setpoints are different and where the inverters have different power ratings. We anticipate the proposed model being useful in analyzing the dynamics of large collections of parallel-connected inverters with minimal computational complexity. The aggregate model is validated against measurements obtained from a multi-inverter experimental setup consisting of three 750-VA paralleled grid-connected inverters, hence establishing robustness of the analytical result to parametric variations seen in practice.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Optimizing DER Participation in Inertial and Primary-Frequency Response

This paper develops an approach to enable the optimal participation of distributed energy resources (DERs) in inertial and primary-frequency response alongside conventional synchronous generators. Leveraging a reduced-order model description of frequency dynamics, DERs' synthetic inertias and droop coefficients are designed to meet time-domain performance objectives of frequency overshoot and steady-state regulation. Furthermore, an optimization-based method centered around classical economic dispatch is developed to ensure that DERs share the power injections for inertial- and primary-frequency response in proportion to their power ratings. Simulations for a modified New England test-case system composed of ten synchronous generators and six instances of the IEEE 37-node test feeder with frequency-responsive DERs validate the design strategy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

TU-AB-BRC-08: Egs-brachy, a Fast and Versatile Monte Carlo Code for Brachytherapy Applications

Purpose: To introduce egs-brachy, a new, fast, and versatile Monte Carlo code for brachytherapy applications. Methods: egs-brachy is an EGSnrc user-code based on the EGSnrc C++ class library (egs++). Complex phantom, applicator, and source model geometries are built using the egs++ geometry module. egs-brachy uses a tracklength estimator to score collision kerma in voxels. Interaction, spectrum, energy fluence, and phase space scoring are also implemented. Phase space sources and particle recycling may be used to improve simulation efficiency. HDR treatments (e.g. stepping source through dwell positions) can be simulated. Standard brachytherapy seeds, as well as electron and miniature x-ray tube sources are fully modelled. Variance reduction techniques for electron source simulations are implemented (Bremsstrahlung cross section enhancement, uniform Bremsstrahlung splitting, and Russian Roulette). TG-43 parameters of seeds are computed and compared to published values. Example simulations of various treatments are carried out on a single 2.5 GHz Intel Xeon E5-2680 v3 processor core. Results: TG-43 parameters calculated with egs-brachy show excellent agreement with published values. Using a phase space source, 2% average statistical uncertainty in the PTV ((2mm){sup 3} voxels) can be achieved in 10 s for 100 {sup 125}I or {sup 103}Pd seeds in a 36.2 cm{sup 3} prostate PTV, 31 s for 64 {sup 103}Pd seeds in a 64 cm{sup 3} breast PTV, and 56 s for a miniature x-ray tube in a 27 cm{sup 3} breast PTV. Comparable uncertainty is reached in 12 s in a (1 mm){sup 3} water voxel 5 mm away from a COMS 16mm eye plaque with 13 {sup 103}Pd seeds. Conclusion: The accuracy of egs-brachy has been demonstrated through benchmarking calculations. Calculation times are sufficiently fast to allow full MC simulations for routine treatment planning for diverse brachytherapy treatments (LDR, HDR, miniature x-ray tube). egs-brachy will be available as free and open-source software to the medical physics research community. This work is partially funded by the Canada Research Chairs program, the Natural Sciences and Engineering Research Council of Canada, and the Ontario Ministry of Research and Innovation (Ontario Early Researcher Award).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the influence of conformational uncertainty in biomolecular solvation

Biomolecules exhibit conformational fluctuations near equilibrium states, inducing uncertainty in various biological properties in a dynamic way. We have developed a general method to quantify the uncertainty of target properties induced by conformational fluctuations. Using a generalized polynomial chaos (gPC) expansion, we construct a surrogate model of the target property with respect to varying conformational states. We also propose a method to increase the sparsity of the gPC expansion by defining a set of conformational “active space” random variables. With the increased sparsity, we employ the compressive sensing method to accurately construct the surrogate model. We demonstrate the performance of the surrogate model by evaluating fluctuation-induced uncertainty in solvent-accessible surface area for the bovine trypsin inhibitor protein system and show that the new approach offers more accurate statistical information than standard Monte Carlo approaches. Further more, the constructed surrogate model also enables us to directly evaluate the target property under various conformational states, yielding a more accurate response surface than standard sparse grid collocation methods. In particular, the new method provides higher accuracy in high-dimensional systems, such as biomolecules, where sparse grid performance is limited by the accuracy of the computed quantity of interest. Our new framework is generalizable and can be used to investigate the uncertainty of a wide variety of target properties in biomolecular systems.

coarse graining↗
DOE OSTI2017

CHROTRAN, 1.0

CHROTRAN is a fork of the widely-used PFLOTRAN flow and reactive transport numerical simulation code. It implements custom physics and chemistry appropriate to the design of in-situ reduction of heavy metals such as Cr(VI) in groundwater. CHROTRAN includes full dynamics for five species: the metal to be remediated, an electron donor, biofilm, a nontoxic conservative bio-inhibitor, and a biocide. Direct abiotic reduction by donor-metal interaction as well as donor-driven biomass growth and bio-reduction are modeled, along with crucial processes such as donor sorption, and biofilm inactivation. The software implementation handles heterogeneous flow fields, arbitrarily many chemical species and amendment injection points, and features full coupling between flow and reactive transport, allowing for assessment of the effect of bio-fouling.

Hansen, ScottK.↗
DOE OSTI2018

Ex Situ Catalytic Fast Pyrolysis of Lignocellulosic Biomass to Hydrocarbon Fuels: 2018 State of Technology and Future Research

This report will update modeled costs for a fixed bed ex situ catalytic fast pyrolysis process for the conversion of biomass to liquid hydrocarbon fuels based on recent bench-scale experimental results. It will also update previously published information for this conversion pathway in the DOE Bioenergy Technologies Office 2016 Multi-Year Program Plan with current understanding of research goals achievable by 2022 and associated modeled cost reductions. All costs will be updated to a consistent set of assumptions in 2016 dollars.

09 BIOMASS FUELS↗
DOE OSTI2018

High-Octane Gasoline from Lignocellulosic Biomass via Syngas and Methanol/Dimethyl Ether Intermediates: 2018 State of Technology and Future Research

This report will update modeled costs for the production of high-octane gasoline via gasification of biomass followed by syngas cleanup and conversion to methanol/dimethyl ether. The 2018 State of Technology will be based on recent bench-scale experimental results. This report will also update previously published information for this conversion pathway in the DOE Bioenergy Technologies Office 2016 Multi-Year Program Plan with current understanding of research goals achievable by 2022 and associated modeled cost reductions. All costs will be updated to a consistent set of assumptions in 2016 dollars.

09 BIOMASS FUELS↗
DOE OSTI2020

Parameterized Neural Ordinary Differential Equations: Applications to Computational Physics Problems

This work proposes an extension of neural ordinary differential equations (NODEs) by introducing an additional set of ODE input parameters to NODEs. This extension allows NODEs to learn multiple dynamics specified by the input parameter instances. Our extension is inspired by the concept of parameterized ordinary differential equations, which are widely investigated in computational science and engineering contexts, where characteristics of the governing equations vary over the input parameters. We apply the proposed parameterized NODEs (PNODEs) for learning latent dynamics of complex dynamical processes that arise in computational physics, which is an essential component for enabling rapid numerical simulations for time-critical physics applications. For this, we propose an encoder-decoder-type framework, which models latent dynamics as PNODEs. We demonstrate the effectiveness of PNODEs with important benchmark problems from computational physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗