Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

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

Kennedy, Christopher A.↗
NASA NTRS1997

Physiological Response of Plants Grown on Porous Ceramic Tubes

This research involves the manipulation of the root-zone water potential for the purposes of discriminating the rate limiting step in the inorganic nutrient uptake mechanism utilized by higher plants. This reaction sequence includes the pathways controlled by the root-zone conditions such as water tension and gradient concentrations. Furthermore, plant based control mechanisms dictated by various protein productions are differentiated as well. For the nutrients limited by the environmental availability, the kinetics were modeled using convection and diffusion equations. Alternatively, for the nutrients dependent upon enzyme manipulations, the uptakes are modeled using Michaelis-Menten kinetics. In order to differentiate between these various mechanistic steps, an experimental apparatus known as the Porous Ceramic Tube - Nutrient Delivery System (PCT-NDS) was used. Manipulation of the applied suction pressure circulating a nutrient solution through this system imposes a change in the matric component of the water potential. This compensates for the different osmotic components of water potential dictated by nutrient concentration. By maintaining this control over the root-zone conditions, the rate limiting steps in the uptake of the essential nutrients into tomato plants (Lycopersicon esculentum cv. Cherry Elite) were differentiated. Results showed that the uptake of some nutrients were mass transfer limited while others were limited by the enzyme kinetics. Each of these were adequately modeled with calculations and discussions of the parameter estimations provided.

Tsao, David↗
NASA NTRS1998

A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy

Stable and accurate interface conditions are derived for the linear advection-diffusion equation. The conditions are functionally independent of the spatial order of accuracy and rely only on the form of the discrete operator. We focus on high-order finite-difference operators that satisfy the summation-by-parts (SBP) property. We prove that stability is a natural consequence of the SBP operators used in conjunction with the new boundary conditions. In addition, we show that the interface treatments are conservative. New finite-difference operators of spatial accuracy up to sixth order are constructed: these operators satisfy the SBP property. Finite-difference operators are shown to admit design accuracy (p(sup th)-order global accuracy) when (p - 1)(sup th)-order stencil closures are used near the boundaries if the physical boundary conditions are implemented to at least p(sup th)-order accuracy. Stability and accuracy are demonstrated on the nonlinear Burgers' equation for an twelve-subdomain problem with randomly distributed interfaces.

Carpenter, Mark H.↗
NASA NTRS1997

A Thermal and Electrical Analysis of Power Semiconductor Devices

The state-of-art power semiconductor devices require a thorough understanding of the thermal behavior for these devices. Traditional thermal analysis have (1) failed to account for the thermo-electrical interaction which is significant for power semiconductor devices operating at high temperature, and (2) failed to account for the thermal interactions among all the levels involved in, from the entire device to the gate micro-structure. Furthermore there is a lack of quantitative studies of the thermal breakdown phenomenon which is one of the major failure mechanisms for power electronics. This research work is directed towards addressing. Using a coupled thermal and electrical simulation, in which the drift-diffusion equations for the semiconductor and the energy equation for temperature are solved simultaneously, the thermo-electrical interactions at the micron scale of various junction structures are thoroughly investigated. The optimization of gate structure designs and doping designs is then addressed. An iterative numerical procedure which incorporates the thermal analysis at the device, chip and junction levels of the power device is proposed for the first time and utilized in a BJT power semiconductor device. In this procedure, interactions of different levels are fully considered. The thermal stability issue is studied both analytically and numerically in this research work in order to understand the mechanism for thermal breakdown.

Vafai, Kambiz↗
NASA NTRS1996

Frequency Response of Pressure Sensitive Paints

An experimental method for measuring the frequency response of Pressure Sensitive Paints (PSP) is presented. These results lead to the development of a dynamic correction technique for PSP measurements which is of great importance to the advancement of PSP as a measurement technique. The ability to design such a dynamic corrector is most easily formed from the frequency response of the given system. An example of this correction technique is shown. In addition to the experimental data, an analytical model for the frequency response is developed from the one dimensional mass diffusion equation.

Winslow, Neal A.↗
NASA NTRS1999

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Advection-Diffusion Problems

Test problems are used to examine the performance of several one-dimensional numerical schemes based on the space-time conservation and solution element (CE/SE) method. Investigated in this paper are the CE/SE schemes constructed previously for solving the linear unsteady advection-diffusion equation and the schemes derived here for solving the nonlinear viscous and inviscid Burgers equations. In comparison with the numerical solutions obtained using several traditional finite-difference schemes with similar accuracy, the CE/SE solutions display much lower numerical dissipation and dispersion errors.

