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 433 records · Page 24

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. Finally, the patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

36 MATERIALS SCIENCE↗

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗

A consistent and conservative Phase-Field method for multiphase incompressible flows

In the present study, a consistent and conservative Phase-Field method, including both the model and scheme, is developed for multiphase flows with an arbitrary number of immiscible and incompressible fluid phases. The consistency of mass conservation and the consistency of mass and momentum transport are implemented to address the issue of physically coupling the Phase-Field equation, which locates different phases, to the hydrodynamics. These two consistency conditions, as illustrated, provide the “optimal” coupling because (i) the new momentum equation resulting from them is Galilean invariant and implies the kinetic energy conservation, regardless of the details of the Phase-Field equation, and (ii) failures of satisfying the second law of thermodynamics or the consistency of reduction of the multiphase flow model only result from the same failures of the Phase-Field equation but are not due to the new momentum equation. Physical interpretation of the consistency conditions and their formulations are first provided, and general formulations that are obtained from the consistency conditions and independent of the interpretation of the velocity are summarized. Then, the present consistent and conservative multiphase flow model is completed by selecting a reduction consistent Phase-Field equation. Several novel techniques are developed to inherit the physical properties of the multiphase flows after discretization, including the gradient-based phase selection procedure, the momentum conservative method for the surface force, and the general theorems to preserve the consistency conditions on the discrete level. Equipped with those novel techniques, a consistent and conservative scheme for the present multiphase flow model is developed and analyzed. The scheme satisfies the consistency conditions, conserves the mass and momentum, and assures the summation of the volume fractions to be unity, on the fully discrete level and for an arbitrary number of phases. All those properties are numerically validated. Finally, numerical applications demonstrate that the present model and scheme are robust and effective in studying complicated multiphase dynamics, especially for those with large-density ratios.

97 MATHEMATICS AND COMPUTING↗

Fluid Dynamics and Thermodynamics of Vapor Phase Crystal Growth

The ground-based research effort under this program is concerned with systematic studies of the effects of variations: (1) of the relative importance of buoyancy-driven convection, and (2) of diffusion and viscosity conditions on crystal properties. These experimental studies are supported by thermodynamic characterizations of the systems, based on which fluid dynamic parameters can be determined. The specific materials under investigation include: the GeSe-GeI4, Ge-GeI4, HgTe-HgI2, and Hg sub (1-x)Cd sub (x) Te-HgI2 systems. Mass transport rate studies of the GeSe-GeI system as a function of orientation of the density gradient relative to the gravity vector demonstrated the validity of flux anomalies observed in earlier space experiments. The investigation of the effects of inert gases on mass flux yielded the first experimental evidence for the existence of a boundary layer in closed ampoules. Combined with a thorough thermodynamic analysis, a transport model for diffusive flow including chemical vapor transport, sublimation, and Stefan flow was developed.

Wiedemeier, H.↗

Application of a self-adaptive grid method to complex flows

A directional-split, modular, user-friendly grid point distribution code is applied to several test problems. The code is self-adaptive in the sense that grid point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and constants relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flow-field gradients and/or geometric features, are used to determine the point distribution. Points are redistributed along grid lines in a specified direction in an elliptic manner over a user-specified subdomain, while orthogonality and smoothness are controlled in a parabolic (marching) manner in the remaining directions. Multidirectional adaption is achieved by sequential application of the method in each coordinate direction. The flow-field solution is redistributed onto the newly distributed grid points after each unidirectional adaption by a simple one-dimensional interpolation scheme. For time-accurate schemes such interpolation is not necessary and time-dependent metrics are carried in the fluid dynamic equations to account for grid movement.

Deiwert, G. S.↗

Development of new flux splitting schemes

