SEARCH · Engineering Papers

Results for “reduced order methods”

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 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Development of ultra-pure gadolinium sulfate for the Super-Kamiokande gadolinium project

This paper reports the development and detailed properties of about 13 metric tons of gadolinium sulfate octahydrate, Gd 2 (SO 4 ) 3 • 8H 2 O⁠, which has been dissolved into Super-Kamiokande (SK) in the summer of 2020. We evaluate the impact of radioactive impurities in Gd 2 (SO 4 ) 3 • 8H 2 O on diffuse supernova neutrino background searches and solar neutrino observation and confirm the need to reduce radioactive and fluorescent impurities by about three orders of magnitude from commercially available high-purity Gd 2 (SO 4 ) 3 • 8H 2 O⁠. In order to produce ultra-high-purity Gd 2 (SO 4 ) 3 • 8H 2 O, we have developed a method to remove impurities from gadolinium oxide, Gd 2 O 3 , consisting of acid dissolution, solvent extraction, and pH control processes, followed by a high-purity sulfation process. All of the produced ultra-high-purity Gd 2 (SO 4 ) 3 • 8H 2 O is assayed by inductively coupled plasma mass spectrometry and high-purity germanium detectors to evaluate its quality. Because of the long measurement time of high-purity germanium detectors, we have employed several underground laboratories for making parallel measurements including the Laboratorio Subterráneo de Canfranc in Spain, Boulby in the UK, and Kamioka in Japan. In the first half of production, the measured batch purities were found to be consistent with the specifications. However, in the latter half, the Gd 2 (SO 4 ) 3 • 8H 2 O contained one order of magnitude more 228Ra than the budgeted mean contamination. This was correlated with the corresponding characteristics of the raw material Gd 2 O 3 , in which an intrinsically large contamination was present. Based on their modest impact on SK physics, they were nevertheless introduced into the detector. To reduce 228 Ra for the next stage of gadolinium loading to SK, a new process has been successfully established.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗

DOE OSTI2023

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Faster Randomized Dynamical Decoupling

We present a randomized dynamical decoupling (DD) protocol that can substantially improve the performance of any given deterministic DD scheme for suppressing coherent noise by using no more than two additional pulses. Our construction is implemented by probabilistically applying sequences of pulses, which, when combined, effectively eliminate the error terms that scale linearly with the system-environment coupling strength. As a result, we show that a randomized protocol using a few pulses can outperform deterministic DD protocols that require considerably more pulses. Furthermore, we prove that the randomized protocol provides an improvement compared to deterministic DD sequences that aim to reduce the error in the system’s Hilbert space, such as Uhrig DD, which had been previously regarded to be optimal. To rigorously evaluate the performance, we introduce new analytical methods suitable for analyzing higher-order DD protocols that might be of independent interest. Here, we also present numerical simulations confirming the significant advantage of using randomized protocols compared to widely used deterministic protocols.

Quantum algorithms & computation↗

DOE OSTI2023

Trapped-ion quantum simulation of collective neutrino oscillations

It is well known that the neutrino flavor in extreme astrophysical environments changes under the effect of three contributions: the vacuum oscillation, the interaction with the surrounding matter, and the collective oscillations due to interactions between different neutrinos. The latter adds a nonlinear contribution to the equations of motion, making the description of their dynamics complex. In this work we study various strategies to simulate the coherent collective oscillations of a system of N neutrinos in the two-flavor approximation using quantum computation. This was achieved by using a pair-neutrino decomposition designed to account for the fact that the flavor Hamiltonian, in the presence of the neutrino-neutrino term, presents an all-to-all interaction that makes the implementation of the evolution dependent on the qubit topology. We analyze the Trotter error caused by the decomposition demonstrating that the complexity of the implementation of time evolution scales polynomially with the number of neutrinos and that the noise from near-term quantum device simulation can be reduced by optimizing the quantum circuit decomposition and exploiting a full-qubit connectivity. We find that the gate complexity using second order Trotter-Suzuki formulas scales better with system size than with other decomposition methods such as quantum signal processing. In conclusion, we finally present the application and the results of our algorithm on a real quantum device based on trapped-ion qubits.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

An evaluation of multi-fidelity methods for quantifying uncertainty in projections of ice-sheet mass change

