Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “integral equations”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Is the Adoption of Electric Vehicles (EVs) and Solar Photovoltaics (PVs) Interdependent or Independent? An Integrated EVs-PVs Modeling Framework

Transportation and residential energy consumption account for more than one-half of the overall energy consumption in the United States. Adoption of electric vehicles (EVs) can play a key role in decarbonizing the transportation sector, while the adoption of renewable energy sources (e.g., solar photovoltaics or PVs) could bring similar benefits to the residential energy sector and in turn support transport electrification. Although the market share for both EVs and PVs continue to grow, both these emerging technologies are deployed rather disjoint without considering the existence of potential similarities among users who own (or aspire to own) EVs and PVs. This might be due to lack of understanding of the behavioral interdependence in consumer preferences towards these technologies. To fill this gap in knowledge, this study utilizes data from the 2018 WholeTraveler Transportation Behavior Study to develop an integrated model system that explores EV and PV adoption behaviors. A structural equations model (SEM) is employed that incorporates direct effects as well as error correlations among the adoption behaviors for EVs and PVs. Model results indicate that the adoption behavior for both these technologies is indeed interconnected and significantly influenced by attitudes, values and personality traits. Findings from this research suggest that incentives (e.g., subsidies) that drive 'bundled' adoption of EVs and PVs could accelerate the adoption of both these sustainable technologies. The results highlights the need to consider transport and building energy efficient technology adoption behavior in a single integrated structure.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2026

Scalable dark matter searches using integrated photonics

Dark matter (DM) with masses of order an electronvolt or below can have a nonzero coupling to electromagnetism while being compatible with cosmological observations. In these models, the ambient DM behaves as a new classical source in Maxwell’s equations, which can excite potentially detectable electromagnetic (EM) fields in the laboratory. We propose a new integrated-photonics–based approach to search for dark matter candidates in the 0.1–few eV mass range. This approach offers a wide range of wavelength-scale devices like resonators and wave guides that are readily fabricated in large quantities, enabling a scalable and novel search. In particular, we demonstrate that refractive index-modulated resonators, such as etched/grooved microrings, or patterned slabs, support EM modes with efficient coupling to DM. When excited by DM, these modes are read out by coupling the resonators to a wave guide that terminates on a micron-scale–sized single photon detector, such as a single pixel of a low-noise charge-coupled device or a superconducting nanowire. We then estimate the sensitivity of this experimental concept in the context of axionlike particle and dark photon models of DM, demonstrating that nanophotonic confinement and scalability can extend dark matter sensitivity into previously unexplored parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Power System Frequency Dynamics Modeling, State Estimation, and Control using Neural Ordinary Differential Equations (NODEs) and Soft Actor-Critic (SAC) Machine Learning Approaches

With the global energy transition of the electric power system, grid control, supervision, and protection is becoming more challenging. With the increasing integration of renewable energy sources (RES), the system dynamics are changing, causing traditional power system dynamic modeling with swing equation-based modeling approaches to fail. Additionally, the converter-dominated power grid is decreasing the system inertia, making the power system more fragile to the frequency swings. This paper first investigates and compares the application of a model-based Kalman filter state estimation approach with (i) a model-free machine learning approach --- neural ordinary differential equations (NODEs) --- and (ii) a data-driven system identification (SysId) approach to model and infer critical state values of the power system frequency dynamics. Then a model predictive control (MPC) framework is compared to a model-free Soft Actor-Critic (SAC) reinforcement learning (RL) control algorithm in providing efficient fast frequency response (FFR) to the power system frequency dynamics. The approaches are compared in terms of their performance goals as well as their per-timestep computational efficiency. Furthermore, the comparative study for state estimation shows that for the model-free requirement, both NODEs and SysId can provide accurate state estimates; however, with increasing model complexity, NODEs can be a better choice for model identification. Similarly, the results from the FFR comparative study show that the SAC RL-based FFR, once trained, outperforms MPC with better control signals and faster computation time, making the SAC RL-based FFR better option for providing FFR to the power system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Quantized State Integrator With Second Order Errors Over Monotonic Segments

