SEARCH · Engineering Papers

Results for “Sensitivity Matrix”

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.

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Interpretation of satellite-measured bidirectional reflectance from Cirrus cloudy atmospheres

The interpretation of the observed bidirectional reflectance from cirrus cloudy atmospheres is presented. A theoretical model was developed for the computation of the transfer of solar radiation in an anisotropic medium with particular applications to oriented ice crystals in cirrus clouds. In this model, the adding principle for radiative transfer was used with modifications to account for the anisotropy of scattering particles and the associated scattering phase matrix. The single-scattering properties, including the phase function, single-scattering albedo, and extinction cross sections, were used for randomly and horizontally oriented hexagonal ice crystals in radiative transfer computations. The radiative transfer model developed for the cirrus clouds was modified to account for the scattering contributions from the atmosphere and the surface. In order to test the relevance and significance of the ice crystal model for the interpretation of observed bidirectional reflectance from satellites, visible radiances collected on the half hour by the GOES series were selected. The data are calibrated and corrected with the proper filter functions, then navigated to match selected landmark data. A number of clear and cloudy cases during the cirrus IFO of the Fire experiment were chosen for theoretical analyses. The cloud particle shape and size distributions that were taken during satellite overpasses are used in radiative transfer calculations. The sensitivities of the shape, orientation, and size distribution of ice crystals on the reflected intensities at the top of the atmosphere are investigated. Finally, the relative importance of these cloud microphysical properties in the interpretation of satellite bidirectional reflectance are assessed and presented.

Takano, Y.↗

NASA NTRS1999

The Density Matrix of H20 - N2 In the Coordinate Representation: A Monte Carlo Calculation of the Far-Wing Line Shape

The far-wing line shape theory within the binary collision and quasistatic framework has been developed using the coordinate representation. Within this formalism, the main computational task is the evaluation of multidimensional integrals whose variables are the orientational angles needed to specify the initial and final positions of the system during transition processes. Using standard methods, one is able to evaluate the 7-dimensional integrations required for linear molecular systems, or the 7-dimensional integrations for more complicated asymmetric-top (or symmetric-top) molecular systems whose interaction potential contains cyclic coordinates. In order to obviate this latter restriction on the form of the interaction potential, a Monte Carlo method is used to evaluate the 9-dimensional integrations required for systems consisting of one asymmetric-top (or symmetric-top) and one linear molecule, such as H20-N2. Combined with techniques developed previously to deal with sophisticated potential models, one is able to implement realistic potentials for these systems and derive accurate, converged results for the far-wing line shapes and the corresponding absorption coefficients. Conversely, comparison of the far-wing absorption with experimental data can serve as a sensitive diagnostic tool in order to obtain detailed information on the short-range anisotropic dependence of interaction potentials.

Ma, Q.↗

NASA NTRS2008

Characterization of a New Phase and Its Effect on the Work Characteristics of a Near-Stoichiometric Ni30Pt20Ti50 High-Temperature Shape Memory Alloy (HTSMA)

A new phase observed in a nominal Ni30Pt20Ti50 (at.%) high temperature shape memory alloy has been characterized using transmission electron microscopy and 3-D atom probe tomography. This phase forms homogeneously in the B2 austenite matrix by a nucleation and growth mechanism and results in a concomitant increase in the martensitic transformation temperature of the base alloy. Although the structure of this phase typically contains a high density of faults making characterization difficult, it appears to be trigonal (-3m point group) with a(sub o) approx. 1.28 nm and c(sub o) approx. 1.4 nm. Precipitation of this phase increases the microhardness of the alloy substantially over that of the solution treated and quenched single-phase material. The effect of precipitation strengthening on the work characteristics of the alloy has been explored through load-biased strain-temperature testing in the solution-treated condition and after aging at 500 C for times ranging from 1 to 256 hours. Work output was found to increase in the aged alloy as a result of an increase in transformation strain, but was not very sensitive to aging time. The amount of permanent deformation that occurred during thermal cycling under load was small but increased with increasing aging time and stress. Nevertheless, the dimensional stability of the alloy at short aging times (1-4 hours) was still very good making it a potentially useful material for high-temperature actuator applications.

