Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural model decomposition”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Computational and Experimental Study of Energetic Materials in a Counterflow Microgravity Environment

Counterflow diffusion flames are studied for various fuels flowing against decomposition products from solid ammonium perchlorate (AP) pellets in order to obtain fundamental understanding of composite propellant flame structure and chemistry. We illustrate this approach through a combined experimental and numerical study of a fuel mixture consisting of C2H4 CO + H2, and C2H2 + C2H4 flowing against solid AP. For these particular AP-fuel systems, the resulting flame zone simulates the various flame structures that are ex+ to exist between reaction products from Ap crystals and a hydrocarbon binder. As in all our experimental studies, quantitative species and temperature profiles have been measured between the fuel exit and AP surface. Species measured included CN, NH, NO, OH, N2, CO2, CO, H2, CO, HCl, and H2O. Temperature was measured using a thermocouple at the exit, spontaneous Raman scattering measurements throughout the flame, OH rotational population distributions, and NO vibrational population distributions. The burning rate of AP was also measured as a function of strain rate, given by the separation distance between the AP surface and the gaseous hydrocarbon fuel tube exit plane. This distance was nominally set at 5 mm, although studies have been performed for variations in separation distance. The measured 12 scalars are compared with predictions from a detailed gas-phase kinetics model consisting of 86 species and 531 reactions. Model predictions are found to be in good agreement with experiment and illustrate the type of kinetic features that may be expected to occur in propellants when AP particle size distributions are varied. Furthermore, the results constitute the continued development of a necessary database and validation of a comprehensive model for studying more complex AP-solid fuel systems in microgravity. Exploratory studies have also been performed with liquid and solid fuels at normal gravity. Because of melting (and hence dripping) and deep thermal wave penetration into the liquid, these experiments were found feasible, but not used for obtaining quantitative data. Microgravity experiments are needed to eliminate the dripping and boiling phenomena of these systems at normal gravity. Microgravity tests in the NASA Glenn 2.2 second drop tower were performed (1) to demonstrate the feasibility of performing propellant experiments using the NASA Glenn microgravity facilities, (2) to develop the operational procedures for safe handing of the energetic materials and disposal of their toxic combustion by-products and (3) to obtain initial measurements of the AP burning rate and flame structure under microgravity conditions. Experiments were conducted on the CH4/AP system previously studied at normal gravity using a modified design of the counterflow burner and a NASA Glenn Pig Rig, i.e., one of the existing drop rigs for general-purpose usage. In these experiments, the AP burning rate was measured directly with a linear variable differential transducer (LVDT) and video imaging of the flame structure was recorded ignition was achieved by hot wires stretched across the AP surfaces. Initial drop tower combustion data show that with the same burner separation distance and flow conditions of the normal gravity experiments, the AP burning rate is approximately a factor of two lower. This difference is likely a result of radiation effects, but further tests with longer test times need to be conducted to verify that steady state conditions were achieved under microgravity conditions.

Takahashi, Fumiaki↗
DOE OSTI2013

Thermal Pretreatment of Wood for Co-gasification/co-firing of Biomass and Coal

Utilization of biomass as a co-feed in coal and biomass co-firing and co-gasification requires size reduction of the biomass. Reducing biomass to below 0.2 mm without pretreatment is difficult and costly because biomass is fibrous and compressible. Torrefaction is a promising thermal pretreatment process and has the advantages of increasing energy density, improving grindability, producing fuels with more homogenous compositions and hydrophobic behavior. Temperature is the most important factor for the torrefaction process. Biomass grindability is related to cell wall structure, thickness and composition. Thermal treatment such as torrefaction can cause chemical changes that significantly affect the strength of biomass. The objectives of this study are to understand the mechanism by which torrefaction improves the grindability of biomass and discuss suitable temperatures for thermal pretreatment for co-gasification/co-firing of biomass and coal. Wild cherry wood was selected as the model for this study. Samples were prepared by sawing a single tangential section from the heartwood and cutting it into eleven pieces. The samples were consecutively heated at 220, 260, 300, 350, 450 and 550⁰C for 0.5 hr under flowing nitrogen in a tube furnace. Untreated and treated samples were characterized for physical properties (color, dimensions and weight), microstructural changes by SEM, and cell wall composition changes and thermal behaviors by TGA and DSC. The morphology of the wood remained intact through the treatment range but the cell walls were thinner. Thermal treatments were observed to decompose the cell wall components. Hemicellulose decomposed over the range of ~200 to 300⁰C and resulted in weakening of the cell walls and subsequently improved grindability. Furthermore, wood samples treated above 300⁰C lost more than 39% in mass. Therefore, thermal pretreatment above the hemicelluloses decomposition temperature but below 300⁰C is probably sufficient to improve grindability and retain energy value.

