Results for “Point Kinetics Equations” — Page 3

NASA NTRS1991

Modeling turbulent compressible flows - The mass fluctuating velocity and squared density

This paper deals with single-point closure theory for compressible turbulent flow, including the effects of compressibility on the turbulence. In particular, the combination of the pressure dilatation and the dilatation dissipation, terms which appear on the turbulent kinetic energy equation, are modeled. Model parameters in these transport equations are determined by comparing predictions with boundary layer measurements. Finally, predictions with a k-epsilon model, including the new formulations, are presented for the compressible shear layer.

Taulbee, D.↗

NASA NTRS1990

Flux estimation of the FIFE planetary boundary layer (PBL) with 10.6 micron Doppler lidar

A method is devised for calculating wind, momentum, and other flux parameters that characterize the planetary boundary layer (PBL) and thereby facilitate the calibration of spaceborne vs. in situ flux estimates. Single Doppler lidar data are used to estimate the variance of the mean wind and the covariance related to the vertically pointing fluxes of horizontal momentum. The skewness of the vertical velocity and the range of kinetic energy dissipation are also estimated, and the surface heat flux is determined by means of a statistical Navier-Stokes equation. The conclusion shows that the PBL structure combines both 'bottom-up' and 'top-down' processes suggesting that the relevant parameters for the atmospheric boundary layer be revised. The conclusions are of significant interest to the modeling techniques used in General Circulation Models as well as to flux estimation.

Gal-Chen, Tzvi↗

NASA NTRS2008

Flow Solver for Incompressible 2-D Drive Cavity

This software solves the Navier-Stokes equations for the incompressible driven cavity flow problem. The code uses second-order finite differencing on a staggered grid using the Chorin projection method. The resulting intermediate Poisson equation is efficiently solved using the fast Fourier transform. Time stepping is done using fourth-order Runge-Kutta for stability at high Reynolds numbers. Features include check-pointing, periodic field snapshots, ongoing reporting of kinetic energy and changes between time steps, time histories at selected points, and optional streakline generation.

Kalb, Virginia L.↗

NASA NTRS2011

The Prediction of Broadband Shock-Associated Noise Including Propagation Effects

An acoustic analogy is developed based on the Euler equations for broadband shock- associated noise (BBSAN) that directly incorporates the vector Green's function of the linearized Euler equations and a steady Reynolds-Averaged Navier-Stokes solution (SRANS) as the mean flow. The vector Green's function allows the BBSAN propagation through the jet shear layer to be determined. The large-scale coherent turbulence is modeled by two-point second order velocity cross-correlations. Turbulent length and time scales are related to the turbulent kinetic energy and dissipation. An adjoint vector Green's function solver is implemented to determine the vector Green's function based on a locally parallel mean flow at streamwise locations of the SRANS solution. However, the developed acoustic analogy could easily be based on any adjoint vector Green's function solver, such as one that makes no assumptions about the mean flow. The newly developed acoustic analogy can be simplified to one that uses the Green's function associated with the Helmholtz equation, which is consistent with the formulation of Morris and Miller (AIAAJ 2010). A large number of predictions are generated using three different nozzles over a wide range of fully expanded Mach numbers and jet stagnation temperatures. These predictions are compared with experimental data from multiple jet noise labs. In addition, two models for the so-called 'fine-scale' mixing noise are included in the comparisons. Improved BBSAN predictions are obtained relative to other models that do not include the propagation effects, especially in the upstream direction of the jet.

Miller, Steven↗

NASA NTRS1998

Jet Aeroacoustics: Noise Generation Mechanism and Prediction

This report covers the third year research effort of the project. The research work focussed on the fine scale mixing noise of both subsonic and supersonic jets and the effects of nozzle geometry and tabs on subsonic jet noise. In publication 1, a new semi-empirical theory of jet mixing noise from fine scale turbulence is developed. By an analogy to gas kinetic theory, it is shown that the source of noise is related to the time fluctuations of the turbulence kinetic theory. On starting with the Reynolds Averaged Navier-Stokes equations, a formula for the radiated noise is derived. An empirical model of the space-time correlation function of the turbulence kinetic energy is adopted. The form of the model is in good agreement with the space-time two-point velocity correlation function measured by Davies and coworkers. The parameters of the correlation are related to the parameters of the k-epsilon turbulence model. Thus the theory is self-contained. Extensive comparisons between the computed noise spectrum of the theory and experimental measured have been carried out. The parameters include jet Mach number from 0.3 to 2.0 and temperature ratio from 1.0 to 4.8. Excellent agreements are found in the spectrum shape, noise intensity and directivity. It is envisaged that the theory would supercede all semi-empirical and totally empirical jet noise prediction methods in current use.

