Search NASA⌕ Search

SEARCH · Search NASA

Results for “Laplace transform”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Monte Carlo Laplace Transform Estimator for Radiation Transport

This work formulates and implements a Laplace transform estimator in a simple Monte Carlo radiation transport code. The estimator maps flux-based quantities of interest, like reaction rates, from a desired phase-space dimension to the complex Laplace domain. This on-the-fly Monte Carlo integration technique enables the spectral analysis of arbitrary nuclear systems via the Laplace transform. A simple code tests the estimator in neutron slowing-down problems across various infinite media, and the results compare well with Ganapol’s uninverted analytical solution of the neutron slowing-down equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Laplace Transform–Based Quantum Eigenvalue Transformation via Linear Combination of Hamiltonian Simulation

Eigenvalue transformations, which include solving time-dependent differential equations as a special case, have a wide range of applications in scientific and engineering computation. While quantum algorithms for singular value transformations are well studied, eigenvalue transformations are distinct, especially for nonnormal matrices. Here, we propose an efficient quantum algorithm for performing a class of eigenvalue transformations that can be expressed as a certain type of matrix Laplace transformation. This allows us to significantly extend the recently developed linear combination of Hamiltonian simulation method [D. An, J.-P. Liu, and L. Lin, Phys. Rev. Lett., 131 (2023), 150603; D. An, A. M. Childs, and L. Lin, Commun. Math. Phys. 407, 19 (2026)] to represent a wider class of eigenvalue transformations, such as powers of the matrix inverse, 𝐴 −𝑘 , and the exponential of the matrix inverse, 𝑒 −𝐴 −1 . The latter can be interpreted as the solution of a mass-matrix differential equation of the form form 𝐴⁢𝑢′⁡⁡(𝑡) =−𝑢⁡(𝑡). We demonstrate that our eigenvalue transformation approach can solve this problem without explicitly inverting 𝐴, thereby reducing the computational complexity.

Laplace transform↗
DOE OSTI2021

Machine-learning-based inversion of nuclear responses

A microscopic description of the interaction of atomic nuclei with external electroweak probes is required for elucidating aspects of short-range nuclear dynamics and for the correct interpretation of neutrino oscillation experiments. Nuclear quantum Monte Carlo methods infer the nuclear electroweak response functions from their Laplace transforms. Inverting the Laplace transform is a notoriously ill-posed problem; and Bayesian techniques, such as maximum entropy, are typically used to reconstruct the original response functions in the quasielastic region. In this work, we present a physics-informed artificial neural network architecture suitable for approximating the inverse of the Laplace transform. Utilizing simulated, albeit realistic, electromagnetic response functions, we show that this physics-informed artificial neural network outperforms maximum entropy in both the low-energy transfer and the quasielastic regions, thereby allowing for robust calculations of electron scattering and neutrino scattering on nuclei and inclusive muon capture rates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Second Order System Study

During my education in mathematics, engineering and physics, I learned transform pairs and their usage mechanics but I never remember seeing the derivations of the solutions to second order ordinary differential equations (ODE) and difference equations. A solution to a question posed in a potential funder meeting put me on a path to solving second order systems using the five principal Fourier based methods: Fourier transform (FT), Z-transform (ZT), discrete time Fourier transform (DTFT), discrete Fourier transform (DFT) and Laplace transform (LT).

42 ENGINEERING↗
DOE OSTI2026

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fault isolation and fault-tolerant control for Takagi-Sugeno fuzzy time-varying delay stochastic distribution systems

A fault isolation, estimation, and fault-tolerant control (FTC) scheme for nonlinear time-varying delay stochastic distribution control systems was presented in this paper. The Takagi-Sugeno fuzzy model was adopted to approach the nonlinear dynamics of time-varying delay systems. According to the output equivalence principle and Laplace transformation, an augmented state vector was given to solve the time-varying delay problem. When multiple actuator faults and interference occur simultaneously, fault detection, isolation and fault estimation was designed to obtained the fault information. To decouple faults and obtain the value and location information of the fault, the system was separated into two parts through the designed multiple conversion matrices, in which one subsystem was only affected by one actuator fault. This has simplified the design of fault isolation and estimation. A adaptive observer for fault estimation was given. Then, fault information such as the time, location, and size was determined. The observer gain matrices were calculated using linear matrix inequality (LMI). When a fault was detected and diagnosed, a FTC algorithm was devised using the proportional-integral control scheme to compensate the fault as much as possible. It has been shown that even if multiple faults actuator occurred simultaneously, the FTC controller still ensured the output probability density function of the system traced the desired probability density function when a fault occurred. Finally, the expected results were obtained through the simulation example, which confirmed the effectiveness of the method.

42 ENGINEERING↗
DOE OSTI2025

Current interrupt method for calculating the electrochemical impedance in a solid oxide electrolysis stack

