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 397 records · Page 22

Numerical Study of Buoyancy and Different Diffusion Effects on the Structure and Dynamics of Triple Flames

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames. Results of numerical simulations show that gravity strongly affects the triple flame speed through its contribution to the overall flow field. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release, is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed. Downward-propagating triple flames at relatively strong gravity effects have exhibited instabilities. These instabilities are generated without any artificial forcing of the flow. Instead disturbances are initiated by minute round-off errors in the numerical simulations, and subsequently amplified by instabilities. A linear stability analysis on mean profiles of stable triple flame configurations have been performed to identify the most amplified frequency in spatially developed flows. The eigenfunction equations obtained from the linearized disturbance equations are solved using the shooting method. The linear stability analysis yields reasonably good agreements with the observed frequencies of the unstable triple flames. The frequencies and amplitudes of disturbances increase with the magnitude of the gravity vector. Moreover, disturbances appear to be most amplified just downstream of the premixed branches. The effects of mixing width and differential diffusion are investigated and their roles on the flame stability are studied.

Chen, Jyh-Yuan↗

Magnetohydrodynamic equilibrium. II - General integrals of the equations with one ignorable coordinate

The steady equations of hydromagnetics for the isentropic or nonisentropic flow of an inviscid magnetofluid of high electrical conductivity, with one ignorable coordinate in a general orthogonal system, are treated. Several integrals of the equations are established thereafter reducing them to a scalar, quasi-linear, second order, partial differential equation for the magnetic potential. Simple solutions of this final equation are presented. The result, together with a similar treatment of helically symmetric hydromagnetic flows presented in a subsequent paper, allows a unified and systematic approach to the solution of problems involving steady hydromagnetic fields with a topological invariance in various curvilinear coordinates.

Tsinganos, K. C.↗

Highly Accurate Quartic Force Field and Rovibrational Spectroscopic Constants for the Azirinyl Cation (c-C2NH2+) and Its Isomers

Azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high resolution experimental or astronomical observations.

Partha P Bera↗

Numerical simulation of a plane turbulent mixing layer, with applications to isothermal, rapid reactions

A hybrid method has been developed for the numerical prediction of turbulent mixing in a spatially-developing, free shear layer. Most significantly, the computation incorporates the effects of large-scale structures, Schmidt number and Reynolds number on mixing, which have been overlooked in the past. In flow field prediction, large-eddy simulation was conducted by a modified 2-D vortex method with subgrid-scale modeling. The predicted mean velocities, shear layer growth rates, Reynolds stresses, and the RMS of longitudinal velocity fluctuations were found to be in good agreement with experiments, although the lateral velocity fluctuations were overpredicted. In scalar transport, the Monte Carlo method was extended to the simulation of the time-dependent pdf transport equation. For the first time, the mixing frequency in Curl's coalescence/dispersion model was estimated by using Broadwell and Breidenthal's theory of micromixing, which involves Schmidt number, Reynolds number and the local vorticity. Numerical tests were performed for a gaseous case and an aqueous case. Evidence that pure freestream fluids are entrained into the layer by large-scale motions was found in the predicted pdf. Mean concentration profiles were found to be insensitive to Schmidt number, while the unmixedness was higher for higher Schmidt number. Applications were made to mixing layers with isothermal, fast reactions. The predicted difference in product thickness of the two cases was in reasonable quantitative agreement with experimental measurements.

Lin, P.↗

Visual System for Browsing, Analysis, and Retrieval of Data (ViSBARD)

ViSBARD software provides a way of visualizing multiple vector and scalar quantities as measured by many spacecraft at once. The data are displayed three-dimensionally along the orbits that may be shown either as connected lines or as points. The data display allows the rapid determination of vector configurations, correlations among many measurements at multiple points, and global relationships. Things such as vector field rotations and dozens of simultaneous variables are very difficult to see in (complementary) panel plot representations. The current and next generations of space physics missions require a means to display from tens to hundreds of time series of data in such a way that the mind can comprehend them for the purposes of browsing data, retrieving them in directly useful form, and analyzing them in a global context. Sets of many spacecraft, each carrying many instruments yielding nearly continuous data at high time resolution, have become one of the most effective ways to make progress in understanding the extended, ionized (plasma) atmosphere of the Earth and the Sun. For large collections of data to be effective, they must be extremely readily accessible, with simple, comprehensible overviews of what is available. ViSBARD provides a means to answer these concerns. The ViSBARD package also acts as a remote repository browser; an interface to a Virtual Observatory. Therefore, data can be pulled directly into the application, as opposed to searching for it and downloading separately.

