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 181 records · Page 10

Efficient Neural Network Approaches for Conditional Optimal Transport with Applications in Bayesian Inference

In this work, we present two neural network approaches that approximate the solutions of static and dynamic conditional optimal transport (COT) problems. Both approaches enable conditional sampling and conditional density estimation, which are core tasks in Bayesian inference—particularly in the simulation-based (“likelihood-free”) setting. Our methods represent the target conditional distribution as a transformation of a tractable reference distribution. Obtaining such a transformation, chosen here to be an approximation of the COT map, is computationally challenging even in moderate dimensions. To improve scalability, our numerical algorithms use neural networks to parameterize candidate maps and further exploit the structure of the COT problem. Our static approach approximates the map as the gradient of a partially input convex neural network. It uses a novel numerical implementation to increase computational efficiency compared to state-of-the-art alternatives. Our dynamic approach approximates the conditional optimal transport via the flow map of a regularized neural ODE; compared to the static approach, it is slower to train but offers more modeling choices and can lead to faster sampling. We demonstrate both algorithms numerically, comparing them with competing state-of-the-art approaches, using benchmark datasets and simulation-based Bayesian inverse problems.

97 MATHEMATICS AND COMPUTING↗

The effect of gravitational modulation on convection in vertical Bridgman growth

During vertical directional solidification of a binary alloy at constant velocity, buoyancy-driven solutal convection may occur due to the solute gradient associated with the solidification process. This problem is further complicated if time-periodic forcing is considered, which is relevant to materials processing in a microgravity environment or as a means of dynamic control of flow instabilities. The effect of time-periodic modulation is studied by introducing a gravitational acceleration which is a sinusoidal function of time. The onset of solutal convection is treated by a stability analysis of the linearized governing equations and boundary conditions. Solutions are obtained numerically employing two distinct computational implementations of Floquet theory. Results are presented for materials with large Schmidt number, and an analysis for large frequency and large Schmidt number yields a more complete description of the behavior in this relevant limit.

Murray, B. T.↗

Tip-Clearance Vortex Characterized With Three-Dimensional Digital Particle Image Velocimetry

An optical measurement technique known as Three-Dimensional Digital Particle Image Velocimetry (3-D DPIV) was used to characterize the tip clearance flow in NASA Glenn Research Center's low-speed axial compressor. 3-D DPIV is a technique in which a stereoscopic imaging system consisting of two cross-correlation cameras is used to record particles entrained in a flow as a laser light sheet is pulsed at two instances in time. Although 3-D DPIV has been used elsewhere, this is the first time it has been used to measure compressor tip clearance flows. In-house modifications of the DPIV system include the use of effective seeding technology and a novel system to perform a priori calibrations at all five measurement planes, greatly reducing facility run time. Computational fluid dynamics predictions, which are used to guide design changes toward improving the efficiency and operating range of turbomachinery, can be verified and improved by comparison with 3-D DPIV measurements of the actual tip clearance flow. This measurement campaign dealt with the characterization of the tip clearance vortex in the first stage of a four-stage axial compressor. The tip clearance vortex is formed in compressors operating with a clearance gap between the moving rotor blade tips and the stationary casing when a leakage flow, forced from the pressure side of the blade over the blade tip, forms a vortical structure on the suction side of the blade. 3-D DPIV is ideally suited to measure the clearance vortex for two reasons: (1) this technique captures the entire blade passage flow at one instant in time, so that wandering of the vortex during the measurement does not smear out velocity gradients in the flow field, and (2) the spanwise component of velocity changes sign across the vortex core, providing a more accurate measurement of the vortex location than was available with previous two-dimensional measurement approaches. These two attributes of the data will enable computational fluid dynamics researchers to validate their predictions to a level of accuracy not previously attainable. In an effort to understand the effects of this flow phenomenon on the operation of the low speed axial compressor, data were acquired at two mass flow coefficients: 0.395 (design operating point) and 0.35 (operating point just above stall). In order to identify the spatial extent, location, and magnitude of the tip clearance vortex for each mass flow condition, data were acquired at five equally spaced spanwise locations, from 90- to 100-percent span. The data presented shows a cross section of the tip clearance vortex at 94-percent span and a mass flow coefficient of 0.395, where the color contours represent the radial velocity component. The blue (flow towards the hub) and red (flow towards the casing) regions illustrate the extent of the tip clearance vortex. The interface between the blue and red regions delineates the location of the vortex core.

John, W. Trevor↗

Potential vorticity structure of Titan's polar vortices from Cassini CIRS observations

