Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Study of Tritium Transport in Nuclear-Grade Graphite and Molten Fluoride Salt Systems

High-temperature salts are proposed as coolant in recent reactor technologies including Fluoride Salt-Cooled High Temperature Reactor (FHR) with fluoride salt as primary and secondary loop coolant, Molten Salt Reactors (MSRs) with fuel salt, and high-magnetic-field fusion machines with salt coolant blanket. The baseline coolant is {sup 7}Li{sub 2}BeF{sub 4} (flibe) due to its good thermal-hydraulic and neutronic properties. Tritium produced in flibe due to neutron irradiation makes the issue of tritium transport an area of great interest. In the temperature range of 427- 704 deg. C, tritium readily diffuses through metallic reactor structural materials, which makes tritium control and capture a key technical challenge. It is anticipated that tritium can be captured by graphite, as graphite material showed a high affinity for tritium in the Molten Salt Reactor Experiment (MSRE) due to the porous and irregular crystallite microstructure. Two types of nuclear grade graphite which used as core components for FHR: matrix graphite, which encapsulates the TRISO (Tri-structural-isotropic) fuel particles in FHR, and high purity, isotropic nuclear grade graphite IG-110 which can be used as moderator and reflector in reactor core for both fission and fusion reactor. The microstructure of matrix graphite and IG-110 is different from that of matrix graphite, due to the raw material and manufacturing processes. There are three phenomena important for predicting tritium transport in graphite-salt systems: tritium species diffusion in molten flibe; tritium species transport at the salt-graphite interface; and tritium species transport in graphite. In this summary microstructure characterization of IG-110 and matrix graphite, characterization of graphite-salt interface as well as techniques for trace tritium detection in graphite are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

PHB and Montmorillonite Clay Composites as KNO{sub 3} and NPK Support for a Controlled Release

In order to improve nutrient uptake by plants and at the same time decrease losses by leaching, polymeric compounds based on biodegradable poly(3-hydroxybutyrate) (PHB) and montmorillonite clay (MMt) were melt processed with KNO{sub 3} and NPK fertilizers (Fert). We present a comprehensive analysis of releasing profile of the potassium (K{sup +}), nitrogen (NH{sub 4}{sup +} and NO{sub 3}{sup −}) and phosphorous (PO{sub 4}{sup 3−}) from PHB/MMt/fertilizer at 90/0/10, 50/0/50, 80/10/10, 40/10/50, 25/25/50 formulation. The amount and type of the fertilizer and the presence of clay delineated the release of the cations and anions from the matrix. Cations released faster than anions through MMt composites due to cations–MMt interactions and its diffusion through MMt phase. The nutrients delivery from composites was lower than pure fertilizer. Composites formulations with smaller amount of PHB and MMt (25 wt%) were enough to encapsulate the fertilizer and delay the ions releasing. Further, the simple method of preparation also attains the sustainability of the whole process.

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

Bio-inspired redox mediated electrolyte for high performance flexible supercapacitor applications over broad temperature domain

For practical supercapacitors, recent advances in electrolytes have diminished long-standing problems such chemical stability, low ionic conductivity and limited electrochemical performance. Herewith, we communicate a facile method to produce anhydrous redox-active Bio-inspired electrolytes for flexible supercapacitors with application over a broad temperature range. The glycerol (Gly) gel electrolyte consists of various fractions of potassium hydroxide (KOH) and redox active unit, i.e., ammonium molybdate, (NH 4 ) 2 MoO 4 (abbreviated as Mo), in the anhydrous state. The doped glycerol matrix acts as a framework that provides a stable gel electrolyte, forming hierarchical pathway for ion diffusion. The developed gel network with Gly-KOH-MoX (X = 3, 5, 10, and 15%, w/w) provides an effective route to carry the ions onto the pores of the electrodes. The supercapacitor with Gly-KOH5-Mo10 gel provides a specific capacitance of 328 F g –1 , and the same device achieved an energy density of 45.6 Wh kg –1 at a power density of 497 W kg –1 . In addition, the supercapacitor maintained extraordinary cyclic robustness and wide temperature tolerance. Thus, the low cost and environmentally benign, cyclic durability of Gly-KOH-Mo makes these gel electrolytes a potential candidate for use in various energy storage devices.

25 ENERGY STORAGE↗
DOE OSTI2024

Correlations of conserved quantities at finite baryon density

