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 1,675 records · Page 93

An adaptive grid with directional control

An adaptive grid generator for adaptive node movement is here derived by combining a variational formulation of Winslow's (1981) variable-diffusion method with a directional control functional. By applying harmonic-function theory, it becomes possible to define conditions under which there exist unique solutions of the resulting elliptic equations. The results obtained for the grid generator's application to the complex problem posed by the fluid instability-driven magnetic field reconnection demonstrate one-tenth the computational cost of either a Eulerian grid or an adaptive grid without directional control.

Brackbill, J. U.↗

Unusual satellite-electron signature within the Uranian magnetosphere and its implications regarding whistler electron loss processes

It has been reported that during the outbound (nightside) portion of the Voyager 2 encounter with the Uranian magnetosphere, intense whistler mode emissions were observed near the magnetic equator (lambda(sub m) approx. 16 deg) and at L shell values between approx. 5.5 and approx. 9 R(sub U). Comprehensive calculations of whistler-driven pitch angle diffusion, in previous work, have yielded strong diffusion electron lifetimes of approx. 1 hour for 20 to 40 keV electrons. In this paper we report on an unusual and sharply defined charged particle feature that: (1) involved electrons between 22 and 35 keV, (2) was observed during the time period of the intense whistler mode observations, (3) was aligned very accurately and sharply with the minimum L shell position (L approx. 7.5) of the satellite Ariel, and (4) has an appearance that suggests that electrons were removed only at and beyond Ariel's minimum-L. On the basis of our conclusion that the signature was caused by electron interactions with either Ariel or materials distributed along Ariel's orbit, the signature could not have been generated for at least 10 hours prior to its observation. Thus the calculated whistler loss times are in apparent conflict with the signature observation. A scenario of events is proposed to explain the data that involves substormlike electron acceleration on the Uranian nightside and a subsequent sculpting of the electron spatial distributions via interactions with Ariel or materials distributed along Ariel's orbit. The possibility exists that the accurate alignment of the sharp electron feature with Ariel's minimum-L, and the absorptionlike character of the feature, are accidental, and that the feature is caused by dynamical processes (e.g., substorms). In this case the dynamical processes must be quite dissimilar to those occurring in the Earth's magnetosphere.

Mauk, B. H.↗

The First National Student Conference: NASA University Research Centers at Minority Institutions

The conference includes contributions from 13 minority universities with NASA University Research Centers. Topics discussed include: leadership, survival strategies, life support systems, food systems, simulated hypergravity, chromium diffusion doping, radiation effects on dc-dc converters, metal oxide glasses, crystal growth of Bil3, science and communication on wheels, semiconductor thin films, numerical solution of random algebraic equations, fuzzy logic control, spatial resolution of satellite images, programming language development, nitric oxide in the thermosphere and mesosphere, high performance polyimides, crossover control in genetic algorithms, hyperthermal ion scattering, etc.

Daso, Endwell O.↗

Overview of Terra and Aqua MODIS Status

This presentation will consist of two one-hour lectures. The first lecture will review the characteristics of the MODIS instrument and how it responds to the performance specifications. The pre-launch and on-orbit procedures and results to characterize and maintain instrument output so as to allow the prescribed science to be done will also be outlined. This will include defining the equations used to account for gain changes and any degradation in the optics of the instrument or the on-board calibration devices themselves. The changes in the solar diffuser and the radiance versus scan angle performance of the scan mirror will also be reviewed. Overall it will be shown that the instrument has and is meeting specifications of 2% reflectance relative to the sun, 5% in the radiance observed in the reflected solar radiation bands, and 0.5-1% in the thermal bands. The second lecture will review and describe highlights in the geophysical products and related science results. There are approximately 40 geophysical products related to observations of land, ocean and atmosphere features. Many of the results are unprecedented and offer considerable advances over those achievable with heritage instruments such as the NOAA Advanced Very High Resolution Radiometer (AVHRR). The literature is showing a steady growth of publications in scientific journals using MODIS data or products. The future is also bright in that a follow-on instrument based on the MODIS will be flown on the National Polar-Orbiting Environmental Satellite Series (NPOESS) starting around 2010.

Salomonson, Vincent V.↗

Magnetic anisotropy and two-dimensional short-range chemical ordering in Ba 1−𝑥 ⁢Na 𝑥 ⁢Fe 2 ⁢As 2

A true understanding of the properties of pnictide superconductors requires the development of high-quality materials and performing measurements designed to unravel their intrinsic properties and short-range nematic correlations which are often obscured by extrinsic effects such as poor crystallinity, inhomogeneity, domain formation, and twinning. Here, in this paper, we report the systematic growth of high-quality Na-substituted BaFe 2 As 2 single crystals and their characterization using pulsed magnetic fields x-ray diffraction and x-ray diffuse scattering. Analysis of the properties and compositions of the highest-quality crystals shows that their actual Na stoichiometry is about 50-60% of the nominal content and that the targeted production of crystals with specific compositions is accessible. We derived a reliable equation to estimate the Na stoichiometry based on the measured superconducting T c of these materials. Attempting to force spin reorientation and induce tetragonality, orthorhombic Ba 1-x Na x Fe 2 As 2 single crystals subjected to out-of-plane magnetic fields up to 31.4T are found to exhibit strong in-plane magnetic anisotropy demonstrated by the insufficiency of such high fields in manipulating the relative population of their twinned domains or in suppressing the orthorhombic order. Broad x-ray diffuse-intensity rods observed at temperatures between 30 and 300 K uncover short-range structural correlations. Local structure modeling together with 3D-Δ pair-distribution function mapping of real-space interatomic vectors show that the diffuse scattering arises from in-plane short-range chemical correlations of the Ba and Na atoms coupled with short-range atomic displacements within the same plane due to an effective size difference between the two atomic species.

