Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scalar fields”

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 379 records · Page 21

Near-wall similarity in three-dimensional turbulent boundary layers. I - Model review. II - Pressure-driven flow results

Eleven proposed models for near-wall similarity for three-dimensional turbulent boundary layer flows are reviewed. Six of these models are comparatively simple scalar models and five are more complex and/or two-component vector models. Ten of the models can be tested as to their validity or predictive capability with the aid of measured mean velocity field, wall pressure field, and direct wall shear stress field (magnitude and direction) data. One of the models cannot be tested owing to its dependence on two parameters that are at present extremely difficult (if not impossible) to measure. Ten three-dimensional near-wall similarity models are then evaluated with direct wall shear, velocity field, and pressure gradient data from a three-dimensional pressure-driven boundary layer flow. In a primary focus of the interval where y+ is between 50 and 300, graphical results suggest that six simpler models and the freestream component of one complex model are adequate for profiles with monotone increasing skew up to about 15 deg.

Pierce, F. J.↗

On the Alignment of Strain, Vorticity and Scalar Gradient in Turbulent, Buoyant, Nonpremixed Flames

The alignment of vorticity and scalar gradient with the eigendirections of the rate of strain tensor is investigated in turbulent buoyant nonpremixed horizontal and vertical flames. The uniqueness of a buoyant nonpremixed flame is that it contains regions with distinct alignment characteristics. The strain-enstrophy angle Psi is used to identify these regions. Examination of the vorticity field and the vorticity production in these different regions indicates that Psi and consequently the alignment properties near the flame surface identified by the mixture fraction band F approximately equals F(sub st) differ from those in the fuel region, F > F(sub st) and the oxidizer region, F < F(sub st). The F approximately equals F(sub st) band shows strain-dominance resulting in vorticity/alpha alignment while F > F(sub st) (and F < F(sub st) for the vertical flame) band(s) show(s) vorticity/beta alignment. The implication of this result is that the scalar dissipation, epsilon(sub F), attains its maximum value always near F approximately equals F(sub st). These results are also discussed within the framework of recent dynamical results [Galanti et al., Nonlinearity 10, 1675 (1997)] suggesting that the Navier-Stokes equations evolved towards an attracting solution. It is shown that the properties of such an attracting solution are also consistent with our results of buoyant turbulent nonpremixed flames.

Boratav, O. N.↗

Spaceborne magnetometry

Low orbit (160 km altitude) satellite magnetometry is discussed. Magnetic field simulations and theoretical analyses were performed. They suggest that at 160 km, resolution value is 10, compared to 0.6 for Magsat, and that the satellite can reveal significant information on crustal tectonics. Vector and scalar magnetometers should be carried.

Taylor, P. T.↗

Transmitter and receiver antenna gain analysis for laser radar and communication systems

A comprehensive and fairly self-contained study of centrally obscured optical transmitting and receiving antennas is presented and is intended for use by the laser radar and communication systems designer. The material is presented in a format which allows the rapid and accurate evaluation of antenna gain. The Fresnel approximation to scalar wave theory is reviewed and the antenna analysis proceeds in terms of the power gain. Conventional range equations may then be used to calculate the power budget. The transmitter calculations, resulting in near and far field antenna gain patterns, assumes the antenna is illuminated by a laser operating in the fundamental cavity mode. A simple equation is derived for matching the incident source distribution to a general antenna configuration for maximum on-axis gain. An interpretation of the resultant gain curves allows a number of auxiliary design curves to be drawn which display the losses in antenna gain due to pointing errors and the cone angle of the outgoing beam as a function of antenna size and central obscuration. The use of telescope defocusing as an approach to spreading the beam for target acquisition is compared to some alternate methods.

Klein, B. J.↗

Vectorization of a penalty function algorithm for well scheduling

In petroleum engineering, the oil production profiles of a reservoir can be simulated by using a finite gridded model. This profile is affected by the number and choice of wells which in turn is a result of various production limits and constraints including, for example, the economic minimum well spacing, the number of drilling rigs available and the time required to drill and complete a well. After a well is available it may be shut in because of excessive water or gas productions. In order to optimize the field performance a penalty function algorithm was developed for scheduling wells. For an example with some 343 wells and 15 different constraints, the scheduling routine vectorized for the CYBER 205 averaged 560 times faster performance than the scalar version.

Absar, I.↗

Eddy viscosity measurements in a rectangular jet

The flow field of a rectangular jet with a 2:1 aspect ratio was studied at a Reynolds number of 100,000 (Mach number 0.09) using three-dimensional laser Doppler velocimetry (LDV). Velocity gradients, Reynolds stress tensor components, and scalar eddy viscosities are presented for the major and minor axis planes of the jet. The eddy viscosity model was found to be applicable only in the direction of maximum mean velocity gradient.

Swan, David H.↗

Modeling Sound Propagation Through Non-Axisymmetric Jets

