SEARCH · Search NASA

Results for “Hydrodynamic Stability Analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 3rd, California State University, Long Beach, CA, January 21-24, 1985, Proceedings

Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.

Source record↗

NASA NTRS1990

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

NASA NTRS1987

Stability of a rigid rotor supported on oil-film journal bearings under dynamic load

Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.

Majumdar, B. C.↗

NASA NTRS1984

The dynamic fission instability and the origin of the Moon

A theory for the formation of the Moon which involves the dynamic fission of a rapidly rotating protoplanet, which might then result in the formation of the Earth and the Moon is discussed. The fission hypothesis was originally based on analytic, linearized models of the growth of asymmetry in homogenous bodies. The fully nonlinear evolution of the dynamic instability in inviscid, compressible bodies was calculated by numerical techniques. It was found that the dynamic instability degenerates into the ejection of a ring of matter with a substantial fraction of the mass, leaving behind a central body with most of the mass. The linearized analytical approach and the numerical approach were used to show that dynamic fission probably does not occur in rocky protoplanets. The numerical calculations are performed with a fully three dimensional hydrodynamical code, which allows the nonlinear, time evolution of the instability to be followed. Sequences of uniformly rotating equilibria were constructed and are used as the initial models for the fission calculations. An initially imposed asymmetry consisting of a 10% binary perturbation in the density was found to disappear on the rotational period time scale. No dynamic instability occurred. This result are verified by including the velocity dissipation terms in the linearized analysis of the stability of a Maclaurin spheroid: the dynamic instability disappears when the simulated viscous dissipation terms are included. It is concluded that any rocky body, even with considerable partial melt or a molten core, should be stable to dynamic fission; any rotational instability that occurs can only result in equatorial mass loss.

Boss, A. P.↗

NASA NTRS1997

The Stability of Radiatively Cooling Jets. 2: Nonlinear Evolution

We use two-dimensional time-dependent hydrodynamical simulations to follow the growth of the Kelvin-Helmholtz (K-H) instability in cooling jets into the nonlinear regime. We focus primarily on asymmetric modes that give rise to transverse displacements of the jet beam. A variety of Mach numbers and two different cooling curves are studied. The growth rates of waves in the linear regime measured from the numerical simulations are in excellent agreement with the predictions of the linear stability analysis presented in the first paper in this series. In the nonlinear regime, the simulations show that asymmetric modes of the K-H instability can affect the structure and evolution of cooling jets in a number of ways. We find that jets in which the growth rate of the sinusoidal surface wave has a maximum at a so-called resonant frequency can be dominated by large-amplitude sinusoidal oscillations near this frequency. Eventually, growth of this wave can disrupt the jet. On the other hand, nonlinear body waves tend to produce low-amplitude wiggles in the shape of the jet but can result in strong shocks in the jet beam. In cooling jets, these shocks can produce dense knots and filaments of cooling gas within the jet. Ripples in the surface of the jet beam caused by both surface and body waves generate oblique shock "spurs" driven into the ambient gas. Our simulations show these shock "spurs" can accelerate ambient gas at large distances from the jet beam to low velocities, which represents a new mechanism by which low-velocity bipolar outflows may be driven by high-velocity jets. Rapid entrainment and acceleration of ambient gas may also occur if the jet is disrupted. For parameters typical of protostellar jets, the frequency at which K-H growth is a maximum (or highest frequency to which the entire jet can respond dynamically) will be associated with perturbations with a period of - 200 yr. Higher frequency (shorter period) perturbations excite waves associated with body modes that produce internal shocks and only small-amplitude wiggles within the jet. The fact that most observed systems show no evidence for large-amplitude sinusoidal oscillation leading to disruption is indicative that the perturbation frequencies are generally large, consistent with the suggestion that pro- tostellar jets arise from the inner regions (r less than 1 AU) of accretion disks.

Stone, James M.↗

NASA NTRS1990

Instability of a supersonic shock free elliptic jet

This paper presents a comparison of the measured and the computed spatial stability properties of an aspect ratio 2 supersonic shock free elliptic jet. The shock free nature of the elliptic jet provides an ideal test of validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with linear hydrodynamic stability theory. Both aerodynamic and acoustic data were measured. The data are used to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. A hybrid numerical scheme is applied to solve the Rayleigh problem governing the inviscid linear spatial stability of the jet. The measured mean velocity profiles are used to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles used in the stability analysis. Computational results are presented for several modes of instability at two jet cross sections. The acoustic measurements show that a varicose instability is the jet's perferred mode of motion. The stability analysis predicts that the Strouhal number varies linearly as a function of axial distance in the jet's initial mixing region, which is in good qualitative agreement with previous measurements.

Baty, Roy S.↗

NASA NTRS1980

Viscous flow drag reduction; Symposium, Dallas, Tex., November 7, 8, 1979, Technical Papers

The symposium focused on laminar boundary layers, boundary layer stability analysis of a natural laminar flow glove on the F-111 TACT aircraft, drag reduction of an oscillating flat plate with an interface film, electromagnetic precipitation and ducting of particles in turbulent boundary layers, large eddy breakup scheme for turbulent viscous drag reduction, blowing and suction, polymer additives, and compliant surfaces. Topics included influence of environment in laminar boundary layer control, generation rate of turbulent patches in the laminar boundary layer of a submersible, drag reduction of small amplitude rigid surface waves, and hydrodynamic drag and surface deformations generated by liquid flows over flexible surfaces.

Hough, G. R.↗

NASA NTRSDate not supplied

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

