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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2025

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump↗

DOE OSTI2021

Randomized Sketching Algorithms for Low-Memory Dynamic Optimization

This paper develops a novel limited-memory method to solve dynamic optimization problems. The memory requirements for such problems often present a major obstacle, particularly for problems with PDE constraints such as optimal flow control, full waveform inversion, and optical tomography. In these problems, PDE constraints uniquely determine the state of a physical system for a given control; the goal is to find the value of the control that minimizes an objective. While the control is often low dimensional, the state is typically more expensive to store. This paper suggests using randomized matrix approximation to compress the state as it is generated and shows how to use the compressed state to reliably solve the original dynamic optimization problem. Concretely, the compressed state is used to compute approximate gradients and to apply the Hessian to vectors. The approximation error in these quantities is controlled by the target rank of the sketch. This approximate first- and second-order information can readily be used in any optimization algorithm. As an example, we develop a sketched trust-region method that adaptively chooses the target rank using a posteriori error information and provably converges to a stationary point of the original problem. Numerical experiments with the sketched trust-region method show promising performance on challenging problems such as the optimal control of an advection-reaction-diffusion equation and the optimal control of fluid flow past a cylinder.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Planar turbulent wakes under pressure gradient: Integral and self-similarity analyses

By using a combination of integral and self-similarity analyses, the generalized analytical solutions for the mean transverse velocity and Reynolds shear stress are rigorously derived for the first time for the far field of planar turbulent wakes under arbitrary pressure gradients. Specifically, by assuming self-similarity for the mean axial velocity, the analytical formulation for the mean transverse velocity is obtained from the integral of the mean continuity equation, and the analytical formulation for the Reynolds shear stress is obtained from the integral of the momentum equation. The generalized analytical formulations for the mean transverse velocity and Reynolds shear stress consist of multiple components, each with its unique scale and physical mechanism. In the zero pressure gradient limit, the generalized formulations recover the single-scale equations reported by Wei, Liu, and Livescu. Furthermore, simpler approximate formulations for the mean transverse velocity and Reynolds shear stress are also obtained, and show excellent agreement with the experimental measurements. Here, the findings provide new insights into the properties of planar turbulent wakes under pressure gradients, filling some long-standing gaps in the existing literature.

42 ENGINEERING↗

NASA NTRS1973

Turbulent shear flow of continuum plasma.

The turbulence structure of weakly ionized, continuum plasma is analyzed with emphasis on the turbulent fluctuation and the transport properties of the ionized species. The mean gradients of the velocity and ionized-species concentrations are considered to be uniform. It is assumed that the turbulence Reynolds number is small and that the turbulence field is homogeneous in order to close the dynamic equations of turbulence. It was found that the ratio of the Debye shielding length to the local integral scale of turbulence must be sufficiently small in order for the detailed quasi-neutrality to prevail. The condition of the detailed quasi-neutrality is necessary, in addition to the equal mean concentrations of electrons and ions, in order for the turbulence properties involving the electrons to be the same as those involving the ions.

Chung, P. M.↗

NASA NTRS1979

Lightman-Eardley instabilities and accretion disk thickening

After reviewing the role of Compton scattering in accretion disks around black holes, it is discussed whether Lightman-Eardley (LE) secular instabilities can trigger and maintain Pringle-Rees (PR) thermal instabilities. The radiative-transfer-equation and equation-of-state criteria for LE stability in alpha-viscosity-law disk models and dynamic viscosity criteria for more general situations is derived. On the basis of these considerations the LE instability is insufficient for inducing PR instabilities and hot thick inner regions important in accretion-disk models of compact hard X-ray sources. The density thinning due to radial velocity gradients in the accretion flow is suggested as a more likely and satisfactory mechanism.

Stoeger, W. R.↗

NASA NTRS1974

Magnetic field and electric current structure in the chromosphere