The Cassini mission has provided the best opportunity to date to extensively study the seasonal variation in Titan’s atmosphere, with observations spanning almost half a Titan year (𝐿𝑠 = 293 − 93°). An important feature in the Titan middle-atmosphere is the formation of a polar vortex. Observations have shown that an initially well-developed northern vortex enriched with trace gas species gradually breaks down after spring equinox as a new vortex emerges in southern winter. Here we use Cassini CIRS observations to derive the temperature and composition of the middle-atmosphere. We use the gradient wind equation to first estimate the mean zonal winds, and then the Potential Vorticity (PV) throughout Titan’s atmosphere over the timespan of the Cassini mission. PV is a useful diagnostic quantity for studying the dynamics of polar vortices because it is materially conserved for adiabatic and frictionless flows, and can be inverted to find all other dynamical fields. Our results show the formation of a strong zonal jet in the winter hemisphere, with wind velocities reaching 220 m/s, which is consistent with previous studies. An annular PV structure is also observed over the winter poles, whereby a ring of PV encircles a local minima over the pole. Such distributions are often found to be unstable without a restoring force, yet they are seen here in numerous observations in both the northern and southern hemispheres. A comparison with the annular Martian vortices shows that latent heat release from condensation or subsidence-induced adiabatic heating may explain the origin and stability of the annulus. Finally, we investigate the evolution of the size of the vortices and the role of strong PV gradients as a dynamical mixing barrier for trace gas species across the vortex edge. We find that longer lived gases are less confined to the vortex than those with shorter photochemical lifetimes.

Conor Andrew Nixon↗

Phase-space dynamics of Alfvén mode chirping

In tokamak discharges, toroidal Alfvén eigenmodes often experience complex semi-periodic frequency modulation known as chirping. These events modify the local high energy particle distribution and are expected to occur in many future fusion devices, which include energetic beams or fusion products. In this paper, we present a study of simulations of mode chirping made in order to better understand its phase-space properties in a realistic tokamak configuration. We find a mechanism that permits rapid repeated chirping with strong amplitude variation in each chirp. Each chirp is then associated with an amplitude crash to low magnitude and local manipulation of the density gradients through a shift of mode phase through π. The chirping produces high density clumps, which propagate down the fast ion density gradient and low density holes that propagate up the density gradient away from the resonance. This flow of particles across the resonance provides an energy source and local gradients for repeated chirping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear dynamics and transport driven by energetic particle instabilities using a gyro-Landau closure model *

Abstract Energetic particle (EP) destabilized Alfvén eigenmode (AE) instabilities are simulated for a DIII-D experimental case with a pulsed neutral beam using a gyro-Landau moments model which introduces EP phase-mixing effects through closure relations. This provides a computationally efficient reduced model which is applied here in the nonlinear regime over timescales that would be difficult to address with more complete models. The long timescale nonlinear evolution and related collective transport losses are examined including the effects of zonal flow/current generation, nonlinear energy cascades, and EP profile flattening. The model predicts frequencies and mode structures that are consistent with experimental observations. These calculations address issues that have not been considered in previous modelling: the EP critical gradient profile evolution in the presence of zonal flows/currents, and the dynamical nature of the saturated state. A strong level of intermittency is present in the predicted instability-driven transport; this is connected to the zonal flow growth and decay cycles and nonlinear energy transfers. Simulation of intermittent AE-enhanced EP transport will be an important issue for the protection of plasma facing components in the next generation of fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing initial baryon stopping and equation of state with rapidity-dependent directed flow of identified particles

In this work, using a (3 + 1)-dimensional hybrid framework with parametric initial conditions, we study the rapidity-dependent directed flow $v_1 (y)$ of identified particles, including pions, kaons, protons, and lambdas in heavy-ion collisions. Cases involving Au + Au collisions are considered, performed at $\sqrt {^SNN}$ ranging from 7.7 to 200 GeV. The dynamics in the beam direction is constrained using the measured pseudo-rapidity distribution of charged particles and the net proton rapidity distribution. Within this framework, the directed flow of mesons is driven by the sideward pressure gradient from the tilted source, and that of baryons mainly due to the initial asymmetric baryon distribution with respect to the beam axis driven by the transverse expansion. Our approach successfully reproduces the rapidity- and beam energy-dependence of $v_1$ for both mesons and baryons. We find that the $v_1 (y)$ of baryons has strong constraining power on the initial baryon stopping, and together with that of mesons, the directed flow probes the equation of state of the dense nuclear matter at finite chemical potentials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fluid dynamics of unsteady separated flow. II - Lifting surfaces

