Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface diffusion”

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 181 records · Page 10

Ignition, Burning, and Extinction of a Strained Fuel Strip

Flame structure and ignition and extinction processes associated with a strained fuel strip are explored numerically using detailed transport and complex kinetics for a propane-air reaction. Ignition modes are identified that are similar to those predicted by one-step activation energy asymptotics, i.e., modes in which diffusion flames can ignite as independent or dependent interfaces, and modes in which single premixed or partially premixed flames ignite. These ignition modes are found to be dependent on critical combinations of strain rate, fuel strip thickness, and initial reactant temperatures. Extinction in this configuration is seen to occur due to fuel consumption by adjacent flames, although viscosity is seen to have the effect of delaying extinction by reducing the effective strain rate and velocity field experienced by the flames.

Selerland, T.↗

Ignition, Burning and Extinction of a Strained Fuel Strip with Complex Kinetics

Flame structure and ignition and extinction processes associated with a strained fuel strip are explored numerically using detailed transport and complex kinetics for a propane-air reaction. Ignition modes are identified that are similar to those predicted by one-step activation energy asymptotics, i.e., modes in which diffusion flames can ignite as independent or dependent interfaces and modes in which single premixed or partially premixed flames ignite and burn. These ignition modes have been found to be dependent on critical combinations of strain rate, fuel strip thickness and initial reactant temperatures. The formation of NO/NO2 is found to be strongly dependent on strain rate and the local molecular mixing of reactants which occurs as a consequence of strain. Extinction in this configuration is seen to occur due to fuel consumption by adjacent flames, although viscosity is seen to have the effect of delaying extinction by reducing the effective strain rate experienced by the flames. Response of the flames to oscillatory strain rates is seen to be strongly dependent on the amplitude and frequency of the oscillation.

Selerland, T.↗

Influences of Models on the Unsteady Pressure Characteristics of the NASA National Transonic Facility

Pressure fluctuations have been measured over the course of several tests in the National Transonic Facility to study unsteady phenomenon both with and without the influence of a model. Broadband spectral analysis will be used to characterize the length scales of the tunnel. Special attention will be given to the large-scale, low frequency data that influences the Mach number and force and moment variability. This paper will also discuss the significance of the vorticity and sound fields that can be related to the Common Research Model and will also highlight the comparisons to an empty tunnel configuration. The effectiveness of vortex generators placed at the interface of the test section and wind tunnel diffuser showed promise in reducing the empty tunnel unsteadiness, however, the vortex generators were ineffective in the presence of a model.

Jones, Gregory↗

Deep Chandra Observation and Numerical Studies of the Nearest Cluster Cold Front in the Sky

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are approximately 10 per cent brighter than the surrounding gas and are separated by approximately 15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches approximately 5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

hydrodynamics↗

Diffusivities of Redox-Sensitive Elements in Basalt vs. Oxygen Fugacity Determined by LA-ICP-MS

Several diffusion experiments were conducted in a piston cylinder device across a range of oxygen fugacities (FMQ-3 FMQ-1.2, FMQ+6) at 1 GPa and 1300 C. This was done to explore the effects of oxygen fugacity (fO2) on diffusivity of redox sensitive trace elements. This allows investigation of how these elements diffuse across the fO2 range encountered in different reservoirs on planets and moons in our solar system. The University of Rochester LA-ICP-MS system was used for analysis of samples. Analyses were conducted using an Agilent 7900 quadrupole mass spectrometer connected to a Photon Machines 193 nm G2 laser ablation (LA) system equipped with a HelEx 2-volume sample chamber. Spots used were 35 micrometers circles spaced at 65 micrometers intervals. Laser fluence was 7.81 J/cm^2 with a rep rate of 10 Hz. The iolite software package was used to reduce data collected from laser ablation analysis of experiments with Si-29 used as the internal standard isotope. Iolite's global fit module was used to simultaneously fit elements' diffusivities in each experiment while keeping the Matano interface constant. Elements analysed include V, Nb, W, Mo, La, Ce, Pr, Sm, Eu, Gd, Ta, and W. Figures

