Search NASA⌕ Search

SEARCH · Search NASA

Results for “Field Theory Hydrodynamics”

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 37 records · Page 2

The interaction of radio frequency electromagnetic fields with atmospheric water droplets and applications to aircraft ice prevention

The feasibility of computerized simulation of the physics of advanced microwave anti-icing systems, which preheat impinging supercooled water droplets prior to impact, was investigated. Theoretical and experimental work performed to create a physically realistic simulation is described. The behavior of the absorption cross section for melting ice particles was measured by a resonant cavity technique and found to agree with theoretical predictions. Values of the dielectric parameters of supercooled water were measured by a similar technique at lambda = 2.82 cm down to -17 C. The hydrodynamic behavior of accelerated water droplets was studied photograhically in a wind tunnel. Droplets were found to initially deform as oblate spheroids and to eventually become unstable and break up in Bessel function modes for large values of acceleration or droplet size. This confirms the theory as to the maximum stable droplet size in the atmosphere. A computer code which predicts droplet trajectories in an arbitrary flow field was written and confirmed experimentally. The results were consolidated into a simulation to study the heating by electromagnetic fields of droplets impinging onto an object such as an airfoil. It was determined that there is sufficient time to heat droplets prior to impact for typical parameter values. Design curves for such a system are presented.

Hansman, R. J., Jr.↗

Thermocapillary-Induced Phase Separation with Coalescence

Research has been undertaken on interactions of two or more deformable drops (or bubbles) in a viscous fluid and subject to a temperature, gravitational, or flow field. An asymptotic theory for nearly spherical drops shows that small deformations reduce the coalescence and phase separation rates. Boundary-integral simulations for large deformations show that bubbles experience alignment and enhanced coalescence, whereas more viscous drops may break as a result of hydrodynamic interactions. Experiments for buoyancy motion confirm these observations. Simulations of the sedimentation of many drops show clustering phenomena due to deformations, which lead to enhanced phase separation rates, and simulations of sheared emulsions show that deformations cause a reduction in the effective viscosity.

Davis, Robert H.↗

Theory of electrohydrodynamic instabilities in electrolytic cells

The paper develops the theory of the hydrodynamic stability of an electrolytic cell as a function of the imposed electric current. A new electrohydrodynamic instability is encountered when the current is forced to exceed the Nernst limit. The convection is driven by the volume force exerted by the electric field on space charges in the electrolyte. This intrinsic instability is found to be easily masked by extrinsic convection sources such as gravity or stirring. A linear stability analysis is performed and a dimensionless number Le is derived whose value determines the convection pattern.

Bruinsma, R.↗

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗

Magnetic dynamo action in two-dimensional turbulent magneto-hydrodynamics

Two-dimensional magnetohydrodynamic turbulence is explored by means of numerical simulation. Previous analytical theory, based on non-dissipative constants of the motion in a truncated Fourier representation, is verified by following the evolution of highly non-equilibrium initial conditions numerically. Dynamo action (conversion of a significant fraction of turbulent kinetic energy into long-wavelength magnetic field energy) is observed. It is conjectured that in the presence of dissipation and external forcing, a dual cascade will be observed for zero-helicity situations. Energy will cascade to higher wave numbers simultaneously with a cascade of mean square vector potential to lower wave numbers, leading to an omni-directional magnetic energy spectrum which varies as 1/k 3 at lower wave numbers, simultaneously with a buildup of magnetic excitation at the lowest wave number of the system. Equipartition of kinetic and magnetic energies is expected at the highest wave numbers in the system.

Fyfe, D.↗

Magnetic dynamo action in two-dimensional turbulent magneto-hydrodynamics

Two-dimensional magnetohydrodynamic turbulence is explored by means of numerical simulation. Previous analytical theory, based on non-dissipative constants of the motion in a truncated Fourier representation, is verified by following the evolution of highly non-equilibrium initial conditions numerically. Dynamo action (conversion of a significant fraction of turbulent kinetic energy into long-wavelength magnetic field energy) is observed. It is conjectured that in the presence of dissipation and external forcing, a dual cascade will be observed for zero-helicity situations. Energy will cascade to higher wavenumbers simultaneously with a cascade of mean square vector potential to lower wavenumbers, leading to an omni-directional magnetic energy spectrum.

Fyfe, D.↗

Shear flow-ballooning instability as a possible mechanism for hydromagnetic fluctuations

