Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Emulator”

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 703 records · Page 39

Thermally induced cracking on the massive concrete structure of the NSLS II synchrotron and its engineering remediation

Synchrotron accelerator facilities such as the NSLS II require extreme stability, both transient (short-term) and quasi-static (long-term) to achieve the desired resolution performance. Consequently, even μm-level movements, particularly differential movements between locations in the concrete structure supporting the accelerator electron beam lattice (storage ring) or high sensitivity experiments (experimental floor) will lead to serious degradation of its performance. Differential settlement in the overall structure or structural movement exceeding anticipated levels will inevitably degrade the performance and will require intervention. Presented in this paper are the design philosophy of the NSLS II ring structure favoring a monolithic ring, the observed cracking behavior of the young NSLS II concrete following casting and in combination with extreme ambient temperature fluctuation, the results of a non-linear, high-fidelity numerical analysis used to emulate the observed cracking and establish the driving mechanism, the numerical analysis-based identification of the crack-arresting solution and finally the implementation of the remediation solution and the long-term performance of the adopted engineering solution. The multi-stage process revealed that computational methods such as non-linear finite element methods have the potential of providing engineering guidance even when complex structures and in combination with non-linear materials, such as steel reinforcement and concrete are involved.

42 ENGINEERING↗

Bayesian model updating with finite element vs surrogate models: Application to a miter gate structural system

Bayesian finite element (FE) model updating using direct model evaluations of large-scale high-fidelity FE models is extremely computationally expensive. Surrogate models can be used as fast emulators of FE models to accelerate the model calibration process. The physics/mechanics-based FE models are still the underpinning behind the surrogate models. Here, this paper evaluates the loss in accuracy and the gain in computational time while performing Bayesian model updating by using surrogate model evaluations compared to using direct FE model evaluations. This evaluation is crucial before entirely relying on surrogate models in model updating for structural health monitoring (SHM) and damage prognosis (DP) purposes. This paper also demonstrates Bayesian updating and surrogate model construction of large-scale high-fidelity FE models of infrastructure systems. In this regard, the miter gate structural system is considered as the testbed structure. Three predominant damage modes (loss of contact between gate and wall, loss of thickness due to corrosion, and loss of tension in the diagonal rods) are considered for model updating purposes. Bayesian model updating is performed using direct FE evaluations by leveraging parallel computing. Two types of surrogates, namely polynomial chaos expansion (PCE) and Gaussian process regression (GPR), are developed for the miter gate. Model updating is performed again using the trained surrogate models, and the updating results are compared with their counterparts obtained using the direct FE evaluation results. The posterior distribution of the FE model parameters obtained using the trained surrogates are sufficiently accurate with respect to the posterior obtained utilizing the direct FE evaluations. In addition, an approximate 4-fold decrease in the computational time was observed when using surrogate model evaluations instead of direct FE evaluations for model updating.

42 ENGINEERING↗

Bayes_Opt-SWMM: A Gaussian process-based Bayesian optimization tool for real-time flood modeling with SWMM

Real-time flood model plays a pivotal role in averting urban flood damage, particularly when there is minimal lead time for preparatory measures. However, urban flood modeling in real-time often contends with inherent uncertainties arising from input data uncertainty and parameter ambiguities. Here this study introduces a real-time calibration (RTC) tool called Bayes_Opt-SWMM, specifically tailored for real-time urban flood modeling and uncertainty optimization. This tool leverages the Gaussian process-based Bayesian optimization algorithm and interfaces seamlessly with the Stormwater Management Model (SWMM). It integrates real-time model forcing data and flood monitoring collected through sensors and gauges which are strategically placed within critical locations of urban drainage systems. Our approach hinges on the Surrogate Model based Uncertainty Optimization (SMUO) concept, providing an avenue for enhancing real-time flood modeling. Bayes_Opt-SWMM runs the optimization process using a surrogate model called Gaussian Process emulator with two inference methods: (1) the Gaussian Process (GP) model and (2) Markov Chain Monte Carlo (MCMC) algorithm in GP model (GP_MCMC). Furthermore, three acquisition functions, namely Expected Improvement (EI), Maximum Probability of Improvement (MPI), and Lower Confidence Bound (LCB), facilitate optimal parameter fitting within the surrogate models. The efficiency of GP-based surrogate models in learning SWMM model parameters, leads to an improved uncertainty quantification and accelerated real-time flood modeling in urban areas. Overall, Bayes_Opt-SWMM emerges as a cost-effective and valuable tool for real-time flood modeling and monitoring, with significant potential for managing intelligent storm water systems in urban environments.