Garg, A.↗

NASA NTRS1992

Closed-form analysis of fiber-matrix interface stresses under thermo-mechanical loadings

Closed form techniques for calculating fiber matrix (FM) interface stresses, using repeating square and diamond regular arrays, were presented for a unidirectional composite under thermo-mechanical loadings. An Airy's stress function micromechanics approach from the literature, developed for calculating overall composite moduli, was extended in the present study to compute FM interface stresses for a unidirectional graphite/epoxy (AS4/3501-6) composite under thermal, longitudinal, transverse, transverse shear, and longitudinal shear loadings. Comparison with finite element results indicate excellent agreement of the FM interface stresses for the square array. Under thermal and longitudinal loading, the square array has the same FM peak stresses as the diamond array. The square array predicted higher stress concentrations under transverse normal and longitudinal shear loadings than the diamond array. Under transverse shear loading, the square array had a higher stress concentration while the diamond array had a higher radial stress concentration. Stress concentration factors under transverse shear and longitudinal shear loadings were very sensitive to fiber volume fraction. The present analysis provides a simple way to calculate accurate FM interface stresses for both the square and diamond array configurations.

Naik, Rajiv A.↗

NASA NTRS1996

Processing and Properties of Fiber Reinforced Polymeric Matrix Composites: Processing Robustness of IM7/PETI Polyimide Composites - Part 2

The processability of a phenylethynyl terminated imide (PETI) resin matrix composite was investigated. Unidirectional prepregs were made by coating an N-methylpyrrolidone solution of the amide acid oligomer onto unsized IM7. Two batches of prepregs were used: one was made by NASA in-house, and the other was from an industrial source. The composite processing robustness was investigated with respect to the effect of B-staging conditions, the prepreg shelf life, and the optimal processing window. Rheological measurements indicated that PETI's processability was only slightly affected over a wide range of B-staging temperatures (from 250 C to 300 C). The open hole compression (OHC) strength values were statistically indistinguishable among specimens consolidated using various B-staging conditions. Prepreg rheology and OHC strengths were also found not to be affected by prolonged (i.e., up to 60 days) ambient storage. An optimal processing window was established using response surface methodology. It was found that IM7/PETI composite is more sensitive to the consolidation temperature than to the consolidation pressure. A good consolidation was achievable at 371 C/100 Psi, which yielded an OHC strength of 62 Ksi at room temperature. However, processability declined dramatically at temperatures below 350 C.

Hou, Tan-Hung↗

DOE OSTI2021

Prompt Gamma Analysis of ARIES Materials and Updates to the 2015 Calibration Equations

Prompt gamma (PG) analysis is a nondestructive, nuclear, elemental analysis technique that uses charged particle reactions to interrogate a sample, and elements present in the sample matrix are identified through the characteristic gamma-rays emitted from the product nuclei in alpha-p and alpha-n nuclear reactions. This technique has been applied to plutonium oxide packaged in over 4,000 individual 3013 containers, and the concentrations of certain light elements were determined from the integrated peak areas based on a calibration that was published previously. This report provides the results for a new population of 3013 containers packaged with oxide materials produced by the conversion of metal by Advanced Recovery and Integrated Extraction System (ARIES) project using the direct metal oxidation (DMO) process and muffle furnaces. New PG and analytical chemistry data collected since 2015 were added to the existing calibration data sets to refine the calibration parameters. Calibration equations were also developed for determining beryllium and fluorine at low concentrations in high-purity ARIES product oxides. The new fluorine calibration provides an order of magnitude greater sensitivity and results in an additional 1,408 containers in the original population identified as having fluorine as an impurity. Additionally, equations for calculating the lower limits of detection (LLDs) as a function of the actual counting time (live time) were obtained using WLS regression technique for samples in the calibration data set. This resulted in changes to the LLDs published previously.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2010

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗

DOE OSTI2021

Modulation of Akt vs Stat3 activity by the focal adhesion kinase in non-neoplastic mouse fibroblasts

Adhesion of cells to each other and to the extracellular matrix (ECM) are both required for cellular functions. Cell-to-cell adhesion is mediated by cadherins and their engagement triggers the activation of Stat3, which offers a potent survival signal. Adhesion to the ECM on the other hand, activates FAK which attracts and activates Src, as well as receptor tyrosine kinases (RTKs), the PI3k/Akt and Ras/Erk pathways. However, the effect of cell density upon FAK and Akt activity has not been examined. We now demonstrate that, interestingly, despite being potent Stat3 activators, Src and RTKs are unable to activate Stat3 in sparsely growing (i.e., without cadherin engagement), non-neoplastic cells attached to the ECM. In contrast, cell aggregation (i.e., cadherin engagement in the absence of adhesion to a solid substratum) was found to activate both Stat3 and Akt. Pharmacologic or genetic reduction of FAK activity abolished Akt activity at low densities, indicating that FAK is an important activator of Akt in this setting. Notably, FAK knockout increased cellular sensitivity to the Stat3 inhibitor CPA7, while FAK reintroduction restored resistance to this drug. These findings suggest a complementary role of integrin/FAK/Akt and cadherin/Stat3-mediated pro-survival pathways, which may be of significance during neoplastic transformation and metastasis.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Modelling nearest neighbour distributions of biased tracers using hybrid effective field theory

ABSTRACT We investigate the application of hybrid effective field theory (HEFT) – which combines a Lagrangian bias expansion with subsequent particle dynamics from N-body simulations – to the modelling of k-nearest neighbour cumulative distribution functions (kNN-CDFs) of biased tracers of the cosmological matter field. The kNN-CDFs are sensitive to all higher order connected N-point functions in the data, but are computationally cheap to compute. We develop the formalism to predict the kNN-CDFs of discrete tracers of a continuous field from the statistics of the continuous field itself. Using this formalism, we demonstrate how kNN-CDF statistics of a set of biased tracers, such as haloes or galaxies, of the cosmological matter field can be modelled given a set of low-redshift HEFT component fields and bias parameter values. These are the same ingredients needed to predict the two-point clustering. For a specific sample of haloes, we show that both the two-point clustering and the kNN-CDFs can be well-fit on quasi-linear scales (≳ 20h−1Mpc) by the second-order HEFT formalism with the same values of the bias parameters, implying that joint modelling of the two is possible. Finally, using a Fisher matrix analysis, we show that including kNN-CDF measurements over the range of allowed scales in the HEFT framework can improve the constraints on σ8 by roughly a factor of 3, compared to the case where only two-point measurements are considered. Combining the statistical power of kNN measurements with the modelling power of HEFT, therefore, represents an exciting prospect for extracting greater information from small-scale cosmological clustering.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Search for Heavy Neutral Leptons at the MINER$\nu$A detector

