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

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A Theory for the Response of Tropical Moist Convection to Mechanical Orographic Forcing

Abstract Spatial patterns of tropical rainfall are strongly influenced by mountains. Although theories for precipitation induced by convectively stable upslope ascent exist for the midlatitudes, these do not represent the interaction of moist convection with orographic forcing. Here, we present a theory for convective precipitation produced by the mechanical interaction of a tropical ridge with a basic-state horizontal wind. Deviations from this basic state are represented as the sum of a “dry” perturbation, due to the stationary orographic gravity wave, and a “moist” perturbation that carries the convective response. The moist component dynamics are subject to the weak temperature gradient approximation; they are forced by the dry mode’s influence on lower-tropospheric moisture and temperature. Analytical solutions provide estimates of the precipitation distribution, including peak precipitation, upstream extent, and rain shadow extent. The theory can be used with several degrees of complexity depending on the technique used to compute the dry mode, which can be drawn from linear mountain wave theory or full numerical simulations. To evaluate the theory, we use a set of convection-permitting simulations with a flow-perpendicular ridge in a long channel. The theory makes a good prediction for the cross-slope precipitation profile, indicating that the organization of convective rain by orography can be quantitatively understood by considering the effect of stationary orographic gravity waves on a lower-tropospheric convective quasi-equilibrium state.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Frontal Maintenance in Submesoscale Flows

Classic deformation theory includes parameters—divergence, total strain, and vorticity—that are invariant to changes in the coordinate system. However, these parameters are sometimes ambiguous with respect to characterizing how fronts are formed and maintained because the presence of a front imposes a reference coordinate system. To help remedy this ambiguity, we propose a framework in frontal coordinates based on along- and cross-front velocity gradients to better characterize frontal maintenance, which can also be used to define divergence and normal strain in frontal coordinates. The framework with these four parameters (along-, cross-front velocity gradients, divergence, and normal strain in frontal coordinate) defines eight characteristic flow types at a front, providing a complete characterization of the flow that strengthens or weakens a front. This framework highlights the importance of the “strain efficiency” concept, which unambiguously defines the contribution of total strain to frontogenesis. Two examples, one based on a realistic simulation of submesoscales in the northern Gulf of Mexico and the other based on an idealized model with similar flow characteristics, are provided to demonstrate how this framework can be used to enhance our understanding of frontal dynamics in submesoscale flows.

frontogenesis↗
DOE OSTI2023

The transport dynamics of tens of micrometer-sized water droplets in RF atmospheric pressure glow discharges

Abstract Charging of particles having diameters of tens of microns has been extensively studied at atmospheric pressure in the context of, for example, electrostatic precipitators where the focus was on unipolar charging. The ambipolar charging of particles in atmospheric pressure plasmas, and of droplets in particular, has received less attention. The plasma activation of droplets is of interest for water purification, fertilizer production and materials synthesis, all of which depend on the transport of the droplets through the plasma, which in turn depends on their charging. In this paper, we report on the transport dynamics of water droplets, tens of microns in diameter, carried by the gas flow through an atmospheric pressure radiofrequency glow discharge sustained in helium. The droplets pass through the plasma with minimal evaporation and without reaching the Rayleigh limit. The droplet trajectory in the presence and absence of the plasma provides insights on the forces acting on the droplet. The measurements were analyzed using results from a three-dimensional fluid model and a two-dimensional plasma hydrodynamics model. We found that the transport dynamics as the droplet enters and leaves the plasma are due to differential charging of the droplet in the plasma gradients of the bounding sheaths to the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Numerical investigations in three-dimensional internal flows

