Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compositional simulation”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Compiled records of carbon isotopes in atmospheric CO 2 for historical simulations in CMIP6

Abstract. The isotopic composition of carbon (Δ 14 C and δ 13 C) in atmospheric CO 2 and in oceanic and terrestrial carbon reservoirs is influenced by anthropogenic emissions and by natural carbon exchanges, which can respond to and drive changes in climate. Simulations of 14 C and 13 C in the ocean and terrestrial components of Earth system models (ESMs) present opportunities for model evaluation and for investigation of carbon cycling, including anthropogenic CO 2 emissions and uptake. The use of carbon isotopes in novel evaluation of the ESMs' component ocean and terrestrial biosphere models and in new analyses of historical changes may improve predictions of future changes in the carbon cycle and climate system. We compile existing data to produce records of Δ 14 C and δ 13 C in atmospheric CO 2 for the historical period 1850–2015. The primary motivation for this compilation is to provide the atmospheric boundary condition for historical simulations in the Coupled Model Intercomparison Project 6 (CMIP6) for models simulating carbon isotopes in the ocean or terrestrial biosphere. The data may also be useful for other carbon cycle modelling activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Where Big Data and Prediction Meet

Our ability to assemble and analyze massive data sets, often referred to under the title of “big data”, is an increasingly important tool for shaping national policy. This in turn has introduced issues from privacy concerns to cyber security. But as IBM’s John Kelly emphasized in the last Innovation, making sense of the vast arrays of data will require radically new computing tools. In the past, technologies and tools for analysis of big data were viewed as quite different from the traditional realm of high performance computing (HPC) with its huge models of phenomena such as global climate or supporting the nuclear test moratorium. Looking ahead, this will change with very positive benefits for both worlds. Societal issues such as global security, economic planning and genetic analysis demand increased understanding that goes beyond existing data analysis and reduction. The modeling world often produces simulations that are complex compositions of mathematical models and experimental data. This has resulted in outstanding successes such as the annual assessment of the state of the US nuclear weapons stockpile without underground nuclear testing. Ironically, while there were historically many test conducted, this body of data provides only modest insight into the underlying physics of the system. A great deal of emphasis was thus placed on the level of confidence we can develop for the predictions. As data analytics and simulation come together, there is a growing need to assess the confidence levels in both data being gathered and the complex models used to make predictions. An example of this is assuring the security or optimizing the performance of critical infrastructure systems such as the power grid. If one wants to understand the vulnerabilities of the system or impacts of predicted threats, full scales tests of the grid against threat scenarios are unlikely. Preventive measures would need to be predicated on well-defined margins of confidence in order to take mitigating actions that could have wide ranging impacts. There is a rich opportunity for interaction and exchange between the HPC simulation and data analytics communities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Polyacrylonitrile Composites of Ag–Al–Si–O Aerogels and Xerogels as Iodine and Iodide Sorbents

Mechanically robust engineered sorbents were constructed by embedding Ag-loaded aerogels or xerogels in a polyacrylonitrile (PAN) matrix at 30 or 50 mass% loadings of gel in PAN. Here, the process for making the base gels as well as the composite materials is described herein, a technology that is very flexible for making hybrid sorbents using various inorganic or organic active sorbents. The PAN-gel composite materials were evaluated for iodine uptake in the form of gaseous iodine [I 2(g) ] at room temperature and at 150°C as well as iodide (I – ) at room temperature from simulated liquid nuclear wastes. The composites show high I 2(g) uptake and high I – uptake (= 93 mass%) in simulated liquid secondary waste spiked with 5 mg L -1 I – . These results demonstrate the utility of such a composite form for iodine removal from contaminated gas and liquid streams.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Experimental Testing and Numerical Modeling of Ballistic Impact on S-2 Glass/SC 15 Composites

