Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gradient flow”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Strong gradient effects on neoclassical electron transport and the bootstrap current in large aspect ratio tokamaks

Standard approaches to neoclassical theory do not extend into regions of strong gradients in tokamaks such as the pedestal and internal transport barriers. Here, we calculate the modifications to neoclassical electron physics inside strong gradient regions of large aspect ratio tokamaks in the banana regime. We show that these modifications are due to the different ion flow and the strong poloidal variation of the potential. We also provide a physical interpretation of the mechanisms that drive poloidal asymmetries and hence a poloidal electric field. We apply our model to two specific example cases of pedestal profiles, calculating the neoclassical electron flux and the bootstrap current. We find that, depending on the ion flow, weak gradient neoclassical theory overestimates or underestimates the neoclassical electron transport and the bootstrap current in regions with strong gradients. We show that the determination of the mean parallel flow is more complex than in weak gradient neoclassical theory. For vanishing turbulence, we can determine the radial electric field for a given flow profile in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗
DOE OSTI2018

Deriving GENERIC from a Generalized Fluctuation Symmetry

Much of the structure of macroscopic evolution equations for relaxation to equilibrium can be derived from symmetries in the dynamical fluctuations around the most typical trajectory. For example, detailed balance as expressed in terms of the Lagrangian for the path-space action leads to gradient zero-cost flow. We expose a new such fluctuation symmetry that implies GENERIC, an extension of gradient flow where a Hamiltonian part is added to the dissipative term in such a way as to retain the free energy as Lyapunov function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗
DOE OSTI2017

Numerical and Experimental Investigation of Turbulent Transport Control via Shaping of Radial Plasma Flow Profiles

Turbulence, and turbulence-driven transport are ubiquitous in magnetically confined plasmas, where there is an intimate relationship between turbulence, transport, instability driving mechanisms (such as gradients), plasma flows, and flow shear. Though many of the detailed physics of the interrelationship between turbulence, transport, drive mechanisms, and flow remain unclear, there have been many demonstrations that transport and/or turbulence can be suppressed or reduced via manipulations of plasma flow profiles. This is well known in magnetic fusion plasmas [e.g., high confinement mode (H-mode) and internal transport barriers (ITB’s)], and has also been demonstrated in laboratory plasmas. However, it may be that the levels of particle transport obtained in such cases [e.g. H-mode, ITB’s] are actually lower than is desirable for a practical fusion device. Ideally, one would be able to actively feedback control the turbulent transport, via manipulation of the flow profiles. The purpose of this research was to investigate the feasibility of using both advanced model-based control algorithms, as well as non-model-based algorithms, to control cross-field turbulence-driven particle transport through appropriate manipulation of radial plasma flow profiles. The University of New Mexico was responsible for the experimental portion of the project, while our collaborators at the University of Montana provided plasma transport modeling, and collaborators at Lehigh University developed and explored control methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Non-modal theory of the kinetic ion temperature gradient driven instability of plasma shear flows across the magnetic field

The temporal evolution of the kinetic ion temperature gradient driven instability and of the related anomalous transport of the ion thermal energy of plasma shear flow across the magnetic field is investigated analytically. This instability develops in a steady plasma due to the inverse ion Landau damping and has the growth rate of the order of the frequency when the ion temperature is equal to or above the electron temperature. The investigation is performed employing the non-modal methodology of the shearing modes which are the waves that have a static spatial structure in the frame of the background flow. The solution of the governing linear integral equation for the perturbed potential displays that the instability experiences the non-modal temporal evolution in the shearing flow during which the unstable perturbation becomes very different from a canonical modal form. It transforms into the non-modal structure with vanishing frequency and growth rate with time. The obtained solution of the nonlinear integral equation, which accounts for the random scattering of the angle of the ion gyro-motion due to the interaction of ions with ensemble of shearing waves, reveals similar but accelerated process of the transformations of the perturbations into the zero frequency structures. It was obtained that in the shear flow the anomalous ion thermal conductivity decays with time. It is a strictly non-modal effect, which originates from the temporal evolution of the shearing modes turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Reconstruction of three-dimensional anisotropic structure from small-angle scattering experiments

