Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maxwell 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Computational Modeling of Vertical-Cavity Surface-Emitting Lasers

Vertical-Cavity Surface-Emitting Lasers (VCSELs) have become increasingly important for many applications, such as optical interconnects and optical data storages. In many applications, critical requirements are very often high beam quality and/or high output powers. To help the design and optimization of VCSELs for such applications, it is essential to have a simulation model that includes the time evolution of the transverse space dependence of the hewn. Such models allow arbitrary transverse profiles of laser intensity or carrier density to develop in time without any a priori assumption about the number and types of transverse modes. In addition, extensive research in the past has shown that many-body effects are important in influencing the optical gain and the refractive index of semiconductor quantum wells. It is thus important to include the key many-body effects in simulation models for semiconductor lasers. A computational model, called the Maxwell-Effective Bloch Equations (MEBE) was developed recently that simulates the above features and effects. This model combines the space-time resolution with microscopic many-body effects in a way that is both accurate enough and computationally manageable. This model has been used to simulate edge-emitting high power lasers recently. A new finite-difference algorithm has been developed for solving these equations and will be shown in the presentation. That algorithm will be used to calculate the detailed time evolution of the transverse mode profiles of VCSELs. The VCSELs are based on InGaAs/GaAs quantum well structures. Effects of pumping current and refractive index profiles will be investigated in detail. The simulation results will be compared with other results the use the rate-equation model.

Goorjian, Peter M.↗
NASA NTRS1998

Simulation of Transverse Modes in Vertical-Cavity Surface-Emitting Lasers

An approximation to the Maxwell-Semiconductor Bloch equations is solved for the transverse modes of vertical-cavity surface-emitting lasers. A finite-difference algorithm will be used to calculate the detailed time evolution of the transverse mode profiles of VCSELs. Effects of pumping current and refractive index profiles will be investigated.

Goorjian, Peter M. P.↗
NASA NTRS1998

Computational Nonlinear of Vertical-Cavity Surface-Emitting Lasers

Vertical-Cavity Surface-Emitting Lasers (VCSELs) have become increasingly important for many applications, such as optical interconnects and optical data storages. In many applications, critical requirements are very often high beam quality and/or high output powers. To help the design and optimization of VCSELs for such applications, it is essential to have a simulation model that includes the time evolution of the transverse space dependence of the beam. Such models allow arbitrary transverse profiles of laser intensity or carrier density to develop in time without any a priori assumption about the number and types of transverse modes. In addition, extensive research in the past has shown that many-body effects are important in influencing the optical gain and the refractive index of semiconductor quantum wells. It is thus important to include the key many-body effects in simulation models for semiconductor lasers. A computational model, called the Maxwell-Effective Bloch Equations (MEBE) was developed recently that simulates the above features and effects. This model combines the space-time resolution with microscopic many-body effects in a way that is both accurate enough and computationally manageable. This model has been used to simulate edge-emitting high power lasers recently. A new finite-difference algorithm has been developed for solving these equations and will be shown in the presentation. That algorithm will be used to calculate the detailed time evolution of the transverse mode profiles of VCSELs. The VCSELs are based on InGaAs/GaAs quantum well structures. Effects of pumping current and refractive index profiles will be investigated in detail. The simulation results will be compared with other results the use the rate-equation model.

Goorjian, Peter M.↗
DOE OSTI2021

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗
DOE OSTI2022

VPIC-Kokkos

VPIC is a particle-in-cell (PIC) code that solves the Maxwell-Boltzmann system of equations and can be used to simulate a wide variety of plasmas, including those that are fully kinetic and relativistic. Utilizing the Kokkos performance-portable framework, VPIC achieves high performance on multiple CPU and GPU architectures and is adaptable to future platforms with minimal developer effort. VPIC features very powerful input decks, allowing insertion of arbitrary C++ code for custom diagnostics, boundary conditions, and additional physics models.

Luedtke, Scott↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electromagnetic Response of Superconducting RF Cavities

Recently reported anomalies in frequency shift of order δƒ ~ 0.1–10 kHz for Niobium SRF cavities in a narrow temperature region near T c ≃ 9 K are sensitive to the surface preparation and Nitrogen doping. We developed methods for calculating the surface current response of Nb SRF cavities, as well as the resonant frequency shift and quality factor, as functions of temperature, frequency and disorder based on the Keldysh formulation of the theory for nonequilibrium superconductivity coupled with Maxwell’s theory and boundary conditions for the response functions and electromagnetic field. We show that the anomaly is sensitive to impurity disorder. Our results for the anomaly in the frequency shift are in good agreement with experimental results for Nb with an anisotropic gap on the Fermi surface and inhomogeneous non-magnetic disorder. We also show that the quality factor as a function of the impurity scattering rate is maximum for in SRF cavities with intermediate disorder, ℏ/τΔ~O(1) with the maximum $Q$ decreasing with increasing frequency.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1981

