Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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.

208 records · Page 12

Kinetic plasma processes occurring in the outer plasmasphere

One area of data analysis work that was begun under this contract is the fitting of the perpendicular velocity distributions of equatorially trapped ions with a Kappa function. This type of characterization of the trapped ions will be very useful for comparison with velocity distributions produced by the model. A second area of data analysis is to study data from consecutive passes when DE 1's apogee was near the magnetic equator and the spacecraft was often skimming along nearly the same L shell. In 1982 three such periods occurred in May, June, and July. For these consecutive events we have Kp histories, density measurements from a number of sources (Whistler data, DE SFR, ISEE SFR) and consecutive samples of ion pitch angle distributions along field lines. It is clear from this data how the pitch angle distributions evolve during a flux tube refilling event. Our modeling of the flow of plasma along closed field lines is following two basic tracks. The first is a study of the basic refilling process without the effect of wave-particle heating near the equator or the effect of large or abrupt field-aligned electric potential drops. This model includes the effects of Coulomb self-collisions and collisions with the O+ ions in the topside ionosphere. The second track is a study of the effects of wave produced pitch-angle scattering and perpendicular heating occurring near the magnetic equator, in connection with the development of large potential drops that result from electron heating and the development of density gradients.

Wilson, Gordon R.↗

The Effects of Ground and Space Processing on the Properties of Organic, Polymeric, and Colloidal Materials

In recent years, a great deal of interest has been directed toward the use of organic materials in the development of high-efficiency optoelectronic and phototonic devices. There is a myriad of possibilities among organic materials which allow flexibility in the design of unique structures with a variety of functional groups. The use of nonlinear optical (NLO) organic materials as thin film wave-guides allows full exploitation of their desirable qualifies by permitting long interaction lengths and large susceptibilities allowing modest power input. There are several methods in use to prepare thin films such as Langmuir-Blodgett (LB) and self-assembly techniques, vapor deposition, growth from sheared solution or melt, and melt growth between glass plates. Organic-based materials have many features that make them desirable for use in optical devices, such as high second-and third-order nonlinearity, flexibility of molecular design, and damage resistance to optical radiation. However, their use in devices has been hindered by processing difficulties for crystals and thin films. We discuss the potential role of microgravity processing of a few organic and polymeric materials. It is of interest to note how materials with second-and third-order NLO behavior may be improved in a diffusion-limited environment and ways in which convection may be detrimental to these materials. We focus our discussion on third-order materials for all-optical switching, and second-order materials for frequency conversion and electrooptics. The goal of minimizing optical loss obviously depends on processing methods. For solution-based processes, such as solution crystal growth and solution photopolymerization, it is well known that thermal and solutal density gradients can initiate buoyancy-driven convection. Resultant fluid flows can affect transport of material to and from growth interfaces and become manifest in the morphology and homogeneity of the growing film or crystal. Likewise, buoyancy-driven convection can hinder production of defect-free, high-quality crystals or films during crystal and film growth by vapor deposition.

Frazier, Donald O.↗

Hydrodynamics of Two-Phase Flows through Porous Media in Microgravity: Packed Bed Reactor Experiment onboard of the International Space Station

