Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 631 records · Page 35

Evolution Models of Helium White Dwarf–Main-sequence Star Merger Remnants

It is predicted that orbital decay by gravitational-wave radiation and tidal interaction will cause some close binary stars to merge within a Hubble time. The merger of a helium-core white dwarf with a main-sequence (MS) star can produce a red giant branch star that has a low-mass hydrogen envelope when helium is ignited and thus become a hot subdwarf. Because detailed calculations have not been made, we compute post-merger models with a stellar evolution code. We find the evolutionary paths available to merger remnants and find the pre-merger conditions that lead to the formation of hot subdwarfs. We find that some such mergers result in the formation of stars with intermediate helium-rich surfaces. These stars later develop helium-poor surfaces owing to diffusion. Combining our results with a model population and comparing to observed stars, we find that some observed intermediate helium-rich hot subdwarfs can be explained as the remnants of the mergers of helium-core white dwarfs with low-mass MS stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Frost Growth and Densification in Laminar Flow Over Flat Surfaces

One-dimensional frost growth and densification in laminar flow over flat surfaces has been theoretically investigated. Improved representations of frost density and effective thermal conductivity applicable to a wide range of frost circumstances have been incorporated. The validity of the proposed model considering heat and mass diffusion in the frost layer is tested by a comparison of the predictions with data from various investigators for frost parameters including frost thickness, frost surface temperature, frost density and heat flux. The test conditions cover a range of wall temperature, air humidity ratio, air velocity, and air temperature, and the effect of these variables on the frost parameters has been exemplified. Satisfactory agreement is achieved between the model predictions and the various test data considered. The prevailing uncertainties concerning the role air velocity and air temperature on frost development have been elucidated. It is concluded that that for flat surfaces increases in air velocity have no appreciable effect on frost thickness but contribute to significant frost densification, while increase in air temperatures results in a slight increase the frost thickness and appreciable frost densification.

Kandula, Max↗

The iterative thermal emission method: A more implicit modification of IMC

For over 40 years, the Implicit Monte Carlo (IMC) method has been used to solve challenging problems in thermal radiative transfer. These problems typically contain regions that are optically thick and diffusive, as a consequence of the high degree of “pseudo-scattering” introduced to model the absorption and reemission of photons from a tightly-coupled, radiating material. IMC has several well-known features that could be improved: a) it can be prohibitively computationally expensive, b) it introduces statistical noise into the material and radiation temperatures, which may be problematic in multiphysics simulations, and c) under certain conditions, solutions can be nonphysical, in that they violate a maximum principle, where IMC-calculated temperatures can be greater than the maximum temperature used to drive the problem. We have developed a variant of IMC called iterative thermal emission IMC, which is designed to have a reduced parameter space in which the maximum principle is violated. ITE IMC is a more implicit version of IMC in that it uses the information obtained from a series of IMC photon histories to improve the estimate for the end of time step material temperature during a time step. A better estimate of the end of time step material temperature allows for a more implicit estimate of other temperature-dependent quantities: opacity, heat capacity, Fleck factor (probability that a photon absorbed during a time step is not reemitted) and the Planckian emission source. We have verified the ITE IMC method against 0-D and 1-D analytic solutions and problems from the literature. These results are compared with traditional IMC. We perform an infinite medium stability analysis of ITE IMC and show that it is slightly more numerically stable than traditional IMC. We find that significantly larger time steps can be used with ITE IMC without violating the maximum principle, especially in problems with non-linear material properties. The ITE IMC method does however yield solutions with larger variance because each sub-step uses a different Fleck factor (even at equilibrium)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

Transport phenomena in the crystallization of lysozyme by osmotic dewatering and liquid-liquid diffusion in low gravity