Here, the work presented in this paper details both an experimental program and an associated numerical modeling effort to characterize and predict the ballistic response of S-2 glass/SC15 epoxy composite panels. The experimental program consisted of ¼ inch diameter soft carbon steel spheres impacting ¼ and ½ inch thick flat composite panels at velocities ranging from 220 to 1570 m/s. High speed cameras were used to capture the impact event and resulting residual velocity of the spheres for each test configuration. After testing, each panel was inspected both visually and with ultrasonic C-scan techniques to determine the extent and depth of damage imparted on the panel by the impactor. The numerical modeling efforts utilized the anisotropic multi-constituent composite model (MCM) within the CTH shock physics hydrocode. The MCM model allows for evaluation of damage at the constituent level through continuum averaged stress and strain fields. The model also accounts for the inherent coupling of the equation of state and strength response that occurs in anisotropic composite materials. Finally, the simulation results are compared against the experimentally measured residual velocity as a quantitative metric and against the measured damage extent and patterns as a qualitative metric. The comparisons show good agreement in residual velocity and damage extent.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Atmospheric-pressure electric discharge as an instrument of chemical activation of water solutions

Results of experimental studies and numerical simulations of physicochemical characteristics of plasmas generated in different types of atmospheric-pressure discharges (pulsed streamer corona, gliding electric arc, dielectric barrier discharge, glow-discharge electrolysis, diaphragmatic discharge, and dc glow discharge) used to initiate various chemical processes in water solutions are analyzed. Typical reactor designs are considered. Data on the power supply characteristics, plasma electron parameters, gas temperatures, and densities of active particles in different types of discharges excited in different gases and their dependences on the external parameters of discharges are presented. The chemical composition of active particles formed in water is described. Possible mechanisms of production and loss of plasma particles are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Effect of feed composition on the production of off-gases during vitrification of simulated low-activity nuclear waste

During the vitrification of nuclear waste, hazardous and radioactive emissions are generated from the feed-to-glass conversion reactions, in addition to discharges from forced air bubbling and air inleakage. Although the major gaseous emissions are water vapor, nitrogen, and carbon dioxide, various monitored environmental pollutants are also released, such as nitrogen oxides or sulfur dioxide. In addition, reactions between organics and nitrates in the feed may also form products of incomplete combustion such as carbon monoxide and acetonitrile. Although off-gas emissions are commonly measured during both laboratory- and pilot-scale melter testing, no predictive tool is currently available to a priori estimate the composition of gaseous emissions during nuclear waste vitrification. This work forms a basis for the development of such predictive tool by measuring gas evolution from a broad range of simulated low-activity waste melter feeds using evolved gas analysis data and developing correlations between the feed and off-gas compositions. Using reaction stoichiometry and regression analysis, we demonstrate that next to the content of nitrogen and organic carbon in the feed, the gaseous emissions are affected by the feed reduction-oxidation conditions – the more the feed is reduced, the less nitrogen monoxide, and more carbon monoxide and acetonitrile evolves. In conclusion, the results presented in this work provide a first step towards reducing the amount of expensive physical melter testing and the regression analysis provides a simple tool for rapid optimization of feed composition with respect to off-gas composition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Modeling magnetic evolution and exchange hardening in disordered magnets: The example of Mn 1 - x Fe x Ru 2 Sn Heusler alloys

We demonstrate how exchange hardening can arise in a chemically disordered solid solution from a first-principles statistical mechanics approach. A general mixed-basis chemical and magnetic cluster expansion has been developed, and applied to the Mn 1-x Fe x Ru 2 Sn Heusler alloy system; single-phase solid solutions between antiferromagnetic MnRu 2 Sn and ferromagnetic FeRu 2 Sn with disorder on the Mn/Fe sublattice that exhibit unexpected exchange hardening. Monte Carlo simulations applied to the cluster expansion are able to reproduce the experimentally measured magnetic transition temperatures and the bulk magnetization as a function of composition. The magnetic ordering around a site is shown to be dependent not only on bulk composition, but also on the identity of the site and the local composition around that site. The simulations predict that local antiferromagnetic orderings form inside a bulk ferromagnetic region at intermediate compositions that drives the exchange hardening. Furthermore, the antiferromagnetic regions disorder at a lower temperature than the ferromagnetic regions, providing an atomistic explanation for the experimentally observed decrease in exchange hardening with increasing temperature. These effects occur on a length scale too small to be resolved with previously used characterization techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2017

