SEARCH · Engineering Papers

Results for “discretization”

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 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Per-Phase and 3-Phase Optimal Coordination of Directional Overcurrent Relays Using Genetic Algorithm

Penetration of the power grid by renewable energy sources, distributed storage, and distributed generators is becoming increasingly common. Increased utilization of these distributed energy resources (DERs) has given rise to additional protection coordination concerns, particularly where they are utilized in an unbalanced manner or where loading among phases is unbalanced. Digital relays such as the SEL-751 (produced by Schweitzer Engineering Laboratories, Pullman, WA, USA) series have the capability of being set on a per-phase basis. This capability is underutilized in common practice. Additionally, in optimization algorithms for determining relay settings, the time-overcurrent characteristics (TOCs) of relays are generally not treated as variables and are assigned before running the optimization algorithm. In this paper, TOC options themselves are treated as discrete variables to be considered in the optimization algorithm. A mixed integer nonlinear programming problem (MINLP) is set up where the goal is to minimize relay operating times. A genetic algorithm (GA) approach is implemented in MATLAB where two cases are considered. In the first case, the TOC and Time dial setting (TDS) of each relay is set on a three-phase basis. In the second case, per-phase settings are considered. Relay TDSs and TOCs are both considered as simultaneous discrete control variables. Despite the stochastic nature of using per-phase settings for unbalanced systems is found to generally allow for shorter operating times. However, for relatively balanced systems, it is best to use three-phase settings if computation time is of importance.

Matthews, Ronald C.↗

DOE OSTI2021

Contaminant Source Identification from Finite Sensor Data: Perron–Frobenius Operator and Bayesian Inference

Sensors in the built environment ensure safety and comfort by tracking contaminants in the occupied space. In the event of contaminant release, it is important to use the limited sensor data to rapidly and accurately identify the release location of the contaminant. Identification of the release location will enable subsequent remediation as well as evacuation decision-making. In previous work, we used an operator theoretic approach—based on the Perron–Frobenius (PF) operator—to estimate the contaminant concentration distribution in the domain given a finite amount of streaming sensor data. In the current work, the approach is extended to identify the most probable contaminant release location. The release location identification is framed as a Bayesian inference problem. The Bayesian inference approach requires considering multiple release location scenarios, which is done efficiently using the discrete PF operator. The discrete PF operator provides a fast, effective and accurate model for contaminant transport modeling. The utility of our PF-based Bayesian inference methodology is illustrated using single-point release scenarios in both two and three-dimensional cases. The method provides a fast, accurate, and efficient framework for real-time identification of contaminant source location.

42 ENGINEERING↗

DOE OSTI2022

Projecting the Thermal Response in a HTGR-Type System during Conduction Cooldown Using Graph-Laplacian Based Machine Learning

Accurate prediction of an off-normal event in a nuclear reactor is dependent upon the availability of sensory data, reactor core physical condition, and understanding of the underlying phenomenon. This work presents a method to project the data from some discrete sensory locations to the overall reactor domain during conduction cooldown scenarios similar to High Temperature Gas-cooled Reactors (HTGRs). The existing models for conductive cooldown in a heterogeneous multi-body system, such as an assembly of prismatic blocks or pebble beds relies on knowledge of the thermal contact conductance, requiring significant knowledge of local thermal contacts and heat transport possibilities across those contacts. With a priori knowledge of bulk geometry features and some discrete sensors, a machine learning approach was devised. The presented work uses an experimental facility to mimic conduction cooldown with an assembly of 68 cylindrical rods initially heated to 1200 K. High-fidelity temperature data were collected using an infrared (IR) camera to provide training data to the model and validate the predicted temperature data. The machine learning approach used here first converts the macroscopic bulk geometry information into Graph-Laplacian, and then uses the eigenvectors of the Graph-Laplacian to develop Kernel functions. Support vector regression (SVR) was implemented on the obtained Kernels and used to predict the thermal response in a packed rod assembly during a conduction cooldown experiment. The usage of SVR modeling differs from most models today because of its representation of thermal coupling between rods in the core. When trained with thermographic data, the average normalized error is less than 2% over 400 s, during which temperatures of the assembly have dropped by more than 500 K. The rod temperature prediction performance was significantly better for rods in the interior of the assembly compared to those near the exterior, likely due to the model simplification of the surroundings.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2026

Multigroup Thermal Radiation Transport with Tensor Trains

