Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scalar field data”

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 127 records · Page 7

Optical interferometry in fluid dynamics research

Optical interferometry techniques were applied to the investigation of transonic airfoil flow fields in large wind tunnels. Holographic interferometry techniques were used to study 2 dimensional symmetric NACA 64A010 and Douglas Aircraft Co. DSMA671 supercritical airfoil performance in the NASA Ames 2 x 2 ft transonic wind tunnel. Quantitative data obtained from the interferograms were compared to the surface pressure data. The agreement obtained verified the accuracy of the flow visualization and demonstrated the potential for acquiring quantitative scalar results. Measurements of the inviscid flow speed and the boundary layer and wake velocity profiles were extracted from the interferograms and compared to laser Doppler velocimeter measurements. These results were also in good agreement. A method for acquiring real time interferometric data in large scale facilities was developed. This method, based on the point diffraction interferometer, was successfully tested in the 2 x 2 ft transonic wind tunnel. The holographic and real time interferometry methods were applied to the investigations of circulation control airfoils utilizing the Coanda effect. These results reveals the details of the jet interacting with the trailing edge boundary layer and the other parameters affecting the lift augmentation.

Bachalo, W. D.↗

Visualization and Quantification of Rotor Tip Vortices in Helicopter Flows

Helicopter aeromechanics encompasses a highly vortical flow field. The vortices generated at each blade tip contain unsteady, complex, three-dimensional structures, which interact with each other, other blades, the fuselage and various components of the helicopter. It is crucial to understand vortex kinematics and their subsequent dynamic evolution. Much research has been devoted to the understanding of helicopter vortex dynamics, including a number of experimental studies.1-6 In May 2010 Particle Image Velocimetry (PIV) measurements of a full-scale UH-60A rotor were acquired in the National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel.1 These measurements were taken at a plane just downstream of the advancing blade in the vicinity of the blade tip—the so-called PIV plane. The resulting PIV data were then processed using an ensemble-average approach to create graphical representations of the vortical wake velocity and vorticity fields, which, in turn, have enhanced the understanding of rotorcraft vortical wake flow field physics and have provided a more detailed validation of vortical wake computer simulations.7 A common approach used to analyze flow field features is to compute and plot color contour maps of various scalar quantities such as pressure, velocity magnitude and vorticity magnitude. For example, the color map of the vorticity magnitude is typically used to determine vortical flow structure. With this approach the vortex core may appear larger or smaller, depending on the contour levels that are selected. Thus, the resulting visualization is sensitive to user-specified contour levels. For vortex core radius measurements, it is more accurate to calculate the vortex core radius using the cross-flow velocity profile across the vortex core. The task of extracting the cross-flow velocity profile can be time consuming with existing tools since the user needs to manually select the core center then specify sampling points along the profile axis. The task becomes even more challenging when the associated grid system uses AMR (Adaptive Mesh Refinement) where the profile axis could span multiple grid blocks. There are a number of existing techniques for profiling of vortex core attributes;8-9 however, these techniques are not fully automatic in that the user still needs to select the vortex core center to compute the cross-flow velocity profile. The present study introduces a new color map scheme that is based on the vortex core radius, which is fully automatic and does not require user intervention. Analysis and visualization of blade tip vortices on the PIV plane using the proposed new color map scheme are described in Section II. The new approach is evaluated using two case studies, which are described in Section III. The paper ends with a summary in Section IV.

Quantification↗

Theoretical frameworks for testing relativistic gravity. IV - A compendium of metric theories of gravity and their post-Newtonian limits.

Metric theories of gravity are compiled and classified according to the types of gravitational fields they contain, and the modes of interaction among those fields. The gravitation theories considered are classified as (1) general relativity, (2) scalar-tensor theories, (3) conformally flat theories, and (4) stratified theories with conformally flat space slices. The post-Newtonian limit of each theory is constructed and its Parametrized Post-Newtonian (PPN) values are obtained by comparing it with Will's version of the formalism. Results obtained here, when combined with experimental data and with recent work by Nordtvedt and Will and by Ni, show that, of all theories thus far examined by our group, the only currently viable ones are general relativity, the Bergmann-Wagoner scalar-tensor theory and its special cases (Nordtvedt; Brans-Dicke-Jordan), and a recent, new vector-tensor theory by Nordtvedt, Hellings, and Will.

Ni, W.-T.↗

The study of PDF turbulence models in combustion