Maximizing both accuracy and efficiency has been the primary objective in designing a numerical algorithm for computational fluid dynamics (CFD). This is especially important for solutions of complex three dimensional systems of Navier-Stokes equations which often include turbulence modeling and chemistry effects. Recently, upwind schemes have been well received for their capability in resolving discontinuities. With this in mind, presented are two new flux splitting techniques for upwind differencing. The first method is based on High-Order Polynomial Expansions (HOPE) of the mass flux vector. The second new flux splitting is based on the Advection Upwind Splitting Method (AUSM). The calculation of the hypersonic conical flow demonstrates the accuracy of the splitting in resolving the flow in the presence of strong gradients. A second series of tests involving the two dimensional inviscid flow over a NACA 0012 airfoil demonstrates the ability of the AUSM to resolve the shock discontinuity at transonic speed. A third case calculates a series of supersonic flows over a circular cylinder. Finally, the fourth case deals with tests of a two dimensional shock wave/boundary layer interaction.

Liou, Meng-Sing↗

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurements of flow phenomena induced by suction through perforated and partially plugged surfaces

Efforts were directed towards completing construction of the windtunnel test section, assembling instrumentation, programming the data acquisition and reduction system, adjusting the streamwise pressure gradient of the test section, calibrating the hot-wire anemometer probe, and constructing and testing a smoke generator. The test section was installed in the wind tunnel and is completely operational. The streamwise pressure gradient was adjusted to be nominally zero at a free-stream velocity of 3.05 m/s (10 ft/s). This was accomplished by adjusting the upper wall of the test section to be slightly divergent. The change in static pressure between any two streamwise locations in the test section was less than one percent of the free-stream dynamic pressure. A suitable means was found for accurately calibrating the hot-wire probe which is used to measure boundary-layer velocity profiles and fluctuating velocities.

Wilkinson, S. P.↗

Raman Temperature and Density Measurements in Supersonic Jets

Prediction of flow-field properties in supersonic jets using computational fluid dynamics (CFD) code predictions has become routine; however, obtaining accurate solutions becomes more challenging when there is a significant temperature difference between the jet core and the ambient air and/or compressibility effects are significant. Benchmark sets of flow field property data are required in order to assess current CFD capabilities and develop better modeling approaches for these turbulent flow fields where accurate calculation of temperatures and turbulent heat flux is important. Particle Image Velocimetry (PIV), spontaneous rotational Raman scattering (SRS) spectroscopy and Background Oriented Schlieren (BOS) have been previously used to acquire measurements of the mean and root mean square (rms) velocities, the mean and rms gas temperatures and density gradients in subsonic jet flows and film cooling flows. In this work, the ability to measure density is added to the list of measurands available using the acquired Raman spectra. The suite of measurement techniques are now applied to supersonic jet flows. The computation of the local gas pressure in the potential core of an over-expanded jet is demonstrated using the Raman measured gas temperature and density. Additionally, a unique density feature in temperature matched, perfectly expanded jet flow shear layers identified using BOS was verified using the Raman measurement technique. These non-intrusive flow measurements are compared against RANS predictions of the supersonic jet flow properties as a means of assessing their prediction accuracy.

Raman Spectroscopy↗

Ink-substrate interactions during 3D printing revealed by time-resolved coherent X-ray scattering

Additive printing techniques are regarded as revolutionary and versatile methods of advanced device manufacturing, stemming from the possibility to pattern materials on a custom-based approach and the potential to create novel microstructures and achieve new functionalities. Despite these advantages, the inherent anisotropy of the printing process is a source of property gradients within the printed materials, often associated with variable and/or poor performance. Up to date, the evolutionary pathways associated with printing have largely remained unaddressed, mainly owing to the difficulty to study the transformations induced in the material during processing. Time-resolved coherent X-ray scattering techniques, such as X-ray photon correlation spectroscopy, enable the in situ study of transient nanoscale and mesoscale states in a large variety of materials, including amorphous ones, by directly accessing the most relevant timescales and length scales of their nanoscale and mesoscale dynamics, self-assembly, and mesostructure evolution. In this work, we conduct in operando studies of continuous-flow direct writing with colloidal inks, focusing on how the ink formulation and ink-substrate interactions affect the processes that determine the macroscopic properties of the printed materials. We find fundamental differences in the ink structural relaxations emerging from the primary colloid properties (monodisperse versus aggregated colloids) and the substrates surface energy and mechanical properties. Our work helps to reveal and quantify the basic science governing the evolution of 3D-printed materials during processing, ultimately improving engineering criteria for the design of printable materials.