Correlations involving the seven conserved quantities, namely energy, baryon number, electric charge, strangeness, and the three components of momentum, give rise to correlations in heavy-ion collisions. Through the utilization of a simple one-dimensional hydrodynamic model, we calculate the evolution of the entire 7×7 matrix of correlations as a function of relative spatial rapidity. This comprehensive analysis accounts for finite baryon density, which results in off-diagonal correlations between the charge-related quantities and the energy-momentum quantities. These correlations in coordinate space are subsequently transformed into correlations in momentum space using statistical weighting. As a result, the entire matrix of correlations is revealed to be highly sensitive to the equation of state, viscosity, and diffusivity.

Equations of state of nuclear matter↗
DOE OSTI2017

FOKKER-PLANCK ANALYSIS OF TRANSVERSE COLLECTIVE INSTABILITIES IN ELECTRON STORAGE RINGS

We analyze single bunch transverse instabilities due to wakefields using a Fokker-Planck model. We expand on the work of Suzuki [1], writing out the linear matrix equation including chromaticity, both dipolar and quadrupolar transverse wakefields, and the effects of damping and diffusion due to the synchrotron radiation. The eigenvalues and eigenvectors determine the collective stability of the beam, and we show that the predicted threshold current for transverse instability and the profile of the unstable agree well with tracking simulations. In particular, we find that predicting collective stability for high energy electron beams at moderate to large values of chromaticity requires the full Fokker-Planck analysis to properly account for the effects of damping and diffusion due to synchrotron radiation.

FOKKER-PLANCK ANALYSIS↗
DOE OSTI2019

Parfume (particle Fuel Model)

PARFUME is an advanced gas-cooled reactor fuel performance modeling and analysis code. It has been developed as an integrated mechanistic code that evaluates the thermal, mechanical, and physio-chemical behavior of fuel particles during irradiation and postulated accident conditions to determine the failure probability of a population of fuel particles given the particle-to-particle statistical variations in physical dimensions and material properties that arise from the fuel fabrication process, accounting for all viable mechanisms that can lead to particle failure. The code also determines the diffusion of fission products from the fuel through the particle coating layers, and through the fuel matrix to the coolant boundary. The subsequent release of fission products is calculated at the fuel element level (release of fission products from the compact or pebble). PARFUME calculates the release fraction as the ratio of the number of atoms released from the fuel element to the amount.

Gougar, Hans D.↗
DOE OSTI2018

Investigation on Explosive Welding of Zr{sub 53}Cu{sub 35}Al{sub 12} Bulk Metallic Glass with Crystalline Copper

A Zr{sub 53}Cu{sub 35}Al{sub 12} bulk metallic glass (BMG) was welded to a crystalline Cu using explosive welding technique. The morphology and the composition of the composite were characterized using optical microscopy, scanning electron microscopy, energy-dispersive x-ray spectroscopy and transmission electron microscopy. The investigation indicated that the BMG and Cu were tightly joined together without visible defects, and a thin diffusion layer appeared at the interface. The captured jet at the end of the welding region mostly comes from the Cu side. Amorphous and partially crystallized structures have been observed within the diffusion layer, but the BMG in close proximity to the interface still retains its amorphous state. Nanoindentation tests reveal that the interface exhibits an increment in hardness compared with the matrix on both sides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Analysis of Silver Release from Furnace-Tested TRISO Particles

High temperature testing of intact TRISO particles previously irradiated in the AGR-5/6/7 experiment was performed in the Furnace for Irradiated TRISO Testing (FITT) to directly confirm diffusive release of silver and europium from intact TRISO particles. Testing was conducted from 1,100–1,600°C for exposure times up to 100 h to directly confirm silver through layer release below safety testing temperatures. The FITT analysis showed highest levels of silver release at 1,300–1,400°C which confirms athermal release behaviors previous observed in step-wise and transient safety tests. Additionally, release was non-uniform with some particles releasing a majority of their inventory while others did not appear to release silver under identical testing conditions, which was consistent with historic observations. An assessment of the effective silver diffusion coefficient in the SiC layer was conducted and indicated maximum values in the 1,300–1,400°C range. The magnitude of the calculated diffusion coefficients also exceeded currently accepted diffusion coefficients. The release behavior of europium from intact particles was also analyzed in FITT for at 1,450–1,550°C for 500 h to 750 h. Direct confirmation of europium release below safety testing temperatures was confirmed absent contributions from matrix release. The analysis indicated europium release follows a general Arrhenius behavior and suggests general uniform release and indicates irradiation conditions influence observed release response. Calculated diffusion kinetics agree well with historic experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗
DOE OSTI2017

