Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Compilers”

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 1,657 records · Page 92

The origin, transport, and evolution of ash in engine particulate filters

Engine particulate filters have been widely applied across the world to control engine exhaust particulate matter (or particulate number) emissions. With increasing vehicle mileage, ash accumulation deteriorates vehicle fuel economy and complicates on-board control. Extending filter service life with ash loading has significant economic and environmental impacts. Many studies have been conducted in characterizing ash accumulation and evaluating its impacts on filter performance. However, comprehensive reviews covering all the key issues in the field are rather rare. This paper reviews the extensive prior research on filter ash, and not only summarizes the experimental observations but also elucidates the fundamental mechanisms. The review covers the areas of ash origin, accumulation, transport, evolution, and artificial acceleration methods. The previously reported data of ash properties is compiled and analyzed. Furthermore, the advantages and disadvantages of ash acceleration approaches are also discussed in detail. Based on the cumulative understanding, a few potential ways to improve ash management are discussed in this paper. In summary, the present work systematically reviews the previous observations and understanding of ash aging in particulate filters and identifies areas that need further research, which can be useful guidance for future studies.

42 ENGINEERING↗

Market analysis for the integration of new power technologies: A case study of the deployment of hybrid fossil-based generator plus energy storage (ES-FE)

This study examines the national landscape of hybridized fossil energy (FE) power plants with energy storage (ES) technologies (“ES-FE”) and presents the compilation of an ES-FE dataset, which includes over 65 ES-FE projects and concepts in the United States, comprising approximately 500 MWh of co-located ES capacity with FE power plants. This study also estimates the economic feasibility of adding ES to existing FE power plants by characterizing the potential revenues that can be generated by the ES component through flexibility and capacity value. The analysis focuses on ES technologies with 2- to 10-h. durations located in four U.S. independent system operators (ISOs): Midcontinent ISO (MISO), Electric Reliability Council of Texas (ERCOT), PJM Interconnection (PJM), and California ISO (CAISO), which have +70,000 MW of combined FE power capacity that could add ES. Annual revenues are estimated for the ES component using a what-if-analysis approach, for capacity value, price arbitrage, or ancillary services provision. The results show that annual revenues depend on the end-use storage service, wholesale electricity and capacity market prices, and ES technology operation parameters such as discharging duration and cycling frequency. When performing a sensitivity analysis, ES accrues $7–178/kW-yr. via price arbitrage and ancillary services provision in the four ISOs, and $13–92/kW-yr. when providing capacity value only in MISO and PJM. A cash flow analysis is performed to estimate the net present value (NPV) of the ES addition using a range of ES costs. The study finds that for most ES technologies considered, these revenues alone are insufficient to achieve economic feasibility. In conclusion, of the 1645 total runs analyzed, 115 had positive NPVs (7%). Therefore, other revenue streams or monetizable benefits are necessary to achieve the break-even point.

20 FOSSIL-FUELED POWER PLANTS↗

A review of energy storage for power system resilience: Functions, metrics, and applications

Aging infrastructure, increasing operational complexity, and surging electricity demand from artificial intelligence and electrification are straining the grid and heightening the risks of disruptions, making resilience a critical priority. Energy storage is increasingly deployed to provide critical power supply, fast grid support, and rapid restoration. However, current practice lacks consistent metrics and systematic methodologies to rigorously quantify the resilience benefits of storage. This paper provides a comprehensive review of energy storage in resilience enhancement, focusing on functional roles, quantification metrics, and integration strategies. A structured resilience metrics library is compiled and categorized to encompass both technical and economic performance aspects. Existing methodologies for resilience-oriented storage planning and operations are critically examined. Key technical and practical challenges are identified, and future research directions are outlined to strengthen storage contributions to grid resilience.

Benefit quantification↗

Natural attenuation of uranium in a fluvial Wetland: Importance of hydrology and speciation

