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 127 records · Page 7

An improved stochastic weighted particle method for boundary driven flows

Here, the stochastic weighted particle method (SWPM) is a generalization of the Direct Simulation Monte Carlo (DSMC) method where particle weights are variable and dynamic. SWPM is backed by a strong theoretical foundation but has not been critically evaluated for problems of practical interest. A thorough assessment of SWPM for boundary-driven flows reveals significant numerical artifacts near the boundary, notably a diverging heat flux. To correct the boundary heat flux, two modifications to SWPM are proposed: separated grouping and a spatially-dependent weight transfer function. To gauge the relative efficiency of SWPM in comparison to DSMC, a high-Mach-number wheel flow which forms a strong density gradient is also simulated.

97 MATHEMATICS AND COMPUTING↗

Predicting Turbulent Convective Heat Transfer in Fully Developed Duct Flows

The performance of an explicit algebraic stress model (EASM) is assessed in predicting the turbulent flow and forced heat transfer in both straight and wavy ducts, with rectangular, trapezoidal and triangular cross-sections, under fully developed conditions. A comparison of secondary flow patterns. including velocity vectors and velocity and temperature contours, are shown in order to study the effect of waviness on flow dynamics, and comparisons between the hydraulic parameters. Fanning friction factor and Nusselt number, are also presented. In all cases. isothermal conditions are imposed on the duct walls, and the turbulent heat fluxes are modeled using gradient-diffusion type models. The formulation is valid for Reynolds numbers up to 10(exp 5) and this minimizes the need for wall functions that have been used with mixed success in previous studies of complex duct flows. In addition, the present formulation imposes minimal demand on the number of grid points without any convergence or stability problems. Criteria in terms of heat transfer and friction factor needed to choose the optimal wavy duct cross-section for industrial applications among the ones considered are discussed.

Rokni, Masoud↗

A similarity rule for compressibility and sidewall-boundary-layer effects in two-dimensional wind tunnels

The effect of the sidewall boundary layer on flow in two-dimensional wind tunnels is determined. The small-disturbance and isentropic approximations are made, and the sidewall-boundary-layer dynamics are modeled with the von Karman momentum-integral equation. The effects of the edge-velocity-gradient term in the sidewall momentum integral, which is usually dominant near the model, and the compressibility term are shown to be similar. It is shown that the effect of sidewall suction around the model is not similar to two-dimensional flow. Comparisons with experiment are made to verify the similarity rule.

Barnwell, R. W.↗

Comparison of Various Supersonic Turbine Tip Designs to Minimize Aerodynamic Loss and Tip Heating

The rotor tips of axial turbines experience high heat flux and are the cause of aerodynamic losses due to tip clearance flows, and in the case of supersonic tips, shocks. As stage loadings increase, the flow in the tip gap approaches and exceeds sonic conditions. This introduces effects such as shock-boundary layer interactions and choked flow that are not observed for subsonic tip flows that have been studied extensively in literature. This work simulates the tip clearance flow for a flat tip, a diverging tip gap and several contoured tips to assess the possibility of minimizing tip heat flux while maintaining a constant massflow from the pressure side to the suction side of the rotor, through the tip clearance. The Computational Fluid Dynamics (CFD) code GlennHT was used for the simulations. Due to the strong favorable pressure gradients the simulations assumed laminar conditions in the tip gap. The nominal tip gap width to height ratio for this study is 6.0. The Reynolds number of the flow is 2.4 x 10(exp 5) based on nominal tip width and exit velocity. A wavy wall design was found to reduce heat flux by 5 percent but suffered from an additional 6 percent in aerodynamic loss coefficient. Conventional tip recesses are found to perform far worse than a flat tip due to severe shock heating. Overall, the baseline flat tip was the second best performer. A diverging converging tip gap with a hole was found to be the best choice. Average tip heat flux was reduced by 37 percent and aerodynamic losses were cut by over 6 percent.

Shyam, Vikram↗

Transport coefficient approach for characterizing nonequilibrium dynamics in soft matter

Nonequilibrium states in soft condensed matter require a systematic approach to characterize and model materials, enhancing predictability and applications. Among the tools, X-ray photon correlation spectroscopy (XPCS) provides exceptional temporal and spatial resolution to extract dynamic insight into the properties of the material. However, existing models might overlook intricate details. We introduce an approach for extracting the transport coefficient, denoted as $J(t)$, from the XPCS studies. This coefficient is a fundamental parameter in nonequilibrium statistical mechanics and is crucial for characterizing transport processes within a system. Our method unifies the Green–Kubo formulas associated with various transport coefficients, including gradient flows, particle–particle interactions, friction matrices, and continuous noise. We achieve this by integrating the collective influence of random and systematic forces acting on the particles within the framework of a Markov chain. We initially validated this method using molecular dynamics simulations of a system subjected to changes in temperatures over time. Subsequently, we conducted further verification using experimental systems reported in the literature and known for their complex nonequilibrium characteristics. The results, including the derived $J(t)$ and other relevant physical parameters, align with the previous observations and reveal detailed dynamical information in nonequilibrium states. This approach represents an advancement in XPCS analysis, addressing the growing demand to extract intricate nonequilibrium dynamics. Further, the methods presented are agnostic to the nature of the material system and can be potentially expanded to hard condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗

Quasi-steady solar wind dynamics

Progress in understanding the large scale dynamics of quasisteady, corotating solar wind structure was reviewed. The nature of the solar wind at large heliocentric distances preliminary calculations from a 2-D MHD model are used to demonstrate theoretical expectations of corotating structure out to 30 AU. It is found that the forward and reverse shocks from adjacent CIR's begin to interact at about 10 AU, producing new shock pairs flanking secondary CIR's. These sawtooth secondary CIR's interact again at about 20 AU and survive as visible entities to 30 AU. The model predicts the velocity jumps at the leading edge of the secondary CIR's at 30 AU should be very small but there should still be sizable variations in the thermodynamic and magnetic parameters. The driving dynamic mechanism in the distant solar wind is the relaxation of pressure gradients. The second topic is the influence of weak, nonimpulsive time dependence in quasisteady dynamics. It is suggested that modest large scale variations in the coronal flow speed on periods of several hours to a day may be responsible for many of the remaining discrepancies between theory and observation. Effects offer a ready explanation for the apparent rounding of stream fronts between 0.3 and 1.0 AU discovered by Helios.

Pizzo, V. J.↗

Numerical computation of space shuttle orbiter flow field

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous chemically reacting flow fields around 3-D bodies. The flow medium is assumed to be a multicomponent mixture of thermally perfect but calorically imperfect gases. The new PNS code solves the gas dynamic and species conservation equations in a coupled manner using a noniterative, implicit, approximately factored, finite difference algorithm. The space-marching method is made well-posed by special treatment of the streamwise pressure gradient term. The code has been used to compute hypersonic laminar flow of chemically reacting air over cones at angle of attack. The results of the computations are compared with the results of reacting boundary-layer computations and show excellent agreement.

Tannehill, John C.↗