Two methods of protein crystallization, osmotic dewatering and liquid-liquid diffusion, like the vapor diffusion (hanging-drop and sessile-drop) methods allow a gradual approach to supersaturation conditions. The crystallization of hen egg-white lysozyme, an extensively characterized protein crystal, in the presence of sodium chloride was used as an experimental model with which to compare these two methods in low gravity and in the laboratory. Comparisons of crystal growth rates by the two methods under the two conditions have, to date, indicated that the rate of crystal growth by osmotic dewatering is nearly the same in low gravity and on the ground, while much faster crystal growth rates can be achieved by the liquid-liquid diffusion method in low gravity.

Todd, Paul↗

Effects of Gravity on the Double-Diffusive Convection during Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe

General 2-D and 3-D finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe in terrestrial and microgravity conditions. Parametric research was undertaken to investigate effects of gravity level, gravity vector orientation and growth velocity on the pattern of melt convection, shape of crystal/melt interface and radial thermal gradient. Verification of the model was undertaken by comparison with previously published results. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures obtained from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. For gravity levels higher then 10(exp -7) of terrestrial one it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by 50% by choosing proper orientation of the ampoule. The predicted interface shape is in agreement with one obtained experimentally by quenching.

Bune, Andris↗

Efficient Reformulation of Solid Phase Diffusion in Electrochemical-Mechanical Coupled Models for Lithium-Ion Batteries: Effect of Intercalation Induced Stresses

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms with an appropriate discretization or approximation for the solid phase diffusion within the electrode particle. One of the major difficulties in simulating Li-ion battery models is the need for simulating solid-phase diffusion in the second radial dimension r within the particle. It increases the complexity of the model as well as the computation time/cost to a great extent. This is Particularly true for the inclusion of pressure induced diffusion inside particles experiencing volume change. A computationally efficient representation for solid-phase diffusion is discussed in this paper. The operating condition has a significant effect on the validity, accuracy, and efficiency of various approximations for the solid-phase transport governed by pressure induced diffusion. This paper introduces efficient methods for solid phase reformulation - (1) parabolic profile approach and (2) a mixed order finite difference method for approximating/representing solid-phase concentration variations within the active materials of porous electrodes for macroscopic models for lithium-ion batteries. (C) 2013 The Electrochemical Society. All rights reserved.

De, S↗

Theory and observation of a dynamically evolving negative ion plasma

The study described here examines the full range of negative-ion chemistry in the upper ionosphere by using current reaction rate data to investigate the many chemical paths SF6-type injections might take in an F region environment. Special attention is given to the conditions required to create heavy negative ions that persist long enough to affect the dynamical properties of the F region. The ambipolar diffusion characteristics of a three-component plasma are described, and estimates are given of the incoherent scatter spectra obtained from such a plasma. Model calculations using a first-order chemical code are defined and tested in order to investigate the actual types of negative-ion plasma capable of being created under nighttime conditions.

Mendillo, M.↗

Analytic Couple Modeling Introducing Device Design Factor, Fin Factor, Thermal Diffusivity Factor, and Inductance Factor

A set of convenient thermoelectric device solutions have been derived in order to capture a number of factors which are previously only resolved with numerical techniques. The concise conversion efficiency equations derived from governing equations provide intuitive and straight-forward design guidelines. These guidelines allow for better device design without requiring detailed numerical modeling. The analytical modeling accounts for factors such as i) variable temperature boundary conditions, ii) lateral heat transfer, iii) temperature variable material properties, and iv) transient operation. New dimensionless parameters, similar to the figure of merit, are introduced including the device design factor, fin factor, thermal diffusivity factor, and inductance factor. These new device factors allow for the straight-forward description of phenomenon generally only captured with numerical work otherwise. As an example a device design factor of 0.38, which accounts for thermal resistance of the hot and cold shoes, can be used to calculate a conversion efficiency of 2.28 while the ideal conversion efficiency based on figure of merit alone would be 6.15. Likewise an ideal couple with efficiency of 6.15 will be reduced to 5.33 when lateral heat is accounted for with a fin factor of 1.0.

Seebeck effect↗

Triso Analysis Tool For Predictive Source Terms