36 MATERIALS SCIENCE↗

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Absorption and Emission Characteristics of Diffuse Spherical Enclosures

An analysis was made to determine the energy absorbed when radiation from an external source enters a spherical cavity with diffusely reflecting walls. It was found that both the overall energy absorbed in the cavity and the local distribution of absorbed energy could be expressed in terms of simple algebraic equations that are valid for any arbitrary spatial and directional distribution of the incoming radiation. In addition, the characteristics of an isothermal spherical cavity as a possible source of near black-body radiation were investigated. This information was also expressed by simple algebraic relations.

ABSORPTION↗

Possible theoretical explanations for occasional days of non-field-aligned diffusion at neutron monitor energies

It has been shown previously (Anath et al., 1973 and Kane, 1974) that 20 to 25% of days, the diffusion component of the cosmic-ray neutron diurnal anisotropy is directed more than 30 degrees away from the ecliptic projection of the interplanetary magnetic field averaged over the same 24 hours. A number of explanations for this deviation are discussed and it is concluded that transverse gradient drifts due to gradients perpendicular to the ecliptic are likely, that diurnal variations in the diffusion component of the neutron anisotropy may affect results from single stations and that the 24 hour mean interplanetary magnetic field may not be the field appropriate to the streaming equation at neutron monitor energies.

Forman, M. A.↗

Repeated-cascade theory of strong turbulence in a magnetized plasma

A two-dimensional Navier-Stokes equation of vorticity in fluid turbulence is used to model drift turbulence in a plasma with a strong constant magnetic field and a constant mean density gradient. The nonlinear eddy diffusivity is described by a time-integrated Lagrangian correlation of velocities, and the repeated-cascade method is employed to choose the rank accounting for nearest-neighbor interactions, to calculate the Lagrangian correlation, and to close the correlation hierarchy. As a result, the diffusivity becomes dependent on the plasma's induced diffusion and is represented by a memory chain that is cut off by similarity and inertial randomization. Spectral laws relating the kinetic-energy spectrum to the -5, -5/2, -3, and -11 powers of wavenumber are derived for the velocity subranges of production, approach to inertia, inertia, and dissipation, respectively. It is found that the diffusivity is proportional to some inverse power of the magnetic field, that power being 1, 2/3, 5/6, and 2, respectively, for the four velocity subranges.

Tchen, C. M.↗

Numerical method for boundary layers with blowing - The exponential box scheme

The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.

El-Mistikawy, T. M.↗

The hydrogenation of metals upon interaction with water

Hydrogen evolution at 600 deg and 5 x 10 to the 7th power - 10 to the 6th power torr from AVOOO Al samples, which were pickled in 10 percent NaOH, is discussed. An H evolution kinetic equation is derived for samples of equal vol. and different surfaces (5 and 20 sq cm). The values of the H evolution coefficient K indicated an agreement with considered H diffusion mechanism through an oxide layer. The activation energy for the H evolution process, obtained from the K-temp. relation, was 13,000 2000 cal/g-atom.

Andreyev, L. A.↗

The equilibrium structure of thin magnetic flux tubes. I

A model atmosphere within a thin magnetic flux tube that is embedded in an arbitrarily stratified medium is presently constructed by solving the radiative transfer equation in the two-stream approximation for gray opacity, under the assumption that the magnetic field is sufficiently strong to warrant the neglect of both thermal conduction and convective diffusion; energy inside the flux tube therefore being transported solely by radiation. The structure of the internal atmosphere is determined on the basis of the hydrostatic and radiative equilibrium conditions of the tube embedded in an external atmosphere. The gas temperature along the axis of the tube is computed, and the geometry of the flux tube is determined on the basis of magnetic flux conservation and total pressure equilibrium.

Ferrari, A.↗

Vapor concentration measurement with photothermal deflectometry

Theoretical and experimental results for using the photothermal deflection technique to measure vapor species concentration, while minimizing the disturbance of the transport (material) parameters due to vapor heating, are developed and described. In contrast to common practice, the above constraints require using a pump-beam duty cycle of less than 50 percent. The theoretical description of the shortened heating time is based on a step-function formulation of the pumping cycle. The results are obtained as closed-form solutions of the energy equation for many chopping cycles until steady state is reached, by use of a Green's-function method. The Euler formulation of the Fermat principle is used to calculate the deflection angle. The equations are expanded to include the effects of vapor velocity on both the temperature and temperature gradient profiles. The effects of finite (unfocused) pump and probe beams and thermal (Soret) diffusion are also accounted for. Excellent agreement between theory and experiment is obtained.

Banish, R. Michael↗

Continuum theories for fluid-particle flows: Some aspects of lift forces and turbulence

A general framework is outlined for the modeling of fluid particle flows. The momentum exchange between the constituents embodies both lift and drag forces, constitutive equations for which can be made explicit with reference to known single particle analysis. Relevant results for lift are reviewed, and invariant representations are posed. The fluid and particle velocities and the particle volume fraction are then decomposed into mean and fluctuating parts to characterize turbulent motions, and the equations of motion are averaged. In addition to the Reynolds stresses, further correlations between concentration and velocity fluctuations appear. These can be identified with turbulent transport processes such as eddy diffusion of the particles. When the drag force is dominant, the classical convection dispersion model for turbulent transport of particles is recovered. When other interaction forces enter, particle segregation effects can arise. This is illustrated qualitatively by consideration of turbulent channel flow with lift effects included.

Mctigue, David F.↗

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic.

Turkel, E.↗

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time-dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that remove the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock-capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional viscous flow over a blunt biconic.

Turkel, E.↗

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