Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The Fornax Deep Survey with VST. I. The extended and diffuse stellar halo of NGC 1399 out to 192 kpc

We have started a new, deep multi-imaging survey of the Fornax cluster, dubbed the Fornax Deep Survey (FDS), at the VLT Survey Telescope (VST). In this paper we present the deep photometry inside two square degrees around the bright galaxy NGC 1399 in the core of the cluster. We found that the core of the Fornax cluster is characterized by a very extended and diffuse envelope surrounding the luminous galaxy NGC 1399: we map the surface brightness out to 33 arcmin (∼192 kpc) from the galaxy center and down to μ{sub g} ∼ 31 mag arcsec{sup −2} in the g band. The deep photometry allows us to detect a faint stellar bridge in the intracluster region on the west side of NGC 1399 and toward NGC 1387. By analyzing the integrated colors of this feature, we argue that it could be due to the ongoing interaction between the two galaxies, where the outer envelope of NGC 1387 on its east side is stripped away. By fitting the light profile, we found that there exists a physical break radius in the total light distribution at R = 10 arcmin (∼58 kpc) that sets the transition region between the bright central galaxy and the outer exponential halo, and that the stellar halo contributes 60% of the total light of the galaxy (Section 3.5). We discuss the main implications of this work on the build-up of the stellar halo at the center of the Fornax cluster. By comparing with the numerical simulations of the stellar halo formation for the most massive bright cluster galaxies (i.e., 13<logM{sub 200}/M{sub ⊙}<14), we find that the observed stellar halo mass fraction is consistent with a halo formed through the multiple accretion of progenitors with stellar mass in the range 10{sup 8}–10{sup 11} M{sub ⊙}. This might suggest that the halo of NGC 1399 has also gone through a major merging event. The absence of a significant number of luminous stellar streams and tidal tails out to 192 kpc suggests that the epoch of this strong interaction goes back to an early formation epoch. Therefore, different from the Virgo cluster, the extended stellar halo around NGC 1399 is characterized by a more diffuse and well-mixed component, including the intracluster light.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

DNS of a turbulent lifted DME jet flame

A three-dimensional direct numerical simulation (DNS) of a turbulent lifted dimethyl ether (DME) slot jet flame was performed at elevated pressure to study interactions between chemical reactions with low-temperature heat release (LTHR), negative temperature coefficient (NTC) reactions and shear generated turbulence in a jet in a heated coflow. By conditioning on mixture fraction, local reaction zones and local heat release rate, the turbulent flame is revealed to exhibit a “pentabrachial” structure that was observed for a laminar DME lifted flame [Krisman et al., (2015)]. The propagation characteristics of the stabilization and triple points are also investigated. Potential stabilization points, spatial locations characterized by preferred temperature and mixture fraction conditions, exhibit autoignition characteristics with large reaction rate and negligible molecular diffusion. The actual stabilization point which coincides with the most upstream samples from the pool of potential stabilization points fovr each spanwise location shows passive flame structure with large diffusion. The propagation speed along the stoichiometric surface near the triple point is compared with the asymptotic value obtained from theory [Ruetsch et al., (1995)]. At stoichiometric conditions, the asymptotic and averaged DNS values of flame displacement speed deviate by a factor of 1.7. However, accounting for the effect of low-temperature species on the local flame speed increase, these two values become comparable. In conclusion, this suggests that the two-stage ignition influences the triple point propagation speed through enhancement of the laminar flame speed in a configuration where abundant low-temperature products from the first stage, low-temperature ignition are transported to the lifted flame by the high-velocity jet.

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

A Fokker–Planck approach to a moment closure for mixing in variable-density turbulence