Initial CEED API (CEED-MS10)

Key components of CEED software development involve fast finite element operator storage and evaluation, architecture optimizations, performant algorithms for all orders, global kernels for finite element operators, efficient use of the memory sub-system, optimal data locality and motion, enhanced scalability and parallelism, and fast tensor contractions. As ultimate goals of the project, CEED is exploring and identifying the best algorithms for the full range of discretizations and applying algorithmic and software development to support ECP applications’ needs. In this milestone, we initiated the development and implementation of a common CEED API consisting of two main components: low-level and high-level APIs. We accomplished the goal of having the low-level component ready for the release of CEED v1.0 in the next quarter, see CEED-MS13. In this report we also provide a highlight from CEED’s integration with the MARBL app and the corresponding Laghos miniapp. A main theme in this work has been the systematic transition of advanced algorithms, developed and tested in Laghos, to the large-scale multi-physics settings of the MARBL code. The low-level CEED API provides a set of FE kernels and components for writing new low-level kernels. Examples include: vector and sparse linear algebra, element matrix assembly over a batch of elements, partial assembly and action for efficient high-order operators like mass, diffusion, advection, etc. The main goal of the low-level API is to establish the basis for the high-level API. Also, identifying such low-level kernels and providing a reference implementation for them serves as the basis for collaboration with vendors and STs. Another goal for the low-level API is to make it easy for applications with their own discretization infrastructure, who do not want to switch directly to using the high-level CEED API, to still evaluate and use the core operations provided by the low-level API. One of the challenges with high-order methods is that a global sparse matrix is no longer a good representation of a high-order linear operator, both with respect to the FLOPs needed for its evaluation, as well as the memory transfer needed for a matvec. Thus, another important question we tried to address while developing the low-level API is: is it possible to abstract out in a natural and efficient way the different implementations and data structures that arise when dealing with various computational device types: CPUs, GPUs, etc. Fully answering this question is beyond the scope of this particular milestone, since that requires longer time as we develop and expand the devices we support. The future high-level API will be the frontend provided by CEED consisting of high-level abstractions such as meshes, finite element spaces and solutions, bilinear forms, solvers, etc. At this high-level, an important aspect of the API will be to provide a uniform abstraction layer for writing code in a portable manner. To achieve that, we plan to utilize the MPI programming model combined with a device abstraction layer (DAL) for finite element computations that we will develop. The DAL will serve as the bridge between the low- and the high-level APIs, allowing a variety of programming models to be used for the low-level kernel implementations. The software artifacts delivered as part of this milestone include an initial CEED API software implementation (libCEED) provided through the CEED website, http://ceed.exascaleproject.org and the CEED GitHub organization, http://github.com/ceed. In addition to details and results from the above R&D efforts, in this document we are also reporting on other project-wide activities performed in Q1 of FY18, including: the Nek5000 v17.0 MFEM v3.3.2 and PETSc v3.8 releases, collaborations with ECP/ST and SciDAC projects, the inaugural Nek5000 hackathon and various outreach activities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Sub-nanometer glass surface dynamics induced by illumination

Illumination is known to induce stress and morphology changes in opaque glasses. Amorphous silicon carbide (a-SiC) has a smaller bandgap than the crystal. Thus, we were able to excite with 532 nm light a 1 μm amorphous surface layer on a SiC crystal while recording time-lapse movies of glass surface dynamics by scanning tunneling microscopy (STM). Photoexcitation of the a-SiC surface layer through the transparent crystal avoids heating the STM tip. Up to 6 × 10{sup 4} s, long movies of surface dynamics with 40 s time resolution and sub-nanometer spatial resolution were obtained. Clusters of ca. 3-5 glass forming units diameter are seen to cooperatively hop between two states at the surface. Photoexcitation with green laser light recruits immobile clusters to hop, rather than increasing the rate at which already mobile clusters hop. No significant laser heating was observed. Thus, we favor an athermal mechanism whereby electronic excitation of a-SiC directly controls glassy surface dynamics. This mechanism is supported by an exciton migration-relaxation-thermal diffusion model. Individual clusters take ∼1 h to populate states differently after the light intensity has changed. We believe the surrounding matrix rearranges slowly when it is stressed by a change in laser intensity, and clusters serve as a diagnostic. Such cluster hopping and matrix rearrangement could underlie the microscopic mechanism of photoinduced aging of opaque glasses.

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

SU-E-T-353: Effects of Time and Temperature On a Potential Reusable 3D Dosimeter