We introduce a novel quantized state integration method that is analogous to the explicit, second order Runge-Kutta technique. A feature of this method is that, for a single differential equation over an interval where the solution is monotonic, it exhibits an error proportional to the square of the integration quantum. We offer a theoretical proof of this behavior and demonstrate it with a numerical example. At the same time, an extension of the method to a system with input causes the errors to become proportional to the integration quantum. The latter observation fits established theory; the former appears to be a new result that might offer new insights into the construction of quantized state integration methods.

Mahawattege, Rasika↗
DOE OSTI2021

Fishnet four-point integrals: integrable representations and thermodynamic limits

In this work, we consider four-point integrals arising in the planar limit of the conformal “fishnet” theory in four dimensions. They define a two-parameter family of higher-loop Feynman integrals, which extend the series of ladder integrals and were argued, based on integrability and analyticity, to admit matrix-model-like integral and determinantal representations. In this paper, we prove the equivalence of all these representations using exact summation and integration techniques. We then analyze the large-order behaviour, corresponding to the thermodynamic limit of a large fishnet graph. The saddle-point equations are found to match known two-cut singular equations arising in matrix models, enabling us to obtain a concise parametric expression for the free-energy density in terms of complete elliptic integrals. Interestingly, the latter depends non-trivially on the fishnet aspect ratio and differs from a scaling formula due to Zamolodchikov for large periodic fishnets, suggesting a strong sensitivity to the boundary conditions. We also find an intriguing connection between the saddle-point equation and the equation describing the Frolov-Tseytlin spinning string in AdS 3 × S 1 , in a generalized scaling combining the thermodynamic and short-distance limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A high fidelity and user-friendly equation-oriented optimization model for carbon capture using a novel water-lean solvent

Research Triangle Institute (RTI) International and SLB have developed a novel water-lean solvent technology for carbon capture, demonstrating low specific reboiler duty (SRD) values at capture rates exceeding 90%. At the Technology Centre Mongstad (TCM) pilot plant, the technology achieved an SRD of 2.55 GJ/t-CO2 at 95% capture, utilizing an intercooler and a 5°C temperature approach in the lean/rich solvent cross exchanger. To meet varying carbon capture targets for Front End Engineering and Design (FEED) studies and to enable real-time optimization and advanced process control, an efficient optimization model is required. This model needs to minimize energy demand for a given capture rate and determine optimal operating parameters in response to fluctuating flue gas conditions. While an existing Aspen Plus simulation model, developed by RTI and SLB, accurately matches TCM plant data, its sequential modular (SM) strategy is too slow for real-time applications due to recycle streams and tight heat integration inherent in solvent-based carbon capture processes. Although an equation-oriented (EO) modeling strategy is more suitable for optimizing these processes, its adoption has been limited by several factors: feature limitations in Aspen Plus EO mode (e.g., lack of balance block support), a less user-friendly interface for variable identification and loop solving, complex troubleshooting of convergence issues, and the necessity for accurate initial values.

carbon capture↗
DOE OSTI2023

Dynamic modelling of slip in a wind turbine spherical roller main bearing

This paper considers the problem of the dynamic modelling of macro slip in spherical roller bearings. By revisiting the fundamental physics which drive these systems, potential issues in existing models have been identified. Furthermore, in pure rolling conditions it was found that governing differential equations become "stiff", requiring the use of implicit methods of time integration. The problem of individual roller macro slip in a wind turbine main bearing is then investigated using a simplified representation of system dynamics. Model results indicate clear links between slip/friction and the operational strategy of the wind turbine, as well as significantly higher frictional effects in the downwind main bearing row. Due to modelling simplifications, these results should not yet be considered conclusive, with further work required.

17 WIND ENERGY↗
DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Matlab implementation of a novel semi-structured kinetic model for methanotroph-photoautotroph cocultures

This paper presents the MatLab implementation details of a novel semi-structured kinetic model for methanotroph-photoautotroph cocultures. This includes the parameterization of the modeling equations, and the initialization of the simulation based on experimental conditions. More importantly, it provides details on how the differential equations governing mass balances in both gas and liquid phases are integrated together to simulate the system dynamics over time. The semi-structured kinetic model for methanotroph-photoautotroph coculture is validated using a wide range of experimental conditions. The model: Accurate predicts both the coculture growth in liquid phase and the gas composition changes in head space over time; Explicitly models the exchange of in situ produced O 2 and CO 2 within the coculture; Considers the self-shading effect on the growth of photoautotroph.