Source term modeling for TRi-structural ISOtropic (TRISO) fuel has been performed for previous reactor designs, but few are available in the open literature. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that can be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. To meet this need, the TRISO Analysis Tool for Predictive Source terms (TRISO-ATOPS) was developed. This model calculates the release of the key safety-important fission products by diffusion through the kernel, silicon carbide and graphite based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system for removing fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. From this initial condition, the model calculates the fission product release for any transient temperature profile, and the fission product releases can then be used to assess radiological dose to the workers and the public using conventional radiological dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy Advanced Gas Reactor TRISO fuel development program. The example cases in this work demonstrate the flexibility of the model

Stoyer, Benjamin [Idaho National Laboratory (INL),↗

The Fundamental Creep Behavior Model of Gr.91 Alloy by ICME Approach (Final Report)

The current project mainly focused on the fundamental creep cracking mechanism of the Grade 91 system during the operation conditions in the advanced technologies FE power plants and build the link among Composition-Processing Parameters-Phase Stability-Microstructure-Creep Resistance. A model based on computational thermodynamics and diffusion kinetics will be developed to provide guidance on how to improve the creep resistance of the alloy system. The PI mainly investigate the Gr.91 base alloy and weldment with the Integrated Computational Materials Engineering (ICME) approach. The long-term goal of the proposed program is to develop a model for creep and fatigue resistant alloys with different elemental systems, compositions, and processing parameters of weldment and heat-treatment. In addition, new creep and fatigue resistant alloys are to be designed based on the predictions from the model. The computational as well as experimental results include three sections: 1) The compositional optimization and secondary phases evaluation regarding the creep resistance in Grade 91 steel through the CALPHAD approach. In this section, the formation of the critical secondary phases, M23C6, MX and Z-phase were well predicted based on many conditions like the temperature as well as compositions. Meanwhile, the role of different alloying elements was also considered as they will directly affect the formation of the secondary phases mentioned above based on the compositional evolution of the alloying elements. 2) An investigation of Creep Resistance in Grade 91 Steel through Computational Thermodynamics. We systematically considered the evolution of the four critical temperatures, Ac3, Ac1, the threshold of M23C6 and Z-phase, on basis of the computational approach. Meanwhile, the equilibrium cooling as well as Scheil simulations were also considered to perform further predictions based on the different cooling rates. 3) Creep lifetime and microstructural evolution of the Grade 91 steel. In this part, mainly experimental approach was performed on the real creep test from the beginning of the commercial Gr.91 alloys. The alloys were normalized, tempered, welded (Gleeble) and PWHT and finally subjected to the creep lifetime analyses. Finally, all the samples after each process were further characterized based on the OM, XRD and SEM technique to draw final conclusions.

01 COAL, LIGNITE, AND PEAT↗

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗

Combined Experimental and Modeling Study of the Interactions of Acid Gas with Common Spacecraft Surfaces for Fire Safety Applications

