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 217 records · Page 12

Finite Element Modeling of the Posterior Eye in Microgravity

Microgravity experienced during spaceflight affects astronauts in various ways, including weakened muscles and loss of bone density. Recently, visual impairment and intracranial pressure (VIIP) syndrome has become a major concern for space missions lasting longer than 30 days. Astronauts suffering from VIIP syndrome have changes in ocular anatomical and visual impairment that persist after returning to earth. It is hypothesized that a cephalad fluid shift in microgravity may increase the intracranial pressure (ICP), which leads to an altered biomechanical environment of the posterior globe and optic nerve sheath (ONS).Currently, there is a lack of knowledge of how elevated ICP may lead to vision impairment and connective tissue changes in VIIP. Our goal was to develop a finite element model to simulate the acute effects of elevated ICP on the posterior eye and optic nerve sheath. We used a finite element (FE) analysis approach to understand the response of the lamina cribrosa and optic nerve to the elevations in ICP thought to occur in microgravity and to identify which tissue components have the greatest impact on strain experienced by optic nerve head tissues.

intraocular pressure↗

Phase diagram of the two-dimensional negative-U Hubbard model

Theoretical arguments and numerical calculations are used to discuss the phase diagram of the two-dimensional negative-U Hubbard model. The results are consistent with (1) a vanishing transition temperature at half-filling but with a ground state having both superconducting and charge-density-wave long-range order, and (2) a Kosterlitz-Thouless transition at a finite temperature into a superconducting state with power-law decay of the pairing correlations away from half-filling.

Scalettar, R. T.↗

Relative efficiency and accuracy of two Navier-Stokes codes for simulating attached transonic flow over wings

In the present study, two codes which solve the three-dimensional Thin-Layer Navier-Stokes (TLNS) equations are used to compute the steady-state flow for two test cases representing typical finite wings at transonic conditions. Several grids of C-O topology and varying point densities are used. After a description of each code and test case, standards for determining code efficiency and accuracy are defined and applied to determine the relative performance of the two codes in predicting turbulent transonic wing flows. Comparisons of computed surface pressure distributions with experimental data are made.

Bonhaus, Daryl L.↗

Relative efficiency and accuracy of two Navier-Stokes codes for simulating attached transonic flow over wings

Two codes which solve the 3-D Thin Layer Navier-Stokes (TLNS) equations are used to compute the steady state flow for two test cases representing typical finite wings at transonic conditions. Several grids of C-O topology and varying point densities are used to determine the effects of grid refinement. After a description of each code and test case, standards for determining code efficiency and accuracy are defined and applied to determine the relative performance of the two codes in predicting turbulent transonic wing flows. Comparisons of computed surface pressure distributions with experimental data are made.

Bonhaus, Daryl L.↗

Low-frequency ionization-driven instability of a discrete auroral arc model

The low-frequency (time scales of tens of seconds) dynamics of the auroral ionosphere differs from that of the nonauroral ionosphere by the presence of strong, unstable space- and time-dependent ionospheric ionization produced by precipitating auroral electrons. If recombination is relatively unimportant (as at high ionospheric heights), it is shown that, in general, transport processes cannot remove this ionization as fast as it is created, and no equilibrium is possible. These nonequilibrium phenomena are investigated in the context of a nonlinear adiabatic auroral model, which has previously been studied in static situations. A linearized local perturbation analysis is given of what amounts to a current-driven E x B gradient-drift instability with an ionization source, as well as some exact nonlinear solutions valid in a finite but limited spatial region. These solutions show continuing motion of auroral potential and plasma density, as the aurora tries to shift its ionization problems from one place to another. The analysis gives clues to the possibility of generation of chaos and of fine-scale spatial structure.

Cornwall, John M.↗

Probability distributions for the magnification of quasars due to microlensing

Gravitational microlensing can magnify the flux of a lensed quasar considerably and therefore possibly influence quasar source counts or the observed quasar luminosity function. A large number of distributions of magnification probabilities due to gravitational microlensing for finite sources are presented, with a reasonable coverage of microlensing parameter space (i.e., surface mass density, external shear, mass spectrum of lensing objects). These probability distributions were obtained from smoothing two-dimensional magnification patterns with Gaussian source profiles. Different source sizes ranging from 10 exp 14 cm to 5 x 10 exp 16 cm were explored. The probability distributions show a large variety of shapes. Coefficients of fitted slopes for large magnifications are presented.

Wambsganss, Joachim↗

