Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient flow dynamics”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A dynamic subgrid-scale model for compressible turbulence and scalar transport

The dynamic subgrid-scale (SGS) model of Germano et al. (1991) is generalized for the large eddy simulation (LES) of compressible flows and transport of a scalar. The model was applied to the LES of decaying isotropic turbulence, and the results are in excellent agreement with experimental data and direct numerical simulations. The expression for the SGS turbulent Prandtl number was evaluated using direct numerical simulation (DNS) data in isotropic turbulence, homogeneous shear flow, and turbulent channel flow. The qualitative behavior of the model for turbulent Prandtl number and its dependence on molecular Prandtl number, direction of scalar gradient, and distance from the wall are in accordance with the total turbulent Prandtl number from the DNS data.

Moin, P.↗
NASA NTRS1995

Domain decomposition for aerodynamic and aeroacoustic analyses, and optimization

The overarching theme was the domain decomposition, which intended to improve the numerical solution technique for the partial differential equations at hand; in the present study, those that governed either the fluid flow, or the aeroacoustic wave propagation, or the sensitivity analysis for a gradient-based optimization. The role of the domain decomposition extended beyond the original impetus of discretizing geometrical complex regions or writing modular software for distributed-hardware computers. It induced function-space decompositions and operator decompositions that offered the valuable property of near independence of operator evaluation tasks. The objectives have gravitated about the extensions and implementations of either the previously developed or concurrently being developed methodologies: (1) aerodynamic sensitivity analysis with domain decomposition (SADD); (2) computational aeroacoustics of cavities; and (3) dynamic, multibody computational fluid dynamics using unstructured meshes.

Baysal, Oktay↗
DOE OSTI2022

Measurements of the imploding plasma sheath in triple-nozzle gas-puff z pinches

Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A regularized-interface method as a unified formulation for simulations of high-pressure multiphase flows

The injection of multi-species fluids into high-pressure and high-temperature environments beyond the species' critical points is commonly found in engineering applications. At these conditions, for immiscible species, both subcritical interfacial dynamics and supercritical mixing can coexist due to variations in temperature around the mixture critical point. The modeling of these complex transcritical phenomena for large-scale configurations is so far not possible. To address this issue, we propose the Regularized-Interface Method (RIM) as a unified formulation that can describe both sub- and supercritical processes as well as the transition between them. The proposed method is derived via filtering of the nanoscale interface-resolving formulation based on van der Waals' linear gradient theory. Thus, this approach allows for the consistent modeling of interfacial dynamics that vanishes at supercritical conditions, while significantly reducing the temporal and spatial resolution constraints of the original nanoscale formulation. The resulting RIM formulation is examined in interface-capturing simulations of sub-, trans-, and supercritical fuel injection processes, involving droplets and jets. Furthermore, these results highlight the importance of resolving spatio-temporal transitions from subcritical interfacial dynamics to supercritical mixing in high-pressure multiphase simulations, in contrast to commonly employed diffused-interface methods, where interfacial dynamics are often neglected.

Interface capturing↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗
DOE OSTI2024

ELM-free enhanced D α H-mode with near zero NBI torque injection in DIII-D tokamak

Enhanced D α H-mode (EDA H-mode), an ELM-free H-mode regime, and the concomitant quasi-coherent mode (QCM) are explored in neutral beam heated, lower single null plasmas with near zero external torque injection. This regime exhibits good energy confinement (H 98y2 ~ 1) with β N ~ 2, high density, regime access at low input power, and no ELMs. New time–resolved measurements of electron and ion density, temperature, plasma rotation, and radial electric field as well as turbulence, flows, thermal and particle transport, and linear instability calculations examining the dynamics of the QCM are presented. Measurements using multiple fluctuation diagnostics reveal the QCM to be a separatrix spanning mode, peaking just inside the separatrix, existing in a wide range of k ⊥ ρ s ~ 0.1–1.2 with multiple harmonics, and propagating with a very small phase velocity in the plasma frame, where k ⊥ is the binormal wavenumber and ρ s is the ion sound radius. Linear gyrokinetic simulations of this EDA H-mode discharge with CGYRO indicate that the trapped electron mode (TEM) and electron temperature gradient mode are the dominant instabilities in the radial region where QCM is also unstable. Qualitative analysis indicates that the properties of TEM are consistent with the experimentally observed characteristics of the QCM. These similarities suggest that the QCM is consistent with a TEM instability existing in the edge region of the EDA H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