42 ENGINEERING↗
DOE OSTI2026

Local Pair Natural Orbital-Based Coupled-Cluster Theory through Full Quadruples (DLPNO–CCSDTQ)

In this work, we implement a local pair natural orbitalbased coupled-cluster method through the full treatment of quadruple excitations (CCSDTQ). The domain-based local pair natural orbital (DLPNO) approach, which has successfully been applied to lower levels of coupled-cluster theory, is utilized in our algorithm, and thus our algorithm is called DLPNO-CCSDTQ. For simplicity in the working equations and in the implementation, we t 1 -dress the twoelectron integrals as well as Fock matrix elements. Our method can recover CCSDTQ-CCSDT and CCSDTQ-CCSDT(Q) energy differences on the order of 0.01−0.05 kcal mol −1 , even at a loose quadruples natural orbital (QNO) occupation number cutoff of 3.33 × 10 −6 . To highlight the capabilities of our code and its potential future applications, we showcase computations that would be intractable with canonical CCSDTQ, such as the benzene dimer, (H 2 O) 17 , and adamantane. With sufficient computing resources, computations up to 15 heavy atoms (40 atoms overall) may be feasible for fully bonded 3D systems.

Cluster chemistry↗
DOE OSTI2020

Absolute stimulated Raman side scatter in direct-drive laser-produced plasmas

Production of suprathermal electrons by stimulated Raman scattering (SRS) is a principal concern for contemporary direct-drive inertial confinement fusion experiments at the National Ignition Facility and similar systems since such electrons penetrate and preheat the target core, preventing efficient implosion. The higher temperatures and longer scale lengths in these experiments favor SRS over two-plasmon decay, which predominated in earlier experiments. In particular, current experiments are expected to exceed the threshold for absolute Raman side scatter, which would then dominate the interaction since it grows temporally until saturated by nonlinear mechanisms such as hot-electron production. Until recently, analyses of SRS side scatter have treated the case of a single laser beam incident on a plasma, but the direct-drive approach to laser fusion employs a multitude of beams to drive the implosion. Here in this paper, we present an analysis that can be applied to an arbitrary number of beams with varied angles of incidence and polarizations. In the case of a single beam, it allows a physically motivated derivation and verification of an analytic threshold formula. In the general case of multiple beams and arbitrary orientation and polarizations, the threshold is found by numerical integration of a set of first-order linear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

MEUMAPPS (C++ Version)

Many materials, metal alloys in particular, have features on the on micrometer or nanometer scale that have a large impact on the properties of the material. These features are known as the microstructure of the material. Understanding why and how the microstructure forms in a material is of fundamental scientific interest as well as of significant technological interest. The capability to predict microstructure evolution in a material allows the intentional design of microstructures and hence the intentional design of material properties. The phase-field method is one of the leading methods for predicting microstructure evolution. One of the most significant problems for phase-field models is their computational expense. Even limited phase-field simulations can easily require thousands of CPU core-hours to complete, which significantly limits their use. This code provides both a general framework for creating scalable, GPU-accelerated phase-field model applications as well as several applications themselves. The code is capable of using hundreds of GPUs efficiently, which greatly reduces the time required to perform simulations. The code is written with an emphasis on performance portability, that is the ability for the code to run efficiently on a number of different computing architectures without modification of the source code. The performance portability of this code is primarily enabled through the use of two libraries, Kokkos (performance portable data structures and execution patterns) and heFFTe (performance portable distributed 3D fast Fourier transforms). The code consists of a core library, applications, and tests. The core library includes shared functionality between applications. This includes interfaces with fast Fourier transform (FFT) libraries such as heFFTe, data structures based on Kokkos, file input and output capabilities, and a solver for infinitesimal strain mechanical equilibrium problems. Five applications are included in the code. The flagship application is the MEUMAPPS-SS application, which implements the Kim-Kim-Suzuki phase-field model for precipitation for an arbitrary number of phases and components in a metal alloy. Five simpler applications are also included that solve the Eshelby inclusion problem, Allen-Cahn equation, the coupled Allen-Cahn and diffusion equations, and the Cahn-Hilliard equation. The code includes two applications to solve the Cahn-Hilliard equation, one with constant-step-size first-order time integration and the second with adaptive high-order time integration.