09 BIOMASS FUELS↗
DOE OSTI2025

Autonomous direct freeform fabrication strategy for multi-axis additive manufacturing

Multi-axis additive manufacturing (M-AM) enables precise material deposition along both planar and curved layers, eliminating the need for support structures through a continuous material deposition approach. In contrast to conventional 2-dimensional planar layers constrained to a fixed building orientation, the deposition on freeform layers demands the specification of guided curves to determine material deposition directions which are no longer to be fixed to a build direction. There are challenges that arise when fabricating components with multiple “build” directions, necessitating the decomposition of geometries and the specification of guided curves for the resulting volumes. Furthermore, multi-axis systems introduce heightened challenges due to an increased degree of freedom in motion. Consequently, the potential risks of collision between the deposited geometry and the motion platform become a notable concern. This research proposes a freeform layering algorithm to address the challenge of seamless transitions between planar and curved layers in the process planning of M-AM. The algorithm computes 3D “printable” layers by leveraging topological information derived from the geometry to be built and integrates collision-free manufacturability considerations into the computational process. These accumulated volumes serve as a “substrate” and support volumes for subsequent deposition, allowing later layers to be built without the need for additional support material. In conclusion, the proposed method successfully devises a freeform layering approach suitable for intricate models that demand substantial support, thus enabling the fabrication of diverse geometries in a manner previously unachievable.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pre-conditioned BFGS-based uncertainty quantification in elastic full-waveform inversion

SUMMARY Full-waveform inversion has become an essential technique for mapping geophysical subsurface structures. However, proper uncertainty quantification is often lacking in current applications. In theory, uncertainty quantification is related to the inverse Hessian (or the posterior covariance matrix). Even for common geophysical inverse problems its calculation is beyond the computational and storage capacities of the largest high-performance computing systems. In this study, we amend the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm to perform uncertainty quantification for large-scale applications. For seismic inverse problems, the limited-memory BFGS (L-BFGS) method prevails as the most efficient quasi-Newton method. We aim to augment it further to obtain an approximate inverse Hessian for uncertainty quantification in FWI. To facilitate retrieval of the inverse Hessian, we combine BFGS (essentially a full-history L-BFGS) with randomized singular value decomposition to determine a low-rank approximation of the inverse Hessian. Setting the rank number equal to the number of iterations makes this solution efficient and memory-affordable even for large-scale problems. Furthermore, based on the Gauss–Newton method, we formulate different initial, diagonal Hessian matrices as pre-conditioners for the inverse scheme and compare their performances in elastic FWI applications. We highlight our approach with the elastic Marmousi benchmark model, demonstrating the applicability of pre-conditioned BFGS for large-scale FWI and uncertainty quantification.

58 GEOSCIENCES↗
DOE OSTI2019

Operando FTIR Studies of Electrolyte Solution Structures and Dynamics at Electrode-Electrolyte Interfaces