We develop a theory for the cascade mixing terms in a moment closure approach to binary active scalar mixing in variable-density turbulence. To address the variable-density complications we apply, as a principle and constraint, the conservation of the probability density function (PDF) through a Fokker–Planck equation with bounded sample space whose attractor is the beta PDF with skewness. Mixing is related to a single-point PDF as a realisability principle to provide mathematically rigorous expressions for the small scale statistics in terms of largescale moments. The problem of the unknown small-scale mixing is replaced with the determination of the drift and diffusion terms of a Fokker–Planck equation in a beta-PDF-convergent stochastic process. We find that realisability of a beta-convergent process requires the mixing time-scale ratio, taken as a constant in passive scalar mixing, to be a function of the mean mass fraction, mean fluid density, the Atwood number, the density-volume correlation and moments of the density field. We develop and compare the new model with direct numerical simulations data of non-stationary homogeneous variable-density turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

SOLAR WIND STRAHL BROADENING BY SELF-GENERATED PLASMA WAVES

This Letter reports on the results of numerical simulations which may provide a possible explanation for the strahl broadening during quiet solar conditions. The relevant processes involved in the broadening are due to kinetic quasi-linear wave-particle interaction. Making use of static analytical electron distribution in an inhomogeneous field, it is found that self-generated electrostatic waves at the plasma frequency, i.e., Langmuir waves, are capable of scattering the strahl component, resulting in energy and pitch-angle diffusion that broadens its velocity distribution significantly. The present theoretical results provide an alternative or complementary explanation to the usual whistler diffusion scenario, suggesting that self-induced electrostatic waves at the plasma frequency might play a key role in broadening the solar wind strahl during quiet solar conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air

Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.

42 ENGINEERING↗
DOE OSTI2012

Suprathermal ion transport in simple magnetized torus configurations

Inspired by suprathermal ion experiments in the basic plasma experiment TORPEX, the transport of suprathermal ions in ideal interchange mode turbulence is theoretically examined in the simple magnetized torus configuration. We follow ion tracer trajectories as specified by ideal interchange mode turbulence imported from a numerical simulation of drift-reduced Braginskii equations. Using the variance of displacements, {sigma}{sup 2}(t){approx}t{sup {gamma}}, we find that {gamma} depends strongly on suprathermal ion injection energy and the relative magnitude of turbulent fluctuations. The value of {gamma} also changes significantly as a function of time after injection, through three distinguishable phases: ballistic, interaction, and asymmetric. During the interaction phase, we find the remarkable presence of three regimes of dispersion: superdiffusive, diffusive, and subdiffusive, depending on the energy of the suprathermal ions and the amplitude of the turbulent fluctuations. We contrast these results with those from a 'slab' magnetic geometry in which subdiffusion does not occur during the interaction phase. Initial results from TORPEX are consistent with data from a new synthetic diagnostic used to interpret our simulation results. The simplicity of the simple magnetized torus makes the present work of interest to analyses of more complicated contexts ranging from fusion devices to astrophysics and space plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Mathematical model of testing of pipeline integrity by thermal fields

Thermal fields testing at the ground surface above a pipeline are considered. One method to obtain and investigate an ideal thermal field in different environments is a direct numerical simulation of heat transfer processes taking into account the most important physical factors. In the paper a mathematical model of heat propagation from an underground source is described with accounting of physical factors such as filtration of water in soil and solar radiation. Thermal processes are considered in 3D origin where the heat source is a pipeline with constant temperature and non-uniform isolated shell (with 'damages'). This problem leads to solution of heat diffusivity equation with nonlinear boundary conditions. Approaches to analysis of thermal fields are considered to detect damages.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Magnetic relaxation, current sheets, and structure formation in an extremely Tenuous fluid medium

The process of relaxation of a unidirectional magnetic field in a highly conducting tenuous fluid medium is considered. Null points of the field play a critical role in this process. During an initial stage of relaxation, variations in magnetic pressure are eliminated, and current sheets build up in the immediate neighborhood of null points. This initial phase is followed by a long diffusive phase of slow algebraic decay of the field, during which fluid is continuously sucked into the current sheets, leading to exponential growth of fluid density and concentration of mass around the null points, which show a tendency to cluster. Ultimately, this second phase of algebraic decay gives way to a final period of exponential decay of the field. The peaks of density at the null points survive as a fossil relic of the decay process. Numerical solution of the governing equations provides convincing confirmation of this three-stage scenario. Generalizations to two- and three-dimensional fields are briefly considered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Li1−xNiO2 Many-body DMC Benchmark Dataset