In combustion computations, it is known that the predictions of chemical reaction rates are poor if conventional turbulence models are used. The probability density function (pdf) method seems to be the only alternative that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus is the only viable approach for more accurate turbulent combustion calculations. The fact that the pdf equation has a very large dimensionality renders finite difference schemes extremely demanding on computer memories and thus impractical. A logical alternative is the Monte Carlo scheme. Since CFD has a certain maturity as well as acceptance, it seems that the use of a combined CFD and Monte Carlo scheme is more beneficial. Therefore, a scheme is chosen that uses a conventional CFD flow solver in calculating the flow field properties such as velocity, pressure, etc., while the chemical reaction part is solved using a Monte Carlo scheme. The discharge of a heated turbulent plane jet into quiescent air was studied. Experimental data for this problem shows that when the temperature difference between the jet and the surrounding air is small, buoyancy effect can be neglected and the temperature can be treated as a passive scalar. The fact that jet flows have a self-similar solution lends convenience in the modeling study. Futhermore, the existence of experimental data for turbulent shear stress and temperature variance make the case ideal for the testing of pdf models wherein these values can be directly evaluated.

Hsu, Andrew T.↗

Effects of spacecraft potential on three-dimensional electron measurements in the solar wind

Using the three-dimensional, low-energy electron spectrometer aboard the Ulysses spacecraft, we have measured the gyrotropicity of electron distributions in the solar wind. In order to make these observations, we have developed a new technique for correcting spacecraft charging effects in three-dimensional, low-energy particle measurements. Comparisons of ion and electron number and current densities, and the alignment of electron temperature anisotropies with the local magnetic field, are presented as evidence of the improvement in the accuracy of the electron moments resulting from the spacecraft charging corrections. The implications of our charging correction technique go beyond simple scalar corrections to the Ulysses measurements. We discuss the effects of our charging correction upon the measurements of temporal and radial gradients in a plasma environment and for two-dimensionally obtained low-energy particle data.

Scime, Earl E.↗

Visualizing Vector Fields Using Line Integral Convolution and Dye Advection

We present local and global techniques to visualize three-dimensional vector field data. Using the Line Integral Convolution (LIC) method to image the global vector field, our new algorithm allows the user to introduce colored 'dye' into the vector field to highlight local flow features. A fast algorithm is proposed that quickly recomputes the dyed LIC images. In addition, we introduce volume rendering methods that can map the LIC texture on any contour surface and/or translucent region defined by additional scalar quantities, and can follow the advection of colored dye throughout the volume.

Shen, Han-Wei↗

Source Parameter Inversion for Recent Great Earthquakes from a Decade-long Observation of Global Gravity Fields

We quantify gravity changes after great earthquakes present within the 10 year long time series of monthly Gravity Recovery and Climate Experiment (GRACE) gravity fields. Using spherical harmonic normal-mode formulation, the respective source parameters of moment tensor and double-couple were estimated. For the 2004 Sumatra-Andaman earthquake, the gravity data indicate a composite moment of 1.2x10(exp 23)Nm with a dip of 10deg, in agreement with the estimate obtained at ultralong seismic periods. For the 2010 Maule earthquake, the GRACE solutions range from 2.0 to 2.7x10(exp 22)Nm for dips of 12deg-24deg and centroid depths within the lower crust. For the 2011 Tohoku-Oki earthquake, the estimated scalar moments range from 4.1 to 6.1x10(exp 22)Nm, with dips of 9deg-19deg and centroid depths within the lower crust. For the 2012 Indian Ocean strike-slip earthquakes, the gravity data delineate a composite moment of 1.9x10(exp 22)Nm regardless of the centroid depth, comparing favorably with the total moment of the main ruptures and aftershocks. The smallest event we successfully analyzed with GRACE was the 2007 Bengkulu earthquake with M(sub 0) approx. 5.0x10(exp 21)Nm. We found that the gravity data constrain the focal mechanism with the centroid only within the upper and lower crustal layers for thrust events. Deeper sources (i.e., in the upper mantle) could not reproduce the gravity observation as the larger rigidity and bulk modulus at mantle depths inhibit the interior from changing its volume, thus reducing the negative gravity component. Focal mechanisms and seismic moments obtained in this study represent the behavior of the sources on temporal and spatial scales exceeding the seismic and geodetic spectrum.

Parameter Inversion↗

A scalar/vector potential formulation of an airfoil in nonuniform stream