We investigate the application of tensor-train (TT) algorithms to multigroup thermal radiation transport (i.e., photon radiation transport). The TT framework enables simulations at discretizations that might otherwise be computationally infeasible on conventional hardware. We show that solutions to certain multigroup problems possess an intrinsic low-rank structure, which the TT representation leverages effectively. This enables us to solve problems where the discretized solution size exceeds a trillion parameters on a single node. The solver is evaluated on a range of test problems with varying levels of complexity, consistently achieving compression factors greater than 100× and speedups exceeding 2×. We also investigate alternative TT topologies by analyzing the low-rank structure of the merged spatio-spectral core to assess the potential for greater compression. This analysis suggests that compression gains could increase by factors as large as 7. Our results indicate that the low-rank structure of the merged spatio-spectral core captures the spatio-spectral complexity of the solution, largely driven by the opacity structure of the medium. Beyond identifying opportunities for improved compression, this analysis highlights the types of errors that may arise in angle-integrated quantities when exploiting this low-rank structure.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

A DG-IMEX Method for Two-moment Neutrino Transport: Nonlinear Solvers for Neutrino–Matter Coupling

Neutrino-matter interactions play an important role in core-collapse supernova (CCSN) explosions as they contribute to both lepton number and/or four-momentum exchange between neutrinos and matter, and thus act as the agent for neutrino-driven explosions. Due to the multiscale nature of neutrino transport in CCSN simulations, an implicit treatment of neutrino-matter interactions is desired, which requires solutions of coupled nonlinear systems in each step of the time integration scheme. In this paper we design and compare nonlinear iterative solvers for implicit systems with energy coupling neutrino-matter interactions commonly used in CCSN simulations. Specifically, we consider electron neutrinos and antineutrinos, which interact with static matter configurations through the Bruenn 85 opacity set. The implicit systems arise from the discretization of a nonrelativistic two-moment model for neutrino transport, which employs the discontinuous Galerkin (DG) method for phase-space discretization and an implicit-explicit (IMEX) time integration scheme. In the context of this DG-IMEX scheme, we propose two approaches to formulate the nonlinear systems — a coupled approach and a nested approach. For each approach, the resulting systems are solved with Anderson-accelerated fixedpoint iteration and Newton’s method. The performance of these four iterative solvers has been compared on relaxation problems with various degree of collisionality, as well as proto-neutron star deleptonization problems with several matter profiles adopted from spherically symmetric CCSN simulations. Here, numerical results suggest that the nested Anderson-accelerated fixed-point solver is more efficient than other tested solvers for solving implicit nonlinear systems with energy coupling neutrino-matter interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Scientific machine learning for closure models in multiscale problems: A review

Here, closure problems are omnipresent when simulating multiscale systems, where some quantities and processes cannot be fully prescribed despite their effects on the simulation's accuracy. Recently, scientific machine learning approaches have been proposed as a way to tackle the closure problem, combining traditional (physics-based) modeling with data-driven (machine-learned) techniques, typically through enriching differential equations with neural networks. This paper reviews the different reduced model forms, distinguished by the degree to which they include known physics, and the different objectives of a priori and a posteriori learning. The importance of adhering to physical laws (such as symmetries and conservation laws) in choosing the reduced model form and choosing the learning method is discussed. The effect of spatial and temporal discretization and recent trends toward discretization-invariant models are reviewed. In addition, we make the connections between closure problems and several other research disciplines: inverse problems, Mori-Zwanzig theory, and multi-fidelity methods. In conclusion, much progress has been made with scientific machine learning approaches for solving closure problems, but many challenges remain. In particular, the generalizability and interpretability of learned models is a major issue that needs to be addressed further.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Manufacturing control systems and logic for prognosis of defects in composite materials

Presented are manufacturing control systems for composite-material structures, methods for assembling/operating such systems, and transfer molding techniques for predicting and ameliorating void conditions in fiber-reinforced polymer panels. A method for forming a composite-material construction includes receiving a start signal indicating a fiber-based preform is inside a mold cavity, and transmitting a command signal to inject pressurized resin into the mold to induce resin flow within the mold cavity and impregnate the fiber-based preform. An electronic controller receives, from a distributed array of sensors attached to the mold, signals indicative of pressure and/or temperature at discrete locations on an interior face of the mold cavity. The controller determines a measurement deviation between a calibrated baseline value and the pressure and/or temperature values for each of the discrete locations. If any one of the respective measurement deviations exceeds a calibrated threshold, a void signal is generated to flag a detected void condition.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Automated Adversary Emulation for Cyber-Physical Systems via Reinforcement Learning