The dataset contains all numerical data generated in support of the manuscript “Many‑body Benchmark of Electronic Charge and Spin Densities for Li1–xNiO2​” (Journal of Chemical Theory and Computation, DOI: 10.1021/acs.jctc.5c02097, URL: https://pubs.acs.org/doi/10.1021/acs.jctc.5c02097). The materials included in this repository are: 1. Data files used to produce all figures and tables in the main manuscript and supporting information. 2. Benchmark density‑functional theory (DFT) datasets used for the charge‑ and spin‑density analyses. 3. Reference many‑body diffusion Monte Carlo (DMC) calculations and associated input/output files.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Experimental Investigation of Forced Convection and Natural Circulation Cooling of a VHTR Core under Normal Operation and Accident Scenarios

This NEUP project #15-8205 entitled “Experimental Investigation of Forced Convection and Natural Circulation Cooling of a VHTR Core under Normal Operation and Accident Scenarios”, has investigated important phenomena in the safety of Very High Temperature Reactors (VHTRs) using the experimental facilities at City College of New York (CCNY) and Kansas State University (KSU), in collaboration with the Idaho State University (ISU) and Idaho National Laboratory (INL). The main objectives of this project were to investigate the forced convection and bypass flow phenomena, natural circulation flow and heat transfer, and graphite oxidation due to air ingress which could occur in a Very High Temperature Reactor (VHTR) with a prismatic core. High pressure/high temperature experiments have been conducted to obtain data that could be used for validation of VHTR design and safety analysis codes. The focus of these experiments was on the generation of benchmark data for design and off-design heat transfer for forced and natural circulation in a VHTR core. This work has built on the previous NEUP project #11-3218 which investigated forced convection and flow laminarization phenomena in a prismatic core of a VHTR and natural circulation of a pure gas under Pressurized Conduction Cooldown (PCC) and Depressurized Conduction Cooldown (DCC) conditions. Under each Task, extensive literature reviews on natural circulation, bypass flow, air ingress and graphite oxidation phenomena were conducted. In Task 1, scaling analyses were conducted by Idaho State University followed by design and construction of experimental facilities needed to conduct both bypass flow and natural circulation experiments at CCNY. Natural circulation experiments were conducted in Task 2 to investigate the air ingress phenomena, and bypass flow phenomena in Task 3, both at CCNY. In Task 4, radiation heat transfer and graphite oxidation tests were performed at Kansas State University. The experimental data were then analyzed to obtain quantitative results on natural circulation of pure gases and gas mixtures under DCC conditions, bypass flow rates, radiation heat transfer and graphite oxidation. In addition, CFD and Direct Numerical Simulations were conducted by INL and CCNY to investigate the flow laminarization phenomenon in strongly heated gas flows. The scaling analyses have shown that the key nondimensional numbers and parameters all have acceptable distortion levels principally because the CCNY facility’s test sections have representative cooling channel dimensions and are using either helium or a scalable gas such as nitrogen at representative state conditions. The natural circulation experiments involving helium and nitrogen gas mixtures at CCNY have elucidated the effect of nitrogen (simulating air) ingress into a lower plenum on the onset of natural circulation, mass flow rates and heat transfer coefficients under different power input levels in the riser section. Also, the transport of nitrogen from the lower plenum to upper plenum due to natural circulation was quantified for the first time. Nusselt number correlations for the heat transfer coefficient under natural circulation as a function of Reynolds and Prandtl numbers was also developed. The experiments at KSU revealed that prior to the onset of natural circulationin a VHTR core. This work has built on the previous NEUP project #11-3218 which investigated forced convection and flow laminarization phenomena in a prismatic core of a VHTR and natural circulation of a pure gas under Pressurized Conduction Cooldown (PCC) and Depressurized Conduction Cooldown (DCC) conditions. Under each Task, extensive literature reviews on natural circulation, bypass flow, air ingress and graphite oxidation phenomena were conducted. In Task 1, scaling analyses were conducted by Idaho State University followed by design and construction of experimental facilities needed to conduct both bypass flow and natural circulation experiments at CCNY. Natural circulation experiments were conducted in Task 2 to investigate the air ingress phenomena, and bypass flow phenomena in Task 3, both at CCNY. In Task 4, radiation heat transfer and graphite oxidation tests were performed at Kansas State University. The experimental data were then analyzed to obtain quantitative results on natural circulation of pure gases and gas mixtures under DCC conditions, bypass flow rates, radiation heat transfer and graphite oxidation. In addition, CFD and Direct Numerical Simulations were conducted by INL and CCNY to investigate the flow laminarization phenomenon in strongly heated gas flows. The scaling analyses have shown that the key nondimensional numbers and parameters all have acceptable distortion levels principally because the CCNY facility’s test sections have representative cooling channel dimensions and are using either helium or a scalable gas such as nitrogen at representative state conditions.The natural circulation experiments involving helium and nitrogen gas mixtures at CCNY have elucidated the effect of nitrogen (simulating air) ingress into a lower plenum on the onset of natural circulation, mass flow rates and heat transfer coefficients under different power input levels in the riser section. Also, the transport of nitrogen from the lower plenum to upper plenum due to natural circulation was quantified for the first time. Nusselt number correlations for the heat transfer coefficient under natural circulation as a function of Reynolds and Prandtl numbers was also developed. The experiments at KSU revealed that prior to the onset of natural circulation(ONC), molecular diffusion plays a significant role in air-ingress, however, at higher temperatures convection currents may be influencing air ingress and in-turn ONC times. Bypass flow experiments faced difficulties in the measurement of gas velocity in the main channel due to the failure of a high temperature Hot Wire Anemometer probe, so an energy balance method was instead employed to estimate the flow distribution between the main and bypass flow channels. In the graphite oxidation experiments conducted at KSU, oxidation was found to increase the thermal diffusivity of nuclear grade graphite. This change is believed to be driven by oxidation decreasing graphite sample density rather than raising the thermal conductivity. Oxidation was also shown to greatly alter graphite emissivity. If large amounts of impurities have built up within the graphite over time and have significantly increased its effective ash content, then air ingress is likely to severely inhibit the radiative heat transfer abilities of the graphite due to the formation of low emissivity, ash layers on its surface. However, if the purity of the graphite at the time of the air ingress is still close to that of its virgin specifications, the oxidation induced roughening of the graphite’s outer surface would have the effect of slightly raising its emissivity. The flow laminarization phenomenon was investigated by INL by examining the results of Direct Numerical Simulations of previous experiments, and by CCNY conducting its own DNS simulations using DOE’s Nek5000 code

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗
DOE OSTI2021

Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times

Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.

42 ENGINEERING↗
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 OSTI2012

DIRECT NUMERICAL SIMULATION OF RADIATION PRESSURE-DRIVEN TURBULENCE AND WINDS IN STAR CLUSTERS AND GALACTIC DISKS

The pressure exerted by the radiation of young stars may be an important feedback mechanism that drives turbulence and winds in forming star clusters and the disks of starburst galaxies. However, there is great uncertainty in how efficiently radiation couples to matter in these high optical depth environments. In particular, it is unclear what levels of turbulence the radiation can produce, and whether the infrared radiation trapped by the dust opacity can give rise to heavily mass-loaded winds. In this paper, we report a series of two-dimensional flux-limited diffusion radiation-hydrodynamics calculations performed with the code ORION in which we drive strong radiation fluxes through columns of dusty matter confined by gravity in order to answer these questions. We consider both systems where the radiation flux is sub-Eddington throughout the gas column, and those where it is super-Eddington at the midplane but sub-Eddington in the atmosphere. In the latter, we find that the radiation-matter interaction gives rise to radiation-driven Rayleigh-Taylor instability, which drives supersonic turbulence at a level sufficient to fully explain the turbulence seen in Galactic protocluster gas clouds, and to make a non-trivial contribution to the turbulence observed in starburst galaxy disks. However, the instability also produces a channel structure in which the radiation-matter interaction is reduced compared to time-steady analytic models because the radiation field is not fully trapped. For astrophysical parameters relevant to forming star clusters and starburst galaxies, we find that this effect reduces the net momentum deposition rate in the dusty gas by a factor of {approx}2-6 compared to simple analytic estimates, and that in steady state the Eddington ratio reaches unity and there are no strong winds. We provide an approximation formula, appropriate for implementation in analytic models and non-radiative simulations, for the force exerted by the infrared radiation field in this regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Transport and optical conductivity in the Hubbard model: A high-temperature expansion perspective

We derive analytical expressions for the spectral moments of the dynamical response functions of the Hubbard model using the high-temperature series expansion. We consider generic dimension d as well as the infinite- d limit, arbitrary electron density n, and both finite and infinite repulsion U. We use moment-reconstruction methods to obtain the one-electron spectral function, the self-energy, and the optical conductivity. They are all smooth functions at high temperature and, at large U, they are featureless with characteristic widths of the order of the lattice hopping parameter t. In the infinite-d limit, we compare the series expansion results with accurate numerical renormalization group and interaction expansion quantum Monte Carlo results. We find excellent agreement down to surprisingly low temperatures, throughout most of the bad-metal regime, which applies for T≳ (1 - n)D, the Brinkman-Rice scale. The resistivity increases linearly in T at high temperature without saturation. This results from the 1/T behavior of the compressibility or kinetic energy, which play the role of the effective carrier number. In contrast, the scattering time (or diffusion constant) saturates at high T. We find that σ( n,T) ≈ (1 - n)σ(n = 0,T) to a very good approximation for all n, with σ(n = 0,T) ∝ t/T at high temperatures. The saturation at small n occurs due to a compensation between the density dependence of the effective number of carriers and that of the scattering time. The T dependence of the resistivity displays a kneelike feature which signals a crossover to the intermediate-temperature regime where the diffusion constant (or scattering time) starts increasing with decreasing T. At high temperatures, the thermopower obeys the Heikes formula, while the Wiedemann-Franz law is violated with the Lorenz number vanishing as 1/T 2 . The relevance of our calculations to experiments probing high-temperature transport in materials with strong electronic correlations or ultracold atomic gases in an optical lattice is briefly discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗
DOE OSTI2020

Multigrid Algorithms with Projection and Prolongation over Elements of the Phase Space for K-Eigenvalue Transport Problems

This paper describes new multilevel acceleration methods for solving the multigroup neu- tron transport eigenvalue problems. These multilevel algorithms use different projection and prolongation operators in the phase space. The Nonlinear Diffusion Acceleration (NDA) method with multiple grids in energy is formulated with the prolongation oper- ator based on multiplication iterative correction and linear-in-energy mapping. Another multilevel NDA method uses the projection operator with coarsening in energy between the high-order transport and low-order NDA equations. The third algorithm is formu- lated with the partial-current based CMFD low-order equations and applies projection operators in space and energy. The numerical results are presented.

Cornejo, Luke↗
DOE OSTI2013

A 2D/1D coupling neutron transport method based on the matrix MOC and NEM methods

A new 2D/1D coupling method based on the matrix MOC method (MMOC) and nodal expansion method (NEM) is proposed for solving the three-dimensional heterogeneous neutron transport problem. The MMOC method, used for radial two-dimensional calculation, constructs a response matrix between source and flux with only one sweep and then solves the linear system by using the restarted GMRES algorithm instead of the traditional trajectory sweeping process during within-group iteration for angular flux update. Long characteristics are generated by using the customization of commercial software AutoCAD. A one-dimensional diffusion calculation is carried out in the axial direction by employing the NEM method. The 2D and ID solutions are coupled through the transverse leakage items. The 3D CMFD method is used to ensure the global neutron balance and adjust the different convergence properties of the radial and axial solvers. A computational code is developed based on these theories. Two benchmarks are calculated to verify the coupling method and the code. It is observed that the corresponding numerical results agree well with references, which indicates that the new method is capable of solving the 3D heterogeneous neutron transport problem directly. (authors)

97 MATHEMATICS AND COMPUTING↗