Tam, Christopher↗

NASA NTRS2009

Device for Measuring Low Flow Speed in a Duct

A multiple-throat venturi system has been invented for measuring laminar flow of air or other gas at low speed (1 to 30 cm/s) in a duct while preserving the laminar nature of the flow and keeping the velocity profile across the duct as nearly flat as possible. While means for measuring flows at higher speeds are well established, heretofore, there have been no reliable means for making consistent, accurate measurements in this speed range. In the original application for which this system was invented, the duct leads into the test section of a low-speed wind tunnel wherein uniform, low-speed, laminar flow is required for scientific experiments. The system could also be used to monitor a slow flow of gas in an industrial process like chemical vapor deposition. In the original application, the multiple- throat venturi system is mounted at the inlet end of the duct having a rectangular cross section of 19 by 14 cm, just upstream of an assembly of inlet screens and flow straighteners that help to suppress undesired flow fluctuations (see Figure 1). The basic venturi measurement principle is well established: One measures the difference in pressure between (1) a point just outside the inlet, where the pressure is highest and the kinetic energy lowest; and (2) the narrowest part (the throat) of the venturi passage, where the kinetic energy is highest and the pressure is lowest. Then by use of Bernoulli s equation for the relationship between pressure and kinetic energy, the volumetric flow speed in the duct can be calculated from the pressure difference and the inlet and throat widths. The design of this system represents a compromise among length, pressure recovery, uniformity of flow, and complexity of assembly. Traditionally, venturis are used to measure faster flows in narrower cross sections, with longer upstream and downstream passages to maintain accuracy. The dimensions of the passages of the present venturi system are sized to provide a readily measurable pressure drop. Multiple throats are used to minimize the length needed to recover internal energy and enable the velocity profile to recover to near flatness.

Quinn, Frank↗

NASA NTRS2009

Measuring Thermal Conductivity of a Small Insulation Sample

A multiple-throat venturi system has been invented for measuring laminar flow of air or other gas at low speed (1 to 30 cm/s) in a duct while preserving the laminar nature of the flow and keeping the velocity profile across the duct as nearly flat as possible. While means for measuring flows at higher speeds are well established, heretofore, there have been no reliable means for making consistent, accurate measurements in this speed range. In the original application for which this system was invented, the duct leads into the test section of a low-speed wind tunnel wherein uniform, low-speed, laminar flow is required for scientific experiments. The system could also be used to monitor a slow flow of gas in an industrial process like chemical vapor deposition. In the original application, the multiple- throat venturi system is mounted at the inlet end of the duct having a rectangular cross section of 19 by 14 cm, just upstream of an assembly of inlet screens and flow straighteners that help to suppress undesired flow fluctuations (see Figure 1). The basic venturi measurement principle is well established: One measures the difference in pressure between (1) a point just outside the inlet, where the pressure is highest and the kinetic energy lowest; and (2) the narrowest part (the throat) of the venturi passage, where the kinetic energy is highest and the pressure is lowest. Then by use of Bernoulli s equation for the relationship between pressure and kinetic energy, the volumetric flow speed in the duct can be calculated from the pressure difference and the inlet and throat widths. The design of this system represents a compromise among length, pressure recovery, uniformity of flow, and complexity of assembly. Traditionally, venturis are used to measure faster flows in narrower cross sections, with longer upstream and downstream passages to maintain accuracy. The dimensions of the passages of the present venturi system are sized to provide a readily measurable pressure drop. Multiple throats are used to minimize the length needed to recover internal energy and enable the velocity profile to recover to near flatness.

Miller, Robert A.↗

NASA NTRS1991

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup - alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state towards which a complete self-preseving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗

NASA NTRS1992

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup -alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state toward which a complete self-preserving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗

NASA NTRS2014

Properties of Desert Sand and CMAS Glass

