A study of longitudinal oscillations of propellant tanks and wave propagations in feed lines. Part II - Wave propagation in elastic pipe filled with incompressible viscous fluid
Fluid wave propagation in liquid propellant feed system
SEARCH · Search NASA
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.
Fluid wave propagation in liquid propellant feed system
The evolution of a disturbance in the presence of a uniform mean velocity gradient is calculated by a power-series solution of the incompressible Navier-Stokes equations. Terms through those in time cubed are retained in the solution. For the initial condition a three-dimensional cosine distribution with two harmonic terms is assumed. The nonlinear interaction of these harmonic terms produces new harmonics which in turn interact. For large velocity gradients the energy of the disturbance can grow with time. The results shed some light on the maintenance or growth of turbulence in a shear flow.
For both the linearized prolate and oblate problems, perturbation series for all orders of eccentricity are obtained. The 'weak separation' method of solving partial-differential-equation boundary-value problems assumes that each linearly independent solution is a sum of products of single-variable functions. The method is applied to a hierarchy of partial differential equations. The weak-separation solutions are constructed by generalizing only the Bessel functions of the homogeneous solution.
A computational method for solving nonlinear problems relating to impact and penetration of a rigid body into a fluid type medium is presented. The numerical techniques, based on the Marker-and-Cell method, gives the pressure and velocity of the flow field. An important feature in this method is that the force and displacement of the rigid body interacting with the fluid during the impact and sinking phases are evaluated from the boundary stresses imposed by the fluid on the rigid body. A sample problem of low velocity penetration of a rigid block into still water is solved by this method. The computed time histories of the acceleration, pressure, and displacement of the block show food agreement with experimental measurements. A sample problem of high velocity impact of a rigid block into soft clay is also presented.
The stability of an electrically conducting fluid subjected to two dimensional disturbance was investigated. A physical system consisting of two parallel infinite vertical plates which are thermally insulated was studied. An external magnetic field of constant strength was applied to normal plates. The fluid was heated from below so that a steady temperature gradient was maintained in the fluid. The governing equations were derived by perturbation technique, and solutions were obtained by a modified Galerkin method. It was found that the presence of the magnetic field increases the stability of the physical system and instability can occur in the form of neutral or oscillatory instability.
Explore the source record for details and available documents.
Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. The resulting models, for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-veocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studied. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SHSV in the vicinity of the low-velocity zone. Near the East Pacific Rise, however, SVSH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of about 50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures A new model for Average Earth is obtained which includes a thin LID. This model extends the fit of a P.R.E.M. type model to shorter period surface waves.
Measurements were taken of the path of a simulated typical tissue segment or 'particle' within a rotating fluid as a function of gravitational strength, fluid rotation rate, particle sedimentation rate, and particle initial position. Parameters were examined within the useful range for tissue culture in the NASA rotating wall culture vessels. The particle moves along a nearly circular path through the fluid (as observed from the rotating reference frame of the fluid) at the same speed as its linear terminal sedimentation speed for the external gravitational field. This gravitationally induced motion causes an increasing deviation of the particle from its original position within the fluid for a decreased rotational rate, for a more rapidly sedimenting particle, and for an increased gravitational strength. Under low gravity conditions (less than 0.1 G), the particle's motion through the fluid and its deviation from its original position become negligible. Under unit gravity conditions, large distortions (greater than 0.25 inch) occur even for particles of slow sedimentation rate (less than 1.0 cm/sec). The particle's motion is nearly independent of the particle's initial position. Comparison with mathematically predicted particle paths show that a significant error in the mathematically predicted path occurs for large particle deviations. This results from a geometric approximation and numerically accumulating error in the mathematical technique.
The Navier-Stokes stress principle is checked in the light of Maxwell's mechanism of friction and in connection herewith the possibility of another theorem is indicated. The Navier-Stokes stress principle is in general predicated upon the conception of the plastic body. Hence the process is a purely phenomenological one, which Newton himself followed with his special theorem for one-dimensional flows. It remained for Maxwell to discover the physical mechanism by which the shear inflow direction is developed: According to it, this shear is only 'fictitious' as it merely represents the substitute for a certain transport on macroscopic motion quantity, as conditioned by Brown's moiecular motion and the diffusion, respectively. It is clear that this mechanism is not bound to the special case of the one-dimensioilal flows, but holds for any flow as expression of the diffusion, by which a fluid differs sharply from a plastic body. If it is remembered, on the other hand, that the cause of the stresses on the plastic body lies in a certain cohesion of the molecules, it appears by no means self evident that this difference in the mechanism of friction between fluid and plastic body should not prevail in the stress principle as well, although it certainly is desirable in any case, at least subsequently, to establish the general theorem in the sense of Maxwell. Actually, a different theorem is suggested which, in contrast to that by Navier-Stokes, has the form of an unsymmetrical matrix. Without anticipating a final decision several reasons are advanced by way of a special flow which seem to affirm this new theorem. To make it clear that the problem involved here still awaits its final solution, is the real purpose behind the present article.
The dynamics of miscible displacements in a cylindrical tube are being investigated experimentally and numerically, with a view to understand the complex processes that occur, for example, in enhanced oil recovery, hydrology, and filtration. We have observed complex shapes of the interface between two liquids that mix with each other when the less viscous liquid is displaced by the more viscous one in a tube. A less viscous fluid that displaces a more viscous fluid is known to propagate in the form of a "finger," and a flight experiment proposed by Maxworthy et al. to investigate the miscible-interface dynamics is currently being developed by NASA. From the current theory of miscible displacements, which was developed for a porous medium satisfying Darcy's law, it can be shown that in the absence of gravity the interface between the fluids is destabilized and thus susceptible to fingering only when a more viscous fluid is displaced by a less viscous one. Therefore, if the interface is initially flat and the more viscous fluid displaces the less viscous fluid, the interface ought to be stable and remain flat. However, numerical simulations by Chen and Meiburg for such displacement in a cylindrical tube show that the interface is unstable and a finger of the more viscous fluid is indeed formed. Preliminary experiments performed at the NASA Glenn Research Center show that not only can fingering occur when the more viscous fluid displaces a less viscous one in a cylindrical tube, but also that under certain conditions the advancing finger achieves a sinuous or snakelike shape. These experiments were performed using silicone oils in a vertical pipette of small diameter. In the initial configuration, the more viscous fluid rested on top of the less viscous one, and the interface was nominally flat. A dye was added to the upper liquid for ease of observation of the interface between the fluids. The flow was initiated by draining the lower fluid from the bottom of the pipette, at speeds less than 0.1 mm/sec.
Fluid-flow simulation over a computer-generated aircraft surface is generated using inviscid and viscous simulations. A fluid-flow mesh of fluid cells is obtained. At least one inviscid fluid property for the fluid cells is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. A set of intersecting fluid cells that intersects the aircraft surface are identified. One surface mesh polygon of the surface mesh is identified for each intersecting fluid cell. A boundary-layer prediction point for each identified surface mesh polygon is determined. At least one boundary-layer fluid property for each boundary-layer prediction point is determined using the at least one inviscid fluid property of the corresponding intersecting fluid cell and a boundary-layer simulation that simulates fluid viscous effects. At least one updated fluid property for at least one fluid cell is determined using the at least one boundary-layer fluid property and the inviscid fluid simulation.
We report on experiments in which a more viscous fluid displaces a less viscous one in a vertical cylindrical tube. These experiments were performed using silicone oils in a vertical pipette of small diameter. The more viscous fluid also had a slightly larger density than the less viscous fluid. In the initial configuration, the fluids were at rest, and the interface was nominally flat. A dye was added to the more viscous fluid for ease of observation of the interface between the fluids. The flow was initiated by pumping the more viscous fluid into the less viscous one. The displacement velocity was such that the Reynolds number was smaller than unity and the Peclet number for mass transfer between the fluids was large compared to unity. For upward displacement of the more viscous fluid from an initially stable configuration, an axisymmetric finger was observed under all conditions. However, a needle-shaped spike was seen to propagate from the main finger in many cases, similar to that observed by Petitjeans and Maxworthy for the displacement of a more viscous fluid by a less viscous one.
Magnetohydrodynamic wakes in a viscous conducting fluid
The case of two incompressible viscous fluids (with different densities) which meet at a planar interface and are subject to an externally imposed oscillating acceleration directed normal to the interface is studied. The linear evolution of perturbations to the basic-state flow is investigated. It is found that the critical Stokes-Reynolds number and the most unstable perturbation wavelengths are insensitive to the degree of density and viscosity differences between the two fluids.
Difference method computations of nonstationary flow of viscous incompressible fluid
Longitudinal wave propagation in elastic pipe filled with incompressible viscous streaming fluid
This report presents experiments of a less viscous liquid displaced by a more viscous, miscible liquid in a vertical tube of small diameter. The liquids used are silicone oils of various viscosities. The more viscous fluid has a higher density than the less viscous fluid. The experiments showed that the interface between the liquids can be unstable and exhibits an asymmetric, sinuous shape for downward displacement. The Reynolds number is negligibly small in these experiments and the Peclet number for mass transfer is O(50)-0(l0(exp 5)). The influence of gravity is characterized by the dimensionless parameter F=g(Delta)(rho)d(sup 2)/(mu)U where g is the acceleration due to gravity, Delta(rho) is the density difference between the fluids, d is the tube diameter, mu is the viscosity of the more viscous fluid, and U is the displacement speed. Visualization and PIV measurement of the unstable linger are reported. Experiments were also conducted to study immiscible displacements in tubes in order to compare with previously obtained results for miscible fluids.
We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.