SEARCH · Engineering Papers

Results for “SCALARS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Cosmic microwave background probes models of inflation

Inflation creates both scalar (density) and tensor (gravity wave) metric perturbations. We find that the tensor-mode contribution to the cosmic microwave background anisotropy on large-angular scales can only exceed that of the scalar mode in models where the spectrum of perturbations deviates significantly from scale invariance. If the tensor mode dominates at large-angular scales, then the value of DeltaT/T predicted on 1 deg is less than if the scalar mode dominates, and, for cold-dark-matter models, bias factors greater than 1 can be made consistent with Cosmic Background Explorer (COBE) DMR results.

Davis, Richard L.↗

NASA NTRS1993

Dynamic localization and second-order subgrid-scale models in large eddy simulations of channel flow

The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Cabot, William H.↗

NASA NTRS1994

Gain weighted eigenspace assignment

This report presents the development of the gain weighted eigenspace assignment methodology. This provides a designer with a systematic methodology for trading off eigenvector placement versus gain magnitudes, while still maintaining desired closed-loop eigenvalue locations. This is accomplished by forming a cost function composed of a scalar measure of error between desired and achievable eigenvectors and a scalar measure of gain magnitude, determining analytical expressions for the gradients, and solving for the optimal solution by numerical iteration. For this development the scalar measure of gain magnitude is chosen to be a weighted sum of the squares of all the individual elements of the feedback gain matrix. An example is presented to demonstrate the method. In this example, solutions yielding achievable eigenvectors close to the desired eigenvectors are obtained with significant reductions in gain magnitude compared to a solution obtained using a previously developed eigenspace (eigenstructure) assignment method.

Davidson, John B.↗

NASA NTRS1994

Analysis and design of numerical schemes for gas dynamics 1: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence

The theory of non-oscillatory scalar schemes is developed in this paper in terms of the local extremum diminishing (LED) principle that maxima should not increase and minima should not decrease. This principle can be used for multi-dimensional problems on both structured and unstructured meshes, while it is equivalent to the total variation diminishing (TVD) principle for one-dimensional problems. A new formulation of symmetric limited positive (SLIP) schemes is presented, which can be generalized to produce schemes with arbitrary high order of accuracy in regions where the solution contains no extrema, and which can also be implemented on multi-dimensional unstructured meshes. Systems of equations lead to waves traveling with distinct speeds and possibly in opposite directions. Alternative treatments using characteristic splitting and scalar diffusive fluxes are examined, together with modification of the scalar diffusion through the addition of pressure differences to the momentum equations to produce full upwinding in supersonic flow. This convective upwind and split pressure (CUSP) scheme exhibits very rapid convergence in multigrid calculations of transonic flow, and provides excellent shock resolution at very high Mach numbers.

Jameson, Antony↗

NASA NTRS1995

Optical implementation of inner product neural associative memory

An optical implementation of an inner-product neural associative memory is realized with a first spatial light modulator for entering an initial two-dimensional N-tuple vector and for entering a thresholded output vector image after each iteration until convergence is reached, and a second spatial light modulator for entering M weighted vectors of inner-product scalars multiplied with each of the M stored vectors, where the inner-product scalars are produced by multiplication of the initial input vector in the first iterative cycle (and thresholded vectors in subsequent iterative cycles) with each of the M stored vectors, and the weighted vectors are produced by multiplication of the scalars with corresponding ones of the stored vectors. A Hughes liquid crystal light valve is used for the dual function of summing the weighted vectors and thresholding the sum vector. The thresholded vector is then entered through the first spatial light modulator for reiteration of the process cycle until convergence is reached.

Liu, Hua-Kuang↗

NASA NTRS1995

Turbulence modeling for non-equilibrium flow

The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Durbin, P. A.↗

NASA NTRS1996

Analysis of Subgrid Eddy Diffusivity in Turbulence Using Direct Numerical Simulation Databases

