Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transport theory”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK Utility Programs for DIF3D

The DIF3D-VARIANT solver is part of the DIF3D code which was developed at ANL for engineering level predictive analysis of nuclear reactor cores. Over the course of the last 20 years, the DIF3D-VARIANT option of DIF3D has become increasingly more used because of improvements in computer memory capabilities and performance. This is beneficial to the users of DIF3D as DIF3D-VARIANT is able to apply both p- and h- spatial mesh refinement strategies in diffusion and transport theory. Given its greater use, utility programs were created to assist in post-processing the DIF3D output. This document discusses the basic execution process of DIF3D and the output files it creates but is primarily focused on being the manual for the utility programs EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Rational sub-nanometer manipulation of polymer morphology for efficient chemical separations

The ability to control polymer morphology on the sub-nanometer length scale has broad implications for chemical separations. To achieve such control on easily processable systems, this proposal focuses on the synthesis and characterization of polymers containing appended labile moieties that are easily detached by thermolysis or UV irradiation deep within the glassy state. Once liberated, these moieties can diffuse from the polymer matrix as gaseous products, leaving behind templated pathways for selective diffusion and sorption of small molecules. With a specific target of creating polymeric membrane materials with unprecedented diffusion and (ad)sorption characteristics for chemical separations, synthesis of new materials will be complemented with advanced metrologies, simulations, and evaluation of thermodynamic and transport theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK Utility Programs for DIF3D (Rev. 1)

The DIF3D-VARIANT solver is part of the DIF3D code which was developed at ANL for engineering level predictive analysis of nuclear reactor cores. Over the course of the last 20 years, the DIF3D-VARIANT option of DIF3D has become increasingly more used because of improvements in computer memory capabilities and performance. This is beneficial to the users of DIF3D as DIF3D-VARIANT is able to apply both p- and h-spatial mesh refinement strategies in diffusion and transport theory. Given its greater use, utility programs were created to assist in post-processing the DIF3D output. This document discusses the basic execution process of DIF3D and the output files it creates but is primarily focused on being the manual for the utility programs EvaluateFlux, VARPOW, VARPEAK, and DIF3DtoVTK.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Drift kinetic theory of alpha transport by tokamak perturbations

We report that upcoming tokamak experiments fuelled with deuterium and tritium are expected to have large alpha particle populations. Such experiments motivate new attention to the theory of alpha particle confinement and transport. A key topic is the interaction of alpha particles with perturbations to the tokamak fields, including those from ripple and magnetohydrodynamic modes like Alfvén eigenmodes. These perturbations can transport alphas, leading to changed localization of alpha heating, loss of alpha power and damage to device walls. Alpha interaction with these perturbations is often studied with single-particle theory. In contrast, we derive a drift kinetic theory to calculate the alpha heat flux resulting from arbitrary perturbation frequency and periodicity (provided these can be studied drift kinetically). Novel features of the theory include the retention of a large effective collision frequency resulting from the resonant alpha collisional boundary layer, correlated interactions over many poloidal transits and finite orbit effects. Heat fluxes are considered for the example cases of ripple and the toroidal Alfvén eigenmode (TAE). The ripple heat flux is small. The TAE heat flux is significant and scales with the square of the perturbation amplitude, allowing the derivation of constraints on mode amplitude for avoidance of significant alpha depletion. A simple saturation condition suggests that TAEs in one upcoming experiment will not cause significant alpha transport via the mechanisms in this theory. However, saturation above the level suggested by the simple condition, but within numerical and experimental experience, which could be accompanied by the onset of stochasticity, could cause significant transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Chiral kinetic theory of anomalous transport induced by torsion

In Weyl semimetals subjected to torsion, there are two different kinds of chirality: (i) The (coordinate-space) shape of the twisted crystal is chiral and (ii) the momentum space contains chiral quasiparticles. Here we construct a general kinetic theory of anomalous transport using the phase space (coordinate and momentum spaces combined) Berry curvature induced by torsion in Weyl systems. We describe how torsion generates the chiral chemical potential and thus leads to the chiral magnetic effect (CME) in the presence of a background magnetic field. We propose to measure the CME current induced by the torsion as a way to detect the anomalous coupling between the coordinate-space and momentum-space chiralities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Perturbation theory and thermal transport in mass-disordered alloys: Insights from Green's function methods

