Search NASA⌕ Search

SEARCH · Search NASA

Results for “LAPLACE TRANSFORMATION”

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 OSTI2018

A ONE-WAY COUPLED FAST BURST REACTOR BENCHMARK THROUGH LAPLACE TRANSFORM

In several previous publications, a slab Fast Burst Reactor (FBR), characterized by four non-linear coupled transient conservation equations for the neutron flux ?, material density and displacement ?, u and temperature T, is considered. In this work, before addressing a simplified benchmark version, we revisit the consistency of the four conservation equations. We then develop a One-Way Coupled (OWC) high precision linear benchmark in preparation for full FBR non-linear multiphysics. As will be shown, the importance of the OWC benchmark enables a direct comparison with a MOOSE? based simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2018

Semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model

We present a semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model. It combines Laplace transform and finite difference schemes. The Laplace transform handles the time variable whereas the finite difference method deals with the spatial coordinate. This semi-numerical scheme is not restricted by space discretization and allows the computation of a solution at any time without compromising numerical stability or the mass conservation principle. Our formulation results in a non-analytically-solvable second-order differential equation whose spatial discretization yields a tridiagonal linear algebraic system. Moreover, we describe comparisons between semi-numerical and semi-analytical solutions for particular cases. Results agree well with those from semi-analytic solutions. Furthermore, we expose a parametric analysis from the coupled model in order to show the effects of relevant parameters on pressure profiles and flow rates for the case where neither analytic nor semi-analytic solutions are available.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2018

On Integral of a Semi-Markov Diffusion Process

Let (X(t)) (t ≥ 0) be a semi-Markov diffusion process. The process (J(T )) (T ≥ 0) equal to the integral of (X(t)) on interval [0, T ) is studied. The relation between one-dimensional differential equation of the second order of elliptical type and asymptotics of a solution to Dirichlet problem on an interval with length tending to zero is established. This relation is used to derive a differential equation for the Laplace transform of the semi-Markov generating function of the process (J(t)).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Boundary node Petrov–Galerkin method in solid structures

Based on the interpolation of the Lagrange series and the Finite Block Method (FBM), the formulations of the Boundary Node Petrov–Galerkin Method (BNPGM) are presented in the weak form in this paper and their applications are demonstrated to the elasticity of functionally graded materials, subjected to static and dynamic loads. By introducing the mapping technique, a block of quadratic type is transformed from the Cartesian coordinate to the normalized coordinate with 8 seeds for two-dimensional problems. The first-order partial differential matrices of boundary nodes are obtained in terms of the nodal values of the boundary node, which can be utilized to determine the tractions on the boundary. Time-dependent partial differential equations are analyzed in the Laplace transformed domain and the Durbin’s inversion method is applied to determine the physical values in the time domain. Illustrative numerical examples are given and comparison has been made with the analytical solutions, the Boundary Element Method (BEM) and the Finite Element Method (FEM).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Effects of Dufour and fractional derivative on unsteady natural convection flow over an infinite vertical plate with constant heat and mass fluxes

In this paper, we analyze the effects of Dufour number and fractional-order derivative on unsteady natural convection flow of a viscous and incompressible fluid over an infinite vertical plate with constant heat and mass fluxes. The fractional constitutive model is obtained using fractional calculus approach. The Caputo fractional derivative operator is used in this problem. The dimensionless system of equations has been solved by employing Laplace transformation technique. Closed form solutions for concentration, temperature and velocity are presented in the form of Wright function and complementary error function. Effects of fractional and physical parameters on temperature and velocity profiles are illustrated graphically.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Numerical approach based on fractional-order Lagrange polynomials for solving a class of fractional differential equations

In this study, we propose a new set of fractional functions based on the Lagrange polynomials to solve a class of fractional differential equations. Fractional differential equations are the best tools for modelling natural phenomenon that are elaborated by fractional calculus. Therefore, we need an accurate and efficient technique for solving them. The main purpose of this article is to generalize new functions based on Lagrange polynomials to the fractional calculus. At first, we present a new representation of Lagrange polynomials and in continue, we propose a new set of fractional-order functions which are called fractional-order Lagrange polynomials (FLPs). Besides, a general formulation for operational matrices of fractional integration and derivative of FLPs on arbitrary nodal points are extracted. These matrices are obtained using Laplace transform. The initial value problems is reduced to the system of algebraic equations using the operational matrix of fractional integration and collocation method. Also, we find the upper bound of error vector for the fractional integration operational matrix and we indicate convergence of approximations of FLPs. Illustrative examples are included to demonstrate the validity and applicability of the proposed technique.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