Roberts, Aaron↗

Direct numerical simulation of buoyantly driven turbulence

Numerical simulations of homogeneous turbulence subject to buoyant forcing were performed. The presence of a mean temperature gradient combined with a gravitational field results in a forcing term in the momentum equations. The development of the turbulence was studied and compared to the decay of similar fields in the absence of gravity. In the buoyantly driven field, the vorticity is preferentially aligned with the intermediate eigenvector of the strain-rate tensor and the local temperature gradient is more likely to be aligned with the most compressive eigenvector. These relationships are qualitatively similar to those observed in previous shear flow results studied by Ashurst (1987). A tensor diffusivity model for passive scalar transport developed from shear flow results in Rogers, Moin, and Reynolds (1986) also predicts this buoyant scalar transport, indicating that the relationship between the scalar flux and the Reynolds stress is similar in both flows.

Ashurst, William T.↗

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗

Secluded WIMPs, Dark QED with Massive Photons, and the Galactic Center Gamma-Ray Excess

We discuss a particular secluded WIMP dark matter model consisting of neutral fermions as the dark matter candidate and a Proca-Wentzel (PW) field as a mediator. In the model that we consider here, dark matter WIMPs interact with standard model (SM) particles only through the PW field of approximately MeV-multi-GeV mass particles. The interactions occur via a U(1)' mediator, V'(sub mu), which couples to the SM by kinetic mixing with U(1) hypercharge bosons, B'(sub mu). One important difference between our model and other such models in the literature is the absence of an extra singlet scalar, so that the parameter with dimension of mass M(sup 2, sub V) is not related to a spontaneous symmetry breaking. This QED based model is also renormalizable. The mass scale of the mediator and the absence of the singlet scalar can lead to interesting astrophysical signatures. The dominant annihilation channels are different from those usually considered in previous work. We show that the GeV energy gamma-ray excess in the galactic center region, as derived from Fermi-LAT Gamma-ray Space Telescope data, can be attributed to such secluded dark matter WIMPs, given parameters of the model that are consistent with both the cosmological dark matter density and the upper limits on WIMP spin-independent elastic scattering. Secluded WIMP models are also consistent with suggested upper limits on a DM contribution to the cosmic-ray antiproton flux.

Gamma-Rays↗

ISOLEV - A level surface cutting plane program for fluid flow data

A computer program, ISOLEV, allows interactive visualization of CFD scalar and vector functions by means of user-specified sections of the data. The application is based on table lookups which govern both isosurface generation on hexahedral grid cells and recursive subdivision of the cells. The program supports Gouraud-shaded color maps of the data, surface-on-surface maps, and deformation surfaces of vector fields. The execution of the code for animated sweeps is improved by resorting to the cells in the database and maintaining an active set of cells to be rendered.

Kerlick, G. David↗

Assessment and development of second order turbulence models

Second order turbulence models describe the effect of mean flow and external agencies (such as buoyancy) on the evolution of turbulence. Therefore, in principle, these models give a more accurate description of complicated flow fields (e.g. flows with large anisotropy in turbulence, such as near the leading edge of a turbine blade) than the two equation models. The objectives of this document are to assess the performance of the various second order turbulence models in benchmark flows and to seek improvements where necessary in models for the pressure correlation term in the scalar flux equation and for the scalar dissipation equation.

Shabbir, Aamir↗

Thermocapillary migration of a gas bubble in an arbitrary direction with respect to a plane surface

The thermocapillary migration of a gas bubble in an unbounded fluid in the presence of a neighboring rigid plane surface is considered in the limit of negligible Reynolds and Marangoni numbers. Results are given for a scalar interaction parameter defined as the ratio of the speed of the bubble in the presence of the plane surface to the speed in its absence. It is suggested that the weaker interaction effects noted for the case of thermocapillary migration relative to the case of motion due to a body force such as that caused by a gravitational field is attributable to the more rapid decay, away from the bubble, of the disturbance velocity and temperature gradient fields. The surface is found to exert the greatest influence in the case of motion normal to it, and the weakest influence in the case of parallel motion.

Meyyappan, M.↗

The effect of sensor sheltering and averaging techniques on wind measurements at the Shuttle Landing Facility