Abstract. This study investigated the computational benefits of using multi-fidelity statistical estimation (MFSE) algorithms to quantify uncertainty in the mass change of Humboldt Glacier, Greenland, between 2007 and 2100 using a single climate change scenario. The goal of this study was to determine whether MFSE can use multiple models of varying cost and accuracy to reduce the computational cost of estimating the mean and variance of the projected mass change of a glacier. The problem size and complexity were chosen to reflect the challenges posed by future continental-scale studies while still facilitating a computationally feasible investigation of MFSE methods. When quantifying uncertainty introduced by a high-dimensional parameterization of the basal friction field, MFSE was able to reduce the mean-squared error in the estimates of the statistics by well over an order of magnitude when compared to a single-fidelity approach that only used the highest-fidelity model. This significant reduction in computational cost was achieved despite the low-fidelity models used being incapable of capturing the local features of the ice-flow fields predicted by the high-fidelity model. The MFSE algorithms were able to effectively leverage the high correlation between each model's predictions of mass change, which all responded similarly to perturbations in the model inputs. Consequently, our results suggest that MFSE could be highly useful for reducing the cost of computing continental-scale probabilistic projections of sea-level rise due to ice-sheet mass change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Impacts of Bulk Microphysics Scheme Structural Choices on Simulations of Rain Initiation Through Drop Coalescence

This study examines how different structural choices in bulk microphysics schemes impact the simulation of warm rain initiation. A single liquid category (SLC) approach prognosing up to four moments of a single drop size distribution (DSD) is compared to the traditional two-category, two-moment approach with separate DSDs for cloud and rain (four total prognostic variables). Different methods for calculating tendencies of the prognostic variables from drop collision-coalescence are also tested: a discretized numerical-integration approach, machine learning via neural networks, lookup tables, and traditional power law fits. Relative to simulations using a bin microphysics model, SLC gives smaller error overall than the two-category approach when numerical integration is used to calculate the collision-coalescence tendencies for both. Replacing the numerical integration with a pre-computed lookup table reduces computational cost with little loss of accuracy. However, using fitted power laws with SLC to represent the collision-coalescence tendencies substantially reduces accuracy and leads to an order of magnitude increase in error. It is also demonstrated that with SLC, reasonably accurate solutions are obtained using only three prognostic moments, while a two-moment SLC scheme leads to substantial error. Overall, both the choice of prognostic moments (e.g., SLC vs. two-category) and method to calculate the collision-coalescence tendencies are important to consider for minimizing errors in bulk schemes. SLC with a sufficiently detailed calculation of the collision-coalescence tendencies provides accurate solutions for a reasonable computational cost, providing a viable alternative to the traditional two-category, two-moment approach for bulk microphysics.

320 (cloud physics and chemistry)↗

DOE OSTI2020

Enhanced dynamic equivalent identification method of large-scale power systems using multiple events

The increasing complexity of the interconnected power system makes high-fidelity dynamic simulation models computationally more intensive. To improve computation efficiency, model reduction techniques have been investigated to only preserve the dynamics in a limited area of interest (study area), while deriving equivalent representation for the external area. For this purpose, a measurement-based reduction approach using system identification techniques has been previously proposed. Accordingly, external areas are represented by dynamic equivalent loads using transfer function estimation. In this paper, in order to enhance the accuracy of the reduced model in preserving dynamics of the study area, multiple “grid” events of different types and at different locations are considered for identifying the parameters of the equivalent loads. Additionally, case studies are carried out in the NPCC 140-bus system. Comparisons are made between multiple events training and single event training, highlighting the advantages of the proposed method in providing a better representation of the grid dynamics under different operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Enhanced Dynamic Equivalent Identification Method of Large-Scale Power Systems Using Multiple Event

The increasing complexity of the interconnected power system makes high-fidelity dynamic simulation models computationally more intensive. To improve computation efficiency, model reduction techniques have been investigated to only preserve the dynamics in a limited area of interest (study area), while deriving equivalent representation for the external area. For this purpose, a measurement-based reduction approach using system identification techniques has been previously proposed. In that case, external areas are represented by dynamic equivalent loads using transfer function estimation. In this paper, in order to enhance the accuracy of the reduced model in preserving dynamics of the study area, multiple “grid” events of different types and at different locations are considered for identifying the parameters of the equivalent loads. Case studies are carried out in the NPCC 140-bus system. Comparison are made between multiple events training and single event training, highlighting the advantages of the proposed method in providing a better representation of the grid dynamics under different operating conditions.

Xue, Yaosuo↗

DOE OSTI2022

Physically rigorous reduced-order flow models of fractured subsurface environments without explosive computational cost

Fractured media models comprise discontinuities of multiple lengths (e.g. fracture lengths and apertures, wellbore area) that fall into the relatively insignificant length scales spanning millimeter-scale fractures to centimeter-scale wellbores in comparison to the extensions of the field of interest, and challenge the conventional discretization methods imposing highly-fine meshing and formidably large numerical cost. By utilizing the recent developments in the finite element analysis of electromagnetics that allow to represent material properties on a hierarchical geometry, this project develops computational capabilities to model fluid flow, heat conduction, transport and induced polarization in large-scale geologic environments that possess geometrically-complex fractures and man-made infrastructures without explosive computational cost. The computational efficiency and robustness of this multi-physics modeling tool are demonstrated by considering various highly-realistic complex geologic environments that are common in many energy and national security related engineering problems.