BeamDyn: a high-fidelity wind turbine blade solver in the FAST modular framework

Here, this paper presents a numerical implementation of the geometrically exact beam theory based on the Legendre-spectral-finite-element (LSFE) method. The displacement-based geometrically exact beam theory is presented, and the special treatment of three-dimensional rotation parameters is reviewed. An LSFE is a high-order finite element with nodes located at the Gauss-Legendre-Lobatto points. These elements can be an order of magnitude more computationally efficient than low-order finite elements for a given accuracy level. The new module, BeamDyn, is implemented in the FAST modularization framework for dynamic simulation of highly flexible composite-material wind turbine blades within the FAST aeroelastic engineering model. The framework allows for fully interactive simulations of turbine blades in operating conditions. Numerical examples are provided to validate BeamDyn and examine the LSFE performance as well as the coupling algorithm in the FAST modularization framework. BeamDyn can also be used as a stand-alone high-fidelity beam tool.

17 WIND ENERGY↗
DOE OSTI2015

Stealth Dark Matter: Dark scalar baryons through the Higgs portal

We present a new model of "Stealth Dark Matter": a composite baryonic scalar of an SU(N D ) strongly coupled theory with even N D ≥ 4. All mass scales are technically natural, and dark matter stability is automatic without imposing an additional discrete or global symmetry. Constituent fermions transform in vectorlike representations of the electroweak group that permit both electroweak-breaking and electroweak-preserving mass terms. This gives a tunable coupling of stealth dark matter to the Higgs boson independent of the dark matter mass itself. We specialize to SU(4), and investigate the constraints on the model from dark meson decay, electroweak precision measurements, basic collider limits, and spin-independent direct detection scattering through Higgs exchange. We exploit our earlier lattice simulations that determined the composite spectrum as well as the effective Higgs coupling of stealth dark matter in order to place bounds from direct detection, excluding constituent fermions with dominantly electroweak-breaking masses. A lower bound on the dark baryon mass m B ≳ 300 GeV is obtained from the indirect requirement that the lightest dark meson not be observable at LEP II. Furthermore, we briefly survey some intriguing properties of stealth dark matter that are worthy of future study, including collider studies of dark meson production and decay; indirect detection signals from annihilation; relic abundance estimates for both symmetric and asymmetric mechanisms; and direct detection through electromagnetic polarizability, a detailed study of which will appear in a companion paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Dosimetric characterization and organ dose assessment in digital breast tomosynthesis: Measurements and Monte Carlo simulations using voxel phantoms

