Search NASA⌕ Search

SEARCH · Search NASA

Results for “Body Force”

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 307 records · Page 17

Ducted fan acoustic radiation including the effects of nonuniform mean flow and acoustic treatment

Forward and aft acoustic propagation and radiation from a ducted fan is modeled using a finite element discretization of the acoustic field equations. The fan noise source is introduced as equivalent body forces representing distributed blade loading. The flow in and around the nacelle is assumed to be nonuniform, reflecting the effects of forward flight and flow into the inlet. Refraction due to the fan exit jet shear layer is not represented. Acoustic treatment on the inlet and exhaust duct surfaces provides a mechanism for attenuation. In a region enclosing the fan a pressure formulation is used with the assumption of locally uniform flow. Away from the fan a velocity potential formulation is used and the flow is assumed nonuniform but irrotational. A procedure is developed for matching the two regions by making use of local duct modal amplitudes as transition state variables and determining the amplitudes by enforcing natural boundary conditions at the interface between adjacent regions in which pressure and velocity potential are used. Simple models of rotor alone and rotor/exit guide vane generated noise are used to demonstrate the calculation of the radiated acoustic field and to show the effect of acoustic treatment. The model has been used to assess the success of four techniques for acoustic lining optimization in reducing far field noise.

Eversman, Walter↗

Transient studies of G-induced capillary dryout and rewet

A transient, one-dimensional numerical code is developed to model the liquid motion in an axial groove with square cross section. Axial variation in liquid level, shear stress and heat transfer between the groove wall and the liquid, evaporation and transient body forces are accounted for in the model. Dryout and rewet of the groove are allowed; the front location is determined numerically using conservation of mass and linear extrapolation. Several numerical test results are presented and discussed.

Reagan, M. K.↗

Thermoviscoplastic analysis of fibrous periodic composites by the use of triangular subvolumes

The non-linear viscoplastic behavior of fibrous periodic composites is analyzed by discretizing the unit cell into triangular subvolumes. A set of these subvolumes can be configured by the analyst to construct a representation for the unit cell of a periodic composite. In each step of the loading history the total strain increment at any point is governed by an integral equation which applies to the entire composite. A Fourier series approximation allows the incremental stresses and strains to be determined within a unit cell of the periodic lattice. The non-linearity arising from the viscoplastic behavior of the constituent materials comprising the composite is treated as a fictitious body force in the governing integral equation. Specific numerical examples showing the stress distributions in the unit cell of a fibrous tungsten/copper metal-matrix composite under viscoplastic loading conditions are given. The stress distribution resulting in the unit cell when the composite material is subjected to an overall transverse stress loading history perpendicular to the fibers is found to be highly heterogeneous, and typical homogenization techniques based on treating the stress and strain distributions within the constituent phases as homogeneous result in large errors under inelastic loading conditions.

Walker, Kevin P.↗

The effect of wave-particle interactions on the polar winds O(+)