Direct numerical simulation data in isotropic turbulence at Taylor-scale Reynolds number 90 are used to compute subgrid transfer and spectral eddy diffusivities for passive scalars at Schmidt numbers 1 and 1/8. The results are consistent with spectral transfer processes dominated by a local forward cascade in wavenumber space which is a consequence of moderately nonlocal interactions in the scalar fields. Data decompositions show that resolvable-subgrid interactions make (in general) substantially stronger contributions than fully-subgrid interactions to the eddy diffusivity, which has a non-negligible inverse cascade, or backscatter, component. The eddy diffusivity has a strong cusp at wavenumbers close to the cutoff value separating resolvable and subgrid ranges. At low wavenumbers the eddy diffusivity is generally small, but there the details depend on Schmidt number. The behavior of the subgrid turbulent Schmidt number (the ratio of effective eddy viscosity to eddy diffusivity) is observed to depend strongly on Schmidt number. In interpreting the results at a given spectral cutoff it is important to note that the spectrum of less diffusive scalars (of higher Schmidt number) has more high-wavenumber content.

Yeung, P. K.↗

NASA NTRS1996

Rayleigh Light Scattering for Concentration Measurements in Turbulent Flows

Despite intensive research over a number of years, an understanding of scalar mixing in turbulent flows remains elusive. An understanding is required because turbulent mixing has a pivotal role in a wide variety of natural and technologically important processes. As an example, the mixing and transport of pollutants in the atmosphere and in bodies of water are often dependent on turbulent mixing processes. Turbulent mixing is also central to turbulent combustion which underlies most hydrocarbon energy use in modern societies as well as in unwanted fire behavior. Development of models for combusting flows is therefore crucial, however, an understanding of scalar mixing is required before useful models of turbulent mixing and, ultimately, turbulent combustion can be developed. An important subset of turbulent flows is axisymmetric turbulent jets and plumes because they are relatively simple to generate, and because the provide an appropriate test bed for the development of general theories of turbulent mixing which can be applied to more complex geometries and flows. This paper focuses on a number of experimental techniques which have been developed at the National Institute of Standards and Development for measuring concentration in binary axisymmetric turbulent jets. In order to demonstrate the value of these diagnostics, some of the more important results from earlier and on-going investigations are summarized. Topics addressed include the similarity behavior of variable density axisymmetric jets, the behavior of absolutely unstable axisymmetric helium jets, and the role of large scale structures and scalar dissipation in these flows.

Pitts, William M.↗

NASA NTRS1997

Bioreactor Mass Transport Studies

The objectives of the proposed research efforts were to develop both a simulation tool and a series of experiments to provide a quantitative assessment of mass transport in the NASA rotating wall perfused vessel (RWPV) bioreactor to be flown on EDU#2. This effort consisted of a literature review of bioreactor mass transport studies, the extension of an existing scalar transport computer simulation to include production and utilization of the scalar, and the evaluation of experimental techniques for determining mass transport in these vessels. Since mass transport at the cell surface is determined primarily by the relative motion of the cell assemblage and the surrounding fluid, a detailed assessment of the relative motion was conducted. Results of the simulations of the motion of spheres in the RWPV under microgravity conditions are compared with flight data from EDU#1 flown on STS-70. The mass transport across the cell membrane depends upon the environment, the cell type, and the biological state of the cell. Results from a literature review of cell requirements of several scalars are presented. As a first approximation, a model with a uniform spatial distribution of utilization or production was developed and results from these simulations are presented. There were two candidate processes considered for the experimental mass transport evaluations. The first was to measure the dissolution rate of solid or gel beads. The second was to measure the induced fluorescence of beads as a stimulant (for example hydrogen peroxide) is infused into the vessel. Either technique would use video taped images of the process for recording the quantitative results. Results of preliminary tests of these techniques are discussed.

Kleis, Stanley J.↗

NASA NTRS1997

Mixing Enhancement in a Lobed Injector

An experimental investigation of the non-reactive mixing processes associated with a lobed fuel injector in a coflowing air stream is presented. The lobed fuel injector is a device which generates streamwise vorticity, producing high strain rates which can enhance the mixing of reactants while delaying ignition in a controlled manner. The lobed injectors examined in the present study consist of two corrugated plates between which a fuel surrogate, CO2, is injected into coflowing air. Acetone is seeded in the CO2 supply as a fuel marker. Comparison of two alternative lobed injector geometries is made with a straight fuel injector to determine net differences in mixing and strain fields due to streamwise vorticity generation. Planar laser-induced fluorescence (PLIF) of the seeded acetone yields two-dimensional images of the scalar concentration field at various downstream locations, from which local mixing and scalar dissipation rates are computed. It is found that the lobed injector geometry can enhance molecular mixing and create a highly strained flowfield, and that the strain rates generated by scalar energy dissipation can potentially delay ignition in a reacting flowfield.

Smith, L. L.↗

NASA NTRS1994

