Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrodynamic stability”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Model for photothermal ionization and molecular recombination during pulsed ablation of polyethylene

The surface erosion of polyethylene is studied under conditions characteristic of extremely high-rate heating expected in laser-driven and high energy density experiments. A coupled model for photoionization and volumetric photothermal degradation of polyethylene is developed that includes finite rate chemistry and molecular recombination of radicals. First, the model is used to explore the evolution of radicals and hydrocarbons during isochoric heating, for which it is found that polyethylene is not in chemical equilibrium for rates >108K/s. Then, the model is used to explore the cooperativity between photoionization and photothermal ionization during one-dimensional ablation from a pulsed heating source, accounting for coupled energy deposition, thermokinetics, thermochemistry, hydrodynamics, mechanics, and thermal conduction in the ALEGRA multi-physics code. The ablation depth per energy pulse is found to agree with an analytical model for instantaneous energy deposition and absent molecular recombination or photoionization. Otherwise, molecular recombination of radicals reduces the depth of ablation per pulse, owing to thermal stabilization on forming branched molecules. Radical formation through photoionization can be compensated for by recombination, but it accelerates ablation at high enough photoionization rates. Finally, the influence of molecular recombination on the ablation depth is found to decrease for long pulses, owing to decreased recombination reaction rates at lower vapor densities near the ablated surface. These and other findings highlight competition between the processes of scission, recombination, and photoionization during pulsed ablation of polyethylene.

Moore, N. W. (ORCID:0000000272259674)↗
NASA NTRS2001

Fluid Physics Under a Stochastic Acceleration Field

The research summarized in this report has involved a combined theoretical and computational study of fluid flow that results from the random acceleration environment present onboard space orbiters, also known as g-jitter. We have focused on a statistical description of the observed g-jitter, on the flows that such an acceleration field can induce in a number of experimental configurations of interest, and on extending previously developed methodology to boundary layer flows. Narrow band noise has been shown to describe many of the features of acceleration data collected during space missions. The scale of baroclinically induced flows when the driving acceleration is random is not given by the Rayleigh number. Spatially uniform g-jitter induces additional hydrodynamic forces among suspended particles in incompressible fluids. Stochastic modulation of the control parameter shifts the location of the onset of an oscillatory instability. Random vibration of solid boundaries leads to separation of boundary layers. Steady streaming ahead of a modulated solid-melt interface enhances solute transport, and modifies the stability boundaries of a planar front.

Vinals, Jorge↗
NASA NTRS1992

On radial oscillations in viscous accretion discs surrounding neutron stars

Radial oscillations resulting from axisymmetric perturbations in viscous accretion disks surrounding neutron stars in X-ray binary systems have been investigated. Within the framework of the alpha-viscosity model a series of hydrodynamic calculations demonstrates that the oscillations are global for alpha of about 1. On the other hand, for alpha of 0.4 or less, the oscillations are local and confined to the disk boundaries. If viscous stresses acting in the radial direction are included, however, it is found that the disk can be stabilized. The application of such instabilities in accretion disks, without reference to the boundary layer region between the neutron star (or magnetosphere) and the inner edge of the disk, to the phenomenology of quasi-periodic oscillations is brought into question.

Chen, Xingming↗
NASA NTRS1986

Convective effects in float-zone and Czochralski melts

The hydrodynamics of crystal-growth melts is a relatively new research area. Numerical modeling of these processes is necessary. The work discussed herein is in two parts: numerical simulations of the flow in a Czochralski melt, and also of that in a float zone. In addition, for the float-zone case, energy stability theory will be used to determine stability bounds for the onset of oscillatory thermo-capillary flow. Convective effects in crystal-growth melts arise from a variety of mechanisms. Temperature gradients both in the direction of gravity and normal to it give rise to convection due to buoyancy effects. Rotation of the crucible and/or crystal causes a forced convection which may augment or oppose the buoyancy-driven flow. Finally, thermo-capillary forces (due to the variation of surface tension with temperature) drive surface motions which in turn generate convection in the bulk fluid. All of these mechanisms are present in either Czochralski or float-zone growth. The objective of the Czochralski modeling is to develop an accurate numerical simulation of the flow in a Czochralski silicon melt and to investigate the effects of various parameters on the flow properties. Like some earlier investigations, the intent is to simulate the effects of buoyancy, forced and thermo-capillary convection, including unsteady effects. Unlike earlier work, the aim is to include the effects of a variable free surface and freezing interface and, possibly incorporate nonaxisymmetric effects.

Neitzel, G. P.↗
NASA NTRS2002

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.↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
NASA NTRS1980

Instability thresholds for flexible rotors in hydrodynamic bearings

Two types of fixed pad hydrodynamic bearings (multilobe and pressure dam) were considered. Optimum and nonoptimum geometric configurations were tested. The optimum geometric configurations were determined by using a theoretical analysis and then the bearings were constructed for a flexible rotor test rig. It was found that optimizing bearings using this technique produces a 100% or greater increase in rotor stability. It is shown that this increase in rotor stability is carried out in the absence of certain types of instability mechanisms such as aerodynamic crosscoupling. However, the increase in rotor stability should greatly improve rotating machinery performance in the presence of such forces as well.