A fire in a spacecraft poses detrimental consequences and risks mission success in addition to crew safety. This is compounded during long-duration missions when the crew has limited options to recover from a fire. A common spacecraft fire concern is the smoldering of wire insulation, typically made from Polyvinyl chloride (PVC) or Polytetrafluoroethylene (PTFE). This creates acid gases such as Hydrogen Chloride (HCl), Hydrogen Fluoride (HF) and Hydrogen Cyanide (HCN). These poisonous gases are hazardous to the crew. They also interact with common surfaces within the spacecraft more than dominant combustion products such as CO2 and H2O. This makes them more difficult to track for potential fire detection techniques, or for postfire clean-up. It is imperative to be able to understand and predict the fate of these poisonous species in a microgravity environment in order to design a safe vehicle. HCl interacts with a number of materials inside a spacecraft. Primary among these materials is aluminum, which is abundantly used due to its strong and light weight nature. Aluminum has a natural oxide layer that protects it from corrosion but is typically treated to enhance this oxide layer. Among these treatments is a chromate conversion coating (CCC), which provides a thin enough protective oxide layer to still conduct electricity, and a traditional anodized material that has a thicker oxide layer that does not conduct electricity. Nomex is another common material found inside a spacecraft. It is a flame-resistant woven polymer that is related to nylon. This commercially available material is used for cargo storage bags and as a fire barrier. Physics-based models were developed to predict the uptake of HCl by these materials. The ultimate objective of these models is to predict the fate of HCl within the spacecraft so that sensors can be placed in meaningful locations in future missions based on the model predictions. To support these modeling efforts, experiments were performed in a cast acrylic test cell that measured the difference between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. Different uptake capacities were realized for each type of sample tested. A computational fluid dynamics model (CFD) model of the reactor was then constructed that used a one-step global reaction rate with calibratable reaction (or kinetic) constants. These constants were calibrated to match the HCl uptake on the CCC aluminum samples, and the same kinetic constants were then tested for the stock and anodized aluminum samples. Model predictions matched the experimental data for the stock aluminum, and to a much lesser extent, the anodized aluminum. The model was additionally validated at different flow rates, sample surface areas, and inlet concentrations, and showed good agreement for all stock and CCC samples. The model did not accurately predict the HCl uptake in the anodized samples compared to the other two types of aluminum. Adjusting the kinetic constants and transport properties did little to improve the prediction. X-Ray Photoelectron Spectroscopy (XPS) was used to determine that the oxide layer thickness of anodized aluminum is approximately 5,000 nm, compared to 250 nm for CCC and 50 nm for stock. XPS also revealed presence of chlorine further down in the aluminum oxide layer in anodized samples than CCC and stock samples after the samples were saturated with HCl, indicating that accounting for diffusion of HCl into the oxide layer is important for accurate prediction of HCl uptake onto anodized aluminum. Consequently, a multi-scale model was developed and tested. First, a single pore inside the anodized aluminum oxide layer was modeled and is referred to as the pore-scale model. In this model, HCl diffused through the pore and reacted with the aluminum oxide pore wall to create aluminum chloride. The sample was then saturated when the mass transfer resistance through the growing aluminum chloride layer became too large for the HCl to reach the aluminum oxide wall and continue the reaction. This pore-scale model was coupled to the reactor-scale model using a concentration-dependent diffusion coefficient, resulting in much more accurate predictions (approximately half the sum square error of the aforementioned reactor-scale model that produced good agreement for stock and CCC) for a variety of operating conditions. The amount of water vapor or relative humidity (RH) in the flow during a reactor experiment was determined to influence HCl uptake. Experiments were performed to understand the interaction of gaseous HCl with aluminum surfaces in the presence of water vapor. The results show that increasing levels of RH increased the capacity of aluminum to adsorb HCl but decreased the capacity of Nomex to uptake HCl. A series of tests were performed on individual aluminum samples after they had been saturated with a fixed concentration of HCl in dry air conditions with the goal of determining how their HCl uptake capacity changes after various treatments with water relative to the original saturation tests. HCl-saturated aluminum samples subjected to a second dry air flow at the same HCl concentration as the original test had an uptake of 23.5% of the original sample with no treatment in between. Saturated aluminum samples subjected to an air flow with a RH of 90% in between tests had an uptake of 35.6% of the original. Saturated aluminum samples submerged in distilled water for 12 hours in between tests had an uptake of 82.2% of the original sample. Previously saturated aluminum tested with HCl and a 50% RH air flow resulted in similar uptake characteristics in multiple repeated tests. The experiments show the profound effect water vapor has on HCl uptake onto aluminum surfaces. In the samples subjected to water vapor or liquid water, capillary condensation and capillary diffusion alters the transport of HCl significantly. A model was proposed that developed a relationship between RH and the coefficient of HCl diffusion in aluminum chloride. This produced an “S-shaped” curve with diffusion coefficient as a function of RH, with 45% RH represented as the point where the diffusion coefficient is halfway between no water saturation and 100% water saturation in the aluminum chloride product layer. No difference in uptake characteristics for the experiment or model were realized between 50% and 62% RH. The results from the large-scale microgravity experiment, Saffire, are discussed as they pertain to the fate of HCl throughout a spacecraft. HCl was released, both as a standalone event, and in concurrence with the burning of a structured cloth. These events only produced a small response in the far field HCl sensor, while a PMMA burn that did not produce HCl had a significantly greater response. A ground-based large-scale facility was constructed to flow acid gas at the scale and configuration realized in the Saffire experiments. A CFD model of this duct was constructed to test kinetic parameters developed in this work at a larger scale and different geometric configuration and to predict the results of the large-scale facility. The models developed in this work were used to interpret the results of the microgravity tests and lead the discussion on what further experiments and models are needed in order to predict the fate of acid gas in a spacecraft environment. To summarize, the major contributions of this work are as follows: the capacity to uptake HCl, with and without the presence of water vapor, was measured for a variety of real spacecraft surfaces. Several different models (single reactor-scale, multiscale, spacecraft-scale) were developed and with the aid of modeling, the rate of uptake for those surfaces was also predicted and validated. The kinetic parameters determined from the small-scale reactor experiments and models were used to predict large-scale and microgravity tests. Conclusions from this research will be used in the design of spacecraft vehicles and large-scale microgravity fire safety experiments. The models built by this work will aid designers in sensor placement and could be used to predict acid gas transport from fires in partial gravity, as would be seen in Lunar and Martian habitats.

