Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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 721 records · Page 40

The electromagnetic ion cyclotron instability in the Io torus

The electromagnetic ion cyclotron waves which are expected to exist in the Jovian magnetosphere are investigated. The temperature anisotropy generated by the inward radial diffusion of hot ions gives rise to an instability of L mode waves in the off-equatorial region of the Io torus. The resulting pitch angle scattering has been suggested as the cause for the precipitation of ions into the loss cone and auroral excitation. The linear wave dispersion is first examined, and the nonlinear wave amplitude for the saturated state is studied. Two estimates of the wave saturation level are checked by performing an electromagnetic hybrid simulation. Estimated nonlinear saturation amplitudes are compared with those resulting from linear amplification in a finite length. The result shows that the waves in the Jovian magnetosphere produced by the hot protons are mostly in a linear regime.

Machida, S.↗

Modeling the Diffuse Cloud-Top Optical Emissions from Ground and Cloud Flashes

A number of studies have indicated that the diffuse cloud-top optical emissions from intra-cloud (IC) lightning are brighter than that from normal negative cloud-to-ground (CG) lightning, and hence would be easier to detect from a space-based sensor. The primary reason provided to substantiate this claim has been that the IC is at a higher altitude within the cloud and therefore is less obscured by the cloud multiple scattering medium. CGs at lower altitudes embedded deep within the cloud are more obscured, so CG detection is thought to be more difficult. However, other authors claim that because the CG source current (and hence luminosity) is typically substantially larger than IC currents, the greater CG source luminosity is large enough to overcome the effects of multiple scattering. These investigators suggest that the diffuse cloud top emissions from CGs are brighter than from ICs, and hence are easier to detect from space. Still other investigators claim that the detection efficiency of CGs and ICs is about the same because modern detector sensitivity is good enough to "see" either flash type no matter which produces a brighter cloud top emission. To better assess which of these opinions should be accepted, we introduce an extension of a Boltzmann lightning radiative transfer model previously developed. It considers characteristics of the cloud (geometry, dimensions, scattering properties) and specific lightning channel properties (length, geometry, location, current, optical wave front propagation speed/direction). As such, it represents the most detailed modeling effort to date. At least in the few cases studied thus far, it was found that IC flashes appear brighter at cloud top than the lower altitude negative ground flashes, but additional model runs are to be examined before finalizing our general conclusions.

Solakiewicz, Richard↗

A study of four soft X-ray enhancements

Four enhanced regions of diffuse X-ray emission are studied in the 0.15-2 KeV band. Two regions are consistent with clouds of plasma with a temperature of 1 million K about 50 pc in diameter at a distance of 150 pc. These clouds require a particle density of only a factor of 3 above an ambient hot component which may fill a substantial portion of the galactic disk in the vicinity of the sun. Another region near Eridanus is consistent with a low integrated amount of atomic hydrogen along the direction of observation, thereby permitting a greater path length of plasma to contribute to the observed flux. Finally, an extended region in Centaurus shows a temperature of about 20 million K and may be an old supernova remnant expanding into a low-density region of the galactic disk.

Long, K. S.↗

Forced Convection and Sedimentation Past a Flat Plate

The steady laminar flow of a well-mixed suspension of monodisperse solid spheres, convected steadily past a horizontal flat plate and sedimenting under the action of gravity, is examined. It is shown that, in the limit as Re approaches infinity and epsilon approaches 0, where Re is the bulk Reynolds number and epsilon is the ratio of the particle radius a to the characteristic length scale L, the analysis for determining the particle concentration profile has several aspects in common with that of obtaining the temperature profile in forced-convection heat transfer from a wall to a fluid stream moving at high Reynolds and Prandtl numbers. Specifically, it is found that the particle concentration remains uniform throughout the O(Re(exp -1/2)) thick Blasius boundary layer except for two O(epsilon(exp 2/3)) thin regions on either side of the plate, where the concentration profile becomes non-uniform owing to the presence of shear-induced particle diffusion which balances the particle flux due to convection and sedimentation. The system of equations within this concentration boundary layer admits a similarity solution near the leading edge of the plate, according to which the particle concentration along the top surface of the plate increases from its value in the free stream by an amount proportional to X(exp 5/6), with X measuring the distance along the plate, and decreases in a similar fashion along the underside. But, unlike the case of gravity settling on an inclined plate in the absence of a bulk flow at infinity considered earlier, here the concentration profile remains continuous everywhere. For values of X beyond the region near the leading edge, the particle concentration profile is obtained through the numerical solution of the relevant equations. It is found that, as predicted from the similarity solution, there exists a value of X at which the particle concentration along the top side of the plate attains its maximum value phi(sub m) and that, beyond this point, a stagnant sediment layer will form that grows steadily in time. This critical value of X is computed as a function of phi(sub s), the particle volume fraction in the free stream. In contrast, but again in conformity with the similarity solution, for values of X sufficiently far removed from the leading edge along the underside of the plate, a particle-free region is predicted to form adjacent to the plate. This model, with minor modifications, can be used to describe particle migration in other shear flows, as, for example, in the case of crossflow microfiltration.

