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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Two-Step Multi-Physics Analysis of an Annular Linear Induction Pump for Fission Power Systems

One of the key technologies associated with fission power systems (FPS) is the annular linear induction pump (ALIP). ALIPs are used to circulate liquid-metal fluid for transporting thermal energy from the nuclear reactor to the power conversion device. ALIPs designed and built to date for FPS project applications have not performed up to expectations. A unique, two-step approach was taken toward the multi-physics examination of an ALIP using ANSYS Maxwell 3D and Fluent. This multi-physics approach was developed so that engineers could investigate design variations that might improve pump performance. Of interest was to determine if simple geometric modifications could be made to the ALIP components with the goal of increasing the Lorentz forces acting on the liquid-metal fluid, which in turn would increase pumping capacity. The multi-physics model first calculates the Lorentz forces acting on the liquid metal fluid in the ALIP annulus. These forces are then used in a computational fluid dynamics simulation as (a) internal boundary conditions and (b) source functions in the momentum equations within the Navier-Stokes equations. The end result of the two-step analysis is a predicted pump pressure rise that can be compared with experimental data.

Fission Power Systems↗

NASA NTRS2016

Two-Step Multi-Physics Analysis of an Annular Linear Induction Pump for Fission Power Systems

Fission Power Systems↗

DOE OSTI2020

The Differential Equations That Model Diseases

In this brief tutorial, I will describe the simplest way of modeling a deadly infectious disease, such as COVID-19. I will show that at early times, the disease grows exponentially, but at late times, falls off like a Bell curve. I will solve the equations numerically with Python. A common mathematical model is the so-called compartmental model, where diseases can move people between "categories" such as healthy, infectious, or truly sick. This is a large, well established field and there are many good resources on this topic. I started with the SIAM review article The mathematics of infectious diseases. However, there are many excellent textbooks and more modern reviews as well.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Bringing Trimmed Serendipity Methods to Computational Practice in Firedrake

We present an implementation of the trimmed serendipity finite element family, using the open-source finite element package Firedrake. The new elements can be used seamlessly within the software suite for problems requiring H 1 , H (curl), or H (div)-conforming elements on meshes of squares or cubes. To test how well trimmed serendipity elements perform in comparison to traditional tensor product elements, we perform a sequence of numerical experiments including the primal Poisson, mixed Poisson, and Maxwell cavity eigenvalue problems. Overall, we find that the trimmed serendipity elements converge, as expected, at the same rate as the respective tensor product elements, while being able to offer significant savings in the time or memory required to solve certain problems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2018

Climates of Warm Earth-Like Planets. I. 3D Model Simulations

We present a large ensemble of simulations of an Earth-like world with increasing insolation and rotation rate. Unlike previous work utilizing idealized aquaplanet congurations we focus our simulations on modern Earth-like topography. The orbital period is the same as modern Earth, but with zero obliquity and eccentricity. The atmosphere is 1 bar N2-dominated with CO2=400 ppmv and CH4=1 ppmv. The simulations include two types of oceans; one without ocean heat transport (OHT) between grid cells as has been commonly used in the exoplanet literature, while the other is a fully coupled dynamic bathtub type ocean. The dynamical regime transitions that occur as day length increases induce climate feedbacks producing cooler temperatures, rst via the reduction of water vapor with increasing rotation period despite decreasing shortwave cooling by clouds, and then via decreasing water vapor and increasing shortwave cloud cooling, except at the highest insolations. Simulations without OHT are more sensitive to insolation changes for fast rotations while slower rotations are relatively insensitive to ocean choice. OHT runs with faster rotations tend to be similar with gyres transporting heat poleward making them warmer than those without OHT. For slower rotations OHT is directed equator-ward and no high latitude gyres are apparent. Uncertainties in cloud parameterization preclude a precise determination of habitability but do not a affect robust aspects of exoplanet climate sensitivity. This is the first paper in a series that will investigate aspects of habitability in the simulations presented herein. The datasets from this study are open source and publicly available.

Planetary systems↗

DOE OSTI2021

Electron-scale temperature gradients in kinetic equilibrium: MMS observations and Vlasov–Maxwell solutions

For an ideal gas, the notion of thermodynamic equilibrium requires that the system’s temperature be spatially uniform. We demonstrate that for a collisionless plasma, kinetic equilibrium can be achieved even in the presence of a temperature gradient. Motivated by recent MMS magnetopause observations of thin current layers exhibiting electron-scale temperature gradients, here we present an exact solution to the Vlasov-Maxwell system which accounts for the observed electron temperature gradient and reproduces much of the plasma’s spatial variation along one dimension. Furthermore, the current layer in the Vlasov-Maxwell model is self-consistently supported by a prominent crescent shaped electron velocity distribution that naturally arises due to the temperature transition, which agrees with the electron distribution structures measured by MMS during two asymmetric current layers encountered at the magnetopause.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

The gabbro-eclogite phase transition and the elevation of mountain belts on Venus

