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 109 records · Page 6

Interstellar scattering of the radio source 2013+370

Results are reported from a six-station 1663-MHz VLBI observation of the radio source 2013+370. The image of this source is dominated by interstellar scattering and is observed to be elliptical, with a major axis (FWHM) of 17 marcsec and an axial ratio of of 0.70. A fit of the observations to a model visibility function yields a value for alpha, the spectral index of the power-law density irregularity spectrum, of alpha=3.79 + or - 0.05. Corrections for the finite angular size of the source would probably produce a downward revision of alpha. The present observations are consistent with the interstellar density spectrum being close to Kolmogoroff for those scales which determine the angular broadening.

Spangler, Steven R.↗

Charging in the ac Conductance of a Double Barrier Resonant Tunneling Structure

There have been many studies of the linear response ac conductance of a double barrier resonant tunneling structure (DBRTS), both at zero and finite dc biases. While these studies are important, they fail to self consistently include the effect of the time dependent charge density in the well. In this paper, we calculate the ac conductance at both zero and finite do biases by including the effect of the time dependent charge density in the well in a self consistent manner. The charge density in the well contributes to both the flow of displacement currents in the contacts and the time dependent potential in the well. We find that including these effects can make a significant difference to the ac conductance and the total ac current is not equal to the simple average of the non-selfconsistently calculated conduction currents in the two contacts. This is illustrated by comparing the results obtained with and without the effect of the time dependent charge density included correctly. Some possible experimental scenarios to observe these effects are suggested.

Anantram, M. P.↗

Geometrically Nonlinear Finite Element Analysis of a Composite Space Reflector

Lightweight aerospace structures, such as low areal density composite space reflectors, are highly flexible and may undergo large deflection under applied loading, especially during the launch phase. Accordingly, geometrically nonlinear analysis that takes into account the effect of finite rotation may be needed to determine the deformed shape for a clearance check and the stress and strain state to ensure structural integrity. In this study, deformation of the space reflector is determined under static conditions using a geometrically nonlinear solid shell finite element model. For the solid shell element formulation, the kinematics of deformation is described by six variables that are purely vector components. Because rotational angles are not used, this approach is free of the limitations of small angle increments. This also allows easy connections between substructures and large load increments with respect to the conventional shell formulation using rotational parameters. Geometrically nonlinear analyses were carried out for three cases of static point loads applied at selected points. A chart shows results for a case when the load is applied at the center point of the reflector dish. The computed results capture the nonlinear behavior of the composite reflector as the applied load increases. Also, they are in good agreement with the data obtained by experiments.

Lee, Kee-Joo↗

A phenomenological pulsar model

Particle injection energies and rates previously calculated for the stellar wind generation by rotating magnetized neutron stars are adopted. It is assumed that the ambient space-charge density being emitted to form this wind is bunched. These considerations immediately place the coherent radio frequency luminosity from such bunches near 10 to the 28th erg/s for typical pulsar parameters. A comparable amount of incoherent radiation is emitted for typical (1 second) pulsars. For very rapid pulsars, however, the latter component grows more rapidly than the available energy sources. The comparatively low radio luminosity of the Crab and Vela pulsars is attributed to both components being limited in the same ratio. The incoherent radiation essentially has a synchotron spectrum and extends to gamma-ray energies; consequently the small part of the total luminosity that is at optical wavelengths is unobservable. Assuming full coherence at all wavelengths short of a critical length gives a spectral index for the flux density of -8/3 at higher frequencies. The finite energy available from the injected particles would force the spectrum to roll over below about 100 MHz, although intrinsic morphological factors probably enter for any specific pulsar as well.

Michel, F. C.↗

Ground-state properties of the periodic Anderson model

The ground-state energy, hybridization matrix element, local moment, and spin-density correlations of a one-dimensional, finite-chain, periodic, symmetric Anderson model are obtained by numerical simulations and compared with perturbation theory and strong-coupling results. It is found that the local f-electron spins are compensated by correlation with other f-electrons as well as band electrons leading to a nonmagnetic ground state.

Blankenbecler, R.↗

MHD surface waves in high- and low-beta plasmas. I - Normal-mode solutions

Since the first paper by Barston (1964) on electrostatic oscillations in inhomogeneous cold plasmas, it has been commonly accepted that all finite layers with a continuous profile in pressure, density and magnetic field cannot support normal surface waves but instead the waves always decay through phase mixing (also called resonant absorption). The problem is reanalyzed by studying a compressible current sheet of a general structure with rotation of the magnetic field included. All inhomogeneous layers considered in the high-beta plasma limit do not support normal modes. However, in the limit of a low-beta plasma there are some cases when normal-mode solutions are recovered. The latter means that the process of resonant absorption is not common for all inhomogeneous layers.

Musielak, Z. E.↗

Pre-mixed flame simulations for non-unity Lewis numbers

A principal effect of turbulence on premixed flames in the flamelet region is to wrinkle the flame fronts. For non-unity Lewis numbers (Le), the local flame structure is altered in curved regions. This effect is examined using direct numerical simulations of the three dimensional, constant density, decaying isotropic turbulence with a single step, finite rate chemical reaction. Simulations of Lewis numbers 0.8, 1.0, and 1.2 are compared. The turbulent flame speed, S(sub T), increases as Le decreases. The correlation between S(sub T) and u prime found in previous Le = 1 simulations has a strong Lewis number dependency. The variance of the pdf of the flame curvature increases as Le decreases, indicating that the flames become more wrinkled. A strong correlation between local flame speed and curvature was found. For Le greater than 1, the flame speed increases in regions concave towards the products and decreases in convex regions. The opposite correlation was found for Le less than 1. The mean temperature of the products was also found to vary with Lewis number. For Le = 0.8, it is less than the adiabatic flame temperature and for Le = 1.2 it is greater.