Molecular Static Third-Order Polarizabilities of Carbon-Cage Fullerene and Their Correlation with Three Geometric Properties: Symmetry, Aromaticity, and Size

The static third-order polarizabilities (gamma) of C60, C70, five isomers of C78 and two isomers of C84 were analyzed in terms of three properties, from a geometric point of view: symmetry, aromaticity and size. The polarizability values were based on the finite field approximation using a semiempirical Hamiltonian (AM1) and applied to molecular structures obtained from density functional theory calculations. Symmetry was characterized by the molecular group order. The selection of 6-member rings as aromatic was determined from an analysis of bond lengths. Maximum interatomic distance and surface area were the parameters considered with respect to size. Based on triple linear regression analysis, it was found that the static linear polarizability (alpha) and gamma in these molecules respond differently to geometrical properties: alpha depends almost exclusively on surface area while gamma is affected by a combination of number of aromatic rings, length and group order, in decreasing importance. In the case of alpha, valence electron contributions provide the same information as all-electron estimates. For gamma, the best correlation coefficients are obtained when all-electron estimates are used and when the dependent parameter is ln(gamma) instead of gamma.

Moore, C. E.↗

Theoretical Studies of Liquid He-4 Near the Superfluid Transition

We performed theoretical studies of liquid helium by applying state of the art simulation and finite-size scaling techniques. We calculated universal scaling functions for the specific heat and superfluid density for various confining geometries relevant for experiments such as the confined helium experiment and other ground based studies. We also studied microscopically how the substrate imposes a boundary condition on the superfluid order parameter as the superfluid film grows layer by layer. Using path-integral Monte Carlo, a quantum Monte Carlo simulation method, we investigated the rich phase diagram of helium monolayer, bilayer and multilayer on a substrate such as graphite. We find excellent agreement with the experimental results using no free parameters. Finally, we carried out preliminary calculations of transport coefficients such as the thermal conductivity for bulk or confined helium systems and of their scaling properties. All our studies provide theoretical support for various experimental studies in microgravity.

Manousakis, Efstratios↗

Strip Yield Model Numerical Application to Different Geometries and Loading Conditions

A new numerical method based on the strip-yield analysis approach was developed for calculating the Crack Tip Opening Displacement (CTOD). This approach can be applied for different crack configurations having infinite and finite geometries, and arbitrary applied loading conditions. The new technique adapts the boundary element / dislocation density method to obtain crack-face opening displacements at any point on a crack, and succeeds by obtaining requisite values as a series of definite integrals, the functional parts of each being evaluated exactly in a closed form.

Hatamleh, Omar↗

High Enthalpy Studies of Capsule Heating in an Expansion Tunnel Facility

Measurements were made on an Orion heat shield model to demonstrate the capability of the new LENS-XX expansion tunnel facility to make high quality measurements of heat transfer distributions at flow velocities from 3 km/s (h(sub 0) = 5 MJ/kg) to 8.4 km/s (h(sub 0) = 36 MJ/kg). Thirty-nine heat transfer gauges, including both thin-film and thermocouple instruments, as well as four pressure gauges, and high-speed Schlieren were used to assess the aerothermal environment on the capsule heat shield. Only results from laminar boundary layer runs are reported. A major finding of this test series is that the high enthalpy, low-density flows displayed surface heating behavior that is observed to be consistent with some finite-rate recombination process occurring on the surface of the model. It is too early to speculate on the nature of the mechanism, but the response of the gages on the surface seems generally repeatable and consistent for a range of conditions. This result is an important milestone in developing and proving a capability to make measurements in a ground test environment and extrapolate them to flight for conditions with extreme non-equilibrium effects. Additionally, no significant, isolated stagnation point augmentation ("bump") was observed in the tests in this facility. Cases at higher Reynolds number seemed to show the greatest amount of overall increase in heating on the windward side of the model, which may in part be due to small-scale particulate.

Dufrene, Aaron↗

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗

An efficient sampling technique for sums of bandpass functions

A well known sampling theorem states that a bandlimited function can be completely determined by its values at a uniformly placed set of points whose density is at least twice the highest frequency component of the function (Nyquist rate). A less familiar but important sampling theorem states that a bandlimited narrowband function can be completely determined by its values at a properly chosen, nonuniformly placed set of points whose density is at least twice the passband width. This allows for efficient digital demodulation of narrowband signals, which are common in sonar, radar and radio interferometry, without the side effect of signal group delay from an analog demodulator. This theorem was extended by developing a technique which allows a finite sum of bandlimited narrowband functions to be determined by its values at a properly chosen, nonuniformly placed set of points whose density can be made arbitrarily close to the sum of the passband widths.

