Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Rotor Vortex Filaments: Living on the Slipstream's Edge

The purpose of this paper is to gain a better understanding of rotor wake evolution in hover and axial flow by deriving an analytical solution for the time dependent behavior of vortex filament circulation and core size. This solution is applicable only for vortex filaments in the rotor far-wake. A primarily inviscid vortex/shear layer interaction (where the slipstream boundary is modeled as a shear layer) has been identified in this analytical treatment. This vortex/shear layer interaction results in decreasing, vortex filament circulation and core size with time. The inviscid vortex/shear layer interaction is shown, in a first-order treatment, to be of greater magnitude than viscous diffusion effects. The rate of contraction, and ultimate collapse, of the vortex filament core is found to be directly proportional to the rotor inflow velocity. This new insight into vortex filament decay promises to help reconcile several disparate observations made in the literature and will, hopefully, promote new advances in theoretical modeling of rotor wakes.

Young, Larry A.↗

Heterogeneous Combustion of Porous Graphite Particles in Microgravity

Recent theoretical investigations on graphite particle combustion have employed several levels of heterogeneous reaction models, ranging from global to elementary models, to describe the oxidation of carbon to gaseous products. Unlike the counterpart homogeneous reaction models, these heterogeneous reaction models are not well developed because of the difficulties associated with decoupling the physical characteristics of the solid (e.g. surface area taking part in combustion) from the chemical kinetic data. This is certainly true for porous graphite particle combustion, where heterogeneous and homogeneous reactions occur within the pores and play an important role in the overall oxidation process. As a result, there are considerable uncertainties of physical phenomena predicted using different heterogeneous kinetic models available in the literature. A good example, discussed later in this paper, is the predicted critical particle size below which the mass burning rate becomes exponentially small. The main goal of this study is to understand the basic mechanism controlling such rapid changes in burning rates, by developing a model where physical contributions are decoupled from chemical rate constants in a consistent manner. Another important goal of the proposed study is to develop a truly intrinsic, detailed heterogeneous reaction model for porous graphite combustion at high-temperatures, and to derive a systematically reduced heterogeneous reaction model in terms of the elementary reaction rate constants of the detailed model. The validation of chemical kinetic models describing the heterogeneous and homogeneous combustion in and around a spherically symmetric porous graphite particle can be considerably simplified by experimental measurements obtained under microgravity conditions. A vital component of this study is to conduct such supporting experiments on particle burning rate and surface temperature using NASA microgravity facilities, in close coordination with the theoretical effort. The basic understanding obtained and models developed as part of this project will be useful for optimal design of coal combustion devices. These models can also be extended to investigate the role of heterogeneous chemistry on pollutant formation pathways in combustion devices. The theoretical approach developed here, with pore diffusion effects decoupled from the chemical effects, can also be extended to understand the heterogeneous combustion of other porous fuels, for example, combustion of magnesium in a CO2 environment for propulsion in the Martian atmosphere.

Chelliah, Harsha K.↗

Tensile Behavior of Al2o3/feal + B and Al2o3/fecraly Composites

The feasibility of Al2O3/FeAl + B and Al2O3/FeCrAlY composites for high-temperature applications was assessed. The major emphasis was on tensile behavior of both the monolithics and composites from 298 to 1100 K. However, the study also included determining the chemical compatibility of the composites, measuring the interfacial shear strengths, and investigating the effect of processing on the strength of the single-crystal Al2O3 fibers. The interfacial shear strengths were low for Al203/FeAl + B and moderate to high for Al203/FeCrAlY. The difference in interfacial bond strengths between the two systems affected the tensile behavior of the composites. The strength of the Al203 fiber was significantly degraded after composite processing for both composite systems and resulted in poor composite tensile properties. The ultimate tensile strength (UTS) values of the composites could generally be predicted with either rule of mixtures (ROM) calculations or existing models when using the strength of the etched-out fiber. The Al2O3/FeAl + B composite system was determined to be unfeasible due to poor interfacial shear strengths and a large mismatch in coefficient of thermal expansion (CTE). Development of the Al2O3/FeCrAlY system would require an effective diffusion barrier to minimize the fiber strength degradation during processing and elevated temperature service.

Draper, S. L.↗

Quantitative Examination of Corrosion Damage by Means of Thermal Response Measurements