The well-stabilized electrode-electrolyte interface (EEI) is critical to battery performance and safety, thus understanding the electrolyte solution structures (e.g., ion solvation and ionic association) and dynamics within EEI is an important area of study. Unfortunately, due to the sensitive nature of the electrodes, ex-situ characterization methods are limited and challenging to interpret as chemical and structural changes occur upon deconstructing the cell. In this work, we present operando characterization of near-surface electrolyte solution structures and dynamics by ATR-FTIR spectroscopy during battery cycling. We studied two selected model electrodes, a lithium nickel manganese cobalt oxide (NMC) or a planar silicon (Si), in a standard Gen2 electrolyte of 1.2 M LiPF6 in ethylene carbonate (EC):ethyl methyl carbonate (EMC) (3:7 wt%). In the former case, the composite cathode on aluminum foil is placed face-down on an ATR crystal. A ring of lithium metal serves as the counter/reference electrode and no separator is required. In the latter case, a lightly doped p-type silicon wafer is phosphorous-doped at the surface to increase the conductivity of the sample without compromising infrared transparency. The ATR crystal is then pressed on the back of the silicon with a minimal air gap to ensure sufficient optical contact between the two materials. The front of wafer is then placed in a modified coin-type cell geometry with a minimal electrolyte, a separator, and a lithium metal counter/reference electrode. In both experiments, we analyze the electrolyte solution structures within the EEI and find similar patterns as a function of potential. Firstly, we find that EMC decomposes faster than EC -suggesting that EMC plays a greater role in the electrochemical reactions on the electrode surface (e.g., formation of a protective layer to prevent further electrolyte decomposition). Secondly, we monitor the increase/decrease of the vibrational absorption peaks attributed to carbonates coordinated with Li+ ions relative to the decrease/increase of their uncoordinated counterparts, as expected when changing the local concentration of Li+ ions near the electrode surface during deintercalation/intercalation into the electrode. Finally, we observe that the features assigned to the Li+-PF6 -/-solvent molecules contact ion pair relative to the uncoordinated PF6 - anion. These observations give important information regarding the electrolyte solution structures and dynamics near the electrode surfaces. We look forward to applying these techniques to additional lithiumion battery systems with different cathode or anode materials and beyond lithium-ion systems to characterize surface chemistry/evolution (on-going research) and solution structures/dynamics (in particular, (de)solvation behaviors) within the EEI, which can be a bottleneck to efficient ion transfer from the electrolyte to the electrode.

14 SOLAR ENERGY↗
DOE OSTI2016

Noise suppression for dual-energy CT via penalized weighted least-square optimization with similarity-based regularization

Purpose: Dual-energy CT (DECT) expands applications of CT imaging in its capability to decompose CT images into material images. However, decomposition via direct matrix inversion leads to large noise amplification and limits quantitative use of DECT. Their group has previously developed a noise suppression algorithm via penalized weighted least-square optimization with edge-preservation regularization (PWLS-EPR). In this paper, the authors improve method performance using the same framework of penalized weighted least-square optimization but with similarity-based regularization (PWLS-SBR), which substantially enhances the quality of decomposed images by retaining a more uniform noise power spectrum (NPS). Methods: The design of PWLS-SBR is based on the fact that averaging pixels of similar materials gives a low-noise image. For each pixel, the authors calculate the similarity to other pixels in its neighborhood by comparing CT values. Using an empirical Gaussian model, the authors assign high/low similarity value to one neighboring pixel if its CT value is close/far to the CT value of the pixel of interest. These similarity values are organized in matrix form, such that multiplication of the similarity matrix to the image vector reduces image noise. The similarity matrices are calculated on both high- and low-energy CT images and averaged. In PWLS-SBR, the authors include a regularization term to minimize the L-2 norm of the difference between the images without and with noise suppression via similarity matrix multiplication. By using all pixel information of the initial CT images rather than just those lying on or near edges, PWLS-SBR is superior to the previously developed PWLS-EPR, as supported by comparison studies on phantoms and a head-and-neck patient. Results: On the line-pair slice of the Catphan{sup ©}600 phantom, PWLS-SBR outperforms PWLS-EPR and retains spatial resolution of 8 lp/cm, comparable to the original CT images, even at 90% reduction in noise standard deviation (STD). Similar performance on spatial resolution is observed on an anthropomorphic head phantom. In addition, results of PWLS-SBR show substantially improved image quality due to preservation of image NPS. On the Catphan{sup ©}600 phantom, NPS using PWLS-SBR has a correlation of 93% with that via direct matrix inversion, while the correlation drops to −52% for PWLS-EPR. Electron density measurement studies indicate high accuracy of PWLS-SBR. On seven different materials, the measured electron densities calculated from the decomposed material images using PWLS-SBR have a root-mean-square error (RMSE) of 1.20%, while the results of PWLS-EPR have a RMSE of 2.21%. In the study on a head-and-neck patient, PWLS-SBR is shown to reduce noise STD by a factor of 3 on material images with image qualities comparable to CT images, whereas fine structures are lost in the PWLS-EPR result. Additionally, PWLS-SBR better preserves low contrast on the tissue image. Conclusions: The authors propose improvements to the regularization term of an optimization framework which performs iterative image-domain decomposition for DECT with noise suppression. The regularization term avoids calculation of image gradient and is based on pixel similarity. The proposed method not only achieves a high decomposition accuracy, but also improves over the previous algorithm on NPS as well as spatial resolution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Raman and infrared spectra of plutonium (IV) oxalate and its thermal degradation products