Purpose: Preliminary studies of a novel, optically-clearing PRESAGE 3D dosimeter formulation (Presage-RU) demonstrated potential reusability. This study investigates the effects of time and temperature on the accuracy and reusability of Presage-RU, and reports on progress toward developing a reusable 3D dosimeter. Methods: Presage-RU was cast as small volume samples (1×1×4.5cm). The effect of dose response sensitivity with reirradiation and time was evaluated by irradiating samples from 0–10Gy, measuring change in optical density (ΔOD), clearing at room temperature (RT) (5–7 days to fully clear), and then repeating for a total of 5 irradiations. Effects of heating on clearing rate were investigated by irradiating samples to 8Gy, then tracking measurements with samples held at RT, 35°C, and 45°C. Two cylindrical dosimeters (11cm diameter, 9.5cm length) were evaluated for dosimetric accuracy when stored at RT and −3°C prior to irradiation. Plans delivered were 2 overlapping AP fields (RT) and VMAT (-3°C). Results: Heating the dosimeters reduced the clearing half-life from 16.3h at RT to 5.8h (35°C) and 5.1h (45°C), but also increased background ΔOD by 1.7x (35°C) and 2.3x (45°C). Reductions in dose response were more closely linked to age than reirradiation, and storage at RT showed pronounced desensitization from dosimeter edges. These results suggest desensitization from oxygen diffusion. It should be noted that atmospheric diffusion into the dosimeter is not seen in standard, single-use PRESAGE, and is likely caused by differences in the Presage-RU polyurethane matrix. The dosimeter kept in cold storage, however, showed no evidence of desensitization and exhibited accuracy on par with standard PRESAGE with a 3%/3mm 3D gamma passing rate of 98.1%. Conclusions: Presage-RU is sensitive to storage temperatures and time, both of which affect oxygen diffusion and subsequent desensitization. Development shows promising progress with further formulation optimization as the next step toward achieving a successful reusable 3D dosimeter. This work was supported by NIH R01CA100835. John Adamovics is the president of Heuris Inc., which commercializes PRESAGE.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Solid-state reaction mediated microstructural evolution in a spark plasma sintered in situ Ti–B 4 C composite

For this study, a novel porosity-free in situ Ti–B 4 C composite was fabricated via spark plasma sintering (SPS). Detailed analysis of the phase evolution, coupling results from XRM, XRD and SEM-EBSD-EDS, indicate that primarily TiB 2 precipitates formed due to the in-situ reactions between the boron-carbide and titanium powders. The precipitation of TiB 2 resulted in the formation of a graphitic C-rich thin layer circumscribing the partially reacted B 4 C particles. Further progression of the in-situ reaction leads to the out-diffusion of the excess carbon from the B 4 C particles, across the graphite and TiB 2 phases, forming TiC upon reacting with the titanium matrix. Therefore, the final microstructure primarily consisted of TiB 2 , TiC, and partially reacted B 4 C phases, with small amounts of TiB, α-Ti, and graphite. Furthermore, the microstructure in these SPS processed in situ composites appeared to be the product of a solid-state Ti–B 4 C diffusion couple, substantially different from their fusion-based additively manufactured counterparts. Nano-indentation tests revealed a remarkably high average hardness of ~25 GPa for this SPS-processed Ti–B 4 C composite and comparable (with literature) phase-specific hardness and modulus values for the constituent TiB 2 , TiC, and B 4 C phases.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Thermal Conductivity and Diffusivity for SiC Fibers for Use in ATF Cladding Composites