E-Invariant Quantized Motion of Valence Quarks

In sub-proton space wave processes are impossible. The analog of the Klein–Gordon equation in sub-proton space is elliptical and describes a stationary system with a constant number of particles. For dynamical processes, separation of variables is used and in each quantum of motion of the quark two states are distinguished: a localization state and a translation state with infinite velocity. Alternation of these states describes the motion of a quark. The mathematical expectations of the lifetimes of the localization states and the spatial extents of the translation states for a free quark and for a quark in a centrally symmetric potential are found. The action after one quantum of motion is equal to the Planck constant. The one-sided Laplace transform is used to determine the Green’s function. Use of path integrals shows that the quantized trajectory of a quark is a broken line enveloping the classical trajectory of oscillation of the quark. Comparison of the calculated electric charge distribution in a proton with its experimental value gives satisfactory results. A hypothesis is formulated, according to which the three Grand Geometries of space correspond to the three main interactions of elementary particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Electron Spectrum of Nonlinear Cold Emission from a Metal under the Action of a Laser Shot

The nonlinear emission of electrons from a metal under the action of a femtosecond moderate-intensity laser pulse (laser shot) has been studied. A theoretical model of the process has been constructed based on the 1D nonstationary Schrödinger equation in the vacuum half-space with given boundary conditions for the electron wavefunction. This equation has been solved using the Laplace transformation. It has been assumed that the states of free electrons in a metal, which are described by the Sommerfeld theory of metals, are insignificantly influenced by the laser field. The energy spectrum of emitted electrons has been obtained, and its dependence on the parameters of the lased shot has been found. The calculated spectrum of nonlinear electron emission from a tungsten nanotip under the action of a 6.5-fs-long laser shot generating a field of 9.26 V/nm agrees with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Vacuum polarization in Coulomb field revisited

Simplified derivation of Wichmann–Kroll term is presented. The derivation uses two formulas for hypergeometric functions, but otherwise is elementary. It is found that Laplace transform of the vacuum charge density diverges at zero momentum transfer. This divergence has nothing to do with known ultraviolet divergence. The latter is related to the large momentum behavior of the pertinent integral, while the former to the small momentum behavior. When these divergences are removed, the energy shift caused by vacuum polarization for an ordinary hydrogen obtained here is in an exact agreement with the result obtained by Wichmann and Kroll. Also, for muonic hydrogen the result obtained here reasonably agrees with that given in literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2018

The Lorenzo–Hartley’s function for fractional calculus and its applications pertaining to fractional order modelling of anomalous relaxation in dielectrics

Considering the applications and a great significance of the study pertaining to generalized functions in applied sciences, the author investigated the Lorenzo–Hartley’s function which is a generalization of classical functions widely used in fractional calculus. The Laplace transform pairs are derived and the requirements for the Lorenzo–Hartley’s function G{sub ν,μ,δ}(a,c,t) to be completely monotonic (for t>0) are investigated. The author has shown the applications of this generalized function to describe the relaxation models, particularly in dielectrics. The Lorenzo–Hartley’s function G{sub ν,μ,δ}(a,c,t) of a real variable t is considered to investigate a computable mathematical framework for standard Debye and non-Debye relaxation processes in dielectric materials. The non-negative spectral distribution function is obtained for the corresponding response function. It is also demonstrated that the classical models like Cole–Cole (C–C), Davidson–Cole (D–C), and Havriliak–Negami (H–N) for non-Debye relaxation and standard Debye relaxation are the particular cases of the Lorenzo–Hartley’s function. Some of the study cases are also worked out to visualize the effects of variations of parameters on the response function and corresponding spectral distribution function. A generalized and unified fractional relaxation differential equation which governs response functions for classical dielectrics models pertaining to non-Debye relaxation (C–C, D–C, and H–N) is also established in the present investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