Among the four mountain belts surrounding Lakshmi Planum, Maxwell Montes is the highest and stands up to 11 km above the mean planetary radius and 7 km above Lakshmi Planum. The bulk composition and radioactive heat production of the crust on Venus, where measured, are similar to those of terrestrial tholeiitic basalt. Because the thickness of the low-density crust may be limited by the gabbro-garnet granulite-eclogite phase transitions, the 7-11 km maximum elevation of Maxwell Montes is difficult to understand except in the unlikely situation that the crust contains a large volume of magma. A possible explanation is that the base of the crust is not in phase equilibrium. It has been suggested that under completely dry conditions, the gabbro-eclogite phase transition takes place by solid-state diffusion and may require a geologically significant time to run to completion. Solid-state diffusion is a strongly temperature-dependent process. In this paper we solve the thermal evolution of the mountain belt to attempt to constrain the depth of the gabbro-eclogite transition and thus to assess this hypothesis quantitatively. The one-dimensional heat equation is solved numerically by a finite difference approximation. The deformation of the horizontally shortening crustal and mantle portions of the thermal boundary layer is assumed to occur by pure shear, and therefore the vertical velocity is given by the product of the horizontal strain rate and depth.

Namiki, Noriyuki↗

DOE OSTI2025

Extracting scattering amplitudes for arbitrary two-particle systems with one-particle left-hand cuts via lattice QCD

We derive a general formalism that relates the spectrum of two-particle systems in a finite volume to physical scattering amplitudes, taking into account the presence of any left-hand branch cuts due to single-particle exchanges. The method first relates the finite-volume spectrum to an infinite-volume short-range quantity, denoted ${\mathcal{M}}_0$, and then relates the latter to the physical scattering amplitudes via known integral equations. The derivation of both relations is performed using all-orders perturbation theory and is exact up to neglected exponentially suppressed volume dependence. The relations hold for arbitrary two-particle systems with any number of coupled channels, non-identical and non-degenerate particles, and any intrinsic spin.

algorithms↗

DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗

DOE OSTI2022

Machine learning of consistent thermodynamic models using automatic differentiation

In this study, we propose a data-driven method to describe consistent equations of state (EOS) for arbitrary systems. Complex EOS are traditionally obtained by fitting suitable analytical expressions to thermophysical data. A key aspect of EOS is that the relationships between state variables are given by derivatives of the system free energy. In this work, we model the free energy with an artificial neural network and utilize automatic differentiation to directly learn the derivatives of the free energy. We demonstrate this approach on two different systems, the analytic van der Waals EOS and published data for the Lennard-Jones fluid, and we show that it is advantageous over direct learning of thermodynamic properties (i.e., not as derivatives of the free energy but as independent properties), in terms of both accuracy and the exact preservation of the Maxwell relations. Furthermore, the method implicitly provides the free energy of a system without explicit integration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Decay amplitudes to three hadrons from finite-volume matrix elements

We derive relations between finite-volume matrix elements and infinite-volume decay amplitudes, for processes with three spinless, degenerate and either identical or non-identical particles in the final state. This generalizes the Lellouch-Lüscher relation for two-particle decays and provides a strategy for extracting three-hadron decay amplitudes using lattice QCD. Unlike for two particles, even in the simplest approximation, one must solve integral equations to obtain the physical decay amplitude, a consequence of the nontrivial finite-state interactions. We first derive the result in a simplified theory with three identical particles, and then present the generalizations needed to study phenomenologically relevant three-pion decays. The specific processes we discuss are the CP-violating K → 3π weak decay, the isospin-breaking η → 3π QCD transition, and the electromagnetic γ* → 3π amplitudes that enter the calculation of the hadronic vacuum polarization contribution to muonic g - 2.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1978

Monte Carlo solution of Boltzmann equation for a simple model of highly nonequilibrium diatomic gases - Translational rotational energy relaxation

Theoretical results pertaining to internally excited translational-rotational energy relaxation in a spatially uniform diatomic gas far removed from solid boundaries are obtained by solving the Boltzmann equation by means of the Monte Carlo direct simulation method. The analysis is based on calculations involving three different types of initial conditions: equilibrium, nonequilibrium-equipartition (i.e., equipartition is satisfied, but the distributions are perturbed), and nonequilibrium-nonequipartition (i.e., both equipartition and the distributions are perturbed). Results of monatomic-gas simulations are also included to facilitate comparisons with the coupled translational-rotational relaxation simulations, and some simulations for a normal shock-wave structure are briefly examined. The results show that: (1) single-step transitions are the significant mechanisms of intermodal energy transfer; (2) translational-rotational transitions are coupled most efficiently for low-lying states of rotationally excited molecules and least efficiently for highly rotationally excited molecules; and (3) relaxation occurs via a successive set of distributions that are not Maxwell-Boltzmann (nonlocal Maxwellian).

Yoshikawa, K. K.↗

DOE OSTI2025

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗

DOE OSTI2022

Supergravity black holes, Love numbers, and harmonic coordinates