A recent innovation in aerodynamic calculation techniques is the formulation based upon scalar and vector potentials (Stokes-Helmholtz decomposition) for the velocity field. This technique is presented: an implementation involving classical panel method is studied. The application of interest is to two-dimensional airfoils moving into a nonuniform approach flow. Comparisons with theoretical and numerical results are included. Furthermore the variations of lift and moment coefficients of quasi-steady simulated flight of an airfoil through microburst wind shear data are studied.

Wang, J. C. A.↗

Accurate Predictions of Mean Geomagnetic Dipole Excursion and Reversal Frequencies, Mean Paleomagnetic Field Intensity, and the Radius of Earth's Core Using McLeod's Rule

The geomagnetic spatial power spectrum R(sub n)(r) is the mean square magnetic induction represented by degree n spherical harmonic coefficients of the internal scalar potential averaged over the geocentric sphere of radius r. McLeod's Rule for the magnetic field generated by Earth's core geodynamo says that the expected core surface power spectrum (R(sub nc)(c)) is inversely proportional to (2n + 1) for 1 less than n less than or equal to N(sub E). McLeod's Rule is verified by locating Earth's core with main field models of Magsat data; the estimated core radius of 3485 kn is close to the seismologic value for c of 3480 km. McLeod's Rule and similar forms are then calibrated with the model values of R(sub n) for 3 less than or = n less than or = 12. Extrapolation to the degree 1 dipole predicts the expectation value of Earth's dipole moment to be about 5.89 x 10(exp 22) Am(exp 2)rms (74.5% of the 1980 value) and the expected geomagnetic intensity to be about 35.6 (mu)T rms at Earth's surface. Archeo- and paleomagnetic field intensity data show these and related predictions to be reasonably accurate. The probability distribution chi(exp 2) with 2n+1 degrees of freedom is assigned to (2n + 1)R(sub nc)/(R(sub nc). Extending this to the dipole implies that an exceptionally weak absolute dipole moment (less than or = 20% of the 1980 value) will exist during 2.5% of geologic time. The mean duration for such major geomagnetic dipole power excursions, one quarter of which feature durable axial dipole reversal, is estimated from the modern dipole power time-scale and the statistical model of excursions. The resulting mean excursion duration of 2767 years forces us to predict an average of 9.04 excursions per million years, 2.26 axial dipole reversals per million years, and a mean reversal duration of 5533 years. Paleomagnetic data show these predictions to be quite accurate. McLeod's Rule led to accurate predictions of Earth's core radius, mean paleomagnetic field intensity, and mean geomagnetic dipole power excursion and axial dipole reversal frequencies. We conclude that McLeod's Rule helps unify geo-paleomagnetism, correctly relates theoretically predictable statistical properties of the core geodynamo to magnetic observation, and provides a priori information required for stochastic inversion of paleo-, archeo-, and/or historical geomagnetic measurements.

Voorhies, Coerte V.↗

The convection electric field model for the magnetosphere based on Explorer 45 observations

Large-scale electric fields in the magnetosphere, represented by the scalar potentials of corotation and convection fields, are examined in order to interpret both plasmapause positions at different local times and energetic particle penetrations inside the plasmapause observed by the Explorer 45 (S3-A) dc E field measurements and particle detectors. The Volland-Stern type convection electric field is assumed. A uniform dawn-dusk convection electric field which corresponds to an exponent (gamma) of unity is not consistent with the observational data, but the Volland-Stern model with gamma = 2 more likely represents the convection electric field during the hours of electric field increases during magnetic activity. In the time-independent case, the shape of the plasmapause does not depend on the intensity of the convection electric field, so that by using two plasmapause distances at different local times along the S3 outbound and inbound orbits the shape factor gamma of the convection field could be determined. The average value of gamma for 42 cases was found to be 2.4.

Ejiri, M.↗

Optical interferometry in fluid dynamics research

Optical interferometry techniques have been applied to the investigation of transonic airfoil flow fields in large-scale wind tunnels. Holographic interferometry techniques were used in the study of two-dimensional symmetric NACA 64A010 and Douglas Aircraft Company DSMA671 supercritical airfoil performance in the NASA Ames 2 ft x 2 ft transonic wind tunnel. Quantitative data obtained from the interferograms were compared to the surface pressure data. The excellent agreement obtained verified the accuracy of the flow visualization and demonstrated the potential for acquiring quantitative scalar results. Measurements of the inviscid flow speed and the boundary layer and wake velocity profiles were extracted from the interferograms and compared to laser Doppler velocimeter measurements. These results were also in good agreement. A method for acquiring real-time interferometric data in large-scale facilities was developed. This method, based on the point diffraction interferometer, was successfully tested in the Ames 2 ft x 2 ft transonic wind tunnel. The holographic and real-time interferometry methods were applied to the investigations of circulation control airfoils utilizing the Coanda effect. These results revealed the details of the jet interaction with the trailing edge boundary layer and the other parameters affecting the lift augmentation.