Wang, Xiao-Yen↗
NASA NTRS1998

A Parallel Non-Overlapping Domain-Decomposition Algorithm for Compressible Fluid Flow Problems on Triangulated Domains

This paper considers an algebraic preconditioning algorithm for hyperbolic-elliptic fluid flow problems. The algorithm is based on a parallel non-overlapping Schur complement domain-decomposition technique for triangulated domains. In the Schur complement technique, the triangulation is first partitioned into a number of non-overlapping subdomains and interfaces. This suggests a reordering of triangulation vertices which separates subdomain and interface solution unknowns. The reordering induces a natural 2 x 2 block partitioning of the discretization matrix. Exact LU factorization of this block system yields a Schur complement matrix which couples subdomains and the interface together. The remaining sections of this paper present a family of approximate techniques for both constructing and applying the Schur complement as a domain-decomposition preconditioner. The approximate Schur complement serves as an algebraic coarse space operator, thus avoiding the known difficulties associated with the direct formation of a coarse space discretization. In developing Schur complement approximations, particular attention has been given to improving sequential and parallel efficiency of implementations without significantly degrading the quality of the preconditioner. A computer code based on these developments has been tested on the IBM SP2 using MPI message passing protocol. A number of 2-D calculations are presented for both scalar advection-diffusion equations as well as the Euler equations governing compressible fluid flow to demonstrate performance of the preconditioning algorithm.

Barth, Timothy J.↗
NASA NTRS1998

Re-examination of a Classic Experiment to Measure the Positronium-Helium Cross Section

In 1975, before the advent of positronium beams, a clever experiment was carried out in an attempt to measure low-energy Ps-atom scattering cross-sections, especially that of Ps-He. A series of parallel thin plates was placed in the vessel containing the sample of gas, and positrons were allowed to enter the chamber from a radioactive source. The time spectrum of annihilation radiation was observed in the traditional way, and the pickoff annihilation due to thermalized Ps striking one of the plates was to be the special feature of the experiment. Using a diffusion analysis the authors derived cross- sections for several gases, and for helium the value obtained was sigma = 0.0166 pi(alpha)(sub o, sup 2). Even then this value was thought to be very much too small, while recent measurements and some calculations suggest a more reasonable value would be sigma = 3 to 13 pi(alpha)(sub o,sup 2). It has always been puzzling that an apparently well-designed experiment should give such an unsatisfactory result. We have now re-examined the analysis of the data with some interesting consequences which may explain the discrepancy. Two main observations support our re-analysis. First, we note that the mean free path for Ps-He scattering would be quite long if the cross-section were as small as that quoted above; the diffusion method itself would seem to be questionable. For the larger values, however, there would be no such problem. Second, it was assumed that when the annihilation rate had settled down and was following an exponential decay curve the asymptotic solution of the diffusion equation had been reached. We find, instead, that a superposition of exponentials can accurately represent the decay curve, consistent with the higher cross-section value.

Drachman, Richard J.↗
NASA NTRS2001

Ar Atmosphere: Implications for Structure and Composition of Mercury's Crust

We examine the possibilities of sustaining an argon atmosphere by diffusion from the upper 10 km of crust, and alternatively by effusion from a molten or previously molten area at great depth . Ar-40 in the atmospheres of the planets is a measure of potassium abundance in the interiors since Ar-40 is a product of radiogenic decay of K-40 by electron capture with the subsequent emission of a 1.46 eV gamma-ray. Although the Ar-40 in the earth's atmosphere is expected to have accumulated since the late bombardment, Ar-40 in surface-bounded exospheres is eroded quickly by photoionization and electron impact ionization. Thus, the argon content in the exospheres of the Moon, Mercury and probably Europa is representative of current effusion rather than accumulation over the lifetime of the body. Argon content will be a function of K content, temperature, grain size distribution, connected pore volume and possible seismic activity. Although Mercury and the Moon differ in many details, we can train the solutions to diffusion equations to predict the average lunar atmosphere. Then these parameters can be varied for Hermean conditions. Assuming a lunar crustal potassium abundance of 300 ppm, the observed argon atmosphere requires equilibrium between the argon production in the upper 9 Km of the moon (1.135 x 10(exp -3) cm(exp -3) s(exp -1)) and its loss. Hodges et al. conclude that this loss rate and the observed time variability requires argon release through seismic activity, tapping a deep argon source. An important observation is that the extreme surface of the Moon is enhanced in argon rather than depleted, as one would expect from outgassing of radiogenic argon. Manka and Michel concluded that ion implantation explains the surface enhancement of Ar-40. About half of the argon ions produced in the lunar atmosphere would return to the surface, where they would become embedded in the rocks. Similarly, at Mercury we expect the surface rocks to be enhanced in Ar-40 wherever the magnetosphere has been open over time. Thus the measurement of surface composition will reveal the long-term effects of solar wind-magnetosphere interaction. Additional information is contained in the original extended abstract.