Here, in this work the time domain response of Solid Oxide Electrolysis Cells (SOEC) to a current interruption was transformed into the frequency domain using a carrier function Laplace transform, which is fit to the experimental data using a MATLAB Complex Nonlinear Least Squares (CNLS) solver. The hardware implementation, consisting principally of a high-speed switch and a fast-logging Analog to Digital Converter (ADC), was assembled and tested using a calibration module to assess the accuracy, repeatability, and speed of acquisition of the prototype device as compared against a calibrated commercial impedance spectrometer. Additionally, the current interrupt device and commercial FRA were used to acquire the impedance spectra of a four cell SOEC stack with a large, 300 cm 2 , active cell area.

SOEC↗
DOE OSTI2024

Uncertainty-quantification-enabled inversion of nuclear responses

Nuclear quantum many-body methods rely on integral transform techniques to infer properties of electroweak response functions from ground-state expectation values. Retrieving the energy dependence of these responses is highly nontrivial, especially for quantum Monte Carlo methods, as it requires inverting the Laplace transform, a notoriously ill-posed problem. Here, in this work, we propose an artificial neural network architecture suitable for accurate response function reconstruction with precise estimation of the uncertainty of the inversion. We demonstrate the capabilities of this new architecture benchmarking it against maximum entropy and previously developed neural network methods designed for a similar task, paying particular attention to its robustness noise in the Euclidean

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software

The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).

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

Dirichlet-type absorbing boundary conditions for peridynamic scalar waves in two-dimensional viscous media

Construction of absorbing boundary conditions (ABCs) for nonlocal models is generally challenging, primarily due to the fact that nonlocal operators are commonly associated with volume constrained boundary conditions. Moreover, application of Fourier and Laplace transforms, which are essential for the majority of available methods for ABCs, to nonlocal models is complicated. In this paper, we propose a simple method to construct accurate ABCs for peridynamic scalar wave-type problems in viscous media. The proposed ABCs are constructed in the time and space domains and are of Dirichlet type. Consequently, their implementation is relatively simple, since no derivatives of the wave field are required. The proposed ABCs are derived at the continuum level, from a semi-analytical solution of the exterior domain using harmonic exponential basis functions in space and time (plane-wave modes). The numerical implementation is done using a meshfree collocation approach employed within a boundary layer adjacent to the interior domain boundary. The modes satisfy the peridynamic numerical dispersion relation, resulting in a compatible solution of the interior region (near-field) with that of the exterior region (far-field). The accuracy and stability of the proposed ABCs are demonstrated with several numerical examples in two-dimensional unbounded domains.

42 ENGINEERING↗
DOE OSTI2023

Analyticity of replica correlators and modular ETH

We study the two point correlation function of a local operator on an n -sheeted replica manifold corresponding to the half-space in the vacuum state of a conformal field theory. In analogy with the inverse Laplace transform, we define the Renyi transform of this correlation function, which is a function of one complex variable w, dual to the Renyi parameter n . Inspired by the inversion formula of Caron-Huot, we argue that if the Renyi transform ƒ(w) has bounded behavior at infinity in the complex w plane, the discontinuity of the Renyi transform disc ƒ(w) provides the unique analytic continuation in n of the original replica correlation function. We check our formula by explicitly calculating the Renyi transform of a particular replica correlator in a large N holographic CFT d in dimensions d > 2. We also discover that the discontinuity of the Renyi transform is related to the matrix element of local operators between two distinct eigenstates of the modular Hamiltonian. We calculate the Renyi transform in 2 d conformal field theories, and use it to extract the off-diagonal elements of (modular) ETH. We argue that in 2 d , this is equivalent to the off-diagonal OPE coefficients of a CFT and show that our technique exactly reproduces recent results in the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Orientation reversal and the Chern-Simons natural boundary

We show that the fundamental property of preservation of relations, underlying resurgent analysis, provides a new perspective on crossing a natural boundary, an important general problem in theoretical and mathematical physics. This reveals a deeper rigidity aspect of resurgence in a quantum field theory path integral. The physical context here is the non-perturbative completion of complex Chern-Simons theory that associates to a 3-manifold a collection of q-series invariants labeled by Spinc structures, for which crossing the natural boundary corresponds to orientation reversal of the 3-manifold. Our new resurgent perspective leads to a practical numerical algorithm that generates q-series which are dual to unary q-series composed of false theta functions. Until recently, these duals were only known in a limited number of cases, essentially based on Ramanujan’s mock theta functions, and the common belief was that the duals might not even exist in the general case. Resurgence analysis identifies as primary objects Mordell integrals: up to changes of variables, they are Laplace transforms of resurgent functions. Their unique Borel summed transseries decomposition on either side of the Stokes line is simply the unique decomposition into real and imaginary parts. In turn, the latter are combinations of unary q-series in terms of q and its modular counterpart $\overset{\sim }{q}$ , and are resurgent by construction. The Mordell integral is analytic across the natural boundary of the q and $\overset{\sim }{q}$ series, and uniqueness of a similar decomposition which preserves algebraic relations on the other side of the boundary defines the unique boundary crossing of the q series. We demonstrate that this continuation can be efficiently implemented numerically. In the cases where unique mock modular identities are known, they are found by this numerical procedure, but the procedure can go well beyond the known list of identities. A particularly interesting feature of the resurgent approach is that it reveals new aspects, and is very different from other known approaches based on indefinite theta series, Appell-Lerch sums, and representation theory of logarithmic vertex operator algebras.

