Search NASA⌕ Search

SEARCH · Search NASA

Results for “Slab geometry”

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 37 records · Page 2

Verification of local electrostatic gyrokinetic simulation of driftwave instability in field-reversed configuration

Local linear simulations of ion temperature gradient (ITG) instabilities using a gyrokinetic particle code GKD1D have been verified by comparing with the analytic dispersion relation in the slab geometry and by benchmarking with a global gyrokinetic code GTC-X in a realistic field-reversed configuration (FRC). Results of ITG simulations using either adiabatic or kinetic electrons exhibit excellent agreements between the two codes for linear frequencies, growth rates, and mode structures. The ITG modes in the FRC scrape-off layer region are mostly dominated by the poloidal m = 1 harmonic. Kinetic electrons enhance the subdominant m = 2 harmonic and the even parity dominates the ITG mode structure. Kinetic electrons increase the growth rates for the low-n modes, but decrease the growth rates for the high-n modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dispersion relation for gauge-free electromagnetic drift kinetics

Recently, a new approach to gyrokinetics, invariant under electromagnetic gauge transformations, was developed. The gyrocenter equations of motion are now expressed in terms of the perturbed fields instead of the potentials, in a form suitable for numerical simulations and analytic studies. In this paper, we verify that the long-wavelength limit, i.e., the drift-kinetic limit of the new gyrokinetic theory, is in line with existing work, providing a solid foundation for simulations. Here, we compute the dispersion relation of the new drift-kinetic theory in slab geometry and find agreement with a long-wavelength limit of the full Vlasov–Maxwell model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electric fields near undulating dielectric membranes

Dielectric interfaces are crucial to the behavior of charged membranes, from graphene to synthetic and biological lipid bilayers. Understanding electrolyte behavior near these interfaces remains a challenge, especially in the case of rough dielectric surfaces. A lack of analytical solutions consigns this problem to numerical treatments. Here, we report an analytic method for determining electrostatic potentials near curved dielectric membranes in a two-dimensional periodic “slab” geometry using a periodic summation of Green’s functions. This method is amenable to simulating arbitrary groups of charges near surfaces with two-dimensional deformations. We concentrate on one-dimensional undulations. We show that increasing membrane undulation increases the asymmetry of interfacial charge distributions due to preferential ionic repulsion from troughs. In the limit of thick membranes, we recover results mimicking those for electrolytes near a single interface. Our work demonstrates that rough surfaces generate charge patterns in electrolytes of charged molecules or mixed-valence ions.

Chemistry↗

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probability of Initiation in Neutron Transport

We discuss the numerical solution of the nonlinear integro-differential equation for the probability of a divergent neutron chain in a stationary system (i.e., the probability of initiation (POI)). We follow the development described in Bell’s classic paper on the stochastic theory of neutron transport. As noted by Bell, the linearized form of this equation resembles the linear adjoint neutron transport equation. A matrix formalism for the discretized steady state (or forward) neutron equation in slab geometry is first developed and is then used to derive the discrete adjoint equation. A main advantage of this discrete development is that the resulting discrete adjoint equation does not depend upon how the multigroup cross sections for the forward problem are obtained. That is, we derive the discrete adjoint directly from the discrete forward equations rather than discretizing directly the adjoint equation. This also guarantees that the discrete adjoint operator is consistent with the inner product used to define the adjoint operator. We discuss three approaches for the numerical solution of the POI equations, and present numerical results on several test problems. The three solution methods are a simple fixed-point iteration, a second approach that is akin to a nonlinear Power iteration, and a third approach which uses a Newton-Krylov nonlinear solver. We also give sufficient conditions to guarantee the existence and uniqueness of nontrivial solutions to our discrete POI equations when the discrete system is supercritical, and that only the trivial solution exists when the discrete system is subcritical. Our approach is modeled after the analysis presented for the continuous POI equations by Mokhtar-Kharroubi and Jarmouni-Idrissi, and by Pazy and Rabinowitz.

42 ENGINEERING↗

Direct prediction of saturated neoclassical tearing modes in slab using an equilibrium approach

We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al 2020 Phys. Plasmas 27 070701; Loizu and Bonfiglio 2023 J. Plasma Phys. 89 905890507), we make use of Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502), an equilibrium solver based on the variational principle of the multi-region relaxed magnetohydrodynamics (MHDs), featuring stepped pressure profiles and arbitrary magnetic topology. We work in slab geometry and employ a simple bootstrap current model J bs = C$\boldsymbol{\nabla}$p to study the bootstrap-driven tearing modes, scanning over the asymptotic matching parameter Δ' and bootstrap current strength. Saturated island widths produced by SPEC agree well with the predictions of an initial value resistive MHDs code (Huang and Bhattacharjee 2016 Astrophys. J. 818 20) while being orders of magnitude faster to calculate. Additionally, we observe good agreement with a simple analytical modified Rutherford equation, without requiring any fitting coefficients. The match is obtained for both linearly unstable classical tearing modes in the presence of bootstrap current, and NTMs, which are linearly stable but nonlinear-unstable due to the effects of the bootstrap current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Recent progress in L-H transition studies at JET: Tritium, Helium, Hydrogen and Deuterium

