A multi-physics analysis capability for engine materials
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
SEARCH · Engineering Papers
Results for “Elemental analysis”
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.
Nonintrusive manufactured solutions for ablation
The study of heat transfer and ablation plays an important role in many problems of scientific and engineering interest. As with the computational simulation of any physical phenomenon, the first step toward establishing credibility in ablation simulations involves code verification. Code verification is typically performed using exact and manufactured solutions. However, manufactured solutions generally require the invasive introduction of an artificial forcing term within the source code such that the code solves a modified problem for which the solution is known. In this paper, we present a nonintrusive method for manufacturing solutions for a non-decomposing ablation code, which does not require the addition of a source term.
42 ENGINEERING↗
A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations
Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.
42 ENGINEERING↗
N‐Doping Effects On Electrocatalytic Water Splitting of Non‐Noble High‐Entropy Alloy Nanoparticles Prepared by Inert Gas Condensation
Abstract The unique catalytic activities of high‐entropy alloys (HEAs) emerge from the complex interaction among different elements in a single‐phase solid solution. As a “green” nanofabrication technique, inert gas condensation (IGC) combined with laser source opens up a highly efficient avenue to develop HEA nanoparticles (NPs) for catalysis and energy storage. In this work, the novel N‐doped non‐noble HEA NPs are designed and successfully prepared by IGC. The N‐doping effects of HEA NPs on oxygen evolution reaction (OER) and hydrogen evolution reaction (HER) are systematically investigated. The results show that N‐doping is conducive to improving the OER, but unfavorable for HER activity. The FeCoNiCrN NPs achieve an overpotential of 269.7 mV for OER at a current density of 10 mA cm −2 in 1.0 M KOH solution, which is among the best reported values for non‐noble HEA catalysts. The effects of the differences in electronegativity, ionization energy and electron affinity energy among mixed elements in N‐doped HEAs are discussed as inducing electron transfer efficiency. Combined with X‐ray photoelectron spectroscopy and the extended X‐ray absorption fine structure analysis, an element‐design strategy in N‐doped HEAs electrocatalysts is proposed to improve the intrinsic activity and ameliorate water splitting performance.
Chemistry↗
Detailed biomass fast pyrolysis kinetics integrated to computational fluid dynamic (CFD) and discrete element modeling framework: Predicting product yields at the bench-scale
Fast pyrolysis is an intricate process due to the variability and anisotropy of lignocellulosic biomass and the complicated chemistry and physics during conversion in a bubbling fluidized bed reactor (BFBR). The complexity of biomass fast pyrolysis lends itself well to computational fluid dynamics (CFD) and discrete element (DEM) analysis, which promises to reduce experimental time and its associated cost. This work investigated switchgrass fast pyrolysis simulated by computational fluid dynamics coupled with a discrete element method to track individual reacting biomass particles throughout a bench-scale BFBR reactor. We accounted for the fast pyrolysis chemistry through a comprehensive reaction scheme with secondary cracking reactions. We performed a three-step reduction for secondary cracking reactions to convert the full cracking scheme into a reduced scheme easily incorporated into our model. We assessed the impact of operational conditions on the steady-state yields of liquid bio-oil, non-condensable gases (NCG), at 550 °C over a range of fluidization numbers (2 – 6 Umf), reported as a ratio to the minimum fluidization velocity (Umf). At steady-state, the volatile bio-oil yield had a range of 49.3–50.4 wt%. Levoglucosan was the primary volatile component present with 21 wt% of the bio-oil while water was the second largest with 20 wt%. The reduction of the secondary reaction schemes did not appreciably affect the overall yields of switchgrass pyrolysis compared to the full secondary scheme.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Task 12 PV Sustainability - Human Health Risk Assessment Methods for PV, Part 3: Module Disposal Risks
End-of-life management is important for addressing large future photovoltaic (PV) waste volumes and conserving raw materials for use in new PV modules. In regions without regulatory mandates for PV recycling, end-of-life PV modules can be disposed in accordance with general waste laws. Given the use of various metals as raw materials in semiconductor compounds and electrical contacts in commercial PV modules, some stakeholders have raised concerns regarding potential environmental impacts if PV modules are subject to improper disposal instead of being recycled or disposed in sanitary landfills, as is required in most world regions. To evaluate these concerns, screening-level risk assessment methods are developed herein that evaluate potential human health risks from groundwater and surface (air, soil, surface water) exposure pathways. The methods estimate potential impacts from disposal of end-of-life (EOL) PV modules in non-sanitary landfills under the following worst-case conditions: no leachate collection or groundwater monitoring, no liner for preventing leachate migration, uncovered waste, and lack of stormwater management. Examining worst-case conditions allows the exploration of maximum potential risk to attempt to ensure disposal does not increase health risk above regulatory thresholds. Specifically, this report presents an analysis of potential human health risks associated with non-sanitary landfill disposal for three PV technologies, focusing on release of the highest-prioritized chemical element for each: lead (Pb) in crystalline-silicon (c-Si) PV modules, cadmium (Cd) in thin film cadmium telluride (CdTe) PV modules, and selenium (Se) in thin film copper indium selenide (CIS) PV modules. The prioritization of these chemical elements for analysis is based on stakeholder interest. Because the methodology is chemical-specific, the risk assessment results for these chemicals cannot be directly generalized to other chemicals, although the risk assessment methodology can be applied to other chemicals. If the chemicals chosen are indeed the ones presenting greatest risk, then the results herein should represent the upper bound of health risk from exposure to a single constituent.
14 SOLAR ENERGY↗
The virtual element method for linear elastodynamics models: Design, analysis, and implementation
We design the conforming virtual element method for the numerical simulation of two dimensional time-dependent elastodynamics problems. We investigate the performance of the method both theoretically and numerically. We prove the stability and the convergence of the semi-discrete approximation in the energy norm and derive optimal error estimates. We also show the convergence in the L 2 norm. The performance of the virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally seen in the p-refinement setting.
97 MATHEMATICS AND COMPUTING↗
Insight into the interfacial microstructure and chemistry of hot isostatically pressed AA6061-AA6061 bonds for U-10Mo fuel cladding application
Here, in the present study, hot isostatic pressing (HIP) bonding was carried out using AA6061-AA6061 plates for U-10Mo monolithic fuel cladding application. HIP was performed at 560°C with an applied isostatic pressure of 103 MPa for 90 min. Interfacial microstructure characterizations were performed using scanning electron microscopy (SEM) combined with electron backscatter diffraction (EBSD) to quantify precipitate fractions and study recrystallization, grain growth, and orientation relationships across the interface. EBSD results confirmed that no interfacial grain growth occurred during the HIP process. Transmission electron microscopy (TEM) was also employed to quantify composition and identify phases of various precipitates and an oxide layer that formed at the interface during HIP bonding. Microstructural analyses revealed the formation of a large fraction of Mg 2 Si precipitate and a ~50 nm thick Mg-oxide layer at the interface. Further TEM analysis confirmed that the oxide layer was not continuous, and Mg 2 Al 2 O 5 oxide particles with an average size of 18 nm were present in the oxide layer. Finite element modeling analysis was performed to study the stress distribution at the interface during HIP processing.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Continuum Damage Mechanics Modeling of High-Temperature Flaw Propagation: Application to Creep Crack Growth In 316H Standardized Specimens and Nuclear Reactor Components
Predicting creep crack growth (CCG) of flaws found during operation in high-temperature alloy components is essential for assessing the remaining lifetime of those components. While defect assessment procedures are available for this purpose in design codes, these are limited in their range of applicability. This study assesses the application of a local damage-based finite-element methodology as a more general technique for the prediction of CCG at high temperatures on a variety of structural configurations. Numerical results for stainless steel 316H, which are validated against experimental data, show the promise of this approach. This integration of continuum damage mechanics (CDM) based methodologies, together with adequate inelastic models, into assessment procedures can therefore inform the characterization of CCG under complex operating conditions, while avoiding excessive conservatism. This article shows that such modeling frameworks can be calibrated to experimental data and used to demonstrate that the degree of tri-axiality ahead of a growing creep crack affects its rate of growth. The framework is also successfully employed in characterizing CCG in a realistic reactor pressure vessel geometry under an arbitrary loading condition. These results are particularly relevant to the nuclear power industry for defect assessment and inspections as part of codified practices of structural components with flaws in high-temperature reactors.
42 ENGINEERING↗
Design, performance, and cost savings of using GF-PC additively manufactured tooling for stamping of HSS 590 sheet metal
Conventional sheet metal forming tooling in the automotive industry is made up of hardened steel and used for mass-production. Prototype tooling made of metal is durable, but it is only used for a small number of parts despite its high cost, contributing heavily to the vehicle development cost. Additive Manufacturing (AM) offers a low-cost alternative for fabricating tooling suitable for low-volume production. Here, this study investigates the feasibility of using AM polymer composite tooling for the stamping of 1.5mm thickness HSS 590 steel sheets through a two-pronged approach – experimental and numerical analysis. Sheet metal stamping experiments were conducted with glass fiber reinforced polycarbonate (GF-PC) AM polymer tooling and their performance was evaluated based on various metrics such as tool deformation and part accuracy. Finite element simulations of the stamping process were performed, which accurately captured the tool deformations seen in experiments with the use of anisotropic material models. The effective use of simulations in optimizing process parameters to achieve the desired final part geometry is also demonstrated. The GF-PC AM polymer composite tooling was found to be well suited for low volume production applications with potential cost savings.
36 MATERIALS SCIENCE↗
Composition Quantification of SiGeSn Alloys Through Time-of-Flight Secondary Ion Mass Spectrometry: Calibration Methodologies and Validation With Atom Probe Tomography
Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) is a powerful technique for elemental compositional analysis and depth profiling of materials. However, it encounters the problem of matrix effects that hinder its application. In this work, we introduce a pioneering ToF-SIMS calibration method tailored for SixGeySnz ternary alloys. SixGe1-x and Ge1-zSnz binary alloys with known compositions are used as calibration reference samples. Through a systematic SIMS quantification study of SiGe and GeSn binary alloys, we unveil a linear correlation between secondary ion intensity ratio and composition ratio for both SiGe and GeSn binary alloys, effectively mitigating the matrix effects. Extracted relative sensitivity factor (RSF) value from SixGe1-x (0.07 < x < 0.83) and Ge1-zSnz (0.066 < z < 0.183) binary alloys are subsequently applied to those of SixGeySnz (0.011 < x < 0.113, 0.863 < y < 0.935 and 0.023 < z < 0.103) ternary alloys for elemental compositions quantification. These values are cross-checked by Atom Probe Tomography (APT) analysis, an indication of the great accuracy and reliability of as-developed ToF-SIMS calibration process. Furthermore, the proposed method and its reference sample selection strategy in this work provide a low-cost as well as simple-to-follow calibration route for SiGeSn composition analysis, thus driving the development of next-generation multifunctional SiGeSn-related semiconductor devices.
36 MATERIALS SCIENCE↗
Graphite-Embedded High-Performance Insulated Metal Substrate for Wide-Bandgap Power Modules
Emerging wide-bandgap (WBG) semiconductor devices such as silicon carbide (SiC) metal–oxide semiconductor field-effect transistors (MOSFETs) and gallium nitride high-electron-mobility transistors can handle high power in reduced semiconductor areas better than conventional Si-based devices owing to superior material properties. With increased power loss density in a WBG-based converter and reduced die size in power modules, thermal management of power devices must be optimized for high performance. Here, we present a graphite-embedded insulated metal substrate (thermally-annealed-pyrolytic-graphite-embedded insulated metal substrate—IMSwTPG) designed for WBG power modules. Theoretical thermal performance analysis of graphite-embedded metal cores is presented, with design details for IMSwTPG with embedded graphite to replace a direct-bonded copper (DBC) substrate. The proposed IMSwTPG is compared with an aluminum nitride-based DBC substrate using finite-element thermal analysis for steady-state and transient thermal performance. The solutions’ thermal performances are compared under different coolant temperature and thermal loading conditions, and the proposed substrate's electrical performance is validated with static and dynamic characterization. Using graphite-embedded substrates, junction-to-case thermal resistance of SiC MOSFETs can be reduced up to 17%, and device current density can be increased by 10%, regardless of the thermal management strategy used to cool the substrate. Reduced transient thermal impedance of up to 40% of dies owing to increased heat capacity is validated in transient thermal simulations and experiments. The half-bridge power module's electrical performance is evaluated for on -state resistance, switching performance, and switching loss at three junction temperature conditions. The proposed substrate solution has minimal impact on conduction and switching performance of SiC MOSFETs.
42 ENGINEERING↗
MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes
Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.
DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗
RUScal : Software for the analysis of resonant ultrasound spectroscopy measurements
Resonant ultrasound spectroscopy is used to nondestructively measure the elastic resonances of small solids to elucidate the material's elastic properties or other qualities like size, shape, or composition. In this work, we introduce the software RUScal for the purpose of determining elastic properties by analyzing the eigenfrequencies of solid specimens with common shapes, such as rectangular parallelepipeds, cylinders (solid and hollow tube), ellipsoids, and octahedrons, as well as irregularly shaped ellipsoids that can be described analytically. All symmetry classes are supported, from isotropic to triclinic, along with the option to add or remove up to three orthogonal mirror planes as well as the ability to reorient the crystal axes with respect the sample edges via Euler angles. Additional features include tools to help find initial sets of elastic constants, including grid exploration and Monte Carlo methods, a tool to analyze frequencies as a function of sample length or crystal orientation, an error analysis tool to assess fit quality, and formatting of the input and output files for batch fitting, e.g., as a function of temperature. This software was validated with published resonant ultrasound spectroscopy data for various materials, shapes, and symmetries with noted improvements in calculation time compared to finite element methods.
47 OTHER INSTRUMENTATION↗
A layered 2D computational framework: Theory and applications to nuclear fuel behavior
Nuclear fuel performance computer codes have been developed over the last 50 years to analyze fuel behavior under various operating conditions. Traditionally, these codes used quasi-two-dimensional (also commonly known as 1.5D) representations of the fuel rod, which model the rod using a set of one-dimensional axisymmetric models that represent the behavior at specific axial positions on the rod. Modern fuel performance codes have the ability to investigate full three-dimensional (3D) effects and couple to other physics-based codes for true multiphysics simulations. However, with increasing complexity comes increasing computational costs. Many phenomena of interest involve azimuthally-varying behavior that cannot be represented using the aforementioned quasi-two-dimensional approach, but do not require the use of a full 3D model. To efficiently address these problems, there is a need for a computational framework that provides a compromise between the quasi-two-dimensional and full 3D models. In this paper, we present a new quasi-three-dimensional approach that represents the fuel as a set of 2D planar models that represent the behavior of the fuel cross-section at various axial positions. Presented here are the theory behind the methodology, test cases to illustrate proper implementation, and practical applications of its use in the BISON fuel performance code for a variety of cases in nuclear fuel analysis, such as fuel fracture, axial fuel relocation, and cladding distension and oxidation.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Chemical Pressure-Derived Assembly Principles for Dodecagonal Quasicrystal Approximants and Other Complex Frank–Kasper Phases
The structures of complex intermetallic compounds can often be interpreted in terms of assemblies of units from simpler parent phases. For example, dodecagonal quasicrystals appear, when viewed down their high-symmetry axes, as plane-filling arrangements of square and triangular tiles corresponding to the Cr 3 Si and Al 3 Zr 4 structure types, respectively. The atomic arrangements and cell-dimensions at the (100) faces of the cells of these structures provide a close geometrical match, which underlies not only dodecagonal quasicrystals and their approximants, but also the much more common σ-phase structure. In this Article, we show that such intergrowth of parent structures can arise from more than just geometrical coincidences, but can be driven by a complementary matching of atomic packing forces. DFT-Chemical Pressure (CP) analysis on elemental versions of the Cr 3 Si and Al 3 Zr 4 types reveal that in both cases arrays of positive interatomic pressures inhibit the formation of optimal contacts elsewhere in the structures. When they are lined up at the potential Cr 3 Si/Al 3 Zr 4 interfaces, however, positive pressures from the two structures interdigitate rather than coincide, providing the opportunity for the relaxation of strained interatomic contacts. That such relief is afforded by the interfaces is confirmed by CP analysis of the σ-phase (FeCr-type) structure. Building on this scheme, we introduce the CP interface function to represent how the CP features of atoms within a structure impact planes or other surfaces that could serve as interfaces between different structures. Using this function, we then explore how the favorability of interfaces between Cr 3 Si and Al 3 Zr 4 -type units is tuned by partial elemental substitution with Si, as well as their potential matches with Laves phase units. Furthermore, the emerging picture provides an account for features of the quasicrystal approximants Mn 7 VSi 2 and Mn 81.5 Si 18.5 , as well as a framework for approaching intermetallic intergrowth structures more broadly.
36 MATERIALS SCIENCE↗
Mach Conference, 2022.
There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiC/SiC) composite tubes due to their superior mechanical properties. The SiC/SiC tubes are known to remain stable even after prolonged exposure to radiation. However, there are many variables in the manufacturing of SiC/SiC composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. The intention is to identify the distribution of elastic constants of SiC/SiC tubes which would help model the performance of the novel cladding material. To do so, a sensitivity analysis will be carried out to determine the dominant variables that influence the elastic constants of SiC/SiC tubes. A chosen number of combinations of dominant variables will be sampled through the Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in the circumferential and longitudinal direction, Poisson's ratio, and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used as a surrogate to FE analysis. Once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The uncertainty can be quantified by calculating the coefficient of variation (CV) based on MCS-generated data which will be presented at the conference.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Identification of Distribution of Elastic Constants of SiC/SiC Braided Tubes
There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiCf/SiCm) composite tubes due to their superior mechanical properties. There are many parameters in the manufacturing of SiCf/SiCm composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. A sensitivity analysis will be carried out to determine the dominant parameters among the abovementioned ones. We intend to determine the variability in elastic constants of SiCf/SiCm composite tubes considering the variability in dominant parameters. A chosen number of combinations of dominant variables will be sampled through Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in circumferential and longitudinal direction, Poisson's ratio and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used a surrogate to FE analysis. once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The variability can be calculated based on MCS-generated data which will be presented at the conference.
42 ENGINEERING↗