A nuclear fuel fabrication facility released 43,500 kg of uranium into a riparian wetland located on the Savannah River Site between 1955 and 1988. Studies were undertaken to evaluate hydrological and geochemical processes influencing uranium accumulation in the wetland. Gamma-radiation-mapping surveys were conducted by systematically walking over the contaminated wetland with backpacks equipped with global positioning systems and NaI gamma detectors. Based on maps compiled from >700,000 gamma spectra and eight sediment uranium depth profiles, it was determined that 94% of the released uranium remained in the wetland. The uranium in the wetland is concentrated in five multi-hectare areas along the stream, accounting for ~11% of the land area adjacent to the stream. While land type (upland or wetland) and topography provided a reasonable first approximation of where much of the uranium was deposited, hydrological watershed modeling revealed that the stream velocity was especially slow through many of the hot spots. Here, using autoradiography combined with SEM/EDX measurements of contaminated sediments, surprisingly few hot particles were detected. Instead, uranium was evenly distributed throughout the sampled sediment, suggesting that dissolved uranium had bound to sediment particles that became suspended and later deposited in low energy (low flow velocity) portions of the stream. EXAFS suggested that U atoms were present as individual ions in disordered complexes within the sediment. Furthermore, linear combination analyses suggested that the predominant component of the U(VI) was adsorbed to sediment minerals (~70%) and a minor component (~30%) was associated with organic matter phases. Furthermore, these studies show that wetlands can be extraordinarily effective at binding and retaining uranium, thereby providing a natural barrier to the transport of uranium out of a watershed. However, significant anthropogenic or climatic changes to wetlands, such as those associated with flooding, forest fires, or land use, may disrupt the complex hydrological and biogeochemical balance necessary to maintain long-term immobilization of uranium.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Universal correlation for falling film evaporation heat transfer coefficients of water and seawater

Horizontal falling film evaporators are widely utilized in desalination industries to increase fresh water supply. However, universal correlations for seawater falling film evaporation under varied operational and geometrical conditions are simply unavailable in open literature. Thus, this study aims to develop such a universal correlation for both plain and enhanced tubes. The detailed heat transfer mechanisms are reviewed, and rational parameters are incorporated to develop the universal correlation. A dataset of 994 data points from 9 sources covering a wide range of conditions was compiled. These conditions include Reynolds numbers from 10 to 7235, heat fluxes from 7.7 to 208 kW/m -2 , saturation temperatures from 278 to 401 K, salinities from 0 to 60 gsalt kg -1 water, tube diameters from 15.8 to 50.8 mm, and liquid feeder height to diameter ratios from 1 to 2.25. Upon analysis, it was found that most of the recommended existing correlations exhibited poor predictive accuracy, as evidenced by larger MADs. The developed correlation in this study demonstrated the best predictive accuracy for the entire dataset, yielding a MAD of 16.8 % and an R 2 of 0.82. Furthermore, the performance of the new correlation was individually assessed across a broader spectrum of operational and design conditions, reflecting the individual conditions’ influences with an overall MAD of 20 %.

42 ENGINEERING↗

Unification of [FeFe]-hydrogenases into three structural and functional groups

[FeFe]-hydrogenases (Hyd) are structurally diverse enzymes that catalyze the reversible oxidation of hydrogen (H 2 ). Recent biochemical data demonstrate new functional roles for these enzymes, including those that function in electron bifurcation where an exergonic reaction is coupled with an endergonic reaction to drive the reversible oxidation/production of H 2 . To identify the structural determinants that underpin differences in enzyme functionality, a total of 714 homologous sequences of the catalytic subunit, HydA, were compiled. Bioinformatics approaches informed by biochemical data were then used to characterize differences in inferred quaternary structure, HydA active site protein environment, accessory iron-sulfur clusters in HydA, and regulatory proteins encoded in HydA gene neighborhoods. HydA homologs were clustered into one of three classification groups, Group 1 (G1), Group 2 (G2), and Group 3 (G3). G1 enzymes were predicted to be monomeric while those in G2 and G3 were predicted to be multimeric and include HydB, HydC (G2/G3) and HydD (G3) subunits. Variation in the HydA active site and accessory iron-sulfur clusters did not vary by group type. Group-specific regulatory genes were identified in the gene neighborhoods of both G2 and G3 Hyd. Analyses of purified G2 and G3 enzymes by mass spectrometry strongly suggests that they are post-translationally modified by phosphorylation. In conclusion, these results suggest that bifurcation capability is dictated primarily by the presence of both HydB and HydC in Hyd complexes, rather than by variation in HydA.

09 BIOMASS FUELS↗

Contribution of citizen science towards international biodiversity monitoring