For over 80 years, plutonium dioxide has been routinely produced via thermal decomposition of hydrated plutonium(IV) oxalate. Despite the longstanding utility of this process, the chemical structures of starting materials and intermediates produced during this thermal conversion remain ill-defined. To help resolve this uncertainty, we measured high-resolution Raman and infrared spectra of Pu(C 2 O 4 ) 2 ·6H 2 O that was heated to 25, 100, 220, 250, 350, and 450 °C in air. Our measurements show that Pu(C 2 O 4 ) 2 ·6H 2 O has a rich vibrational spectrum with at least 15 Raman bands between 180 cm -1 and 1900 cm -1 and 9 infrared bands between 800 cm -1 and 4000 cm -1 . As Pu(C 2 O 4 ) 2 ·6H 2 O is heated, water is liberated, and the oxalate ligand decomposes to produce plutonium oxycarbide species. When heated to 350 °C or higher, vibrational spectra are consistent with PuO 2 with some residual carbon-containing species. Full vibrational spectra, powder X-ray diffraction, and scanning electron microscopy measurements of Pu(C 2 O 4 ) 2 ·6H 2 O and its thermal degradation products are presented herein along with approximate assignments for observed spectral bands. These data can be used to validate and potentially improve existing computational models that describe the chemical structure of compounds produced during thermal degradation of plutonium (IV) oxalate. Given the utility of plutonium (IV) oxalate in synthesizing plutonium dioxide, these results are expected to provide value in the fields of nuclear fuel processing, nuclear nonproliferation, and nuclear forensics.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Rank-reduced coupled-cluster. III. Tensor hypercontraction of the doubles amplitudes

In this work, we develop a quartic-scaling implementation of coupled-cluster singles and doubles (CCSD) based on low-rank tensor hypercontraction (THC) factorizations of both the electron repulsion integrals (ERIs) and the doubles amplitudes. This extends our rank-reduced (RR) coupled-cluster method to incorporate higher-order tensor factorizations. The THC factorization of the doubles amplitudes accounts for most of the gain in computational efficiency as it is sufficient, in conjunction with a Cholesky decomposition of the ERIs, to reduce the computational complexity of most contributions to the CCSD amplitude equations. Further THC factorization of the ERIs reduces the complexity of certain terms arising from nested commutators between the doubles excitation operator and the two-electron operator. We implement this new algorithm using graphical processing units and demonstrate that it enables CCSD calculations for molecules with 250 atoms and 2500 basis functions using a single computer node. Furthermore, we show that the new method computes correlation energies with comparable accuracy to the underlying RR-CCSD method.

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