A unified linear electromagnetic analysis of both the Kelvin-Helmholtz (shear flow) instability and of the ballooning (interchange) instability is carried out on the basis of MHD theory. In the analysis, the concept of the Richardson instability of hydrodynamic flows is extended into the hydromagnetic context by unifying both the shear flow and the ballooning instability. As essential concept of the analysis is the role played by the magnetic buoyancy due to an effective gravity produced by the curvature of the field lines which provides the basic step by which both instabilities could be coupled. The results of the study are applied to the plasmapause to explain the excitation of hydromagnetic waves in that region, including the effect of the hot particles from the plasma sheet.

Vinas, A. F.↗

Energy transport and dynamics

We report findings concerning energy transport and dynamics in flares during the impulsive and gradual phases based on new ground-based and space observations (notably from Yohkoh). A preheating sometimes occurs during the impulsive phase. Ca XIX line shifts are confirmed to be good tracers of bulk plasma motions, although strong blue shifts are not as frequent as previously claimed. They often appear correlated with hard X-rays but, for some events, the concept that electron beams provide the whole energy input to the thermal component seems not to apply. Theory now yields: new diagnostics of low-energy proton and electric beams; accurate hydrodynamical modeling of pulse beam heating of the atmosphere; possible diagnostics of microflares (based on X-ray line ratio or on loop variability); and simulated images of chromospheric evaporation fronts. For the gradual phase, the continual reorganization of magnetic field lines over active regions determines where and when magnetic reconnection, the mechanism favored for energy release, will occur. Spatial and temporal fragmentation of the energy release, observed at different wavelengths, is considered to be a factor as well in energy transport and plasma dynamics.

Schmieder, Brigitte↗

Trailing edge noise prediction from measured surface pressures

Trailing edge (TE) noise is investigated for the case of a two-dimensional airfoil embedded in a uniform low Mach number flow, and the usefulness of several TE noise theories is examined by applying them to the measured data. The TE noise spectra and directivity are quantitatively determined for the case of a high Reynolds number and a fully turbulent boundary layer. Parameters include angle of attack, flow velocity and TE bluntness. Evanescent wave theories by Chase (1975) and Chandiramani (1974) are compared to the present results and show good agreement. Agreement of the near field pressure scatter phenomenon analysis with measurements implies that the basic assumptions used in the analysis are correct, i.e., the turbulent boundary layer (TBL) flow passes the trailing edge into the wake region. No hydrodynamic wake shedding activity is confirmed for the two-sided TBL flow, and a method incorporating the principles of the coherent output power method is used to determine the sound field. The near field edge scatter model is found to not only establish optimum sizing of edge treatment for noise control, but also to separate and identify the scattered field from the incident hydrodynamic field.

Brooks, T. F.↗

Inverse energy cascades in three-dimensional turbulence

Fully three-dimensional magnetohydrodynamic (MHD) turbulence at large kinetic and low magnetic Reynolds numbers is considered in the presence of a strong uniform magnetic field. It is shown by numerical simulation of a model of MHD that the energy inverse cascades to longer length scales when the interaction parameter is large. While the steady-state dynamics of the driven problem is three-dimensional in character, the behavior has resemblance to two-dimensional hydrodynamics. These results have implications in turbulence theory, MHD power generator, planetary dynamos, and fusion reactor blanket design.

Hossain, Murshed↗

A First-Principles Model Based on Saturation of the Electron Cyclotron Drift Instability for Electron Transport in Hydrodynamics Simulations of Hall Thruster Plasmas

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Lopez Ortega, Alejandro↗

An Experimental Investigation of Incompressible Richtmyer-Meshkov Instability