High-fidelity Multidisciplinary Sensitivity Analysis and Design Optimization for Rotorcraft Applications

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. A sensitivity-enabled fluid dynamics solver and a nonlinear flexible multibody dynamics solver are coupled to predict aerodynamic loads and structural responses of helicopter rotor blades. A discretely consistent adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system is validated by conducting simulations for a benchmark rotorcraft model and comparing solutions with established analyses and experimental data. Sensitivities of lift computed by the multidisciplinary sensitivity analysis are verified by comparison with the sensitivities obtained by complex-variable simulations. Finally the multidisciplinary sensitivity analysis is applied to a constrained gradient-based design optimization for a HART-II rotorcraft configuration.

Wang, Li↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗
DOE OSTI2020

Steady detonation propagation in thin channels with strong confinement

We examine asymptotically the dynamics of two-dimensional, steady detonation wave propagation and failure for a strongly confined high explosive (HE), in which the width of the explosive is small relative to the reaction zone length. An energy balance equation is derived, which shows how the longitudinal acceleration of subsonic flow behind the detonation shock is influenced both by chemical reaction and by the effects of HE boundary streamline deflection, specifically via the induced rate of change of mass flux through the detonation wave. The latter serves to either counteract or reinforce the acceleration of longitudinal flow, depending on the sign of the gradient of the boundary streamline deflection at the detonation shock. The analysis is valid for general equations of state and chemical reaction rates in the HE. The asymptotically derived form of the energy equation represents an eigenvalue problem for the determination of the steady detonation propagation speed, solved via a shooting method. We explore specific results for ideal and stiffened equations of state, along with a pressure-dependent reaction rate for which changes in the pressure exponent and reaction order are also studied. We consider the influences of both straight and curved HE boundary streamline shapes. The asymptotic analysis reveals significant physical insights into how detonation propagation and failure are affected by strong confinement.

42 ENGINEERING↗
DOE OSTI2022

Examination of stiff ion temperature gradient mode physics in simulations of DIII-D H-mode transport

A systematic evaluation of gyrokinetic and gyrofluid model predictions of ion temperature gradient (ITG) stability and transport using parameters from DIII-D high confinement mode (H-mode) plasmas has been performed. The nonlinear CGYRO code is used to make the gyrokinetic predictions, and the quasilinear TGLF model for the corresponding gyrofluid predictions. The assessments are made at three radii (normalized toroidal flux ρtor = 0.4, 0.55, and 0.7) in three different plasma scenarios with varying levels of neutral beam heating and torque. For each of the nine cases (3 radii × 3 scenarios) considered, ITG turbulence is found to be the dominant long-wavelength instability and transport mechanism. The inclusions of both transverse magnetic fluctuations and dynamic fast beam ions are stabilizing for all cases considered, with strongest effects seen at ρor = 0.4 where the fast ion population and normalized plasma pressure β = 2μ0nT/B2 are highest. The further inclusion of parallel magnetic fluctuations does not have a meaningful impact on the ITG turbulence in these scenarios, but does destabilize (in combination with fast ions) new high-frequency instabilities at ρtor = 0.4 in the high power scenarios. In each case the linear and nonlinear ITG critical gradients are predicted to be lower than the measured ITG scale lengths and their associated uncertainties. Inclusion of equilibrium flow shear in the transport predictions generally leads to an upshift in effective critical gradient rather than a qualitative change in the predicted stiffness, with stronger responses typically seen in the gyrokinetic predictions than in the gyrofluid results. However, in most cases these upshifted gradients still remain below the measured values and their uncertainties. Although the predicted critical gradients are below the measured gradients, both models predicted flux-matching gradients consistent with measured values in six of the nine cases considered, with no clear systematic over- or underprediction. Thus, while the experimental ion temperature profiles do not appear to be closely pinned to the ITG critical gradient, both gyrokinetic and gyrofluid models are able to accurately match the measured gradients reasonably well in most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Impact-generated winds on Venus: Causes and effects