42 ENGINEERING↗

DOE OSTI2022

Effect of bacterial cell addition on Fe(III) reduction and soil organic matter transformation in a farmland soil

The coupled biogeochemical processes of Fe(III) reduction and organic matter transformation profoundly impact terrestrial carbon cycling. However, little is known about the microbial role in soil organic matter (SOM) transformation during Fe(III) bio-reduction. Here we investigated the bio-reduction behavior of a black farmland soil and corresponding SOM transformation under circumneutral and anoxic conditions. A model dissimilatory Fe reducing bacterium, Geobacter sulfurreducens, was added in either live or dead form in order to enhance Fe(III) reduction in soil and to evaluate the accompanying transformation of SOM. The progress of Fe reduction was monitored and SOM transformation was characterized by various spectroscopy methods. Results showed that addition of either dead or live G. sulfurreducens cells increased the Fe(III) reduction rate and extent. Without cell addition, SOM transformation was insignificant within 13 days of incubation, only with some consumption of aliphatic/protein compounds, apparently due to their higher bio-degradability. Addition of dead or live cells resulted in more drastic SOM transformation, but through different mechanisms. With dead cell amendment, cell necromass and debris stimulated the activity of indigenous soil microbes by serving as extra carbon/energy sources. During Fe(III) reduction, the aliphatic/protein compounds were preferentially consumed by indigenous microbial communities, similar to the treatment without cell addition but with a greater extent of consumption. In comparison, addition of live G. sulfurreducens cells stimulated degradation of less bioavailable compounds, including more saturated and higher molecular weight SOM. It is possible that fast depletion of labile SOM by live G. sulfurreducens cells favored utilization of less bioavailable molecules (such as those with more aromatic structures) by the originally dormant species in native microbial community, suggesting an active role of live Fe(III)-reducing bacteria in affecting SOM transformation. In addition, fast assimilation of microbial related carbon, such as aromatic proteins and microbial byproducts, into SOM pools was also observed. Overall, our results suggest that addition of Fe(III)-reducing bacteria to a native soil can enhance Fe(III) reduction and accelerate the turnover of SOM through various mechanisms. The study provides new insights into coupled Fe(III) reduction and SOM transformation, as well as the “priming effect” of SOM using microbial cells as substrate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Mapping the space of quasisymmetric stellarators using optimized near-axis expansion

A method is demonstrated to rapidly calculate the shapes and properties of quasi-axisymmetric and quasi-helically symmetric stellarators. In this approach, optimization is applied to the equations of magnetohydrodynamic equilibrium and quasisymmetry, expanded in the small distance from the magnetic axis, as formulated by Garren & Boozer [ Phys. Fluids B, vol. 3, 1991, p. 2805]. Due to the reduction of the equations by the expansion, the computational cost is significantly reduced, to times of the order of 1 cpu second, enabling wide and high-resolution scans over parameter space. In contrast to traditional stellarator optimization, here, the cost function serves to maximize the volume in which the expansion is accurate. A key term in the cost function is $\| \boldsymbol {\nabla }\boldsymbol B \|$ , the norm of the magnetic field gradient, to maximize scale lengths in the field. Using this method, a database of $5\times 10^5$ optimized configurations is calculated and presented. Quasisymmetric configurations are observed to exist in continuous bands, varying in the ratio of the magnetic axis length to average major radius. Several qualitatively new types of configuration are found, including quasi-helically symmetric fields in which the number of field periods is two or more than six.

Landreman, Matt (ORCID:000000027233577X)↗

DOE OSTI2026

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

DOE OSTI2019

Memory embedded non-intrusive reduced order modeling of non-ergodic flows

Generating a digital twin of any complex system requires modeling and computational approaches that are efficient, accurate, and modular. Traditional reduced order modeling techniques are targeted at only the first two, but the novel nonintrusive approach we present in this study is an attempt at taking all three into account effectively compared to their traditional counterparts. Based on dimensionality reduction using proper orthogonal decomposition (POD), we introduce a long short-term memory neural network architecture together with a principal interval decomposition (PID) framework as an enabler to account for localized modal deformation. As an effective partitioning tool for breaking the Kolmogorov barrier, our PID framework, therefore, can be considered a key element in the accurate reduced order modeling of convective flows. Our applications for convection-dominated systems governed by Burgers, Navier-Stokes, and Boussinesq equations demonstrate that the proposed approach yields significantly more accurate predictions than the POD-Galerkin method and could be a key enabler toward near real-time predictions of unsteady flows.

42 ENGINEERING↗

DOE OSTI2026

Dynamic aperture studies at the Fermilab recycler ring