A method for computing the far-field adjoint Green's function of the generalized acoustic analogy equations under a locally parallel mean flow approximation is presented. The method is based on expanding the mean-flow-dependent coefficients in the governing equation and the scalar Green's function in truncated Fourier series in the azimuthal direction and a finite difference approximation in the radial direction in circular cylindrical coordinates. The combined spectral/finite difference method yields a highly banded system of algebraic equations that can be efficiently solved using a standard sparse system solver. The method is applied to test cases, with mean flow specified by analytical functions, corresponding to two noise reduction concepts of current interest: the offset jet and the fluid shield. Sample results for the Green's function are given for these two test cases and recommendations made as to the use of the method as part of a RANS-based jet noise prediction code.

Leib, Stewart J.↗

Integrals of motion for the classical two-body problem with drag

Integrals of motion for the two-body problem with drag are obtained by operating on the second-order vector differential equation describing the motion. The force field consists of an inverse-square gravitational attraction and a drag force proportional to the velocity vector and inversely proportional to the square of the distance to the attracting center. The developed integrals are the analogs of the Keplerian scalar energy, the vector angular momentum, and the Laplace vector.

Jezewski, D. J.↗

Relativistic astrophysics

The coordinate independence of gravitational radiation and the parameterized post-Newtonian approximation from which it is extended are described. The general consistency of the field equations with Bianchi identities, gauge conditions, and the Newtonian limit of the perfect fluid equations of hydrodynamics are studied. A technique of modification is indicated for application to vector-metric or double metric theories, as well as to scalar-tensor theories.

Smalley, L. L.↗

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↗

Three-Dimensional Multiscale MHD Model of Cometary Plasma Environments

First results of a three-dimensional multiscale MHD model of the interaction of an expanding cometary atmosphere with the magnetized solar wind are presented. The model starts with a supersonic and super-Alfvenic solar wind far upstream of the comet (25 Gm upstream of the nucleus) with arbitrary interplanetary magnetic field orientation. The solar wind is continuously mass loaded with cometary ions originating from a 10-km size nucleus. The effects of photoionization, electron impact ionization, recombination, and ion-neutral frictional drag are taken into account in the model. The governing equations are solved on an adaptively refined unstructured Cartesian grid using our new multiscale upwind scalar conservation laws-type numerical technique (MUSCL). We have named this the multiscale adaptive upwind scheme for MHD (MAUS-MHD). The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the diamagnetic cavity of the comet. The main findings are the following: (1) Mass loading decelerates the solar wind flow upstream of the weak cometary shock wave (M approximately equals 2, M(sub A) approximately equals 2), which forms at a subsolar standoff distance of about 0.35 Gm. (2) A cometary plasma cavity is formed at around 3 x 10(exp 3) km from the nucleus. Inside this cavity the plasma expands outward due to the frictional interaction between ions and neutrals. On the nightside this plasma cavity considerably narrows and a relatively fast and dense cometary plasma beam is ejected into the tail. (3) Inside the plasma cavity a teardrop-shaped inner shock is formed, which is terminated by a Mach disk on the nightside. Only the region inside the inner shock is the 'true' diamagnetic cavity. (4) The model predicts four distinct current systems in the inner coma: the density peak current, the cavity boundary current, the inner shock current, and finally the cross-tail current. (5) The calculated plasma parameters (magnetic field, plasma density, speed, and temperature) are in very good agreement with published Giotto observations.

Gombosi, Tamas I.↗

An implicit marching procedure for the treatment of supersonic flow fields using the conservative full potential equation

A new aerodynamic prediction technique based on the conservative full potential equation is developed for the treatment of supersonic flow fields. This new technique bridges the gap between simplistic linear theory methods and complex Euler solvers. A novel local density linearization concept and a second order accurate retarded density scheme, both producing the correct artificial viscosity, are introduced in developing an implicit marching scheme for solving the scalar phi. The method produces results that compare well with Euler solvers and requires an order of magnitude less computer time and significantly less computer memory over existing nonlinear codes. The scalar phi formulation can be extended to handle subsonic pockets in the marching direction and also is suitable for developing inverse procedures where the shape corresponding to a prescribed loading is sought.

Shankar, V.↗

Beyond the Kirchhoff approximation. II - Electromagnetic scattering

In a paper by Rodriguez (1981), the momentum transfer expansion was introduced for scalar wave scattering. It was shown that this expansion can be used to obtain wavelength-dependent curvature corrections to the Kirchhoff approximation. This paper extends the momentum transfer perturbation expansion to electromagnetic waves. Curvature corrections to the surface current are obtained. Using these results, the specular field and the backscatter cross section are calculated.

Rodriguez, Ernesto↗

Fractal interfaces and product generation in the two-dimensional mixing layer

The dependence of product generation on Peclet and Reynolds numbers in a numerically simulated, reacting, two-dimensional, temporally growing mixing layer is related theoretically to the fractal dimension of the passive scalar interfaces. This reaction is verified using product generation measurements and fractal dimensions derived from the box counting technique. A transition from a low initial dimension to a higher one of approximately 5/3 is identified and shown to be associated to the kinematic distortion of the flow field during the first pairing interaction. It is suggested that the structures reponsible for this transition are nondeterministic, nonrandom, inhomogeneous fractals. In the range of Schmidt numbers investigated (0.25-4), only the large scales are involved. No further transitions, either in the spectra of the vorticity field or in the mixing behavior, are found for Reynolds numbers up to 90,000.

