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 55 records · Page 3

Do massive Brans-Dicke theories of gravitation imitate Brans-Dicke theories with nonzero divergence of energy-momentum tensor

Acharya and Hogan (1973) have introduced a massive scalar field into the usual Brans-Dicke (1961) theory of gravitation. Formally they obtain certain field equations. The assumption of a nonzero divergence for T mu nu (or equivalently the introduction of sources) formally imitates a massive Brans-Dicke and satisfies the condition of Acharya and Hogan that the theory be indistinguishable with the classical test of the Einstein theory. Although it was shown elsewhere that the modified Brans-Dicke theory agrees with the classical test under certain conditions, there were no specified limits on omega, a similar circumstance discovered by Acharya and Hogan for the massive scalar field.

Smalley, L. L.↗

Thermal relics: Do we know their abundances

The relic abundance of a particle species that was once in thermal equilibrium in the expanding Universe depends upon a competition between the annihilation rate of the species and the expansion rate of the Universe. Assuming that the Universe is radiation dominated at early times the relic abundance is easy to compute and well known. At times earlier than about 1 sec after the bang there is little or no evidence that the Universe had to be radiation dominated, although that is the simplest and standard assumption. Because early-Universe relics are of such importance both to particle physics and to cosmology, three nonstandard possibilities are considered in detail for the Universe at the time a species' abundance froze in: energy density dominated by shear (i.e., anisotropic expansion), energy density dominated by some other nonrelativistic species, and energy density dominated by the kinetic energy of the scalar field that sets the gravitational constant in a Brans-Dicke-Jordan cosmological mode. In the second case the relic abundance is less than the standard value, while in the other two cases it can be enhanced by a significant factor. Two other more exotic possibilities for enhancing the relic abundance of a species are also mentioned--a larger value of Newton's constant at early times (e.g., as might occur in superstring or Kaluza-Klein theories) or a component of the energy density at early times with a very stiff equation of state (p greater than rho/3), e.g., a scalar field phi with potential V(phi) = Beta /phi/ (exp n) with n greater than 4. Results have implications for dark matter searches and searches for particle relics in general.

Kamionkowski, Marc↗

On thick domain walls in general relativity

Planar scalar field configurations in general relativity differ considerably from those in flat space. It is shown that static domain walls of finite thickness in curved space-time do not possess a reflection symmetry. At infinity, the space-time tends to the Taub vacuum on one side of the wall and to the Minkowski vacuum (Rindler space-time) on the other. Massive test particles are always accelerated towards the Minkowski side, i.e., domain walls are attractive on the Taub side, but repulsive on the Minkowski side (Taub-vacuum cleaner). It is also proved that the pressure in all directions is always negative. Finally, a brief comment is made concerning the possibility of infinite, i.e., bigger than horizon size, domain walls in our universe. All of the results are independent of the form of the potential V(phi) greater than or equal to 0 of the scalar field phi.

Goetz, Guenter↗

A transport equation for the scalar dissipation in reacting flows with variable density: First results

Although the different regimes of premixed combustion are not well defined, most of the recent developments in turbulent combustion modeling are led in the so-called flamelet regime. The goal of these models is to give a realistic expression to the mean reaction rate (w). Several methods can be used to estimate (w). Bray and coworkers (Libby & Bray 1980, Bray 1985, Bray & Libby 1986) express the instantaneous reaction rate by means of a flamelet library and a frequency which describes the local interaction between the laminar flamelets and the turbulent flowfield. In another way, the mean reaction rate can be directly connected to the flame surface density (Sigma). This quantity can be given by the transport equation of the coherent flame model initially proposed by Marble & Broadwell 1977 and developed elsewhere. The mean reaction rate, (w), can also be estimated thanks to the evolution of an arbitrary scalar field G(x, t) = G(sub O) which represents the flame sheet. G(x, t) is obtained from the G-equation proposed by Williams 1985, Kerstein et al. 1988 and Peters 1993. Another possibility proposed in a recent study by Mantel & Borghi 1991, where a transport equation for the mean dissipation rate (epsilon(sub c)) of the progress variable c is used to determine (w). In their model, Mantel & Borghi 1991 considered a medium with constant density and constant diffusivity in the determination of the transport equation for (epsilon(sub c)). A comparison of different flamelet models made by Duclos et al. 1993 shows the realistic behavior of this model even in the case of constant density. Our objective in this present report is to present preliminary results on the study of this equation in the case of variable density and variable diffusivity. Assumptions of constant pressure and a Lewis number equal to unity allow us to significantly simplify the equation. A systematic order of magnitude analysis based on adequate scale relations is performed on each term of the equation. As in the case of constant density and constant diffusivity, the effects of stretching of the scalar field by the turbulent strain field, of local curvature, and of chemical reactions are predominant. In this preliminary work, we suggest closure models for certain terms, which will be validated after comparisons with DNS data.

Mantel, T.↗

An algebraic model for the turbulent flux of a passive scalar