Purpose: Due to its capability to more accurately detect deep lesions inside the breast by removing the effect of overlying anatomy, digital breast tomosynthesis (DBT) has the potential to replace the standard mammography technique in clinical screening exams. However, the European Guidelines for DBT dosimetry are still a work in progress and there are little data available on organ doses other than to the breast. It is, therefore, of great importance to assess the dosimetric performance of DBT with respect to the one obtained with standard digital mammography (DM) systems. The aim of this work is twofold: (i) to study the dosimetric properties of a combined DBT/DM system (MAMMOMAT Inspiration Siemens{sup ®}) for a tungsten/rhodium (W/Rh) anode/filter combination and (ii) to evaluate organs doses during a DBT examination. Methods: For the first task, measurements were performed in manual and automatic exposure control (AEC) modes, using two homogeneous breast phantoms: a PMMA slab phantom and a 4 cm thick breast-shaped rigid phantom, with 50% of glandular tissue in its composition. Monte Carlo (MC) simulations were performed using Monte Carlo N-Particle eXtended v.2.7.0. A MC model was implemented to mimic DM and DBT acquisitions for a wide range of x-ray spectra (24 –34 kV). This was used to calculate mean glandular dose (MGD) and to compute series of backscatter factors (BSFs) that could be inserted into the DBT dosimetric formalism proposed by Dance et al. Regarding the second aim of the study, the implemented MC model of the clinical equipment, together with a female voxel phantom (“Laura”), was used to calculate organ doses considering a typical DBT acquisition. Results were compared with a standard two-view mammography craniocaudal (CC) acquisition. Results: Considering the AEC mode, the acquisition of a single CC view results in a MGD ranging from 0.53 ± 0.07 mGy to 2.41 ± 0.31 mGy in DM mode and from 0.77 ± 0.11 mGy to 2.28 ± 0.32 mGy in DBT mode. Regarding the BSF, the results achieved may lead to a MGD correction of about 6%, contributing to the improvement of the current guidelines used in these applications. Finally, considering the MC results obtained for the organ dose study, the radiation doses found for the tissues of the body other than the breast were in the range of tens of μSv, and are in part comparable to the ones obtained in standard DM. Nevertheless, in a single DBT examination, some organs (such as lung and thyroid) receive higher doses (of about 9% and 21%, respectively) with respect to the CC DM acquisition. Conclusions: Taking into account an average breast with a thickness of 4.5 cm, the MGDs for DM and DBT acquisitions were below the achievable value (2.0 mGy) defined by the European protocol. Additionally, in the case of a fusion imaging study (DM + DBT), the MGD for a 4.5 cm thick breast is of the order of 1.88 ± 0.36 mGy. Finally, organ dose evaluations underline the need to improve awareness concerning dose estimation of DBT exams for some organs, especially when radiation risk is assessed by using the effective dose.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Research in Lattice Gauge Theory and in the Phenomenology of Neutrinos and Dark Matter

Research in theoretical elementary particle physics was performed by the PI Yannick Meurice and co-PI Mary Hall Reno. New techniques designed for precision calculations of strong interaction physics were developed using the tensor renormalization group method. Large-scale Monte Carlo simulations with dynamical quarks were performed for candidate models for Higgs compositeness. Ab-initio lattice gauge theory calculations of semileptonic decays of B-mesons observed in collider experiments and relevant to test the validity of the standard model were performed with the Fermilab/MILC collaboration. The phenomenology of strong interaction physics was applied to new predictions for physics processes in accelerator physics experiments and to cosmic ray production and interactions. A research focus has been on heavy quark production and their decays to neutrinos. The heavy quark contributions to atmospheric neutrino and muon fluxes have been evaluated, as have the neutrino fluxes from accelerator beams incident on heavy targets. Results are applicable to current and future particle physics experiments and to astrophysical neutrino detectors such as the IceCube Neutrino Observatory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Impact of dry intrusion events on the composition and mixing state of particles during the winter Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA)

Abstract. Long-range transport of continental emissions has a far-reaching influence over remote regions, resulting in substantial change in the size, morphology, and composition of the local aerosol population and cloud condensation nuclei (CCN) budget. Here, we investigate the physicochemical properties of atmospheric particles collected on board a research aircraft flown over the Azores during the winter 2018 Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA) campaign. Particles were collected within the marine boundary layer (MBL) and free troposphere (FT) after long-range atmospheric transport episodes facilitated by dry intrusion (DI) events. Chemical and physical properties of individual particles were investigated using complementary capabilities of computer-controlled scanning electron microscopy and X-ray spectromicroscopy to probe particle external and internal mixing state characteristics. Furthermore, real-time measurements of aerosol size distribution, cloud condensation nuclei (CCN) concentration, and back-trajectory calculations were utilized to help bring into context the findings from offline spectromicroscopy analysis. While carbonaceous particles were found to be the dominant particle type in the region, changes in the percent contribution of organics across the particle population (i.e., external mixing) shifted from 68 % to 43 % in the MBL and from 92 % to 46 % in FT samples during DI events. This change in carbonaceous contribution is counterbalanced by the increase in inorganics from 32 % to 57 % in the MBL and 8 % to 55 % in FT. The quantification of the organic volume fraction (OVF) of individual particles derived from X-ray spectromicroscopy, which relates to the multi-component internal composition of individual particles, showed a factor of 2.06 ± 0.16 and 1.11 ± 0.04 increase in the MBL and FT, respectively, among DI samples. We show that supplying particle OVF into the κ-Köhler equation can be used as a good approximation of field-measured in situ CCN concentrations. We also report changes in the κ values in the MBL from κMBL, non-DI=0.48 to κMBL, DI=0.41, while changes in the FT result in κFT, non-DI=0.36 to κFT, DI=0.33, which is consistent with enhancements in OVF followed by the DI episodes. Our observations suggest that entrainment of particles from long-range continental sources alters the mixing state population and CCN properties of aerosol in the region. The work presented here provides field observation data that can inform atmospheric models that simulate sources and particle composition in the eastern North Atlantic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