The flow in the transonic test facility was investigated using the three dimensional computational fluid dynamics techniques. The application of the full Navier-Stokes three dimensional code to the flow qualities in the contraction section of transonic wind tunnel is discussed. Initially, two dimensional solutions indicated the possibility for large secondary flow to exist as a result of the asymmetries involved in the contraction section as it is constructed. The results of a full three dimensional solution indicate that only minor pressure variations actually occur in the contraction section within any given cross flow plane. Further analysis of the three dimensional solution indicated that these slight lateral pressure gradients lead to negligible secondary flows, except within a small region in the corners within the boundary layer. On the basis of present solution, it would not be expected that any flow asymmetries and/or secondary flow present within contraction section are associated with the methods by which the contraction is implemented in its present configuration.

Rose, William C.↗
NASA NTRS1992

Dynamics Explorer measurements of particles, fields, and plasma drifts over a horse-collar auroral pattern

As shown from ground-based measurements and satellite-borne imagers, one type of global auroral pattern characteristic of quiet (usually northward IMF) intervals is that of a contracted but thickened emission region of a pattern referred to as 'horse-collar' aurora (Hones et al., 1989). In this report we use the Dynamics Explorer data set to examine a case in which this horse-collar pattern was observed by the DE-1 auroral imager, while at the same time DE-2, at lower altitude, measured precipitating particles, electric and magnetic fields, and plasma drifts. Our analysis shows that, in general, there is close agreement between the optical signatures and the particle precipitation patterns. In many instances, over scales ranging from tens to a few hundred kilometers, electron precipitation features and upward field-aligned currents are observed at locations where the plasma flow gradients indicate negative V-average x E. The particle, plasma, and field measurements made along the satellite track and the 2D perspective of the imager provide a means of determining the configuration of convective flows in the high-latitude ionosphere during this interval of northward IMF. Recent mapping studies are used to relate the low-altitude observations to possible magnetospheric source regions.

Sharber, J. R.↗
NASA NTRS2009

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗
NASA NTRS1994

Thermocapillary flow with evaporation and condensation and its effect on liquid retention in low-G fluid acquisition devices

The steady motion, thermal and free surface behavior of a volatile, wetting liquid in microgravity are studied using scaling and numerical techniques. The objective is to determine whether the thermocapillary and two-phase convection arising from thermodynamic nonequilibrium along the porous surfaces of spacecraft liquid acquisition devices could cause the retention failures observed with liquid hydrogen and heated vapor pressurant. Why these devices seem immune to retention loss when pressurized with heated helium or heated directly through the porous structure was also examined. Results show that highly wetting fluids exhibit large negative and positive dynamic pressure gradients towards the meniscus interline when superheated and subcooled, respectively. With superheating, the pressure variation and recoil force arising from liquid/vapor phase change exert the same influence on surface morphology and promote retention. With subcooling, however, the pressure distribution produces a suction that degrades mechanical equilibrium of the surface. This result indicates that thermocapillary-induced deformation arising from subcooling and condensation is the likely cause for retention loss. In addition, increasing the level of nonequilibrium by reducing accommodation coefficient suppresses deformation and explains why this failure mode does not occur in instances of direct screen heating or pressurization with a heated inert gas.

Schmidt, George R.↗
NASA NTRS1995

Euler Technology Assessment program for preliminary aircraft design employing SPLITFLOW code with Cartesian unstructured grid method

This report documents results from the Euler Technology Assessment program. The objective was to evaluate the efficacy of Euler computational fluid dynamics (CFD) codes for use in preliminary aircraft design. Both the accuracy of the predictions and the rapidity of calculations were to be assessed. This portion of the study was conducted by Lockheed Fort Worth Company, using a recently developed in-house Cartesian-grid code called SPLITFLOW. The Cartesian grid technique offers several advantages for this study, including ease of volume grid generation and reduced number of cells compared to other grid schemes. SPLITFLOW also includes grid adaptation of the volume grid during the solution convergence to resolve high-gradient flow regions. This proved beneficial in resolving the large vortical structures in the flow for several configurations examined in the present study. The SPLITFLOW code predictions of the configuration forces and moments are shown to be adequate for preliminary design analysis, including predictions of sideslip effects and the effects of geometry variations at low and high angles of attack. The time required to generate the results from initial surface definition is on the order of several hours, including grid generation, which is compatible with the needs of the design environment.