Pelekasis, Nikolaos A.↗

Flow-induced gelation of living (micellar) polymers

The effect of shear velocity gradients on the size (L) of rodlike micelles in dilute and semidilute solution is considered. A kinetic equation is introduced for the time-dependent concentration of aggregates of length L, consisting of 'bimolecular' combination processes L + L-prime yield (L + L-prime) and unimolecular fragmentations L yield L + (L - L-prime). The former are described by a generalization (from spheres to rods) of the Smoluchowski mechanism for shear-induced coalesence of emulsions, and the latter by incorporating the tension-deformation effects due to flow. Steady-state solutions to the kinetic equation are obtained, with the corresponding mean micellar size evaluated as a function of the Peclet number P (i.e., the dimensionless ratio of the flow rate and the rotational diffusion coefficient). For sufficiently dilute solutions, only a weak dependence of the micellar size on P is found. In the semidilute regime, however, an apparent divergence in the micellar size at P of about 1 suggests a flow-induced first-order gelation phenomenon.

Bruinsma, Robijn↗

Transitional Gas Jet Diffusion Flames in Microgravity

Drop tower experiments were performed to identify buoyancy effects in transitional hydrogen gas jet diffusion flames. Quantitative rainbow schlieren deflectometry was utilized to optically visualize the flame and to measure oxygen concentration in the laminar portion of the flame. Test conditions consisted of atmospheric pressure flames burning in quiescent air. Fuel from a 0.3mm inside diameter tube injector was issued at jet exit Reynolds numbers (Re) of 1300 to 1700. Helium mole percentage in the fuel was varied from 0 to 40%. Significant effects of buoyancy were observed in near field of the flame even-though the fuel jets were momentum-dominated. Results show an increase of breakpoint length in microgravity. Data suggest that transitional flames in earth-gravity at Re<1300 might become laminar in microgravity.

Agrawal, Ajay K.↗

Catalytic ignition model in a monolithic reactor with in-depth reaction

Two transient models have been developed to study the catalytic ignition in a monolithic catalytic reactor. The special feature in these models is the inclusion of thermal and species structures in the porous catalytic layer. There are many time scales involved in the catalytic ignition problem, and these two models are developed with different time scales. In the full transient model, the equations are non-dimensionalized by the shortest time scale (mass diffusion across the catalytic layer). It is therefore accurate but is computationally costly. In the energy-integral model, only the slowest process (solid heat-up) is taken as nonsteady. It is approximate but computationally efficient. In the computations performed, the catalyst is platinum and the reactants are rich mixtures of hydrogen and oxygen. One-step global chemical reaction rates are used for both gas-phase homogeneous reaction and catalytic heterogeneous reaction. The computed results reveal the transient ignition processes in detail, including the structure variation with time in the reactive catalytic layer. An ignition map using reactor length and catalyst loading is constructed. The comparison of computed results between the two transient models verifies the applicability of the energy-integral model when the time is greater than the second largest time scale of the system. It also suggests that a proper combined use of the two models can catch all the transient phenomena while minimizing the computational cost.

Tien, Ta-Ching↗

Measurement of Flaw Size From Thermographic Data

Simple methods for reducing the pulsed thermographic responses of delaminations tend to overestimate the size of the delamination, since the heat diffuses in the plane parallel to the surface. The result is a temperature profile over the delamination which is larger than the delamination size. A variational approach is presented for reducing the thermographic data to produce an estimated size for a flaw that is much closer to the true size of the delamination. The method is based on an estimate for the thermal response that is a convolution of a Gaussian kernel with the shape of the flaw. The size is determined from both the temporal and spatial thermal response of the exterior surface above the delamination and constraints on the length of the contour surrounding the delamination. Examples of the application of the technique to simulation and experimental data are presented to investigate the limitations of the technique.

Winfree, William P.↗

Analysis of localized fringes in the holographic optical Schlieren system

The relation between localization of interference fringes in classical and holographic interferometry is reviewed and an application of holographic interferometry is considered for which the object is a transparent medium with nonhomogeneous refractive index. The technique is based on the analysis of the optical path length change of the object wave as it propagates through a transparent medium. Phase shifts due to variations of the speed of light within the medium give rise to an interference pattern. The resulting interferogram can be used to determine the physical properties of the medium or transparent object. Such properties include the mass density of fluids, electron densities of plasmas, the temperature of fluids, the chemical species concentration of fluids, and the state of stress in solids. The optical wave used can be either a simple plane or spherical wave, or it may be a complicated spatial wave scattered by a diffusing screen. The mathematical theory on the formation and analysis of localized fringes, the general theoretical concepts used, and a computer code for analysis are included along with the inversion of fringe order data.

