Poincare's hydrodynamic analogy in celestial mechanics
Poincare hydrodynamic analogy in celestial mechanics, relating differential equations for dynamic systems with two degrees of freedom and two and three dimensional flow
Szebehely, V.↗
SEARCH · Search NASA
Results for “HYDRODYNAMIC EQUATION”
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.
NASTRAN solutions of problems described by simultaneous parabolic differential equations
NASTRAN solution techniques are shown for a numerical analysis of a class of coupled vector flow processes described by simultaneous parabolic differential equations. To define one physical problem type where equations of this form arise, the differential equations describing the coupled transfers of heat and mass in mechanical equilibrium with negligible mass average velocity are presented and discussed. Also shown are the equations describing seepage when both electrokinetic and hydrodynamic forces occur. Based on a variational statement of the general problem type, the concepts of scalar transfer elements and parallel element systems are introduced. It is shown that adoptation of these concepts allows the direct use of NASTRAN's existing Laplace type elements for uncoupled flow (the heat transfer elements) for treating multicomponent coupled transfer. Sample problems are included which demonstrate the application of these techniques for both steady-state and transient problems.
Mason, J. B.↗
Rapid expansion of polytropes
An extremely massive unevolved star may be idealized as an n = 3 polytropic sphere supported by radiation pressure. Such a polytrope is subject to explosive, supersonic expansion, whose later stages can be described by similarity solution of the dynamical equations. This picture is confirmed by numerical hydrodynamic simulation, which shows that an n = 3 polytrope, initially at rest near equilibrium, undergoes an explosive expansion which approaches the similarity solution at large times.
Barnes, A.↗
Hydrodynamics of binary coalescence. 2: Polytropes with gamma = 5/3
We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.
Rasio, Frederic A.↗
Nonlinearity Effects of Lateral Density Diffusion Coefficient on Gain-Guided VCSEL Performance
Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. In this paper, we apply a hydrodynamic model developed from the semiconductor Bloch equations to numerically study the effects of nonlinearity in the diffusion coefficient on single mode operation and direct modulation of a gain-guided InGaAs/GaAs multiple quantum well laser, operating not too far from threshold. We found that a small diffusion coefficient is advantageous for lowering the threshold current and increasing the modulation bandwidth. Most importantly, the effects of nonlinearity in the coefficient can be approximately reproduced by replacing the coefficient with an effective constant diffusion coefficient, which corresponds roughly to the half height density of the density distribution.
Li, Jian-Zhong↗
Protoplanetary Nebula Evolution using the Beta Viscosity Model
The evolutionary dynamics of a protoplanetary disk is an important component of the planet formation process. In particular, the dynamic and thermodynamic field plays a critical role in chemical evolution, the migration of dust particles in the nebula, and the radial transport of meteoritic components. The dynamic evolution is investigated using analytical solutions of the surface density transport equations using a turbulence model based on hydrodynamic generation of turbulence. It captures the major properties of the disk including region of separation between radial inflow and-outflow and the evolution of the central plane temperature. The analytical formulas are compared with available numerical solutions based on the alpha viscosity model. The beta viscosity model, heretofore used for steady-state disks, is shown to be a useful approximation for unsteady problems.
Davis, Sanford S.↗
The steady global corona
The formation of the steady coronal structure which consists of coronal streamers and holes is seen in terms of a model employing numerical solutions of time-dependent, dissipationless, magnetohydrodynamic equations of motion applicable to the meridional plane. A coronal streamer consists of closed magnetic field lines near the solar surface with overlying and adjacent open field lines, and the atmosphere, which is stationary in the closed region, flows outward in the open region, or coronal hole. The steady coronal structure is obtained by starting the numerical calculation with a state comprising a polytropic, hydrodynamical solution to the steady-state radial equation of motion coupled with a dipole magnetic field. Global coronal configurations are calculated for values of the plasma beta which vary from 0.1 to 100. It is found that the height and lateral extension of the closed region are only weakly dependent on the data.
Steinolfson, R. S.↗
Complete Galilean-Invariant Lattice BGK Models for the Navier-Stokes Equation
Galilean invariance has been an important issue in lattice-based hydrodynamics models. Previous models concentrated on the nonlinear advection term. In this paper, we take into account the nonlinear response effect in a systematic way. Using the Chapman-Enskog expansion up to second order, complete Galilean invariant lattice BGK models in one dimension (theta = 3) and two dimensions (theta = 1) for the Navier-Stokes equation have been obtained.
Qian, Yue-Hong↗
Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines
Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.
Singh, Nagendra↗
Face-seal lubrication. 2: Theory of response to angular misalignement
A theoretical analysis was made of a hypothetical seal operating mode. The hypothetical seal model provides for three degrees of primary ring motion and includes the force and moments induced by primary ring response to seat angular misalignment. This ring response causes a relative angular misalignment between the faces of the primary seal. Hydrodynamic pressure generation is produced by this misalignment. The analysis is based on the Reynolds equation in short bearing form and on a balance of forces and moments that arise from hydrodynamic and secondary seal friction effects. A closed form solution was obtained that can be solved for film thickness and relative angular misalignment.
Ludwig, L. P.↗
The oscillations of superfluid neutron stars
The effects of superfluid hydrodynamics on the oscillations of neutron stars are investigated. The equations describing the small-amplitude pulsations of a neutron star's neutron-proton superfluid mixture are reduced to a system of three second-order equations for three scalar potentials. A variational principle is developed from which the frequencies of the modes of these superfluid oscillations may be estimated. These pulsations are studied by finding analytical solutions to the equations for simple uniform models of nonrotating neutron stars and numerical solutions for realistic models (including rotation). These solutions reveal that the lowest frequency modes are almost indistinguishable from those based on ordinary-fluid hydrodynamics. The analytical solutions also reveal the existence of a new set of modes having no ordinary-fluid counterpart.
Lindblom, Lee↗
Determination of boundaries of cloudiness
Cloud boundary determinations based upon adiabatic flow
CLOUD COVER↗
Dynamic simulation of coronal mass ejections
A model is developed for the formation and propagation through the lower corona of the loop-like coronal transients in which mass is ejected from near the solar surface to the outer corona. It is assumed that the initial state for the transient is a coronal streamer. The initial state for the streamer is a polytropic, hydrodynamic solution to the steady-state radial equation of motion coupled with a force-free dipole magnetic field. The numerical solution of the complete time-dependent equations then gradually approaches a stationary coronal streamer configuration. The streamer configuration becomes the initial state for the coronal transient. The streamer and transient simulations are performed completely independent of each other. The transient is created by a sudden increase in the pressure at the base of the closed-field region in the streamer configuration. Both coronal streamers and coronal transients are calculated for values of the plasma beta (the ratio of thermal to magnetic pressure) varying from 0.1 to 100.
Steinolfson, R. S.↗
An explicit algorithm for a fluid approach to nonlinear optics propagation using splitting and rezoning techniques
The paper presents a three-dimensional analysis of the nonlinear light matter interaction in a hydrodynamic context. It is reported that the resulting equations are a generalization of the Navier-Stokes equations subjected to an internal potential which depends solely upon the fluid density. In addition, three numerical approaches are presented to solve the governing equations using an extension of McCormack predict-corrector scheme. These are a uniform grid, a dynamic rezoned grid, and a splitting technique. It is concluded that the use of adaptive mapping and splitting techniques with McCormack two-level predictor-corrector scheme results in an efficient and reliable code whose storage requirements are modest compared with other second order methods of equal accuracy.
Mattar, F. P.↗
Soliton turbulence
Theoretical and numerical works in atmospheric turbulence have used the Navier-Stokes fluid equations exclusively for describing large-scale motions. Controversy over the existence of an average temperature gradient for the very large eddies in the atmosphere suggested that a new theoretical basis for describing large-scale turbulence was necessary. A new soliton formalism as a fluid analogue that generalizes the Schrodinger equation and the Zakharov equations has been developed. This formalism, processing all the nonlinearities including those from modulation provided by the density fluctuations and from convection due to the emission of finite sound waves by velocity fluctuations, treats large-scale turbulence as coalescing and colliding solitons. The new soliton system describes large-scale instabilities more explicitly than the Navier-Stokes system because it has a nonlinearity of the gradient type, while the Navier-Stokes has a nonlinearity of the non-gradient type. The forced Schrodinger equation for strong fluctuations describes the micro-hydrodynamical state of soliton turbulence and is valid for large-scale turbulence in fluids and plasmas where internal waves can interact with velocity fluctuations.
Tchen, C. M.↗
Radiative shocks and hydrogen molecules in pregalactic gas - The effects of postshock radiation
Previous calculations of radiative shocks in pregalactic gas are generalized and substantially improved. The hydrodynamical conservations are solved, along with the rate equations for nonequilibrium ionization, recombination, and molecule formation, and the equation of radiative transfer, for steady-state shocks in a gas of primordial compositions. This calculation self-consistently includes the effects of the diffuse postshock emission as well as a possible external radiation flux on the postshock flow and the preshock ionization levels. A previous result is confirmed that the shocked gas cools faster than it can recombine and, as a result, is able to form an H2 concentration as high as 0.003 via the formation of H(-) and H2(+) intermediaries due to the enhanced nonequilibrium ionization fraction at 10,000 K. It is shown that, as the level of the external radiation flux is increased, the formation of and cooling by H2 molecules can be inhibited and delayed. The detailed radiation spectra and ionizing photon number fluxes emergent from these shocks are presented.
Kang, Hyesung↗
A multiscale numerical study into the cascade of kinetic energy leading to severe local storms
The cascade of kinetic energy from macro- through mesoscales is studied on the basis of a nested grid system used to solve a set of nonlinear differential equations. The kinetic energy cascade and the concentration of vorticity through the hydrodynamic spectrum provide a means for predicting the location and intensity of severe weather from large-scale data sets. A mechanism described by the surface pressure tendency equation proves to be important in explaining how initial middle-tropospheric mass-momentum imbalances alter the low-level pressure field.
Paine, D. A.↗
Hydrodynamics of Spherical Flows and Geometry of Premixed Flames near the Stagnation Point of Axisymmetric Viscous Counterflows
Counterflow premixed flames play a significant role in the modeling of laminar flames. This is in part motivated by the fact that stretched premixed flames simulate local flamelet dynamics within turbulent premixed flames. In the present study, the modified form of the Navier-Stokes equation for reactive fields introduced earlier is employed to investigate the hydrodynamics of spherical flows embedded within counterflows. The geometry of premixed flames near the stagnation point is also determined. The predictions are in favorable agreement with the experimental observations and prior numerical studies.
Sohrab, Siavash H.↗