54 ENVIRONMENTAL SCIENCES↗

A generic method for interfacing IEDs using low voltage interfaces to real-time simulators with hardware in the loop

This study presents an interface method (IM) to interconnect real-time simulators (RTSs) without amplifiers for different intelligent electronic devices (IEDs), such as protective relays, smart meters, and other devices. It is a significant topic for the power systems protection community because there are many protection devices with different requirements. These protection devices have pinouts that measure bus voltages/currents and breaker pole states and generate trip/close pulse signals. Most of the emulation test beds with RTSs and HIL need amplifiers to connect IEDs from different product manufacturers because the low-voltage interface levels are not available in some instruction manuals. Additionally, commercial amplifiers have cutoff frequency limitations and have an excessive cost when several IEDs are wired onto test beds. Then, the validation of protection and control systems with RTSs and IEDs in the loop from different vendors without amplifiers results in an unfeasible alternative. In this paper, the IM was based on using an interface box, power source, and RTS to identify IED inputs/outputs. Here, the sequence to recognize IED pinouts and calculate the current/voltage scaling factors is described. The IM was validated satisfactorily with RTS and a protection device (with unknown pinouts) in the loop.

42 ENGINEERING↗

Accelerating Laboratory-Based, All-Solid-State Battery Research and Development: Industrially Relevant Small-Batch Dry-Processing, Small and Low-Cost Test Fixtures, and Short-Tolerant Separators

This work presents improvements to mixing (for dry-processed electrodes), separator robustness, and cell testing jigs for all-solid-state batteries (ASSBs). These developments combine synergistically to enable rapid research and development of ASSB catholytes. A mechanical "kneader" is designed and built which emulates the mixing and fibrillation of poly-tetrafluoroethylene binder that occurs in industrially relevant twin-screw extruders. A modest improvement in the mixing and dispersion catholyte materials is observed using micro-resolution X-ray computed tomography (X-ray CT). Catholytes are paired with an In-metal anode and separated by a novel, dual-layer, polyaramid-fiber-supported, sulfide-solid electrolyte (Li6 PS5 Cl) separator. This separator achieves high short tolerance which enabled a 95% success rate across 20 attempted cells and to date over 100 successful cells have been built. ASSBs were tested in a 2032 coin-cell format. Simultaneous cycling of a large number of cells is further enabled by a compact and low cost (~ 38 USD per jig) pressure application jig of which 100 have been constructed to date. The utility of these advances is demonstrated through a brief study on the effect cathode-side conductive carbon and current collector type has on capacity and rate performance.

25 ENERGY STORAGE↗

Synergistic foam stabilization and transport improvement in simulated fractures with polyelectrolyte complex nanoparticles: Microscale observation using laser etched glass micromodels

Inaccessibility to direct pore scale observation in hydrocarbon recovery of tight shale formations poses a great challenge to water-energy nexus initiatives and necessitates the use of high throughput technologies to emulate environmentally friendly processes. Herein, we employ a precise glass micromodel fabrication and visualization method to isolate the supercritical CO 2 bubbles surrounded by CO 2 -water lamella prepared in saline produced water stabilized with molecular complexation of zwitterionic surfactants (ZS) and polyelectrolyte complex nanoparticles (PECNP). The Selective Laser Enhanced Etching (SLE) technique was selected for micromodel simulation of high-pressure flow. Two representative designs, (1) fracture/micro-crack network 28 and (2) fracture/matrix were etched on fused silica glass with a laser printing machine and scCO 2 foam was injected to study the foamability, propagation, stability, and fluid loss properties. The highly monodispersed and uniformly distributed array of scCO 2 bubbles were detected in flow of scCO 2 foam in highly saline brine containing ionic complexes of positively charged PECNPs and ZS, whereas foam flow with the lamella containing ZS in fractures offered a noticeably large and polydisperse array of scCO 2 bubbles. scCO 2 bubble motion and deformation were traced, and local description of foam flow was visually examined. The confined array of scCO 2 bubbles stabilized by ZS in microcracks was affected by bubble growth and coalescence, whereas the super-populated array of monodispersed scC O2 bubbles with lamella containing complexes of PECNP and ZS were able to fill the channels with stable configurations within the timeframe of comparative stability measurements. The ability of complex fluid to prevent the formation damage was evaluated through fluid loss visualization in micromodels. Probing scCO 2 foam transport in homogenous porous media revealed smaller volume leak-off for scCO 2 foam containing PECNP-ZS ionic complexes.

