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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗
NASA NTRS1994

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗
NASA NTRS1973

Apollo 15 measurement of lunar surface brightness temperatures - Thermal conductivity of the upper 1.5 meters of regolith.

In situ measurements of lunar surface brightness temperatures made as a part of the Apollo Lunar Surface Experiments Package at the Apollo 15 Hadley Rille landing site are reported. Data derived from five thermocouples of the Heat Flow Experiment, which are lying on or just above the surface, are used to examine the thermal properties of the upper 15 cm of the lunar regolith using eclipse and nighttime cool-down temperatures. Application of finite-difference techniques in modeling the lunar soil shows that the thermocouple data are best fit by a model consisting of a low-density and low-thermal conductivity surface layer approximately 2 cm thick overlying a region increasing in conductivity and density with depth. Conductivities on the order of 0.00001 W per cm per deg K are postulated for the upper layer, with conductivity increasing to the order of 0.0001 W per cm per deg K at depths exceeding 20 cm. An increase in mean temperature with depth indicates that the ratio of radiative to conductive transfer at 350 K is 2.7 for at least the upper few centimeters of lunar soil; this value is nearly twice that measured for returned lunar fines.

Keihm, S. J.↗
NASA NTRS1985

Modeling of dense sprays from LOX/H2 coaxial injectors under supercritical conditions

Assumptions commonly employed in the analytical description of a thin spray become unacceptable as the droplet density and droplet volume displacement become large. Changes in the semiempirical models describing drop-drop or drop-gas interactions as well as in the basic modeling approach itself may be necessary. A number of efforts undertaken to improve such models of the atomization-evaporation process to render them applicable in the dense-spray, supercritical regime are described. The improvements include those for variable drop density, supercritical properties, stripping evaporation and drop reincorporation into the liquid stream. The models are developed for use with the ARICC three-phase finite difference combustion code. Selected results from a sample calculation are presented, demonstrating the enhanced realism achievable with the upgraded models.

Liang, P. Y.↗
NASA NTRS1988

Viscous shock-layer solutions for the low-density hypersonic flow past long slender bodies

Results are obtained for the surface pressure, drag, heat-transfer, and skin-friction coefficients for hyperboloids and sphere cones. Body half angles from 5 to 22.5 degrees are considered for various low-density flow conditions. Recently obtained surface-slip and shock-slip equations are employed to account for the low-density effects. The method of solution employed for the viscous shock-layer (VSL) equations is a partially coupled spatial-marching implicit finite-difference technique. The flow cases analyzed include highly cooled long slender bodies in high Mach number flows. The present perfect-gas VSL calculations compare quite well with available experimental data. Results have also been obtained from the steady-state Navier-Stokes (NS) equations by successive approximations. Comparison between the NS and VSL results indicates that VSL equations even with body and shock-slip boundary conditions may not be adequate in the stagnation region at altitudes greater than about 75 km for the cases analyzed here.

Gupta, R. N.↗
NASA NTRS1994

Experimental verification of a progressive damage model for composite laminates based on continuum damage mechanics

Progressive failure is a crucial concern when using laminated composites in structural design. Therefore the ability to model damage and predict the life of laminated composites is vital. The purpose of this research was to experimentally verify the application of the continuum damage model, a progressive failure theory utilizing continuum damage mechanics, to a toughened material system. Damage due to tension-tension fatigue was documented for the IM7/5260 composite laminates. Crack density and delamination surface area were used to calculate matrix cracking and delamination internal state variables, respectively, to predict stiffness loss. A damage dependent finite element code qualitatively predicted trends in transverse matrix cracking, axial splits and local stress-strain distributions for notched quasi-isotropic laminates. The predictions were similar to the experimental data and it was concluded that the continuum damage model provided a good prediction of stiffness loss while qualitatively predicting damage growth in notched laminates.

Coats, Timothy William↗
NASA NTRSDate not supplied

Preparing Quantum Many-Body Scar States on Quantum Computers