Two computational methods are presented that enable a characterization of corrosion damage to be performed from thermal response measurements derived from a standard flash thermographic inspection. The first is based upon a one dimensional analytical solution to the heat diffusion equation and presumes the lateral extent of damage is large compared to the residual structural thickness, such that lateral heat diffusion effects can be considered insignificant. The second proposed method, based on a finite element optimization scheme, addresses the more general case where these conditions are not met. Results from an experimental application are given to illustrate the precision, robustness and practical efficacy of both methods.

Rajic, Nik↗

Ion Layer Separation and Equilibrium Zonal Winds in Midlatitude Sporadic E

In-situ observations of a moderately strong mid-latitude sporadic-E layer show a separation in altitude between distinct sublayers composed of Fe(+), Mg(+), and NO(+). From these observations it is possible to estimate the zonal wind field consistent with diffusive equilibrium near the altitude of the layer. The amplitude of the zonal wind necessary to sustain the layer against diffusive effects is less than 10 meters per second, and the vertical wavelength is less than 10 km.

Earle, G. D.↗

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗

Turbulent Jet Flames Into a Vitiated Coflow

Examined is the vitiated coflow flame, an experimental condition that decouples the combustion processes of flows found in practical combustors from the associated recirculating fluid mechanics. The configuration consists of a 4.57 mm diameter fuel jet into a coaxial flow of hot combustion products from a lean premixed flame. The 210 mm diameter coflow isolates the jet flame from the cool ambient, providing a hot environment similar to the operating conditions of advanced combustors; this important high temperature element is lacking in the traditional laboratory experiments of jet flames into cool (room) air. A family of flows of increasing complexity is presented: 1) nonreacting flow, 2) all hydrogen flame (fuel jet and premixed coflow), and 3) set of methane flames. This sequence of experiments provides a convenient ordering of validation data for combustion models. Laser Raman-Rayleigh-LIF diagnostics at the Turbulent Diffusion Flame laboratory of Sandia National Laboratories produced instantaneous multiscalar point measurements. These results attest to the attractive features of the vitiated coflow burner and the well-defined boundary conditions provided by the coflow. The coflow is uniform and steady, isolating the jet flame from the laboratory air for a downstream distance ranging from z/d = 50-70. The statistical results show that differential diffusion effects in this highly turbulent flow are negligible. Complementing the comprehensive set of multiscalar measurements is a parametric study of lifted methane flames that was conducted to analyze flame sensitivity to jet and coflow velocity, as well as coflow temperature. The linear relationship found between the lift-off height and the jet velocity is consistent with previous experiments. New linear sensitivities were found correlating the lift-off height to coflow velocity and temperature. A blow-off study revealed that the methane flame blows off at a common coflow temperature (1260 K), regardless of coflow or jet velocity. An explanation for this phenomenon is that entrainment of ambient air at the high lift-off heights prevents autoignition. Analysis of the results suggests that flame stabilization occurs through a combination of flame propagation, autoignition, and localized extinction processes. Proposed is an expanded view of distributed reaction combustion based on analysis of the distributions of probe volume conditions at the stabilization region of the lifted hydrogen and methane flames. Turbulent eddies the size of the flame thickness mix fuel and hot coflow across the flame front, thereby enhancing the reaction zone with autoignition of reactants at elevated temperatures; this is the reverse effect of turbulent flames in ambient air, where intense turbulence in cool mixtures result in localized extinction. Each of the three processes (i.e., flame propagation, autoignition and localized extinction) contributes to flame stabilization in varying degrees, depending on flow conditions.

Holdeman, James D.↗

Neuroprotection trek--the next generation: neuromodulation I. Techniques--deep brain stimulation, vagus nerve stimulation, and transcranial magnetic stimulation

Neuromodulation denotes controlled electrical stimulation of the central or peripheral nervous system. The three forms of neuromodulation described in this paper-deep brain stimulation, vagus nerve stimulation, and transcranial magnetic stimulation-were chosen primarily for their demonstrated or potential clinical usefulness. Deep brain stimulation is a completely implanted technique for improving movement disorders, such as Parkinson's disease, by very focal electrical stimulation of the brain-a technique that employs well-established hardware (electrode and pulse generator/battery). Vagus nerve stimulation is similar to deep brain stimulation in being well-established (for the treatment of refractory epilepsy), completely implanted, and having hardware that can be considered standard at the present time. Vagus nerve stimulation differs from deep brain stimulation, however, in that afferent stimulation of the vagus nerve results in diffuse effects on many regions throughout the brain. Although use of deep brain stimulation for applications beyond movement disorders will no doubt involve placing the stimulating electrode(s) in regions other than the thalamus, subthalamus, or globus pallidus, the use of vagus nerve stimulation for applications beyond epilepsy-for example, depression and eating disorders-is unlikely to require altering the hardware significantly (although stimulation protocols may differ). Transcranial magnetic stimulation is an example of an external or non-implanted, intermittent (at least given the current state of the hardware) stimulation technique, the clinical value of which for neuromodulation and neuroprotection remains to be determined.