The Surface Chemistry of Methanol on Cu 3 Pd(111): Effects of Metal Alloying and Reaction with Hydrogen

Synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS) was used to study the adsorption and surface chemistry of methanol on a Cu 3 Pd(111) model surface. The composition and morphological properties of the pristine Cu 3 Pd(111) substrate were analyzed using a combination of low-energy electron diffraction (LEED), scanning tunneling microscopy (STM), and low-energy ion scattering (LEIS). The results of ion scattering showed segregation of Pd toward the surface, with a Pd/Cu ratio close to 0.5, a larger value than the ratio of 0.33 expected for a Cu 3 Pd bimetallic structure. The surface of the alloy exhibited a good LEED pattern with three-fold long-range periodicities. In STM, clusters of palladium with hexagonal arrays and Pd- Pd distances of 2.7-2.8 Å were detected. Bonding to copper perturbed the electronic properties of the atoms in the Pd clusters, shifting their 4d states toward higher binding energy with respect to the Fermi level. Further, the valence band spectrum of Cu 3 Pd(111) exhibited a line shape that was very different from those displayed by Cu(111) or Pd(111). At low pressures, the adsorption of methanol on Cu 3 Pd(111) at 300 K mainly produced CH 3 O, CO and CH x species. AP-XPS showed that most Pd atoms in the surface of the Cu 3 Pd(111) alloy interacted with the decomposition products of methanol. No significant changes were observed in the core levels of copper upon the adsorption and dissociation of methanol, suggesting that the molecule mainly interacted with Pd sites of the alloy. Reaction with hydrogen led to fast removal of CH x , C and PdC x species from Cu 3 Pd(111) at moderate (< 450 K) temperatures and prevented a CH x → C transformation. If the stability of adsorbed CH 3 O is used as a descriptor for the hydrogenation of CO 2 to methanol, Cu 3 Pd should be a much better catalyst than monometallic palladium. This may be a consequence of electronic and ensemble effects in the alloy that moderate the reactivity of Pd sites.

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

Leveraging High-resolution Molecular Composition of Soil Organic Matter to Enhance Carbon Cycling Modeling

Soils store more carbon than the atmosphere and vegetation combined, yet Earth system models still struggle to predict how this vast reservoir will respond to environmental change. A central limitation is that most soil biogeochemical models represent organic matter using bulk conceptual pools or chemically homogeneous fractions, preventing direct use of rapidly expanding molecular-scale datasets. Here we develop and test a new soil decomposition framework that explicitly integrates high-resolution information on organic matter composition. First, we construct a molecularly informed litter decomposition module in which plant inputs are partitioned into five functional compound classes—carbohydrates, proteins, lignin-like aromatics, lipids, and carbonyls—using a molecular mixing model calibrated to solid-state 13 C Nuclear Magnetic Resonance (NMR) spectra. Class-specific kinetics, lignin-dependent physical protection, and substrate-driven microbial carbon use efficiency allow the module to capture metabolic tradeoffs associated with enzyme production and nutrient limitation. We then embed this litter module within a microbially explicit whole-soil model that tracks the transformation of these compound classes through particulate organic matter, dissolved organic matter, mineral-associated organic matter, and microbial biomass. High-resolution Fourier Transform Ion Cyclotron Resonance mass spectrometry (FTICR-MS) data are used to link internal pools to measurable soil organic matter fractions and to constrain key process parameters. Applications at soil-core and ecosystem scales demonstrate that the new model reproduces observed soil respiration dynamics while providing mechanistic attribution of CO 2 fluxes to specific chemical classes and pools. Compared to existing frameworks such as the Community Land Model soil biogeochemistry module and the Millennial model, our approach maintains competitive predictive skill while substantially improving interpretability and opportunities for data–model integration. This work illustrates a viable pathway for leveraging molecular-scale observations to reduce structural uncertainty in soil carbon–climate feedback projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