An analytic method is described which uses static experimental data to predict the separated flow effect on rigid and elastic vehicle dynamics. Key parameters in the analytic relationship between steady and nonsteady aerodynamics are: the time lag occurring before a change of flow conditions can affect the separation-induced aerodynamic loads; the accelerated flow effect (i.e., the pressure gradient lag relative to the static aerodynamic characteristics); and the moving wall effect (i.e., the effect of the nonsteady boundary condition at the vehicle surface). Using the existing experimental data base, an analytic theory is formulated that can predict the separation-induced unsteady aerodynamics if the static characteristics are known from theory or experiment. Reference is made to increased-maneuverability advanced aircraft and to Space Shuttle Orbiter aerodynamics.

Ericsson, L. E.↗

Mean Flow Augmented Acoustics in Rocket Systems

Combustion instability in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. Recent advances in energy based modeling of combustion instabilities require accurate determination of acoustic frequencies and mode shapes. Of particular interest is the acoustic mean flow interactions within the converging section of a rocket nozzle, where gradients of pressure, density, and velocity become large. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The present study aims to implement the French model within the COMSOL Multiphysiscs framework and analyzes one of the author's presented test cases.

Fischbach, Sean↗

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Vibration isolation technology: Sensitivity of selected classes of experiments to residual accelerations

The solution was sought of a 2-D axisymmetric moving boundary problem for the sensitivity of isothermal and nonisothermal liquid columns and the sensitivity of thermo-capillary flows to buoyancy driven convection caused by residual accelerations. The sensitivity of a variety of space experiments to residual accelerations are examined. In all the cases discussed, the sensitivity is related to the dynamic response of a fluid. In some cases the sensitivity can be defined by the magnitude of the response of the velocity field. This response may involve motion of the fluid associated with internal density gradients, or the motion of a free liquid surface. For fluids with internal density gradients, the type of acceleration to which the experiment is sensitive will depend on whether buoyancy driven convection must be small in comparison to other types of fluid motion (such as thermocapillary flow), or fluid motion must be suppressed or eliminated (such as in diffusion studies, or directional solidification experiments). The effect of the velocity on the composition and temperature field must be considered, particularly in the vicinity of the melt crystal interface. As far as the response to transient disturbances is concerned the sensitivity is determined by both the magnitude and frequency the acceleration and the characteristic momentum and solute diffusion times.

Alexander, J. Iwan D.↗

Model-Free Primal-Dual Methods for Network Optimization with Application to Real-Time Optimal Power Flow

This paper examines the problem of real-time optimization of networked systems and develops online algorithms that steer the system towards the optimal trajectory without explicit knowledge of the system model. The problem is modeled as a dynamic optimization problem with time-varying performance objectives and engineering constraints. The design of the algorithms leverages the online zero-order primal-dual projected-gradient method. In particular, the primal step that involves the gradient of the objective function (and hence requires a networked systems model) is replaced by its zero-order approximation with two function evaluations using a deterministic perturbation signal. The evaluations are performed using the measurements of the system output, hence giving rise to a feedback interconnection, with the optimization algorithm serving as a feedback controller. The paper provides some insights on the stability and tracking properties of this interconnection. Finally, the paper applies this methodology to a real-time optimal power flow problem in power systems, and shows its efficacy on the IEEE 37-node distribution test feeder for reference power tracking and voltage regulation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Dynamic behavior of an unsteady trubulent boundary layer

Experiments on an unsteady turbulent boundary layer are reported in which the upstream portion of the flow is steady (in the mean) and in the downstream region, the boundary layer sees a linearly decreasing free stream velocity. This velocity gradient oscillates in time, at frequencies ranging from zero to approximately the bursting frequency. For the small amplitude, the mean velocity and mean turbulence intensity profiles are unaffected by the oscillations. The amplitude of the periodic velocity component, although as much as 70% greater than that in the free stream for very low frequencies, becomes equal to that in the free stream at higher frequencies. At high frequencies, both the boundary layer thickness and the Reynolds stress distribution across the boundary layer become frozen. The behavior at higher amplitude is quite similar. At sufficiently high frequencies, the boundary layer thickness remains frozen at the mean value over the oscillation cycle, even though flow reverses near the wall during a part of the cycle.

Parikh, P. G.↗

Dynamic behavior of an unsteady turbulent boundary layer

Experiments on an unsteady turbulent boundary layer are reported in which the upstream portion of the flow is steady (in the mean) and in the downstream region, the boundary layer sees a linearly decreasing free stream velocity. This velocity gradient oscillates in time, at frequencies ranging from zero to approximately the bursting frequency. For the small amplitude, the mean velocity and mean turbulence intensity profiles are unaffected by the oscillations. The amplitude of the periodic velocity component, although as much as 70 percent greater than that in the free stream for very low frequencies, becomes equal to that in the free stream at higher frequencies. At high frequencies, both the boundary layer thickness and the Reynolds stress distribution across the boundary layer become frozen. The behavior at higher amplitude is quite similar. At sufficiently high frequencies, the boundary layer thickness remains frozen at the mean value over the oscillation cycle, even though flow reverses near the wall during a part of the cycle. Previously announced in STAR as N81-28392