Kurtz, R. L.↗

The Effects of pH on the Growth and Aspect Ratio of Chicken Egg White Lysozyme Crystals Prepared in Different Buffers

We have measured the nucleation and aspect ratio of CEWL crystals grown by vapor diffusion in acetate, butyrate, carbonate, succinate, and phosphate buffers in a range of pH spanning the pK(sub a) of these buffers. The nucleation numbers drop off significantly in the vicinity of pK(sub a) for each of the buffers except the phosphate system, in which we used only the pH range around the second titration point(pK2). There is a concomitant increase in the sizes of the crystals. Some typical nucleation number results are shown. These data support and extend other observations. In addition, we have examined changes in aspect ratio which accompany the suppression of nucleation within each buffer system. The length of the face in the [001] direction was measured, and compared to the width of the (110) face in the [110] type directions. We find that while the aspect ratio of the crystals is affected by pH, it is dominated by a correlation with the size of the crystals. Small crystals are longer in the [0011 direction than crystals that are larger (higher pH within a buffer system). This relationship is found to hold independent of the choice of buffer. These results are consistent with those of Judge et al, who used a batch process which resulted in uniform sizing of crystals at each pH. In these experiments, we specifically avoid agitating the protein/salt buffer mixture when combining the two. This permits the formation of a range of sizes at a given pH. The results for a .05 M acetate 5% NaCl buffer are also shown. We will discuss these results in light of a growth model.

Gibson, U. J.↗

A Tabletop Tool for Modeling Life Support Systems

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations.

Ramachandran, N.↗

A turbulence module for the NPARC code

A turbulence module is developed for the 2D version of the NPARC code which is currently restricted to planar or axisymmetric flows without swirling. Four turbulence models have been built into the module: Baldwin-Lomax, Chien, Shih-Lumley and CMOTT models. The first is a mixing-length eddy-viscosity model which is mainly used for initialization of computational fields and the last three are the low Reynolds number two-equation models. Unlike chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +), an advantage for separated flow calculations. Contrary to the NPARC and most other compressible codes, the non-delta form of transport equations is used which leads to a simpler linearization and is more effective than using the delta form in ensuring the positiveness of the turbulent kinetic energy and its dissipation rate. To reduce numerical diffusion while maintaining necessary stability, a second-order accurate and bounded scheme is used for the convective terms of the turbulent transport equations. This scheme is implemented in a deferred correction manner so that the main coefficients of the resulting difference equations are always positive, thus making the numerical solutions process unconditionally stable. The system of equations is solved via a decoupled method and by the alternating direction TDMA of Thomas. The module can be easily linked to the NPARC code for turbulent flow calculations.

Zhu, J.↗

Combustion of Metals in Reduced-Gravity and Extraterrestrial Environments

As a result of the ongoing exploration of Mars and the several unmanned and possibly manned missions planned for the near future, increased attention has been given to the use of the natural resources of the planet for rocket propellant production and energy generation. Since the atmosphere of Mars consists of approximately 95% carbon dioxide (CO2), this gas is the resource of choice to be employed for these purposes. Since many metals burn vigorously with CO2, these may be used as an energy source or as propellants for a research vehicle on the surface of Mars. Shafirovich and Goldshleger conducted experiments with spherical particles up to 2.5 mm in diameter and found that the burning process was controlled by diffusion and that the particles exhibited pulsating combustion due to superheating of the Mg vapor trapped inside a protective oxide shell. They also proposed a reaction mechanism based on the gas-phase reaction, Mg + CO2 yields MgO + CO and the heterogeneous reaction Mg + CO yields MgO + C occurring on the sample surface. In all the above studies with large Mg particles, the burning process is invariably influenced by strong convective currents that accelerate the combustion reaction and shorten the burning times. Although these currents are nearly absent in the burning of small particles, the high emissivity of the flames, rapid reaction, and small length scales make the gathering of any useful information on burning rates and flame structure very difficult. The goal of this investigation is to provide a detailed study of flame structure by taking advantage of large, free-floating spherical metal samples and their corresponding long burning times available in a weightless environment. The use of reduced gravity is essential to eliminate the intrusive buoyant flows that plague high temperature metal reactions, to remove the destructive effect of gravity on the shape of molten metal samples, and to study the combustion behavior of metals in the presence of solid oxides undisturbed by natural convection. This work presents the most complete modeling of metal particle burning to date for Mg with CO2 and O2.