Heavy Neutral Leptons (HNL) are particles hypothesised to provide a mass generation mechanism for the active (observed) neutrino species, which are known to have nonzero mass from the definitive observation of neutrino oscillations. HNL are eigenstates of mass of the order $\mathcal{O}(0.1 − 1\,\,\mathrm{GeV}/c^{2})$, which mix into the active flavour eigenstates through the extended leptonic mixing matrix. Apart from neutrino mass, they could provide a natural dark matter candidate and a mechanism for matter-antimatter asymmetry in the early Universe, giving rise to today’s matter-dominated cosmos. Searches for HNL typically attempt to either confirm the existence of HNL through an excess in data that is most compatible with an HNL hypothesis, or by setting limits on the HNL parameter space $\left(M_{N4}, \left|U_{\alpha 4}\right|^{2}\right)$ in the case of statistically insignificant excess. Such searches are intensifying as part of a global research programme at both colliders and accelerator / atmospheric neutrino experiments. One such setting is the MINER$\nu$A experiment, located in the NuMI beamline at Fermilab. MINER$\nu$A has collected a large amount of data over seven years of operation to measure the cross-sections of neutrino-nucleus interactions necessary to drive systematic uncertainties down, in order for neutrino oscillation experiments to achieve sensitivity to CP violation in the neutrino sector. With the high-energy, high-intensity NuMI beam, and with good timing, position, and angular resolution leading to sensitivity to HNL decays, MINER$\nu$A is fertile ground for an HNL search, which is performed in this thesis. Specifically, a novel, experiment-agnostic and general simulation of HNL production and decay is presented and deployed; this simulation has been incorporated in the ubiquitous GENIE neutrino event generator for use with neutrino experiments in the future. The event selection and background characterisation is discussed in detail, including the primary background coming from charged-current coherent and diffractive pion production from neutrino-nucleus interactions, and constraints on the background using control regions from MINER$\nu$A data are derived. Finally, by means of a fake-data study given a background of $\mathcal{O}(500)$ events, the discovery potential and limit-setting capacity of MINER$\nu$A is demonstrated, and ways to improve this capacity are expounded upon.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Environmental matrix and moisture are key determinants of microbial phenotypes expressed in a reduced complexity soil-analog

Soil moisture and porosity regulate microbial metabolism by influencing factors such as redox conditions, substrate availability, and soil connectivity. However, the inherent biological, chemical, and physical heterogeneity of soil complicates laboratory investigations into microbial phenotypes that mediate community metabolism. This difficulty arises from challenges in accurately representing the soil environment and in establishing a tractable microbial community that limits confounding variables. To address these challenges in our investigation of community metabolism, we use a reduced-complexity microbial consortium grown in a soil analog using a glass-bead matrix amended with chitin. Long-read and short-read metagenomes, metatranscriptomes, metaproteomes, and metabolomes were analyzed to test the effects of soil structure and moisture on chitin degradation. Our soil structure analog system greatly altered microbial expression profiles compared to the liquid-only incubations, emphasizing the importance of incorporating environmental parameters, like pores and surfaces, for understanding microbial phenotypes relevant to soil ecosystems. These changes were mainly driven by differences in overall expression of chitin-degrading Streptomyces species and stress-tolerant Ensifer. Our findings suggest that the success of Ensifer in a structured environment is likely related to its ability to repurpose carbon via the glyoxylate shunt while potentially using polyhydroxyalkanoate granules as a C source. We also identified traits like motility, stress resistance, and biofilm formation that underlie the degradation of chitin across our treatments and inform how they may ultimately alter carbon use efficiency. Together our results demonstrate that community functions like decomposition are sensitive to environmental conditions and more complex than the multi-enzyme pathways involved in depolymerization.

Rodriguez-Ramos, Josue A↗

NASA NTRS2002

Development of Three-Dimensional Multicellular Tissue-Like Constructs for Mutational Analysis Using Macroporous Microcarriers

