Search NASA⌕ Search

SEARCH · Search NASA

Results for “constitutive flow model”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analyzing Hardware and Software Common Cause Failures in Digital Instrumentation and Control Systems using Dual Error Propagation Method

This paper develops a methodology for quantifying software common cause failures (CCFs) in digital instrumentation and control (I&C) systems of nuclear power plants. To support the transition of analog I&C systems to digital in nuclear power plants, probabilistic risk assessment (PRA) techniques are used. The hardware components of the I&C systems have reliability databases that can be used in the PRA studies. However, the failure data for redundant software components of the systems is sparse. Failure of components constitutes a CCF, wherein two or more components or systems fail due to a single shared cause and coupling mechanism. This paper proposes a quantification approach that can simultaneously model hardware and software components, incorporate the CCFs of software systems in the models, and bridge the gap between the failure quantification of models and the development of CCF parametric databases. We demonstrate the dual error propagation method (DEPM) by developing I&C systems failure models for a representative digital reactor trip system. The DEPM models are built to simulate the control and data flows within the systems and can accommodate failure states. By expanding DEPM to software CCFs, we generated alpha factor parameter estimates for each of the modeled error propagation mechanisms.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2020

Projecting Antarctica's contribution to future sea level rise from basal ice shelf melt using linear response functions of 16 ice sheet models (LARMIP-2)

Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Measurement of Wind Loading on Heliostats at the Crescent Dunes Power Plant: An Overview

The cost of solar collectors constitutes almost one third of the total cost of a CSP plant. One of the ongoing challenges in the design of these collectors is wind loading on mirrors, support structures, and drives. A particular challenge is dynamic wind loading, caused by the turbulent wind flow. To date, the design of solar collector structures has relied on wind tunnel experiments and numerical simulations that do not entirely capture the dynamic effects observed at scale. The CSP industry has shown increased interest in validating the idealized assumptions with measurements obtained in operational settings to improve wind load assumptions and increase reliability and cost-efficiency of the collector design. Further, performance models need realistic assumptions about wind loading and its impact on optical performance. In a parabolic trough field campaign, NREL successfully collected a wealth of long-term, high-resolution wind and loads data [2], that are publicly available and that can be used for the above-mentioned purposes. Heliostats are impacted differently by wind than parabolic troughs, due to their different shape, size, and field layout. To study the impact of wind and turbulence on heliostats, we initiated another field campaign in an operational power- tower plant, Crescent Dunes, in Nevada, USA. In this work, we present an overview of the measurements, first results, and potential implications of wind driven loads on Heliostats.

concentrated solar power↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2021

Accounting for the effect of dislocation climb-mediated flow on the anisotropy and texture evolution of Mg alloy, AZ31B

The current study explores the role that dislocation climb has in mediating plasticity a Mg alloy at moderately elevated temperatures. Interrupted tensile tests were performed on samples of Mg alloy, AZ31B, sheet in the O temper condition over a range of strain rates (10 -5 to 10 -1 s -1 ) and temperatures (20–350°C) along the rolling and transverse directions. Experimental measurements of the resulting strain anisotropy and texture evolution were used as constraints during a parametric study employing a new crystal plasticity model (VPSC-CLIMB), which explicitly accounts for the kinematics of dislocation climb. The results reveal that the climb of basal (a) dislocations is not only important for dislocation recovery, but also demonstrate that climb accommodates a significant fraction of the strain in conditions where a power-law creep-type constitutive response prevails. This work does not discredit the notion that non-basal slip of (a) and (c+a) dislocations is important over a wide range of temperatures and strain rates. However, it demonstrates that the activation of dislocation climb as the mechanistic change within the power law regime provides an explanation for a wide range of observations, including the simultaneous reduction in strain anisotropy, slowed texture evolution, and rapid increase in strain rate sensitivity. Finally, it is hypothesized that these conclusions may even apply to cases in which grain boundary sliding and/or dynamic recrystallization are observed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Simulations of biomass compression-screw feeding using a compressible non-Newtonian constitutive model

There is global interest in the conversion of biomass into sustainable low-carbon-footprint fuels and chemicals as an alternative to non-renewable fossil feedstocks. Feeding biomass solids into pressurized reactors is one of the key steps in biomass conversion. Predicting mechanical failure and energy requirements for this step helps avoid upstream processing bottlenecks and enables efficient operation of a biorefinery. Here, in this work, we developed a predictive computational model for biomass screw feeders that capture the highly viscous, non-Newtonian and compressible behavior of biomass slurries. Biomass compressible behavior is formulated by an equation of state and the non-Newtonian rheology is represented by a density-dependent viscosity model. Experimental data from two compression screw-feeder systems are presented as a validation for our model. Our model successfully predicted the location of the compressed biomass “plug”, biomass flow rate, and the required torque at different operating conditions for the experimental conditions studied in this work.

09 BIOMASS FUELS↗
DOE OSTI2021

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2022

Amplification-free nucleic acid detection with a fluorescence-based waveguide biosensor

Early detection of pathogens using nucleic acids in clinical samples often requires sensitivity at the single-copy level, which currently necessitates time-consuming and expensive nucleic acid amplification. Here, we describe 1) a redesigned flow cell in the shape of a trapezoid-subtracted geometric stadium, and 2) modified experimental procedures that allow for the measurement of sub-attomolar analytes in microliter quantities on a fluorescence-based waveguide biosensor. We verified our instrumental sensitivity with a 200-μL sample of a fluorescent streptavidin conjugate at 100 zM (100 zeptomolar, or 100·10 −21 mol L −1 ) and theoretically explored the applicability of this modified sensing platform in a sandwich immunoassay format using a Langmuir adsorption model. We present assays that demonstrate specific detection of synthetic influenza A DNA (in buffer) and RNA (in saliva) oligonucleotides at the single-copy level (200 μL at 10 zM) using a fluorescent molecular beacon. Lastly, we demonstrate detection of isolated genomic influenza A RNA at a clinically relevant concentration. This work constitutes a sensitivity improvement of over twelve orders of magnitude compared to our previous nucleic acid detection work, illustrating the significant enhancements that can be gained with optimized experimental design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