COMPRESSIBLE RELATIVISTIC MAGNETOHYDRODYNAMIC TURBULENCE IN MAGNETICALLY DOMINATED PLASMAS AND IMPLICATIONS FOR A STRONG-COUPLING REGIME

In this Letter, we report compressible mode effects on relativistic magnetohydrodynamic (RMHD) turbulence in Poynting-dominated plasmas using three-dimensional numerical simulations. We decomposed fluctuations in the turbulence into 3 MHD modes (fast, slow, and Alfvén) following the procedure of mode decomposition in Cho and Lazarian, and analyzed their energy spectra and structure functions separately. We also analyzed the ratio of compressible mode to Alfvén mode energy with respect to its Mach number. We found the ratio of compressible mode increases not only with the Alfvén Mach number, but also with the background magnetization, which indicates a strong coupling between the fast and Alfvén modes. It also signifies the appearance of a new regime of RMHD turbulence in Poynting-dominated plasmas where the fast and Alfvén modes are strongly coupled and, unlike the non-relativistic MHD regime, cannot be treated separately. This finding will affect particle acceleration efficiency obtained by assuming Alfvénic critical-balance turbulence and can change the resulting photon spectra emitted by non-thermal electrons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Warming and Increased Respiration Have Transformed an Alpine Steppe Ecosystem on the Tibetan Plateau From a Carbon Dioxide Sink Into a Source

Cold region ecosystems store vast amounts of soil organic carbon (C), which upon warming and decomposition can affect the C balance and potentially change these ecosystems from C sinks to carbon dioxide (CO 2 ) sources. We quantified the decadal year-round CO 2 flux from an alpine steppe-ecosystem on the Tibetan Plateau using eddy covariance and automatic chamber approaches during a period of significant warming (0.13°C per 10 years; and 0.18°C in the non-growing season alone: 1st October to next 30th April). The results showed that ongoing climate change, mainly warming within the topsoil layers, is the main reason for the site’s change from a sink for to a source of CO 2 in the atmosphere. Non-growing-season ecosystem respiration accounted for 51% of the annual ecosystem respiration and has increased significantly. The growing seasons (1st May to 30th September) were consistent CO 2 sink periods without significant changes over the study period. Observations revealed high-emission events from the end of the non-growing season to early in the growing season (1st March to fifteenth May), which significantly (p < 0.01) increased at a rate of 22.6 g C m -2 decade -1 , ranging from 14.6 ± 10.7 g C m -2 yr -1 in 2012 to 35.3 ± 12.1 g C m -2 yr -1 in 2017. Further, structural equation modeling suggested that active layer warming was the key factor in explaining changes in ecosystem respiration, leading to significant changes in net ecosystem exchange over the period 2011–2020 and indicated that these changes have already transformed the ecosystem from a CO 2 sink into a source. These results can be used to improve our understanding of the sensitivity of ecosystem respiration to increased warming during the non-growing period.

58 GEOSCIENCES↗
DOE OSTI2021

Diol isomer revealed as a source of methyl ketene from propionic acid unimolecular decomposition

Carboxylic acids are important combustion intermediates, especially regarding the combustion of oxygenates such as ethyl esters. The purpose of this study is to discern the unimolecular decomposition pathways of one such carboxylic acid, propionic acid, using microreactor flow experiments with line-tunable photoionization mass spectrometry and infrared spectroscopy, along with high-level electronic structure calculations (CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory) and master equation theory. Here, microreactor experiments were performed at 300–1500 K, pressures decreasing from roughly 300 Torr to vacuum pressure, and around 100 μs residence times. Primary products are methyl ketene, ketene, ethylene, and methyl radical. Theory suggests the two lowest energy bond fissions are active at these conditions and responsible for the production of ketene, methyl radical, and some ethylene. Importantly, theory revealed the significance of the isomerization of propionic acid to propene-1,1-diol, and the subsequent dehydrogenation reaction of the diol as an alternative explanation for the unimolecular formation of methyl ketene. This is predicted to be the dominant thermal decomposition pathway at lower temperatures (up to ~850 K). Line-tunable VUV photons allowed for isomer resolution of the products and showed no evidence of the diol of propionic acid surviving until detection, suggesting propene-1,1-diol decomposes either to methyl ketene and water rapidly, fragments upon ionization to a distonic ion inseparable from methyl ketene through our detection methods, or never stabilizes as propionic acid well-skips directly to methyl ketene. Infrared spectroscopy showed no evidence of the decarboxylation of propionic acid at these conditions. Updated rate constants and branching ratios for propionic acid decomposition are calculated and provided for future modeling studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Second Generation RLV Space Vehicle Concept