To meet collective obligations towards biodiversity conservation and monitoring, it is essential that the world's governments and non-governmental organisations as well as the research community tap all possible sources of data and information, including new, fast-growing sources such as citizen science (CS), in which volunteers participate in some or all aspects of environmental assessments. Through compilation of a database on CS and community-based monitoring (CBM, a subset of CS) programs, we assess where contributions from CS and CBM are significant and where opportunities for growth exist. We use the Essential Biodiversity Variable framework to describe the range of biodiversity data needed to track progress towards global biodiversity targets, and we assess strengths and gaps in geographical and taxonomic coverage. Our results show that existing CS and CBM data particularly provide large-scale data on species distribution and population abundance, species traits such as phenology, and ecosystem function variables such as primary and secondary productivity. Only birds, Lepidoptera and plants are monitored at scale. Most CS schemes are found in Europe, North America, South Africa, India, and Australia. We then explore what can be learned from successful CS/CBM programs that would facilitate the scaling up of current efforts, how existing strengths in data coverage can be better exploited, and the strategies that could maximise the synergies between CS/CBM and other approaches for monitoring biodiversity, in particular from remote sensing. More and better targeted funding will be needed, if CS/CBM programs are to contribute further to international biodiversity monitoring.

59 BASIC BIOLOGICAL SCIENCES↗

Biomass supply chain equipment for renewable fuels production: A review

The production of renewable fuels is a critical component of global strategies to reduce greenhouse gas (GHG) emissions. Moreover, the collection of raw materials for its production can provide added benefits such as reduction of wildfire risk, additional income for farmers, and decreased disposal costs. Although there is substantial literature on design and modeling of supply chains, the authors were unable to find a single reference with the information needed for the selection and cost estimation of each type of equipment involved in the supply chain. Therefore, the goal of this research is to gather information necessary for the construction and utilization of models that might drive the identification of a feasible supply chain to produce renewable fuels at a commercial scale. The primary objectives are to 1) understand the supply chain of critical feedstocks for renewable fuels production; 2) identify the equipment commercially available for collection and adequation of feedstock; and 3) consolidate information regarding equipment cost, energy consumption, and efficiency, as well as feedstock storage and transportation systems. This paper provides a compilation for five feedstock types studied for sustainable aviation fuel production: 1) agricultural residues and grasses, 2) forest residues, 3) urban wood waste, 4) oilseeds, 5) fats, oils & greases. All the technologies involved from the field to the gate of the preprocessing or conversion unit were reviewed. The information on fats, oils & greases supply chains and equipment purposely designed for forest thinning and pruning was very limited.

: Feedstock, Collection and Adequation, Renewable ↗

Italian prototype building models for urban scale building performance simulation

Urban building energy modeling (UBEM) seeks to evaluate strategies to optimize building energy use at urban scale to support a city's building energy goals. Prototype building models are usually developed to represent typical urban building characteristics of a specific use type, construction year, and climate zone, as detailed characteristics of individual buildings at urban scale are difficult to obtain. This study investigated the Italian building stock, developing 46 building prototypes, based on construction year, for residential and office buildings. The study included 16 single-family buildings, 16 multi-family buildings, and 14 office buildings. Building envelope properties and heating, ventilation, and air conditioning system characteristics were defined according to existing building energy codes and standards for climatic zone E, which covers about half the Italian municipalities. Novel contributions of this study include (1) detailed specifications of prototype building energy models for Italian residential and office buildings that can be adopted by UBEM tools, and (2) a dataset in GeoJSON format of Italian urban buildings compiled from diverse data sources and national standards. Finally, the developed prototype building specifications, the building dataset, and the workflow can be applied to create other building prototypes and to support Italian national building energy efficiency and environmental goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analytical models of hydrogen transport in graphite

The importance of graphite-hydrogen chemical reactions to fusion, fission, and hydrogen storage applications, combined with the rapidly evolving knowledge on the underlying mechanisms, has led to the development of multiple models to describe hydrogen transport in graphite. Significant differences exist among these models, resulting from discrepancies in the modeling assumptions, intended degree of fidelity, and conditions of applicability. This paper attempts at reconciling these apparent differences by providing a comprehensive description of the constitutive equations governing hydrogen transport in graphite at high-temperature, identifying outstanding gaps in knowledge, illustrating how these different models approach them, and proposing alternative analytical formulations grounded on experimental results from hydrogen-graphite studies. Governing equations, closing relations, and simplifying assumptions are discussed for hydrogen transport at the inter-granular and intra-granular level, accompanied by compiled experimental data and illustrated energy diagrams associated to the proposed transport mechanisms. Analytical formulations are provided to reproduce competing hypotheses on the mechanisms, supporting the development of a range of computational models that can enable resolution of outstanding knowledge gaps through comparative testing against experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting concrete compressive strength using hybrid ensembling of surrogate machine learning models

