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

Patterned dielectric back contact design for GaAs thermophotovoltaic devices

-Patterned-dielectric back contact structures in optoelectronic devices are designed to boost the reflectance of light from the device back surface while retaining a low-resistance pathway for electrical conductance. Their reduced light absorption at near- and sub-bandgap photon energies leads to improved luminescence in light-emitting diodes, greater photon recycling, voltage, and efficiency in photovoltaic cells, and greater recuperation of unabsorbed sub-bandgap light in thermophotovoltaic (TPV) systems. However, diffraction from the patterned features can deflect incident light in propagation directions that lead to light trapping and parasitic absorption in the cell. Here, in this article, we use rigorous coupled-wave analysis (RCWA) to study three-dimensional diffractive scattering of electromagnetic waves by periodic metal point-contact gratings on 1.42-eV GaAs TPV cells, to analyze their effect on unwanted sub-bandgap absorption in order to achieve higher TPV system efficiency. Solutions of Maxwell's equations calculated using RCWA are compared to measured sub-bandgap reflectance in experimental GaAs TPV devices with varying metal point-contact diameters and spacing. Modeling and experiments indicate decreased total reflectance due to these diffractive effects for a small point contact diameter of 1 μm, and this effect is much stronger at higher contact coverage fractions.

14 SOLAR ENERGY↗
DOE OSTI2021

Scattering of radiofrequency waves by randomly modulated density interfaces in the edge of fusion plasmas

In the scrape-off layer and the edge region of a tokamak, the plasma is strongly turbulent and scatters the radiofrequency (RF) electromagnetic waves that propagate through this region. It is important to know the spectral properties of these scattered RF waves, whether used for diagnostics or for heating and current drive. The spectral changes influence the interpretation of the obtained diagnostic data, and the current and heating profiles. A full-wave, three-dimensional (3-D) electromagnetic code ScaRF has been developed for studying the RF wave propagation through turbulent plasma. ScaRF is a finite-difference frequency-domain (FDFD) method used for solving Maxwell's equations. The magnetized plasma is defined through the cold plasma by the anisotropic permittivity tensor. As a result, ScaRF can be used to study the scattering of any cold plasma RF wave. It can also be used for the study of the scattering of electron cyclotron waves in ITER-type and medium-sized tokamaks such as TCV, ASDEX-U and DIII-D. For the case of medium-sized tokamaks, there is experimental evidence that drift waves and rippling modes are present in the edge region. Hence, we have studied the scattering of RF waves by periodic density interfaces (plasma gratings) in the form of a superposition of spatial modes with varying periodicity and random amplitudes. The power reflection coefficient (a random variable) is calculated for different realizations of the density interface. In this work, the uncertainty of the power reflection coefficient is rigorously quantified by use of the Polynomial Chaos Expansion method in conjunction with the Smolyak sparse-grid integration, which is known as the PCE-SG method. The PCE-SG method is proven to be accurate and more efficient compared with alternative methods such as the Monte Carlo (MC) approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗
DOE OSTI2020

Quadrupole topological photonic crystals

Quadrupole topological phases, exhibiting protected boundary states that are themselves topological insulators of lower dimensions, have recently been of great interest. Extensions of these ideas from current tight binding models to continuum theories for realistic materials require the identification of quantized invariants describing the bulk quadrupole order. Here we identify the analog of quadrupole order in Maxwell’s equations for a gyromagnetic photonic crystal (PhC) through a double-band-inversion process. The quadrupole moment is quantized by the simultaneous presence of crystalline symmetry and broken time-reversal symmetry, which is confirmed using three independent methods: analysis of symmetry eigenvalues, numerical calculations of the nested Wannier bands and the expectation value of the quadrupole operator. Furthermore, we reveal the boundary manifestations of quadrupole phases as quantized edge polarizations and fractional corner charges. The latter are the consequence of a filling anomaly of energy bands as first predicted in electronic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Optical analogues to the equatorial Kerr–Newman black hole

Optical analogues to black holes allow the investigation of general relativity in a laboratory setting. Previous works have considered analogues to Schwarzschild black holes in an isotropic coordinate system; the major drawback is that required material properties diverge at the horizon. We present the dielectric permittivity and permeability tensors that exactly reproduce the equatorial Kerr–Newman metric, as well as the gradient-index material that reproduces equatorial Kerr–Newman null geodesics. Importantly, the radial profile of the scalar refractive index is finite along all trajectories except at the point of rotation reversal for counter-rotating geodesics. Construction of these analogues is feasible with available ordinary materials. A finite-difference frequency-domain solver of Maxwell’s equations is used to simulate light trajectories around a variety of Kerr–Newman black holes. For reasonably sized experimental systems, ray tracing confirms that null geodesics can be well-approximated in the lab, even when allowing for imperfect construction and experimental error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas

A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.

Solmaz, Evrim (ORCID:0000000211830881)↗
DOE OSTI2021

Simulations of wavelength-multiplexed holography for single-shot spatiotemporal characterization of NIF’s advanced radiographic capability (ARC) laser

We simulate the use of a newly developed single-shot wavelength-multiplexed holography-based diagnostic, STRIPED FISH, to fully characterize the as-delivered laser pulses of the National Ignition Facility’s Advanced Radiographic Capability (NIF-ARC) laser. To that end, we have performed simulations of the NIF-ARC pulse incorporating (a) a time-integrated spatial-profile measurement and a complete temporal-intensity-and-phase measurement using a frequency resolved optical gating, but without any spatiotemporal pulse characterizations, and (b) simulated first-order spatiotemporal distortions, which could be measured on a single shot if a STRIPED FISH device were deployed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Effects of Coulomb collisions on lower hybrid drift waves inside a laboratory reconnection current sheet