A problem which has now come to a head in general relativity, primarily as a result of the many recent detections of gravitational waves from the coalescence of compact binaries, is that there are very few viable theories against which to “test” it. To perform realistic tests of theories of gravity, we need to be able to look beyond general relativity and evaluate the consistency of a parametrized, physically acceptable, family of black hole metric alternatives with observational data from, especially, gravitational wave detections using, for example, an agnostic Bayesian approach. In this paper we further examine properties of one class of such metrics, which in fact arise as solutions of ungauged supergravity. In particular, we examine the massless, neutral, minimally coupled scalar wave equation in a general stationary, axisymmetric background metric such as that of a charged rotating black hole, when the scalar field is either time independent or in the low-frequency, near-zone limit, with a view to calculating the Love numbers of tidal perturbations, and of obtaining harmonic coordinates for the background metric. For a four-parameter family of charged asymptotically flat rotating black hole solutions of ungauged supergravity theory known as STU black holes, which includes Kaluza-Klein black holes and the Kerr-Sen black hole as special cases, we find that all time-independent solutions, and hence the harmonic coordinates of the metrics, are identical to those of the Kerr solution. In the low-frequency limit we find the scalar fields exhibit the same SL(2,R) symmetry as holds in the case of the Kerr solution. We point out extensions of our results to a wider class of metrics, which includes solutions of Einstein-Maxwell-dilaton theory.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Einstein-Maxwell theory and the Weyl double copy

The Weyl double copy relates vacuum solutions in general relativity to Abelian gauge fields in Minkowski spacetime. In a previous work, we showed how the Weyl double copy can be extended to provide a treatment of external gravitational sources consistent with the classical Kerr-Schild double copy. Using this generalization, here we provide a complete double copy analysis of electrovacuum Petrov type D spacetimes. This includes the first analysis of the charged C-metric, whose single copy interpretation invokes the two-potential formalism of electrodynamics. Here, we also present the first double copy prescription for the Ricci spinor, which for non-accelerating spacetimes, takes a form similar to the original double copy relation for the Weyl spinor.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2013

Modeling Insights into Deuterium Excess as an Indicator of Water Vapor Source Conditions

Deuterium excess (d) is interpreted in conventional paleoclimate reconstructions as a tracer of oceanic source region conditions, such as temperature, where precipitation originates. Previous studies have adopted co-isotopic approaches to estimate past changes in both site and oceanic source temperatures for ice core sites using empirical relationships derived from conceptual distillation models, particularly Mixed Cloud Isotopic Models (MCIMs). However, the relationship between d and oceanic surface conditions remains unclear in past contexts. We investigate this climate-isotope relationship for sites in Greenland and Antarctica using multiple simulations of the water isotope-enabled Goddard Institute for Space Studies (GISS) ModelE-R general circulation model and apply a novel suite of model vapor source distribution (VSD) tracers to assess d as a proxy for source temperature variability under a range of climatic conditions. Simulated average source temperatures determined by the VSDs are compared to synthetic source temperature estimates calculated using MCIM equations linking d to source region conditions. We show that although deuterium excess is generally a faithful tracer of source temperatures as estimated by the MCIM approach, large discrepancies in the isotope-climate relationship occur around Greenland during the Last Glacial Maximum simulation, when precipitation seasonality and moisture source regions were notably different from present. This identified sensitivity in d as a source temperature proxy suggests that quantitative climate reconstructions from deuterium excess should be treated with caution for some sites when boundary conditions are significantly different from the present day. Also, the exclusion of the influence of humidity and other evaporative source changes in MCIM regressions may be a limitation of quantifying source temperature fluctuations from deuterium excess in some instances.

simulation↗

DOE OSTI2020

Generalizing the relativistic quantization condition to include all three-pion isospin channels

We present a generalization of the relativistic, finite-volume, three-particle quantization condition for non-identical pions in isosymmetric QCD. The resulting formalism allows one to use discrete finite-volume energies, determined using lattice QCD, to constrain scattering amplitudes for all possible values of two- and three-pion isospin. As for the case of identical pions considered previously, the result splits into two steps: the first defines a non-perturbative function with roots equal to the allowed energies, E n (L), in a given cubic volume with side-length L. This function depends on an intermediate three-body quantity, denoted K d f , 3 , which can thus be constrained from lattice QCD in- put. The second step is a set of integral equations relating K d f , 3 to the physical scattering amplitude, M 3 . Both of the key relations, En(L) ↔ K d f , 3 and K d f , 3 ↔ M 3 , are shown to be block-diagonal in the basis of definite three-pion isospin, Iπππ , so that one in fact recovers four independent relations, corresponding to Iπππ = 0, 1, 2, 3. We also provide the generalized threshold expansion of K d f , 3 for all channels, as well as parameterizations for all three-pion resonances present for Iπππ = 0 and Iπππ = 1. As an example of the utility of the generalized formalism, we present a toy implementation of the quantization condition for Iπππ = 0, focusing on the quantum numbers of the ω and h 1 resonances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