Bachalo, W. D.↗

A Model for Jet-Surface Interaction Noise Using Physically Realizable Upstream Turbulence Conditions

This paper is a continuation of previous work in which a generalized Rapid Distortion Theory (RDT) formulation was used to model low-frequency trailing-edge noise. The research was motivated by proposed next-generation aircraft configurations where the exhaust system is tightly integrated with the airframe. Data from recent experiments at NASA on the interaction between high-Reynolds-number subsonic jet flows and an external flat plate showed that the power spectral density (PSD) of the far-field pressure underwent considerable amplification at low frequencies. For example, at the 900 observation angle, the low-frequency noise could be as much as 10dB greater than the jet noise itself. In this paper, we present predictions of the noise generated by the interaction of a rectangular jet with the trailing edge of a semi-infinite flat plate. The calculations are based on a formula for the acoustic spectrum of this noise source derived from an exact formal solution of the linearized Euler equations involving (in this case) one arbitrary convected scalar quantity and a Rayleigh equation Green's function. A low-frequency asymptotic approximation for the Green's function based on a two-dimensional mean flow is used in the calculations along with a physically realizable upstream turbulence spectrum, which includes a finite de-correlation region. Numerical predictions, based on three-dimensional RANS solutions for a range of subsonic acoustic Mach number jets and nozzle aspect ratios are compared with experimental data. Comparisons of the RANS results with flow data are also presented for selected cases. We find that a finite decorrelation region increases the low-frequency algebraic decay (the low frequency "rolloff") of the acoustic spectrum with angular frequency thereby producing much closer agreement with noise data for Strouhal numbers less than 0.1. Secondly, the large-aspectratio theory is able to predict the low-frequency amplification due to the jet-edge interaction reasonably well, even for moderate aspect ratio nozzles. We show also that the noise predictions for smaller aspect ratio jets can be fine-tuned using the appropriate RANS-based mean flow and turbulence properties.

Jet↗

A Model for Jet-Surface Interaction Noise Using Physically Realizable Upstream Turbulence Conditions

This paper is a continuation of previous work in which a generalized Rapid Distortion Theory (RDT) formulation was used to model low-frequency trailing-edge noise. The research was motivated by proposed next-generation aircraft configurations where the exhaust system is tightly integrated with the airframe. Data from recent experiments at NASA on the interaction between high-Reynolds-number subsonic jet flows and an external flat plate showed that the power spectral density (PSD) of the far-field pressure underwent considerable amplification at low frequencies. For example, at the 900 observation angle, the low-frequency noise could be as much as 10dB greater than the jet noise itself. In this paper, we present predictions of the noise generated by the interaction of a rectangular jet with the trailing edge of a semi-infinite flat plate. The calculations are based on a formula for the acoustic spectrum of this noise source derived from an exact formal solution of the linearized Euler equations involving (in this case) one arbitrary convected scalar quantity and a Rayleigh equation Green's function. A low-frequency asymptotic approximation for the Green's function based on a two-dimensional mean flow is used in the calculations along with a physically realizable upstream turbulence spectrum, which includes a finite de-correlation region. Numerical predictions, based on three-dimensional RANS solutions for a range of subsonic acoustic Mach number jets and nozzle aspect ratios are compared with experimental data. Comparisons of the RANS results with flow data are also presented for selected cases. We find that a finite decorrelation region increases the low-frequency algebraic decay (the low frequency "rolloff") of the acoustic spectrum with angular frequency thereby producing much closer agreement with noise data for Strouhal numbers less than 0.1. Secondly, the large-aspectratio theory is able to predict the low-frequency amplification due to the jet-edge interaction reasonably well, even for moderate aspect ratio nozzles. We show also that the noise predictions for smaller aspect ratio jets can be fine-tuned using the appropriate RANS-based mean flow and turbulence properties.

Aircraft↗

A transport model of the turbulent scalar-velocity