We have developed a local, linear theoretical model for lower hybrid drift waves that can be used for plasmas in the weakly collisional regime. Two cases with typical plasma and field parameters for the current sheet of the magnetic reconnection experiment have been studied. For a case with a low electron beta (β e =0.25, high guide field case), the quasi-electrostatic lower hybrid drift wave is unstable, while the electromagnetic lower hybrid drift wave has a positive growth rate for a high-β e case (β e =8.9, low guide field case). For both cases, including the effects of Coulomb collisions reduces the growth rate but collisional impacts on the dispersion and growth rate are limited (≲20%).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Using extended MHD to explore lasers as a trigger for x-pinches

X-pinches have been shown to be a source of extremely intense x-ray emissions useful for diagnosing plasma dynamics and imaging biological objects. The most striking feature of an x-pinch is the hotspot, the point source from where all the x rays come from. Unfortunately, the exact timing and location of the hotspot are still unpredictable. Since an x-pinch hotspot forms from instabilities (like an m = 0 mode), we will computationally explore whether we can use a high-power laser to control the timing and location of these instabilities. Our goal is to reduce the temporal and spatial jitter associated with the x-ray burst. Using an extended MHD (XMHD) code, we explore the non-relativistic instability generation using a current profile of a 250 kA linear transformer driver system and laser characteristics of SLAC's Matter in Extreme Conditions laboratory. Our results include both laser-penetration results using a boundary-defined electromagnetic (EM) wave and instability results from a power-deposition method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The physical meaning of the magnetic scalar potential and its use in the design of hermetic electromagnetic coils

The magnetic scalar potential U(r⃗) has an intuitive physical interpretation, representing the electric current I = ΔU, which must be directed between any pair of isocontours along the boundary of a region, differing in potential by ΔU, in order to generate the corresponding magnetic intensity H⃗=−∇⃗U inside the region with no tangential component outside the boundary. This physical significance is exploited to invert the design process of hermetic (fringeless) electromagnetic coils from standard iterative techniques (calculating the magnetic field of refined current distributions) to a procedure for calculating the exact surface currents required to generate a given field configuration. A practical construction algorithm is given to produce the prescribed field in a “Target” region (constrained only by Maxwell’s equations) while confining the fringe field to a specified hermetic “Return” region. Example coils are analyzed along with a discussion of the limits of precision of the constructed field.

Crawford, C. B. (ORCID:0000000219324334)↗
DOE OSTI2021

Improved prescription for winding an electromagnet

We describe an improvement on the magnetic scalar potential approach to design an electromagnet, which incorporates the need to wind the coil as a helix. Any magnetic field that can be described by a magnetic scalar potential is produced with high fidelity within a Target region; all fields are confined within a larger Return. The helical winding only affects the field in the Return.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Driving rotamak currents with minimal power dissipation

The Rotamak is a proposed thermonuclear fusion device which employs rotating magnetic fields (RMF) to generate an azimuthal current to produce a field-reversed configuration. The efficiency of the currents that produce the field reversal by RMFs was debated some 40 years ago. The debate revolved around whether the currents would incur dissipation by the conventional Spitzer perpendicular resistivity, or whether some other relation between current and dissipation would be more appropriate. By employing an electron–ion pitch-angle scattering model, we find that the dissipation is non-Spitzer in nature. However, curiously, there appears to exist a regime where the power dissipated to maintain the current becomes vanishingly small.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Local analysis of fast magnetic reconnection

Fast magnetic reconnection is defined by the topology of the magnetic field lines changing on a timescale that is approximately an order of magnitude longer than the topology-conserving ideal-evolution timescale. Fast reconnection is an intrinsic property of Faraday's law when the evolving magnetic field depends non-trivially on all three spatial coordinates and is commonly observed—even when the effects that allow topology breaking are arbitrarily small. The associated current density need only be enhanced by a factor of approximately ten and flows in thin but broad ribbons along the magnetic field. These results follow from the variation in the separation of neighboring pairs of magnetic field lines, which in an ideal evolution typically increases exponentially with time, and the existence of a spatial scale below which magnetic field lines freely change their identities due to non-ideal effects, such as resistivity. Traditional reconnection theory ignores exponentially large variations and relies on the current density reaching a magnitude that is exponentially larger than is actually required. Here, an analysis of the behavior of magnetic field lines in the neighborhood of an arbitrarily chosen line is used to obtain more precise and rigorous results on intrinsic reconnection. The maximum parallel kinetic energy of collisionless charged particles is shown to have an exponential increase in time during a generic magnetic evolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Metriplectic foundations of gyrokinetic Vlasov–Maxwell–Landau theory

Here, this Letter reports on a metriplectic formulation of a collisional, nonlinear full-f electromagnetic gyrokinetic theory compliant with energy conservation and monotonic entropy production. In an axisymmetric background magnetic field, the toroidal angular momentum is also conserved. Notably, a new collisional current, contributing to the gyrokinetic Maxwell–Ampère equation and the gyrokinetic charge conservation law, is discovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Dissociation slowdown by collective optical response under strong coupling conditions

We consider an ensemble of diatomic molecules resonantly coupled to an optical cavity under strong coupling conditions at normal incidence. Photodissociation dynamics is examined via direct numerical integration of the coupled Maxwell–Schrödinger equations with molecular rovibrational degrees of freedom explicitly taken into account. It is shown that the dissociation is significantly affected (slowed down) when the system is driven at its polaritonic frequencies. The observed effect is demonstrated to be of transient nature and has no classical analog. An intuitive explanation of the dissociation slowdown at polaritonic frequencies is proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