Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maxwell”

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 415 records · Page 23

Novel fractal-textured solar absorber surfaces for concentrated solar power

The solar absorptance of a thermal receiver surface significantly affects the photothermal conversion efficiency of concentrated solar power (CSP) plants. The development of low-cost Gen3 CSP systems calls for increasing the solar absorptance of the thermal receivers at operating temperatures above 750 °C. This study presents an innovative approach to the fractal, multiscale texturing of absorber surfaces to significantly enhance solar energy absorption for Gen3 CSP applications. The absorber surface is described in terms of its fractal parameters that are uniquely determined from surface profile measurements. The interaction of solar radiation with the fractal surface is numerically simulated by solving the governing Maxwell’s equations for electromagnetic wave propagation to investigate the effect of texturing on the absorptance of the surfaces. It is shown that absorption of the solar spectrum increases with an increase in the fractal dimension and the multiscale asperity height of the surface texturing. Here, the theoretical model is demonstrated to be in close agreement with experimental measurements of spectral absorptance of electrodeposited copper (Cu), copper mono-oxide (CuO), and copper-manganese oxide (CuMnO) surfaces that are textured to produce a range of fractal parameters by tailoring the deposition parameters. Fractal surface texturing is shown to reduce reflectance by over one order of magnitude, yielding an absorptance of greater than 0.985 for CuMnO. For the first time, the study presents a significant effective means of increasing solar absorptance and a fundamental theoretical description of the underlying physics.

14 SOLAR ENERGY↗

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↗

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Interaction of radio frequency waves with cylindrical density filaments: scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas , vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci. , vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

Physics↗

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)↗

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗

Theoretical Pathway toward Improved Reverse Osmosis Membrane Selectivity for Neutral Solutes: Inspiration from Gas Separations

We report highly rigid membrane materials with tailored structures have exhibited permeabilities and selectivities that exceed the polymer upper bound in gaseous and organic solvent separations for challenging mixtures containing species that are similar in size and shape. One potential question is whether such membrane materials can maintain meaningful guest diffusivities in situations where the microporous spaces are essentially full of guest molecules. Here, we use a simplified transition state theory approach to estimate the diffusivity of water and small organics within a microporous membrane. The transition state theory model is parameterized using experimental values from zeolites and carbon molecular sieve materials found in the literature. We demonstrate the differences in transport and Maxwell–Stefan diffusivities based on guest species loading with different isotherm behaviors. These calculations theorize a path forward for highly selective reverse osmosis membranes for aqueous phase separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Harmonic and Subharmonic RF Injection Locking of THz Metasurface Quantum-Cascade VECSEL

Harmonic and subharmonic RF injection locking is demonstrated in a terahertz (THz) quantum-cascade vertical-external-cavity surface-emitting laser (QC-VECSEL). By tuning the RF injection frequency around integer multiples and submultiples of the cavity round-trip frequency, different harmonic and subharmonic orders can be excited in the same device. Modulation-dependent behavior of the device has been studied with recorded lasing spectral broadening and locking bandwidths in each case. In particular, harmonic injection locking results in the observation of harmonic spectra with bandwidths over 200 GHz. A semiclassical Maxwell-density matrix formalism has been applied to interpret QC-VECSEL dynamics, which aligns well with experimental observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Substituent Effects in Mechanochemical Allowed and Forbidden Cyclobutene Ring-Opening Reactions

Woodward and Hoffman once jested that a very powerful Maxwell demon could seize a molecule of cyclobutene at its methylene groups and tear it open in a disrotatory fashion to obtain butadiene (Woodward, R. B.; Hoffmann, R. The Conservation of Orbital Symmetry. Angew. Chem., Int. Ed. 1969, 8, 781–853). Nearly 40 years later, that demon was discovered, and the field of covalent polymer mechanochemistry was born. In the decade since our demon was befriended, many fundamental investigations have been undertaken to build up our understanding of force-modified pathways for electrocyclic ring-opening reactions. In this work, we seek to extend that fundamental understanding by exploring substituent effects in allowed and forbidden ring-opening reactions of cyclobutene (CBE) and benzocyclobutene (BCB) using a combination of single-molecule force spectroscopy (SMFS) and computation. We show that, while the forbidden ring-opening of cis-BCB occurs at a lower force than the allowed ring-opening of trans-BCB on the time scale of the SMFS experiment, the opposite is true for cis- and trans-CBE. Such a reactivity flip is explained through computational analysis and discussion of the so-called allowed/forbidden gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields

As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

54 ENVIRONMENTAL SCIENCES↗

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strong Sensitivity of Simulated Biomass Burning Aerosol Transport and Radiative Effects Over the South Atlantic to Carbonaceous Aerosol Aging and Particle Density