Computational Study of Axisymmetric Off-Design Nozzle Flows

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles operating off-design at transonic Mach numbers have been completed. These computations span the very difficult transonic flight regime with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined, including the Explicit Algebraic Stress model. Computations have been completed at freestream Mach numbers of 0.9 and 1.2, and nozzle pressure ratios (NPR) of 4 and 6. Calculations completed with variable time-stepping (steady-state) did not converge to a true steady-state solution. Calculations obtained using constant timestepping (timeaccurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was the result of using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show reasonable agreement with experimental data. The SST turbulence model demonstrates the best overall agreement with experimental data.

DalBello, Teryn↗

Damping inter-area oscillation using reinforcement learning controlled TCSC

Appropriate damping control of the inter-area oscillation would ensure the grid stability and maintain the tie-line power transfer capability. In this paper, we proposed a novel reinforcement learning (RL) based power oscillation damping (POD) controller that uses Thyristor Controlled Series Compensators (TCSC) to damp inter-area oscillations. By leveraging the unbiased gradient direction estimation of the natural evolution strategy (NES), we successfully regulated the power flows on the tie-lines and damped inter-area oscillations through dynamically modulating the inserted reactance of the TCSC. The proposed RL-based POD controller has been tested on both two-area four-machine system and North American Minni-WECC system. Extensive studies have demonstrated the excellent performance of the proposed RL-based TCSC POD controller in damping inter-area oscillations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Corrosion Testing Trade Study on Testing Methods in Molten Salts

In assessing molten salt compatibility with structural materials, static testing is often reported because it is easy to conduct. However, in order to assess mass transfer, where material dissolves in the salt in the hot section and may precipitate out of the liquid in colder regions, a flowing salt experiment in a temperature gradient is needed. While a forced convection or pumped loop produces prototypic flowing conditions, it also is a multi-million dollar investment. Among the less expensive dynamic testing options, the thermal convection loop (TCL) was favored during the 1960’s development of molten salt reactor technology. Because the workforce and core facilities are currently available to conduct TCL experiments, it is the least expensive strategy to evaluate candidate structural alloys in flowing molten salt conditions.

36 MATERIALS SCIENCE↗

Numerical Simulation of Low-Density Shock-Wave Interactions

Computational Fluid Dynamics (CFD) numerical simulations of low-density shock-wave interactions for an incident shock impinging on a cylinder have been performed. Flow-field density gradient and surface pressure and heating define the type of interference pattern and corresponding perturbations. The maximum pressure and heat transfer level and location for various interaction types (i.e., shock-wave incidence with respect to the cylinder) are presented. A time-accurate solution of the Type IV interference is employed to demonstrate the establishment and the steadiness of the low-density flow interaction.

Glass, Christopher E.↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validation of Cryogenic Propellant Tank Self-Pressurization by Leveraging Reduced Order Modeling within Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank self-pressurization was performed using a hybrid Computational Fluid Dynamics (CFD) and reduced order modeling methodology. Data from a liquid hydrogen ground test conducted at the K-site facility at the National Aeronautics and Space Administration (NASA) Glenn Research Center was used for the validation effort. Liquid phase dynamics were explicitly resolved with a CFD tool. Vapor phase dynamics were modeled as a point mass that communicated heat from the tank wall to the liquid phase via a boundary condition used at the gas-liquid interface. The method proved to be more accurate, robust, and efficient than explicit resolution of the dynamics using a standard Volume of Fluid (VOF) methodology. The subject pressurization process was found to be heavily dependent upon both the relatively high liquid temperature gradient near the gas-liquid interface and the natural convection flow path. Modeling the gas-liquid interface as an immovable surface eliminated temperature gradient destroying gas-liquid interface velocities observed in VOF simulations, and correspondingly enabled more rapid simulation since interface advection was not allowed. The single phase computational domain also facilitated the ability to demonstrate spatial resolution convergence of natural convection cells within the liquid which significantly impacted the tank pressurization rate. This work was used to demonstrate the critical physics for tank self-pressurization and numerical methodologies that may be used to best resolve those physics. The findings informed development and operation of production level CFD tools used in the Fluid Dynamics Branch at NASA Marshall Space Flight Center.

J. M. Brodnick↗

Large Eddy Simulation Study for Fluid Disintegration and Mixing

A new modeling approach is based on the concept of large eddy simulation (LES) within which the large scales are computed and the small scales are modeled. The new approach is expected to retain the fidelity of the physics while also being computationally efficient. Typically, only models for the small-scale fluxes of momentum, species, and enthalpy are used to reintroduce in the simulation the physics lost because the computation only resolves the large scales. These models are called subgrid (SGS) models because they operate at a scale smaller than the LES grid. In a previous study of thermodynamically supercritical fluid disintegration and mixing, additional small-scale terms, one in the momentum and one in the energy conservation equations, were identified as requiring modeling. These additional terms were due to the tight coupling between dynamics and real-gas thermodynamics. It was inferred that if these terms would not be modeled, the high density-gradient magnitude regions, experimentally identified as a characteristic feature of these flows, would not be accurately predicted without the additional term in the momentum equation; these high density-gradient magnitude regions were experimentally shown to redistribute turbulence in the flow. And it was also inferred that without the additional term in the energy equation, the heat flux magnitude could not be accurately predicted; the heat flux to the wall of combustion devices is a crucial quantity that determined necessary wall material properties. The present work involves situations where only the term in the momentum equation is important. Without this additional term in the momentum equation, neither the SGS-flux constant-coefficient Smagorinsky model nor the SGS-flux constant-coefficient Gradient model could reproduce in LES the pressure field or the high density-gradient magnitude regions; the SGS-flux constant- coefficient Scale-Similarity model was the most successful in this endeavor although not totally satisfactory. With a model for the additional term in the momentum equation, the predictions of the constant-coefficient Smagorinsky and constant-coefficient Scale-Similarity models were improved to a certain extent; however, most of the improvement was obtained for the Gradient model. The previously derived model and a newly developed model for the additional term in the momentum equation were both tested, with the new model proving even more successful than the previous model at reproducing the high density-gradient magnitude regions. Several dynamic SGS-flux models, in which the SGS-flux model coefficient is computed as part of the simulation, were tested in conjunction with the new model for this additional term in the momentum equation. The most successful dynamic model was a "mixed" model combining the Smagorinsky and Gradient models. This work is directly applicable to simulations of gas turbine engines (aeronautics) and rocket engines (astronautics).

Bellan, Josette↗

Interaction between step fuel injectors on opposite walls in a supersonic combustor model

Results are presented from an experimental investigation of perpendicular, hydrogen fuel injection and combustion from opposing walls in a scramjet combustor model using a longitudinally staged laterally inline step-injection configuration. The model represents a portion of the flow in the Langley integrated modular scramjet engine combustor operating at a flight Mach number of 7. When operating at a ratio of jet pressure to free-stream dynamic pressure of 3, the injectors produce a bulk equivalence ratio of unity. This investigation represents part of a continuing study of the modular engine fuel injectors and is specifically designed to eliminate the adverse lateral pressure gradient observed at the injector location in a previous test. Flow survey contours at three axial locations, ranging from one-third of the engine combustor length to the total engine combustor length, are used to determine mixing efficiency and fuel distribution. Wall static pressures are analyzed by using one-dimensional theory to determine the combustion efficiency. Results show a significant improvement over previous injector designs tested in this duct geometry.

Mcclinton, C. R.↗