Lowest-order quantum perturbation theory (Fermi's golden rule) for phonon-disorder scattering has been used to predict thermal conductivities in several semiconducting alloys with surprising success given its underlying hypothesis of weak and dilute disorder. In this paper, we explain how this is possible by focusing on the case of maximally mass-disordered Mg 2 Si 1-x Sn x . In this work, we use a Chebyshev polynomials Green's function method, which allows a full treatment of disorder on very large systems (tens of millions of atoms) to probe individual phonon linewidths and frequency-resolved thermal transport. We demonstrate that the success of perturbation theory originates from the specific form of mass disorder terms in the phonon Green's function and from the interplay between anharmonic and disorder scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Many-Body Theory of Energy Transport and Conversion at the Nanoscale

This project addressed the fundamental challenge of understanding interacting quantum systems far from equilibrium, while simultaneously exploring the potential of nanostructured materials for applications in energy-conversion technologies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Phonon Thermal Transport in UO 2 via Self-Consistent Perturbation Theory

Computing thermal transport from first-principles in UO 2 is complicated due to the challenges associated with Mott physics. Here, we use irreducible derivative approaches to compute the cubic and quartic phonon interactions in UO 2 from first principles, and we perform enhanced thermal transport computations by evaluating the phonon Green’s function via self-consistent diagrammatic perturbation theory. Our predicted phonon lifetimes at T = 600 K agree well with our inelastic neutron scattering measurements across the entire Brillouin zone, and our thermal conductivity predictions agree well with previous measurements. Finally, both the changes due to thermal expansion and self-consistent contributions are nontrivial at high temperatures, though the effects tend to cancel, and interband transitions yield a substantial contribution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quasilinear theory and modelling of gyrokinetic turbulent transport in tokamaks

The theory, development, and validation of reduced quasilinear models of gyrokinetic turbulent transport in the closed flux surface core of tokamaks is reviewed. In combination with neoclassical collisional transport, these models are successful in accurately predicting core tokamak plasma temperature, density, rotation, and impurity profiles in a variety of confinement regimes. Refined experimental tests have been performed to validate the predictions of the quasilinear models, probing changes in the dominant gyrokinetic instabilities, as reflected in fluctuation measurements, cross-phases, and transport properties. These tests continue to produce a deeper understanding of the complex mix of instabilities at both electron and ion gyroradius scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytic insights into nonlocal energy transport. III. steady state Fokker Planck theory in spherical and planar geometry

This paper develops an approximate steady state Fokker Planck theory for energetic electron transport in a spherical laser fusion target. First, we apply the theory to only a small population of electrons, which is specified at the outset. Then, one determines the nonlocal electron energy flux of these electrons at a particular position. These energetic electrons may either come from the tail of the thermal distribution function or be generated by an instability. This paper develops two approximate methods of solution, which we call the “characteristic method” and “sparse eigenfunction.” The former works only in planar geometry and the latter in both planar and spherical geometry. Comparison of the two methods in planar geometry shows that even though the approximations are very different, they give about the same result, increasing their credibility. For the example we have chosen, spherical effects are not important for electrons from the tail of the distribution function but may well be for instability generated electrons, which have much higher energy. Comparing planar to spherical, one finds an additional spherical barrier protecting the fuel. It turns out that the associated fuel preheat a Fokker Planck model predicts is considerably less than that predicted by the Krook models as developed at both NRL and other places.

Manheimer, Wallace (ORCID:0000000163342591)↗
DOE OSTI2020

Vibrational hierarchy leads to dual-phonon transport in low thermal conductivity crystals

Many low-thermal-conductivity ( κ L ) crystals show intriguing temperature ( T ) dependence of κ L : κ L ∝ T -1 (crystal-like) at intermediate temperatures whereas weak T -dependence (glass-like) at high temperatures. It has been in debate whether thermal transport can still be described by phonons at the Ioffe-Regel limit. In this work, we propose that most phonons are still well defined for thermal transport, whereas they carry heat via dual channels: normal phonons described by the Boltzmann transport equation theory, and diffuson-like phonons described by the diffusion theory. Three physics-based criteria are incorporated into first-principles calculations to judge mode-by-mode between the two phonon channels. Case studies on La 2 Zr 2 O 7 and Tl 3 VSe 4 show that normal phonons dominate low temperatures while diffuson-like phonons dominate high temperatures. Our present dual-phonon theory enlightens the physics of hierarchical phonon transport as approaching the Ioffe-Regel limit and provides a numerical method that should be practically applicable to many materials with vibrational hierarchy.

42 ENGINEERING↗
DOE OSTI2020

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗
DOE OSTI2021

Exponential Time Differencing Schemes for Fuel Depletion and Transport in Molten Salt Reactors: Theory and Implementation

A numerical framework for modeling depletion and mass transport in liquid-fueled molten salt reactions is presented based on exponential time differencing. The solution method involves using the finite volume method to transform the system of partial differential equations (PDEs) into a much larger system of ordinary differential equations. The key part of this method involves solving for the exponential of a matrix. We explore six different algorithms to compute the exponential in a series of progression problems that explore physical transport phenomena in molten salt reactors. This framework shows good results for solving linear parabolic PDEs with each of the six matrix exponential algorithms. For large problems, the series solvers such as Padé and Taylor have large run times, which can be mitigated by using the Krylov subspace.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Transport model comparison studies of intermediate-energy heavy-ion collisions

Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, which test separately selected ingredients of a transport code, and more realistic calculations of heavy-ion collisions. Over the years, six studies have been performed within this project. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. After a compact description of the underlying transport approaches, we review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann–Uehling–Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. Furthermore, our evaluation studies often indicate improved strategies in performing transport simulations and thus can provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