Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid instabilities”

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 541 records · Page 30

Effects of Inlet Distortion on Aeromechanical Stability of a Forward-Swept High-Speed Fan

Concerns regarding noise, propulsive efficiency, and fuel burn are inspiring aircraft designs wherein the propulsive turbomachines are partially (or fully) embedded within the airframe; such designs present serious concerns with regard to aerodynamic and aeromechanic performance of the compression system in response to inlet distortion. Separately, a forward-swept high-speed fan was developed to address noise concerns of modern podded turbofans; however this fan encounters aeroelastic instability (flutter) as it approaches stall. A three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is applied to analyze and corroborate fan performance with clean inlet flow. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions, is modified and then applied in a computational study to preliminarily assess the effects of inlet distortion on aeroelastic stability of the fan. Computational engineering application and implementation issues are discussed, followed by an investigation into the aeroelastic behavior of the fan with clean and distorted inlets.

Herrick, Gregory P.↗

Extensional Flow of a Polystyrene Boger Fluid Through a 4:1:4 Axisymmetric Contraction/Expansion

The creeping flow of a dilute (0.025 wt%) monodisperse polystyrene/polystyrene Boger fluid through a 4:1:4 axisymmetric contraction/expansion is experimentally observed for a wide range of Deborah numbers. Pressure drop measurements across the orifice plate show a large extra pressure drop that increases monotonically with Deborah number above the value observed for a similar Newtonian fluid at the same flow rate. This enhancement in the dimensionless pressure drop is not associated with the onset of a flow instability, yet it is not predicted by existing steady-state or transient numerical computations with simple dumbbell models. It is conjectured that this extra pressure drop is the result of an additional dissipative contribution to the polymeric stress arising from a stress-conformation hysteresis in the strong non-homogeneous extensional flow near the contraction plane. Such a hysteresis has been independently measured and computed in recent studies of homogeneous transient uniaxial stretching of PS/PS Boger fluids. Flow visualization and velocity field measurements using digital particle image velocimetry (DPIV) show large upstream growth of the corner vortex with increasing Deborah number. At large Deborah numbers, the onset of an elastic instability is observed, first locally as small amplitude fluctuations in the pressure measurements, and then globally as an azimuthal precessing of the upstream corner vortex accompanied by periodic oscillations in the pressure drop across the orifice.

Rothstein, Jonathan P.↗

Final Report for Tech-X Corporations’s contribution to Enhancing Understanding of High Energy Density Plasmas from Wire Array and Solid Liner Implosions Using Fluid Modeling with Kinetic Closures

Recent results from experiments and simulation of magnetically driven pulsed power liners have explored the role of the early-time electrothermal instability in the evolution of the magneto-Rayleigh-Taylor instability. Our focus will be on understanding the development of such instabilities and the potential stabilization mechanisms via electron thermal conduction, viscosity, and large magnetic fields which we expect could play a significant role in supporting the success of the MagLIF program. Experiments have shown that studies of high-energy density plasmas from wire-array implosions require physics modeling that goes well beyond simple models such as ideal magnetohydrodynamics. The goal of this work is to provide increased understanding of these experiments by employing simulations with a multi-fluid extended-MHD model which uses kinetic closures for thermal conduction, resistivity and viscosity. We will use codes easily available to the wider research community, including university students, with a secondary goal of providing the community with well-benchmarked tools capable of advanced modeling of high-energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The proton cyclotron instability and the anisotropy/B inverse correlation

