Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “momentum ordering”

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 19 records

Closure models of turbulent third-order momentum and temperature fluctuations

A study is made for the development and computations of the separating and reattaching shear flows. The highlight of the study is an attempt to predict the third-moments of turbulent velocity which is responsible for the diffusion transport of the Reynolds stresses. The present computations show that the third-moments obtained by employing the law-Reynolds number model of transport equations improve the prediction of the third-moments. The modelling for scalar variables is also performed for the heat transfer computations. Since the transport equations for (u sub j theta)* and (u sub i u sub j theta) have been given in the last NASA CR, the study has been extended further to the modeling of (theta sup 2), (u sub i theta sup 2) and epsilon sub theta (dissipation rate of (theta sup 2)). The formulations are shown in this report.

Amano, R. S.↗

Consistent boundary conditions for reduced Navier-Stokes (RNS) scheme applied to three-dimensional internal viscous flows

Three-dimensional internal viscous flow problems are presently addressed by a consistent and efficient set of boundary conditions for the multisweep space-marching, pressure-elliptic Reduced Navier-Stokes scheme. Since continuity is by this means directly satisfied at all points in the flow domain, the first-order momentum equations are directly solvable for pressure without the requirement for a Poisson pressure-correction equation. Incompressible flow solutions are obtained for straight and curved ducts of square cross section, in order to validate the procedure. Usefulness is demonstrated for internal flows with strong interactions, as would be found in turbomachine geometries.

Reddy, D. R.↗

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗

RIPPLE - A new model for incompressible flows with free surfaces

A new free surface flow model, RIPPLE, is summarized. RIPPLE obtains finite difference solutions for incompressible flow problems having strong surface tension forces at free surfaces of arbitrarily complex topology. The key innovation is the continuum surface force model which represents surface tension as a (strongly) localized volume force. Other features include a higher-order momentum advection model, a volume-of-fluid free surface treatment, and an efficient two-step projection solution method. RIPPLE's unique capabilities are illustrated with two example problems: low-gravity jet-induced tank flow, and the collision and coalescence of two cylindrical rods.

Kothe, D. B.↗

A second-order theory for transverse ion heating and momentum coupling due to electrostatic ion cyclotron waves

A second-order theory for electrostatic instabilities driven by counterstreaming ion beams is developed which describes momentum coupling and heating of the plasma via wave-particle interactions. Exchange rates between the waves and particles are derived, which are suitable for the fluid equations simulating microscopic effects on macroscopic scales. Using a fully kinetic simulation, the electrostatic ion cyclotron instability due to counterstreaming H(+) beams has been simulated. A power spectrum from the kinetic simulation is used to evaluate second-order exchange rates. The calculated heating and momentum loss from second-order theory is compared to the numerical simulation.

Miller, Ronald H.↗

Ab initio calculation of a global potential, vibrational energies, and wave functions for HCN/HNC, and a simulation of the (A-tilde)-(X-tilde) emission spectrum

A potential energy surface for the HCN/HNC system which is a fit to extensive, high-quality ab initio, coupled-cluster calculations is presented. All HCN and HNC states with energies below the energy of the first delocalized state are reported and characterized. Vibrational transition energies are compared with all available experimental data on HCN and HNC, including high CH-overtone states up to 23,063/cm. A simulation of the (A-tilde)-(X-tilde) stimulated emission pumping (SEP) spectrum is also reported, and the results are compared to experiment. Franck-Condon factors are reported for odd bending states of HCN, with one quantum of vibrational angular momentum, in order to compare with the recent assignment by Jonas et al. (1992), on the basis of axis-switching arguments of a number of previously unassigned states in the SEP spectrum.

Bowman, Joel M.↗

The radio source around Eta Carinae

We present high spatial resolution radio observations of the peculiar southern star Eta Carinae, made with the Australian Telescope. The images, at 8 and 9 GHz with a resolution of 1.0 arcsec show a source of dimension 10 arcsec and total flux of 0.7 Jy dominated by a strong central peak. The radio emission is unpolarized and offers no support to models which invoke degenerate stars or more exotic objects within the core of Eta Car. In these data we find no evidence for more than one energy source in the core with arcsecond separations as some infrared observations have suggested. Several levels of structure are evident in the radio image, which shows symmetry on the larger scales. Conventional formulae for stellar wind radio sources give a mass loss rate of order 3 x 10(exp -4) Solar Mass/yr based on the radio flux in the central peak, which yields a wind momentum flux of order 20% of the momentum flux available from the star's radiation field. The radio emission at these frequencies is consistent with thermal emission from gas flowing away from a 'luminous blue variable' star (LBV) Eta Car is probably the brightest thermal stellar wind radio source in the sky.

White, S. M.↗