36 MATERIALS SCIENCE↗

Models of thermal/chemical boundary layer convection: Potential application to Venus

The upper boundary layer of Venus is comprised of at least two distinct chemical components, mantle and crust. Fluid dynamical models of convection within Venus' mantle were primarily of the thermal boundary layer type. Models assessing the ability of convective mantle flows to deform the crust were undertaken, but models exploring the effects of a variable thickness crust on mantle convection were largely lacking. A Venusian crust of variable thickness could couple back into, and alter, the mantle flow patterns that helped create it, leading to deformation mechanisms not predicted by purely thermal boundary layer convection models. This possibility is explored through a finite element model of thermal/chemical boundary layer convection. Model results suggest that a crust of variable thickness can serve as a mantle flow driver by perturbing lateral temperature gradients in the upper mantle. Resulting mantle flow is driven by the combination of free convective and nonuniform crustal distribution. This combination can lead to a flow instability manifest in the occurrence of episodic mantle lithosphere subduction initiated at the periphery of a crustal plateau. The ability of a light, near surface, chemical layer to potentially alter mantle flow patterns suggest that mantle convection and the creation and/or deformation of such a chemical layer may be highly nonseparable problems on time scales of 10(exp 8) years.

Lenardic, A.↗

Lidar-Based Evaluation of HRRR Performance in California’s Diablo Range

The performance of the NOAA High-Resolution Rapid Refresh (HRRR) model for capturing low-level winds near a wind energy production site during summer 2019 is evaluated. This study catalogs the ability of HRRR to predict boundary layer dynamics relevant to wind energy interests over complex terrain, which has presented challenges for weather and energy forecasting. Performance is evaluated by comparing HRRR output to wind-profiling Doppler lidars at Lawrence Livermore National Laboratory Site 300. HRRR captured the diurnal profile of horizontal winds in the observed 150-m layer, despite strong underpredictions (∼4 m s −1 ) during evening and nighttime hours. These underpredictions may be a result of local speedup flows observed by the lidars, which were unresolved in HRRR due to their small spatial extent. HRRR bias magnitude relative to observations was found to be minimal during days with synoptic-scale troughs and strong 850-hPa geopotential gradients, while bias magnitude was maximal during days with synoptic ridging and weak 850-hPa geopotential gradients. To translate wind speed predictions to energy forecasting, generic turbine models were used to estimate power generation for turbines characteristic of the nearby Altamont Pass Wind Resource Area. Results show that HRRR-based energy estimates predicted daytime power generation adequately relative to lidar-based estimates with an 18-h lead time (bias magnitude < 0.4 MW from 0900 to 1400 LT) but overpredicted power during the rest of the diurnal cycle (bias > 1 MW). These results demonstrate conditions under which HRRR performs well for wind energy applications in complex terrain, while highlighting biases that require further investigation to support usage of a high-resolution model for wind energy forecasts.

Boundary layer↗

The Evolution of Finite Amplitude Wavetrains in Plane Channel Flow

We consider a viscous incompressible fluid flow driven between two parallel plates by a constant pressure gradient. The flow is at a finite Reynolds number, with an 0(l) disturbance in the form of a traveling wave. A phase equation approach is used to discuss the evolution of slowly varying fully nonlinear two dimensional wavetrains. We consider uniform wavetrains in detail, showing that the development of a wavenumber perturbation is governed by Burgers equation in most cases. The wavenumber perturbation theory, constructed using the phase equation approach for a uniform wavetrain, is shown to be distinct from an amplitude perturbation expansion about the periodic flow. In fact we show that the amplitude equation contains only linear terms and is simply the heat equation. We review, briefly, the well known dynamics of Burgers equation, which imply that both shock structures and finite time singularities of the wavenumber perturbation can occur with respect to the slow scales. Numerical computations have been performed to identify areas of the (wavenumber, Reynolds number, energy) neutral surface for which each of these possibilities can occur. We note that the evolution equations will breakdown under certain circumstances, in particular for a weakly nonlinear secondary flow. Finally we extend the theory to three dimensions and discuss the limit of a weak spanwise dependence for uniform wavetrains, showing that two functions are required to describe the evolution. These unknowns are a phase and a pressure function which satisfy a pair of linearly coupled partial differential equations. The results obtained from applying the same analysis to the fully three dimensional problem are included as an appendix.