Allaire, P. E.↗
NASA NTRS2017

Assessing Climate Change Impacts on the Stability of Small Tidal Inlets: Part 2- Data Rich Environments

Climate change (CC) is likely to affect the thousands of bar-built or barrier estuaries (here referred to as Small tidal inlets - STIs) around the world. Any such CC impacts on the stability of STIs, which governs the dynamics of STIs as well as that of the inlet-adjacent coastline, can result in significant socio-economic consequences due to the heavy human utilisation of these systems and their surrounds. This article demonstrates the application of a process based snap-shot modelling approach, using the coastal morphodynamic model Delft3D, to 3 case study sites representing the 3 main STI types; Permanently open, locationally stable inlets (Type 1), Permanently open, alongshore migrating inlets (Type 2) and Seasonally/Intermittently open, locationally stable inlets (Type 3). The 3 case study sites (Negombo lagoon - Type 1, Kalutara lagoon - Type 2, and Maha Oya river - Type 3) are all located along the southwest coast of Sri Lanka. After successful hydrodynamic and morphodynamic model validation at the 3 case study sites, CC impact assessment are undertaken for a high end greenhouse gas emission scenario. Future CC modified wave and riverflow conditions are derived from a regional scale application of spectral wave models (WaveWatch III and SWAN) and catchment scale applications of a hydrologic model (CLSM) respectively, both of which are forced with IPCC Global Climate Model output dynamically downscaled to approximately 50 km resolution over the study area with the stretched grid Conformal Cubic Atmospheric Model CCAM. Results show that while all 3 case study STIs will experience significant CC driven variations in their level of stability, none of them will change Type by the year 2100. Specifically, the level of stability of the Type 1 inlet will decrease from 'Good' to 'Fair to poor' by 2100, while the level of (locational) stability of the Type 2 inlet will also decrease with a doubling of the annual migration distance. Conversely, the stability of the Type 3 inlet will increase, with the time till inlet closure increasing by approximately 75%. The main contributor to the overall CC effect on the stability of all 3 STIs is CC driven variations in wave conditions and resulting changes in longshore sediment transport, not Sea level rise as commonly believed.

IPC↗
DOE OSTI2021

Magnetic-field generation and its effect on ablative Rayleigh–Taylor instability in diffusive ablation fronts

The effects of self-generated magnetic fields on the ablative Rayleigh–Taylor (RT) instability are investigated in the linear regime. The main governing parameters are the Froude number (Fr), which stands for the ratio between ablative convection and acceleration of the target, and the Mach number at the ablation front (Ma), assumed to be small (isobaricity). During the development of the RT instability, magnetic fields are generated due to misalignment between pressure and density gradients (Biermann-battery effect). They accumulate at the section of the ablation front where the Nernst and the plasma velocities cancel each other. The magnetic field modifies the dynamics of the instability through the Righi–Leduc term, which acts as a heat source in the energy equation. It is found that the B fields affect perturbations with short wavelengths up to the most unstable wave in the spectrum. Here, the B field plays a destabilizing role for moderate Froude numbers and becomes stabilizing for large Froude numbers. For plastic ablators, the Fr threshold is found to be Fr = 5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Hydrodynamics of rotating stars and close binary interactions: Compressible ellipsoid models

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

Lai, Dong↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multi-frequency General Relativistic Radiation-hydrodynamics with M 1 Closure

We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++ , which expands the two-moment, M 1 , radiation treatment from gray to multi-frequency transport, including Doppler and gravitational frequency shifts. The solver accommodates either photon (Bose–Einstein) or neutrino (Fermi–Dirac) statistical distribution functions with absorption, emission, and elastic scattering processes. An implicit scheme is implemented to simultaneously solve the primitive inversion problem together with the radiation–matter coupling source terms, providing stability over a broad range of opacities and optical depths where the interaction terms can be stiff. We discuss our formulations and numerical methods, and validate our methods against a wide variety of test problems spanning optically thin to thick regimes in flat, weakly curved, and strongly curved spacetimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Experiments on excitation of Alfvén eigenmodes by alpha-particles with bump-on-tail distribution in JET DTE2 plasmas