This study aims to implement a hybrid ensemble surrogate machine learning technique in predicting the compressive strength (CS) of concrete, an important parameter used for durability design and service life prediction of concrete structures in civil engineering projects. For this purpose, an experimental database consisting of 1030 records has been compiled from the machine learning repository of the University of California, Irvine. The database was used to train and validate four conventional machine learning (CML) models, namely Artificial Neural Network (ANN), Linear and Non-Linear Multivariate Adaptive Regression Splines (MARS-L and MARS-C), Gaussian Process Regression (GPR), and Minimax Probability Machine Regression (MPMR). Subsequently, the predicted outputs of CML models were combined and trained using ANN to construct the Hybrid Ensemble Model (HENSM). It is observed that the proposed HENSM produces higher predictive accuracy compared to the CML models used in the present study. The predictive performance of all models for CS prediction was compared using the testing dataset and it is found that the HENSM model attained the highest predictive accuracy in both phases. Based on the experimental results, the newly constructed HENSM model is very potential to be a new alternative in handling the overfitting issues of CML models and hence, can be used to predict the concrete CS, including the design of less polluting and more sustainable concrete constructions.

36 MATERIALS SCIENCE↗

The multiple depleted mantle components in the Hawaiian-Emperor chain

Oceanic island basalts are targeted for geochemical study because they provide a direct window into mantle composition and a wealth of information on the dynamics and timescales associated with Earth mixing. Previous studies mainly focused on the shield volcanic stage of oceanic islands and the more fusible, enriched mantle components that are easily distinguished in those basalts. Mantle depleted compositions are typically more difficult to resolve unless large amounts of this material participated in mantle melting (e.g., mid-ocean ridges), or unique processes allow for their compositions to be erupted undiluted, such as very small degrees of melting of a source with minimal fusible enriched components (e.g., rejuvenated basalts) or as xenoliths (e.g., abyssal peridotites). Mantle depleted components, defined here as material with low time-integrated Rb/Sr (low 87 Sr/ 86 Sr) and high time-integrated Sm/Nd and Lu/Hf ratios (high 143 Nd/ 144 Nd and 176 Hf/ 177 Hf) relative to primitive mantle, derive from a potentially very large volume reservoir (up to 80% of the mantle), and therefore need adequate characterization in order estimate the composition of the Earth and mantle-derived melts. Here, this review focuses on mantle depleted compositions in oceanic island basalts using the Hawaiian-Emperor chain as a case study. The Hawaiian-Emperor chain is the ~6000 km long geological record of the deeply sourced Hawaiian mantle plume, active for >81 Myr. Hawaiian volcanism evolves through four volcanic stages as a volcano traverses the Hawaiian plume: alkalic preshield, tholeiitic shield (80–90% volcano volume), alkalic postshield (~1%), and silica undersaturated rejuvenated (<0.1%). We report Pb-Sr-Nd-Hf isotope compositions and trace element concentrations of three rejuvenated Northwest Hawaiian Ridge basalts and compare them to an exhaustive compiled dataset of basalts from the Hawaiian Islands to the Emperor Seamounts. The Northwest Hawaiian Ridge (NWHR) includes 51 volcanoes spanning ~42 m.y. between the bend in the Hawaiian-Emperor chain and the Hawaiian Islands where there is no high-precision isotopic data published on the rejuvenated-stage over ~47% of the chain. NWHR and Hawaiian Island rejuvenated basalts are geochemically similar, indicating a consistent source for rejuvenated volcanism over ~12.5 million years. In contrast, shield-stage basalts from the oldest Emperor Seamounts are more depleted in isotopic composition (i.e., higher 176 Hf/ 177 Hf, and 143 Nd/ 144 Nd with lower 87 Sr/ 86 Sr and 208 Pb*/ 206 Pb*) and trace element concentrations (i.e., much lower concentrations of highly incompatible elements) than all other depleted Hawaiian basalts younger than the bend, including NWHR rejuvenated basalts. The strongly depleted source for the oldest Emperor Seamounts (>70 Ma) was likely related to interaction with the Kula-Pacific-Izanagi mid-ocean ridge spreading system active near the Hawaiian plume in the Late Cretaceous. In contrast, the incompatible trace element ratios of NWHR rejuvenated basalts require a distinct source in the Hawaiian mantle plume that was imprinted by ancient (>1 Ga) partial melting, likely ancient recycled oceanic lithosphere. This review of the geochemistry of Hawaiian depleted components documents the need for the sampling of multiple distinctive depleted compositions, each preferentially melted during specific periods of Hawaiian plume activity. This suggests that the composition of depleted components can evolve during the lifetime of the mantle plume, as observed for enriched components in the Hawaiian mantle plume. Changes in the composition of depleted components are dominantly controlled by the upper mantle tectonic configurations at the time of eruption (i.e., proximity to a mid-ocean ridge), as this effect overwhelms the signal imparted by potentially sampling different lower mantle components through time.