Spacecraft observations in the strongly compressed subsolar magnetosheath show an inverse correlation between the proton temperature anistrophy (T(sub perpendicular to p)/T(sub parallel to p) greater than 1 where perpendicular to and parallel to denote directions perpendicular and parallel to the background magnetic field) and the parallel proton beta (beta(sub parallel to p). This manuscript uses one-dimensional hybrid simulations of the proton cyclotron anisotropy instability in homogeneous electron-proton plasmas to study this correlation which may represent a limited closure relation for fluid theories of anisotropic space plasmas. The emphasis is on driven simulation which increase the temperature anisotropy by periodically reducing the magnetic-field-aligned velocities of the protons. The late-time states from ensembles of both initial value and driven simulations yield very similar expressions for the proton anisotropy/beta(sub parallel to P) inverse correlation, and provide a basis for explaining differences between sheath observations from different spacecraft. The driven simulations also yield expressions for the maximum instability growth rate and the fluctuating field energy as functions of beta (sub parallel to p) and a parameter characterizing the anisotropy driver.

Gary, S. Peter↗

Numerical investigation of the control of the secondary instability process in boundary layers

A numerical model has been developed for investigating boundary layer transition control for a flat plate boundary layer. Active control of a periodically forced boundary layer in an incompressible fluid is studied using surface heating techniques. The spatially evolving boundary layer is simulated. The Navier-Stokes and energy equations are integrated using a fully implicit finite difference/spectral method. Temperature perturbations are introduced locally along finite heater strips to directly attenuate the instability waves in the flow. A feedback control loop is employed in which a downstream sensor is used to monitor wall shear stress fluctuations. Active control of small amplitude two-dimensional and three-dimensional disturbances is numerically simulated. With proper phase control, in-phase reinforcement and out-of-phase attenuation are demonstrated. A receptivity study of the localized temperature perturbations is made. It is shown that narrow heater strips are more receptive in that they maximize the amplitude level of the disturbances in the flow. Active control of the early stages of the fundamental breakdown process is also numerically simulated. Control is achieved with either two-dimensional or three-dimensional control inputs.

Kral, L. D.↗

The fluid mechanics of continuous flow electrophoresis in perspective

Buoyancy alters the flow in continuous flow electrophoresis chambers through the mechanism of hydrodynamic instability and, when the instability is supressed by careful cooling of the chamber boundaries, by restructuring the axial flow. The expanded roles of buoyancy follow upon adapting the size of the chamber and the electric field so as to fractionate certain sorts of cell populations. Scale-up problems, hydrodynamic stability and the altered flow fields are discussed to show how phenomena overlooked in the design and operations of narrow-gap devices take on an overwhelming importance in wide-gap chambers

Saville, D. A.↗

Instabilities in uranium plasma.

The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P equal to about .00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.

Tidman, D. A.↗

Beat, modulational, and decay instabilities of a circularly polarized Alfven wave

A circulary polarized low-frequency electomagnetic pump wave propoagating along an ambient magnetic field is known to be unstable to the growth of several parallel-propagating parametric instabilties. If ion-cyclotron effects are retained in a two-fluid description, the dispersion relation is a sixth-order polynomial. We present a series of new analytical approximations to this dispersion relation. We emphasize new results for the beat instability that occurs as an interaction of the forward prpagating upper sideband with the backward propagating lower sideband. The nature of the beat instabitlity depends on beta = (v(sub sound)/v(sub A)(exp 2) and on the sense of polarization of the pump wave. The beat and decay instabilities can occur together if the pump is left-handed (i.e., ion resonant) and if beta is less than or approximately 1, but they cannot occur together if the pump is right-handed. For a left-handed pump the beat mode is the only instability if beta is greater than or approximately 1. If the pump is right-handed and beta is greater than or approximately 1, then the beat instability exists only when the pump amplitude exceeds a threshold value, and the beat will be the only instability if the pump amplitude is large enough to stabilize the modulational instability. If the pump is left-handed and beta is less than or approximately 1, then the beat mode is stabilized when the pump amplitude becomes sufficiently large. The beat instability primarily produces a forward propagating transverse wave in the upper sideband. Thus if beta is greater than or approximately 1, the instabilities considered here do not produce the backward propagating waves which are thought to affect turbulence and the evolution of cross helicity in the solar wind. New analytical results are presented also for the decay and modulational instabilites when beta is approximately equal to 1.

Hollweg, Joseph V.↗

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗

Magnetic Susceptibility Effects and Lorentz Damping in Diamagnetic Fluids

A great number of crystals (semi-conductor and protein) grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity and g-jitter. Both static and dynamic (rotating or travelling wave) magnetic fields can be used to reduce the effects of convection in materials processing. In semi-conductor melts, due to their relatively high electrical conductivity, the induced Lorentz force can be effectively used to curtail convective effects. In melts/solutions with reduced electrical conductivity, such as aqueous solutions used in solution crystal growth, protein crystal growth and/or model fluid experiments for simulating melt growth, however, the variation of the magnetic susceptibility with temperature and/or concentration can be utilized to better damp fluid convection than the Lorentz force method. This paper presents a comprehensive, comparative numerical study of the relative damping effects using static magnetic fields and gradients in a simple geometry subjected to a thermal gradient. The governing equations are formulated in general terms and then simplified for the numerical calculations. Operational regimes, based on the best damping technique for different melts/solutions are identified based on fluid properties. Comparisons are provided between the numerical results and available results from experiments in surveyed literature.

Ramachandran, Narayanan↗

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

2024 NSF CAREER Research Summary

My current research goals have adjusted from a sole focus on Rayleigh-Taylor (RT) driven Biermann Battery magnetic field growth to RT-driven magnetic reconnection caused by distortion of an imposed magnetic field or self-generated fields. The majority of work from 2024 focused on single-fluid extended-MHD studies of Biermann battery magnetic field generation during the Rayleigh-Taylor instability and its impacts on turbulent growth. R-T driven magnetic field growth can produce notable effects during the deceleration-phase of ICF experiments, though theoretical analysis indicates the impacts will be inconsequential for astrophysical turbulence dominated by high strength external fields. This research pivot provides greater application to a wider range of conditions, which now include astrophysics along with ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of solar wind shock models with tensor plasma pressure for data analysis

The development of solar wind shock models with tensor plasma pressure and the comparison of some of the shock models with the satellite data from Pioneer 6 through Pioneer 9 are reported. Theoretically, difficulties were found in non-turbulent fluid shock models for tensor pressure plasmas. For microscopic shock theories nonlinear growth caused by plasma instabilities was frequently not clearly demonstrated to lead to the formation of a shock. As a result no clear choice for a shock model for the bow shock or interplanetary tensor pressure shocks emerged.

Abraham-Shrauner, B.↗

Hydrodynamic stability of jets produced by mass accreting systems

The existing model for pulsed X-ray emission from the source Hercules X-1 is reviewed. A necessary part of this model is a processing accretion disk which turns the source on and off with 35 day cycle. It is usually assumed that precession of the primary star in this binary system, Hz Hercules, slaves the disk to its precession rate. This model can account for the system behavior in a qualitative manner. Precession of Hz Hercules with 35 day period requires precession of the binary orbit. Pulse arrival times from Herc X-1 have been analyzed for orbital precession. The inclusion of precession does not significantly improve the results obtained assuming a non-precessing orbit. The fluid dynamical stability of extra-galactic jets and the possible consequences of Kelvin-Helmholtz instability at the jet surface external medium interface are considered.

Hardee, P. E.↗

Thermocapillary Flows and Their Stability: Effects of Surface Layers and Contamination

The objective is to study quantitatively thermocapillary flows in two-dimensional and axisymmetric geometries in order to learn the impact of such flows on float-zone refining configurations. The work begins with the analytical study of steady flows in slots and zones to obtain predictions of surface curvature, flow and heat transport. It then examines the instability of such flows. Three new instability mechanisms have been identified and the results show that thermocapillary flows of high Prandtl number fluid become convectively unstable while those of low Prandtl number fluid become hydrodynamically unstable. The work includes numerical simulations of steady, high Marangoni number flows in unit aspect ratio boxes and asymptotic solutions in large boxes. Furthermore, the effect of surface contamination on the steady thermocapillary flows was examined. This fundamental study focused on fluid dynamical systems used to understand events in the flow of melts in float-zone geometries both on Earth and in microgravity environments.

Davis, S. H.↗