Quantum many-body scar states are highly excited eigenstates of many-body systems that exhibit atypical entanglement and correlation properties relative to other eigenstates at the same energy density. Scar states also give rise to infinitely long-lived coherent dynamics when the system is prepared in a special initial state having finite overlap with them. Many models with exact scar states have been constructed, but the fate of scarred eigenstates and dynamics when these models are perturbed is difficult to study with classical computational techniques. In this work, we propose state preparation protocols for individual scar states in a particular model, as well as superpositions of them that give rise to coherent dynamics. For superpositions of scar states, we propose both a linear depth unitary and a finite-depth nonunitary state preparation protocol, the latter of which uses measurement and postselection to reduce the circuit depth. For individual scarred eigenstates, we propose a circuit architecture with polynomial depth. We also provide proof of principle implementations of these protocols on superconducting quantum hardware.

Quantum Computing↗
NASA NTRS2002

Coupling Between Microstrip Lines and Finite Ground Coplanar Lines Embedded in Polyimide Layers for 3D-MMICs on Silicon

Three-dimensional circuits built upon multiple layers of polyimide are required for constructing Si/SiGe monolithic microwave/mm-wave integrated circuits on CMOS (low resistivity) Si wafers. It is expected that these circuits will replace the ones fabricated on GaAs and reduce the overall system cost. However, the closely spaced transmission lines that are required for a high-density circuit environment are susceptible to high levels of cross-coupling, which degrades the overall circuit performance. In this paper, theoretical and experimental results on coupling and ways to reduce it are presented for two types of transmission lines: a) the microstrip line and b) the Finite Ground Coplanar (FGC) line. For microstrip lines it is shown that a fence of metalized via-holes can significantly reduce coupling, especially in the case when both lines are on the same polyimide layer or when the shielding structure extends through several polyimide layers. For closely spaced microstrip lines, coupling is lower for a metal filled trench shield than a via-hole fence. Coupling amongst microstrip lines is dependent on the ratio of line separation to polyimide thickness and is primarily due to magnetic fields. For FGC lines it is shown that they have in general low coupling that can be reduced significantly when there is even a small gap between the ground planes of each line. FGC lines have approximately 8 dB lower coupling than coupled coplanar waveguides (CPW). In addition, forward and backward characteristics of the FGC lines do not resemble those of other transmission lines such as microstrip. Therefore, the coupling mechanism of the FGC lines is different compared to thin film microstrip lines.

Ponchak, G. E.↗
NASA NTRS2003

Structural Reliability Using Probability Density Estimation Methods Within NESSUS

A reliability analysis studies a mathematical model of a physical system taking into account uncertainties of design variables and common results are estimations of a response density, which also implies estimations of its parameters. Some common density parameters include the mean value, the standard deviation, and specific percentile(s) of the response, which are measures of central tendency, variation, and probability regions, respectively. Reliability analyses are important since the results can lead to different designs by calculating the probability of observing safe responses in each of the proposed designs. All of this is done at the expense of added computational time as compared to a single deterministic analysis which will result in one value of the response out of many that make up the density of the response. Sampling methods, such as monte carlo (MC) and latin hypercube sampling (LHS), can be used to perform reliability analyses and can compute nonlinear response density parameters even if the response is dependent on many random variables. Hence, both methods are very robust; however, they are computationally expensive to use in the estimation of the response density parameters. Both methods are 2 of 13 stochastic methods that are contained within the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) program. NESSUS is a probabilistic finite element analysis (FEA) program that was developed through funding from NASA Glenn Research Center (GRC). It has the additional capability of being linked to other analysis programs; therefore, probabilistic fluid dynamics, fracture mechanics, and heat transfer are only a few of what is possible with this software. The LHS method is the newest addition to the stochastic methods within NESSUS. Part of this work was to enhance NESSUS with the LHS method. The new LHS module is complete, has been successfully integrated with NESSUS, and been used to study four different test cases that have been proposed by the Society of Automotive Engineers (SAE). The test cases compare different probabilistic methods within NESSUS because it is important that a user can have confidence that estimates of stochastic parameters of a response will be within an acceptable error limit. For each response, the mean, standard deviation, and 0.99 percentile, are repeatedly estimated which allows confidence statements to be made for each parameter estimated, and for each method. Thus, the ability of several stochastic methods to efficiently and accurately estimate density parameters is compared using four valid test cases. While all of the reliability methods used performed quite well, for the new LHS module within NESSUS it was found that it had a lower estimation error than MC when they were used to estimate the mean, standard deviation, and 0.99 percentile of the four different stochastic responses. Also, LHS required a smaller amount of calculations to obtain low error answers with a high amount of confidence than MC. It can therefore be stated that NESSUS is an important reliability tool that has a variety of sound probabilistic methods a user can employ and the newest LHS module is a valuable new enhancement of the program.