Adversary emulation is an offensive exercise that provides a comprehensive assessment of a system’s resilience against cyber attacks. However, adversary emulation is typically a manual process, making it costly and hard to deploy in cyber-physical systems (CPS) with complex dynamics, vulnerabilities, and operational uncertainties. In this paper, we develop an automated, domain-aware approach to adversary emulation for CPS. We formulate a Markov Decision Process (MDP) model to determine an optimal attack sequence over a hybrid attack graph with cyber (discrete) and physical (continuous) components and related physical dynamics. We apply model-based and model-free reinforcement learning (RL) methods to solve the discrete-continuous MDP in a tractable fashion. As a baseline, we also develop a greedy attack algorithm and compare it with the RL procedures. We summarize our findings through a numerical study on sensor deception attacks in buildings to compare the performance and solution quality of the proposed algorithms.

Bhattacharya, Arnab↗

DOE OSTI2020

Efficient Distance-based Global Sensitivity Analysis for Terrestrial Ecosystem Modeling

Sensitivity analysis in terrestrial ecosystem modeling is important for understanding controlling processes, guiding model development, and targeting new observations to reduce parameter and prediction uncertainty. Complex and computationally expensive terrestrial ecosystem models (TEM) limit the number of ensemble simulations, requiring sophisticated and efficient methods to analyze sensitivities of multiple model responses to different types of parameter uncertainties. In this study, we propose a distance-based global sensitivity analysis (DGSA) method. DGSA first classifies model response samples into a small set of discrete classes and then calculates the distance between parameter frequency distributions in different classes to measure the parameter sensitivity. The principle is that, if the parameter distribution is the same in each class, then the model response is insensitive to the parameter, while a large difference in the distributions indicates the parameter is influential to the response. Built on this idea, DGSA can be applied to analyze sensitivity of a single and a group of responses to different kinds of parameter uncertainties including continuous, discrete and even stochastic. Besides the main-effect sensitivity from a single parameter, DGSA can also quantify the sensitivity from parameter interactions. Additionally, DGSA is computationally efficient which can use a small number of model evaluations to obtain an accurate and statistically significant result. We applied DGSA to two TEMs, one having eight parameters and three kinds of model responses, and the other having 47 parameters and a long-period response. We demonstrated that DGSA can be used for sensitivity problems with multiple responses and high-dimensional parameters efficiently.

Lu, Dan↗

DOE OSTI2021

Coupling Finite Element and Finite Volume Simulation Within MOOSE

This work focuses on adding capability to couple finite element (FE) variables into finite volume (FV) physics within INL’s Multiphysics Object-Oriented Simulation Environment (MOOSE). This coupling can allow for improvement of multiphysics simulations where one set of physics is best suited for an FE discretization, while another set is best suited for an FV discretization. Electrohydrodynamics, which govern liquid metal reactor concepts and plasma dynamics, is a fitting example where the electromagnetic field equations are solved with FE and the fluid flow is solved with FV. The new FE to FV coupling method can be summarized as taking the element or face average of the FE variable value or gradient and applying that quantity directly in the FV equation objects.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Performance of a Hybrid HVAC-Integrated Thermal Storage Device

Thermal equipment in buildings is a primary contributor to peak loads on the electrical grid. Thermal energy storage is a cost-effective strategy to decouple electric use from thermal loads, thus reducing grid peak costs for building owners. One method for storing thermal energy in a building is to integrate a phase change material (PCM) directly into the heating, ventilation, and air conditioning system. These systems often require additional glycol loops, pumps, valves, and heat exchangers to couple the storage to the cooling system and building space, which increases the complexity and cost. This work will discuss an alternate approach where the storage is added directly into the heat pump evaporator. A detailed two-dimensional finite difference heat transfer model of a PCM-refrigerant-glycol heat exchanger was developed to simulate the performance of this component. The fluid stream was discretized along the flow direction to capture changes in the fluid properties and local heat transfer rates, and the phase change material was discretized in both the x and y directions to capture the movement of the melt front. The model was used to understand the impact of different material and geometric properties on the charge and discharge characteristics of the device. Finally, a Ragone framework analogous to that used for electrochemical batteries was used to maximize the energy density and round-trip efficiency of the device while supplying loads appropriate for space cooling in buildings.

buildings↗

DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗

DOE OSTI2020

The Influence of Processing Parameters on the Transition Zone for Blended Material 3D Printing