04 OIL SHALES AND TAR SANDS↗

Software-Defined Network for End-to-end Networked Science at the Exascale

Domain science applications and workflow processes are currently forced to view the network as an opaque infrastructure into which they inject data and hope that it emerges at the destination with an acceptable Quality of Experience. There is little ability for applications to interact with the network to exchange information, negotiate performance parameters, discover expected performance metrics, or receive status/troubleshooting information in real time. The work we presen here is motivated by a vision for a new smart network and smart application ecosystem that will provide a more deterministic and interactive environment for domain science workflows. The Software-Defined Network for End-to-end Networked Science at Exascale (SENSE) system includes a model-based architecture, implementation, and deployment which enables automated end-to-end network service instantiation across administrative domains. An intent based interface allows applications to express their high-level service requirements, an intelligent orchestrator and resource control systems allow for custom tailoring of scalability and real-time responsiveness based on individual application and infrastructure operator requirements. This allows the science applications to manage the network as a first-class schedulable resource as is the current practice for instruments, compute, and storage systems. Deployment and experiments on production networks and testbeds have validated SENSE functions and performance. Emulation based testing verified the scalability needed to support research and education infrastructures. Key contributions of this work include an architecture definition, reference implementation, and deployment. This provides the basis for further innovation of smart network services to accelerate scientific discovery in the era of big data, cloud computing, machine learning and artificial intelligence.

97 MATHEMATICS AND COMPUTING↗

Deep Learning Explicit Differentiable Predictive Control Laws for Buildings

We present a differentiable predictive control (DPC) methodology for learning constrained control laws for unknown nonlinear systems. DPC poses an approximate solution to multiparametric programming problems emerging from explicit nonlinear model predictive control (MPC). Contrary to approximate MPC, DPC does not require supervision by an expert controller. Instead, a system dynamics model is learned from a small dataset of recorded observations of the perturbed system's dynamics and the control law is optimized offline by interaction with the learned system model. The DPC method is based on two sequential steps, i) system identification using a constrained neural state-space model, and ii) optimization of an explicit control law parametrized by another neural network in closed-loop simulation with the identified neural state-space model. The combination of a differentiable closed-loop system and penalty methods for constraint handling of system outputs and inputs allows us to optimize the control law's parameters directly by backpropagating economic MPC loss through the learned system model. By incorporating domain knowledge and leveraging established techniques from optimal control, our method leverages deep neural networks as nonlinear function approximators for system identification and control while avoiding concomitant costs of intractably large datasets, and computationally expensive over-parametrized models. The scalability, data efficiency, and constrained optimal control capability of the proposed DPC method are demonstrated in simulation using a multi-zone building emulator.

Drgona, Jan↗

Coordinated frequency regulation among asynchronous AC grids with an MTDC system

With the increasing penetration level of renewable power, the power systems’ inertias are decreasing. To ensure that the frequency is above the threshold of the under frequency load shedding, one effective way is to share the spinning reserves among asynchronous AC grids through HVDC systems. In this paper, a coordinated frequency regulation scheme that enables spinning reserves sharing is proposed for such hybrid AC/DC grids. This scheme is a communication-free method that consists of the droop control, the inertia emulation control and the frequency safety control. The combined response of the three controls can increase not only the frequency nadir but also the settled frequency of the disturbed grid. In addition, the hysteresis comparators and the power injection limits are adopted in the frequency regulation scheme to reduce the coupling effect among AC grids. Thus, under small and harmless disturbance, the frequency regulation scheme will be inactive; the disturbance in strong AC grids would not affect the weak AC systems. The frequency regulation scheme coordinates well with the automatic generation control of the grids. Furthermore, the feasibility and effectiveness of the proposed frequency regulation scheme are verified through simulation studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Adaptive critic design-based reinforcement learning approach in controlling virtual inertia-based grid-connected inverters