As-received desert sand from a Middle East country has been characterized for its phase composition and thermal stability. X-ray diffraction analysis showed the presence of quartz (SiO2), calcite (CaCO3), gypsum (CaSO4.2H2O), and NaAlSi3O8 phases in as-received desert sand and showed weight loss of approx. 35 percent due to decomposition of CaCO3 and CaSO4.2H2O when heated to 1400 C. A batch of as-received desert sand was melted into calcium magnesium aluminosilicate (CMAS) glass at approx. 1500 C. From inductively coupled plasma-atomic emission spectrometry, chemical composition of the CMAS glass was analyzed to be 27.8CaO-4MgO-5Al2O3-61.6SiO2-0.6Fe2O3-1K2O (mole percent). Various physical, thermal and mechanical properties of the glass have been evaluated. Bulk density of CMAS glass was 2.69 g/cc, Young's modulus 92 GPa, Shear modulus 36 GPa, Poisson's ratio 0.28, dilatometric glass transition temperature (T (sub g)) 706 C, softening point (T (sub d)) 764 C, Vickers microhardness 6.3 +/- 0.4 GPa, indentation fracture toughness 0.75 +/- 0.15 MPa.m (sup 1/2), and coefficient of thermal expansion (CTE) 9.8 x 10 (exp -6)/degC in the temperature range 25 to 700 C. Temperature dependence of viscosity has also been estimated from various reference points of the CMAS glass using the Vogel-Fulcher-Tamman (VFT) equation. The glass remained amorphous after heat treating at 850 C for 10 hr but crystallized into CaSiO3 and Ca-Mg-Al silicate phases at 900 C or higher temperatures. Crystallization kinetics of the CMAS glass has also been investigated by differential thermal analysis (DTA). Activation energies for the crystallization of two different phases in the glass were calculated to be 403 and 483 kJ/mol, respectively.

Desert sand↗

NASA NTRS1986

Analysis of Doppler lidar wind measurements

Doppler lidar and multiple Doppler radar data were obtained in a convectively mixed planetary boundary layer. The lidar measurements were possible due to scattering from existing aerosols; radar reflecting chaff was released in the atmosphere to make it visible to the multiple Doppler radar network. The data were analyzed to obtain detailed horizontal wind structures. The divergence of the horizontal wind was calculated and the anelastic continuity equation integrated to obtain vertical air motions. Differences between the areally averaged quantities and the grid point values provided a measure of the fluctuations in the wind components or the turbulent wind fluctuations. Vertical profiles of the mean winds and quantities related to the turbulent kinetic energy components and the turbulent momentum transfers were also calculated.

Srivastava, R. C.↗

DOE OSTI2025

Mean-field dynamo as a quantum-like modulational instability

Presented here is a novel formulation of the mean-field dynamo as a modulational instability of magnetohydrodynamic (MHD) turbulence. This formulation, termed mean-field wave kinetics (MFWK), is based on the Weyl symbol calculus and allows describing the interaction between the mean fields (magnetic field and fluid velocity) and turbulence without requiring scale separation that is commonly assumed in the literature. The turbulence is described by the Wigner–Moyal equation for the spectrum of the two-point correlation matrix (Wigner matrix) of magnetic-field and velocity fluctuations and depicts the turbulence as an effective plasma of quantum-like particles that interact via the mean fields. Eddy–eddy interactions, which serve as ‘collisions’ in this effective plasma, are modelled within the standard minimal tau approximation to aid comparison with existing theories. Using MFWK, the non-local electromotive force is calculated for generic turbulence from first principles, modulo the limitations of MFWK. This result is then used to study, both analytically and numerically, the modulational modes of MHD turbulence, which appear as linear instabilities of the said effective quantum-like plasma of fluctuations. The standard α 2 -dynamo and other known results are reproduced as special cases. A new dynamo effect is predicted that is driven by correlations between the turbulent flow velocity and the turbulent current.

astrophysical plasmas↗

NASA NTRS1988

Correlations of velocity and temperature fluctuations in the stagnation-point flow of circular cylinder in turbulent flow

Boundary layer flow and turbulence transport analyses to study the influence of the free-stream turbulence on the surface heat transfer rate and the skin friction around the stagnation point of a circular cylinder in a turbulent flow are presented. The analyses are formulated with the turbulent boundary layer equations, the Reynolds stress transport equations and the k - epsilon two-equation turbulence modeling. The analyses are used to calculate the time-averaged turbulence double correlations, the mean flow properties, the surface heat transfer rate and the skin friction with an isotropic turbulence in the freestream. The analytical results are described and compared with the existing experimental measurements. Depending on the free-stream turbulence properties, the turbulence kinetic energy can increase or decrease as the flow moves toward the surface. However, the turbulence kinetic energy induces large Reynolds normal stresses at the boundary layer edge. The Reynolds normal stresses change the boundary layer profiles of the time-averaged double correlations of the velocity and temperature fluctuations, the surface heat transfer rate and the skin friction. The free-stream turbulence dissipation rate can affect the stagnation-point heat transfer rate but the influence of the free-stream temperature fluctuation on the heat transfer rate is insignificant.

Wang, Chi R.↗

NASA NTRS1977

A special purpose symbolic manipulation program for the derivation of the equations of motion for large flexible structures