The pressure of the dense atmosphere of Venus significantly changes the appearance of ejecta deposits relative to craters on the Moon and Mercury. Conversely, specific styles and sequences of ejecta emplacement can be inferred to represent different intensities of atmospheric response winds acting over different timescales. Three characteristic timescales can be inferred from the geologic record: surface scouring and impactor-controlled (angle and direction) initiation of the long fluidized run-out flows; nonballistic emplacement of inner, radar-bright ejecta facies and radar-dark outer facies; and very late reworking of surface materials. These three timescales roughly correspond to processes observed in laboratory experiments that can be scaled to conditions on Venus (with appropriate assumptions): coupling between the atmosphere and earlytime vapor/melt (target and impactor) that produces an intense shock that subsequently evolves into blast/response winds; less energetic dynamic response of the atmosphere to the outward-moving ballistic ejecta curtain that generates nonthermal turbulent eddies; and late recovery of the atmosphere to impact-generated thermal and pressure gradients expressed as low-energy but long-lived winds. These different timescales and processes can be viewed as the atmosphere equivalent of shock melting, material motion, and far-field seismic response in the target. The three processes (early Processes, Atmospheric Processes, and Late Recovery Winds) are discussed at length.

Schultz, Pater H.↗
NASA NTRS1996

Aerodynamic shape optimization using control theory

Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Reuther, James↗
DOE OSTI2023

A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach

Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.

42 ENGINEERING↗
DOE OSTI2022

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗
NASA NTRS2001

Flame Synthesis of Single- and Multi-Walled Carbon Nanotubes and Nanofibers

Metal-catalyzed carbon nanotubes are highly sought for a diverse range of applications that include nanoelectronics, battery electrode material, catalysis, hydrogen storage media and reinforcing agents in polymer composites. These latter applications will require vast quantities of nanotubes at competitive prices to be economically feasible. Moreover, reinforcing applications may not require ultrahigh purity nanotubes. Indeed, functionalization of nanotubes to facilitate interfacial bonding within composites will naturally introduce defects into the tube walls, lessening their tensile strength. Current methods of aerosol synthesis of carbon nanotubes include laser ablation of composite targets of carbon and catalyst metal within high temperature furnaces and decomposition of a organometallics in hydrocarbons mixtures within a tube furnace. Common to each approach is the generation of particles in the presence of the reactive hydrocarbon species at elevated temperatures. In the laser-ablation approach, the situation is even more dynamic in that particles and nanotubes are borne during the transient cooling phase of the laser-induced plasma for which the temperature far exceeds that of the surrounding hot gases within the furnace process tube. A shared limitation is that more efficient methods of nanoparticle synthesis are not readily incorporated into these approaches. In contrast, combustion can quite naturally create nanomaterials such as carbon black. Flame synthesis is well known for its commercial scalability and energy efficiency. However, flames do present a complex chemical environment with steep gradients in temperature and species concentrations. Moreover, reaction times are limited within buoyant driven flows to tens of milliseconds. Therein microgravity can greatly lessen temperature and spatial gradients while allowing independent control of flame residence times. In preparation for defining the microgravity experiments, the work presented here focuses on the effect of catalyst particle size and reactant gas in 1g.

VanderWal, R. L.↗
NASA NTRS1994

Ignition dynamics of a laminar diffusion flame in the field of a vortex embedded in a shear flow

The role of streamwise-spanwise vorticity interactions that occur in turbulent shear flows on flame/vortex interactions is examined by means of asymptotic analysis and numerical simulation in the limit of small Mach number. An idealized model is employed to describe the interaction process. The model consists of a one-step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes which are then allowed to mix and react in the presence of a streamwise vortex embedded in a shear flow. It is found that the interaction of the streamwise vortex with shear gives rise to small-scale velocity oscillations which increase in magnitude with shear strength. These oscillations give rise to regions of strong temperature gradients via viscous heating, which can lead to multiple ignition points and substantially decrease ignition times. The evolution in time of the temperature and mass-fraction fields is followed, and emphasis is placed on the ignition time and structure as a function of vortex and shear strength.

Macaraeg, Michele G.↗