A Navier-Strokes Chimera Code on the Connection Machine CM-5: Design and Performance

We have implemented a three-dimensional compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the 'chimera' approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. A parallel machine like the CM-5 is well-suited for finite-difference methods on structured grids. The regular pattern of connections of a structured mesh maps well onto the architecture of the machine. So the first design choice, finite differences on a structured mesh, is natural. We use centered differences in space, with added artificial dissipation terms. When numerically solving the Navier-Stokes equations, there are liable to be some mesh cells near a solid body that are small in at least one direction. This mesh cell geometry can impose a very severe CFL (Courant-Friedrichs-Lewy) condition on the time step for explicit time-stepping methods. Thus, though explicit time-stepping is well-suited to the architecture of the machine, we have adopted implicit time-stepping. We have further taken the approximate factorization approach. This creates the need to solve large banded linear systems and creates the first possible barrier to an efficient algorithm. To overcome this first possible barrier we have considered two options. The first is just to solve the banded linear systems with data spread over the whole machine, using whatever fast method is available. This option is adequate for solving scalar tridiagonal systems, but for scalar pentadiagonal or block tridiagonal systems it is somewhat slower than desired. The second option is to 'transpose' the flow and geometry variables as part of the time-stepping process: Start with x-lines of data in-processor. Form explicit terms in x, then transpose so y-lines of data are in-processor. Form explicit terms in y, then transpose so z-lines are in processor. Form explicit terms in z, then solve linear systems in the z-direction. Transpose to the y-direction, then solve linear systems in the y-direction. Finally transpose to the x direction and solve linear systems in the x-direction. This strategy avoids inter-processor communication when differencing and solving linear systems, but requires a large amount of communication when doing the transposes. The transpose method is more efficient than the non-transpose strategy when dealing with scalar pentadiagonal or block tridiagonal systems. For handling geometrically complex problems the chimera strategy was adopted. For multiple zone cases we compute on each zone sequentially (using the whole parallel machine), then send the chimera interpolation data to a distributed data structure (array) laid out over the whole machine. This information transfer implies an irregular communication pattern, and is the second possible barrier to an efficient algorithm. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran. We make use of the Connection Machine Scientific Software Library (CMSSL) for the linear solver and array transpose operations.

Jespersen, Dennis C.↗

NASA NTRS2003

Subgrid Combustion Modeling for the Next Generation National Combustion Code

In the first year of this research, a subgrid turbulent mixing and combustion methodology developed earlier at Georgia Tech has been provided to researchers at NASA/GRC for incorporation into the next generation National Combustion Code (called NCCLES hereafter). A key feature of this approach is that scalar mixing and combustion processes are simulated within the LES grid using a stochastic 1D model. The subgrid simulation approach recovers locally molecular diffusion and reaction kinetics exactly without requiring closure and thus, provides an attractive feature to simulate complex, highly turbulent reacting flows of interest. Data acquisition algorithms and statistical analysis strategies and routines to analyze NCCLES results have also been provided to NASA/GRC. The overall goal of this research is to systematically develop and implement LES capability into the current NCC. For this purpose, issues regarding initialization and running LES are also addressed in the collaborative effort. In parallel to this technology transfer effort (that is continuously on going), research has also been underway at Georgia Tech to enhance the LES capability to tackle more complex flows. In particular, subgrid scalar mixing and combustion method has been evaluated in three distinctly different flow field in order to demonstrate its generality: (a) Flame-Turbulence Interactions using premixed combustion, (b) Spatially evolving supersonic mixing layers, and (c) Temporal single and two-phase mixing layers. The configurations chosen are such that they can be implemented in NCCLES and used to evaluate the ability of the new code. Future development and validation will be in spray combustion in gas turbine engine and supersonic scalar mixing.

Menon, Suresh↗

NASA NTRS2003

Stochastic Mixing Model with Power Law Decay of Variance

Here we present a simple stochastic mixing model based on the law of large numbers (LLN). The reason why the LLN is involved in our formulation of the mixing problem is that the random conserved scalar c = c(t,x(t)) appears to behave as a sample mean. It converges to the mean value mu, while the variance sigma(sup 2)(sub c) (t) decays approximately as t(exp -1). Since the variance of the scalar decays faster than a sample mean (typically is greater than unity), we will introduce some non-linear modifications into the corresponding pdf-equation. The main idea is to develop a robust model which is independent from restrictive assumptions about the shape of the pdf. The remainder of this paper is organized as follows. In Section 2 we derive the integral equation from a stochastic difference equation describing the evolution of the pdf of a passive scalar in time. The stochastic difference equation introduces an exchange rate gamma(sub n) which we model in a first step as a deterministic function. In a second step, we generalize gamma(sub n) as a stochastic variable taking fluctuations in the inhomogeneous environment into account. In Section 3 we solve the non-linear integral equation numerically and analyze the influence of the different parameters on the decay rate. The paper finishes with a conclusion.