Killen, R. M.↗
NASA NTRS1995

On the Calculation of Anisotropic Extinction Coefficients for Rigid Fibrous Ceramic Insulations

The specific anisotropic extinction coefficient e* couples the effective Radiative properties of a fibrous insulation into the radiation diffusion equation. This coefficient can be calculated using Mie scattering theory if fiber diameters, refractive indices and fiber orientation distributions are known. In general, fiber orientation distributions are not readily accessible and past calculations have considered fibers as either randomly distributed or normal to the heat flow direction. In certain rigid fibrous ceramic insulations neither of these cases apply well and a simple procedure is described for approximating e* from values calculated for the random and normal orientation cases. The intrinsic error associated with this scaling procedure is investigated. Numerical computations for several test structures and fiber materials show the average error to be less than 5% for net heat flux and radiation conductivity calculations.

Marshall, Jochen↗
NASA NTRS1995

A Hybrid Parallel Preconditioning Algorithm For CFD

A new hybrid preconditioning algorithm will be presented which combines the favorable attributes of incomplete lower-upper (ILU) factorization with the favorable attributes of the approximate inverse method recently advocated by numerous researchers. The quality of the preconditioner is adjustable and can be increased at the cost of additional computation while at the same time the storage required is roughly constant and approximately equal to the storage required for the original matrix. In addition, the preconditioning algorithm suggests an efficient and natural parallel implementation with reduced communication. Sample calculations will be presented for the numerical solution of multi-dimensional advection-diffusion equations. The matrix solver has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The full paper will show numerous examples in CFD to demonstrate the efficiency and robustness of the method.

Barth,Timothy J.↗
NASA NTRS2002

Development of a New Membrane Casting Apparatus for Studying Macrovoid Defects in Low-G

A new membrane-casting apparatus is developed for studying macrovoid defects in polymeric membranes made by the wet- and dry-casting process in low-gravity. Macrovoids are large (10-50 micron), open cavities interspersed among the smaller pores in the substructure under the gelled skin surface layer of the cast membrane. Although their occurrence is considered endemic to the wet- and dry-casting process since they can lead to compaction or skin rupture in the membrane process, recent studies suggest several useful applications such as transdermal and osmotic drug delivery systems, miniature bioreactors, etc. However, lack of knowledge about the macrovoid formation mechanism is an obstacle to further development of applications using them. An on-going debate is the role of the surface-tension-driven solutocapillary convection during macrovoid formation. The rapid growth of macrovoids within 1-5 seconds and the high polymer concentration in and near macrovoids make it difficult to explain the mechanism of macrovoid growth by diffusion alone, which is the widely accepted hypothesis proposed by Reuvers et al. The hypothesis advanced by our research group can explain this rapid growth via a mechanism that involves diffusion from the casting solution in the meta-stable region to the macrovoid enhanced by solutocapillary convection induced by the steep nonsolvent concentration gradient in the vicinity of the macrovoid. Since macrovoid growth is hypothesized to be the interplay of a solutocapillary-induced driving force counteracted by viscous drag and buoyancy, eliminate the latter provides a means for testing this hypothesis. Moreover, free convection mass transfer in the nonsolvent immersion bath used to cause phase-separation in membrane casting complicates developing a model for both the wet-casting process and macrovoid growth. The low-g environment minimizes gravitationally induced free convection thereby permitting a tractable solution to the ternary diffusion equations that characterize membrane formation. NASA's Parabolic Flight Research Aircraft provides a small window of low-g (approximately 25 s) that can be used to study macrovoid development in both wet- and dry-cast membranes if an appropriate casting apparatus is used. This casting apparatus should be able to cast the membrane in both low- and high-g in a manner so that essential one-dimensional mass transfer conditions are achieved to insure lateral uniformity in the membrane. The apparatus used in previous research on membrane casting in low-gravity was operated with the plunger driven mechanism. The spring-loaded plunger pushes the bottom block containing the polymer casting solution well directly under the absorbent chamber located in the upper stationary block. However, membranes made via this casting apparatus often displayed lateral nonuniformities that precluded obtaining quantitative information on the macrovoid growth process. Thus, it was necessary to determine the reason for these structural irregularities observed in the low-g casting apparatus. Both experimental as well as computer simulation studies of the low-g casting apparatus established that the impulsive action of the plunger caused the undesired structural nonuniformities. The simulation results showed that the width-to-depth aspect ratio of the shallow well that contains the casting solution in this apparatus was not an important factor in minimizing this problem. Even for a 40:1 (width : depth) aspect ratio, any convection induced by the horizontal motion of the interface of the casting solution will be damped out within 6.25x10(exp 4) seconds. However, the experimental studies revealed that the impulsive motion of the plunger caused a 'sloshing' of the casting solution that had to be eliminated. Therefore, the plungerdriven mechanism was changed to a cam-driven mechanism that did not cause any impulsive motion of the casting solution. Other refinements to this new membrane-casting apparatus include provision for removing the membranes from the casting wells in a less destructive manner. This was accomplished by using a slit geometry for the casting well that permitted disassembly for removal of the cast membrane. The materials used in the construction of this casting apparatus were chosen to insure wetting at the side walls and to maintain precise control of the thickness of the polymer solution in the casting well. An additional provision in this new casting apparatus is the ability to carry out both wet- as well as dry-casting. As such, this apparatus permitted the first studies of the wet-casting of polymeric membranes in low-g. Both wet- and dry-casting experiments on NASA's KC-135 research aircraft employing this new membrane-casting apparatus are scheduled in July 2002. The morphology of the resulting membranes will be characterized using an environmental scanning electron microscope (ESEM). The results of these low-g studies will be reported later.