Branch, M. C.↗

Space Laboratory on a Table Top: A Next Generative ECLSS design and diagnostic tool

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale-time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations. Preliminary experimental and modeling work in this area will be presented. This involves testing of a single inlet-exit model with detailed 3-D flow visualization and quantitative diagnostics and computational modeling of the system.

Ramachandran, N.↗

Chemical application of diffusion quantum Monte Carlo

The diffusion quantum Monte Carlo (QMC) method gives a stochastic solution to the Schroedinger equation. This approach is receiving increasing attention in chemical applications as a result of its high accuracy. However, reducing statistical uncertainty remains a priority because chemical effects are often obtained as small differences of large numbers. As an example, the single-triplet splitting of the energy of the methylene molecule CH sub 2 is given. The QMC algorithm was implemented on the CYBER 205, first as a direct transcription of the algorithm running on the VAX 11/780, and second by explicitly writing vector code for all loops longer than a crossover length C. The speed of the codes relative to one another as a function of C, and relative to the VAX, are discussed. The computational time dependence obtained versus the number of basis functions is discussed and this is compared with that obtained from traditional quantum chemistry codes and that obtained from traditional computer architectures.

Reynolds, P. J.↗

The Chemical Composition and Gas-to-Dust Mass Ratio of Nearby Interstellar Matter

We use recent results on interstellar gas toward nearby stars and interstellar by-products within the solar system to select among the equilibrium radiative transfer models of the nearest interstellar material presented in Slavin & Frisch. For the assumption that O/H - 400 parts per million, models 2 and 8 are found to yield good fits to available data on interstellar material inside and outside of the heliosphere, with the exception of the Ne abundance in the pickup ion and anomalous cosmic-ray populations. For these models, the interstellar medium (ISM) at the entry point to the heliosphere has n(H(sup 0)) = 0.202-0.208/cu cm, n(He(sup 0) = 0.0137-0.0152/cu cm, and ionizations X(H) = 0.29-0.30, X(He) = 0.47-0.51. These best models suggest that the chemical composition of the nearby ISM is approx.60%-70% subsolar if S is undepleted. Both H(0) and H(+) need to be included when evaluating abundances of ions found in warm diffuse clouds. Models 2 and 8 yield an H filtration factor of approx.0.46. Gas-to-dust mass ratios for the ISM toward epsilon CMa are R(sub gd) = 178-183 for solar abundances of Holweger or R(sub gd) = 611-657 for an interstellar abundance standard 70% solar. Direct observations of dust grains in the solar system by Ulysses and Galileo yield R(sub gd) appr0x. 115 for models 2 and 8, supporting earlier results (Frisch and coworkers). If the local ISM abundances are subsolar, then gas and dust are decoupled over small spatial scales. The inferred variation in R(sub gd) over parsec length scales is consistent with the fact that the ISM near the Sun is part of a dynamically active cluster of cloudlets flowing away from the Sco-Cen association. Observations toward stars within approx.500 pc show that R(sub gd) correlates with the percentage of the dust mass that is carried by iron, suggesting that an Fe-rich grain core (by mass) remains after grain destruction. Evidently large dust grains (>10(exp -13) g) and small dust grains (<10(exp -13) g) are not well mixed over parsec length spatial scales in the ISM. It also appears that very small C-dominated dust grains have been destroyed in the ISM within several parsecs of the Sun, since C appears to be essentially undepleted. However, if gas-dust coupling breaks down over the cloud lifetime, the missing mass arguments applied here to determine R(sub gd) and dust grain mineralogy are not appropriate.

Frisch, Priscilla C.↗

Quantum Mechanical Study of Nanoscale MOSFET

The steady state characteristics of MOSFETS that are of practical Interest are the drive current, off-current, dope of drain current versus drain voltage, and threshold voltage. In this section, we show that quantum mechanical simulations yield significantly different results from drift-diffusion based methods. These differences arise because of the following quantum mechanical features: (I) polysilicon gate depletion in a manner opposite to the classical case (II) dependence of the resonant levels in the channel on the gate voltage, (III) tunneling of charge across the gate oxide and from source to drain, (IV) quasi-ballistic flow of electrons. Conclusions dI/dV versus V does not increase in a manner commensurate with the increase in number of subbands. - The increase in dI/dV with bias is much smaller then the increase in the number of subbands - a consequence of bragg reflection. Our calculations show an increase in transmission with length of contact, as seen in experiments. It is desirable for molecular electronics applications to have a small contact area, yet large coupling. In this case, the circumferential dependence of the nanotube wave function dictates: - Transmission in armchair tubes saturates around unity - Transmission in zigzag tubes saturates at two.

Svizhenko, Alexei↗

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