X-ray tomography of feed-to-glass transition of simulated borosilicate waste glasses

The feed composition of a high level nuclear waste (HLW) glass melter affects the overall melting rate by influencing the chemical, thermophysical, and morphological properties of a relatively insulating cold cap layer over the molten phase where the primary feed vitrification reactions occur. Data from X ray computed tomography imaging of melting pellets comprised of a simulated high-aluminum HLW feed heated at a rate of 10°C/min reveal the distribution and morphology of bubbles, collectively known as primary foam, within this layer for various SiO 2 /(Li 2 CO 3 +H 3 BO 3 +Na 2 CO 3 ) mass fractions at temperatures between 600°C and 1040°C. To track melting dynamics, cross-sections obtained through the central profile of the pellet were digitally segmented into primary foam and a condensed phase. Pellet dimensions were extracted using Photoshop CS6 tools while the DREAM.3D software package was used to calculate pellet profile area, average and maximum bubble areas, and two-dimensional void fraction. The measured linear increase in the pellet area expansion rates – and therefore the increase in batch gas evolution rates – with SiO 2 /(Li 2 CO 3 +H 3 BO 3 +Na 2 CO 3 ) mass fraction despite an exponential increase in viscosity of the final waste glass at 1050°C and a lower total amount of gas-evolving species suggest that the retention of primary foam with large average bubble size at higher temperatures results from faster reaction kinetics rather than increased viscosity. However, viscosity does affect the initial foam collapse temperature by supporting the growth of larger bubbles. Because the maximum bubble size is limited by the pellet dimensions, larger scale studies are needed to understand primary foam morphology at high temperatures. This temperature-dependent morphological data can be used in future investigations to synthetically generate cold cap structures for use in models of heat transfer within a HLW glass melter.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Three-dimensional FIB/EBSD characterization of irradiated HfAl{sub 3}-Al composite

A thermal neutron-absorbing composite comprised of 28.4 vol.% HfAl{sub 3} intermetallic particles in an Al matrix, developed for use in a conduction-cooled thermal neutron filter, will enable fast flux irradiation of materials and fuels in an existing pressurized water test reactor. In this paper, an EBSD-FIB characterization approach is used to observe the microstructural changes of the HfAl{sub 3}-Al composite in response to simulated long-term neutron irradiation. Using a focused ion beam (FIB), the sample was fabricated to 25 μm x 25 μm x 20 μm and mounted on a transmission electron microscopy (TEM) grid. A series of operations were carried out repeatedly on the sample's top surface to prepare it for scanning electron microscopy (SEM). First, a ∼ 100-nm layer was removed by high voltage FIB milling. Then, several cleaning passes were performed on the newly exposed surface using low voltage FIB milling to improve the SEM image quality. The surface was then scanned by Electron Backscattering Diffraction (EBSD) to obtain a two-dimensional image. After 50 to 100 two-dimensional images were collected, the images were stacked to reconstruct a three-dimensional model using DREAM.3D software. Two such reconstructed three dimensional models were obtained from samples of the original and post-irradiation HfAl{sub 3}-Al composite, from which the most significant microstructural change caused by neutron irradiation appears to be the size reduction of both HfAl{sub 3} and Al grains. Possible causes include thermal expansion and related thermal strain from thermal neutron capture. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Coupling reduced‐order modeling and coarse‐grained CFD‐DEM to accelerate coal gasifier simulation and optimization