Dedicated experiments were performed in JET DTE2 plasmas for obtaining an α-particle bump-on-tail (BOT) distribution aiming at exciting Alfvén eigenmodes (AEs). Neutral beam injection-only heating with modulated power was used so that fusion-born α-particles were the only ions present in the MeV energy range in these DT plasmas. The beam power modulation on a time scale shorter than the α-particle slowing down time was chosen for modulating the α-particle source and thus sustaining a BOT in the α-particle distribution. High-frequency modes in the toroidicity-induced Alfven eigenmode (TAE) frequency range and multiple short-lived modes in a wider frequency range have been detected in these DT discharges with interferometry, soft x-ray cameras, and reflectometry. The modes observed were localised close to the magnetic axis and were not seen in the Mirnov coils. Analysis with the TRANSP and Fokker-Planck FIDIT codes confirms that α-particle distributions with BOT in energy were achieved during some time intervals in these discharges though no clear correlation was found between the times of the high-frequency mode excitation and the BOT time intervals. The combined magneto-hydrodynamic (MHD) and kinetic modelling studies show that the high-frequency mode in the TAE frequency range is best fitted with a TAE of toroidal mode number n= 9. This mode is driven mostly by the on-axis beam ions while the smaller drive due to the pressure gradient of α-particles allows overcoming the marginal stability and exciting the mode. The observed multiple short-lived modes in a wider frequency range are identified as the on-axis kinetic AEs predicted in Rosenbluth and Rutherford.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Modeling of vortex generated sound in solid propellant rocket motors

There is considerable evidence based on both full scale firings and cold flow simulations that hydrodynamically unstable shear flows in solid propellant rocket motors can lead to acoustic pressure fluctuations of significant amplitude. Although a comprehensive theoretical understanding of this problem does not yet exist, procedures were explored for generating useful analytical models describing the vortex shedding phenomenon and the mechanisms of coupling to the acoustic field in a rocket combustion chamber. Since combustion stability prediction procedures cannot be successful without incorporation of all acoustic gains and losses, it is clear that a vortex driving model comparable in quality to the analytical models currently employed to represent linear combustion instability must be formulated.

Flandro, G. A.↗
NASA NTRS1995

Review of Combustion-acoustic Instabilities

Combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system. The problem of combustion-acoustic instability is a concern in many devices for various reasons, as each device may have a unique mechanism causing unsteady heat release rates and many have unique boundary conditions. To accurately predict and quantify combustion-acoustic stabilities, the unsteady heat release rate and boundary conditions need to be accurately determined. The present review brings together work performed on a variety of practical combustion devices. Many theoretical and experimental investigations of the unsteady heat release rate have been performed, some based on perturbations in the fuel delivery system particularly for rocket instabilities, while others are based on hydrodynamic processes as in ramjet dump combustors. The boundary conditions for rocket engines have been analyzed and measured extensively. However, less work has been done to measure acoustic boundary conditions in many other combustion systems.

Oyediran, Ayo↗
DOE OSTI2023

Overcoming Barriers Associated with Oral Delivery of Differently Sized Fluorescent Core‐Shell Silica Nanoparticles

Oral delivery, while a highly desirable form of nanoparticle-drug administration, is limited by challenges associated with overcoming several biological barriers. Here, in this study, the authors study how fluorescent and poly(ethylene glycol)-coated (PEGylated) core-shell silica nanoparticles sized 5 to 50 nm interact with major barriers including intestinal mucus, intestinal epithelium, and stomach acid. From imaging fluorescence correlation spectroscopy studies using quasi-total internal reflection fluorescence microscopy, diffusion of nanoparticles through highly scattering mucus is progressively hindered above a critical hydrodynamic size around 20 nm. By studying Caco-2 cell monolayers mimicking the intestinal epithelia, it is observed that ultrasmall nanoparticles below 10 nm diameter (Cornell prime dots, [C’ dots]) show permeabilities correlated with high absorption in humans from primarily enhanced passive passage through tight junctions. Particles above 20 nm diameter exclusively show active transport through cells. After establishing C’ dot stability in artificial gastric juice, in vivo oral gavage experiments in mice demonstrate successful passage through the body followed by renal clearance without protein corona formation. Results suggest C’ dots as viable candidates for oral administration to patients with a proven pathway towards clinical translation and may generate renewed interest in examining silica as a food additive and its effects on nutrition and health.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analysis into Asymptotic Convergence to Full Nonlinear Solutions and Exploration of the Implication of Numerical Operator Mutation of Differential Systems

A robust, sufficiently accurate and practical hydrodynamic simulation toolset is required as a key component of the modeling and simulation of air-gap electrostatic discharge events. This work was performed to complement these ongoing efforts. In particular, hydrodynamic simulations must be vetted to ensure they are robust and sufficiently accurate over relevant characteristic scales. Verification models were generated in order to cultivate the technical knowledge and expertise needed to properly create, implement and execute numerical simulations. Furthermore, this effort was utilized extensively to educate students on the mathematical and numerical principles underlying hydrodynamic simulations. This education opportunity, provided in a holistic and rigorous manner, has greatly benefited developing scientists and engineers with the necessary understandings and toolsets required to excel at accomplishing the task at hand, and, more generally, it has enabled them to generate key programmatic deliverables. This report articulates several subtilties; specifically, how perturbations, nonlinear behavior, and dissipative mechanisms influence numerical stability, how to properly structure mathematical and numerical solutions, and how to properly generate error estimation/assignment. A more rigorous discussion of the consequences of such topics can be found in the body of this report in Chapters 2 and 3 with qualitative findings discussed in Chapter 4.

97 MATHEMATICS AND COMPUTING↗