The use of Multiple Materials (MM) in Additive Manufacturing (AM) is increasingly important for expanding the range of applications in the manufacturing industry, particularly for large-format processes. Typically, polymer-based AM incorporates MM transitions through discrete interfaces between layers. This arrangement significantly increases the occurrence of layer delamination failures due to decreased bonding between dissimilar polymers. Elimination of discrete material interfaces by continuously transitioning from Material A to B provides a possible solution. Such continuous gradients could be used to create functionally graded structures that take full advantage of AM’s capability to deliberately impart site-specific properties. Cincinnati’s Big Area Additive Manufacturing (BAAM) system at Oak Ridge National Lab has been equipped with a dual-hopper system that enables in-situ material switching specifically intended for functionally graded and MM printing. The resulting material transition exhibits varied behavior based on printing conditions, which can have an impact on part design and resulting mechanical properties. In this work, the transition zone is characterized as a function of the printing screw speed (related to volumetric flow) and the screw geometry.

Brackett, James↗

DOE OSTI2022

Fiber-optic diagnostic system for future accelerator magnets

The next generation high energy physics accelerators will require magnetic fields at ~20 T. HTS coils will be an essential component of future accelerator magnets and several efforts are currently dedicated on designing 20 T HTS- LTS hybrid magnets. Among the existing challenges, there is the lack of a robust quench detection system for hybrid magnet technology. Another big challenge is represented by the high number of training quenches required by Nb3Sn magnets to reach performance level. In this framework it is important to find a tool that allow local real-time monitoring of magnet strain and temperature. In this paper, we propose the use of fiber optics sensors for diagnostic and quench detection in future accelerator superconducting magnets. Discrete and distributed fiber optic sensors have demonstrated to be a promising tool. The goal is to instrument hundreds of accelerator superconducting magnets and to move beyond the proof-of-concept level. Significant developments are still needed. Here, we are going to present the most recent results and discuss the most urgent technical developments in order to make those sensors a robust and reliable diagnostic tool for accelerator superconducting magnets over the next 10 year. We foresee that discrete fiber sensors will be a stable diagnostic probe for superconducting magnets over the next 3 to 5 years. More R&D work will be necessary for distributed fibers. The most urgent needs are the increase of sample rate and sensitivity. Close collaboration with vendors will be necessary to improve mechanical properties and fabrication processes in order to produce hundreds of meters of fiber and instrument a large number of accelerator superconducting magnets. Those R&D efforts will last up to 10 years with a founding level that spans between 5-10 M$.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

Meshfree Methods for Modeling Chemo-Mechanical Cathode Cracking in Li-Ion Batteries

Through repeated charging and discharging, the electrodes of a Li-ion battery experience cyclic loading causing swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries. Chemo-mechanical models are simulated mainly by two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a chemo-mechanical crack propagation model that achieves both intermediate accuracy and computational time between a CDM and CZM. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

cathode cracking↗

DOE OSTI2023

Meshfree Multiphysics Damage Modeling of Li-ion Battery Materials

Through repeated charging and discharging cycles, the electrodes of a Li-ion battery experience significant swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries and material interfaces. Chemo-mechanical damage models are generally simulated by one of two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a coupled chemo-mechanical crack propagation model that enhances accuracy while maintaining high computational efficiency. This model features a fully coupled, iterative electrochemistry solution, which informs the meshfree damage model of impending crack formation. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

damage modeling↗

DOE OSTI2023

THE IMPACT OF VISCOSITY ON MATERIAL TRANSITIONS IN LARGE FORMAT ADDITIVE MANUFACTURING OF POLYMER COMPOSITES

The need to produce complex geometries incompatible with traditional manufacturing techniques has fueled rapid growth in Large-Format Additive Manufacturing (LFAM). Printing of polymer composite materials have generated significant interest, but the production of Multi-Material (MM) structures with location-based material properties continues to be a challenge. Extrusion-based techniques have utilized multiple deposition heads to successfully print MM structures with both stiff and flexible regions, but these techniques often result in discrete material boundaries that concentrate stress and act as failure points. To avoid discrete interfaces, a novel dual-hopper configuration was developed for the Big Area Additive Manufacturing (BAAM) system that creates a blended material region within the structure. The ability to blend and freely switch between stiff polymer composites and flexible polymers enables printing of robust MM structures with site-specific properties. This study characterizes the influence of viscosity on a blended material transition between carbon fiber-filled acrylonitrile butadiene styrene (CF-ABS) and unfilled thermoplastic polyurethane (TPU), which have significantly different viscoelastic behaviors

Brackett, James↗