In this report we present an overview of results from a series of L-H transition experiments undertaken at JET since the installation of the ITER-like-wall (JET-ILW), with Beryllium wall tiles and a Tungsten divertor. Tritium, Helium and Deuterium plasmas have been investigated. Initial results in Tritium show ohmic L-H transitions at low density and the power threshold for the L-H transition (P LH ) is lower in Tritium plasmas than in Deuterium ones at low densities, while we still lack contrasted data to provide a scaling at high densities. In Helium plasmas there is a notable shift of the density at which the power threshold is minimum (n e,min ) to higher values relative to Deuterium and Hydrogen references. Above n e,min (He) the L-H power threshold at high densities is similar for D and He plasmas. Transport modelling in slab geometry shows that in Helium neoclassical transport competes with interchange-driven transport, unlike in Hydrogen isotopes. Measurements of the radial electric field in Deuterium plasmas show that E r shear is not a good indicator of proximity to the L-H transition. Transport analysis of ion heat flux in Deuterium plasmas show a non-linearity as density is decreased below n e,min . Lastly, a regression of the JET-ILW Deuterium data is compared to the 2008 ITPA scaling law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation study of particle transport by weakly coherent mode in the Alcator C-Mod tokamak

A simulation study has been conducted of the physical mechanisms behind the weakly coherent mode (WCM) and its produced particle transport in the I-mode edge plasmas by using the BOUT++ code. The WCM is identified in our simulations by its poloidal and radial distributions as well as its frequency and wavenumber spectra. Its produced radial particle flux is calculated and compared with the experimental value. The good agreement indicates that the WCM is an important particle transport channel in the I-mode pedestal. It is found that the WCM can transport particles across the strong outer shear layer of the Er well established in the formation of I-mode, based on which a possible explanation is provided why I-mode does not feature a density pedestal. The key point lies in the change of the cross-phase between the electric potential and density fluctuations induced by the E × B Doppler shift. In the strong shear layer, although the electric potential fluctuation is significantly suppressed, the cross-phase is close to π/2, resulting in a strong drive of the density fluctuation and particle transport. To identify the physical nature of the WCM, a linear dispersion relation for drift Alfvén modes is derived in the slab geometry. A drift Alfvén wave instability is found to have similar dependence to the simulated linear instability behind the WCM on the resistivity and the parallel electron pressure gradient and thermal force terms in the parallel Ohm's law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics of beam-driven ion cyclotron emission in the large plasma device

Abstract Ion cyclotron emission (ICE) is widely observed from spatially localised minority energetic ion populations in toroidal magnetically confined fusion (MCF) plasmas, both tokamaks and stellarators. Its spectral structure is typically regular with narrow suprathermal peaks, whose frequency separation matches a local energetic ion cyclotron frequency. Here we report the first computational (fully nonlinear self-consistent kinetic particle-in-cell code) and analytical (linear magnetoacoustic cyclotron instability (MCI)) studies of ICE observations from cylindrical plasmas contained in the Large Plasma Device (LAPD). Because LAPD is cylindrical, the plasma physics giving rise to the observed ICE spectrum necessarily excludes toroidal effects. Our approach, previously successful for toroidal plasma ICE, assumes slab geometry and hence is well adapted to LAPD. ICE from LAPD is strongly electrostatic, as distinct from electromagnetic, and is driven by 15 keV beam ions for which the ratio of perpendicular speed to the local Alfven speed, v ⊥ / v A , is 0.15, lower than in MCF plasmas from which beam-driven ICE has previously been observed. Our results are in good agreement with these observations. There is congruence between simulated ICE spectra, obtained in the saturated nonlinear regime of our computations, and observed ICE spectra; the underlying physics is essentially the same as in toroidal plasmas; and there is alignment with linear analytical theory where appropriate. The present work establishes a mapping from the cylindrical LAPD ICE observations to toroidal MCF ICE observations. The LAPD spectra are instances of beam-driven MCI-type ICE in its sub-Alfvenic, predominantly electrostatic manifestation, which has precedents in MCF stretching back to the 1990s. An interesting corollary is that, for many purposes, ICE in toroidal MCF plasmas ‘might as well’ be occurring in a cylinder.