Chern-Simons theories↗
DOE OSTI2023

Peridynamic elastic waves in two-dimensional unbounded domains: Construction of nonlocal Dirichlet-type absorbing boundary conditions

The focus of this paper is on application of peridynamics (PD) to propagation of elastic waves in unbounded domains. We construct absorbing boundary conditions (ABCs) derived from a semi-analytical solution of the PD governing equation at the exterior region. This solution is made up of a finite series of plane waves, as fundamental solutions (modes), which satisfy the PD dispersion relations. The modes are adjusted to transmit the energy from the interior region (near field) to the exterior region (far field). The corresponding unknown coefficients of the series are found in terms of the displacement field at a layer of points adjacent to the absorbing boundary. This is accomplished through a collocation procedure at subregions (clouds) around each absorbing point. The proposed ABCs offer appealing advantages, which facilitate their application to PD. They are of Dirichlet-type, hence their implementation is relatively simple as no derivatives of the field variables are required. They are constructed in the time and space domains and thus application of Fourier and Laplace transforms, cumbersome for nonlocal models, is not required. At the discrete level, the modes satisfy the same numerical dispersion relations of the near field, which makes the far-field solution compatible with that of the near field. We scrutinize the performance of the proposed ABCs through several examples. So our investigation shows that the proposed ABCs perform stably in time with an appropriate level of accuracy even in problems characterized by highly-dispersive propagating waves, including crack propagation in semi-unbounded brittle solids.

42 ENGINEERING↗
DOE OSTI2021

Evaluating dynamic thermal performance of building envelope components using small-scale calibrated hot box tests

The hot box test method has been applied to evaluate both the steady-state (U-value) and dynamic thermal properties of building envelopes. However, the high construction cost of full-scale hot box apparatus and the testing time required (usually several days) may prevent its wider adoption. To overcome the limitations of full-scale hot box tests, here we propose a novel method to evaluate the dynamic thermal performance of building envelope components using a small-scale calibrated hot box and scaled-down specimen. In this paper, the scaling relationships of thermal properties evaluated using a full-size specimen and a scaled-down specimen are established based on the Laplace transform of the heat transfer equations. In addition, dynamic thermal properties obtained from scaled-down experimental tests are compared to the values calculated by the EN ISO13768 (ISO) method. A small-scale hot box with a 355 mm × 355 mm metering area was constructed and calibrated to validate the correlations. Three scaled-down concrete sandwich wall panels were then tested and the scaling relationship was cross-validated using the experimental results, finite difference (FD) simulations, and the ISO method. The results indicate that the dynamic thermal properties obtained from a scaled-down hot box test can be correlated to its full-size counterpart when certain conditions are met. The scaled-down hot box test is demonstrated to be an effective yet economical alternative to a full-scale test with significantly reduced experimentation cost and turn-around time.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Pressure falloff testing to characterize CO 2 plume and dry-out zone during CO 2 injection in saline aquifers

The fluid mobilities and the lateral extent of the CO 2 plume and the dry-out zone are among important unknowns required to effectively manage CO 2 storage in saline aquifers projects. These unknowns can be derived from pressure falloff testing which has been widely used to characterize the subsurface. This paper provides a three-region analytical solution for falloff pressure during the infinite-acting period of CO 2 injection in saline aquifers. The solution is derived using Laplace transformation technique of the governing diffusivity equations and the accompanying initial and boundary conditions. The exact solution is validated against numerical results obtained from a black-oil reservoir simulator for two cases. The solution is approximated into simple forms in real-time domain considering different time periods. Then, a graphical semi-log analysis technique based on the approximate solutions is provided as an interpretation tool of falloff tests where the fluid mobilities and the extents of dry-out zone and CO 2 plume can be inferred. The interpretation technique is applied to synthetic falloff test data to investigate its potential and reliability.

02 PETROLEUM↗
DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗
DOE OSTI2021

Emergence of spin singlets with inhomogeneous gaps in the kagome lattice Heisenberg antiferromagnets Zn-barlowite and herbertsmithite

The kagome Heisenberg antiferromagnet formed by frustrated spins arranged in a lattice of corner-sharing triangles is a prime candidate for hosting a quantum spin liquid (QSL) ground state consisting of entangled spin singlets. However, the existence of various competing states makes a convincing theoretical prediction of the QSL ground state difficult, calling for experimental clues from model materials. The kagome lattice materials Zn-barlowite (ZnCu 3 (OD) 6 FBr) and herbertsmithite (ZnCu 3 (OD) 6 Cl 2 ) do not exhibit long-range order and are considered the best realizations of the kagome Heisenberg antiferromagnet known so far. In this study, we use 63 Cu nuclear quadrupole resonance combined with the inverse Laplace transform to locally probe the inhomogeneity of delicate quantum ground states affected by disorder. We present direct evidence for the gradual emergence of spin singlets with spatially varying excitation gaps, but even at temperatures far below the super-exchange energy scale their fraction is limited to ~60% of the total spins. Theoretical models need to incorporate the role of disorder to account for the observed inhomogeneously gapped behaviour.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.

Tao, Xuecheng↗