It has been well-known that the Ingard-Myers impedance condition, while simple to apply, is subject to the hydrodynamic Kelvin-Helmholtz instability due to its use of a vortex sheet in modeling the flow at the liner boundary. Recently, in the development of a time domain boundary element method for acoustic scattering by treated surfaces, it was found that by neglecting a certain second-order spatial derivative term in the Ingard-Myers formulation, the hydrodynamic instability can be avoided. The present paper aims to provide further analysis of this modified condition, hereby referred to as Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC). It will be shown, based on the dispersion relations of linear waves, that the instability intrinsic to the Ingard-Myers condition is eliminated in the proposed new formulation. Quantitative assessments on the accuracy of TIMIBC for scattering of acoustic waves by lined surfaces will be carried out, and its effectiveness will be demonstrated by numerical examples. Specifically, the accuracy is assessed by comparing solutions obtained by the Ingard-Myers condition with that by the proposed TIMIBC where theoretical reflection coefficients at a lined surface are derived for cases of plane and spherical incident waves. It is found that the TIMIBC provides a good approximation to the original Ingard-Myers condition for flows of low to mid subsonic Mach numbers. Time domain implementations of TIMIBC are also discussed and illustrated with a numerical example using a finite difference scheme. As many studies have shown the Ingard-Myers condition to be the correct limit of boundary layer thickness going to zero, the proposed TIMIBC can offer a practical solution for overcoming the intrinsic instability associated with the Ingard-Myers condition.

Time Domain BEM↗

NASA NTRSDate not supplied

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

Time Domain BEM↗

DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗

NASA NTRS1985

The utilization of satellite data and dynamics in understanding and predicting global weather phenomena

A two layer spectral quasi-geostrophic model is used to simulate the effects of topography on the equilibria, the stability, and the long term evaluation of incipient unstable waves. The flow is forced by latitudinally dependent radiational heating. The nature of the form drag instability of high index equilibria is investigated. The proximity of the equilibrium shear to a resonant value is essential for the instability, provided the equilibrium occurs at a slightly stronger shear than resonance. The properties of the steady Hadley and Rossby required for a thermally forced rotating fluid on a sphere are further explained. An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. Also documented is a steady, axisymmetric model of the general circulation developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. Included is further research on cloud street phenomena. Orientation angles and horizontal wavelengths of boundary layer rolls and cloud streets are determined from an analysis of a truncated spectral model of three dimensional shallow moist Boussinesq convection in a shearing environment is further explained. Relatively broadly spaced roll clouds have orientations for which the Fourier component of the roll perpendicular shear is nearly zero, but the second corresponds to narrowly spaced rolls having orientations for which the Fourier coefficients of both the perpendicular and the parallel components of the shear are nearly equal.

Shirer, H. N.↗

NASA NTRS2021

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Boundary layer transition↗

NASA NTRS1991

Equilibria, stability and bifurcations of rotating columns of fluid subjected to planar disturbances

Long gyrostatically rotating drops bonded by surface tension are amenable to conventional bifurcation analysis and newer, computer-aided analytical methods, and therefore are useful prototypes of three-dimensional drops. A study is conducted by setting aside instability to Rayleigh's axisymmetric mode and investigating the effects of translationally symmetric (planar) disturbances. The disadvantage of employing single-coordinate representation of drop shapes close to break-up is brought out. It is shown that a family of symmetric two-lobed shapes bifurcates from the main family of perfectly cylindrical shapes when the rotation rate reaches a critical value, in accord with the linearized hydrodynamic analysis of Hocking.

Benner, R. E., Jr.↗

NASA NTRS2021

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗

NASA NTRS2021

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

Boundary Layer Transition↗

NASA NTRS1991

Three-Dimensional Navier-Stokes Simulation of Space Shuttle Main Propulsion 17-inch Disconnect Valves

A steady incompressible three-dimensional viscous flow analysis has been conducted for the Space Shuttle external tank/orbiter propellant feed line disconnect flapper valves with upstream elbows. The Navier-Stokes code, INS3D, is modified to handle interior obstacles and a simple turbulence model. The flow solver is tested for stability and convergence in the presence of interior flappers. An under-relaxation scheme has been incorporated to improve the solution stability. Important flow characteristics such as secondary flows, recirculation, vortex and wake regions, and separated flows are observed. Computed values for forces, moments, and pressure drop are in satisfactory agreement with water flow test data covering a maximum tube Reynolds number of 3.5 million. The predicted hydrodynamical stability of the flappers correlates well with the measurements.

Kandula, M.↗

DOE OSTI2025

Assessing the performance of MagLIF with 3D MHD simulations

The Magnetized Liner Inertial Fusion (MagLIF) experimental platform at Sandia National Laboratories has realized a number of recent improvements in electrical current and laser preheat coupling but fusion yields are below expectations, based on modeling with clean (that is, simulations with no interfacial instabilities and no mix) two-dimensional (2D) magneto-hydrodynamics (MHD) simulations, by a factor of 5 or more. However, magnetized liner implosions are known to produce helical magneto-Rayleigh–Taylor (MRT) modes and complex stagnation structures that cannot be directly modeled in 2D. This paper presents the results of 3D HYDRA MHD simulations, including helical MRT, showing that degradation from these instabilities can readily reproduce experimental yields, but it is difficult to simultaneously match the yield along with inferred fuel temperature, pressure, and burn history (assessed with x-rays), which are larger in the experiments considered here. Additional analysis methods and future experiments are proposed to help address the discrepancies. The 3D simulations also show, without substantive improvements to stability, changing the applied Bz or preheat is unlikely to increase performance of MagLIF to the same degree as clean 2D simulations. Finally, the first 3D HYDRA simulations with the Hall term show the ab initio production of helical MRT and produce a clear change in stagnation morphology compared to pre-seeded simulations without Hall. In spite of the differences, the simulations still produce comparable fusion performance but with longer, lower power, neutron yield history.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