Review, Tutorial↗

Correspondence between discrete and continuous models of excitable media: trigger waves

We present a theoretical framework for relating continuous partial differential equation (PDE) models of excitable media to discrete cellular automata (CA) models on a randomized lattice. These relations establish a quantitative link between the CA model and the specific physical system under study. We derive expressions for the CA model's plane wave speed, critical curvature, and effective diffusion constant in terms of the model's internal parameters (the interaction radius, excitation threshold, and time step). We then equate these expressions to the corresponding quantities obtained from solution of the PDEs (for a fixed excitability). This yields a set of coupled equations with a unique solution for the required CA parameter values. Here we restrict our analysis to "trigger" wave solutions obtained in the limiting case of a two-dimensional excitable medium with no recovery processes. We tested the correspondence between our CA model and two PDE models (the FitzHugh-Nagumo medium and a medium with a "sawtooth" nonlinear reaction source) and found good agreement with the numerical solutions of the PDEs. Our results suggest that the behavior of trigger waves is actually controlled by a small number of parameters.

Non-NASA Center↗

Linear Temporal Stability Analysis of a Low-Density Round Gas Jet Injected into a High-Density Gas

It has been observed in previous experimental studies that round helium jets injected into air display a repetitive structure for a long distance, somewhat similar to the buoyancy-induced flickering observed in diffusion flames. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis of a round helium jet injected into air was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. The temporal growth rates of the disturbances increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity.

Lawson, Anthony L.↗

Towards a 3D Space Radiation Transport Code

High-speed computational procedures for space radiation shielding have relied on asymptotic expansions in terms of the off-axis scatter and replacement of the general geometry problem by a collection of flat plates. This type of solution was derived for application to human rated systems in which the radius of the shielded volume is large compared to the off-axis diffusion limiting leakage at lateral boundaries. Over the decades these computational codes are relatively complete and lateral diffusion effects are now being added. The analysis for developing a practical full 3D space shielding code is presented.

Wilson, J. W.↗

Large-scale Atmospheric Transport Processes

Continuing earlier work, we continued an investigation of the seasonal behavior of the edges of the stratospheric surf zone. These edges form a barrier between the rapidly mixed surf zone and the relatively isolated tropics. In collaboration with Dr Lynn Sparling at GSFC, we used a statistical analysis of HALOE and CLAES trace gas data from UARS to identify and locate these edges during each UARS observing period. We found that the edges on both sides of the equator are present all year (a fact that is important for conceptual models of stratospheric transport), though that on the summer side of the equator is much less sharp than the winter edge. The edges migrate seasonally into the summer hemisphere. Their location also shows influence of the QBO, together with the SAO at higher altitudes. Comparisons with effective diffusivities, and the edge locations, suggest that the edge is sustained by surf zone entrainment during winter, but by the residual circulation during summer.

Plumb, R. Alan↗

Apparent I-Xe Cooling Rates of Chondrules Compared with Silicates from the Colomera Iron Meteorite

In I-Xe dating, a regular pattern of increasing Xe-129/128Xe-128 ratio with increasing extraction temperature is often observed. If one makes the crude assumption that the temperatures at which the Xe is extracted in the laboratory is approximately the same as the temperature at which those sites closed 4.6 Ga ago, a (zeroth order) model cooling rate can be found. In order to test and refine this model we can apply the cooling theory of Dodson to those extraction steps approaching the I-Xe isochrons. Using an Arrhenius plot for these temperature fractions, and assuming that an only single phase is involved, the effective diffusion parameters can be estimated (frequency factor and activation energy). From the apparent (zeroth order) cooling rate, the closure temperature can be estimated from the Dodson equation. This model closure temperature can then be compared with the actual laboratory temperature at which the isochron begins. The ratio of the closure temperature and the temperature corresponding to the start of the isochron provides the ratio of the two temperature scales, incorporation and extraction. The actual cooling rate is then given by the apparent (zeroth order) cooling rate times the temperature scale factor. Figure 1 shows Arrhenius plots for I-Xe data

Hohenberg, C. M.↗

Mechanically Induced g-Jitter from Space Station Rotary Joints