The objective of the Packed Bed Reactor Experiment was to develop a fundamental understanding of the hydrodynamics of two-phase flow with no phase transition through porous media in microgravity. This work describes two experimental campaigns flown on the International Space Station (ISS) that were designed to achieve the objective. This work presents the results on flow patterns observed, two-phase flow pressure drop in porous media, and the impact of capillary effects on pressure drop at various gas and liquid flow rates. This work is the first to present predictive correlations of the two-phase friction factor for the different regimes identified based on the superficial liquid and gas velocity in microgravity, the first to hypothesize the different regimes based on the change in slope of the pressure gradient versus flow rate plots, the first to address the effects of the capillarity on the pressure gradient and the first to measure column holdup and assess the effective porosity based on the packing and the trapped gas bubbles. Experimental data on pressure drop in gas (N2) – liquid (water) flow show the presence of four different flow regimes in microgravity. Dispersed bubble flow (DB) and pulse (P) flow regimes are detected at high liquid flow rates, whereas at low liquid flow rates “large bubble” or elongated bubble (LB) and “gas channeling” (GC) regimes are observed. For these different flow regimes, different two-phase flow friction factor and pressure gradient correlations are presented for the first time as a function of the gas and liquid modified Reynolds numbers and Suratman number. Within the viscous-capillary (V-C) regime, it is found that the capillary contribution is the dominant force that contributes to the pressure drop for the wetting case (glass). However, for the non-wetting packing (Teflon), the viscous contribution dominates. It was found that the gas hold-up and pressure drop are functions of bed history at low liquid and gas flow rates with the magnitude of the hysteresis decreasing with increasing flow rates. PBRE-2 results show that the capillary force is a strong function of the superficial liquid velocity but is a much weaker function of the superficial gas velocity and varies inversely with the particle diameter. Within the Viscous-Capillary (V-C) regime, over 90% of the pressure gradient is attributed to the capillary contribution in the gas continuous regime. However, in the large bubble regime, the viscous and capillary contributions were comparable. After the completion of PBRE-2 with glass spherical beads, another bed packed with alumina was installed. Pressure gradient data for the alumina packed bed (PBRE-Water Recovery) were transferred to the PI and were not analyzed by the authors of this work. The Packed Bed Reactor Experiment concluded its on-orbit operation after two successful campaigns of testing in 2017 and in 2021. The flight hardware was brought back from orbit and is being reconditioned for a series of future experiments referred to as PBRE-WRS in support of water recovery (WRS). The WRS series consists of testing different two-phase fluid system components and packed beds to assess their pressure gradient characteristics, which will be used for designing packed bed reactors for various applications relevant to life support. Although there is so much relevance of the experimental results obtained from the two PBRE campaign in microgravity to life support, these findings also apply to other applications that involve two phase flow in porous medium such as fuel cells, transport of nutrient to plants in space and other chemical and materials processing that involve two phase flows. These systems operate differently in microgravity because, due to the lack of buoyancy, the density difference between the phases becomes irrelevant and no longer leads to phase separation.

Packed Bed↗

Responses of Precipitation and Hydrologic Processes to Tropical SST

The goal of the research is to identify the mechanisms in the response of tropical precipitation and atmospheric hydrologic cycle to sea surface temperature (SST) variability at seasonal-to-interannual time scales, and to utilize the knowledge for better understanding of climate feedback processes relevant to global change. As a first step to achieve the goal, we characterize the inter-relationship among convective/stratiform rain, ice/water clouds water vapor, and SST using TRMM satellite data and a cloud-resolving model. We examined the daily hydrologic variables [column water vapor (PW), cloud liquid water (CW), rainfall rates (RR)] as a function of SST using high-resolution data (0.25 x 0.25, daily) derived from TRMM satellite measurements. Comparing the winter of 97/98 (El Nino condition) against the winter of 99/00 (La Nina condition), area-mean values of all four hydrologic variables in cloudy areas within the tropical Pacific are higher in the El Nino winter than in the La Nina winter. This is consistent with previous observational analyses and SST warming experiments (idealized or ENSO-like) that showed the interaction between hydrologic cyclic and radiation at the seasonal to interannual time scales leads to intensified tropical circulation and hydrologic cycle. However, there is evidence that the enhanced hydrologic cycle over the warm pool is accompanied by an expansion of radiatively -driven subsidence in response to a stronger SST gradient between warm pool and surrounding cold pool. The expanding subsidence effectively reduces cloud amounts over the warm pool. Our analysis of daily variability further indicates a more vigorous water cycle characterized by higher PW, CW, and RR in response to overall warming. This is expected from the Clausius Clapeyron relation as a thermodynamic response to warming. However cloudy areas decrease in response to overall warming. This may be due to factors that are fundamentally different. One possibility is that in a warmer climate, there may be alteration in the microphysical processes in clouds, e.g., increased coalescence of cloud size droplets and increased entrainment of dry air from above due to the invigorated convective motions in response to overall warming. As a result, cloud amounts are less. Another interpretation is that, in colder air masses, more ice particles will be produced and, because of their lower densities, are less likely to precipitate, hence the higher cloud amount. This suggests that, when averaged over a sufficiently large area covering both the rising and sinking branches of the tropical circulation, there maybe an increase in the precipitation efficiency in stratiform clouds, in conjunction with increased subsidence in an enhanced water cycle. This is supported by the residence time at different rainrates.