07 ISOTOPE AND RADIATION SOURCES↗

Optimized thermodynamic properties of REE aqueous species (REE 3+ and REEOH 2+ ) and experimental database for modeling the solubility of REE phosphate minerals (monazite, xenotime, and rhabdophane) from 25 to 300 °C

Rare earth elements (REE) are critical elements found in monazite, xenotime, and hydrated REE phosphates which typically form in hydrothermal mineral deposits. Accurate predictions of the solubility of these REE phosphates and the speciation of REE in aqueous fluids are both key to understanding the controls on the transport, fractionation, and deposition of REE in natural systems. Previous monazite and xenotime solubility experiments indicate the presence of large discrepancies between experimentally derived solubility constants versus calculated solubilities by combining different data sources for the thermodynamic properties of minerals and aqueous species at hydrothermal conditions. In this study, these discrepancies were resolved by using the program GEMSFITS to optimize the standard partial molal Gibbs energy of formation (Δ f G° 298 ) of REE aqueous species (REE 3+ and REE hydroxyl complexes) at 298.15 K and 1 bar while keeping the thermodynamic properties fixed for the REE phosphates. A comprehensive experimental database was compiled using solubility data available between 25 and 300 °C. The latter permits conducting thermodynamic parameter optimization of Δ f G° 298 for REE aqueous species. Optimal matching of the rhabdophane solubility data between 25 and 100 °C requires modifying the Δ f G° 298 values of REE 3+ by 1–6 kJ/mol, whereas matching of the monazite solubility data between 100 and 300 °C requires modifying the Δ f G° 298 values of both REE 3+ and REEOH 2+ by ~15–31 kJ/mol and ~2–10 kJ/mol, respectively. For xenotime, adjustments of Δ f G° 298 values by 1–26 kJ/mol are only necessary for the REE 3+ species. The optimizations indicate that the solubility of monazite in acidic solutions is controlled by the light (L)REE 3+ species at <150 °C and the LREEOH 2+ species at >150 °C, whereas the solubility of xenotime is controlled by the heavy (H)REE 3+ species between 25 and 300 °C. Based on the optimization results, we conclude that the revised Helgeson-Kirkham-Flowers equation of state does not reliably predict the thermodynamic properties of REE 3+ , REEOH 2+ , and likely other REE hydroxyl species at hydrothermal conditions. We therefore provide an experimental database (ThermoExp_REE) as a basic framework for future updates, extensions with other ligands, and optimizations as new experimental REE data become available. As a result, the optimized thermodynamic properties of aqueous species and minerals are available open access to accurately predict the solubility of REE phosphates in fluid-rock systems.

58 GEOSCIENCES↗

Plastic additives in the ocean: Use of a comprehensive dataset for meta-analysis and method development

In excess of 13,000 chemicals are added to plastics (‘additives’) to improve performance, durability, and production of plastic products. They are categorized into numerous chemical classes including flame retardants, light stabilizers, antioxidants, and plasticizers. While research on plastic additives in the marine environment has increased over the past decade, there is a lack of methodological standardization. To direct future measurement of plastic additives, we compiled a first-of-its-kind dataset of literature assessing plastic additives in marine environments, delineated by sample type (plastic debris, seawater, sediment, biota). Using this dataset, we performed a meta-analysis to summarize the state of the science. Currently, our dataset includes 217 publications published between 1978 and May 2023. The majority of publications analyzed plastic additives in biota collected from Europe and Asia. Analyses concentrated on plasticizers, brominated flame retardants, and bisphenols. Common sample preparation techniques included Solvent - Agitation extraction for plastic, sediment, and biota samples, and Solid Phase Extraction for seawater samples with dichloromethane and solvent mixtures including dichloromethane as the organic extraction solvent. Finally, most analyses were performed utilizing gas chromatography/mass spectrometry. There are a variety of data gaps illuminated by this meta-analysis, most notably the small number of compounds that have been targeted for detection compared to the large number of additives used in plastic production. The provided dataset facilitates future investigation of trends in plastic additive concentration data in the marine environment (allowing for comparison to toxicity thresholds) and acts as a starting point for optimizing and harmonizing plastic additive analytical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical estimates of flammability bounds for thin condensed fuel diffusion flames in microgravity using detailed models of chemistry and radiation