fire safety↗

Cosmic-ray gradients in the outer heliosphere

We report on a study of the expected spatial variation of the galactic cosmic-ray intensity in the vicinity of the solar-wind termination shock. Model simulations predict that the radial gradients change abruptly at the shock, and that the nature of the effect varies significantly with particle energy. At low energies, the radial gradient changes abruptly from a lower value inside the shock to a higher value outside, whereas at high energies, the higher value of the gradient is inside the shock. This effect, which is a consequence of the matching conditions at the shock and is closely related to diffusive shock acceleration, is qualitatively the same for both heliospheric magnetic polarity states and remains much the same in one-dimensional, two-dimensional and three-dimensional models. Hence drifts do not change the nature of this phenomenon, although they change it quantitatively. The effect may prove to be an important diagnostic tool for the study of the termination of the solar wind.

Kota, J.↗

Benchmark analytical solutions for steady state high frequency broadband sound fields in three rectangular enclosures

Analytical solutions are presented for broadband sound fields in three rectangular enclosures with absorption applied on the floor and ceiling, rigid sidewalls, and a vertically oriented dipole source. The solutions are intended to serve as benchmarks that can be used to assess the performance of broadband techniques, particularly energy-based methods, in a relatively straightforward configuration with precisely specified boundary conditions. A broadband Helmholtz solution is developed using a frequency-by-frequency modal approach to determine the exact band averaged mean-square pressures along spatial trajectories within each enclosure. Due to the specific choice of enclosure configuration and absorption distribution, an approximate specular solution can be obtained through a summation of uncorrelated image sources. Comparisons between the band averaged Helmholtz solution and the uncorrelated image solution reveal excellent agreement for a wide range of absorption levels and improve the understanding of correlation effects in broadband sound fields. In conclusion, a boundary element solution with diffuse boundaries is also presented, which produces consistently higher mean-square pressures in comparison with the specular solution, emphasizing the careful attention that must be placed on correctly modeling reflecting boundary conditions and demonstrating the errors that can result from assuming a Lambertian surface.

42 ENGINEERING↗

Parametric studies with an atmospheric diffusion model that assesses toxic fuel hazards due to the ground clouds generated by rocket launches

Parametric studies were made with a multilayer atmospheric diffusion model to place quantitative limits on the uncertainty of predicting ground-level toxic rocket-fuel concentrations. Exhaust distributions in the ground cloud, cloud stabilized geometry, atmospheric coefficients, the effects of exhaust plume afterburning of carbon monoxide CO, assumed surface mixing-layer division in the model, and model sensitivity to different meteorological regimes were studied. Large-scale differences in ground-level predictions are quantitatively described. Cloud alongwind growth for several meteorological conditions is shown to be in error because of incorrect application of previous diffusion theory. In addition, rocket-plume calculations indicate that almost all of the rocket-motor carbon monoxide is afterburned to carbon dioxide CO2, thus reducing toxic hazards due to CO. The afterburning is also shown to have a significant effect on cloud stabilization height and on ground-level concentrations of exhaust products.