Direct numerical simulations of the unsteady incompressible Navier-Stokes equations have been performed in order to investigate the behavior of passive-scalar fields resulting from mean scalar gradients in each of three orthogonal directions in homogeneous turbulent shear flow. For all orientations of the mean scalar gradient, the sum of the pressure-scalar gradient and velocity gradient-scalar gradient terms in the turbulent scalar flux balance equation are found to be approximately aligned with the scalar flux vector itself. The simulation results are used to obtain dimensionless model coefficients as a function of the turbulence Reynolds and Peclet numbers.

Rogers, Michael M.↗

Near-field structures in three-dimensional spatially-developing wakes

The effects of three-dimensional perturbations on a spatially-developing plane Gaussian wake are studied and quantified using a new numerical code. Results are presented for a elatively simple case since it appears that no 3-D spatially-developing simulations of plane wakes have been performed in the past. The inlet perturbation consists of a 2-D time-periodic forcing plus a spanwise-periodic array of steady streamwise vortices. The distortion of the 2-D vortices caused by the streamwise vortices will be discussed. The vorticity fields are compared to passive scalar fields which simulate the injection of smoke or low heat into the wake.

Buell, Jeffrey C.↗

Stochastic modeling of turbulent reacting flows

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

Fox, R. O.↗

Higher-order derivative correlations and the alignment of small-scale structures in isotropic numerical turbulence

In a three dimensional simulation higher order derivative correlations, including skewness and flatness factors, are calculated for velocity and passive scalar fields and are compared with structures in the flow. The equations are forced to maintain steady state turbulence and collect statistics. It is found that the scalar derivative flatness increases much faster with Reynolds number than the velocity derivative flatness, and the velocity and mixed derivative skewness do not increase with Reynolds number. Separate exponents are found for the various fourth order velocity derivative correlations, with the vorticity flatness exponent the largest. Three dimensional graphics show strong alignment between the vorticity, rate of strain, and scalar-gradient fields. The vorticity is concentrated in tubes with the scalar gradient and the largest principal rate of strain aligned perpendicular to the tubes. Velocity spectra, in Kolmogorov variables, collapse to a single curve and a short minus 5/3 spectral regime is observed.

Kerr, R. A.↗

Turbulence spectrum of strongly conductive temperature field in a rapidly stirred fluid

A simulation is performed of a passive scalar field convected by a rapidly fluctuating velocity field whose correlation time approaches zero. By using a code proposed in a previous study (Chasnov et al., 1988), the turbulence spectrum of the passive temperature field in the conductive subrange is determined. A theoretical model is proposed which explains the result obtained by representing the transfer of scalar variance by an eddy conductivity, whose correlation time is limited by the correlation time of the velocity field.

Chasnov, J.↗

Initial geomagnetic field model from MAGSAT

Magsat data from magnetically quiet days were used to derive a thirteenth degree and order spherical harmonic geomagnetic field model, MGST(3/80). The model utilized both scalar and vector data and fit that data with standard deviations of 8, 52, 55 and 97 nT for the scalar magnitude, B sub r, B sub theta and B sub phi respectively. When compared with earlier models, the Earth's dipole moment continues to decrease at a rate of about 26 nT/year. Evaluation of earlier models with Magsat data shows that the scalar field at the Magsat epoch is best predicted by the POGO(2/72) model but that the AWC/75 and IGS/75 are better for predicting vector fields.

Langel, R. A.↗

Development of Techniques for Visualization of Scalar and Vector Fields in the Immersive Environment

Visualization of scalar and vector fields in the immersive environment (CAVE - Cave Automated Virtual Environment) is important for its application to radiation shielding research at NASA Langley Research Center. A complete methodology and the underlying software for this purpose have been developed. The developed software has been put to use for the visualization of the earth s magnetic field, and in particular for the study of the South Atlantic Anomaly. The methodology has also been put to use for the visualization of geomagnetically trapped protons and electrons within Earth's magnetosphere.

Bidasaria, Hari B.↗

Electromagnetic and scalar diffraction by a right-angled wedge with a uniform surface impedance

The diffraction of an electromagnetic wave by a perfectly-conducting right-angled wedge with one surface covered by a dielectric slab or absorber is considered. The effect of the coated surface is approximated by a uniform surface impedance. The solution of the normally incident electromagnetic problem is facilitated by introducing two scalar fields which satisfy a mixed boundary condition on one surface of the wedge and a Neumann of Dirichlet boundary condition on the other. A functional transformation is employed to simplify the boundary conditions so that eigenfunction expansions can be obtained for the resulting Green's functions. The eigenfunction expansions are transformed into the integral representations which then are evaluated asymptotically by the modified Pauli-Clemmow method of steepest descent. A far zone approximation is made to obtain the scattered field from which the diffraction coefficient is found for scalar plane, cylindrical or sperical wave incident on the edge. With the introduction of a ray-fixed coordinate system, the dyadic diffraction coefficient for plane or cylindrical EM waves normally indicent on the edge is reduced to the sum of two dyads which can be written alternatively as a 2 X 2 diagonal matrix.

Hwang, Y. M.↗

The reduction, verification and interpretation of MAGSAT magnetic data over Canada