When subjected to flow, the structures of many soft-matter systems become anisotropic due to the symmetry breaking of the spatial arrangements of constituent particles at the microscopic level. At present, it is common practice to use various small-angle scattering techniques to explore flow-induced microstructural distortion. However, there has not been a thorough discussion in the literature on how a three-dimensional anisotropic structure can be faithfully reconstructed from two-dimensional small-angle scattering spectra. Here, in this work, we address this issue rigorously from a mathematical perspective by using real spherical harmonic expansion analysis. We first show that, except for cases in which mechanical perturbation is sufficiently small, the existing small-angle scattering techniques generally do not provide complete information on structural distortion. This limitation is caused by the linear dependence of certain real spherical harmonic basis vectors on the flow-vorticity and flow-velocity gradient planes in the Couette shear cell. To circumvent the constraint imposed by this geometry, an alternative approach is proposed in which a parallel sliding plate shear cell is used with a central rotary axis along the flow direction. Finally, from the calculation of rotation of the reference frame, we demonstrate the feasibility of this experimental implementation for a fully resolved three-dimensional anisotropic structure via a case study of sheared polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Particle-in-cell study of the ion-to-electron sheath transition

The form of a sheath near a small electrode, with bias changing from below to above the plasma potential, is studied using 2D particle-in-cell simulations. When the electrode is biased within T e /2 e below the plasma potential, the electron velocity distribution functions (EVDFs) exhibit a loss-cone type truncation due to fast electrons overcoming the small potential difference between the electrode and plasma. No sheath is present in this regime, and the plasma remains quasineutral up to the electrode. The EVDF truncation leads to a presheath-like density and flow velocity gradients. Once the bias exceeds the plasma potential, an electron sheath is present. In this case, the truncation driven behavior persists, but is accompanied by a shift in the maximum value of the EVDF that is not present in the negative bias cases. In conclusion, the flow moment has significant contributions from both the flow shift of the EVDF maximum, and the loss-cone truncation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Density Effects on Lock Exchange Flow During a DCC

The High Temperature Test Facility (HTTF) is a scaled, electric model at Oregon State University (OSU) of the Modular High Temperature Gas Reactor (MHTGR) by General Atomics. There are a few major differences between the design and the OSU facility. In the HTTF the reactor pressure vessel and the reactor cavity are in two separate tanks connected by a crossover duct similar to the crossover duct that would be seen between the reactor core and the steam generator in the MHTGR. This feature allows for the facility to open two main break valves simulating a Double-Ended Guillotine Break (DEGB) of the cross over duct, synonymous to a crossover duct break in the MHTGR. The density gradient driven flow phenomena in the HTTF duct can then be studied further. This break would cause a depressurization of the reactor and a loss of flow. It is called the Depressurized Loss of Forced Convection Accident (D-LOFC), also known as the Depressurized Conduction Cool-Down Event for High Temperature Gas Reactors (HTGRs) because heat removal is by conduction through the core and radiation out of the vessel. The HTTF was built and scaled to maintain the highest accuracy for the DCC. Prototypical temperatures with an inlet temperature of 259 deg. C and an outlet temperature of 687 deg. C are maintained in the HTTF. The temperature profile is also maintained in the smaller volume by using a ceramic core. The HTTF uses graphite heaters because the resistive heat produced does not change significantly with temperature . In order to produce a facility of similar temperatures within the allotted budget, sacrifices had to be made in other areas. The HTTF is scaled one-fourth in diameter, one-fourth in height and one-eighth in pressure of the MHTGR. There are also a few limitations for the HTTF. While in the MHTGR, the cavity is filled with air, the HTTF uses highly pure Nitrogen and does not have the ability to simulate a full depressurization. The DCC event is traditionally comprised of multiple stages: depressurization, molecular diffusion and natural circulation. However, this traditional interpretation ignores lock exchange flow which will occur immediately following depressurization and could have significant consequences in the reactor. It was found that the traditional model underestimates the amount of air ingress by at least 660 times according to Oh and Kim in their theoretical research of air ingress. Molecular diffusion also decreases with stratification. Any model that specifies the gases are stratified has further underestimated the air ingress of the DCC event. In addition, molecular diffusion is a much slower process, allowing the reactor to shut down and cool off considerably from operating temperatures. Thus, understanding lock exchange flow, the air ingress concentrations and the time scale is critical to MHTGR research and licensing. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗
DOE OSTI2017

Turbulent transport and mixing in transitional Rayleigh-Taylor unstable flow: A priori assessment of gradient-diffusion and similarity modeling