Rutland, C. J.↗

Computational Design of Materials: Planetary Entry to Electric Aircraft and Beyond

NASA's projects and missions push the bounds of what is possible. To support the agency's work, materials development must stay on the cutting edge in order to keep pace. Today, researchers at NASA Ames Research Center perform multiscale modeling to aid the development of new materials and provide insight into existing ones. Multiscale modeling enables researchers to determine micro- and macroscale properties by connecting computational methods ranging from the atomic level (density functional theory, molecular dynamics) to the macroscale (finite element method). The output of one level is passed on as input to the next level, creating a powerful predictive model.

Materials Design↗

Dynamic fracture mechanics

Dynamic fracture and crack propagation concepts for ductile materials are reviewed. The equations for calculating dynamic stress integrity and the dynamic energy release rate in order to study dynamic crack propagation are provided. The stress intensity factor versus crack velocity relation is investigated. The uses of optical experimental techniques and finite element methods for fracture analyses are described. The fracture criteria for a rapidly propagating crack under mixed mode conditions are discussed; crack extension and fracture criteria under combined tension and shear loading are based on maximum circumferential stress or energy criteria such as strain energy density. The development and use of a Dugdale model and finite element models to represent crack and fracture dynamics are examined.

Kobayashi, A. S.↗

The calculation of electron density profiles from topside ionograms: Method and applications

A method for converting topside sounder ionograms into topside electron density profiles is discussed. The lamination method used is modified to take into account the variation of electron density and magnetic field within each lamination. Also included is a change of variable to produce a finite integrand of the integral involved, an iteration scheme that permits convergence on an initially unknown density, a second iteration scheme to overcome the problem of an uncertainty in the electron density at the satellite, and a modification to compensate for the changing satellite altitude over the duration of the sounding. Two applications of the technique are discussed: field-aligned traces for computing field-aligned profiles and computer-aided systems for scaling ionograms.

Lockwood, G. E. K.↗

Coupling finite element and integral equation solutions using decoupled boundary meshes

A method is outlined for calculating scattered fields from inhomogeneous penetrable objects using a coupled finite element-integral equation solution. The finite element equation can efficiently model fields in penetrable and inhomogeneous regions, while the integral equation exactly models fields on the finite element mesh boundary and in the exterior region. By decoupling the interior finite element and exterior integral equation meshes, considerable flexibility is found in both the number of field expansion points as well as their density. Only the nonmetal portions of the object need be modeled using a finite element expansion; exterior perfect conducting surfaces are modeled using an integral equation with a single unknown field since E(tan) is identically zero on these surfaces. Numerical convergence, accuracy, and stability at interior resonant frequencies are studied in detail.

Cwik, Tom↗

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗

Dislocation dynamics during the growth of silicon ribbon

The thermal viscoplastic stresses and the dislocation densities in silicon ribbon are computed for an axially changing thermal profile by using an iterative finite difference method. A material constitutive equation (Haasen-Sumino model) which involves an internal variable (mobile dislocation density) is used. The results are interpreted as showing that there is a maximum width of silicon ribbon that can be grown when viscoplasticity and dislocations are considered. This maximum width limitation does not exist if the material behavior is elastic.

Dillon, O. W., Jr.↗

Exploring Li-Air Batteries for High Specific Energy and High Power Applications: A Simulation Study

Commercialization of lithium-air batteries face many challenges, such as electrolyte decomposition, short cycle life, low energy efficiency, low power density, etc. However, commercialization of Li-air batteries for mass sensitive applications such as electric vehicles, portable power source, and drones is more challenging due to additional constraints of safety, electrolyte evaporation, high specific energy requirements, and reliable discharge times. In this presentation, we will present our finite element simulation results comparing Li-O2 and Li-air batteries using power density, energy density, specific power, and discharge times as metrics to evaluate different electrolytes and electrode geometry to reduce total mass and maximize discharge current. We use a finite element model and a discharge product model developed in which is based on porous electrode and concentrated electrolyte theories and the discharge product is modeled using quantum tunneling; for reaction kinetics and oxygen diffusion, an improved model was used. The electrolyte properties such as ion conductivity and ion diffusion were obtained from Molecular Dynamics (MD) simulations while the other parameters for the finite element model were calibrated to match experiments at high discharge current densities (>1.5 mA/cm2). The mass densities of different electrolytes were computed using MD simulations as well. For this presentation, we examine the practical electrochemical mass of a system at different current ratings, the sensitivity of mass to the use of ambient air as compared to pure oxygen as well as the electrolyte, which affects maximum current density and the total mass associated with the electrolyte (which includes the mass of additional components), and, the optimization of battery geometry for total discharge time, average discharge voltage, maximum discharge current, and minimum electrochemical mass.

Mehta, M.↗

Numerical investigations of low-density nozzle flow by solving the Boltzmann equation

A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.

Deng, Zheng-Tao↗

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗

Fluid dynamical description of relativistic nuclear collisions

On the basis of both a conventional relativistic nuclear fluid dynamic model and a two fluid generalization that takes into account the interpenetration of the target and projectile upon contact, collisions between heavy nuclei moving at relativistic speeds are calculated. This is done by solving the relevant equations of motion numerically in three spatial dimensions by use of particle in cell finite difference computing techniques. The effect of incorporating a density isomer, or quasistable state, in the nuclear equation of state at three times normal nuclear density, and the effect of doubling the nuclear compressibility coefficient are studied. For the reaction 20Ne + 238U at a laboratory bombarding energy per nucleon of 393 MeV, the calculated distributions in energy and angle of outgoing charged particles are compared with recent experimental data both integrated over all impact parameters and for nearly central collisions.

Nix, J. R.↗