Lee, Hanyong↗
NASA NTRS2001

Diffusion Time-Scale of Porous Pressure-Sensitive Paint

Pressure-sensitive paint (PSP) is an optical pressure sensor that utilizes the oxygen quenching of luminescence. PSP measurements in unsteady aerodynamic flows require fast time response of the paint. There are two characteristic time-scales that are related to the time response of PSP. One is the luminescent lifetime representing an intrinsic physical limit for the achievable temporal resolution of PSP. Another is the time-scale of oxygen diffusion across the PSP layer. When the time-scale of oxygen diffusion is much larger than the luminescent lifetime, the time response of PSP is controlled by oxygen diffusion. In a thin homogenous polymer layer where diffusion is Fickian, the oxygen concentration 1021 can be described by the diffusion equation in one-dimension.

Liu, Tianshu↗
NASA NTRS2003

Transient Interfacial Phenomena in Miscible Polymer Systems (TIPMPS)

Almost one hundred years ago Korteweg published a theory of how stresses could be induced in miscible fluids by concentration gradients, causing phenomena that would appear to be the same as with immiscible fluids. Miscible fluids could manifest a transient or effective interfacial tension (EIT). To this day, there has been no definitive experiment to confirm Korteweg's model but numerous fascinating and suggestive experiments have been reported. The goal of TIPMPS is to answer the question: Can concentration and temperature gradients in miscible materials induce stresses that cause convection? Many polymer processes involving miscible monomer and polymer systems could be affected by fluid flow and so this work could help understand miscible polymer processing, not only in microgravity, but also on earth. Demonstrating the existence of this phenomenon in miscible fluids will open up a new area of study for materials science. The science objectives of TIPMPS are: (1) Determine if convection can be induced by variation of the width of a miscible interface; (2) Determine if convection can be induced by variation of temperature along a miscible interface; (3) Determine if convection can be induced by variation of conversion along a miscible interface An interface between two miscible fluids can best be created via a spatially-selective photopolymerization of dodecyl acrylate with a photoinitiator, which allows the creation of precise and accurate concentration gradients between polymer and monomer. Optical techniques will be used to measure the refractive index variation caused by the resultant temperature and concentration fields. The viscosity of the polymer will be measured from the increase in the fluorescence of pyrene. Because the large concentration and temperature gradients cause buoyancy-driven convection that prevents the observation of the predicted flows, the experiment must be done in microgravity. In this report, we will consider our efforts to estimate the square gradient parameter, k, and our use of the estimates in modeling of the planned TIPMPS experiments. We developed a model consisting of the heat and diffusion equations with convective terms and of the Navier-Stokes equations with an additional volume force written in the form of the Korteweg stresses arising from nonlocal interaction in the fluid. The fluid's viscosity dependence on polymer conversion and temperature was taken from measurements of poly(dodecyl acrylate). Numerical modeling demonstrated that significant flows would arise for conditions corresponding to the planned experiments.

Pojman, John A.↗
NASA NTRS2000