The mission of the International Space Station is to provide a working laboratory in orbit for research in engineering, life sciences, and microgravity. Among the microgravity disciplines that are preparing to utilize this international resource are materials processing, combustion, fluid dynamics, biotechnology, and fundamental physics. The Station promises to enable significant advances in each of these areas by making available a research facility in which gravitational and other accelerations, and their corresponding buoyancy and diffusion effects on various physical processes, are orders of magnitude lower than they are on Earth. In order to fulfill this promise, it is not enough for the Space Station to simply replicate a typical terrestrial scientific laboratory in orbit. Although an orbiting laboratory is free of most of the effects of gravitational acceleration by virtue of its free fall condition, it also produces structural vibration or jitter that can interfere with the processes under study. To ensure the quality of the acceleration environment and enable a successful mission, the Space Station Program has limited potential disturbances in two ways: first, by isolating the most sensitive payloads from the vehicle structure, and second, by quieting major disturbances at their sources. The first area, payload isolation, is implemented inside the pressurized modules at the rack level. Sub-rack level isolators have also been developed. This paper addresses the second area, disturbance source limits, for one of the major sources of mechanical noise on the Space Station: the Solar Alpha Rotary Joints. Due to the potential for large disturbances to the microgravity environment, an initial analytical prediction of rotary joint vibration output was made. Key components were identified and tested to validate the analytical predictions. Based on the component test results, the final vibration output of the joints was verified by a test on each fully assembled flight unit. This paper describes the Space Station microgravity requirements, the rotary joint hardware, and the disturbance producing aspects of joint operation. The test setup, instrumentation, test conditions, and results for the component level and system level measurements are described. An overall forcing function that describes the maximum torque imparted to the Station is created based on the test results, and these disturbances are shown to meet the applicable Space Station microgravity requirements.

Boucher, Robert L.↗

Radiation Belt Environment Model: Application to Space Weather and Beyond

Understanding the dynamics and variability of the radiation belts are of great scientific and space weather significance. A physics-based Radiation Belt Environment (RBE) model has been developed to simulate and predict the radiation particle intensities. The RBE model considers the influences from the solar wind, ring current and plasmasphere. It takes into account the particle drift in realistic, time-varying magnetic and electric field, and includes diffusive effects of wave-particle interactions with various wave modes in the magnetosphere. The RBE model has been used to perform event studies and real-time prediction of energetic electron fluxes. In this talk, we will describe the RBE model equation, inputs and capabilities. Recent advancement in space weather application and artificial radiation belt study will be discussed as well.

Fok, Mei-Ching H.↗

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

Toward 2D Dynamo Models Calibrated by Global 3D Relativistic Accretion Disk Simulations

Two-dimensional models assuming axisymmetry are an economical way to explore the long-term evolution of black hole accretion disks, but they are only realistic if the feedback of the nonaxisymmetric turbulence on the mean momentum and magnetic fields is incorporated. Dynamo terms added to the 2D induction equation should be calibrated to 3D MHD simulations. For generality, the dynamo tensors should be calibrated as functions of local variables rather than explicit functions of spatial coordinates in a particular basis. In this paper, we study the feedback of non-axisymmetric features on the 2D mean fields using a global 3D, relativistic, Cartesian simulation from the Illinois GRMHD code. We introduce new methods for estimating overall dynamo alpha and turbulent diffusivity effects as well as measures of the dominance of non-axisymmetric components of energies and fluxes within the disk interior. We attempt closure models of the dynamo EMF using least squares fitting, considering both models where coefficient tensors are functions of space and more global, covariant models. None of these models are judged satisfactory, but we are able to draw conclusions on what sorts of generalizations are and are not promising.

Matthew D Duez↗

End Effects and Load Diffusion in Composite Structures

The research carried out here builds on our previous NASA supported research on the general topic of edge effects and load diffusion in composite structures. Further fundamental solid mechanics studies were carried out to provide a basis for assessing the complicated modeling necessary for large scale structures used by NASA. An understanding of the fundamental mechanisms of load diffusion in composite subcomponents is essential in developing primary composite structures. Specific problems recently considered were focussed on end effects in sandwich structures and for functionally graded materials. Both linear and nonlinear (geometric and material) problems have been addressed. Our goal is the development of readily applicable design formulas for the decay lengths in terms of non-dimensional material and geometric parameters. Analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and assessing results from finite element analyses. The decay behavior of stresses and other field quantities provides a significant aid towards this process. The analysis is also amenable to parameter study with a large parameter space and should be useful in structural tailoring studies.

Horgan, Cornelius O.↗