Data from a 1152×760×1280 direct numerical simulation of a Rayleigh-Taylor mixing layer modeled after a small-Atwood-number water-channel experiment is used to investigate the validity of gradient diffusion and similarity closures a priori. The budgets of the mean flow, turbulent kinetic energy, turbulent kinetic energy dissipation rate, heavy-fluid mass fraction variance, and heavy-fluid mass fraction variance dissipation rate transport equations across the mixing layer were previously analyzed at different evolution times to identify the most important transport and mixing mechanisms. Here a methodology is introduced to systematically estimate model coefficients as a function of time in the closures of the dynamically significant terms in the transport equations by minimizing the L2 norm of the difference between the model and correlations constructed using the simulation data. It is shown that gradient-diffusion and similarity closures used for the turbulent kinetic energy K, turbulent kinetic energy dissipation rate ε, heavy-fluid mass fraction variance S, and heavy-fluid mass fraction variance dissipation rate χ equations capture the shape of the exact, unclosed profiles well over the nonlinear and turbulent evolution regimes. Using order-of-magnitude estimates for the terms in the exact transport equations and their closure models, it is shown that several of the standard closures for the turbulent production and dissipation (destruction) must be modified to include Reynolds-number scalings appropriate for Rayleigh-Taylor flow at small to intermediate Reynolds numbers. The late-time, large Reynolds number coefficients are determined to be different from those used in shear flow applications and largely adopted in two-equation Reynolds-averaged Navier-Stokes (RANS) models of Rayleigh-Taylor turbulent mixing. In addition, it is shown that the predictions of the Boussinesq model for the Reynolds stress agree better with the data when additional buoyancy-related terms are included. It is shown that an unsteady RANS paradigm is needed to predict the transitional flow dynamics from early evolution times, analogous to the small Reynolds number modifications in RANS models of wall-bounded flows in which the production-to-dissipation ratio is far from equilibrium. Although the present study is specific to one particular flow and one set of initial conditions, the methodology could be applied to calibrations of other Rayleigh-Taylor flows with different initial conditions (which may give different results during the early-time, transitional flow stages, and perhaps asymptotic stage). Finally, the implications of these findings for developing high-fidelity eddy viscosity-based turbulent transport and mixing models of Rayleigh-Taylor turbulence are discussed.

42 ENGINEERING↗
DOE OSTI2024

Large eddy simulation of lean blow-off in swirl-stabilized flame with the subgrid dissipation concept

The lean blow-off mechanism of the premixed swirl flame is numerically investigated by large eddy simulation (LES) with the subgrid dissipation concept (SDC) combustion model. Three simulated cases cover stable, near blow-off, and transient conditions. Compared with the experiment, the LES-SDC approach captures the flow and combustion features for stable and near blow-off conditions. More importantly, the predictions of the blowoff procedure and duration agree satisfactorily with the experiment, indicating that the LES-SDC approach is a promising tool for predicting strong, unsteady turbulent combustion processes. Further, the numerical results are used to investigate the blow-off mechanism. Two stages in the blow-off procedure are specified. The first is the necking and extinction of the downstream flame surface, and the second is the shrinking of the upstream flame surface. The blow-off mechanism is well explained by the theory of stretched flame extinction. At the end of the recirculation zone, the large negative radial velocity pushes the flame to the central line. The combustion process here can be abstracted as the stretched counter-flame of the reactant-reactant configuration. The excessive flame stretch dominates the flame extinction and triggers the blow-off event. The upstream flame resists the intense stretch with the help of hot product recirculation, and the combustion here can be idealized as the counter-flame of the reactant-product configuration. The alignment of the temperature gradient and flow velocity, together with the excessive stretch, clearly indicates the tendency of flame local extinction. A Damköhler number-based blow-off criterion is raised from the mean flow strain rate and laminar flame bulk extinction strain rate.

42 ENGINEERING↗
DOE OSTI2021

Optimizing $\mathrm{ACRT}$ to reduce inclusion formation during the $\mathrm{VGF}$ growth of cadmium zinc telluride: I. Computational approach

Theoretical computations reveal insight into parameters that govern the accelerated crucible rotation technique (ACRT) applied to the gradient freeze (VGF) growth of cadmium zinc telluride (CZT). A metric based on the classic Mullins and Sekerka instability criterion is put forth that provides a quantitative means of assessing the impact of different ACRT rotation schedules on morphological stability and inclusion formation during CZT growth. Here this metric is employed in a 2-k factorial design of experiments that identifies the maximum rotation rate as the most important ACRT parameter, followed by the acceleration period and the rest period. Subsequent optimization of rotation schedule is performed for a specific VGF growth system, and mechanistic explanations are provided for the preferred ACRT schedule, which improves interface stability via slower rotation rates and longer acceleration periods than would be applied using classical ACRT schedules.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optimizing $\mathrm{ACRT}$ to reduce inclusion formation during the $\mathrm{VGF}$ growth of cadmium zinc telluride: II. Application to experiments

