Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIRAC EQUATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

The 3 K blackbody radiation, Dirac's Large Numbers Hypothesis, and scale-covariant cosmology

A program is described which is intended to derive a generalized system of gravitational equations that allow (but do not require) G to vary, to use the 3-K blackbody radiation to fix the relation between G and the gauge function, and to employ Dirac's (1937) Large Numbers Hypothesis to derive the geometry of the universe. Einstein's equations are retained in their total integrity, but the specification is made that they are valid only when gravitational units are used. A scale-invariant form of Einstein's equations is obtained, and from this are derived the energy conservation law, the baryon-number conservation law, and the appropriate cosmological equations. Dirac's proposals of 1937 and 1973 are incorporated into the formalism, and a gauge based on consolidation of the 3-K blackbody radiation is presented. A unique solution for the geometry of the universe is determined for zero curvature solely from the 3-K radiation and the Large Numbers Hypothesis; this solution predicts a deceleration parameter exactly equal to unity.

Canuto, V.↗
NASA NTRS1995

Pseudoclassical description of the Dirac Oscillator

In this paper we discuss the Dirac Oscillator wave equation in terms of pseudoclassical language, using Grassmann variables to describe the internal degrees of freedom of the oscillator. Regarding the original wave equation as a classical constraint, we use the theory of constrained systems, to develop a reparameterization invariant lagrangian, which is the pseudoclassical equivalent of the quantum case. The consistency of the Hamiltonian formalism and the quantization procedure are also analyzed.

Delsolmesa, Antonio↗
NASA NTRS2007

A Fast Method of Deriving the Kirchhoff Formula for Moving Surfaces

The Kirchhoff formula for a moving surface is very useful in many wave propagation problems, particularly in the prediction of noise from rotating machinery. Several publications in the last two decades have presented derivations of the Kirchhoff formula for moving surfaces in both time and frequency domains. Here we present a method originally developed by Farassat and Myers in time domain that is both simple and direct. It is based on generalized function theory and the useful concept of imbedding the problem in the unbounded three-dimensional space. We derive an inhomogeneous wave equation with the source terms that involve Dirac delta functions with their supports on the moving data surface. This wave equation is then solved using the simple free space Green's function of the wave equation resulting in the Kirchhoff formula. The algebraic manipulations are minimal and simple. We do not need the Green's theorem in four dimensions and there is no ambiguity in the interpretation of any terms in the final formulas. Furthermore, this method also gives the simplest derivation of the classical Kirchhoff formula which has a fairly lengthy derivation in physics and applied mathematics books. The Farassat-Myers method can be used easily in frequency domain.

Farassat, F.↗
DOE OSTI2024

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗
NASA NTRS1992

Wave diffraction in weak cosmic-ray-modified shocks

Weakly multidirectional, long-wavelength cosmic-ray-modified shocks are studied via multiple scale perturbation techniques. The effects of diffraction are discussed in terms of Green's function solutions of the linearized 1 + 3D Burgers and 1 + 3D KdVB equations, and also in terms of solutions with singular Dirac delta initial distributions. The solutions show a monotonic decrease of the wave-front curvature with increasing time owing to the effects of wave diffraction. The shape of the wave surface is discussed in terms of solutions S to the wave eikonal equation corresponding to singular initial conditions. For the fast magnetosonic wave propagating in the positive x-direction, the wave phase surface S = 0 has elliptic cross sections with the planes x = constant and has a convex paraboloidal shape. Plane-wave solutions of the 1 + 3D KdVB equation are discussed.

Webb, G. M.↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

The Duffin-Kemmer-Petiau oscillator

In view of current interest in relativistic spin-one systems and the recent work on the Dirac Oscillator, we introduce the Duffin-Kemmer-Petiau (DKP) equation obtained by using an external potential linear in r. Since, in the non-relativistic limit, the spin 1 representation leads to a harmonic oscillator with a spin-orbit coupling of the Thomas form, we call the equation the DKP oscillator. This oscillator is a relativistic generalization of the quantum harmonic oscillator for scalar and vector bosons. We show that it conserves total angular momentum and that it is exactly solvable. We calculate and discuss the eigenspectrum of the DKP oscillator in the spin 1 representation.