DeWitt, Stephen [Oak Ridge National Lab. (ORNL), O↗
DOE OSTI2022

Statistics of Green's functions on a disordered Cayley tree and the validity of forward scattering approximation

The accuracy of the forward scattering approximation for two-point Green's functions of the Anderson localization model on the Cayley tree is studied. A relationship between the moments of the Green's function and the largest eigenvalue of the linearized transfer-matrix equation is proved in the framework of the supersymmetric functional-integral method. The new large-disorder approximation for this eigenvalue is derived and its accuracy is established. Using this approximation the probability distribution of the two-point Green's function is found and compared with that in the forward scattering approximation (FSA). It is shown that FSA overestimates the role of resonances and thus the probability for the Green's function to be significantly larger than its typical value. The error of FSA increases with increasing the distance between points in a two-point Green's function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Passive Ultrasonic Deterrents to Reduce Bat Mortality in Wind Farms (Final Report)

The overall objective of this project was to develop active and passive ultrasonic deterrent concepts to mitigate bat mortality at wind turbines. The research was supported by the US Department of Energy under the contract DE-EE0008731. Novel ultrasonic bat deterrents are investigated in this study. The deterrents are based on the idea of aerodynamic whistles wherein flow-acoustic resonance is used to produce high-amplitude tonal sound at desired ultrasonic frequencies. The deterrents can be grouped into active and passive. The active deterrents/whistles are driven by pressurized air supplied by an external source (e.g., an air compressor). The passive deterrents/whistles are designed to be mounted on wind turbine blades and are “powered” by the energy in the air moving past the rotor blades. Multiple active and passive deterrent ideas are proposed and investigated numerically and experimentally. The computational analysis involves solving the compressible unsteady Reynolds-averaged Navier-Stokes (URANS) equations and coupling the near-field aerodynamic solution with an integral method based on the Ffowcs Williams-Hawkings acoustic analogy to predict ultrasound in the farfield. Experiments are performed in the anechoic chamber at ISU (for the active whistles) and in the Stability wind tunnel at Virginia Tech (for the passive whistles).

17 WIND ENERGY↗
DOE OSTI2026

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗
DOE OSTI2022

Universal Kardar-Parisi-Zhang Dynamics in Integrable Quantum Systems

Although the Bethe ansatz solution of the spin-1/2 Heisenberg model dates back nearly a century, the anomalous nature of its high-temperature transport dynamics has only recently been uncovered. Indeed, numerical and experimental observations have demonstrated that spin transport in this paradigmatic model falls into the Kardar-Parisi-Zhang (KPZ) universality class. This has inspired the significantly stronger conjecture that KPZ dynamics, in fact, occur in all integrable spin chains with non-Abelian symmetry. Here, in this letter, we provide extensive numerical evidence affirming this conjecture. Moreover, we observe that KPZ transport is even more generic, arising in both supersymmetric and periodically driven models. Motivated by recent advances in the realization of SU(N)-symmetric spin models in alkaline-earth-based optical lattice experiments, we propose and analyze a protocol to directly investigate the KPZ scaling function in such systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Multigrid-Reduction-in-Time for the Rotating Shallow Water Equations

We consider multilevel time-parallel methods for the numerical solution of the rotating shallow water equations. In particular, the multigrid-reduction-in-time (MGRIT) algorithm is used for the parallel time integration. An asymptotic model is used at the coarse levels while the full model is employed at the finer levels. The asymptotic model is well-suited for highly oscillatory partial differential equations like the rotating shallow water equations because it can accurately and stably take the required large time-steps on coarse levels. Our work exploits the flexibility of the MGRIT algorithm in terms of the number of levels and relaxation schemes to show some computational benefits, especially with respect to FCF-relaxation and data reuse.

97 MATHEMATICS AND COMPUTING↗