The three-dimensional vector magnetic field structure in the chromosphere above an active region is deduced by using high-resolution H-alpha filtergrams together with a simultaneous digital magnetogram. An analog model of the field is made with 400 metal wires representing field lines that outline the H-alpha structure. The height extent of the field is determined from vertical field-gradient observations around sunspots, from observed fibril heights, and from an assumption that the sources of the field are largely local. The computed electric currents (typically 10 mA/sq m) are found to flow in patterns not similar to observed features and not parallel to magnetic fields. Force structures correspond to observed solar features; the dynamics to be expected include: downward motion in bipolar areas in the lower chromosphere, an outflow of the outer chromosphere into the corona with radially outward flow above bipolar plage regions, and motion of arch filament systems.

Dravins, D.↗

NASA NTRS1990

Modeling the pressure-strain correlation of turbulence: An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗

NASA NTRS1991

Modelling the pressure-strain correlation of turbulence - An invariant dynamical systems approach

Speziale, Charles G.↗

DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1994

Large-eddy simulation of flow through a plane, asymmetric diffuser

A challenge for traditional turbulence modeling, based on the Reynolds averaged Navier-Stokes equations, remains the accurate prediction of 'mild', adverse pressure-gradient driven separation from a smooth surface. With this study we want to explore the capability of large-eddy simulation to predict the separation which occurs on the deflected wall of an asymmetric, plane diffuser with opening angle of 10 deg. The flow through the plane diffuser exhibits some additional interesting physical phenomena which make it a challenging test case. In addition to 'mild' separation about halfway down the deflected ramp, the flow is characterized by a small backflow zone with stalled fluid in the rear part of the expanding section. The turbulent flow entering the diffuser is subject to combined adverse and radial pressure gradients stemming from the convex curvature. Finally the flow recovers into a developed, turbulent channel flow in the outlet section. Obi et al. provide measurements of mean flow, Reynolds stresses, and pressure recovery, which were obtained by means of LDV in a wind tunnel. The objective of this study is to investigate whether LES with the standard dynamic model is able to accurately predict the flow in the one-sided diffuser and to explore the resolution requirements and associated costs.

Kaltenbach, Hans-Jakob↗

NASA NTRS2003

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗

NASA NTRS1990

A three-dimensional dynamic solution-adaptive mesh algorithm

A solution-adaptive grid algorithm has been developed for use in two and three dimensions. The algorithm uses a transformation from the cartesian coordinate system to a general coordinate space, which will be defined as a parallelepiped. A weighting function for adaption of the grid is developed that will allow adaption to the gradients of any combination of dependent variables in the flow. The adaption is carried out in the parametric space and a simple inverse mapping to return the new parametric space to the physical space is derived. The concept used to relocate the grid-points in the parametric space is based on the center of mass of distributed weights. Solution-adaptive results are presented for various laminar flows in two dimensions and for mathematical weighting functions in three dimensions.

Benson, Rusty A.↗

DOE OSTI2024

Osmotic and phoretic competition explains chemotaxic assembly and sorting

Microscale objects responding to chemical gradients by migrating toward or away from a preferred species is a simple yet constitutive mechanism by which transport occurs in biological organisms. Synthetic chemotaxis provides key physical descriptions of simplified systems that can be used in biological models, or in the creation of advanced responsive material systems. In this article, we provide a quantitative framework for understanding synthetic chemotaxis of microparticles which involves a competition between phoresis and osmosis. We present separate quantitative measurements of phoresis and osmosis acting on individual taxing particles, finding that phoresis follows the long-predicted v ∼ 1 / r 2 scaling while the osmotic contribution depends on the geometry and details of the system, and must be solved on a case-by-case basis. Through this, we are able to develop a more accurate picture of particle transport at the single particle level. Equipped with this approach, we go on to describe how high concentrations of particles in a symmetric chemical gradient grow close-packed hives that reach a steady-state size tunable through light intensity or particle size. Last, we demonstrate that mixed particles experiencing the same chemical gradient will selectively migrate toward or away depending on the nature of the particle surface, thereby locally sorting out a particular species. We anticipate these results will be important in describing both biological and synthetic chemotaxis in phoretic systems and should bring a wealth of studies that take advantage of competing osmotic flows to illicit unexpected dynamic active behavior.