In this report, an adaptive critic design (ACD) approach is proposed to control the phase and voltage of a grid-connected virtual synchronous generator (VSG). The penetration of fast responding inertia-less power converters significantly affect the stability of the power system, especially weak systems such as micro grids. The concept of virtual inertia addresses this concern by virtually emulating the behavior of a synchronous generator. However, the conventional VSG is designed based on two conditions: (i) fixed operating point and (ii) inductive grid connections. The performance of VSGs in low-voltage semi-resistive microgrids is far from optimal. To overcome the aforementioned concerns, a heuristic dynamic programing (HDP) approach is proposed to optimally control grid-connected VSGs. The neural-network-based inherence of the HDP enables the proposed technique to adapt to any impedance angle. The HDP controller includes two subnetworks: (i) the action network that controls the system optimally and (ii) the critic network, which evaluates the effectiveness of the action network. The simulation and experimental results are provided to evaluate the effectiveness of the proposed technique. As shown, the HDP-based approach illustrates a better performance in comparison with the conventional PI-based VSG in various operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A numerical study of an atomizing jet in a resonant acoustic field

By applying state-of-the-art, high-fidelity methods for compressible multiphase flows, this work features a parametric study of a turbulent liquid jet atomizing under the influence of standing acoustic waves with five different amplitude values. We perform these simulations in a semi-periodic domain to isolate interactions between acoustic and atomization processes and to emulate a pure liquid fuel spray placed at the pressure node of a resonant acoustic field. Characterization of the flow includes measurements of the mean dimensions, droplet number, and surface area of the liquid distribution, alongside qualitative depictions of the interface evolution, breakup events, and hydrodynamic fields. As acoustic forcing is provided to the flow, the added energy provokes atomization when shear would otherwise be insufficient for breakup. The magnitude of the forcing determines the rate of droplet production, which accelerates as the evolving liquid surface modifies the induced motion of the gas phase. At the nodal plane, the acoustic field produces an organizing effect on the liquid structures.

42 ENGINEERING↗

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗

CAMERA: A method for cost-aware, adaptive, multifidelity, efficient reliability analysis

Estimating probability of failure in aerospace systems is a critical requirement for flight certification and qualification. Failure probability estimation involves resolving tails of probability distributions, and Monte Carlo sampling methods are intractable when expensive high-fidelity simulations have to be queried. Here, we propose a method to use models of multiple fidelities that trade accuracy for computational efficiency. Specifically, we propose the use of multifidelity Gaussian process models to efficiently fuse models at multiple fidelity, thereby offering a cheap surrogate model that emulates the original model at all fidelities. Furthermore, we propose a novel sequential acquisition function based experiment design framework that can automatically select samples from appropriate fidelity models to make predictions about quantities of interest at the highest fidelity. We use our proposed approach in an importance sampling setting and demonstrate our method on the failure level set and probability estimation on synthetic test functions and two real-world applications, namely, the reliability analysis of a gas turbine engine blade using a finite element method and a transonic aerodynamic wing test case using Reynolds-averaged Navier-Stokes equations. We show that our method predicts the failure boundary and probability more accurately and at a fraction of the computational cost compared with using just a single expensive high-fidelity model. Finally, we show that our sequential approach is guaranteed to asymptotically converge to the true failure boundary with high probability.

97 MATHEMATICS AND COMPUTING↗

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗

Estimating QSVT angles for matrix inversion with large condition numbers

Quantum Singular Value Transformation (QSVT) is a state-of-the-art, near-optimal quantum algorithm that can be used for matrix inversion. The QSVT circuit is parameterized by a sequence of angles that must be pre-calculated classically, with the number of angles increasing as the matrix condition number grows. Computing QSVT angles for ill-conditioned problems is a numerically challenging task. Here, we propose a numerical technique for estimating QSVT angles for large condition numbers. This technique allows one to avoid expensive numerical computations of QSVT angles and to emulate QSVT circuits for solving ill-conditioned problems.

97 MATHEMATICS AND COMPUTING↗

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