Fedotov, S.↗

NASA NTRS2003

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Raman, Venkatramanan↗

NASA NTRS1997

Computation of LFT Uncertainty Bounds with Repeated Parametric Uncertainties

A new methodology in which linear fractional transformation uncertainty bounds are directly constructed for use in robust control design and analysis is proposed. Existence conditions for model validating solutions with or without repeated scalar uncertainty are given. The approach is based on minimax formulation to deal with multiple non-repeated structured uncertainty components subject to fixed levels of repeated scalar uncertainties. Input directional dependence and variations with different experiments are addressed by maximizing uncertainty levels over multiple experimental data sets. Preliminary results show that reasonable uncertainty bounds on structured non-repeated uncertainties can be identified directly from measurement data by assuming reasonable levels of repeated scalar uncertainties.

Lim, K. B.↗

NASA NTRS2005

A Spreadsheet for the Mixing of a Row of Jets with a Confined Crossflow

An interactive computer code, written with a readily available software program, Microsoft Excel (Microsoft Corporation, Redmond, WA) is presented which displays 3 D oblique plots of a conserved scalar distribution downstream of jets mixing with a confined crossflow, for a single row, double rows, or opposed rows of jets with or without flow area convergence and/or a non-uniform crossflow scalar distribution. This project used a previously developed empirical model of jets mixing in a confined crossflow to create an Microsoft Excel spreadsheet that can output the profiles of a conserved scalar for jets injected into a confined crossflow given several input variables. The program uses multiple spreadsheets in a single Microsoft Excel notebook to carry out the modeling. The first sheet contains the main program, controls for the type of problem to be solved, and convergence criteria. The first sheet also provides for input of the specific geometry and flow conditions. The second sheet presents the results calculated with this routine to show the effects on the mixing of varying flow and geometric parameters. Comparisons are also made between results from the version of the empirical correlations implemented in the spreadsheet and the versions originally written in Applesoft BASIC (Apple Computer, Cupertino, CA) in the 1980's.

J D Holderman↗

NASA NTRS2008

Effects of C/O Ratio and Temperature on Sooting Limits of Spherical Diffusion Flames

Limiting conditions for soot particle inception in spherical diffusion flames were investigated numerically. The flames were modeled using a one-dimensional, time accurate diffusion flame code with detailed chemistry and transport and an optically thick radiation model. Seventeen normal and inverse flames were considered, covering a wide range of stoichiometric mixture fraction, adiabatic flame temperature, residence time and scalar dissipation rate. These flames were previously observed to reach their sooting limits after 2 s of microgravity. Sooting-limit diffusion flames with scalar dissipation rate lower than 2/s were found to have temperatures near 1400 K where C/O = 0.51, whereas flames with greater scalar dissipation rate required increased temperatures. This finding was valid across a broad range of fuel and oxidizer compositions and convection directions.

Lecoustre, V. R.↗

NASA NTRS2007

Space Propulsion Design and Analysis

This software provides an improved methodology for predicting launcher base pressure and heat loads for RSRM (Reusable Solid Rocket Motor) launchers by accounting for complex anisotropic stress/strains and variable turbulent Prandtl and Schmidt numbers. A "building block" approach to turbulence model development, and validation has been applied for improved missile/launcher base region analysis. Modifications to existing kappa - epsilon turbulence models and application of scalar variance models are incorporated into a RANS-based method for aeropropulsive flow modeling, directly related to base flow methodology. (RANS stands for Reynolds-averaged Navier-Stokes.) The models are applied in a RANS solver framework and can improve analysis of other complex flow fields. The enhanced models provide a more accurate predictive capability for improving the design and analysis of RSRM launcher configuration. The kappa - epsilon model enhancements have been shown to improve the capability for predicting turbulence effects in base blow environments. The scalar variance models have been assessed over a wide range of flow configurations to improve prediction of turbulent scalar mixing.

Sinha, Neeraj↗