NASA has a long history of conducting development programs and projects in a consistent fashion. Systems Engineering within those programs and projects has also followed a given method outlined by such documents as the NASA Systems Engineering Handbook. The relatively new NASA Space Launch Initiative (SLI) is taking a new approach to developing a space vehicle, with innovative management methods as well as new Systems Engineering processes. With the program less than a year into its life cycle, the efficacy of these new processes has yet to be proven or disproven. At $776M for phase 1, SLI represents a major portion of the NASA focus; however, the new processes being incorporated are not reflected in the training provided by NASA to its engineers. The NASA Academy of Program and Project Leadership (APPL) offers core classes in program and project management and systems engineering to NASA employees with the purpose of creating a "knowledge community where ideas, skills, and experiences are exchanged to increase each other's capacity for strong leadership". The SLI program is, in one sense, a combination of a conceptual design program and a technology program. The program as a whole doesn't map into the generic systems engineering project cycle as currently, and for some time, taught. For example, the NASA APPL Systems Engineering training course teaches that the "first step in developing an architecture is to define the external boundaries of the system", which will require definition of the interfaces with other systems and the next step will be to "define all the components that make up the next lower level of the system hierarchy" where fundamental requirements are allocated to each component. Whereas, the SLI technology risk reduction approach develops architecture subsystem technologies prior to developing architectures. The higher level architecture requirements are not allowed to fully develop and undergo decomposition and allocation down to the subsystems before the subsystems must develop allocated requirements based on the highest level of requirements. In the vernacular of the project cycles prior to the mid 1990's, the architecture definition portion of the program appears to be at a generic Phase A stage, while the subsystems are operating at Phase B. Even the management structure of the SLI program is innovative in its approach to Systems Engineering and is not reflected in the APPL training modules. The SLI program has established a Systems Engineering office as an office separate from the architecture development or the subsystem technology development, while that office does have representatives within these other offices. The distributed resources of the Systems Engineering Office are co-located with the respective Project Offices. This template is intended to provide systems engineering as an integrated function at the Program Level. the program management of SLI and the MAT agree that "program/project managers and the systems engineering team must work closely together towards the single objective of delivering quality products that meet the customer needs". This paper will explore the differences between the methods being taught by NASA, which represent decades of ideas, and those currently in practice in SLI. Time will tell if the innovation employed by SLI will prove to be the model of the future. For now, it is suggested that the training of the present exercise the flexibility of recognizing the new processes employed by a major new NASA program.

Bailey, Michelle↗
DOE OSTI2026

The evolution of coal porosity during pyrolysis

Gasification of coal, municipal waste, or other organic materials is a potential hydrogen source that entails complex thermal decomposition and transport processes. This study provides a multiscale analysis of these processes for sub-bituminous (Usibelli, Healy, Alaska) and lignite (Center, North Dakota) coals and provides data useful for process design. The chemistry, mineralogy, and pore structures of pyrolyzed coal and their evolution with thermal decomposition are discussed. Samples pyrolyzed at 200–1000 °C were analyzed by small-angle neutron scattering; ultra-small, small-, and wide-angle X-ray scattering; and other complementary techniques. Scanning electron microscopy showed new pores in the high-temperature-pyrolyzed material. Upon heating, the coals became progressively denser, and the concentration of hydrogen decreased. Changes in pore volume fell into three temperature ranges: an initial, low-temperature range that, for the Usibelli coal, involved an increase in overall porosity; a mid-temperature range associated with pore volume loss; and a high-temperature range associated with significant porosity increase and char formation. This transformation was paralleled by changes in fractal dimension and correlation length. The higher the pyrolysis temperature the greater the small-pore-volume fraction and overall surface area became. Pyrolysis increased the lateral size of coal crystallites, decreased the amorphous fraction, and increased the aromatics fraction and overall coal rank. Comparisons of neutron and X-ray scattering data and subsequent water uptake studies showed that pre-dried coals can re-hydrate relatively rapidly upon exposure to air, which can significantly affect the porosity calculated from small-angle-scattering data. Fits to the cumulative porosity curves provide a method for modeling the physical and chemical transformation of hydrogen-containing feedstock during gasification.