Hewitt, R. E.↗

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗

Accretion onto magnetized neutron stars - Structure and interchange instability of a model magnetosphere

A self-consistent model is analyzed for the spherical infall of weakly magnetized plasma into the magnetosphere of a slowly rotating, strongly magnetized neutron star. It is shown that spherical infall is probably a good approximation for X-ray sources which accrete from a stellar wind. The location of the standoff shock which halts the hypersonic infall is estimated along with the emission from the shocked layer. The location of the equilibrium magnetopause and the structure of the magnetic field within it are calculated; it is found that the magnetic poles are true cusps and that the entry of gas due to equilibrium flow across a cusp is almost certainly dominated by the interchange instability near the magnetic equator. The energy principle is applied to derive necessary conditions for the occurrence of this instability. The results indicate that the strong magnetic-pressure gradient stabilizes the gas unless moderately strong radiative cooling takes place and that the cooled plasma enters the magnetosphere as long filaments capable of moving between field lines. The rate at which the equilibrium magnetopause can 'absorb' mass and momentum is derived, the validity of the approximations employed is discussed, and the likely evolution of the sinking filaments is outlined to show that the spatial distribution of the plasma is determined mainly by the dynamics and thermodynamics of the filaments rather than the magnetic-field structure.

Arons, J.↗

A lifting surface theory in rotational flow

The partial differential equation for small disturbance steady rotational flow in three dimensions is solved through an integral equation approach. The solution is obtained by using the method of weighted residuals. Specific applications are directed to wings in nonuniform subsonic parallel streams with velocity varying in vertical and spanwise directions and to airfoils in nonuniform freestream. Comparison with limited known results indicates that the present method is reasonably accurate. Numerical results for the lifting pressure of airfoil, lift, induced drag, and pitching moments of airfoil, lift, induced drag, and pitching moments of elliptic, rectangular, and delta wings in a jet, wake, or monotonic sheared stream are presented. It is shown that, in addition to the effect of local dynamic pressures, a positive velocity gradient tends to enhance the lift.

Shiau, M. J.↗

Low-Temperature Oxidation Reactions and Cool Flames at Earth and Reduced Gravity

Non-isothermal studies of cool flames and low temperature oxidation reactions in unstirred closed vessels are complicated by the perturbing effects of natural convection at earth gravity. Buoyant convection due to self-heating during the course of slow reaction produces spatio-temporal variations in the thermal and thus specie concentration fields due to the Arrhenius temperature dependence of the reaction rates. Such complexities have never been quantitatively modeled and were the primary impetus for the development of CSTR's (continuously stirred tank reactors) 30 years ago. While CSTR's have been widely adopted since they offer the advantage of spatial uniformity in temperature and concentration, all gradients are necessarily destroyed along with any structure that may otherwise develop. Microgravity offers a unique environment where buoyant convection can be effectively minimized and the need for stirring eliminated. Moreover, eliminating buoyancy and the need for stirring eliminates complications associated with the induced hydrodynamic field whose influence on heat transport and hot spot formation, hence explosion limits, is not fully realized. The objective of this research is to quantitatively determine and understand the fundamental mechanisms that control the onset and evolution of low temperature reactions and cool flames in both static and flow reactors. Microgravity experiments will be conducted to obtain benchmark data on the structure (spatio-temporal temperature, concentration, flow fields), the dynamics of the chemical fronts, and the ignition diagrams (pressure vs. temperature). Ground-based experiments will be conducted to ascertain the role of buoyancy. Numerical simulations including detailed kinetics will be conducted and compared to experiment.

Pearlman, Howard↗