A visco-plastic framework for interface processes in sedimentary reservoir rocks at HPHT conditions

As energy operations face the challenge of reservoirs at ever-increasing depths, modelling the response of reservoir rocks at high pressure and high temperature (HPHT) conditions is a crucial step for successfully unlocking new resources. At these conditions, rocks can experience thermal and pressure sensitivity, as well as admit internal phase changes at the pore-solid interfaces that determine their macroscopic response to external loading. It is therefore expected that constitutive laws of Geomechanics should be enriched to accommodate such effects. In this work, a multi-physics constitutive theory for sedimentary rocks is proposed within the general framework of viscoplasticity. The viscosity of the material is assumed to be a function of the temperature, pore-pressure and energy required to alter the inter-granular interfaces. This energy is expressed through the chemical potentials of the phases involved during processes like chemical dissolution/precipitation of the interfaces or mechanical debonding. The resulting flow law and corresponding stress equilibrium are coupled to the energy and mass conservation laws, constituting a closed system of equations, which is solved using the Finite Element simulator REDBACK. A series of numerical tests for performance and calibration is then successfully performed against different types of reservoir rocks, confining pressures and temperatures (from room temperature to over 800K). Here, we show that the mechanical response of sedimentary porous rocks at strains usually achieved in laboratory testing can be explained by accounting for the interface processes taking place at the cementitious material bonding the grains, a process that can also determine the brittle-to-ductile transition of sedimentary rocks.

58 GEOSCIENCES↗
DOE OSTI2021

Learning viscoelasticity models from indirect data using deep neural networks

In this study, we propose a novel approach to model viscoelasticity materials, where rate-dependent and non-linear constitutive relationships are approximated with deep neural networks. We assume that inputs and outputs of the neural networks are not directly observable, and therefore common training techniques with input–output pairs for the neural networks are inapplicable. To that end, we develop a novel computational approach to both calibrate parametric and learn neural-network-based constitutive relations of viscoelasticity materials from indirect displacement data in the context of multiple-physics systems. We show that limited displacement data holds sufficient information to quantify the viscoelasticity behavior. We formulate the inverse computation – modeling viscoelasticity properties from observed displacement data – as a PDE-constrained optimization problem and minimize the error functional using a gradient-based optimization method. The gradients are computed by a combination of automatic differentiation and implicit function differentiation rules. The effectiveness of our method is demonstrated through numerous benchmark problems in geomechanics and porous media transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Phase field-volumetric lattice Boltzmann model of ion uptake in porous nuclear waste form materials under continuous flow

The flow field within the mesopores of sorbent particles plays a crucial role in radionuclide diffusion and ion uptake kinetics, thus, impacting the overall performance of porous nuclear waste form materials. To fundamentally understand the influence of microstructures and material properties on the radionuclide absorption and retention processes requires a coupled multi-physics model that considers the advection and diffusion within the flow field, the reaction at liquid-solid interfaces, and finally, the solid-state diffusion within a complex nanoporous medium. Here, this study employs the volumetric lattice Boltzmann method (VLBM) to accurately and efficiently calculate the steady state velocity field inside the mesopores of sorbent particles. The obtained velocity field is then utilized to calculate the advection of ions in the steady flow. A phase field (PF) model of ion uptake is used to describe the reaction occurring at the solid-liquid interface and diffusion inside the porous medium. The integrated PF-VLBM model is verified in terms of the mass conservation and numerical efficiency and validated qualitatively with experimental observation data. Then, it is applied to study the influence of thermodynamic and kinetic properties, as well as flow field conditions on the ion uptake kinetics. The numerical results demonstrate that the ion uptake kinetics in porous particles has three distinct stages, which is in agreement with the observations in continuous flow experiments. In the first stage, the kinetics is predominantly controlled by the flow field and ion diffusivity in the liquid phase. The kinetics in the second stage is primarily governed by ion diffusivity in the solid phase. In the third stage the system reaches a dynamic equilibrium with a net zero uptake flux at the interface. It is also found that porous structures significantly affect the efficiency and capacity of ion uptake. The simulation results can help to understand the physics behind the observed ion uptake kinetics in experiments and to facilitate the development of constitutive equations that can account for heterogeneous microstructures in engineering performance codes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

The weakly nonlinear response and nonaffine interpretation of the Johnson–Segalman/Gordon–Schowalter model

We derive new analytical solutions for the non-affine Johnson-Segalman/Gordon-Schowalter (JS/GS) constitutive equation with a general relaxation kernel in medium-amplitude oscillatory shear (MAOS) deformation. The results show time-strain separable (TSS) nonlinearity, therefore providing new physically-meaningful interpretation to the heuristic TSS nonlinear parameter in MAOS. The upper-convected, lower-convected, and corotational Maxwell models are all subsets of the results presented here. The model assumes that the microscale elements causing stress in the material slip compared to the continuum deformation. We introduce a visualization of the non-affine deformation field that acts on stress-generating elements to reinforce the hysical interpretation of the JS/GS class of models. Lastly, a case study is presented where previously published results, from fitting TSS models to MAOS data, can be re-interpreted based on the concept of non-affine motion of the JS/GS framework.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