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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the nature of two-photon transitions for a collection of molecules in a Fabry–Perot cavity

Here, we investigate the effect of a cavity on nonlinear two-photon transitions of a molecular system and we analyze how such an effect depends on the cavity quality factor, the field enhancement, and the possibility of dephasing. We find that the molecular response to strong light fields in a cavity with a variable quality factor can be understood as arising from a balance between (i) the ability of the cavity to enhance the field of an external probe and promote multiphoton transitions more easily and (ii) the fact that the strict selection rules on multiphoton transitions in a cavity support only one resonant frequency within the excitation range. Although our simulations use a classical level description of the radiation field (i.e., we solve Maxwell–Bloch or Maxwell–Liouville equations within the Ehrenfest approximation for the field–molecule interaction), based on experience with this level of approximation in the past studies of plasmonic and polaritonic systems, we believe that our results are valid over a wide range of external probing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electromagnetic formalism

This is a companion paper to the previous work [P. H. Yoon, Phys. Plasmas 31, 032309 (2024)] in which the nonlinear susceptibilities of weakly turbulent magnetized plasma are derived under a simplifying assumption of electrostatic interaction. The present paper extends the analysis to a general situation of electromagnetic interaction. The main novelty of the previous and present papers is that by employing the Bessel function addition theorem, the mathematical definitions for the susceptibilities are substantially simplified, a procedure that has not been discussed in the existing literature. In the present paper, a full set of Maxwell’s equations are considered in conjunction with the nonlinear Vlasov equation, which is solved by a perturbative method. The result is a fully general nonlinear susceptibility, given in tensorial form, which is applicable for weakly turbulent magnetized plasmas.

Materials Science↗
DOE OSTI2024

Toward an integrated multi-gigahertz ionizing particle diagnostic

Needs arising at both current and future accelerator facilities call for the development of radiation-hardened position-sensing diagnostics that can operate with multi-GHz repetition rates. Such instruments are likely to also have applications in the diagnosis of rapid plasma behavior. Here, building on the recent work of our Advanced Accelerator Diagnostics Collaboration, we are exploring the development of integrated multi-GHz ionizing particle detection systems based on chemical-vapor deposition diamond sensors, with the initial goal of producing a quadrant detector that can determine the intensity and centroid position of a particle beam at a repetition rate between 5 and 10 GHz. Results from our initial high-speed characterization work are presented, including those from a single-channel sensor with a GHz response. Approaches to achieving multi-GHz (5–10 GHz) rate capability, including the design of a dedicated Application Specific Integrated Circuit and the use of 3D RF-solver computer aided design software, are presented and discussed in more detail. 3D RF simulations suggest clean pulses of duration less than 250 ps (FWHM less than 125 ps) can be achieved with the approaches developed by this work.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dependent scattering and fractal microstructure determine the transparency of aerogel monoliths

This study reveals how dependent scattering and microstructure significantly affect electromagnetic wave propagation through aerogel monoliths, contributing to their transparency. Light scattering by particle ensembles is considered “dependent” when the scattering properties rely not only on particle size and optical constants but also on their spatial distribution, typically occurring when the average interparticle distance is small in comparison with the wavelength of incident radiation. Addressing dependent scattering requires solving Maxwell’s equations for complex heterogeneous structures, which is computationally demanding and usually limited to sample thicknesses on the same scale as the wavelength. This study combines computer-generated ambigel microstructures of fractal aggregates of polydisperse nanoparticles and the radiative transfer with reciprocal transaction method to predict the transmittance of thick ambigel slabs. Transmittance measurements of ambiently dried aerogel monoliths (ambigels) with porosities from about 50% to 90% closely matched the predicted values for their digital twins. However, ignoring dependent scattering or particle aggregation led to inaccurate predictions. This study validated the computational framework, and its findings offer insights for designing photonic metamaterials and analyzing their interactions with electromagnetic waves.

Yalcin, Refet A. (ORCID:0000000339973494)↗
DOE OSTI2025

Algebra of invariants for the Vlasov–Maxwell system

The algebra of invariants for both the relativistic and nonrelativistic multispecies Vlasov–Maxwell system is examined, including the case with a fixed ion background. Invariants and their associated fluxes are obtained directly from the Vlasov–Maxwell system. The invariants are shown to Poisson commute with the Hamiltonian and the rest of the Poisson bracket algebra of invariants is identified. Special attention is given to the role played by the monopole condition, ∇ · B.

Fundamental invariants↗
DOE OSTI2026

Gyrotropy-induced symmetry breaking in RF-specific plasma wall interaction: Mechanisms and mitigation strategies

Ion cyclotron range of frequencies (ICRF) heating is essential for creating plasma in next-generation fusion devices. ICRF antennas often produce localized wall erosion, reducing reliability, material survivability, and overall plasma performance. We show that the wave physics' gyrotropy is among the causes of toroidal and poloidal asymmetries in these erosion spots, and that this source of asymmetry can be compensated for either by controlling the poloidal phasing or by modifying the limiter shapes.

Fusion reactors↗
DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Unipolar light: existence, generation, propagation, and impact on microobjects