The escape of the polar wind plasma is an important element in the ionosphere-magnetosphere coupling. Both theory and observations indicate that the wave-particle interactions (WPI) play a significant role in the dynamics of ion outflow along open geomagnetic field lines. A Monte Carlo simulation was developed in order to include the effect of the WPI in addition to the factors that are traditionally included in the 'classical' polar wind (i.e. gravity, electrostatic field, and divergence of geomagnetic field lines). The ion distribution function (f(sub j)), as well as the profiles of its moments (density, drift velocity, temperature, etc.) were found for different levels of WPI, that is, for different values of the normalized diffusion rate in the velocity space (D(tilde) (sub j perpendicular). Although the model included O(=), H(+) and eletrons, we presented only the results related to the O(+) ion. We found that (1) both the density and drift velocity of O(+) increased with the WPI strength, and consequently, the O(+) escape flux was enhanced by a factor of up to 10(exp 5), (2) The O(+) ions could be energized up to a few electron volts; (3) for moderate and high levels of WPI D(tilde) (sub perpendicular) (O(+) greater than (tilde) 1, the distribution function f(O(+)) displayed very pronounced conic features at altitudes around 3 R(sub e). Finally, the interplay between the downward body force, the upward mirror force, and the perpendicular heating resulted in the formation of the 'pressure cooker' effect. This phenomena explained some interesting features of our solution, such as, the peak in the O(=) temperature, and the formation of 'ears' and conics for f(O(+)) around 2.5 R(sub e).

Barakat, A. R.↗

A spectral formalism for computing three-dimensional deformations due to surface loads. 1: Theory

We outline a complete spectral formalism for computing high spatial resolution three-dimensional deformations arising from the surface mass loading of a spherically symmetric planet. The main advantages of the formalism are that all surface mass loads are always described using a consistent mathematical representation and that calculations of deformation fields for various spatial resolutions can be performed by simpley altering the spherical harmonic degree truncation level of the procedure. The latter may be important when incorporating improved observational constraints on a particular surface mass load, when considering potential errors in the computed field associated with mass loading having a spatial scale unresolved by the observational constraints, or when treating a number of global surface mass loads constrained with different spatial resolutions. The advantages do not extend to traditional 'Green's function' approaches which involve surface element discretizations of the global mass loads. Another advantage of the spectral formalism, over the Green's function approach, is that a posteriori analyses of the computed deformation fields are easily performed. In developing the spectral formalism, we consider specific cases where the Earth's mantle is assumed to respond as an elastic, slightly anelastic, or linear viscoelastic medium. In the case of an elastic or slightly anelastic mantle rheology the spectral response equations incorporate frequency dependent Love numbers. The formalism can therefore be used, for example, to compute the potentially resonant deformational response associated with the free core nutation and Chandler wobble eigenfunctions. For completeness, the spectral response equations include both body forces, as arise from the gravitational attraction of the Sun and the Moon, and surface mass loads. In either case, and for both elastic and anelastic mantle rheologies, we outline a pseudo-spectral technique for computing the ocean adjustment associated with the total gravitational perturbation induced by the external forcing. Three-dimensional deformations computed using the usual Love number approach are generally referenced to an origin at the center of mass of the undeformed planet. We derive a spectral technique for transforming the results to an origin located at the center of mass of the deformed planet.

Mitrovica, J. X.↗

Marangoni Effects in the Boiling of Binary Fluid Mixtures

Results of very recent experimental studies indicate that during nucleate boiling in some binary mixture, Marangoni effects augment the gravity driven flow of liquid towards the heated surface. With gravity present, it is impossible to separate the two effects. The reduced gravity environment gives an unique opportunity to explore th role of Marangoni effects on the boiling mechanisms free of gravitational body forces that obscure the role of such effects. However, recent experimental results suggest that under reduced gravity conditions, Marangoni effects is the dominant mechanism of vapor-liquid exchange at the surface for some binary mixture. To further explore such effects, experiments have been conducted with water/2-propanol mixtures at three different concentrations under normal gravity with different orientations of the heater surface and under reduce gravity aboard the DC-9 aircraft at NASA Lewis Research Center. The system pressure was sub atmospheric (approx. 8 kP at 1g(n)) and the bulk liquid temperature varied from low subcooling to near saturation. The molar concentrations of 2-propanol tested were 0.015, 0.025, and 0.1. Boiling curves were obtained both for high gravity (approx. 2g(n)) and reduce gravity (approx. 0.01g(n)). For each concentration of 2-propanol, the critical heat flux has been determined in the flight experiments only for reduced gravity conditions. Comparison of boiling curves and CHF obtained under l-g(n) an reduced gravity indicates that boiling mechanism in this mixtures is nearly independent of gravity. The results also indicate that the Marangoni mechanism is strong enough in these mixtures to sustain the boiling under reduced gravity conditions.

Ahmed, Sayeed↗

Microgravity Particle Dynamics

This research seeks to identify the experiment design parameters for future flight experiments to better resolve the effects of thermal and velocity gradients on gas-solid flows. By exploiting the reduced body forces and minimized thermal convection current of reduced gravity experiments, features of gas-solid flow normally masked by gravitationally induced effects can be studied using flow regimes unattainable under unigravity. This paper assesses the physical scales of velocity, length, time, thermal gradient magnitude, and velocity gradient magnitude likely to be involved in laminar gas-solid multiphase flight experiments for 1-100 micro-m particles.

Clark, Ivan O.↗

On the Boundary Conditions at an Oscillating Contact Line: A Physical/Numerical Experimental Program

We will pursue an improved physical understanding and mathematical model for the boundary condition at an oscillating contact line at high Reynolds number. We expect that the body force is locally unimportant for earth-based systems, and that the local behavior may dominate the mechanics of partially-filled reservoirs in the microgravity environment. One important space-based application for this contact-line study is for Faraday-waves. Oscillations in the direction of gravity (or acceleration) can dominate the fluid motion during take-off and reentry with large steady-state accelerations and in orbit, where fluctuations on the order of 10(exp -4)g occur about a zero mean. Our experience with Faraday waves has shown them to be 'cleaner' than those produced by vertical or horizontal oscillation of walls. They are easier to model analytically or computationally, and they do not have strong vortex formation at the bottom of the plate. Hence many, if not most, of the experiments will be performed in this manner. The importance of contact lines in the microgravity environment is well established. We will compare high resolution measurements of the velocity field (lO micro-m resolution) using particle-tracking and particle-image velocimetry as the fluid/fluid interface is approached from the lower fluid. The spatial gradients in the deviation provide additional means to determine an improved boundary condition and a measure of the slip region. Dissipation, the size of the eddy near the contact line, and hysteresis will be measured and compare to linear and nonlinear models of viscous and irrotational but dissipative models.

Perlin, Marc↗

Waves in Radial Gravity Using Magnetic Fluid

We are beginning laboratory experiments using magnetically active ferrofluids to study surface waves in novel geometries. Terrestrial gravity is eliminated from the dynamics, and the magnetic body force felt by ferrofluid in the presence of a magnetic field gradient is used to create a geopotential field which is a section of or an entire sphere or cylinder. New optical, electromagnetic and ultrasonic diagnostic techniques are under development to initially study capillary-gravity wave propagation and interaction in such geometries.

Ohlsen, Daniel R.↗

The numerical simulation of multistage turbomachinery flows

The effect of the unsteady flow field in a multistage compressor on the time-averaged performance was assessed. The energy transport by the unsteady deterministic flow field was taken into account. The magnitude of the body force resulting from the aerodynamic loading on a blade row was compared to the gradient of the stress tensor associated with the unsteady time-resolved flow field generated by the blade row. The magnitude of the work performed by these forces was compared to the divergence of the energy correlations produced by the unsteady time-resolved flow field. The stress tensor and the energy correlations are non-negligible in the end wall regions. The results suggest that the turbulence is the primary source of flow mixing in the midspan region.

Adamczyk, John J.↗

Accuracy of the Generalized Self-Consistent Method in Modelling the Elastic Behaviour of Periodic Composites

Local stress and strain fields in the unit cell of an infinite, two-dimensional, periodic fibrous lattice have been determined by an integral equation approach. The effect of the fibres is assimilated to an infinite two-dimensional array of fictitious body forces in the matrix constituent phase of the unit cell. By subtracting a volume averaged strain polarization term from the integral equation we effectively embed a finite number of unit cells in a homogenized medium in which the overall stress and strain correspond to the volume averaged stress and strain of the constrained unit cell. This paper demonstrates that the zeroth term in the governing integral equation expansion, which embeds one unit cell in the homogenized medium, corresponds to the generalized self-consistent approximation. By comparing the zeroth term approximation with higher order approximations to the integral equation summation, both the accuracy of the generalized self-consistent composite model and the rate of convergence of the integral summation can be assessed. Two example composites are studied. For a tungsten/copper elastic fibrous composite the generalized self-consistent model is shown to provide accurate, effective, elastic moduli and local field representations. The local elastic transverse stress field within the representative volume element of the generalized self-consistent method is shown to be in error by much larger amounts for a composite with periodically distributed voids, but homogenization leads to a cancelling of errors, and the effective transverse Young's modulus of the voided composite is shown to be in error by only 23% at a void volume fraction of 75%.

Walker, Kevin P.↗

A Near-Wall Reynolds-Stress Closure Without Wall Normals

Turbulent wall-bounded complex flows are commonly encountered in engineering practice and are of considerable interest in a variety of industrial applications. The presence of a wall significantly affects turbulence characteristics. In addition to the wall effects, turbulent wall-bounded flows become more complicated by the presence of additional body forces (e.g. centrifugal force and Coriolis force) and complex geometry. Most near-wall Reynolds stress models are developed from a high-Reynolds-number model which assumes turbulence is homogenous (or quasi-homogenous). Near-wall modifications are proposed to include wall effects in near-wall regions. In this process, wall normals are introduced. Good predictions could be obtained by Reynolds stress models with wall normals. However, ambiguity arises when the models are applied in flows with multiple walls. Many models have been proposed to model turbulent flows. Among them, Reynolds stress models, in which turbulent stresses are obtained by solving the Reynolds stress transport equations, have been proved to be the most successful ones. To apply the Reynolds stress models to wall-bounded flows, near-wall corrections accounting for the wall effects are needed, and the resulting models are called near-wall Reynolds stress models. In most of the existing near-wall models, the near-wall corrections invoke wall normals. These wall-dependent near-wall models are difficult to implement for turbulent flows with complex geometry and may give inaccurate predictions due to the ambiguity of wall normals at corners connecting multiple walls. The objective of this study is to develop a more general and flexible near-wall Reynolds stress model without using any wall-dependent variable for wall-bounded turbulent flows. With the aid of near-wall asymptotic analysis and results of direct numerical simulation, a new near-wall Reynolds stress model (NNWRS) is formulated based on Speziale et al.'s high-Reynolds-stress model with wall-independent near-wall corrections. Moreover, only one damping function is used for flows with a wide range of Reynolds numbers to ensure that the near-wall modifications diminish away from the walls.

Yuan, S. P.↗

Boundary Layer Flow Control with a One Atmosphere Uniform Glow Discharge Surface Plasma

Low speed wind tunnel data have been acquired for planar panels covered by a uniform, glow-discharge surface plasma in atmospheric pressure air known as the One Atmosphere Uniform Glow Discharge Plasma (OAUGDP). Streamwise and spanwise arrays of flush, plasma-generating surface electrodes have been studied in laminar, transitional, and fully turbulent boundary layer flow. Plasma between symmetric streamwise electrode strips caused large increases in panel drag, whereas asymmetric spanwise electrode configurations produced a significant thrust. Smoke wire flow visualization and mean velocity diagnostics show the primary cause of the phenomena to be a combination of mass transport and vortical structures induced by strong paraelectric ElectroHydroDynamic (EHD) body forces on the flow.

J Reece Roth↗

A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks

This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Majumdar, Alok↗

Microgravity Experiments and Numerical Modeling of Rotating Buoyant Convection in a Spherical Shell with Latitudinal Thermal Gradients

The results of numerical model calculations are compared with space laboratory experiments for rotating, thermally driven flow in a hemispherical shell, and various flow regimes are described. A spherically symmetric body force, analogous to gravity, is imposed in the radial direction through the use of a dielectric fluid and an electrostatic potential difference across the gap. The spherical boundaries are maintained at constant temperature profiles (with the pole and the inner sphere being the warmer surfaces), and the equatorial wall is an insulator. Typical parameter combinations result in highly nonlinear, but laminar, flow. For weak enough buoyant forcing, the flow is axisymmetric. It consists of a single meridional cell, rising in warm latitudes (in this case, near the pole) and sinking in cool latitudes, with prograde flow in the equatorial region and near the inner hemisphere for other latitudes, and retrograde motion elsewhere. For fast rotation, the first transition due to the instability of this simple flow is also axisymmetric, consisting of rings of convection in the polar region, which propagate poleward. The first nonaxisymmetric convection occurs at Rayleigh numbers which increase with rotation rate. The form of the convection near the transition also depends upon the rate of rotation. Selected flow patterns near the transition as well as those beyond it are studied numerically. For those cases where there exist laboratory experiments with which to compare, the numerical and experimental results agree very well.

Leslie, Fred W.↗

On the Induced Flow of an Electrically Conducting Liquid in a Rectangular Duct by Electric and Magnetic Fields of Finite Extent

Reported here are the results of a systematic study of a model of the direct-current electromagnetic pump. Of particular interest is the motion imparted to the electrically conducting fluid in the rectangular duct by the body forces that result from applied electric and magnetic fields. The purpose of the investigation is to associate the observed fluid motion with the characteristics of the electric and magnetic fields which cause them. The experiments were carried out with electromagnetic fields that moved a stream of copper sulphate solution through a clear plastic channel. Ink filaments injected into the stream ahead of the region where the fields were applied identify the motion of the fluid elements as they passed through the test channel. Several magnetic field configurations were employed with a two-dimensional electric current distribution in order to study and identify the magnitude of some of the effects on the fluid motion brought about by nonuniformities in the electromagnetic fields. A theoretical analysis was used to guide and evaluate the identification of the several fluid motions observed. The agreement of the experimental data with the theoretical predictions is satisfactory. It is found that sizable variations in the velocity profile and pressure head of the output stream are produced by the shape of the electric and magnetic fields.

Rossow, Vernon J.↗

Interfacing a General Purpose Fluid Network Flow Program with the SINDA/G Thermal Analysis Program

A general purpose, one dimensional fluid flow code is currently being interfaced with the thermal analysis program Systems Improved Numerical Differencing Analyzer/Gaski (SINDA/G). The flow code, Generalized Fluid System Simulation Program (GFSSP), is capable of analyzing steady state and transient flow in a complex network. The flow code is capable of modeling several physical phenomena including compressibility effects, phase changes, body forces (such as gravity and centrifugal) and mixture thermodynamics for multiple species. The addition of GFSSP to SINDA/G provides a significant improvement in convective heat transfer modeling for SINDA/G. The interface development is conducted in multiple phases. This paper describes the first phase of the interface which allows for steady and quasi-steady (unsteady solid, steady fluid) conjugate heat transfer modeling.

Schallhorn, Paul↗

Constitutive and Stability Behavior of Granular Materials Under Very Low Effective Stresses

The study of the behavior of granular materials such as sand is of great importance to the understanding of the response of foundations resting on surface and near-surface soil deposits. A clear understanding of the behavior of such materials can also provide insights into other related issues where low effective stresses are encountered such as liquefaction. The main sources of the constitutive and stability properties of cohesionless granular materials is interparticle friction, which in turn under low confinement stress levels is strongly affected by gravitational body forces under terrestrial (1 gravity) conditions. Under moderate-to-high stress levels, the influence of gravity on the behavior of experiments may not be pronounced and therefore the test results in a terrestrial environment may be acceptable for engineering purposes. However, the conduct of experiments on granular materials under very low stress levels and under quasi-static conditions can only be performed in a microgravity environment. A series of displacement-controlled cyclic triaxial compression experiments were performed in a SPACEHAB module on the Space Shuttle during the STS-79 mission to Mir in September, 1996, and the STS-89 mission in January, 1998. The experiments were conducted on six right cylindrical specimens 75 mm in diameter and 150 mm long at effective confining pressures of 0.05, 0.52 and 1.30 kPa. The results show very high peak strength friction angles in the range of 47.6 to 70.0 degrees, which are mainly due to overconsolidation and grain interlocking effects. It was observed that the residual strength levels were in the same range as that observed at higher confining stress levels. The dilatancy angles were unusually high in the range of 30 to 31 degrees. All specimens display substantial initial stiffnesses and elastic moduli during unloading and reloading events, which are nearly an order of magnitude higher than conventional theories predict. A periodic instability phenomenon which appears to result from buckling of multiple internal arches and columnar systems, augmented by stick-slips was observed in the experiments. Computed Tomography (CT) measurements revealed valuable data about the internal fabric and the specimens deformation patterns. Uniform diffuse bifurcation with multiple radial shear bands was observed in the specimens tested in a microgravity environment. In the axial direction, two major conical surfaces were developed. Spatial nonsymmetrical deformations were observed in specimens tested in terrestrial laboratory.

Alshibli, Khalid A.↗