Jimenez, Javier↗

Fractal interfaces and product generation in the two dimensional mixing layer

The dependence of product generation on Peclet and Reynolds numbers in a numerically simulated, reacting, two dimensional, temporally growing mixing layer is related theoretically to the fractal dimension of the passive scalar interfaces. This relation is verified using product generation measurements and dimensions derived from a standard box counting technique. A transition from a low initial dimension to a higher one of approximately 5/3 is identified and shown to be associated to the kinematic distortion on the flow field during the first pairing interaction. It is suggested that the structures responsible for this transition are non-deterministic, non-random, inhomogeneous fractals. Only the large scales are involved. No further transitions, either in the spectra of the vorticity field or in the mixing behavior, are found for Reynolds numbers up to 90,000.

Jimenez, Javier↗

Variance of phase fluctuations of waves propagating through a random medium

As an electromagnetic wave propagates through a random scattering medium, such as a forest, its energy is attenuated and random phase fluctuations are induced. The magnitude of the random phase fluctuations induced is important in estimating how well a Synthetic Aperture Radar (SAR) can image objects within the scattering medium. The two-layer random medium model, consisting of a scattering layer between free space and ground, is used to calculate the variance of the phase fluctuations induced between a transmitter located above the random medium and a receiver located below the random medium. The scattering properties of the random medium are characterized by a correlation function of the random permittivity fluctuations. The effective permittivity of the random medium is first calculated using the strong fluctuation theory, which accounts for large permittivity fluctuations of the scatterers. The distorted Born approximation is used to calculate the first-order scattered field. A perturbation series for the phase of the received field in the Rytov approximation is then introduced and the variance of the phase fluctuations is also calculated assuming that the transmitter and receiver are in the paraxial limit of the random medium, which allows an analytic solution to be obtained. Results are compared using the paraxial approximation, scalar Green's function formulation, and dyadic Green's function formulation. The effects studied are the dependence of the variance of the phase fluctuations on receiver location in lossy and lossless regions, medium thickness, correlation length and fractional volume of scatterers, depolarization of the incident wave, ground layer permittivity, angle of incidence, and polarization.

Chu, Nelson C.↗

Toward Immersed Boundary Simulation of High Reynolds Number Flows

In the immersed boundary (IB) method, the surface of an object is reconstructed with forcing terms in the underlying flow field equations. The surface may split a computational cell removing the constraint of the near wall gridlines to be aligned with the surface. This feature greatly simplifies the grid generation process which is cumbersome and expensive in particular for structured grids and complex geometries. The IB method is ideally suited for Cartesian flow solvers. The flow equations written in Cartesian coordinates appear in a very simple form and several numerical algorithms can be used for an efficient solution of the equations. In addition, the accuracy of numerical algorithms is dependent on the underlying grid and it usually deteriorates when the grid deviates from a Cartesian mesh. The challenge for the IB method lies in the representation of the wall boundaries and in providing an adequate near wall flow field resolution. The issue of enforcing no-slip boundary conditions at the immersed surface has been addressed by several authors by imposing a local reconstruction of the solution. Initial work by Verzicco et al. was based on a simple linear, one-dimensional operator and this approach proved to be accurate for boundaries largely aligned with the grid lines. Majumdar et al. used various multidimensional and high order polynomial interpolations schemes. These high order schemes, however, are keen to introduce wiggles and spurious extrema. Iaccarino & Verzicco and Kalitzin & Iaccarino proposed a tri-linear reconstruction for the velocity components and the turbulent scalars. A modified implementation that has proven to be more robust is reported in this paper. The issue of adequate near wall resolution in a Cartesian framework can initially be addressed by using a non-uniform mesh which is stretched near the surface. In this paper, we investigate an unstructured approach for local grid refinement that utilizes Cartesian mesh features. The computation of high Reynolds number wall bounded flows is particularly challenging as it requires the consideration of thin turbulent boundary layers, i.e. near wall regions with large gradients of the flow field variables. For such flows, the representation of the wall boundary has a large impact on the accuracy of the computation. It is also critical for the robustness and convergence of the flow solver.

Kalitzin, Georgi↗

Electromagnetic scattering theory

Electromagnetic scattering theory is discussed with emphasis on the general stochastic variational principle (SVP) and its applications. The stochastic version of the Schwinger-type variational principle is presented, and explicit expressions for its integrals are considered. Results are summarized for scalar wave scattering from a classic rough-surface model and for vector wave scattering from a random dielectric-body model. Also considered are the selection of trial functions and the variational improvement of the Kirchhoff short-wave approximation appropriate to large size-parameters. Other applications of vector field theory discussed include a general vision theory and the analysis of hydromagnetism induced by ocean motion across the geomagnetic field. Levitational force-torque in the magnetic suspension of the disturbance compensation system (DISCOS), now deployed in NOVA satellites, is also analyzed using the developed theory.

Bird, J. F.↗