The plasma wake of mesosonic conducting bodies. I - An experimental parametric study of ion focusing by the plasma sheath

The experimental investigation considered is concerned with the deflection of ion streams resulting from the interaction of conducting test bodies with an unmagnetized, mesosonic (supersonic with respect to ions but subsonic with respect to electrons) plasma stream. The investigation is, therefore, limited to plasma-electrostatic interactions. The experimental conditions are similar to those of the spacecraft-ionospheric interaction in that the ionic mass and number density, the electron temperature, and the plasma drift velocity ranges include values appropriate for small satellites or diagnostic probes at 200 to 400 km altitude. The study provides direct observations of deflected ion streams for cylindrical test bodies and gives a detailed description of the effects of the governing, dimensionless parameter ratios obtained from the steady-state, nonmagnetic Maxwell-Vlasov system of equations.

Stone, N. H.↗
NASA NTRS1993

Kinetic relaxation of a non-Maxwellian monatomic gas in a state of gross rest

The kinetic relaxation of a monatomic gas in a state of gross rest is studied by use of Maxwell's second-order moment equation. The molecular models considered range from the hard sphere to the Maxwell molecule. The development is exact for an ellipsoidal distribution and becomes a good approximation for more general distributions. The results show that the relaxation of second-order moments is nearly exponential for power-law molecules and that the characteristic time mu/p appears as the principal controlling parameter while the power-law constant plays a secondary role. Comparisons between results of numerical simulations using Bird's direct simulation Monte Carlo method and theory show excellent agreement for both hard-sphere and Maxwell molecules.

Baganoff, D.↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

The solution of the relaxation problem for the Boltzmann equation by the integral iteration method

The Boltzmann equation is considered in terms of the problem of relaxation of some initial distribution function which depends only on velocities, to Maxwell's distribution function. The Boltzmann equation is given for the relaxation problem in which the distribution function f(t, u, v) is time dependent and is also dependent on two other variables u and v (the velocities of rigid spherical molecules). An iteration process is discussed in which the velocity space u, v is subdivided into squares, the distribution function in each square being approximated by the second-order surface from the values of the distribution function at nine points. The set of all of these points forms a network of u, v values at the nodes of which the distribution function can be found.

Limar, Y. F.↗
NASA NTRS1965

Numerical integration of kinetic equations.

Relaxation to Maxwellian of Balescu-Lenard equation solved numerically as initial value problem for isotropic distributions of electrons in background of positive charge

MAXWELL DISTRIBUTION↗
DOE OSTI2020

Parameter space mapping of the Princeton magnetorotational instability experiment

Extensive simulations of the Princeton Magnetorotational Instability (MRI) Experiment with the Spectral/Finite Element code for Maxwell and Navier-Stokes Equations (SFEMaNS) have been performed to map the MRI-unstable region as a function of inner cylinder angular velocity and applied vertical magnetic field. The angular velocities of the outer cylinder and the end-cap rings follow the inner cylinder in fixed ratios optimized for MRI. Here, we first confirm the exponential growth of the MRI linear phase using idealized conducting vertical boundaries (end caps) rotating differentially with a Taylor-Couette profile. Subsequently, we run a multitude of simulations to scan the experimental parameter space and find that the normalized volume-averaged mean-square radial magnetic field, our main instability indicator, rises significantly where MRI is expected. At various locations, the local radial components of fluid velocity and generated magnetic field are well correlated with the volume-averaged indicator. Based on this correlation, a diagnostic system that will measure the radial magnetic field at several locations on the inner cylinder is proposed as the main comparison between simulation and experiment. A detailed analysis of poloidal mode structures in the SFEMaNS code indicates that MRI, rather than Ekman circulation or Rayleigh instability, dominates the fluid behavior in the region where MRI is expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Finite difference time domain modeling of steady state scattering from jet engines with moving turbine blades

The approach chosen to model steady state scattering from jet engines with moving turbine blades is based upon the Finite Difference Time Domain (FDTD) method. The FDTD method is a numerical electromagnetic program based upon the direct solution in the time domain of Maxwell's time dependent curl equations throughout a volume. One of the strengths of this method is the ability to model objects with complicated shape and/or material composition. General time domain functions may be used as source excitations. For example, a plane wave excitation may be specified as a pulse containing many frequencies and at any incidence angle to the scatterer. A best fit to the scatterer is accomplished using cubical cells in the standard cartesian implementation of the FDTD method. The material composition of the scatterer is determined by specifying its electrical properties at each cell on the scatterer. Thus, the FDTD method is a suitable choice for problems with complex geometries evaluated at multiple frequencies. It is assumed that the reader is familiar with the FDTD method.

Ryan, Deirdre A.↗