Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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

Reduced-model scrape-off layer turbulence (nSOLT) simulations comparing three fueling scenarios

The 2D scrape-off-layer turbulence code (nSOLT) includes 1D Boltzmann neutral–plasma interactions, a model of divertor recycling (introduced here), and a fixed source of plasma concentrated at the core-side boundary. Three fueling methods are considered herein: (1) neutral injection in the far-SOL is accomplished by specifying the density of Franck–Condon distributed neutrals streaming in from the boundary. (2) Divertor recycling is modeled by injecting a fraction of the particle parallel flux in the scrape-off layer (SOL) back into the edge region as a source of plasma. (3) A constant source fuels the edge plasma from the core-side boundary to model pellet injection. For machine parameters (B, R m , and L // ) illustrative of the MAST-U device, and for a deuterium plasma, turbulent equilibria are obtained that share the same plasma fueling rate for each of the three fueling methods, with only one of the sources on in each case. In the presence of self-consistent turbulence, quasi-steady plasma and neutral (deuterium) profiles, fueling efficiencies, SOL transparencies, and heat flux widths are compared. Characteristics of the turbulent fluctuations, including skewness, cross-phases, and power spectra, are described. The calculated fueling efficiencies, SOL transparencies to neutral penetration, and many of the turbulent properties are remarkably similar for all three fueling methods despite significant differences in the plasma profiles. The nonlinear states of the three cases are dominated by separatrix-spanning vortex cells that control particle and heat losses into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗

Magnetic reconnection during magnetospheric substorms

The near earth reconnection model of substorms represents an attempt to place a broad range of observations into a consistent framework. The roles and requirements of reconnection are discussed. High speed plasma sheet flows, thin current sheet instability, substorm triggering, plasmoids and flux ropes in the distant tail, and magnetohydrodynamic simulations are discussed. Substorms are global, coherent sequences of processes involving solar wind/magnetosphere/ionosphere interaction. Magnetic reconnection is required to explain different dayside and polar cap phenomena, which required nightside reconnection. The modification and expansion of the standard near earth neutral line (NENL) model can integrate breakup arcs, current disruption, current wedge features, and localized plasma flows into the magnetic reconnection framework.

Baker, Daniel N.↗

The ionospheric peak on the Venus dayside

The behavior of the ionospheric peak on the dayside of Venus is described and interpreted by combining radio occultation measurements with theoretical calculations. The theoretical models are shown to be able to reproduce the measured electron densities very accurately when careful consideration is given to such parameters as the level of solar activity, the electron temperature, and the neutral density. What is more, the models are able to provide a check on the accuracy of neutral atmospheric models in the vicinity of 140 km. Chemical equilibrium is assumed for the calculation of ion and electron densities. A table giving Pioneer Venus radio occultation measurements of the Venus dayside ionosphere is included.

Cravens, T. E.↗

Lyman alpha emissions in the outer heliosphere

Since 1993, the Ultraviolet Spectrometers aboard both Voyager 1 and 2, were used to study the Ly-alpha glow pattern backscattered by neutral hydrogen atoms at great distance from the sun. When compared to a radiative transfer calculation for an unperturbed interstellar hydrogen flow into the solar system, the data present an excess of intensity seen in the direction of the incoming interstellar wind. An absolute estimate of this Lyman alpha intensity yields 10 to 15 Rayleighs for each spacecraft. We discuss the possibility that this excess is caused by a gradient of hydrogen at large distance from the Sun due to the filtration of interstellar hydrogen when crossing the heliopause. Comparisons with results from the Baranov-Malama model of neutral-plasma coupling at the heliopause are presented. The data yield a range of possible values for the hydrogen gradient which may be used to infer the location of the heliopause.

Quemerais, E.↗

Time-dependent MHD modeling of the global structure of the heliosphere

We present results from time-dependent modeling of the global structure of the heliosphere with neutral and magnetic field effects included. The magnetic field is assumed parallel to the interstellar flow in this two-dimensional axisymmetric model; the neutrals are treated as a fluid. The effects of interstellar neutrals and the interplanetary magnetic field on the location of the termination shock are studied using the most recent estimate of the interstellar medium parameters, results will be compared to those of Baranov and Zaitsev. The effect of the solar wind - VLISM interaction on the density and velocity of interstellar neutrals within the heliosphere will also be presented and related to observations. The response of the termination shock to the solar cycle variation in the solar wind will be compared to the response found previously using an axisymmetric hydrodynamic model without neutrals.

Liewer, P. C.↗

Theoretical spectral distributions and total cross sections for neutral subthreshold pion production in carbon-carbon collisions

A coherent isobar formalism is employed to model subthreshold production of neutral pions in carbon-carbon collisions at energies below 100 MeV/nucleon. No arbitrary scale factors or adjustable free parameters are used in calculation of the Lorentz-invariant cross sections for pion production in the projectile, which produces an excited state that goes to M1 resonance in the target by conservation of spin and isospin. Pion production is also modeled for the projectile, which also reaches M1 resonance. The overall pion spectral distribution in the center of mass system is then integrated over the energy range 35-84 MeV/nucleon. The results expose an energy loss in the incident ions, as observed experimentally, and indicate that an isobar mechanism is responsible for higher energy pion production. Lower energy pions are a result of thermal processes.

Norbury, J. W.↗

Satellite tori at Saturn

The inner satellites of Saturn are icy bodies imbedded in a plasma environment in which they are continuously bombarded by energetic ions, corotating plasma, and solar radiation. Laboratory sputtering experiments indicate that this should result in the injection of substantial amounts of neutral H, H2, OH, H2O, and O2 into the magnetosphere. The atomic processes affecting these neutrals and the neutrals and ions formed from them are modeled, and the steady state neutral and ion densities expected in the plasma tori of Enceladus, Dione-Tethys, and Rhea are calculated. Comparison with observations shows that recombination can limit the Enceladus and Dione-Tethys tori to the observed densities, but that transport rates of at least 4 x 10 to the -8th Saturn radii squared/s are required to limit torus densities at Rhea to the observed values.

Richardson, J. D.↗

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A design methodology for neutral buoyancy simulation of space operations

Neutral buoyancy has often been used in the past for EVA development activities, but little has been done to provide an analytical understanding of the environment and its correlation with space. This paper covers a set of related research topics at the MIT Space Systems Laboratory, dealing with the modeling of the space and underwater environments, validation of the models through testing in neutral buoyancy, parabolic flight, and space flight experiments, and applications of the models to gain a better design methodology for creating meaningful neutral buoyancy simulations. Examples covered include simulation validation criteria for human body dynamics, and for applied torques in a beam rotation task, which is the pacing crew operation for EVA structural assembly. Extensions of the dynamics models are presented for powered vehicles in the underwater environment, and examples given from the MIT Space Telerobotics Research Program, including the Beam Assembly Teleoperator and the Multimode Proximity Operations Device. Future expansions of the modeling theory are also presented, leading to remote vehicles which behave in neutral buoyancy exactly as the modeled system would in space.

Akin, David L.↗

Effects of a wavy neutral sheet on cosmic ray anisotropies

The first results of a three-dimensional numerical code calculating cosmic ray anisotropies is presented. The code includes diffusion, convection, adiabatic cooling, and drift in an interplanetary magnetic field model containing a wavy neutral sheet. The 3-D model can reproduce all the principal observations for a reasonable set of parameters.

Kota, J.↗

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗

Nonlinear Alfvén instability simulation and EP transport for ITER reversed shear (steady-state) and monotonic q-profile regimes *

The nonlinear evolution and dynamics of energetic particle (EP) driven instabilities in both steady-state (reversed shear q-profile) and monotonic q-profile regimes of ITER are examined using the global gyro-Landau closure model FAR3d. Both neutral beam and alpha components are included, leading to synergistic effects between the two populations. In order to maintain computational feasibility, the present nonlinear simulation includes all toroidal mode numbers ranging from n = 0 to 15. It cannot be excluded that n’s above this range could introduce larger linear growth rates; however, as indicated from the results, the n = 0 to 15 range is sufficient to provide a strong linear drive and significant nonlinear transport effects. While the toroidal mode numbers n = 6,13 dominate the linear instability growth phase, nonlinear energy transfers in the saturated phase reverse this trend, leading to the dominance of lower n’s in the saturated phase of the simulation. Zonal flow structure generation and EP density flattening are further consequences of the nonlinear phase. The density profile flattening is caused by collectively driven EP radial transport fluxes, which have non-local characteristics. Instantaneous alpha particle transport fluxes are significant (Γ α ~ 2 x 10 20 m –2 s –1 ) for the reversed shear regime and lower (Γ α ~ 1.1 x 10 20 m –2 s –1 ) for the monotonic q-profile case.

AE↗

Strong-field ionization of the triplet ground state of O 2

Using strong-field ionization as a probe, we observe highly nonperiodic evolution of the spin-rotation wave packet launched by a nonionizing femtosecond pulse in oxygen. The nonperiodicity is readily apparent only in rotationally cold molecules that are pumped with a weak alignment pulse. We show that this behavior is a consequence of the spin-rotation and the spin-spin couplings in the triplet ground state of the neutral molecule. A model that includes these couplings in the field-free Hamiltonian but in neither the alignment nor the ionization step explains most of the observed dynamics, suggesting that neither process depends explicitly on the electronic spin. We also show that the angle dependence of strong-field ionization can be retrieved from the delay-dependent signal even when coupling to spin complicates the rotational dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Freedom of the mixer rotation axis improves performance in the quantum approximate optimization algorithm

Variational quantum algorithms such as the quantum approximate optimization algorithm (QAOA) are particularly attractive candidates for implementation on near-term quantum processors. As hardware realities such as error and qubit connectivity will constrain achievable circuit depth in the near future, new ways to achieve high performance at low depth are of great interest. In this work, we present a modification to QAOA that adds additional variational parameters in the form of freedom of the rotation axis in the XY plane of the mixer Hamiltonian. Here, via numerical simulation, we show that this leads to a drastic performance improvement over standard QAOA at finding solutions to the MaxCut problem on graphs of up to seven qubits. Furthermore, we explore the Z-phase error mitigation properties of our modified Ansatz, its performance under a realistic error model for a neutral atom quantum processor, and the class of problems it can solve in a single round.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Axion baryogenesis puts a new spin on the Hubble tension

We show that a rotating axion field that makes a transition from a matterlike equation of state to a kinationlike equation of state around the epoch of recombination can significantly ameliorate the Hubble tension, i.e., the discrepancy between the determinations of the present-day expansion rate H 0 from observations of the cosmic microwave background on one hand and type Ia supernovae on the other. We consider a specific, UV-complete model of such a rotating axion and find that it can relax the Hubble tension without exacerbating tensions in determinations of other cosmological parameters, in particular the amplitude of matter fluctuations S 8 . We subsequently demonstrate how this rotating axion model can also generate the baryon asymmetry of our Universe, by introducing a coupling of the axion field to right-handed neutrinos. This baryogenesis model predicts heavy neutral leptons that are most naturally within reach of future lepton colliders, but in finely tuned regions of parameter space may also be accessible at the high-luminosity LHC and the beam dump experiment SHiP. Published by the American Physical Society 2024

Co, Raymond T. (ORCID:0000000283957056)↗