Sui, C.-H.↗

Step Bunching: Influence of Impurities and Solution Flow

Step bunching results in striations even at relatively early stages of its development and in inclusions of mother liquor at the later stages. Therefore, eliminating step bunching is crucial for high crystal perfection. At least 5 major effects causing and influencing step bunching are known: (1) Basic morphological instability of stepped interfaces. It is caused by concentration gradient in the solution normal to the face and by the redistribution of solute tangentially to the interface which redistribution enhances occasional perturbations in step density due to various types of noise; (2) Aggravation of the above basic instability by solution flowing tangentially to the face in the same directions as the steps or stabilization of equidistant step train if these flows are antiparallel; (3) Enhanced bunching at supersaturation where step velocity v increases with relative supersaturation s much faster than linear. This v(s) dependence is believed to be associated with impurities. The impurities of which adsorption time is comparable with the time needed to deposit one lattice layer may also be responsible for bunching; (4) Very intensive solution flow stabilizes growing interface even at parallel solution and step flows; (5) Macrosteps were observed to nucleate at crystal corners and edges. Numerical simulation, assuming step-step interactions via surface diffusion also show that step bunching may be induced by random step nucleation at the facet edge and by discontinuity in the step density (a ridge) somewhere in the middle of a face. The corresponding bunching patterns produce the ones observed in experiment. The nature of step bunching generated at the corners and edges and by dislocation step sources, as well as the also relative importance and interrelations between mechanisms 1-5 is not clear, both from experimental and theoretical standpoints. Furthermore, several laws controlling the evolution of existing step bunches have been suggested, though unambiguous conclusions are still missing. Addressing these issues is the major goal of the present project. The theory addressing the above problem, experimental methods, several figures which include: (1) the spatial wave numbers at which the system is neutrally stable as a function of growth velocity for linear kinetics and supersaturation for nonlinear kinetics; (2) a schematic of the experiment of lysozyme crystal growing under conditions of natural convection; (3) fluctuations in time, t, of the normal growth rate, R(t), vicinal slope, p(t) and Fourier Spectra of R(t), discussions and conclusions are presented.

Chernov, A. A.↗

Paramagnetic Liquid Bridge in a Gravity-Compensating Magnetic Field

Magnetic levitation was used to stabilize cylindrical columns of a paramagnetic liquid in air between two solid supports. The maximum achievable length to diameter ratio R(sub max) was approx. (3.10 +/- 0.07), very close to the Rayleigh-Plateau limit of pi. For smaller R, the stability of the column was measured as a function of the Bond number, which could be continuously varied by adjusting the strength of the magnetic field. Liquid bridges supported by two solid surfaces have been attracting scientific attention since the time of Rayleigh and Plateau. For a cylindrical bridge of length L and diameter d, it was shown theoretically that in zero gravity the maximum slenderness ratio R (identically = L/d) is pi. The stability and ultimate collapse of such bridges is of interest because of their importance in a number of industrial processes and their potential for low gravity applications. In the presence of gravity, however, the cylindrical shape of an axisymmetric bridge tends to deform, limiting its stability and decreasing the maximum achievable value of R. Theoretical studies have discussed the stability and possible shapes of axisymmetric bridges. Experiments typically are performed in either a Plateau tank, in which the bridge is surrounded by a density-matched immiscible fluid, or in a space-borne microgravity environment. It has been shown, for example, that the stability limit R can be pushed beyond pi by using flow stabilization, by acoustic radiation pressure, or by forming columns in the presence of an axial electric field. In this work, magnetic levitation was used to simulate a low gravity environment and create quasi-cylindrical liquid columns in air. Use of a magnetic field permits us to continuously vary the Bond number B identically equal to (g)(rho)d(exp 2)/4(sigma), where g is the gravitational acceleration, rho is the density of the liquid, and sigma is the surface tension of the liquid in air. The dimensionless Bond number represents the relative importance of external forces acting on the liquid column to those due to surface tension. Our central result is that in a large magnetic field gradient we could create and stabilize columns of mixtures of water and paramagnetic manganese chloride tetrahydrate (MnCl2.4H2O), achieving a length to diameter ratio very close to pi.