Szumila, Ian↗

J-OCTA Applications to in Space Manufacturing and Nanofluidic Technologies

We consider examples of J-OCTA software application to physics-based analysis of in-space manufacturing technologies. Firstly, we consider an atomistic model of thermal welding at the polymer-polymer interface for applications to 3D manufacturing by following diffusion of polymer chains. We show that each component of the polymer blend has its own characteristic time of diffusion and that both strain–stress and shear viscosity curves agree with experimental data. Next, nonlinear shrinkage of the metal part during manufacturing by bound metal deposition under microgravity, is considered using a multi-scale physics-based approach that spans from atomistic dynamics to full-part shrinkage. At the smallest scale we use atomistic dynamics to simulated sintering of a few Ti6Al4V nanoparticles, study sintering mechanisms and estimate important parameters including grain boundary width and diffusion rates of the atoms at the surface that are used at the other scales of modelling. Finally, we analyse selective ionic conduction through an artificial nanopore in a graphene sheet motivated by applications to desalination and energy harvesting. We build distributions of the number and orientation of water molecules in the ion’s hydration shells providing novel insight into the ion-water-carbon interaction at the nanopore.

molecular dynamics↗

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗

Dendritic growth in the presence of convection

The motion of the freezing front between a dendritic crystal and a supercooled liquid is studied using an interface evolution equation derived from a boundary integral transformation of the transient convective-diffusion equation. A new steady-state theory is introduced that incorporates the effects of convection in dendritic growth. It is shown that in the absence of capillary effects the shape of the crystal-melt interface is a paraboloid of revolution, similar to that found in situations where diffusion is the sole heat transfer mechanism. A relation between the supercooling, the product of the tip velocity and tip radius, and the strength of the flow is derived which reduces to the well-known Ivantsov theory in the absence of convection. A non-linear interface-tracking algorithm is developed and used to study the temporal and spatial evolution of the dendritic interface. The important role of capillarity and convection on the interface dynamics is established and the response of the interface to finite amplitude disturbances is examined for the first time. Tip splitting is identified as the dominant destabilization mechanism in the limit of zero surface tension. Finite surface tension leads to interface stabilization, irrespective of the magnitude and structure of the external perturbations. Finally, convection significantly decreases the magnitude of the freezing velocity.

Beaghton, Pantelis John↗

Growth of Megaspherulites In a Rhyolitic Vitrophyre