A three-dimensional (3-D), tissue-like model was developed for the genotoxic assessment of space environment. In previous experiments, we found that culturing mammalian cells in a NASA-designed bioreactor, using Cytodex-3 beads as a scaffold, generated 3-D multicellular spheroids. In an effort to generate scaffold-free spheroids, we developed a new 3-D tissue-like model by coculturing fibroblast and epithelial cell in a NASA bioreactor using macroporous Cultispher-S(TradeMark) microcarriers. Big Blue(Registered Trademark) Rat 2(Lambda) fibroblasts, genetically engineered to contain multiple copies (>60 copies/cell) of the Lac I target gene, were cocultured with radio-sensitive human epithelial cells, H184F5. Over an 8-day period, samples were periodically examined by microscopy and histology to confirm cell attachment, growth, and viability. Immunohistochemistry and western analysis were used to evaluate the expression of specific cytoskeletal and adhesion proteins. Key cell culture parameters (glucose, pH, and lactate concentrations) were monitored daily. Controls were two-dimensional mono layers of fibroblast or epithelial cells cultured in T-flasks. Analysis of 3-D spheroids from the bioreactor suggests fibroblast cells attached to and completely covered the bead surface and inner channels by day 3 in the bioreactor. Treatment of the 3-day spheroids with dispase II dissolved the Cultisphers(TradeMark) and produced multicellular, bead-less constructs. Immunohistochemistry confirmed the presence of vi.mentin, cytokeratin and E-cadherin in treated spheroids. Examination of the dispase II treated spheroids with transmission electron microscopy (TEM) also showed the presence of desmosomes. These results suggest that the controlled enzymatic degradation of an artificial matrix in the low shear environment of the NASA-designed bioreactor can produce 3-D tissue-like spheroids. 2

Jordan, Jacqueline A.↗

DOE OSTI2022

Exploring the Use of Novel Spatial Accelerators in Scientific Applications

Driven by the need to find alternative accelerators which can viably replace GPUs in next-generation Supercomputing systems, this paper proposes a methodology to enable agile application/hardware co-design. The application-first methodology provides the ability to come up with design of accelerators while working with real-world workloads, available accelerators, and system software. The iterative design process targets a set of kernels in a workload for performance estimates that can prune the design space for later phases of detailed architectural evaluations. To this effect, in this paper, a novel data-parallel device model is introduced that simulates the latency of performance-sensitive operations in an accelerator including data transfers and kernel computation using multi-core CPUs. The use of off-the-shelf simulators, such as pre-RTL simulator Aladdin or multiple tools available for exploring the design of deep neural network accelerators (e.g., Timeloop) is demonstrated for evaluation of various accelerator designs using applications with realistic inputs. Examples of multiple device configurations that are instantiable in a system are explored to evaluate the performance benefit of deploying novel accelerators. The proposed device is integrated with a programming model and system software to potentially explore the impacts of high-level programming languages/compilers and low-level effects such as task scheduling on multiple accelerators. We analyze our methodology for a set of applications that represent high-performance computing (HPC) and graph analytics. The applications include a computational chemistry kernel realized using tensor contractions, triangle counting, GraphSAGE and Breadth-first Search. These applications include kernels such as dense matrix-dense matrix multiplication, sparse matrix-spare matrix multiplication, and sparse matrix-dense vector multiplication. Our results indicate potential performance benefits and insights for system design by including accelerators that realize these kernels along-side general purpose accelerators.

AI, codesign, Accelerated Computing, Modeling and ↗

DOE OSTI2022

Fundamental Investigations into the Properties and Performance of Advanced Materials