Performance tests of the third-order turbulence closure for predictions of separating and recirculating flows in backward-facing steps were studied. Computations of the momentum and temperature fields in the flow domain being considered entail the solution of time-averaged transport equations containing the second-order turbulent fluctuating products. The triple products, which are responsible for the diffusive transport of the second-order products, attain greater significance in separating and reattaching flows. The computations are compared with several algebraic models and with the experimental data. The prediction was improved considerably, particularly in the separated shear layer. Computations are further made for the temperature-velocity double products and triple products. Finally, several advantages were observed in the usage of the transport equations for the evaluation of the turbulence triple products; one of the most important features is that the transport model can always take the effects of convection and diffusion into account in strong convective shear flows such as reattaching separated layers while conventional algebraic models cannot account for these effects in the evaluation of turbulence variables.

Amano, R. S.↗

Extra high speed modified Lundell alternator parameters and open/short-circuit characteristics from global 3D-FE magnetic field solutions

The combined magnetic vector potential - magnetic scalar potential method of computation of 3D magnetic fields by finite elements, introduced in a companion paper, is used for global 3D field analysis and machine performance computations under open-circuit and short-circuit conditions for an example 14.3 kVA modified Lundell alternator, whose magnetic field is of intrinsic 3D nature. The computed voltages and currents under these machine test conditions were verified and found to be in very good agreement with corresponding test data. Results of use of this modelling and computation method in the study of a design alteration example, in which the stator stack length of the example alternator is stretched in order to increase voltage and volt-ampere rating, are given here. These results demonstrate the inadequacy of conventional 2D-based design concepts and the imperative of use of this type of 3D magnetic field modelling in the design and investigation of such machines.

Wang, R.↗

Very Large Eddy Simulations of a Jet-A Spray Reacting Flow in a Single Element LDI Injector With and Without Invoking an Eulerian Scalar DWFDF Method

This paper presents the very large eddy simulations (VLES) of a Jet-A spray reacting flow in a single element lean direct injection (LDI) injector by using the National Combustion Code (NCC) with and without invoking the Eulerian scalar DWFDF method, in which DWFDF is defined as the density weighted time filtered fine grained probability density function. The flow field is calculated by using the time filtered compressible Navier-Stokes equations (TFNS) with nonlinear subscale turbulence models, and when the Eulerian scalar DWFDF method is invoked, the energy and species mass fractions are calculated by solving the equation of DWFDF. A nonlinear subscale model for closing the convection term of the Eulerian scalar DWFDF equation is used and will be briefly described in this paper. Detailed comparisons between the results and available experimental data are carried out. Some positive findings of invoking the Eulerian scalar DWFDF method in both improving the simulation quality and maintaining economic computing cost are observed.

Shih, Tsan-Hsing↗

Geologic analysis of averaged magnetic satellite anomalies

To investigate relative advantages and limitations for quantitative geologic analysis of magnetic satellite scalar anomalies derived from arithmetic averaging of orbital profiles within equal-angle or equal-area parallelograms, the anomaly averaging process was simulated by orbital profiles computed from spherical-earth crustal magnetic anomaly modeling experiments using Gauss-Legendre quadrature integration. The results indicate that averaging can provide reasonable values at satellite elevations, where contributing error factors within a given parallelogram include the elevation distribution of the data, and orbital noise and geomagnetic field attributes. Various inversion schemes including the use of equivalent point dipoles are also investigated as an alternative to arithmetic averaging. Although inversion can provide improved spherical grid anomaly estimates, these procedures are problematic in practice where computer scaling difficulties frequently arise due to a combination of factors including large source-to-observation distances ( 400 km), high geographic latitudes, and low geomagnetic field inclinations.

Goyal, H. K.↗

Numerical experiments in geomagnetic modeling

Numerical tests were made, using least squares fitting of a spherical harmonic model, to a selection of Magsat data to determine the practical limits of this technique with modern computers. The resulting (M102189) model, whose coefficients were adjusted up to n = 50, was compared with M07AV6, a previous model which used least squares (on vector data) for coefficients up to n = 29, and Gauss-Legendre quadrature (on Z residuals) to adjust the coefficients up to n = 63. For the new least squares adjustment to n = 50 a condition number of 115 was obtained for the solution matrix, with a resulting precision of 11 significant figures. The M102189 model shows a lower and more Gaussian residual distribution than did M07AV6, though the Gaussian envelope fits to the residual distributions, even for the scalar field, gives "standard deviations' never lower than 6 nT, a factor of three higher than the estimated Magsat observational errors. Ionospheric currents are noted to have a significant effect on the coefficients of the internal potential functions.

Cain, Joseph C.↗