Megaspherulites occur in the middle zone of a thick sequence of rhyolitic vitrophyre that occupies a small, late Eocene to early Oligocene volcanic-tectonic basin near Silver Cliff, Custer County, Colorado. Diameters of the megaspherulites range from 0.3 m to over 3.66 m, including a clay envelope. The megaspherulites are compound spherulites. consisting of an extremely large number (3.8 x 10(exp 9) to 9.9 x 10(exp 9)) of individual growth cones averaging 3 mm long by 1.25 mm wide at their termination. They are holocrystalline, very fine- to fine-grained, composed of disordered to ordered sanidine (orthoclase) and quartz, and surrounded by a thin K-feldspar, quartz rich rind, an inner clay layer with mordenite, and an outer clay layer composed wholly of 15 A montmorillonite. Whole rock analyses of the megaspherulites show a restricted composition from their core to their outer edge, with an average analyses of 76.3% SiO2, 0.34% CaO, 2.17% Na2O, 6.92% K2O, 0.83% H2O+ compared to the rhyolitic vitrophyre from which they crystallize with 71.07% SiO2, 0.57% CaO, 4.06% Na2O,4.l0% K2O, and 6.40% H2O+. The remaining oxides of Fe2O3 (total Fe), A12O3, MnO,MgO, TiO2, P2O5, Cr2O3, and trace elements show uniform distribution between the megaspherulites and the rhyolitic vitrophyre. Megaspherulite crystallization began soon after the rhyolitic lava ceased to flow as the result of sparse heterogeneous nucleation, under nonequilibrium conditions, due to a high degree of undercooling, delta T. The crystals grow with a fibrous habit which is favored by a large delta T ranging between 245 C and 295 C, despite lowered viscosity, and enhanced diffusion due to the high H2O content, ranging between 5% and 7%. Therefore, megaspherulite growth proceeded in a diffusion controlled manner, where the diffusion, rate lags behind the crystal growth rate at the crystal-liquid interface, restricting fibril lengths and diameters to the 10 micron to 15 micron and 3 micron and 8 micron ranges respectively. Once diffusion reestablishes itself at the crystallization front, a new nucleation event occurs at the terminated tips of the fibril cones and a new cone begins to develop with a similar orientation (small angle branching) to the earlier cones. During crystallization, these fibril cones impinge upon each other, resulting in fibril cone-free areas. These cone-free areas consist of coarser, fine-grained phases, dominated by quartz, which crystallized from the melt as it accumulated between the crystallizing K-feldspar fibrils of the cones. The anhydrous nature of the disordered to ordered sanidine (orthoclase) and quartz, suggests that water in the vitrophyre moved ahead of the crystallization front, resulting in a water rich fluid being enriched in Si, K, Na, Mg, Ca, Sr, Ba, and Y. The clay layers associated with the megaspherulites are therefore, the result of the deuteric alteration between the fractionated water and the vitrophyre, as indicated by the presence of the minerals mordenite and montmorillonite. This silica rich fluid also resulted in the total silicification of the megaspherulites within the upper 3 m of the vitrophyre.

Smith, Robert K.↗

Designing, Modeling, Constructing, and Testing a Flat Panel Speaker and Sound Diffuser for a Simulator

The goal of this project was to design, model, build, and test a flat panel speaker and frame for a spherical dome structure being made into a simulator. The simulator will be a test bed for evaluating an immersive environment for human interfaces. This project focused on the loud speakers and a sound diffuser for the dome. The rest of the team worked on an Ambisonics 3D sound system, video projection system, and multi-direction treadmill to create the most realistic scene possible. The main programs utilized in this project, were Pro-E and COMSOL. Pro-E was used for creating detailed figures for the fabrication of a frame that held a flat panel loud speaker. The loud speaker was made from a thin sheet of Plexiglas and 4 acoustic exciters. COMSOL, a multiphysics finite analysis simulator, was used to model and evaluate all stages of the loud speaker, frame, and sound diffuser. Acoustical testing measurements were utilized to create polar plots from the working prototype which were then compared to the COMSOL simulations to select the optimal design for the dome. The final goal of the project was to install the flat panel loud speaker design in addition to a sound diffuser on to the wall of the dome. After running tests in COMSOL on various speaker configurations, including a warped Plexiglas version, the optimal speaker design included a flat piece of Plexiglas with a rounded frame to match the curvature of the dome. Eight of these loud speakers will be mounted into an inch and a half of high performance acoustic insulation, or Thinsulate, that will cover the inside of the dome. The following technical paper discusses these projects and explains the engineering processes used, knowledge gained, and the projected future goals of this project

Dillon, Christina↗

Growth morphology with anisotropic surface kinetics

The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.

Xiao, Rong-Fu↗

Domain decomposition algorithms and computational fluid dynamics

Some of the new domain decomposition algorithms are applied to two model problems in computational fluid dynamics: the two-dimensional convection-diffusion problem and the incompressible driven cavity flow problem. First, a brief introduction to the various approaches of domain decomposition is given, and a survey of domain decomposition preconditioners for the operator on the interface separating the subdomains is then presented. For the convection-diffusion problem, the effect of the convection term and its discretization on the performance of some of the preconditioners is discussed. For the driven cavity problem, the effectiveness of a class of boundary probe preconditioners is examined.

