Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Three-dimensional plasma diffusion in a very strong magnetic field.

The thermal equilibrium coefficient of spatial diffusion transverse to a strong uniform dc magnetic field is calculated for a fully ionized plasma. The particles are assumed to move transverse to the field only as a consequence of the E x B drift, but to move freely parallel to it. The particles interact only electrostatically. The calculation is done at large, but fixed and finite, plasma volume. It is shown that, as magnetic field density approaches infinity, the leading term of the coefficient of transverse spatial diffusion falls off in a prescribed manner, but contains a multiplicative factor which goes to zero as the plasma volume becomes infinite. The method of calculation fails for unbounded plasmas.

Montgomery, D.↗
NASA NTRS1976

Characterizations of linear sufficient statistics

A necessary and sufficient condition is developed such that there exists a continous linear sufficient statistic T for a dominated collection of totally finite measures defined on the Borel field generated by the open sets of a Banach space X. In particular, corollary necessary and sufficient conditions are given so that there exists a rank K linear sufficient statistic T for any finite collection of probability measures having n-variate normal densities. In this case a simple calculation, involving only the population means and covariances, determines the smallest integer K for which there exists a rank K linear sufficient statistic T (as well as an associated statistic T itself).

Peters, B. C., Jr.↗
NASA NTRS1977

Numerical solution of the compressible Navier-Stokes equations using density gradients as additional dependent variables

Numerical solution of two dimensional, time dependent, compressible viscous Navier-Stokes equations about arbitrary bodies was treated using density gradients as additional dependent variables. Thus, six dependent variables were computed with the SOR iteration method. Besides formulation for pressure gradient terms, a formulation for computing the body density was presented. To approximate the governing equations, an implicit finite difference method was employed. In computing the solution for the flow about a circular cylinder, a problem arose near the wall at both stagnation points. Thus, computations with various conditions were tried to examine the problem. Also, computations with and without formulations are compared. The flow variables were computed on 37 by 40 field first, then on an 81 by 40 field.

Kwon, J. H.↗
NASA NTRS1986

Soliton turbulence

Theoretical and numerical works in atmospheric turbulence have used the Navier-Stokes fluid equations exclusively for describing large-scale motions. Controversy over the existence of an average temperature gradient for the very large eddies in the atmosphere suggested that a new theoretical basis for describing large-scale turbulence was necessary. A new soliton formalism as a fluid analogue that generalizes the Schrodinger equation and the Zakharov equations has been developed. This formalism, processing all the nonlinearities including those from modulation provided by the density fluctuations and from convection due to the emission of finite sound waves by velocity fluctuations, treats large-scale turbulence as coalescing and colliding solitons. The new soliton system describes large-scale instabilities more explicitly than the Navier-Stokes system because it has a nonlinearity of the gradient type, while the Navier-Stokes has a nonlinearity of the non-gradient type. The forced Schrodinger equation for strong fluctuations describes the micro-hydrodynamical state of soliton turbulence and is valid for large-scale turbulence in fluids and plasmas where internal waves can interact with velocity fluctuations.

Tchen, C. M.↗
NASA NTRS1989

A numerical study of scalar gradients in Kelvin-Helmholtz billows

A high resolution numerical technique is used to model the development of a periodically perturbed shear layer imbedded in an initially vertical gradient of a passive scalar. The technique follows the development of the vorticity through an initial linear growth state and well into the nonlinear development of Kelvin-Helmholtz billows, in the zero-viscosity, zero-diffusion limit. The resulting scalar distribution rapidly develops regions of extremely sharp scalar gradients, which wind around the periodically spaced vortical low gradient cores. Vertical cross sections through different parts of the billow structure are presented and compared with rocket measurements of electron density fine structure in the mesosphere. Gradient limits imposed by finite diffusion are calculated, and implications for atmospheric radar observations are discussed.

Parker, J. W.↗
NASA NTRS1989

A Simple Algebraic Grid Adaptation Scheme with Applications to Two- and Three-dimensional Flow Problems

An algebraic adaptive grid scheme based on the concept of arc equidistribution is presented. The scheme locally adjusts the grid density based on gradients of selected flow variables from either finite difference or finite volume calculations. A user-prescribed grid stretching can be specified such that control of the grid spacing can be maintained in areas of known flowfield behavior. For example, the grid can be clustered near a wall for boundary layer resolution and made coarse near the outer boundary of an external flow. A grid smoothing technique is incorporated into the adaptive grid routine, which is found to be more robust and efficient than the weight function filtering technique employed by other researchers. Since the present algebraic scheme requires no iteration or solution of differential equations, the computer time needed for grid adaptation is trivial, making the scheme useful for three-dimensional flow problems. Applications to two- and three-dimensional flow problems show that a considerable improvement in flowfield resolution can be achieved by using the proposed adaptive grid scheme. Although the scheme was developed with steady flow in mind, it is a good candidate for unsteady flow computations because of its efficiency.

Hsu, Andrew T.↗
NASA NTRS1990

Simulation of on-orbit modal tests of large space structures

This paper describes a procedure to analytically simulate a modal test of an on-orbit large space structure (LSS), extract the modal properties, and evaluate the success of the modal test. This procedure addresses some of the major challenges to performing an on-orbit modal test of an LSS including high modal density, low frequency modes, and limitations in excitation capabilities. A finite element model of the orbiting structure is used to predict acceleration responses due to thruster excitations, time-domain modal extraction methods are used to estimate the modal properties, and comparison of frequencies and cross-orthogonality values is used to evaluate the success of the modal test. Several alternative excitation patterns and sensor arrangements were evaluated using a space station model as an example. Results of the simulations indicate that the choice of excitation functions is critical to the success of the test.