Parikh, P. G.↗

Model-Free Primal-Dual Methods for Network Optimization with Application to Real-Time Optimal Power Flow: Preprint

This paper examines the problem of real-time optimization of networked systems and develops online algorithms that steer the system towards the optimal trajectory without explicit knowledge of the system model. The problem is modeled as a dynamic optimization problem with time-varying performance objectives and engineering constraints. The design of the algorithms leverages the online zero-order primal-dual projected-gradient method. In particular, the primal step that involves the gradient of the objective function (and hence requires networked systems model) is replaced by its zero-order approximation with two function evaluations using a deterministic perturbation signal. The evaluations are performed using the measurements of the system output, hence giving rise to a feedback interconnection, with the optimization algorithm serving as a feedback controller. The paper provides some insights on the stability and tracking properties of this interconnection. Finally, the paper applies this methodology to a real-time optimal power flow problem in power systems, and shows its efficacy on the IEEE 37-node distribution test feeder for reference power tracking and voltage regulation.

61 RADIATION PROTECTION AND DOSIMETRY↗

Three-Dimensional Viscous Alternating Direction Implicit Algorithm and Strategies for Shape Optimization

A gradient-based shape optimization based on quasi-analytical sensitivities has been extended for practical three-dimensional aerodynamic applications. The flow analysis has been rendered by a fully implicit, finite-volume formulation of the Euler and Thin-Layer Navier-Stokes (TLNS) equations. Initially, the viscous laminar flow analysis for a wing has been compared with an independent computational fluid dynamics (CFD) code which has been extensively validated. The new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4 with coarse- and fine-grid based computations performed with Euler and TLNS equations. The influence of the initial constraints on the geometry and aerodynamics of the optimized shape has been explored. Various final shapes generated for an identical initial problem formulation but with different optimization path options (coarse or fine grid, Euler or TLNS), have been aerodynamically evaluated via a common fine-grid TLNS-based analysis. The initial constraint conditions show significant bearing on the optimization results. Also, the results demonstrate that to produce an aerodynamically efficient design, it is imperative to include the viscous physics in the optimization procedure with the proper resolution. Based upon the present results, to better utilize the scarce computational resources, it is recommended that, a number of viscous coarse grid cases using either a preconditioned bi-conjugate gradient (PbCG) or an alternating-direction-implicit (ADI) method, should initially be employed to improve the optimization problem definition, the design space and initial shape. Optimized shapes should subsequently be analyzed using a high fidelity (viscous with fine-grid resolution) flow analysis to evaluate their true performance potential. Finally, a viscous fine-grid-based shape optimization should be conducted, using an ADI method, to accurately obtain the final optimized shape.

Pandya, Mohagna J.↗

Electric field and turbulence in global Braginskii simulations across the ASDEX Upgrade edge and scrape-off layer

Turbulence simulations in diverted geometry across the edge and scrape-off layer (SOL) of ASDEX Upgrade are performed with the GRILLIX code (Stegmeir et al 2019 Phys. Plasmas 26 052517). The underlying global (full-f) drift-reduced Braginskii model allows to concurrently study the self-consistent dynamics of the turbulence and the background as well as the evolution of toroidal and zonal flows. Different contributions to the radial electric field are identified. The dominant contribution on closed flux surfaces comes from the ion pressure gradient, due to the diamagnetic drift in the curved magnetic field. Large deviations can be induced, in particular, by the polarization particle flux, leading to zonal flows. The latter are driven by small-scale eddies, but do not exhibit much impact on the overall transport which is driven by ballooning modes at larger scales. Ion viscosity is found to be important in damping poloidal rotation through adjusting of the parallel velocity profile, but not via direct vorticity damping. The zonal flow drive peaks at the separatrix, where a strong shear layer forms due to the sheath-induced counter-propagating SOL flow, allowing for the formation of a transport barrier. The temperature profile across the separatrix is determined by the competition between cross-field transport and outflow in the SOL, the latter being largely controlled by the parallel heat conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Double diffusive convection during directional solidification - Nonlinear analysis

The dynamics of doubly diffusive convective flow during the constant-velocity vertical solidification of binary melts to form single-phase solids is examined, summarizing the results of recent experimental and analytical investigations. The effects of horizontal temperature gradients are considered; the nonlinearities of the system and their numerical treatment are discussed; and numerical results are presented graphically for the solidification of a Pb melt containing 0.0015 wt pct Sn, with growth velocity 0.0002 cm/s and a temperature gradient of 200 K/cm in the liquid.

Coriell, S. R.↗