As part of the Proton Improvement Plan II (PIP-II), Fermilab aims to increase beam intensity delivered to neutrino experiments. In this context, higher intensity injection into the Recycler Ring enhances space charge effects, pushing operations closer to third-order resonances. These resonances reduce the Dynamic Aperture (DA), leading to increased beam loss. This study presents simulations of DA as a function of tune in the Recycler Ring, incorporating chaos indicators such as the Reversibility Error Method (REM) and Frequency Map Analysis (FMA). The effectiveness of existing resonance mitigation strategies is evaluated by quantifying their ability to delay DA degradation. Additionally, the study examines how space charge detuning and DA limitations dictate viable operational tune points for the Recycler Ring.

Gonzalez-Ortiz, Cristhian [Fermilab]↗

DOE OSTI2025

Dynamic Aperture Studies at the Fermilab Recycler Ring

As part of the Proton Improvement Plan II (PIP-II), Fermilab aims to increase beam intensity delivered to neutrino experiments. In this context, higher intensity injection into the Recycler Ring enhances space charge effects, pushing operations closer to third-order resonances. These resonances reduce the Dynamic Aperture (DA), leading to increased beam loss. This study presents simulations of DA as a function of tune in the Recycler Ring, incorporating chaos indicators such as the Reversibility Error Method (REM). The effectiveness of existing resonance mitigation strategies is evaluated by quantifying their ability to delay DA degradation. Additionally, the study examines how space charge detuning and DA limitations dictate viable operational tune points for the Recycler Ring.

Gonzalez-Ortiz, Cristhian Eduardo [Fermilab]↗

DOE OSTI2021

Key Constructs for Reasonable Confidence in System Validation Programs

This paper presents four key constructs developed in preparing IEEE Std P2411TM, Human Factors Engineering Guide for the Validation of System Designs and Integrated Systems Operations at Nuclear Facilities for submission. These constructs reflect judgments that responsible parties make in order to bring validation to closure. Such judgments weigh competing concerns and cost-benefits for stakeholders. Improved processes and methods may support such judgments but cannot alone eliminate uncertainty. Stating and clarifying these constructs aims to reduce uncertainty in the validation process, to better prepare the implementers and reviewers of future validation programs, both in the nuclear industry and beyond. Reasonable Confidence – A proof standard, as used to reach conclusions in a legal case. Legal proof standards are bounded between “preponderance of evidence” and “beyond reasonable doubt”. Reasonable confidence is comparable to an intermediate legal proof standard of “clear and convincing evidence.” Representative Test Set – A collection of scenarios of sufficient variety to represent the anticipated range of operational conditions, events, evolutions, and activities for validating the system(s) under test. Dispositive vs. Diagnostic Criteria – Relevant performance criteria are placed in one of two categories. Dispositive criteria are assessed, without exception, to determine whether a validation test passes or fails overall. Diagnostic criteria are assessed, subject to justified exception, to evaluate the quality or degree of some particular aspect of performance. Repeatability – Consistency in achieving passing results is to be demonstrated for each scenario in the representative test set; thus, the minimum number of formal repetitions for each scenario is two.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

On the characteristics of the wake of a wind turbine undergoing large motions caused by a floating structure: an insight based on experiments and multi-fidelity simulations from the OC6 project Phase III

This study reports the results of the second round of analyses of the Offshore Code Comparison, Collaboration, Continued, with Correlation and unCertainty (OC6) project Phase III. While the first round investigated rotor aerodynamic loading, here, focus is given to the wake behavior of a floating wind turbine under large motion. Wind tunnel experimental data from the UNsteady Aerodynamics for FLOating Wind (UNAFLOW) project are compared with the results of simulations provided by participants with methods and codes of different levels of fidelity. The effect of platform motion on both the near and the far wake is investigated. More specifically, the behavior of tip vortices in the near wake is evaluated through multiple metrics, such as streamwise position, core radius, convection velocity, and circulation. Additionally, the onset of velocity oscillations in the far wake is analyzed because this can have a negative effect on stability and loading of downstream rotors. Results in the near wake for unsteady cases confirm that simulations and experiments tend to diverge from the expected linearized quasi-steady behavior when the rotor reduced frequency increases over 0.5. Additionally, differences across the simulations become significant, suggesting that further efforts are required to tune the currently available methodologies in order to correctly evaluate the aerodynamic response of a floating wind turbine in unsteady conditions. Regarding the far wake, it is seen that, in some conditions, numerical methods overpredict the impact of platform motion on the velocity fluctuations. Moreover, results suggest that the effect of platform motion on the far wake, differently from original expectations about a faster wake recovery in a floating wind turbine, seems to be limited or even oriented to the generation of a wake less prone to dissipation.

17 WIND ENERGY↗