Blelloch, Paul↗
NASA NTRS1993

Quark matter droplets in neutron stars

We show that, for physically reasonable bulk and surface properties, the lowest energy state of dense matter consists of quark matter coexisting with nuclear matter in the presence of an essentially uniform background of electrons. We estimate the size and nature of spatial structure in this phase, and show that at the lowest densities the quark matter forms droplets embedded in nuclear matter, whereas at higher densities it can exhibit a variety of different topologies. A finite fraction of the interior of neutron stars could consist of matter in this new phase, which would provide new mechanisms for glitches and cooling.

Heiselberg, H.↗
NASA NTRS1996

Analysis of Flowfields over Four-Engine DC-X Rockets

The objective of this study is to validate a computational methodology for the aerodynamic performance of an advanced conical launch vehicle configuration. The computational methodology is based on a three-dimensional, viscous flow, pressure-based computational fluid dynamics formulation. Both wind-tunnel and ascent flight-test data are used for validation. Emphasis is placed on multiple-engine power-on effects. Computational characterization of the base drag in the critical subsonic regime is the focus of the validation effort; until recently, almost no multiple-engine data existed for a conical launch vehicle configuration. Parametric studies using high-order difference schemes are performed for the cold-flow tests, whereas grid studies are conducted for the flight tests. The computed vehicle axial force coefficients, forebody, aftbody, and base surface pressures compare favorably with those of tests. The results demonstrate that with adequate grid density and proper distribution, a high-order difference scheme, finite rate afterburning kinetics to model the plume chemistry, and a suitable turbulence model to describe separated flows, plume/air mixing, and boundary layers, computational fluid dynamics is a tool that can be used to predict the low-speed aerodynamic performance for rocket design and operations.

Wang, Ten-See↗
NASA NTRS1994

On the wavelet optimized finite difference method

When one considers the effect in the physical space, Daubechies-based wavelet methods are equivalent to finite difference methods with grid refinement in regions of the domain where small scale structure exists. Adding a wavelet basis function at a given scale and location where one has a correspondingly large wavelet coefficient is, essentially, equivalent to adding a grid point, or two, at the same location and at a grid density which corresponds to the wavelet scale. This paper introduces a wavelet optimized finite difference method which is equivalent to a wavelet method in its multiresolution approach but which does not suffer from difficulties with nonlinear terms and boundary conditions, since all calculations are done in the physical space. With this method one can obtain an arbitrarily good approximation to a conservative difference method for solving nonlinear conservation laws.

Jameson, Leland↗
NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗
NASA NTRS1981

Full potential solution of transonic quasi-3-D flow through a cascade using artificial compressability

The three-dimensional flow in a turbomachinery blade row was approximated by correcting for streamtube convergence and radius change in the throughflow direction. The method is a fully conservative solution of the full potential equation incorporating the finite volume technique on body fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. Comparison of results for several supercritical blades shows good agreement with their hodograph solutions. Other calculations for these profiles as well as standard NACA blade sections indicate that this is a useful scheme analyzing both the design and off-design performance of turbomachinery blading.

Farrell, C.↗
NASA NTRS1993

A C-band backscatter model for lake ice in Alaska

ERS-1 SAR imagery of lake ice growing on shallow tundra lakes in northern Alaska shows interesting radar backscatter variations. Based on the analysis of ice cores from these lakes, a multi-layer backscatter model comprised of the following elements has been developed: (1) specular air-ice; ice-water and ice-frozen soil boundaries; (2) an ice layer of variable thickness; (3) ice sub-layers with air inclusions of variable density, size and shape including spheres, prolate spheroids, and cylinders of finite length. Preliminary model results confirm that backscatter is a sensitive function of greater reflectivity than from an ice-frozen soil interface. The model has also been tested using bubble data derived from ice cores in April 1992. The modelled backscatter is compared with backscatter derived from ERS-1 SAR images obtained at the same time as the fieldwork.

Wakabayashi, H.↗
NASA NTRS1999

Charge and Spin Dynamics of the Hubbard Chains

We calculate the local correlation functions of charge and spin for the one-chain and two-chain Hubbard model using density matrix renormalization group method and the recursion technique. Keeping only finite number of states we get good accuracy for the low energy excitations. We study the charge and spin gaps, bandwidths and weights of the spectra for various values of the on-site Coulomb interaction U and the electron filling. In the low energy part, the local correlation functions are different for the charge and spin. The bandwidths are proportional to t for the charge and J for the spin respectively.

Park, Youngho↗
NASA NTRS1985

Automatic crack propagation tracking

A finite element based approach to fully automatic crack propagation tracking is presented. The procedure presented combines fully automatic mesh generation with linear fracture mechanics techniques in a geometrically based finite element code capable of automatically tracking cracks in two-dimensional domains. The automatic mesh generator employs the modified-quadtree technique. Crack propagation increment and direction are predicted using a modified maximum dilatational strain energy density criterion employing the numerical results obtained by meshes of quadratic displacement and singular crack tip finite elements. Example problems are included to demonstrate the procedure.

Shephard, M. S.↗
NASA NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