Anovitz, Lawrence {Larry} [ORNL] (ORCID:0000000226↗
DOE OSTI2022

Energetics of hydroxylbastnäsite solid solutions, La 1-$\chi$ Nd $\chi$ CO 3 OH

Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.

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

Sooting tendencies of 20 bio-derived fuels for advanced spark-ignition engines

The sooting tendencies of 20 bio-derived fuels indicated as potential blendstocks for spark-ignition engines by the Co-Optimization of Fuels & Engines (Co-Optima) Initiative are studied. The Yield Sooting Index (YSI) is used to quantify the sooting tendencies. The Co-Optima Chemical Model is used to predict the numerical YSIs for all test fuels, and these results are compared with measurements. Experimental YSIs are newly measured for 2 furans and are taken from our previous work for the other compounds. Overall, the predicted YSIs agree well with measurements, within the experimental uncertainties, except for 2,5-dimethylfuran (2,5-DMF). It is found that Polycyclic Aromatic Hydrocarbons (PAH) growth reactions have little influence on the relative soot production from each fuel and therefore on YSI predictions. A modified sensitivity coefficient formulation is proposed to evaluate reaction sensitivities specifically for YSI. This formulation is applied to the 2,5-DMF-doped flame, and six 2,5-DMF initial decomposition reactions, which lead to different major soot precursors, are identified as the most sensitive reactions. The impact of the chemical kinetic uncertainties embedded in these reactions is quantified by randomly perturbing their reaction rates within a factor of two. The resulting prediction uncertainty in the 2,5-DMF-doped flame is found to be ± 8 YSI units (±11%), which demonstrates that the test-fuel-specific decomposition reactions indeed have considerable influence on its YSI prediction. Our research suggests that more accurate kinetic parameters for fuel initial decomposition reactions can potentially improve YSI predictions effectively, without altering the YSI predictions of other fuels with significantly different molecular structures.

09 BIOMASS FUELS↗
DOE OSTI2018

Sesame-style decomposition of KS-DFT molecular dynamics for direct interrogation of nuclear models

A common paradigm used in the construction of equations of state is to decompose the thermodynamics into a superposition of three terms: a static-lattice cold curve, a contribution from the thermal motion of the nuclei, and a contribution from the thermal excitation of the electrons. While statistical mechanical models for crystals provide tractable framework for the nuclear contribution in the solid phase, much less is understood about the nuclear contribution above the % melt temperature $(C^{(nuc)}_ν≈3R)$ and how it should transition to the high-temperature limit $(C^{(nuc)}_ν\sim\frac{3}{2}R)$. Here, we describe an algorithm for extracting both the thermal nuclear and thermal electronic contributions from quantum molecular dynamics (QMD) simulations. We then use the VASP QMD package to probe thermal nuclear behavior of liquid aluminum at normal density to compare the results to semi-empirical models – the Johnson generic model, the Chisolm high-temperature liquid model, and the CRIS model. The thermal electronic energy is also calculated from the VASP results and compared to the Kohn-Sham calculations from QE-DFT simulations, Green’s function average atom calculations, and the Thomas-Fermi atomic model.

74 ATOMIC AND MOLECULAR PHYSICS↗