Intelligent engineering designs require an accurate understanding of material behavior, since any uncertainties or gaps in knowledge must be counterbalanced with heightened factors of safety, leading to overdesign. Therefore, building better structures and pushing the performance of new components requires an improved understanding of the thermomechanical response of advanced materials under service conditions. This dissertation provides fundamental investigations of several advanced materials: thermoset polymers, a common matrix material for fiber-reinforced composites and nanocomposites; aluminum alloy 7075-T6 (AA7075-T6), a high-performance aerospace material; and ceramic matrix composites (CMCs), an advanced composite for extreme-temperature applications. To understand matrix interactions with various interfaces and nanoinclusions at their fundamental scale, the properties of thermoset polymers are studied at the atomistic scale. An improved proximity-based molecular dynamics (MD) technique for modeling the crosslinking of thermoset polymers is carefully established, enabling realistic curing simulations through its ability to dynamically and probabilistically perform complex topology transformations. The proximity-based MD curing methodology is then used to explore damage initiation and the local anisotropic evolution of mechanical properties in thermoset polymers under uniaxial tension with an emphasis on changes in stiffness through a series of tensile loading, unloading, and reloading experiments. Aluminum alloys in aerospace applications often require a fatigue life of over 109 cycles, which is well over ii the number of cycles that can be practically tested using conventional fatigue testing equipment. In order to study these high-life regimes, a detailed ultrasonic cycle fatigue study is presented for AA7075-T6 under fully reversed tension-compression loading. The geometric sensitivity, frequency effects, size effects, surface roughness effects, and the corresponding failure mechanisms for ultrasonic fatigue across different fatigue regimes are investigated. Finally, because CMCs are utilized in extreme environments, oxidation plays an important role in their degradation. A multiphysics modeling methodology is thus developed to address the complex coupling between oxidation, mechanical stress, and oxygen diffusion in heterogeneous carbon fiber-reinforced CMC microstructures.

Schichtel, Jacob J.↗

DOE OSTI2024

Globally optimal interferometry with lossy twin Fock probes

Parity or quadratic spin (e.g., J z 2 ) readouts of a Mach–Zehnder (MZ) interferometer probed with a twin Fock (TF) input state allow saturating the optimal sensitivity attainable among all mode-separable states with a fixed total number of particles but only when the interferometer phase θ is near zero. When more general Dicke state probes are used, the parity readout saturates the quantum Fisher information (QFI) at θ = 0, whereas better-than-standard quantum limit performance of the J z 2 readout is restricted to an o ( N ) occupation imbalance. We show that a method of moments readout of two quadratic spin observables J z 2 and J + 2 + J − 2 is globally optimal for Dicke state probes; i.e., the error saturates the QFI for all θ . In the lossy setting, we derive the time-inhomogeneous Markov process describing the effect of particle loss on TF states, showing that the method of moments readout of four at-most-quadratic spin observables is sufficient for globally optimal estimation of θ when two or more particles are lost. The analysis culminates in a numerical calculation of the QFI matrix for distributed MZ interferometry on the four-mode state | N 4 , N 4 , N 4 , N 4 〉 and its lossy counterparts, showing that an advantage for the estimation of any linear function of the local MZ phases θ 1 and θ 2 (compared to independent probing of the MZ phases by two copies of | N 4 , N 4 〉 ) appears when more than one particle is lost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2004

Processing Robustness for A Phenylethynyl Terminated Polyimide Composite

The processability of a phenylethynyl terminated imide resin matrix (designated as PETI-5) composite is investigated. Unidirectional prepregs are made by coating an N-methylpyrrolidone solution of the amide acid oligomer (designated as PETAA-5/NMP) onto unsized IM7 fibers. Two batches of prepregs are used: one is made by NASA in-house, and the other is from an industrial source. The composite processing robustness is investigated with respect to the prepreg shelf life, the effect of B-staging conditions, and the optimal processing window. Prepreg rheology and open hole compression (OHC) strengths are found not to be affected by prolonged (i.e., up to 60 days) ambient storage. Rheological measurements indicate that the PETAA-5/NMP processability is only slightly affected over a wide range of B-stage temperatures from 250 deg C to 300 deg C. The OHC strength values are statistically indistinguishable among laminates consolidated using various B-staging conditions. An optimal processing window is established by means of the response surface methodology. IM7/PETAA-5/NMP prepreg is more sensitive to consolidation temperature than to pressure. A good consolidation is achievable at 371 deg C (700 deg F)/100 Psi, which yields an RT OHC strength of 62 Ksi. However, processability declines dramatically at temperatures below 350 deg C (662 deg F), as evidenced by the OHC strength values. The processability of the IM7/LARC(TM) PETI-5 prepreg was found to be robust.

Hou, Tan-Hung↗

DOE OSTI2025

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗

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↗