As part of the accident tolerant fuels (ATF) initiative, there is a need to understand the thermophysical properties of cladding for these fuels. One potential cladding material is silicon carbide (SiC) ceramic matrix composite reinforced with SiC fibers because of the corrosion resistance, resistance to hydrogen formation, high bulk thermal conductivity, and favorable mechanical properties of SiC. However, data about the thermal properties of individual SiC fibers in these composites is absent in the literature. This paper presents the measured axial thermal conductivity and thermal diffusivity of a leading SiC fiber (HiNicalon Type S) via a thermal measurement technique (transient electrothermal technique, TET) optimized for thin fibers. The results show that the measured properties are significantly lower than bulk values for SiC. Because of the low neutron absorption, as well as excellent thermal and mechanical properties, zircaloy has primarily been used for the fuel cladding in many nuclear reactors. A major drawback of the zircaloy is its hydrogen generation vulnerability. This issue has been known for some time, but the explosion of hydrogen accident in the Fukushima Daiichi nuclear disaster has motivated efforts for accident tolerant fuels (ATF) and their cladding as suitable substitutes for zircaloy. Fuel cladding made of silicon carbide (SiC) ceramic matrix composite (CMC) is an attractive alternative to the zircaloy. Investigation on the tubular monolithic a-SiC polytropes indicates that it has superior oxidation (hydrogen production) resistance (three orders of magnitude slower) than zircaloy. Additionally, it has well-documented, favorable mechanical properties, along with high temperature corrosion resistance. However, its bulk thermal properties can vary significantly depending on polytrope, the grain size, grain boundary, or impurities. At room temperature, SiC has a thermal conductivity as high as 490 Wm{sup -1}K{sup -1}, where the {sup 6}H polytrope has the highest conductivity. Thermal measurements of the SiC CMC are usually accomplished by laser flash method for thermal diffusivity (a) and thermal conductivity is indirectly calculated through the relationship k=aρcp. Depending on the manufacturing method and the reinforced fiber, the composite thermal conductivity varies from <10 Wm{sup -1}K{sup -1} to ∼70 Wm{sup -1}K{sup -1}. While the flash method is convenient for the measurement of the composite, it is difficult to apply to the measurement of individual fibers. Because of this, the thermal properties of individual SiC fibers is not present in the literature. Because the SiC fiber manufacturing process is different than that of the SiC matrix in the CMC, the properties of the bulk material might be erroneous when applied to the fiber. Being able to measure the fiber properties independently and directly allows for an understanding of the effect of the manufacturing process on the fiber and the opportunity to potentially improve or control the thermal properties of the composite. This improved accuracy can be achieved by measuring the fiber with a technique optimized for individual fiber measurement, known as the Transient Electrothermal Technique or TET. This paper presents the thermal conductivity and diffusivity of a leading SiC fiber (HiNicalon Type S). Measurements are performed with the TET, with the resulting properties being significantly lower than the bulk values of SiC, which is one possible explanation for relatively low CMC values. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Quantum tunneling splittings from path-integral molecular dynamics

We illustrate how path-integral molecular dynamics can be used to calculate ground-state tunnelling splittings in molecules or clusters. The method obtains the splittings from ratios of density matrix elements between the degenerate wells connected by the tunnelling. We propose a simple thermodynamic integration scheme for evaluating these elements. Numerical tests on fully dimensional malonaldehyde yield tunnelling splittings in good overall agreement with the results of diffusion Monte Carlo calculations.

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

Liquid phase epitaxy of binary III–V nanocrystals in thin Si layers triggered by ion implantation and flash lamp annealing

The integration of III–V compound semiconductors in Si is a crucial step towards faster and smaller devices in future technologies. In this work, we investigate the formation process of III–V compound semiconductor nanocrystals, namely, GaAs, GaSb, and InP, by ion implantation and sub-second flash lamp annealing in a SiO{sub 2}/Si/SiO{sub 2} layer stack on Si grown by plasma-enhanced chemical vapor deposition. Raman spectroscopy, Rutherford Backscattering spectrometry, and transmission electron microscopy were performed to identify the structural and optical properties of these structures. Raman spectra of the nanocomposites show typical phonon modes of the compound semiconductors. The formation process of the III–V compounds is found to be based on liquid phase epitaxy, and the model is extended to the case of an amorphous matrix without an epitaxial template from a Si substrate. It is shown that the particular segregation and diffusion coefficients of the implanted group-III and group-V ions in molten Si significantly determine the final appearance of the nanostructure and thus their suitability for potential applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Dynamics in the Plastic Crystalline Phases of Cyclohexanol and Cyclooctanol Studied by Quasielastic Neutron Scattering

Plastic crystals are a promising candidate for solid state ionic conductors. In this work, quasielastic neutron scattering is employed to investigate the center of mass diffusive motions in two types of plastic crystalline cyclic alcohols: cyclohexanol and cyclooctanol. Two separate motions are observed which are attributed to long-range translational diffusion (α-process) and cage rattling (fast β-process). Residence times and diffusion coefficients are calculated for both processes, along with the confinement distances for the cage rattling. In addition, a binary mixture of these two materials is measured to understand how the dynamics change when a second type of molecule is added to the matrix. Furthermore, it is observed that, upon the addition of the larger cyclooctanol molecules into the cyclohexanol solution, the cage size decreases, which causes a decrease in the observed diffusion rates for both the α- and fast β-processes.

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

A computationally-efficient method for flamelet calculations

A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.

42 ENGINEERING↗