Abstract High fidelity three‐dimensional (3‐D) numerical simulations of commercial‐scale coal gasifiers are very time consuming and expensive. This article proposes a reduced‐order modeling approach that uses quasi one‐dimensional (1‐D) CFD–DEM simulation results to serve as an accurate initial condition for the 3‐D simulation, which significantly shortens the physical time from several hours to minutes to achieve steady states. The 3‐D simulation employs a previously validated coarse‐grained CFD‐DEM approach, which further accelerates the simulation. Both 1‐D and 3‐D simulations adopt a particle shrinkage scheme to correctly account for the particle volume change due to chemical reaction and thus provides an accurate bed height. Comparison with available experimental measurements of the bed temperature and synthesis gas composition is used to validate the simulation. The final syngas composition and flowrates are strongly affected by the fuel or gasification agent rate and process operating conditions as expected. This approach is then used to optimize two important operating conditions (air and steam flows) to achieve specific gasifier performance goals for the production of Fischer–Tropsch synthesis gas and production of fuel gas for gas turbines.

Yu, Jia↗
DOE OSTI2017

Compositional evolution of particle-phase reaction products and water in the heterogeneous OH oxidation of model aqueous organic aerosols

Organic compounds present at or near the surface of aqueous droplets can be efficiently oxidized by gas-phase OH radicals, which alter the molecular distribution of the reaction products within the droplet. A change in aerosol composition affects the hygroscopicity and leads to a concomitant response in the equilibrium amount of particle-phase water. The variation in the aerosol water content affects the aerosol size and physicochemical properties, which in turn governs the oxidation kinetics and chemistry. To attain better knowledge of the compositional evolution of aqueous organic droplets during oxidation, this work investigates the heterogeneous OH-radical-initiated oxidation of aqueous methylsuccinic acid (C 5 H 8 O 4 ) droplets, a model compound for small branched dicarboxylic acids found in atmospheric aerosols, at a high relative humidity of 85 % through experimental and modeling approaches. Aerosol mass spectra measured by a soft atmospheric pressure ionization source (Direct Analysis in Real Time, DART) coupled with a high-resolution mass spectrometer reveal two major products: a five carbon atom (C 5 ) hydroxyl functionalization product (C 5 H 8 O 5 ) and a C 4 fragmentation product (C 4 H 6 O 3 ). These two products likely originate from the formation and subsequent reactions (intermolecular hydrogen abstraction and carbon–carbon bond scission) of tertiary alkoxy radicals resulting from the OH abstraction occurring at the methyl-substituted carbon site. Based on the identification of the reaction products, a kinetic model of oxidation (a two-product model) coupled with the Aerosol Inorganic–Organic Mixtures Functional groups Activity Coefficients (AIOMFAC) model is built to simulate the size and compositional changes of aqueous methylsuccinic acid droplets during oxidation. Model results show that at the maximum OH exposure, the droplets become slightly more hygroscopic after oxidation, as the mass fraction of water is predicted to increase from 0.362 to 0.424; however, the diameter of the droplets decreases by 6.1 %. This can be attributed to the formation of volatile fragmentation products that partition to the gas phase, leading to a net loss of organic species and associated particle-phase water, and thus a smaller droplet size. Overall, fragmentation and volatilization processes play a larger role than the functionalization process in determining the evolution of aerosol water content and droplet size at high-oxidation stages.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Structural and Na-ion conduction characteristics of Na 3 PS x Se 4–x