Stewart, R. B.↗

Strain-Rate-Free Diffusion Flames: Initiation, Properties, and Quenching

For about a half century, the stabilization of a steady planar deflagration on a heat-sink-type flat-flame burner has been of extraordinary service for the theoretical modeling and diagnostic probing of combusting gaseous mixtures. However, most engineering devices and most unwanted fire involve the burning of initially unmixed reactants. The most vigorous burning of initially separated gaseous fuel and oxidizer is the diffusion flame. In this useful idealization (limiting case), the reactants are converted to product at a mathematically thin interface, so no interpenetration of fuel and oxidizer occurs. This limit is of practical importance because it often characterizes the condition of optimal performance (and sometimes environmentally objectionable operation) of a combustor. A steady planar diffusion flame is most closely approached in the laboratory in the counterflow apparatus. The utility of this simple-strain-rate flow for the modeling and probing of diffusion flames was noted by Pandya and Weinberg 35 years ago, though only in the last decade or so has its use become internationally common place. However, typically, as the strain rate a is reduced below about 20 cm(exp -1), and the diffusion-flame limit (reaction rate much faster than the flow rate) is approached, the burning is observed to become unstable in earth gravity. The advantageous steady planar flow is not available in the diffusion-flame limit in earth gravity. This is unfortunate because the typical spatial scale in a counterflow is (k/a)(sup 1/2), where k denotes a characteristic diffusion coefficient; thus, the length scale becomes large, and the reacting flow is particularly amenable to diagnostic probing, as the diffusion-flame limit is approached. The disruption of planar symmetry is owing the fact that, as the strain rate a decreases, the residence time (l/a) of the throughput in the counterflow burner increases. Observationally, when the residence time exceeds about 50 msec, the inevitably present convective (Rayleigh-Benard) instabilities, associated with hot-under-cold (flame-under-fresh-reactant) stratification of fluid in a gravitational field, have time to grow to finite amplitude during transit of the burner.

Fendell, Francis↗

Elucidating the role of Cr migration in Ni-Cr exposed to molten FLiNaK via multiscale characterization

Structural materials used in nuclear reactor environments are subjected to coupled extremes such as high temperature, irradiation, and corrosion which act in concert to degrade their functional performance. Connecting alloy microstructure such as grain boundaries, and accumulating point defects with corrosion attack and pore morphology is critical to understanding underlying mechanisms. We uncover the compositional variations and morphology at multiple length scales in corrosion-damaged Ni-Cr alloy after exposure to oxidants in molten fluoride salts. A complex network of dense corrosion pores is detected by surface-level SEM observations. The corrosion pores take on a 1-dimensional morphology and are enriched with Ni and depleted of Cr 1–5 µm from the pore surface. STEM-EDS and 4D-STEM strain maps acquired simultaneously highlight the local variations in composition and structure of a ≤ 200 nm Cr-rich layer identified from a cross-section taken at the bottom of an isolated corrosion pore between the Ni-Cr alloy matrix and the embedded salt. However, the absence of an observed interface between the Ni-Cr alloy matrix and the FCC-structured Cr-rich layer suggests that Cr plating from the salt did not transpire. These findings support a proposed Cr lattice diffusion mechanism rather than Cr-precipitation from the salt to accommodate temperature transient conditions during sample cooling. Through the development of a 1D phase field model, these results are rationalized by formation energies for the Ni- and Cr-oxidation into the molten salt. This study reveals the locally altered microstructure caused by high temperature corrosion in non-steady-state molten salt nuclear reactor environments.

36 MATERIALS SCIENCE↗