Time Series Vegetation Aerodynamic Roughness Fields Estimated from MODIS Observations

Most land surface models used today require estimates of aerodynamic roughness length in order to characterize momentum transfer between the surface and atmosphere. The most common method of prescribing roughness is through the use of empirical look-up tables based solely on land cover class. Theoretical approaches that employ satellite-based estimates of canopy density present an attractive alternative to current look-up table approaches based on vegetation cover type that do not account for within-class variability and are oftentimes simplistic with respect to temporal variability. The current research applies Raupach s formulation of momentum aerodynamic roughness to MODIS data on a regional scale in order to estimate seasonally variable roughness and zero-plane displacement height fields using bulk land cover parameters estimated by [Jasinski, M.F., Borak, J., Crago, R., 2005. Bulk surface momentum parameters for satellite-derived vegetation fields. Agric. For. Meteorol. 133, 55-68]. Results indicate promising advances over look-up approaches with respect to characterization of vegetation roughness variability in land surface and atmospheric circulation models.

Borak, Jordan S.↗

Assessment of an Unstructured-Grid Method for Predicting Aerodynamic Performance of Jet Flaps

The application of a Computational Fluid Dynamics tool to a jet flap control effector on an elliptical airfoil-section wing was investigated. The study utilized the Tetrahedral Unstructured Software System developed at NASA Langley Research Center. The Reynolds-averaged Navier-Stokes flow solver code used was USM3D. The CFD-based jet flap simulations were compared to experimental results from a wind tunnel test conducted at the NASA Langley Transonic Dynamics Tunnel. The wind tunnel model consisted of a six percent thick elliptical airfoil with a modified trailing edge. The jet flap was located at 95% chord and exited at 90 degrees to the lower surface. The experimental model was designed to promote two-dimensional flow across the wing. It was found that the CFD simulation had to model the three-dimensional geometry of the experiment in order to obtain good agreement. Tests were performed at two Mach numbers at several different jet momentum coefficients. In order to be consistent with the experimental method, the CFD lift and pitching moment values were determined by integrating the pressures over the wing.

Cruz, Josue↗

Momentum transfer from oblique impacts

A completely satisfactory experiment would be in a low gravity environment where the effect of momentum imparted by ejecta impacting the surface can be removed or controlled from momentum transfer during impact. Preliminary estimates can be made using a ballistic pendulum. Such experiments were initiated at the NASA-Ames Vertical Gun Range in order to examine momentum transfer due to impact vaporization for oblique impacts. The preliminary results indicate that momentum from oblique impacts is very inefficient: decreasing with increasing impact velocity and perhaps size; increasing with decreasing density; and increasing with increasing impact angle. At face value, such results minimize the effect of momentum transfer by grazing impact; the more probable impact angles of 30 deg would have a greater effect, contrary to the commonly held impression.

Schultz, Peter H.↗

A discrete momentum-conserving explicit algorithm for multibody dynamics analysis

A discrete momentum-conserving, explicit time integration is presented. The accurate feature and simplicity of the present algorithm are realized by a mid-point implicit formula for integrating the Euler parameters and a second-order discrete momentum-conserving form of the central difference algorithm, respectively. The accuracy and robustness of the algorithm is demonstrated by example problems which exhibit large overall rigid motions under holonomic constraints.

Park, K. C.↗

A discrete momentum-conserving explicit algorithm for rigid body dynamics analysis

A discrete momentum-conserving explicit time integration is presented. The accurate feature and simplicity of the present algorithm are realized by a mid-point implicit formula for integrating the Euler parameters and a second-order discrete momentum-conserving form of the central difference algorithm, respectively. The accuracy and robustness of the algorithm is demonstrated by example problems which exhibit large overall rigid motions under holonomic constraints.

Park, K. C.↗

A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing

Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.

Bellan, Josette↗

Mature thunderstorm cloud top structure - Three-dimensional numerical simulation versus satellite observations

Preliminary results of four runs with a three-dimensional model of the effects of vertical wind shear on cloud top height/temperature structure and the internal properties of isolate midlatitude thunderstorms are reported. The model is being developed as an aid to analyses of GEO remote sensing satellite data. The grid is a 27 x 27 x 20 mesh with 2 km horizontal resolution and 0.9 vertical resolution. The total grid is 54 km on a side and 18 km deep. A second-order Crowley scheme for advecting momentum is extended with a third-order correction for spatial truncation error, and the earth-relative horizontal surface wind components are decreased to 50 percent of their values at 0.45 km. A temperature increase with height is included, together with an initial impulse consisting of a nonrotating cylindrical weak buoyant updraft 10 km in radius. The results of the runs are discussed in terms of the time variation of the vertical velocity extrema, the effects of strong and weak shear on a storm, the cloud top height, the Lagrangian dynamics of a thermal couplet, and data from a real storm.