Richtmyer-Meshkov (RM) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and eventually become a turbulent flow. It is closely related to Rayleigh-Taylor instability, which is the instability of a planar interface undergoing constant acceleration, such as caused by the suspension of a heavy fluid over a lighter one in the earth's gravitational field. Like the well-known Kelvin-Helmholtz instability, RM instability is a fundamental hydrodynamic instability which exhibits many of the nonlinear complexities that transform simple initial conditions into a complex turbulent flow. Furthermore, the simplicity of RM instability (in that it requires very few defining parameters), and the fact that it can be generated in a closed container, makes it an excellent test bed to study nonlinear stability theory as well as turbulent transport in a heterogeneous system. However, the fact that RM instability involves fluids of unequal densities which experience negligible gravitational force, except during the impulsive acceleration, requires RM instability experiments to be carried out under conditions of microgravity. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The impulsive acceleration is generated by bouncing a rectangular tank containing two different density liquids off a retractable vertical spring. The initial perturbation is produced prior to release by oscillating the tank in the horizontal direction to produce a standing wave. The instability evolves in microgravity as the tank travels up and then down the vertical rails of a drop tower until hitting a shock absorber at the bottom. Planar Laser Induced Fluorescence (PLIF) is employed to visualize the flow. PLIF images are captured by a video camera that travels with the tank. Figure 1 is as sequence of images showing the development of the instability from the initial sinusoidal disturbance far into the nonlinear regime which is characterized by the appearance of mushroom structures resulting from the coalescence of baroclinic vorticity produced by the impulsive acceleration. At later times in this sequence the vortex cores are observed to become unstable showing the beginnings of the transition to turbulence in this flow. The amplitude of the growing disturbance after the impulsive acceleration is measured and found to agree well with theoretical predictions. The effects of Reynolds number (based on circulation) on the development of the vortices and the transition to turbulence are also determined.

Jacobs, J. W.↗

A High-Resolution X-Ray and Optical Study of SN1006: Asymmetric Expansion and Small-Scale Structure in a Type Ia Supernova Remnant

We introduce a deep (670 ks) X-ray survey of the entire SN 1006 remnant from the Chandra X-Ray Observatory, together with a deep Hα image of SN 1006 from the 4 m Blanco telescope at CTIO. Comparison with Chandra images from 2003 gives the first measurement of the X-ray proper motions around the entire periphery, carried out over a 9 yr baseline. We find that the expansion velocity varies significantly with azimuth. The highest velocity of approx.7400 km/s (almost 2.5 times that in the northwest (NW)) is found along the southeast (SE) periphery, where both the kinematics and the spectra indicate that most of the X-ray emission stems from ejecta that have been decelerated little, if at all. Asymmetries in the distribution of ejecta are seen on a variety of spatial scales. Si-rich ejecta are especially prominent in the SE quadrant, while O and Mg are more uniformly distributed, indicating large-scale asymmetries arising from the explosion itself. Neon emission is strongest in a sharp filament just behind the primary shock along the NWrim, where the pre-shock density is highest. Here the Ne is likely interstellar, while Ne within the shell may include a contribution from ejecta. Within the interior of the projected shell we find a few isolated "bullets" of what appear to be supernova ejecta that are immediately preceded by bowshocks seen in Hα, features that we interpret as ejecta knots that have reached relatively dense regions of the surrounding interstellar medium, but that appear in the interior in projection. Recent three-dimensional hydrodynamic models for Type Ia supernovae display small-scale features that strongly resemble the ones seen in X-rays in SN 1006; an origin in the explosion itself or from subsequent hydrodynamic instabilities both remain viable options. We have expanded the search for precursor X-ray emission ahead of a synchrotron-dominated shock front, as expected from diffusive shock acceleration theory, to numerous regions along both the northeast and southwest rims of the shell. Our data require that a precursor be thinner than about 3, and fainter than about 5% of the post-shock peak. These limits suggest that the magnetic field is amplified by a factor of seven or more in a narrow precursor region, promoting diffusive particle acceleration.

galaxies: active↗

Inflow/outflow boundary conditions and global dynamics of spatial mixing layers

The numerical simulation of incompressible spatially-developing shear flows poses a special challenge to computational fluid dynamicists. The Navier-Stokes equations are elliptic and boundary equations need to be specified at the inflow and outflow boundaries in order to compute the fluid properties within the region of interest. It is, however, difficult to choose inflow and outflow conditions corresponding to a given experimental situation. Furthermore the effects that changes in the boundary conditions or in the size of the computational domain may induce on the global dynamics of the flow are presently unknown. These issues are examined in light of recent developments in hydrodynamic stability theory. The particular flow considered is the spatial mixing layer but it was expected that similar phenomena were bound to occur in other cases such as channel flow, the boundary layer, etc. A short summary of local/global and absolute/convective instability concepts is given. The results of numerical simulations are presented which strongly suggest that global resonances may be triggered in domains of finite streamwise extent although the evolution of the perturbation vorticity field is everywhere locally convective. A relationship between finite domains and pressure sources which might help in devising a scheme to eliminate these difficulties is discussed.

Buell, J. C.↗

On the relation between 'mixing length' and 'direct interaction approximation' theories of turbulence