Finley, Dennis B.↗
NASA NTRS1999

Gravitational Effects on the Morphology and Kinetics of Photodeposition of Polydiacetylene Thin Films From Monomer Solutions

The goal of this proposed work is to study gravitational effects on the photodeposition of polydiacetylene thin films from monomer solutions onto transparent substrates. Polydiacetylenes have been an extensively studied class of organic polymers because they exhibit many unusual and interesting properties, including electrical conductivity and optical nonlinearity. Their long polymeric chains render polydiacetylenes readily conducive to thin film formation, which is necessary for many applications. These applications require thin polydiacetylene films possessing uniform thicknesses, high purity, minimal inhomogeneities and defects (such as scattering centers), etc. Also, understanding and controlling the microstructure and morphology of the films is important for optimizing their electronic and optical properties. The lack of techniques for processing polydiacetylenes into such films has been the primary limitation to their commercial use. We have recently discovered a novel method for the formation of polydiacetylene thin films using photo-deposition from monomer solutions onto transparent substrates with UV light. This technique is very simple to carry out, and can yield films with superior quality to those produced by conventional methods. Furthermore, these films exhibit good third-order properties and are capable of waveguiding. We have been actively studying the chemistry of diacetylene polymerization in solution and the photo-deposition of polydiacetylene thin films from solution. It is well-known that gravitational factors such as buoyancy-driven convection and sedimentation can affect chemical and mass transport processes in solution. One important aspect of polydiacetylene thin film photodeposition in solution, relevant to microgravity science, is that heat generated by absorption of UV radiation induces thermal density gradients that under the influence of gravity, can cause fluid flows (buoyancy-driven convection). Additionally, changes in the chemical composition of the solution during polymerization may cause solutal convection. These fluid flows affect transport of material to and from the film surface and thereby affect the kinetics of the growth process. This manifests itself in the morphology of the resulting films; films grown under the influence of convection tend to have less uniform thicknesses, and can possess greater inhomogeneities and defects. Specifically, polydiacetylene films photodeposited from solution, when viewed under a microscope, exhibit very small particles of solid polymer which get transported by convection from the bulk solution to the surface of the growing film and become embedded. Even when carried out under conditions designed to minimize unstable density gradients (i.e., irradiating the solution from the top), some fluid flow still takes place (particles remain present in the films). It is also possible that defect nucleation may be occurring within the films or on the surface of the substrate; this, too, can be affected by convection (as is the case with crystal growth). Hence films grown in 1-g will, at best, still possess some defects. The objective of this proposal is to investigate, both in 1-g and in low-g, the effects of gravitational factors (primarily convection) on the dynamics of these processes, and on the quality, morphology, and properties of the films obtained.

Paley, Mark S.↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
NASA NTRS1992

Results of a zonally truncated three-dimensional model of the Venus middle atmosphere