Schlesinger, R. E.↗

Geostrophic adjustment in a shallow-water numerical model as it relates to thermospheric dynamics

The theory of geostrophic adjustment and its application to the dynamics of the high latitude thermosphere have been discussed in previous papers based on a linearized treatment of the fluid dynamical equations. However, a linearized treatment is only valid for small Rossby numbers given by Ro = V/fL, where V is the wind speed, f is the local value of the Coriolis parameter, and L is a characteristic horizontal scale for the flow. For typical values in the auroral zone, the approximation is not reasonable for wind speeds greater than 25 m/s or so. A shallow-water (one layer) model was developed that includes the spherical geometry and full nonlinear dynamics in the momentum equations in order to isolate the effects of the nonlinearities on the adjustment process. A belt of accelerated winds between 60 deg and 70 deg latitude was used as the initial condition. The adjustment process was found to proceed as expected from the linear formulation, but that an asymmetry between the response for an eastward and westward flow results from the nonlineawr curvature (centrifugal) terms. In general, the amplitude of an eastward flowing wind will be less after adjustment than a westward wind. For instance, if the initial wind velocity is 300 m/s, the linearized theory predicts a final wind speed of 240 m/s, regardless of the flow direction. However, the nonlinear curvature terms modify the response and produce a final wind speed of only 200 m/s for an initial eastward wind and a final wind speed of almost 300 m/s for an initial westward flow direction. Also, less gravity wave energy is produced by the adjustment of the westward flow than by the adjustment of the eastward flow. The implications are that the response of the thermosphere should be significantly different on the dawn and dusk sides of the auroral oval. Larger flow velocities would be expected on the dusk side since the plasma will accelerate the flow in a westward direction in that sector.

Larsen, M. F.↗

Control of Cavity Resonance Using Steady and Oscillatory Blowing

An experimental study to investigate the effect of steady and oscillatory (with zero net mass flux) blowing on cavity resonance is undertaken. The objective is to study the basic mechanisms of the control of cavity resonance. An actuator is designed and calibrated to generate either steady blowing or oscillatory blowing with A zero net mass flux. The results of the experiment show that both steady and oscillatory blowing are effective, and reduce the amplitude of the dominant resonant mode by 1OdB. The oscillatory blowing is however found to be more superior in that the same effectiveness could be accomplished with a momentum coefficient an order of magnitude smaller than for steady blowing. The experiment also confirms the results of previous computations that suggest the forcing frequency for oscillatory blowing must not be at harmonic frequencies of the cavity resonant modes.

Lamp, Alison M.↗

Towards Automating Spacecraft Attitude Sensor Calibration

With a view towards reducing cost and complexity for spacecraft early mission support at the NASA Goddard Space Flight Center (GSFC), efforts are being made to automate the attitude sensor calibration process. This paper addresses one of the major components needed by such a system. The beneficiaries of an improved calibration process are missions that demand moderate to high precision attitude knowledge or that need to perform accurate attitude slews. Improved slew accuracy reduces the time needed for re-acquisition of fine-pointing after each attitude maneuver, Rapid target acquisition can be very important for astronomical targeting or for off-nadir surface feature targeting by Earth-oriented spacecraft. The normal sequence of on-orbit calibration starts with alignment calibration of the star trackers and possibly the Sun sensor. Their relative alignment needs to be determined using a sufficiently large data set so their fields of view are adequately sampled. Next, the inertial reference unit (IRU) is calibrated for corrections to its alignment and scale factors. The IRU biases are estimated continuously by the onboard attitude control system, but the IRU alignment and scale factors are usually determined on the ground using a batch-processing method on a data set that includes several slews sufficient to give full observability of all the IRU calibration parameters. Finally, magnetometer biases, alignment, and its coupling to the magnetic torquers are determined in order io improve momentum management and occasionally for use in the attitude determination system. The detailed approach used for automating calibrations will depend on whether the automated system resides on the ground or on the spacecraft with an ultimate goal of autonomous calibration. Current efforts focus on a ground-based system driving subsystems that could run either on the ground or onboard. The distinction is that onboard calibration should process the data sequentially rather than in a single large batch since onboard computer data storage is limited. Very good batch- processing calibration utilities have been developed and used extensively at NASA/GSFC for mission support but no sequential calibration utilities are available. To meet this need, this paper presents the mathematical description of a sequential IRU calibration system. The system has been tested using flight data from the Rossi X-ray Timing Explorer (RXTE) during a series of attitude slews. The paper also discusses the current state of the overall automated system and describes plans for adding sequential alignment calibration and other additions that will reduce the amount of analyst time and input.

Sedlak, Joseph↗

Wave Driven Exothermic Heating in the Mesopause Region

A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.

Hickey, Michael P.↗