Chamis, Chrisos C.↗
NASA NTRS1982

Propagation of wide bandwidth signals through strongly turbulent ionized media

Analytic and numerical techniques are presented which directly address the problem of propagation of wide bandwidth signals through random ionized media. This work is applicable to the problems of satellite communication and space based radar observation through a disturbed ionospheric propagation channel that would result from a high altitude chemical release or nuclear detonation. An analytic solution is presented for the two-position, two-frequency mutual coherence function for spherical wave propagation in the strong scattering limit. This solution is used to derive simple expressions for the power impulse response function and to determine its relationship to the important parameters of decorrelation distance and coherence bandwidth which describe the disturbed propagation channel. Results for mean time delay and time delay jitter are presented and compared to direct simulation results and to other theoretical calculations. A numerical/analytical solution to the parabolic wave equation is presented in the form of a multiple phase-screen (MPS) propagation simulation. In this very general solution technique, the random medium is divided into a finite number of layers. The field fluctuation through each layer is obtained by replacing the layer by a centrally located thin phase-screen, whose statistical properties are determined from the statistics of the electron-density irregularities. The waveform then propagates from phase-screen to phase-screen via an exact solution to the Fresnel-Kirchnoff equation. For wide bandwidth waveforms, numerical solutions are obtained at a number of discrete frequencies centered about the carrier and then time-domain solutions are obtained by Fourier transform techniques. Detailed results are presented for a number of interesting cases including propagation of a 20 MHz bandwidth signal through a finite barium cloud at a carrier frequency of 100 MHz. One of the major uses of the MPS propagation simulation has been to provide realizations of the received signal after propagation through a disturbed channel. The MPS simulation obtains a general solution of the parabolic wave equation under both weak and strong scattering conditions. A second technique to directly obtain realizations of wide bandwidth waveforms is presented here. This technique is limited to the case of strong scattering but requires only a fraction of the computer resources needed for MPS signal generation. Detailed comparisons of the two signal generation techniques are presented.

Random Media↗
NASA NTRS2023

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗
NASA NTRS1974

Models of Saturn's rings which satisfy the optical observations

A theoretical model of Saturn's rings is investigated which includes the shadowing effect and realistic anisotropic phase functions for the ring particles. The effects of multiple scattering and the finite size of the sun, including the penumbra, are rigorously included. The permissible range of the relevant parameters, including optical thickness, single scattering albedo, volume density, and phase function are investigated by comparing the theoretical results to observations of the ring brightness vs phase angle, wavelength, and elevation of the sun and earth. Anisotropic scattering by the ring particles is necessary in order to match the observations. The color dependence of the opposition effect is interpreted in terms of the albedo spectrum of the ring particles.

Kawata, Y.↗
NASA NTRS1986

Unsteady design-point flow phenomena in transonic compressors