Lawton, W. M.↗

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗

AC impedance analysis of polypyrrole thin films

The AC impedance spectra of thin polypyrrole films were obtained at open circuit potentials from -0.4 to 0.4 V vs SCE. Two limiting cases are discussed for which simplified equivalent circuits are applicable. At very positive potentials, the predominantly nonfaradaic AC impedance of polypyrrole is very similar to that observed previously for finite porous metallic films. Modeling of the data with the appropriate equivalent circuit permits effective pore diameter and pore number densities of the oxidized film to be estimated. At potentials from -0.4 to -0.3 V, the polypyrrole film is essentially nonelectronically conductive and diffusion of polymer oxidized sites with their associated counterions can be assumed to be linear from the film/substrate electrode interface. The equivalent circuit for the polypyrrole film at these potentials is that previously described for metal oxide, lithium intercalation thin films. Using this model, counterion diffusion coefficients are determined for both semi-infinite and finite diffusion domains. In addition, the limiting low frequency resistance and capacitance of the polypyrrole thin fims was determined and compared to that obtained previously for thicker films of the polymer. The origin of the observed potential dependence of these low frequency circuit components is discussed.

Penner, Reginald M.↗

Analysis of redox additive-based overcharge protection for rechargeable lithium batteries

The overcharge condition in secondary lithium batteries employing redox additives for overcharge protection, has been theoretically analyzed in terms of a finite linear diffusion model. The analysis leads to expressions relating the steady-state overcharge current density and cell voltage to the concentration, diffusion coefficient, standard reduction potential of the redox couple, and interelectrode distance. The model permits the estimation of the maximum permissible overcharge rate for any chosen set of system conditions. Digital simulation of the overcharge experiment leads to numerical representation of the potential transients, and estimate of the influence of diffusion coefficient and interelectrode distance on the transient attainment of the steady state during overcharge. The model has been experimentally verified using 1,1-prime-dimethyl ferrocene as a redox additive. The analysis of the experimental results in terms of the theory allows the calculation of the diffusion coefficient and the formal potential of the redox couple. The model and the theoretical results may be exploited in the design and optimization of overcharge protection by the redox additive approach.

Narayanan, S. R.↗

Source term evaluation for combustion modeling

A modification is developed for application to the source terms used in combustion modeling. The modification accounts for the error of the finite difference scheme in regions where chain-branching chemical reactions produce exponential growth of species densities. The modification is first applied to a one-dimensional scalar model problem. It is then generalized to multiple chemical species, and used in quasi-one-dimensional computations of shock-induced combustion in a channel. Grid refinement studies demonstrate the improved accuracy of the method using this modification. The algorithm is applied in two spatial dimensions and used in simulations of steady and unsteady shock-induced combustion. Comparisons with ballistic range experiments give confidence in the numerical technique and the 9-species hydrogen-air chemistry model.

Sussman, Myles A.↗

Finite linear diffusion model for design of overcharge protection for rechargeable lithium batteries

The overcharge condition in secondary lithium batteries employing redox additives for overcharge protection has been theoretically analyzed in terms of a finite linear diffusion model. The analysis leads to expressions relating the steady-state overcharge current density and cell voltage to the concentration, diffusion coefficient, standard reduction potential of the redox couple, and interelectrode distance. The model permits the estimation of the maximum permissible overcharge rate for any chosen set of system conditions. The model has been experimentally verified using 1,1-prime-dimethylferrocene as a redox additive. The theoretical results may be exploited in the design and optimization of overcharge protection by the redox additive approach.

Narayanan, S. R.↗

Plasma dynamics driven by finite-width current filament and KV potential drops in ionosphere-magnetosphere coupling

We have simulated the plasma dynamics of a quasi-two-dimensional current filament in ionosphere-magnetosphere coupling. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The dynamical interactions of the individual current filaments with the ionospheric and magnetospheric plasma generate kV field-aligned potential drops along the field lines, and the side-by-side display exhibits the evolution of the implied potential structure in the horizontal direction. The primary cross-field motion produced by the horizontal E-field (ExB drift) is perpendicular to both of the significant spatial directions, and is thus ignorable in this geometry. The effects of other cross-field drict processes are discussed. The simulation thus provides insight into the dynamical evolution of 2D potential structures driven by an imposed finite-width field-aligned current profile. Spatial and temporal variation of other plasma parameters (temperatures, density, etc.) that are determined primarily by parallel transport effects are also displayed.

Ganguli, Supriya B.↗