Samant, O. (ORCID:0000000226055363)↗

Uniform convergence of an upwind discontinuous Galerkin method for solving scaled discrete-ordinate radiative transfer equations with isotropic scattering

Here we present an error analysis for the discontinuous Galerkin (DG) method applied to the discrete-ordinate discretization of the steady-state radiative transfer equation with isotropic scattering. Under some mild assumptions, we show that the DG method converges uniformly with respect to a scaling parameter which characterizes the strength of scattering in the system. However, the rate is not optimal and can be polluted by the presence of boundary layers. In one-dimensional slab geometries, we demonstrate optimal convergence when boundary layers are not present and analyze a simple strategy for balance interior and boundary layer errors. Some numerical tests are also provided in this reduced setting.

97 MATHEMATICS AND COMPUTING↗

Extended interaction length laser-driven acceleration in a tunable dielectric structure

The development of long, tunable structures is critical to increasing energy gain in laser-driven dielectric accelerators (DLAs). Here we combine pulse-front-tilt illumination with slab-geometry structures assembled by precisely aligning off-the-shelf 4 mm long transmission gratings to achieve up to 200 keV energy modulation for 6 MeV injected electrons. The effective interaction length is longer than 1 mm, limited by the dephasing of the accelerated particles in the structure. The piezo-based independent mounting system for the gratings allows tuning of the gap and field distribution inside the structure. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Dielectric screening and electric field control of ferromagnetism at the CaMnO 3 /CaRuO 3 interface

Control of magnetism by an applied electric field is a desirable technique for the functionalization of magnetic materials. Motivated by recent experiments, we study the electric field control of the interfacial magnetism of CaRuO 3 /CaMnO 3 (CRO/CMO) (001), a prototype interface between a non-magnetic metal and an antiferromagnetic insulator. Even without the electric field, the interfacial CMO layer acquires a ferromagnetic moment due to a spin-canted state, caused by the Anderson-Hasegawa double exchange (DEX) between the Mn moments and the leaked electrons from the CRO side. An electric field would alter the carrier density at the interface, leading to the possibility of controlling the magnetism, since DEX is sensitive to the carrier density. We study this effect quantitatively using density-functional calculations in the slab geometry. We find a text-book like dielectric screening of the electric field, which introduces polarization charges at the interfaces and the surfaces. The extra charge at the interface enhances the ferromagnetism via the DEX interaction, while away from the interface the original AFM state of the Mn layers remains unchanged. The effect could have potential application in spintronics devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

AK112: Full Waveform Inversion Tomography of Alaska Improves Waveform Fits While Imaging Crustal, Mantle, and Slab Structure

We report a full waveform inversion tomography model of Alaska and the surrounding regions, inferring radially anisotropic shear and isotropic compressional wavespeeds by fitting complete waveforms from 120 regional earthquakes. Our multiscale approach inverted time–frequency phase misfits (maximum period of 100 s), starting with a minimum period of 40 s and ending at 20 s in 7 stages and 112 total iterations. The model (AK112) was evaluated by computing the misfits for 36 independent validation events. We find that misfit reductions were large and equal (∼55%) for both the inversion and validation data sets, providing confidence in the model. AK112 also provides much better waveform fits compared to other reported models for the region, including an isotropic version of itself, highlighting the importance of anisotropy. The model resolves known crustal, upper mantle, and slab structure to depths of 100 km with new detail: sedimentary basins in the Alaskan Shelf, Cook Inlet, and Colville basins, among others; discontinuous lithospheric structure across major terrane boundaries; and subducting slab geometry and back‐arc volcanic sources. In addition to tectonic interpretations, the model enables full waveform simulations for long‐period earthquake ground motions and source characterization (e.g., moment tensor and finite‐fault inversion).

Rodgers, Arthur [Lawrence Livermore National Labor↗

Eigenmode Analysis of Pulsed Neutron Transport Simulations

We discuss the time dependent behavior of some simple pulsed neutron simulations of subcritical problems in slab geometry. Our intent is to investigate the eigenvalue structure of the discretized neutron transport equation and to show under some reasonable assumptions that a dominant time eigenvalue exists that has a nonnegative eigenvector that determines the long time dependent behavior of the solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗

Cylinder block casting slab core cast geometry for sawcut entrance enhancement

An automobile vehicle engine includes multiple water jackets individually formed in a cast engine block proximate to successive ones of multiple cylinder bores. Multiple cast-in place transition regions are individually formed during a casting operation of the cast engine block at entrances to individual ones of the multiple water jackets. Individual ones of multiple sawcuts open into individual ones of the multiple cast-in place transition regions.

Petrus, Ronald John↗