A review is presented of recent works on optical unipolar pulses, whose electric area (integral of the electric field strength over time) is nonzero, which determines the predominant direction of the electric field strength. It is shown that the existence of unipolar pulses does not contradict Maxwell’s equations, and that unipolar pulses can propagate in waveguides. It is emphasised that, along with the spectral, energy, and polarisation parameters, the electric area of short light pulses is also an important characteristic. The unidirectional action of these pulses on microobjects indicates that it is promising to develop methods for generating radiation with such properties. We disciuss methods for the generation, propagation, and interaction of unipolar light with classical and quantum systems, as well as methods for recording the electric area of light pulses. (invited paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Roadmap on spatiotemporal light fields

Spatiotemporal sculpturing of light pulse with ultimately sophisticated structures represents a major goal of the everlasting pursue of ultra-fast information transmission and processing as well as ultra-intense energy concentration and extraction. It also holds the key to unlock new extraordinary fundamental physical effects. Traditionally, spatiotemporal light pulses are always treated as spatiotemporally separable wave packet as solution of the Maxwell’s equations. In the past decade, however, more generalized forms of spatiotemporally nonseparable solution started to emerge with growing importance for their striking physical effects. This roadmap intends to highlight the recent advances in the creation and control of increasingly complex spatiotemporally sculptured pulses, from spatiotemporally separable to complex nonseparable states, with diverse geometric and topological structures, presenting a bird’s eye viewpoint on the zoology of spatiotemporal light fields and the outlook of future trends and open challenges.

42 ENGINEERING↗
DOE OSTI2023

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

ABSTRACT Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell’s equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly rotating magnetized neutron stars, taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. Rs ∼ 0.5 [R*], we show that the new terms lead to a per cent order correction in the magnetic field orientation and strength compared to the case, with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable, when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

Torres, R. (ORCID:0000000291820228)↗
DOE OSTI2024

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Charge Conservation beyond Uniformity: Spatially Inhomogeneous Electromagnetic Response in Periodic Solids

Nonlinear electromagnetic response functions have reemerged as a crucial tool for studying quantum materials, due to recently appreciated connections between optical response functions, quantum geometry, and band topology. Most attention has been paid to responses to spatially uniform electric fields, relevant to low-energy optical experiments in conventional solid state materials. However, magnetic and magnetoelectric phenomena are naturally connected by responses to spatially varying electric fields due to Maxwell’s equations. Furthermore, in the emerging field of moiré materials, characteristic lattice scales are much longer, allowing spatial variation of optical electric fields to potentially have a measurable effect in experiments. In order to address these issues, we develop a formalism for computing linear and nonlinear responses to spatially inhomogeneous electromagnetic fields. Starting with the continuity equation, we derive an expression for the second-quantized current operator that is manifestly conserved and model independent. Crucially, our formalism makes no assumptions on the form of the microscopic Hamiltonian and so is applicable to model Hamiltonians derived from tight-binding or calculations. We then develop a diagrammatic Kubo formalism for computing the wave vector dependence of linear and nonlinear conductivities, using Ward identities to fix the value of the diamagnetic current order by order in the vector potential. We apply our formula to compute the magnitude of the Kerr effect at oblique incidence for a model of a moiré-Chern insulator and demonstrate the experimental relevance of spatially inhomogeneous fields in these systems. We further show how our formalism allows us to compute the (orbital) magnetic multipole moments and magnetic susceptibilities in insulators. Turning to nonlinear response, we use our formalism to compute the second-order transverse response to spatially varying transverse electric fields in our moiré-Chern insulator model, with an eye toward the next generation of experiments in these systems. Published by the American Physical Society 2024

Physics↗
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 OSTI2022

A massively parallel time-domain coupled electrodynamics–micromagnetics solver

We present a high-performance coupled electrodynamics–micromagnetics solver for full physical modeling of signals in microelectronic circuitry. The overall strategy couples a finite-difference time-domain approach for Maxwell’s equations to a magnetization model described by the Landau–Lifshitz–Gilbert equation. The algorithm is implemented in the Exascale Computing Project software framework, AMReX, which provides effective scalability on manycore and GPU-based supercomputing architectures. Furthermore, the code leverages ongoing developments of the Exascale Application Code, WarpX, which is primarily being developed for plasma wakefield accelerator modeling. Our temporal coupling scheme provides second-order accuracy in space and time by combining the integration steps for the magnetic field and magnetization into an iterative sub-step that includes a trapezoidal temporal discretization for the magnetization. The performance of the algorithm is demonstrated by the excellent scaling results on NERSC multicore and GPU systems, with a significant (59×) speedup on the GPU using a node-by-node comparison. We demonstrate the utility of our code by performing simulations of an electromagnetic waveguide and a magnetically tunable filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Hydrodynamics of higher-rank gauge theories

We extend recent work on hydrodynamics with global multipolar symmetries — known as “fracton hydrodynamics” — to systems in which the multipolar symmetries are gauged. We refer to the latter as “fracton magnetohydrodynamics”, in analogy to conventional magnetohydrodynamics (MHD), which governs systems with gauged charge conservation. We show that fracton MHD arises naturally from higher-rank Maxwell’s equations and in systems with one-form symmetries obeying certain constraints; while we focus on “minimal” higher-rank generalizations of MHD that realize diffusion, our methods may also be used to identify other, more exotic hydrodynamic theories (e.g., with magnetic subdiffusion). In contrast to semi-microscopic derivations of MHD, our approach elucidates the origin of the hydrodynamic modes by identifying the corresponding higher-form symmetries. Being rooted in symmetries, the hydrodynamic modes may persist even when the semi-microscopic equations no longer provide an accurate description of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A Virtual Element Method for the Full MHD system

In this manuscript we present a novel discretization for the incompressible MHD system. Our approach follows the framework of the Virtual Element Method and offers two main advantages. The method can be implemented in unstructured meshes making it highly versatile and capable of handling a wide array of problems involving interfaces, free-boundaries or adaptive refinements on the mesh. The second advantage involves the divergence of the magnetic field, our approach guarantees that it remains solenoidal. We include a theoretical proof of the condition on the magnetic field as well as energy estimates and a well-posedness study. The latter sheds light as to the stability properties of the method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