SEARCH · Engineering Papers
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.
The code is a python package that implements a one-step model to estimate quantum efficiency of photoemission from solid state surfaces. The analytical formulation is based on the approach proposed in Phys Rev B 95, 075439 (2017). This code extends that formalism from a model Hamiltonian to the calculation of quantum efficiency from the electronic states derived from solid-state quantum mechanical calculations obtained using the VASP code.
This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.
Triangular metal chalcogenide clusters of the form [M 3 Q 7 L 3 ]An (M = Mo or W; Q = S or Se; L = i Bu 2 NCS 2 – , (CF 3 CH 2 ) 2 NCS 2 – , i Bu 2 NCSe 2 – , or i Bu 2 PS 2 – ; An = Cl – or I – ) have been investigated as molecular analogues of layered metal dichalcogenide (MX 2 ) H 2 -evolution catalysts. These clusters have been evaluated for their relative H 2 -evolving ability under a common photolysis protocol implementing [Ru(bpy) 3 ] 2+ as chromophore and Et 3 N as sacrificial electron donor. With M constant as Mo and with constant supporting ligand, clusters with an all-sulfide core enable greater H 2 -TON than clusters with an all-selenide core. A more active catalyst is produced by [Mo 3 S 7 (S 2 CN i Bu 2 ) 3 ] + I – than its W 3 analogue with the same core sulfide composition and supporting dithiocarbamate ligands. Dichalcogenocarbamate ligands provide more active catalysts than dialkyldithiophosphate ligated clusters, and within the dichalcogenocarbamate set, greater H 2 -turnovers correlate with more-electron-donating ligands (i.e., i Bu 2 NCS 2 – > (CF 3 CH 2 ) 2 NCS 2 – > i Bu 2 NCSe 2 – ). Cluster cations with Cl – as counteranion are very similar in activity H 2 -evolving levels to identical clusters with I – , ruling out any significant interfering effect by I – upon the electron transfer relay between Et 3 N and catalyst. In the aggregate, observations are consistent with a mechanism for H 2 evolution that involves reductive extrusion of H 2 from a metal hydride intermediate.
This work proposes an analytical model considering the effects of hydrodynamic drag and energy barriers induced by liquid solvation forces to predict the in-plane translational diffusivity of a nanoparticle physically adsorbed on a wetted surface.
17th International Symposium on Solid Oxide Fuel Cells (SOFC-XVII), Virtual, July 18-23, 2021
This report describes the integration of new solid and interface-cohesive mechanics systems into MOOSE. The purpose of these new systems is to support the ability of MOOSE to run full-field crystal plasticity finite element method simulations of key material processes in high temperature metallic materials. These simulations could be used to help accurately predict the performance of key high temperature structural materials in future advanced nuclear reactor components. Previous work implemented preliminary versions of many of these systems in MOOSE Apps. The current work reports on their integration into the main MOOSE tensor mechanics module along with associated improvements to the basic formulations and numerical implementations. Finally, the report provides an example of the full-field crystal plasticity simulations now possible in MOOSE, including examples of realistic geometries requiring millions of degrees of freedom to resolve the microstructural features and macroscale geometry.
The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.
The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.
This is the user manual for the SOC Manufacturing Cost Tool spreadsheet. The manual details the use and meaning of each tab, color scheme, and spreadsheet operation. Also detailed are specific instructions for end-user modification and inputs to tailor the tool to their specific technology. To access the cost tool, please visit: <a href="https://netl.doe.gov/energy-analysis/details?id=d224ad08-6a38-402a-9cad-907db19e394f" rel="noopener noreferrer">Energy Analysis | netl.doe.gov</a>.
The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.
In order to ensure that stationary battery energy storage systems (BESSs) provide reliable energy buffering, both for power quality and economic purposes, the degradation must be considered. Cell degradation involves various side reactions and is highly dependent on its operating conditions. To accurately track cell degradation and predict its impact on battery behavior, a comprehensive physics-based degradation model based on an electrolyte phase-enhanced single particle (SP) model is developed. Key degradation physics, namely solid electrolyte interphase (SEI) layer formation and growth, Li plating on the graphite anode, and Mn dissolution on cathode of nickel-cobalt-manganese oxide (NMC622) are considered. The model is validated against experimental capacity data. The results reveal that the deposition rate of both SEI layer and Li metal increase as the charge voltage increases. At the cathode side, the solvent oxidation rate determines the Mn dissolution rate. As a result, the volume fraction of NMC622 in the cathode continually decreases at a gradually increasing rate.
Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.
We perform sign-problem-free determinant quantum Monte Carlo simulations of the optical Su- Schrieffer-Heeger model on a half-filled honeycomb lattice. In particular, we investigate the model’s semi-metal (SM) to Kekulé Valence Bond Solid (KVBS) phase transition at zero and finite temper- atures as a function of phonon energy and interaction strength. Using hybrid Monte Carlo sampling methods we can simulate the model near the adiabatic regime, allowing us to access regions of parameter space relevant to graphene. Our simulations suggest that the SM-KVBS transition is weakly first-order at all temperatures, with graphene situated close to the phase boundary in the SM region of the phase diagram. Furthermore, our results highlight the important role bond-stretching phonon modes play in the formation of KVBS order in strained graphene-derived systems.
Solid-state lithium ion batteries have enhanced thermal stability compared to batteries using liquid electrolytes. Although their practical implementation is limited by undesired side reactions between the electrode and electrolyte, different modification strategies have been proposed to stabilize these reactive interfaces. These approaches have been primarily based upon bulk materials properties, however, which may not necessarily be representative of the chemistry at the solid/solid interface. In this work, first principles calculations and X-ray scattering experiments are used to elucidate molecular-level reactivity and develop design principles for tailored interfaces based on surface and interfacial properties. Because of its well-known instability, the interface between the Li metal anode and the lithium lanthanum titanate (LLTO) solid electrolyte is applied as a model system. Density functional theory calculations are used to describe bulk, surface, and interfacial thermochemistry of the Li/LLTO system, and interfacial reconstruction is probed using ab initio molecular dynamics. These simulations of the Li/LLTO interface demonstrate facile decomposition concomitant with reduction of Ti 4+ ions. Based on further insights from computational analysis of the surface band edge positions, La 2 O 3 is proposed as an interlayer coating and is shown to provide an energetic barrier for interfacial charge transfer and decomposition reactions from X-ray reflectivity measurements and theoretical calculations. The findings in this work suggest that design strategies based on surface and interfacial properties can be used to kinetically stabilize electrode/electrolyte interfaces in solid-state batteries.
Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.
Crystalline solids assembled from superatomic building blocks are attractive functional materials due to their hierarchical structure, multifunctionality, and tunability. An interesting example is Re 6 Se 8 Cl 2 , in which the Re 6 Se 8 building blocks are covalently linked into two-dimensional (2D) sheets that are stacked into a layered van der Waals solid. It is an indirect gap semiconductor that, when heavily doped, becomes a superconductor at low temperatures. Given the finite electronic bandwidths (300–400 meV), carrier properties in this material are expected to be strongly influenced by coupling to phonons. Here, we apply angle-resolved photoemission spectroscopy to probe the valence band edge (VBE) of Re 6 Se 8 Cl 2 . We find that dispersion of the VBE is a strong function of temperature. The bandwidth is W = 120 ± 30 meV at 70 K and decreases by one order of magnitude to W ~ 10 ± 20 meV as temperature is increased to 300 K. This observation reveals the dominant polaronic effects in Re 6 Se 8 Cl 2 , consistent with the Holstein polaron model commonly used to describe molecular solids.
Not provided.