This paper outlines a computer program especially tailored to the task of deriving explicit equations of motion for structures with point-connected substructures. The special purpose program is written in FORTRAN and is designed for performing the specific algebraic operations encountered in the derivation of explicit equations of motion. The derivation is by the Lagrangian approach. Using an orderly kinematical procedure and a discretization and/or truncation scheme, it is possible to write the kinetic and potential energy of each substructure in a compact vector-matrix form. Then, if each element of the matrices and vectors encountered in the kinetic and potential energy is a known algebraic expression, the computer program performs the necessary operations to evaluate the kinetic and potential energy of the system explicitly. Lagrange's equations for small motions about equilibrium can be deduced directly from the explicit form of the system kinetic and potential energy.

Hale, A. L.↗

NASA NTRS1978

Derivation of the equations of motion for complex structures by symbolic manipulation

This paper outlines a computer program especially tailored to the task of deriving explicit equations of motion for structures with point-connected substructures. The special purpose program is written in FORTRAN and is designed for performing the specific algebraic operations encountered in the derivation of explicit equations of motion. The derivation is by the Lagrangian approach. Using an orderly kinematical procedure and a discretization and/or truncation scheme, it is possible to write the kinetic and potential energy of each substructure in a compact vector-matrix form. Then, if each element of the matrices and vectors encountered in the kinetic and potential energy is a known algebraic expression, the computer program performs the necessary operations to evaluate the kinetic and potential energy of the system explicitly. Lagrange's equations for small motions about equilibrium can be deduced directly from the explicit form of the system kinetic and potential energy.

Hale, A. L.↗

NASA NTRS2003

Use of GTE, Mozaic, Sonex, and UARS-MLS Data in Understanding Tropospheric Processes Critical to Model Development

During the first year we focused on the analysis of data collected on over 7600 commercial aircraft flights (the MOZAIC program). The aim was to further our understanding of the fundamental dynamical processes that drive mesoscale phenomena in the upper troposphere and lower stratosphere, and their effects on the advection of passive scalars. Through these studies we made the following findings. 2001]: We derived the Kolmogorov equation for the third-order velocity structure function on an f-plane. We showed how the sign of the function yields the direction of the energy cascade. The remarkable linearity of the measured off diagonal third-order structure function was studied. We suggested that the Coriolis term, which appears explicitly in this equation, may be crucial in understanding the observed kinetic energy spectra at scales larger than 100 km, instead of the nonlinear advection term as previously assumed. Also, we showed that decreases with latitude for mid- to high-latitudes, and tabulated the values. Ozone concentration structure functions were calculated by restricting the data points to approximately isentropic pairings. In this way we were able to make comparisons with scalar advection theories. We found that, at even the largest scales, there was no evidence for the simple scaling predicted for smooth advection/diffusion.

Newell, Reginald E.↗

NASA NTRS1983

A linear analysis of the transition curve for the baroclinic annulus

A primitive equation linear model is shown to be capable of quantitatively predicting the neutral curve for a baroclinic annulus. Flows unstable to small nonaxisymmetric perturbations are assumed to lie on the nonaxisymmetric side of a transition curve, while the stable flows are on the axisymmetric side. The stability in the presence of perturbations is tested with linearized hydrostatic primitive equations. Attention is given to the fastest growing eigenmodes, demonstrating that the linear and nonlinear wave have similar structures at a particular point in the nonaxisymmetric regime. All the waves are similar to the Eddy wave except at the extreme to the transition curve, where the waves maintain little structure with height. The eddy kinetic energy is baroclinic with a potential energy suppresed by the large static stability of the basic state. Eddy dissipation is noted to be significant near all boundaries.

Miller, T. L.↗

NASA NTRS1984

Shear deformation plate continua of large double-layered space structures

The energy equivalence to construct equivalent continua for the actual lattice structure is described. An energy equivalent continuum is defined as that which has the same amount of strain and kinetic energies stored in it as the original lattice structure when both are subjected to the same loading conditions. The equivalent continuum is characterized by its strain and kinetic energies from which the constitutive relations and the equations of motion can be derived. A simple method to model large rigid jointed lattices as continuous media with couple stresses is presented. The transition from the discrete system to the continuous medium is achieved by expanding the displacements and the rotations of the nodal points in a Taylor series about a suitable chosen origin. Basic kinematic assumptions are introduced to insure that the assumptions used in deriving the governing equations of the modeled continuum are satisfied. The number of terms retained in the Taylor series expansion will depend upon the properties to be evaluated. This implies that the kind of continuum needed to model from the discrete lattice, before the actual properties are derived was to predetermined.

Hefzy, M. S.↗