High-frequency response probes which had previously been used exclusively in the MIT Blowndown Facility were successfully employed in two conventional steady state axial flow compressor facilities to investigate the unsteady flowfields of highly loaded transonic compressors at design point operation. Laser anemometry measurements taken simultaneously with the high response data were also analyzed. The time averaged high response data of static and total pressure agreed quite well with the conventional steady state instrumentation except for flow angle which showed a large spread in values at all radii regardless of the type of instrumentation used. In addition, the time resolved measurements confirmed earlier test results obtained in the MIT Blowdown Facility for the same compressor. The results of these tests have further revealed that the flowfields of highly loaded transonic compressors are heavily influenced by unsteady flow phenomena. The high response measurements exhibited large variations in the blade to blade flow and in the blade passage flow. The observed unsteadiness in the blade wakes is explained in terms of the rotor blades' shed vorticity in periodic vortex streets. The wakes were modeled as two-dimensional vortex streets with finite size cores. The model fit the data quite well as it was able to reproduce the average wake shape and bi-modal probability density distributions seen in the laser anemometry data. The presence of vortex streets in the blade wakes also explains the large blade to blade fluctuations seen by the high response probes which is simply due to the intermittent sampling of the vortex street as it is swept past a stationary probe.

Gertz, J. B.↗
NASA NTRS1988

Deformational models of rifting and folding on Venus

Features of presumed tectonic origin on Venus are reviewed, and lithospheric strength envelopes are derived based on laboratory measurements of the deformational properties of crustal and subcrustal rocks, extrapolated to conditions appropriate to Venus. Models for rifting and folding are developed that use this lithospheric structure and take into account both brittle and ductile yielding as well as finite elastic strength. For both rifting and folding, structures with characteristic widths and spacings are predicted whose size depends on the thickness of the lithosphere, density contrast, and elastic properties of the layer. Finally, the model predictions are compared with the widths and spacings of observed tectonic features, and it is concluded that they are consistent with a relatively strong mantle layer separated from a thin brittle surface layer by a ductile lower crust. These results allow constraints to be placed on the crustal thickness and thermal gradient on Venus.

Banerdt, W. Bruce↗
NASA NTRS1990

Comparative dynamics of the ionospheres of Venus and Mars at large solar zenith angles

Simple three-ion spectral models of the dynamics and ion chemistry in the upper ionospheres of Venus and Mars have been constructed by including two minor ions, H(+) and O2(+) for Venus and O(+) and CO2(+) for Mars, along with the major ion, O(+) for Venus and O2(+) for Mars. Horizontal flow velocities and ion densities have been calculated in each case and compared with available measured data. For Venus, the present calculations yield results in agreement with earlier finite-difference model studies and with experimental observations. Calculations for Mars are in agreement with the limited observational data obtained by the Viking landers.

Singhal, R. P.↗
NASA NTRS2008

Comparison of ALE and SPH Simulations of Vertical Drop Tests of a Composite Fuselage Section into Water

Simulation of multi-terrain impact has been identified as an important research area for improved prediction of rotorcraft crashworthiness within the NASA Subsonic Rotary Wing Aeronautics Program on Rotorcraft Crashworthiness. As part of this effort, two vertical drop tests were conducted of a 5-ft-diameter composite fuselage section into water. For the first test, the fuselage section was impacted in a baseline configuration without energy absorbers. For the second test, the fuselage section was retrofitted with a composite honeycomb energy absorber. Both tests were conducted at a nominal velocity of 25-ft/s. A detailed finite element model was developed to represent each test article and water impact was simulated using both Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) approaches in LS-DYNA, a nonlinear, explicit transient dynamic finite element code. Analytical predictions were correlated with experimental data for both test configurations. In addition, studies were performed to evaluate the influence of mesh density on test-analysis correlation.

Jackson, Karen E.↗
NASA NTRS1992

FDDO and DSMC analyses of rarefied gas flow through 2D nozzles

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas expanding through a two-dimensional nozzle and into a surrounding low-density environment. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, the variable hard sphere model is used as a molecular model and the no time counter method is employed as a collision sampling technique. The results of both the FDDO and the DSMC methods show good agreement. The FDDO method requires less computational effort than the DSMC method by factors of 10 to 40 in CPU time, depending on the degree of rarefaction.

Chung, Chan-Hong↗
NASA NTRS1997

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