Recently, U-shape flammability maps have been constructed showing minimal oxygen vs. flame strain for opposed flame spread in micro-gravity by Olson and Ferkul. The U-shape defines the limiting flammability bounds from radiative extinction and flame blow-off. Here, the minimum of the U corresponds to the minimum possible oxidizer concentration where burning can occur, and is an important quantity of interest for fire safety. While high strain extinction bounds have been well analyzed, low strain radiative extinction has not. To estimate low strain extinction, in this study an analytical theory is developed based on thin flame theory coupled with a heat and mass transfer model for solid fuels. A reaction progress variable based on the Damköhler number is adapted in the theory to account for incomplete combustion at high strain rates and enable the capturing of the full flammability map. The analytical model is compared to a one dimensional numerical model w/ detailed chemical kinetics and coupled radiation heat transfer in planar and spherical geometries. The flammability maps are then qualitatively compared to experimental extinguishment data compiled by Olson and Ferkul for cylindrical rods of PMMA showing similar trends. The results show the newly developed analytics capture the radiative extinction bound compared to the numerical model and qualitatively agrees with microgravity data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Rapid multiphase-field model development using a modular free energy based approach with automatic differentiation in MOOSE/MARMOT

In this paper, we present a novel rapid phase-field model development capability using a modular free energy based approach in the finite element based MOOSE framework. The phase-field equations are implemented in a general form as separate code form the model specific data such as the thermodynamic free energy density and mobility functions. This isolates researchers from complicated code and allows them to focus on implementing new material models. In addition, using a native symbolic differentiation algorithm, all required derivatives of the free energy density are automatically determined, removing sources for potential human error, and guaranteeing optimally preconditioned equation systems through fully populated Jacobian matrices. Free energy terms contributing to a phase-field model are abstracted into self contained objects that can be combined at will at simulation run time. Combining multiple chemical and mechanical free energy contributions allows the construction of coupled phase-field, mechanics, and multiphase models. Through just-in-time compilation, we greatly reduce the computational time required to evaluate the differentiated expressions. Lastly, we demonstrate the new capability on a variety of example simulations.

97 MATHEMATICS AND COMPUTING↗

9Cr steel visualization and predictive modeling

The goals of this work were to develop tools to visualize and characterize materials information, to explore new alloy compositions and/or processing, to better predict tensile strength and other mechanical properties. The 9–12 wt% Cr (9Cr) steel test data were compiled from several sources. About 3000 records with 30 predictor variables represent 82 unique steel alloy compositions, variations of thermo-mechanical processing steps and temperatures (homogenization, normalization and tempering cycles), test conditions and outcomes. Detailed data processing steps such as visualization and exploratory analysis, including univariate and bivariate analysis, different clustering techniques used to segment the data, statistical post-hoc analysis to verify significance of the findings, feature engineering to identify predictors of importance, and predictive modeling with cross-validation were performed. The outcome of analysis at each step was reviewed in the context of the domain knowledge of this class of steel, to see if there were underlying physical mechanisms that explain statistical relationships. Data analytics techniques and their parameters were fine-tuned to facilitate interpretation of the results as aligned with the insights from domain experts. Here, the ensemble predictive modeling using Random Forest regressor and post-hoc means comparison corroborated domain knowledge on the role of Co which is known to increase strength of steel alloys through solid-solution strengthening and to affect diffusion of alloying elements and precipitates. Alloys with Co content of 0.7–8 wt% had significantly higher mean strength than the ones without (while having Cr locked within 10–11 wt% range in either group and keeping other compositional element ratios fixed). End products of the computational techniques exploration are presented as the tools that can be used in iterative workflow of materials development and testing.

36 MATERIALS SCIENCE↗