The capabilities of direct interaction approximations (DIA) and weak-coupling approximations (WCA) in modeling turbulence in plasma physics studies are examined. DIA and WCA formalisms are reviewed and a DIA expression for homogeneous incompressible hydrodynamics is simplified to derive an energy spectrum in an arbitrary wavenumber space. A solution is defined for the expression when applied to long wavelength phenomena generated by an instability. The equation is further generalized to yield a transport equation for modeling the spectrum generated by homogeneous turbulence in a linearly unstable system undergoing nonlinear interaction. Finally, the technique is applied to two-dimensional turbulence in a low pressure, weakly ionized plasma embedded in a homogeneous magnetic field and experiencing instabilities caused by electron density and electric potential gradients, such as found in the ionospheric E region.

Sudan, R. N.↗

The Cool Flames Experiment

A space-based experiment is currently under development to study diffusion-controlled, gas-phase, low temperature oxidation reactions, cool flames and auto-ignition in an unstirred, static reactor. At Earth's gravity (1g), natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles via the Arrhenius temperature dependence of the reaction rates. Natural convection is important in all terrestrial cool flame and auto-ignition studies, except for select low pressure, highly dilute (small temperature excess) studies in small vessels (i.e., small Rayleigh number). On Earth, natural convection occurs when the Rayleigh number (Ra) exceeds a critical value of approximately 600. Typical values of the Ra, associated with cool flames and auto-ignitions, range from 104-105 (or larger), a regime where both natural convection and conduction heat transport are important. When natural convection occurs, it alters the temperature, hydrodynamic, and species concentration fields, thus generating a multi-dimensional field that is extremely difficult, if not impossible, to be modeled analytically. This point has been emphasized recently by Kagan and co-workers who have shown that explosion limits can shift depending on the characteristic length scale associated with the natural convection. Moreover, natural convection in unstirred reactors is never "sufficiently strong to generate a spatially uniform temperature distribution throughout the reacting gas." Thus, an unstirred, nonisothermal reaction on Earth does not reduce to that generated in a mechanically, well-stirred system. Interestingly, however, thermal ignition theories and thermokinetic models neglect natural convection and assume a heat transfer correlation of the form: q=h(S/V)(T(bar) - Tw) where q is the heat loss per unit volume, h is the heat transfer coefficient, S/V is the surface to volume ratio, and (T(bar) - Tw ) is the spatially averaged temperature excess. This Newtonian form has been validated in spatially-uniform, well-stirred reactors, provided the effective heat transfer coefficient associated with the unsteady process is properly evaluated. Unfortunately, it is not a valid assumption for spatially-nonuniform temperature distributions induced by natural convection in unstirred reactors. "This is why the analysis of such a system is so difficult." Historically, the complexities associated with natural convection were perhaps recognized as early as 1938 when thermal ignition theory was first developed. In the 1955 text "Diffusion and Heat Exchange in Chemical Kinetics", Frank-Kamenetskii recognized that "the purely conductive theory can be applied at sufficiently low pressure and small dimensions of the vessel when the influence of natural convection can be disregarded." This was reiterated by Tyler in 1966 and further emphasized by Barnard and Harwood in 1974. Specifically, they state: "It is generally assumed that heat losses are purely conductive. While this may be valid for certain low pressure slow combustion regimes, it is unlikely to be true for the cool flame and ignition regimes." While this statement is true for terrestrial experiments, the purely conductive heat transport assumption is valid at microgravity (mu-g). Specifically, buoyant complexities are suppressed at mu-g and the reaction-diffusion structure associated with low temperature oxidation reactions, cool flames and auto-ignitions can be studied. Without natural convection, the system is simpler, does not require determination of the effective heat transfer coefficient, and is a testbed for analytic and numerical models that assume pure diffusive transport. In addition, mu-g experiments will provide baseline data that will improve our understanding of the effects of natural convection on Earth.

Pearlman, Howard↗

On spatial instability of a swirling jet - Theory and experiment

The spatial instability of a swirling jet is investigated both experimentally and theoretically. A hydrodynamic stability analysis is applied to an inviscid incompressible top-hat jet, with a swirl distribution of solid-body rotation and free vortex in and outside the vortex core, respectively. Both plane and helical instability modes are examined. It is found that the top-hat jet with swirl distribution of Rankine vortex type is unstable in all the modes studied. The higher the positive helicity, the less spatially unstable the jet behavior; the higher the negative helicity, the more spatially unstable this behavior becomes. A comparison is made between theoretical results and experimental data on a low-intensity swirling jet. The trend of the initial growth of the instability waves in the near field is captured.

Wu, C.↗

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