Biomass burning aerosol (BBA) impacts climate through aerosol‐cloud‐radiation interactions, but models disagree on the sign and magnitude of BBA radiative effects. We quantify the sensitivity of BBA radiative effects and transport to three BBA‐relevant processes and properties: parameterized oxidative aging of organic aerosol (OA), a combined change to black carbon (BC) density and the method for calculating aerosol refractive index, and reduction in OA density. We evaluate Unified Model simulations against two aircraft campaigns from summer 2017 over the Southeast Atlantic. The model generally performs well, such that discrepancies between the observational data sets may sometimes limit the precision of the evaluation. Our newly developed aging parameterization reproduces observed OA:BC mass ratios well and allows modeled OA:BC to decrease with smoke age, but increases bias in aerosol extinction and changes the BBA radiative effect little (+0.12 W m -2 ). We calculate aerosol refractive index using either a volume‐weighted component average or the Maxwell‐Garnett (MG) mixing assumption, which represents BC as small inclusions in a host material. Compared to MG mixing, the volume‐weighted average refractive index and reduced BC density increase aerosol absorption, substantially increasing the total BBA radiative effect (+2.66 W m -2 ) and amount of BBA transported across the ocean through BC self‐lofting. Reducing OA density to better match literature values changes the total BBA radiative effect by −1.96 W m -2 . Changes to direct radiative effects exceed changes to cloud radiative effects. Our findings emphasize the sensitivity of aerosol radiative effects and transport to these processes and properties, which we suggest could be improved in climate models.

Southeast Atlantic↗

Improving Projections of Antarctic Ice Sheet Contribution to Sea‐Level Change Through 2300 by Capturing Gravitational, Rotational, and Deformational Effects

Ice-mass change induces regionally varying patterns of sea-level change due to gravitational, rotational, and deformational (GRD) effects, which in turn influence marine-based ice stability in Antarctica. For improved projection of the Antarctic Ice Sheet (AIS), there is a need for including GRD effects in modeling and improving understanding of basin-by-basin sensitivity of ice evolution to GRD effects under a range of climate scenarios. We couple a high-resolution, higher-order ice-sheet model with a 1D global sea-level model that fully captures GRD effects, and simulate ice evolution in Antarctica under the Ice Sheet Model Intercomparison Project for CMIP6 experiments. We perform two sets of coupled simulations incorporating 1D Maxwell solid Earth structure suitable for West and East Antarctica and show that the Amundsen Sea Embayment (ASE) in West Antarctica has the highest sensitivity to GRD effects—in high-emission scenarios, grounding-line retreat accelerates by hundreds of kilometers by 2300 without GRD effects, but GRD effects delay this retreat on a timescale of decades. However, we find that delay times do not show a clear relationship to the strength of climate forcing alone. Furthermore, GRD effects can influence ice-sheet dynamics more than the choice of climate model for a given emissions scenario. In contrast, East Antarctica exhibits minimal sensitivity to GRD effects throughout the study period. These findings underscore the critical role of GRD effects in shaping future West AIS evolution, highlighting the importance of constraining the regional 3D Earth structure and bed topography in West Antarctica, particularly the ASE.

58 GEOSCIENCES↗

Nonlinear ionization dynamics of hot dense plasma observed in a laser-plasma amplifier

Understanding the behaviour of matter under conditions of extreme temperature, pressure, density and electromagnetic fields has profound effects on our understanding of cosmologic objects and the formation of the universe. Lacking direct access to such objects, our interpretation of observed data mainly relies on theoretical models. However, such models, which need to encompass nuclear physics, atomic physics and plasma physics over a huge dynamic range in the dimensions of energy and time, can only provide reliable information if we can benchmark them to experiments under well-defined laboratory conditions. Due to the plethora of effects occurring in this kind of highly excited matter, characterizing isolated dynamics or obtaining direct insight remains challenging. High-density plasmas are turbulent and opaque for radiation below the plasma frequency and allow only near-surface insight into ionization processes with visible wavelengths. Here, the output of a high-harmonic seeded laser-plasma amplifier using eight-fold ionized krypton as the gain medium operating at a 32.8 nm wavelength is ptychographically imaged. A complex-valued wavefront is observed in the extreme ultraviolet (XUV) beam with high resolution. Ab initio spatio-temporal Maxwell–Bloch simulations show excellent agreement with the experimental observations, revealing overionization of krypton in the plasma channel due to nonlinear laser-plasma interactions, successfully validating this four-dimensional multiscale model. This constitutes the first experimental observation of the laser ion abundance reshaping a laser-plasma amplifier. The presented approach shows the possibility of directly modelling light-plasma interactions in extreme conditions, such as those present during the early times of the universe, with direct experimental verification.

36 MATERIALS SCIENCE↗

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↗

Optical neural engine for solving scientific partial differential equations

Abstract Solving partial differential equations (PDEs) is the cornerstone of scientific research and development. Data-driven machine learning (ML) approaches are emerging to accelerate time-consuming and computation-intensive numerical simulations of PDEs. Although optical systems offer high-throughput and energy-efficient ML hardware, their demonstration for solving PDEs is limited. Here, we present an optical neural engine (ONE) architecture combining diffractive optical neural networks for Fourier space processing and optical crossbar structures for real space processing to solve time-dependent and time-independent PDEs in diverse disciplines, including Darcy flow equation, the magnetostatic Poisson’s equation in demagnetization, the Navier-Stokes equation in incompressible fluid, Maxwell’s equations in nanophotonic metasurfaces, and coupled PDEs in a multiphysics system. We numerically and experimentally demonstrate the capability of the ONE architecture, which not only leverages the advantages of high-performance dual-space processing for outperforming traditional PDE solvers and being comparable with state-of-the-art ML models but also can be implemented using optical computing hardware with unique features of low-energy and highly parallel constant-time processing irrespective of model scales and real-time reconfigurability for tackling multiple tasks with the same architecture. The demonstrated architecture offers a versatile and powerful platform for large-scale scientific and engineering computations.

Tang, Yingheng (ORCID:0009000153622546)↗