The recent discovery of the isostructrual cubic Na 3 PS 4 and Na 3 PSe 4 as fast Na-ion conductors provided a general structural framework for the exploration of new sodium superionic conductors. In this work, we systematically investigated the structures and ionic conduction characteristics of a series of compounds with the general chemical formula of Na 3 PS x Se 4–x . Synthesis of Na 3 PS 4 under different conditions (e.g., temperature, reaction vessel, mass of the precursors) reveals the reactivity of the precursors with the reaction tubes, producing different polymorphs. X-ray diffraction studies on the solid solution phases Na 3 PS x Se 4–x identified a tetragonal-to-cubic phase transition with increasing Se concentration. This observation is consistent with the computed stability of the tetragonal and cubic polymorphs, where the energy difference between the two polymorphs becomes very close to zero in Se-rich compositions. Furthermore, ab initio molecular dynamic simulations suggest that the fast Na-ion conduction in Na 3 PS x Se 4–x may not be causally related with the symmetry or the composition of these phases. The formation of defects, instead, enables fast Na-ion conduction in this class of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Development of Computational Capabilities to Predict the Corrosion Wastage of Boiler Tubes in Advanced Combustion Systems

A comprehensive corrosion research project consisting of pilot-scale combustion testing and long-term laboratory corrosion study has been successfully performed. A pilot-scale combustion facility available at Brigham Young University was selected and modified to enable burning of pulverized coals under the operating conditions typical for advanced coal-fired utility boilers. Eight United States (U.S.) coals were selected for this investigation, with the test conditions for all coals set to have the same heat input to the combustor. In addition, the air/fuel stoichiometric ratio was controlled so that staged combustion was established, with the stoichiometric ratio maintained at 0.85 in the burner zone and 1.15 in the burnout zone. The burner zone represented the lower furnace of utility boilers, while the burnout zone mimicked the upper furnace areas adjacent to the superheaters and reheaters. From this staged combustion, approximately 3% excess oxygen was attained in the combustion gas at the furnace outlet. During each of the pilot-scale combustion tests, extensive online measurements of the flue gas compositions were performed. In addition, deposit samples were collected at the same location for chemical analyses. Such extensive gas and deposit analyses enabled detailed characterization of the actual combustion environments existing at the lower furnace walls under reducing conditions and those adjacent to the superheaters and reheaters under oxidizing conditions in advanced U.S. coal-fired utility boilers. The gas and deposit compositions were then carefully simulated in a series of 1000-hour laboratory corrosion tests, in which the corrosion performances of different commercial candidate alloys and weld overlays were evaluated at various temperatures for advanced boiler systems. Results of this laboratory study led to significant improvement in understanding of the corrosion mechanisms operating on the furnace walls as well as superheaters and reheaters in coal-fired boilers resulting from the coexistence of sulfur and chlorine in the fuel. A new corrosion mechanism, i.e., “Active Sulfidation Corrosion Mechanism,” has been proposed to account for the accelerated corrosion wastage observed on the furnace walls of utility boilers burning coals containing sulfur and chlorine. In addition, a second corrosion mechanism, i.e., “Active Sulfide-to-Oxide Corrosion Mechanism,” has been identified to account for the rapid corrosion attack on superheaters and reheaters. Both of the newly discovered corrosion mechanisms involve the formation of iron chloride (FeCl2) vapor from iron sulfide (FeS) and HCl, followed by the decomposition of FeCl2 via self-sustaining cycling reactions. For higher alloys containing sufficient chromium, the attack on superheaters and reheaters is dominated by Hot Corrosion in the presence of a fused salt. Furthermore, two stages of the hot corrosion mechanism have been identified and characterized in detail. The initiation of hot corrosion attack induced by molten sulfate leads to Stage 1 “acidic” fluxing and re-precipitation of the protective scale formed initially on the deposit-covered alloy surfaces. Once the protective scale is penetrated, Stage 2 Hot Corrosion is initiated, which is dominated by “basic” fluxing and re-precipitation of the scale in the fused salt. Based on the extensive corrosion information generated from this project, corrosion modeling was performed using non-linear regression analysis. As a result of the modeling efforts, two predictive equations have been formulated, one for furnace walls and the other for superheaters and reheaters. These first-of-the-kind equations can be used to estimate the corrosion rates of boiler tubes based on coal chemistry, alloy compositions, and boiler operating conditions for advanced boiler systems.

Coal-Environmental Processes↗