Science & Technology - Other Topics↗

NASA NTRS1981

Energetic ion beams at the edge of the plasma sheet - ISEE 1 observations plus a simple explanatory model

It has been found that energetic ions with energies in the range from 10 keV to 2 MeV are an important part of the geomagnetic tail particle population. A description is presented of energetic ion, three-dimensional distribution observations from ISEE 1 during 0400-2000 UT, May 3, 1978, a time when several encounters with the geomagnetic tail plasma sheet occurred. Pronounced energetic ion streaming is observed at each entry to and exit from the plasma sheet. Azimuthal gradients observed during these times indicate that this streaming occurs within two gyroradii of the plasma sheet edge. Streaming characteristics show both narrow and broad angular flows at arbitrary angles to B. These flows are observed simultaneously and individually. The time, energy, and pitch angle evolution of the ion fluxes show that ISEE 1 is encountering a dynamic and not a static equilibrium situation at the edge of the plasma sheet.

Williams, D. J.↗

NASA NTRS1999

SAGE: The Self-Adaptive Grid Code

The multi-dimensional self-adaptive grid code, SAGE, is an important tool in the field of computational fluid dynamics (CFD). It provides an efficient method to improve the accuracy of flow solutions while simultaneously reducing computer processing time. Briefly, SAGE enhances an initial computational grid by redistributing the mesh points into more appropriate locations. The movement of these points is driven by an equal-error-distribution algorithm that utilizes the relationship between high flow gradients and excessive solution errors. The method also provides a balance between clustering points in the high gradient regions and maintaining the smoothness and continuity of the adapted grid, The latest version, Version 3, includes the ability to change the boundaries of a given grid to more efficiently enclose flow structures and provides alternative redistribution algorithms.

Davies, Carol B.↗

NASA NTRS2003

Marangoni Convection and Deviations from Maxwells' Evaporation Model

We investigate the convective dynamics of evaporating pools of volatile liquids using an ultra-sensitive thermal imaging camera. During evaporation, there are significant convective flows inside the liquid due to Marangoni forces. We find that Marangoni convection during evaporation can dramatically affect the evaporation rates of volatile liquids. A simple heat balance model connects the convective velocities and temperature gradients to the evaporation rates.

Segre, P. N.↗

DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗

NASA NTRS2011

An Intercomparison of Tropospheric Ozone Retrievals Derived from Two Aura Instruments and Measurements in Western North America in 2006

Two recently developed methods for quantifying tropospheric ozone abundances based on Aura data, the Trajectory‐enhanced Tropospheric Ozone Residual (TTOR) and an assimilation of Aura data into Goddard Earth Observing System Version 4 (ASM), are compared to ozone measurements from ozonesonde data collected in April-May 2006 during the INTEX Ozonesonde Network Study 2006 (IONS‐06) campaign. Both techniques use Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) observations. Statistics on column ozone amounts for both products are presented. In general, the assimilation compares better to sonde integrated ozone to 200 hPa (28.6% difference for TTOR versus 2.7% difference for ASM), and both products are biased low. To better characterize the performance of ASM, ozone profiles based on the assimilation are compared to those from ozonesondes. We noted slight negative biases in the lower troposphere, and slight positive biases in the upper troposphere/lower stratosphere (UT/ LS), where we observed the greatest variability. Case studies were used to further understand ASM performance. We examine one case from 17 April 2006 at Bratt's Lake, Saskatchewan, where geopotential height gradients appear to be related to an underestimation in the ASM in the UT/LS region. A second case, from 21 April 2006 at Trinidad Head, California, is a situation where the overprediction of ozone in the UT/LS region does not appear to be due to current dynamic conditions but seems to be related to uncertainty in the flow pattern and large differences in MLS observations upstream.

UTLS↗