Chan, Tony F.↗

Sulfur transformations at the hydrogen sulfide/oxygen interface in stratified waters and in cyanobacterial mats

Stratified water bodies allow the development of several microbial plates along the water column. The microbial plates develop in relation to nutrient availability, light penetration, and the distribution of oxygen and sulfide. Sulfide is initially produced in the sediment by sulfate-reducing bacteria. It diffuses along the water column creating a zone of hydrogen sulfide/oxygen interface. In the chemocline of Solar Lake oxygen and sulfide coexist in a 0 to 10 cm layer that moves up and down during a diurnal cycle. The microbial plate at the chemocline is exposed to oxygen and hydrogen sulfide, alternating on a diurnal basis. The cyanobacteria occupying the interface switch from anoxygenic photosynthesis in the morning to oxygenic photosynthesis during the rest of the day which results in a temporal build up of elemental sulfur during the day and disappears at night due to both oxidation to thiosulfate and sulfate by thiobacilli, and reduction to hydrogen sulfide by Desulfuromonas sp. and anaerobically respiring cyanobacteria. Sulfate reduction was enhanced in the light at the surface of the cyanobacterial mats. Microsulfate reduction measurements showed enhanced activity of sulfate reduction even under high oxygen concentrations of 300 to 800 micrometer. Apparent aerobic SO sub 4 reduction activity is explained by the co-occurrence of H sub 2. The physiology of this apparent sulfate reduction activity is studied.

Cohen, Y.↗

Directional solidification in an imperfect furnace

The upward directional solidification of a binary alloy in which the rejected solute is more dense than the solvent is considered. The system is stably stratified both thermally and solutally. If heat is lost through the sidewalls, horizontal thermal gradients drive a steady double-diffusive convection. Solidification cells with small aspect ratios (width/length) A are considered, and asymptotic solutions to the steady double-diffusive system are obtained as A tends to 0 under the assumption heat losses are also small. When the thermal Grashof number is large, buoyancy boundary layers are at the sidewalls. The shape of the crystal-melt interface is determined by the coupled effects of local thermodynamic equilibrium, diffusion and convection.

Brattkus, K.↗

Solidification and crystal growth of solid solution semiconducting alloys

Problems associated with the solidification and crytal growth of solid-solution semiconducting alloy crystals in a terrestrial environment are described. A detailed description is given of the results for the growth of mercury cadmium telluride (HgCdTe) alloy crystals by directional solidification, because of their considerable technological importance. A series of HgCdTe alloy crystals are grown from pseudobinary melts by a vertical Bridgman method using a wide range of growth rates and thermal conditions. Precision measurements are performed to establish compositional profiles for the crystals. The compositional variations are related to compositional variations in the melts that can result from two-dimensional diffusion or density gradient driven flow effects ahead of the growth interface. These effects are discussed in terms of the alloy phase equilibrium properties, the recent high temperature thermophysical data for the alloys and the highly unusual heat transfer characteristics of the alloy/ampule/furnace system that may readily lead to double diffusive convective flows in a gravitational environment.

Lehoczky, S. L.↗

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗

The crack problem in bonded nonhomogeneous materials

The plane elasticity problem for two bonded half planes containing a crack perpendicular to the interface was considered. The effect of very steep variations in the material properties near the diffusion plane on the singular behavior of the stresses and stress intensity factors were studied. The two materials were thus, assumed to have the shear moduli mu(o) and mu(o) exp (Beta x), x=0 being the diffusion plane. Of particular interest was the examination of the nature of stress singularity near a crack tip terminating at the interface where the shear modulus has a discontinuous derivative. The results show that, unlike the crack problem in piecewise homogeneous materials for which the singularity is of the form r/alpha, 0 less than alpha less than 1, in this problem the stresses have a standard square-root singularity regardless of the location of the crack tip. The nonhomogeneity constant Beta has, however, considerable influence on the stress intensity factors.

Erdogan, Fazil↗