Here two different accelerated crucible rotation technique (ACRT) rotation schedules are assessed, via theoretical computations and growth experiments, according to their ability to reduce tellurium-rich inclusions during the vertical gradient freeze (VGF) of cadmium zinc telluride (CZT). A vigorous rotation schedule based on classical ACRT guidelines produces a well-mixed melt that is expected to reduce tellurium inclusions compared to growth without crucible rotation. In contrast, a new ACRT schedule, derived from model-based optimization for interface stability, employs slower rotation with longer periods and is predicted to further reduce inclusions. Growth experiments corroborate these expectations. Namely, both ACRT schedules result in crystals with inclusion size and volume significantly decreased from levels found in material grown with no rotation, and material grown using the computationally optimized ACRT parameters exhibits a median inclusion size that is smaller and with a sharper distribution than in material grown via classical ACRT.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermo-economic analysis of a solid oxide fuel cell-gas turbine hybrid with commercial off-the-shelf gas turbine

The thermodynamic and economic performance of a natural gas-fueled solid oxide fuel cell (SOFC)-gas turbine (GT) hybrid system with a commercial off-the-shelf GT is investigated. Until today, commercial GTs are primarily engineered for the direct use of energy dense fuels, such as natural gas, and the integration of commercial GTs into a SOFC-GT hybrid remains challenging. In this work technically feasible and economically viable GT operating modes are identified to gauge the technology’s competitiveness on the free market. Steady state, full load operating conditions were established for various GT operating modes: I) constant spool speed operation, II) variable spool speed operation and III) partially closed compressor inlet guide vanes. For each GT operating mode, the performance was investigated over a range of SOFC fuel utilization factors, while considering physical constraints inside the SOFC, such as local temperature gradients in flow direction as well as overall cell temperature differences. The results of the thermodynamic evaluation served as inputs for the economic analysis of the SOFC-GT hybrid power plant. The results show that the integration of an off-the-shelf GT not only results in a significant derating of the GT, but also substantially impacts the SOFC operation, the main power producer in this SOFC-GT hybrid. To maximize the SOFC power output it is desirable operate the GT in a region of high air mass flow rates and low pressure ratios, which increases the number of stacks that can be accommodated and reduces the SOFC cooling requirement. The lowest costs of electricity are obtained at constant spool speed operation, while the highest efficiencies are reached at variable spool speed operation. Critical for the integration of off-the-shelf GTs, which historically have been designed for natural gas, is the surge margin. The largest surge margins are obtained by closing the compressor inlet guide vanes. Furthermore, higher SOFC fuel utilization factors are shown to increase the surge margin as the turbine firing temperature is decreased.

25 ENERGY STORAGE↗
DOE OSTI2021

Deep Kronecker neural networks: A general framework for neural networks with adaptive activation functions

Here we propose a new type of neural networks, Kronecker neural networks (KNNs), that form a general framework for neural networks with adaptive activation functions. KNNs employ the Kronecker product, which provides an efficient way of constructing a very wide network while keeping the number of parameters low. Our theoretical analysis reveals that under suitable conditions, KNNs induce a faster decay of the loss than that by the feed-forward networks. This is also empirically verified through a set of computational examples. Furthermore, under certain technical assumptions, we establish global convergence of gradient descent for KNNs. As a specific case, we propose the Rowdy activation function that is designed to get rid of any saturation region by injecting sinusoidal fluctuations, which include trainable parameters. The proposed Rowdy activation function can be employed in any neural network architecture like feed-forward neural networks, Recurrent neural networks, Convolutional neural networks etc. The effectiveness of KNNs with Rowdy activation is demonstrated through various computational experiments including function approximation using feed-forward neural networks, solution inference of partial differential equations using the physics-informed neural networks, and standard deep learning benchmark problems using convolutional and fully-connected neural networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Ion parallel closures

Ion parallel closures are obtained for arbitrary atomic weights and charge numbers. For arbitrary collisionality, the heat flow and viscosity are expressed as kernel-weighted integrals of the temperature and flow-velocity gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