Thermographic Imaging of Defects in Anisotropic Composites

Composite materials are of increasing interest to the aerospace industry as a result of their weight versus performance characteristics. One of the disadvantages of composites is the high cost of fabrication and post inspection with conventional ultrasonic scanning systems. The high cost of inspection is driven by the need for scanning systems which can follow large curve surfaces. Additionally, either large water tanks or water squirters are required to couple the ultrasonics into the part. Thermographic techniques offer significant advantages over conventional ultrasonics by not requiring physical coupling between the part and sensor. The thermographic system can easily inspect large curved surface without requiring a surface following scanner. However, implementation of Thermal Nondestructive Evaluations (TNDE) for flaw detection in composite materials and structures requires determining its limit. Advanced algorithms have been developed to enable locating and sizing defects in carbon fiber reinforced plastic (CFRP). Thermal Tomography is a very promising method for visualizing the size and location of defects in materials such as CFRP. However, further investigations are required to determine its capabilities for inspection of thick composites. In present work we have studied influence of the anisotropy on the reconstructed image of a defect generated by an inversion technique. The composite material is considered as homogeneous with macro properties: thermal conductivity K, specific heat c, and density rho. The simulation process involves two sequential steps: solving the three dimensional transient heat diffusion equation for a sample with a defect, then estimating the defect location and size from the surface spatial and temporal thermal distributions (inverse problem), calculated from the simulations.

Plotnikov, Y. A.↗
NASA NTRS2002

Basal melt beneath whillans ice stream and ice streams A and C

We have used a recently derived map of the velocity of Whillans Ice Stream and Ice Streams A and C to help estimate basal melt. Temperature was modeled with a simple vertical advection-diffusion equation, 'tuned' to match temperature profiles. We find that most of the melt occurs beneath the tributaries where larger basal shear stresses and thicker ice favors greater melt (e.g., 10-20 mm/yr). The occurrence of basal freezing is predicted beneath much of the ice plains of Ice Stream C and Whillans Ice Stream. Modelled melt rates for when Ice Stream C was active suggest there was just enough melt water generated in its tributaries to balance basal freezing on its ice plain. Net basal melt for Whillans Ice Stream is positive due to smaller basal temperature gradients. Modelled temperatures on Whillans Ice Stream, however, were constrained by a single temperature profile at UpB. Basal temperature gradients for Whillans B1 and Ice Stream A may have conditions more similar to those beneath Ice Streams C and D, in which case, there may not be sufficient melt to sustain motion. This would be consistent with the steady deceleration of Whillans stream over the last few decades.

glaciers ice sheets↗
NASA NTRS2006

A New Empirical Thermospheric Density Model JB2006 using New Solar Indices

A new empirical atmospheric density model is developed using the CIRA72 (Jacchia 71) model as the basis for the diffusion equations. New solar indices based on orbit based sensor data are used for the solar irradiances in the extreme and far ultraviolet wavelengths. New exospheric temperature and semiannual density equations are employed to represent the major thermospheric density variations. Temperature correction equations are also developed for diurnal and latitudinal effects, and finally density correction factors are used for model corrections required at high altitude (1500-4000 km). The new model, Jacchia-Bowman 2006, is validated through comparisons of accurate daily density drag data previously computed for numerous satellite. For 400 km altitude the standard deviation of 16% for the standard Jacchia model is reduced to 10% for the new JB2006 model for periods of low geomagnetic storm activity.

Bowman, Bruce R.↗
NASA NTRS2007

Short-Term Forecasting of Radiation Belt and Ring Current

A computer program implements a mathematical model of the radiation-belt and ring-current plasmas resulting from interactions between the solar wind and the Earth s magnetic field, for the purpose of predicting fluxes of energetic electrons (10 keV to 5 MeV) and protons (10 keV to 1 MeV), which are hazardous to humans and spacecraft. Given solar-wind and interplanetary-magnetic-field data as inputs, the program solves the convection-diffusion equations of plasma distribution functions in the range of 2 to 10 Earth radii. Phenomena represented in the model include particle drifts resulting from the gradient and curvature of the magnetic field; electric fields associated with the rotation of the Earth, convection, and temporal variation of the magnetic field; and losses along particle-drift paths. The model can readily accommodate new magnetic- and electric-field submodels and new information regarding physical processes that drive the radiation-belt and ring-current plasmas. Despite the complexity of the model, the program can be run in real time on ordinary computers. At present, the program can calculate present electron and proton fluxes; after further development, it should be able to predict the fluxes 24 hours in advance

Fok, Mei-Ching↗