Mahajan, Milind P.↗

Hierarchical Statistical 3D ' Atomistic' Simulation of Decanano MOSFETs: Drift-Diffusion, Hydrodynamic and Quantum Mechanical Approaches

When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain voltage only single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. A pilot version of a hydrodynamic 'atomistic' simulator has been developed in order to study the effect of the nonequilibrium, non local transport in decanano MOSFETs on the random dopant induced current fluctuations. For the first time we have also applied the density gradient approach in 3D to investigate the effect of the quantum confinement on the threshold voltage fluctuations. The developed 'atomistic' simulation techniques have been applied to study various fluctuation resistant MOSFET architectures including epitaxial and delta doped devices.

Asenov, Asen↗

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗

Feasibility of Actively Cooled Silicon Nitride Airfoil for Turbine Applications Demonstrated

Nickel-base superalloys currently limit gas turbine engine performance. Active cooling has extended the temperature range of service of nickel-base superalloys in current gas turbine engines, but the margin for further improvement appears modest. Therefore, significant advancements in materials technology are needed to raise turbine inlet temperatures above 2400 F to increase engine specific thrust and operating efficiency. Because of their low density and high-temperature strength and thermal conductivity, in situ toughened silicon nitride ceramics have received a great deal of attention for cooled structures. However, the high processing costs and low impact resistance of silicon nitride ceramics have proven to be major obstacles for widespread applications. Advanced rapid prototyping technology in combination with conventional gel casting and sintering can reduce high processing costs and may offer an affordable manufacturing approach. Researchers at the NASA Glenn Research Center, in cooperation with a local university and an aerospace company, are developing actively cooled and functionally graded ceramic structures. The objective of this program is to develop cost-effective manufacturing technology and experimental and analytical capabilities for environmentally stable, aerodynamically efficient, foreign-object-damage-resistant, in situ toughened silicon nitride turbine nozzle vanes, and to test these vanes under simulated engine conditions. Starting with computer aided design (CAD) files of an airfoil and a flat plate with internal cooling passages, the permanent and removable mold components for gel casting ceramic slips were made by stereolithography and Sanders machines, respectively. The gel-cast part was dried and sintered to final shape. Several in situ toughened silicon nitride generic airfoils with internal cooling passages have been fabricated. The uncoated and thermal barrier coated airfoils and flat plates were burner rig tested for 30 min without and with air cooling. Without cooling, the surface temperature of the flat plate reached approximately 2350 F. Starting with computer aided design (CAD) files of an airfoil and a flat plate with internal cooling passages, the permanent and removable mold components for gel casting ceramic slips were made by stereolithography and Sanders machines, respectively. The gel-cast part was dried and sintered to final shape. Several in situ toughened silicon nitride generic airfoils with internal cooling passages have been fabricated. The uncoated and thermal barrier coated airfoils and flat plates were burner rig tested for 30 min without and with air cooling. Without cooling, the surface temperature of the flat plate reached approximately 2350 F. With cooling, the surface temperature decreased to approximately 1910 F--a drop of approximately 440 F. This preliminary study demonstrates that a near-net-shape silicon nitride airfoil can be fabricated and that silicon nitride can sustain severe thermal shock and the thermal gradients induced by cooling and, thus, is a viable candidate for cooled components.

Bhatt, Ramakrishna T.↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