This document presents results of a field study of the effect of sheltering of wind sensors by nearby foliage on the validity of wind measurements at the Space Shuttle Landing Facility (SLF). Standard measurements are made at one second intervals from 30-feet (9.1-m) towers located 500 feet (152 m) from the SLF centerline. The centerline winds are not exactly the same as those measured by the towers. A companion study, Merceret (1995), quantifies the differences as a function of statistics of the observed winds and distance between the measurements and points of interest. This work examines the effect of nearby foliage on the accuracy of the measurements made by any one sensor, and the effects of averaging on interpretation of the measurements. The field program used logarithmically spaced portable wind towers to measure wind speed and direction over a range of conditions as a function of distance from the obstructing foliage. Appropriate statistics were computed. The results suggest that accurate measurements require foliage be cut back to OFCM standards. Analysis of averaging techniques showed that there is no significant difference between vector and scalar averages. Longer averaging periods reduce measurement error but do not otherwise change the measurement in reasonably steady flow regimes. In rapidly changing conditions, shorter averaging periods may be required to capture trends.

Merceret, Francis J.↗

A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations

Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

Winckelmans, G. S.↗

Unseeded Scalar Velocity Measurements for Propulsion Flows

Unseeded molecular tagging methods based on single-photon processes that produce long tag lines (>50 mm) have been recently developed and demonstrated by the Combustion Laser Diagnostics Group (Mechanical Engineering Department) at Vanderbilt University [1,2]. In Ozone Tagging Velocimetry (OTV) a line of ozone (O3) is produced by a single photon from a pulsed narrowband argon fluoride (ArF) excimer laser operating at - 193 nm. After a known time delay, t, the position of the displaced (convected in the flow field) O3 tag line is revealed by photodissociation of O3 and subsequent fluorescence of O2, caused by a pulsed laser sheet from a krypton fluoride (KrF) excimer laser operating at - 248 nm. Intensified CCD camera images of the fluorescence are taken from the initial and final tag line locations thus providing unobtrusive means of establishing a velocity profile in the interrogated flow field. The O3 lines are "written" and subsequently "read" by the following reactions:

Pitz, Robert W.↗

Enhanced damping of Alfven waves in the solar corona by a turbulent wave spectrum

The effect of a background spectrum of Alfven waves on the rate of dissipation of a test shear Alfven wave is numerically calculated. The results demonstrate that as the classical resistivity eta and classical viscosity mu become small, the damping rate of the Alfven wave remains large and depends only on the amplitude for the scalar potential of the wave spectrum and the wavenumber of the Alfven wave. The damping rate is virtually independent of eta and mu. The wave spectrum need not be turbulent or stochastic to affect the damping rate. The dissipation rate is nonlinear enhanced by nonstochastic spectra as well as by stochastic spectra if two conditions are met. First, the perpendicular magnetic field associated with Alfven wave spectrum must exceed a certain collision-frequency threshold and second, for nonstochastic spectra only, the magnetic field must exceed a threshold proportional to the parallel wavenumber of the shear Alfven wave. These conditions can be easily satisfied in the solar corona.

Kleva, Robert G.↗

Errors incurred in a plane-wave-type expansion of a Gaussian beam

The multipole expansion obtained by Morita et al. (1968) of the Gaussian laser beam used to levitate an aerosol particle in order that its complete phase matrix may be measured is compared with that of Tsai and Pogorzelski (1975) in order to demonstrate the effect of the incorrect expansion used by Morita. Errors incurred by the use of an equation in which one side satisfies the scalar wave equation while the other side does not and can be reduced to a plane wave amplitude are calculated as functions of the inverse of the wave number times the beam waist, the wave number times the radial spherical coordinate and the angular spherical coordinate. Errors on the order of a few percent, considered undetectable are obtained in the squared-field amplitudes due to the expansion, however, they are found to become significant (several tens of percent) when the angle is zero. It is concluded that the expansion of Morita should only be used in the regions where the spherical angle is less than 0.01 and its product with the wave number and the radial spherical coordinate is less than unity.

Kattawar, G. W.↗

The CMC/3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 3: Programmers' manual

The CMC fluid mechanics program system was developed to transmit the theoretical evolution of finite element numerical solution methodology, applied to nonlinear field problems into a versatile computer code for comprehensive flow field analysis. A detailed view of the code from the standpoint of a computer programmer's use is presented. A system macroflow chart and detailed flow charts of several routines necessary to interact with a theoretican/user to modify the operation of this program are presented. All subroutines and details of usage, primarily for input and output routines are described. Integer and real scalars and a cross reference list denoting subroutine usage for these scalars are outlined. Entry points in dynamic storage vector IZ; the lengths of each vector accompanying the scalar definitions are described. A listing of the routines peculiar to the standard test case and a listing of the input deck and printout for this case are included.

Orzechowski, J. A.↗

On the modelling of non-reactive and reactive turbulent combustor flows

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Nikjooy, Mohammad↗