Although the equatorial rotational speed of the solid surface of Venus is only 4 m s(exp-1), the atmospheric rotational speed reaches a maximum of approximately 100 m s(exp-1) near the equatorial cloud top level (65 to 70 km). This phenomenon, known as superrotation, is the central dynamical problem of the Venus atmosphere. We report here the results of numerical simulations aimed at clarifying the mechanism for maintaining the equatorial cloud top rotation. Maintenance of an equatorial rotational speed maximum above the surface requires waves or eddies that systematically transport angular momentum against its zonal mean gradient. The zonally symmetric Hadley circulation is driven thermally and acts to reduce the rotational speed at the equatorial cloud top level; thus wave or eddy transport must counter this tendency as well as friction. Planetary waves arising from horizontal shear instability of the zonal flow (barotropic instability) could maintain the equatorial rotation by transporting angular momentum horizontally from midlatitudes toward the equator. Alternatively, vertically propagating waves could provide the required momentum source. The relative motion between the rotating atmosphere and the pattern of solar heating, which as a maximum where solar radiation is absorbed near the cloud tops, drives diurnal and semidiurnal thermal tides that propagate vertically away from the cloud top level. The effect of this wave propagation is to transport momentum toward the cloud top level at low latitudes and accelerate the mean zonal flow there. We employ a semispectral primitive equation model with a zonal mean flow and zonal wavenumbers 1 and 2. These waves correspond to the diurnal and semidiurnal tides, but they can also be excited by barotropic or baroclinic instability. Waves of higher wavenumbers and interactions between the waves are neglected. Symmetry about the equator is assumed, so the model applies to one hemisphere and covers the altitude range 30 to 110 km. Horizontal resolution is 1.5 deg latitude, and vertical resolution is 1.5 km. Solar and thermal infrared heating, based on Venus observations and calculations drive the model flow. Dissipation is accomplished mainly by Rayleigh friction, chosen to produce strong dissipation above 85 km in order to absorb upward propagating waves and limit extreme flow velocities there, yet to give very weak Rayleigh friction below 70 km; results in the cloud layer do not appear to be sensitive to the Rayleigh friction. The model also has weak vertical diffusion, and very weak horizontal diffusion, which has a smoothing effect on the flow only at the two grid points nearest the pole.

Newman, M.↗
NASA NTRS1994

A New K-epsilon Eddy Viscosity Model for High Reynolds Number Turbulent Flows: Model Development and Validation

A new k-epsilon eddy viscosity model, which consists of a new model dissipation rate equation and a new realizable eddy viscosity formulation, is proposed. The new model dissipation rate equation is based on the dynamic equation of the mean-square vorticity fluctuation at large turbulent Reynolds number. The new eddy viscosity formulation is based on the realizability constraints: the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. We find that the present model with a set of unified model coefficients can perform well for a variety of flows. The flows that are examined include: (1) rotating homogeneous shear flows; (2) boundary-free shear flows including a mixing layer, planar and round jets; (3) a channel flow, and flat plate boundary layers with and without a pressure gradient; and (4) backward facing step separated flows. The model predictions are compared with available experimental data. The results from the standard k-epsilon eddy viscosity model are also included for comparison. It is shown that the present model is a significant improvement over the standard k-epsilon eddy viscosity model.

Shih, T.-H.↗
NASA NTRS2017

Reconfigurable Model Execution in the OpenMDAO Framework

NASA's OpenMDAO framework facilitates constructing complex models and computing their derivatives for multidisciplinary design optimization. Decomposing a model into components that follow a prescribed interface enables OpenMDAO to assemble multidisciplinary derivatives from the component derivatives using what amounts to the adjoint method, direct method, chain rule, global sensitivity equations, or any combination thereof, using the MAUD architecture. OpenMDAO also handles the distribution of processors among the disciplines by hierarchically grouping the components, and it automates the data transfer between components that are on different processors. These features have made OpenMDAO useful for applications in aircraft design, satellite design, wind turbine design, and aircraft engine design, among others. This paper presents new algorithms for OpenMDAO that enable reconfigurable model execution. This concept refers to dynamically changing, during execution, one or more of: the variable sizes, solution algorithm, parallel load balancing, or set of variables-i.e., adding and removing components, perhaps to switch to a higher-fidelity sub-model. Any component can reconfigure at any point, even when running in parallel with other components, and the reconfiguration algorithm presented here performs the synchronized updates to all other components that are affected. A reconfigurable software framework for multidisciplinary design optimization enables new adaptive solvers, adaptive parallelization, and new applications such as gradient-based optimization with overset flow solvers and adaptive mesh refinement. Benchmarking results demonstrate the time savings for reconfiguration compared to setting up the model again from scratch, which can be significant in large-scale problems. Additionally, the new reconfigurability feature is applied to a mission profile optimization problem for commercial aircraft where both the parametrization of the mission profile and the time discretization are adaptively refined, resulting in computational savings of roughly 10% and the elimination of oscillations in the optimized altitude profile.

Hwang, John T.↗
DOE OSTI2021

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