Nedjadi, Youcef↗
DOE OSTI2025

Relativistic spin precession in homogeneous background fields

The Bargmann-Michel-Telegdi equation, which describes the precession of the spin of a charged Dirac particle moving in a homogeneous electromagnetic field, is generalized to include also other homogeneous background fields. The treatment incorporates observable coefficients that govern operators of mass dimensions three through six in the underlying Dirac effective field theory. A relativistic formulation valid in arbitrary inertial frames is obtained. The results are applicable to searches for new physics beyond the Standard Model, including searches for Lorentz and CPT violation.

Ding, Yunhua [Ohio Wesleyan Univ., Delaware, OH (U↗
NASA NTRS1992

Error analysis of finite element results on plates with nonuniform grids

Discrete Dirac delta functions in two dimensions and the governing plate differential equations are used to produce continuous approximations from discrete finite-element data on nonuniform grids. The resulting continuous approximation when compared to the discrete data provides an error analysis. Moreover, the continuous approximation can be differentiated to compute continuous stresses needed for component failure prediction. Stresses can then be calculated accurately at boundaries and at interfaces where finite-element analysis alone may not provide reliable results. The continuous solution along with the differential equations also provides a residual error which is used in the error analysis. This is a unique feature of this approach. Two application examples of the procedure are provided. One example involves the linear response of a clamped square plate under uniform transverse load. Results from finite-element analysis and the continuous solution from the present analysis are compared to a closed-form elasticity solution. The continuous solution is used in a Zienkiewicz-Zhu error estimator to assess the performance of the present analysis. The continuous solution is observed to dramatically correct the finite-element results. In a second example, results for the nonlinear response of a stiffened panel are compared to the finite-element solution.

Sistla, Rajaram↗
NASA NTRS2012

Positron Interactions with Atoms and Ions

Dirac, in 1928, combining the ideas of quantum mechanics and the ideas of relativity invented the well-known relativistic wave equation. In his formulation, he predicted an antiparticle of the electron of spin n-bar/2. He thought that this particle must be a proton. Dirac published his interpretation in a paper 'A theory of electrons and protons.' It was shown later by the mathematician Hermann Weyl that the Dirac theory was completely symmetric between negative and positive particles and the positive particle must have the same mass as that of the electron. In his J. Robert Oppenheimer Memorial Prize Acceptance Speech, Dirac notes that 'Blackett was really the first person to obtain hard evidence for the existence of a positron but he was afraid to publish it. He wanted confirmation, he was really over cautious.' Positron, produced by the collision of cosmic rays in a cloud chamber, was detected experimentally by Anderson in 1932. His paper was published in Physical Review in 1933. The concept of the positron and its detection were the important discoveries of the 20th century. I have tried to discuss various processes involving interactions of positrons with atoms and ions. This includes scattering, bound states and resonances. It has not been possible to include the enormous work which has been carried out during the last 40 or 50 years in theory and measurements.

Bhatia, Anand K.↗
NASA NTRS1983

Transient difference solutions of the inhomogeneous wave equation: Simulation of the Green's function

A time-dependent finite difference formulation to the inhomogeneous wave equation is derived for plane wave propagation with harmonic noise sources. The difference equation and boundary conditions are developed along with the techniques to simulate the Dirac delta function associated with a concentrated noise source. Example calculations are presented for the Green's function and distributed noise sources. For the example considered, the desired Fourier transformed acoustic pressures are determined from the transient pressures by use of a ramping function and an integration technique, both of which eliminates the nonharmonic pressure associated with the initial transient.

Baumeiste, K. J.↗
NASA NTRS1983

Transient difference solutions of the inhomogeneous wave equation - Simulation of the Green's function

A time-dependent finite difference formulation to the inhomogeneous wave equation is derived for plane wave propagation with harmonic noise sources. The difference equation and boundary conditions are developed along with the techniques to simulate the Dirac delta function associated with a concentrated noise source. Example calculations are presented for the Green's function and distributed noise sources. For the example considered, the desired Fourier transformed acoustic pressures are determined from the transient pressures by use of a ramping function and an integration technique, both of which eliminates the nonharmonic pressure associated with the initial transient.

Baumeister, K. J.↗
NASA NTRS1989

Surface structure of neutron stars with high magnetic fields

The equation of state of cold dense matter in strong magnetic fields is calculated in the Thomas-Fermi and Thomas-Fermi-Dirac approximations. For use in the latter calculation, a new expression is derived for the exchange energy of the uniform electron gas in a strong magnetic field. Detailed calculations of the density profile in the surface region of a neutron star are described for a variety of equations of state, and these show that the surface density profile is strongly affected by the magnetic field, irrespective of whether or not matter in a magnetic field has a condensed state bound with respect to isolated atoms. It is also shown that, as a consequence of the field dependence of the screening potential, magnetic fields can significantly increase nuclear reaction rates.

Fushiki, I.↗
NASA NTRS2008

Nonstandard Analysis and Jump Conditions for Converging Shock Waves

Nonstandard analysis is an area of modern mathematics which studies abstract number systems containing both infinitesimal and infinite numbers. This article applies nonstandard analysis to derive jump conditions for one-dimensional, converging shock waves in a compressible, inviscid, perfect gas. It is assumed that the shock thickness occurs on an infinitesimal interval and the jump functions in the thermodynamic and fluid dynamic parameters occur smoothly across this interval. Predistributions of the Heaviside function and the Dirac delta measure are introduced to model the flow parameters across a shock wave. The equations of motion expressed in nonconservative form are then applied to derive unambiguous relationships between the jump functions for the flow parameters.

Baty, Roy S.↗
NASA NTRS2000

Modeling of Optoelectronic Devices

Ultrafast modulation of semiconductor quantum well (QW) laser is of technological importance for information technology. Improvement by order(s) of magnitude in data transfer rate is possible as terahertz (THz) radiation is available for heating the laser at picosecond time scale. Optical gain modulation in the QW is achieved via temperature modulation of electron-hole plasma (EHP). Applications include free-space THz communication, optical switching, and pulse generation. The EHP in the semiconductor QW is described with a two-band model. Semiconductor Bloch equations with many-body effects are used to derive a hydrodynamical model for the active QW region. Because of ultrafast carrier-carrier scatterings in the order of 50 fs, EHP follows quasiequilibrium Fermi-Dirac distributions and THz field interacts incoherently with it. Carrier-longitudinal optical (LO) phonon scatterings and coherent laser-EHP interaction are treated microscopically in our physical model. A set of hydrodynamical equations for plasma density, temperature, and laser envelop amplitude are derived and Runge-Kutta method is adopted for numerical simulation. A typical 8 nm GaAs/Al(0.3)Ga(0.7) As single QW at 300 K is used. Additional information is contained in the original extended abstract.

Li, Jian-Zhong↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effect of interparticle fields and radiation reaction on beam dynamics

The dynamics of relativistic particles in an intense electromagnetic field can be described by the Landau-Lifshitz (LL) equation, where the radiation reaction (RR) is accounted for via a self-force, and interparticle fields are often neglected as an approximation. However, the inclusion of interparticle fields is necessary to ensure energy-momentum conservation, particularly during coherent emission. Here we present (i) an analytical proof showing that the energy-momentum conservation law of the Hamilton-Rohrlich-Dirac action, which is divergence free and describes a generic system of interacting charges, respects causality and provides physically sensible results; (ii) a simple generalization of the LL equation for many particles evaluated as a function of the total field, i.e., the sum of the external and interparticle fields. By performing first-principles numerical simulations of a neutral, relativistic bunch of electrons and positrons (e − /e + ) colliding with a laser pulse, this theory is shown to satisfy energy-momentum conservation when interparticle fields and RR are simultaneously taken into account; and (iii) the combined effect of interparticle fields and RR primarily affects the tail of the particle energy distribution. Additionally, our first-principles simulations show that the effect of interparticle fields on beam energy loss becomes smaller when most of the radiated energy is incoherent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