Procedures were established to select quiet, least-disturbed MAGSAT data in an effort to study the magnetic field originating in the solid Earth. A preliminary map of scalar field anomalies over Canada and adjacent regions was produced. Data from Fort Churchill magnetic observatory (an auroral zone station) were used in classifying MAGSAT passes as quiet, unsettled, or disturbed. Determinations of long-term quiet levels at Canadian magnetic observatories are being made. For the time interval of each MAGSAT pass, the deviations of the mean hourly values of the magnetic field components are determined from these quiet levels, in particular for Fort Churchill. Primary processing all CHRONINT tapes received to date, involved decoding orbit records, decimating data records, selecting time intervals and latitude intervals, and producing output tapes carrying records with data point time, position and magnetic information. Subsequent processing converts the coordinate system from geocentric to geodetic and removes a reference field.

Coles, R. L.↗

Three-dimensional singular points in aerodynamics

When three-dimensional separation occurs on a body immersed in a flow governed by the incompressible Navier-Stokes equations, the geometrical surfaces formed by the three vector fields (velocity, vorticity and the skin-friction) and a scalar field (pressure) become interrelated through topological maps containing their respective singular points and extremal points. A mathematically consistent description of these singular points becomes inevitable when we want to study the geometry of the separation. A separated stream surface requires, for example, the existence of a saddle-type singular point on the skin-friction surface. This singular point is actually, in the proper language of mathematics, a saddle of index two. The index is a measure of the dimension of the outset (set leaving the singular point). Hence, when a saddle of index two is specified, a two dimensional surface that becomes separated from the osculating plane of the saddle is implied. The three-dimensional singular point is interpreted mathematically and the most common aerodynamical singular points are discussed through this perspective.

Unal, Aynur↗

Field Encapsulation Library The FEL 2.2 User Guide

This document describes version 2.2 of the Field Encapsulation Library (FEL), a library of mesh and field classes. FEL is a library for programmers - it is a "building block" enabling the rapid development of applications by a user. Since FEL is a library intended for code development, it is essential that enough technical detail be provided so that one can make full use of the code. Providing such detail requires some assumptions with respect to the reader's familiarity with the library implementation language, C++, particularly C++ with templates. We have done our best to make the explanations accessible to those who may not be completely C++ literate. Nevertheless, familiarity with the language will certainly help one's understanding of how and why things work the way they do. One consolation is that the level of understanding essential for using the library is significantly less than the level that one should have in order to modify or extend the library. One more remark on C++ templates: Templates have been a source of both joy and frustration for us. The frustration stems from the lack of mature or complete implementations that one has to work with. Template problems rear their ugly head particularly when porting. When porting C code, successfully compiling to a set of object files typically means that one is almost done. With templated C++ and the current state of the compilers and linkers, generating the object files is often only the beginning of the fun. On the other hand, templates are quite powerful. Used judiciously, templates enable more succinct designs and more efficient code. Templates also help with code maintenance. Designers can avoid creating objects that are the same in many respects, but not exactly the same. For example, FEL fields are templated by node type, thus the code for scalar fields and vector fields is shared. Furthermore, node type templating allows the library user to instantiate fields with data types not provided by the FEL authors. This type of flexibility would be difficult to offer without the support of the language. For users who may be having template-related problems, we offer the consolation that support for C++ templates is destined to improve with time. Efforts such as the Standard Template Library (STL) will inevitably drive vendors to provide more thorough, optimized tools for template code development. Furthermore, the benefits will become harder to resist for those who currently subscribe to the least-common-denominator "code it all in C" strategy. May FEL bring you both increased productivity and aesthetic satisfaction.

Moran, Patrick J.↗

Numerical simulation of low Prandtl number turbulent mixing

Numerical simulations of turbulent mixing of strongly diffusive scalar fields were carried out with and without subgrid-scale modeling of the small-scale strain field. For low Reynolds number flows, when the rate of strain field (determined primarily by the small scales) is fully resolved, the scalar microstructure was found to collapse under Batchelor rate-of-strain scaling even for small Prandtl numbers, in agreement with Kerr. For high Reynolds number flows, when small-scale straining is modeled with a subgrid-scale model, the scalar microstructure follows the Batchelor, Howells, and Townsend predictions that the small-scale rate-of-strain is irrelevant.

Gibson, C.↗

Vorticity and divergence in the solar photosphere

We have studied an outstanding sequence of continuum images of the solar granulation from Pic du Midi Observatory. We have calculated the horizontal vector flow field using a correlation tracking algorithm, and from this determined three scalar field: the vertical component of the curl; the horizontal divergence; and the horizontal flow speed. The divergence field has substantially longer coherence time and more power than does the curl field. Statistically, curl is better correlated with regions of negative divergence - that is, the vertical vorticity is higher in downflow regions, suggesting excess vorticity in intergranular lanes. The average value of the divergence is largest (i.e., outflow is largest) where the horizontal speed is large; we associate these regions with exploding granules. A numerical simulation of general convection also shows similar statistical differences between curl and divergence. Some individual small bright points in the granulation pattern show large local vorticities.

Wang, YI↗