Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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,387 records · Page 77

Hierarchical morphologies in co-sputter deposited thin films

Co-depositions of immiscible alloy films at specific processing conditions have yielded hierarchical microstructures which consist of distinct features at multiple length scales, often agglomerates and concentration modulations on the sub-micrometer-scale and fine nanoprecipitates in a matrix on the nanoscale. The present work examined a series of immiscible alloy systems: Cu-Mo, Cu-Ag, Cu-Fe, Cu-Ta, Mo-Ag, Cu-Mo-Ag, to determine the kinetic conditions favorable for hierarchical organization and the formation mechanism of such structures. Thin films of six immiscible systems were sputter co-deposited over a range of deposition rates from 0.12 to 2 nm/s and various temperatures from 400 to 800°C. The resulting microstructures indicate that hierarchical structures form with sufficient disparity in kinetic energy between the constituent atoms, one species being highly mobile (A) and the other relatively immobile (B). This condition arises typically at elevated deposition temperatures and reduced deposition rates but is also alloy dependent. The hierarchical structures form during deposition via phase separation and self-organization processes across the multiple length scales. The adatoms diffuse on the film surface with the highly mobile species swiftly agglomerating into A-rich domains within which B-rich nanoprecipitates form, often self-organizing into periodic arrays. The smallest B-rich nanoprecipitates in the A-rich domains are found to be coherent and in a metastable crystal structure (B taking the structure of A), but coarser precipitates that exhibit the equilibrium structure of B element. The A-rich domains are surrounded by a B-rich matrix that phase-separates into a concentration modulated structure. In conclusion, the observations are interpreted via a model incorporating material properties and process parameters.

36 MATERIALS SCIENCE↗

AGR-1, AGR-2 and AGR-3/4 Release-to-Birth Ratio Data Analysis

A series of Advanced Gas Reactor (AGR) irradiation tests is being conducted in the Advanced Test Reactor at Idaho National Laboratory in support of development and qualification of tristructural isotropic (TRISO) low enriched fuel used in the High Temperature Gas-cooled Reactor (HTGR). Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder shrouded by a steel shell. These capsules are instrumented with thermocouples embedded in the graphite enabling temperature control. AGR configuration and irradiation conditions are based on prismatic HTGR technology distinguished primarily by the use of helium coolant, a low-power-density ceramic core capable of withstanding very high temperatures, and TRISO coated particle fuel. Thus, these tests provide valuable irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, and support development and validation of fuel performance and fission product transport models and codes. The release-rate-to-birth-rate ratio (R/B) for each of fission product isotopes (i.e., krypton and xenon) is calculated from release rates in the sweep gas flow measured by the germanium detectors used in the AGR Fission Product Monitoring (FPM) System installed downstream from each irradiated capsule. Birth rates are calculated based on the fission power in the experiment and fission product generation models. Thus, this R/B is a measure of the ability of fuel kernel, particle coating layers, and compact matrix to retain fission gas atoms preventing their release into the sweep gas flow, especially in the event of particle coating failures that occurred during AGR-2 and AGR-3/4 irradiations. The major factors that govern gaseous radioactive decay, diffusion, and release processes are found to be material diffusion coefficient, temperature, and isotopic decay constant. For each of all AGR capsules, ABAQUS-based three-dimensional finite-element thermal models are created to predict daily averages of fuel compact temperatures for the entire irradiation period, which are used in establishing the R/B correlation with temperature and decay constant. This correlation can be used by reactor designers to estimate fission gas release from postulated failed fuel particles in HTGR cores, which is the key safety factor for fuel performance assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

AGR-2 and AGR-3/4 Release-to-Birth Ratio Data Analysis

A series of Advanced Gas Reactor (AGR) irradiation tests is being conducted in the Advanced Test Reactor at Idaho National Laboratory in support of development and qualification of tristructural isotropic (TRISO) low enriched fuel used in the High Temperature Gas-cooled Reactor (HTGR). Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder shrouded by a steel shell. These capsules are instrumented with thermocouples embedded in the graphite enabling temperature control. AGR configuration and irradiation conditions are based on prismatic HTGR technology distinguished primarily by the use of helium coolant, a low-power-density ceramic core capable of withstanding very high temperatures, and TRISO coated particle fuel. Thus, these tests provide valuable irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, and support development and validation of fuel performance and fission product transport models and codes. The release-rate-to-birth-rate ratio (R/B) for each of fission product isotopes (i.e., krypton and xenon) is calculated from release rates in the sweep gas flow measured by the germanium detectors used in the AGR Fission Product Monitoring (FPM) System installed downstream from each irradiated capsule. Birth rates are calculated based on the fission power in the experiment and fission product generation models. Thus, this R/B is a measure of the ability of fuel kernel, particle coating layers, and compact matrix to retain fission gas atoms preventing their release into the sweep gas flow, especially in the event of particle coating failures that occurred during AGR-2 and AGR-3/4 irradiations. The major factors that govern gaseous radioactive decay, diffusion, and release processes are found to be material diffusion coefficient, temperature, and isotopic decay constant. For each of all AGR capsules, ABAQUS-based three-dimensional finite-element thermal models are created to predict daily averages of fuel compact temperatures for the entire irradiation period, which are used in establishing the R/B correlation with temperature and decay constant. This correlation can be used by reactor designers to estimate fission gas release from postulated failed fuel particles in HTGR cores, which is the key safety factor for fuel performance assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Preheating of Cold, High Moisture Particulate Biomaterials to Reduce Drying Time

Forest Concepts led a two-year project to reduce fossil fuel consumption and improve the energy efficiency of dryers through the application of radio frequency (RF) energy impartation to preheat cold-wet bulk particulate biomass materials. A comprehensive biomass heating and drying mathematical model was developed that enables simulation of drying systems for biomass as well as most other bulk porous particulate materials. Experimental and modeling results conclude that energy consumption for both RF preheating and hot-air preheating are essentially equal for cold, high-moisture feedstocks that are above fiber saturation and therefore not economically beneficial for those conditions. Further results suggest that RF heating may have a direct energy consumption benefit in the late stages of drying and in case of preheating biomass that has a moisture content below fiber saturation. This work was funded in part by the US Department of Energy under contact DE-EE0009130.

09 BIOMASS FUELS↗

Defect-mediated metastability and carrier lifetimes in polycrystalline (Ag,Cu)(In,Ga)Se 2 absorber materials

Using a combination of optical and electrical measurements, we develop a model for metastable defects in Ag-alloyed Cu(In,Ga)Se 2 , one of the leading thin film photovoltaic materials. By controlling the pre-selenization conditions of the back contact prior to the growth of polycrystalline (Ag,Cu)(In,Ga)Se 2 absorbers and subsequently exposing them to various stresses (light soaking and dark-heat), we explore in this paper the nature and role of metastable defects on the electro-optical and photovoltaic performance of high-efficiency solar cell materials and devices. Positron annihilation spectroscopy indicates that dark-heat exposure results in an increase in the concentration of the selenium–copper divacancy complex (V Se –V Cu ), attributed to depassivation of donor defects. Deep-level optical spectroscopy finds a corresponding increase of a defect at Ev+0.98 eV, and deep-level transient spectroscopy suggests that this increase is accompanied by a decrease in the concentration of mid-bandgap recombination centers. Time-resolved photoluminescence excitation spectroscopy data are consistent with the presence of the V Se –V Cu divacancy complex, which may act as a shallow trap for the minority carriers. Light-soaking experiments are consistent with the V Se –V Cu optical cycle proposed by Lany and Zunger, resulting in the conversion of shallow traps into recombination states that limit the effective minority carrier recombination time (and the associated carrier diffusion length) and an increase in the doping density that limits carrier extraction in photovoltaic devices.

14 SOLAR ENERGY↗

Ultrafast (1‐5 sec) Lamination of Perovskite Solar Cells With Self‐Encapsulation Using Rapid Joule Heating

Perovskite solar cells (PSCs) are traditionally fabricated using sequential layer‐by‐layer deposition, in which each layer of the device is processed on top of the preceding layer. This constrains the processing techniques and selection of transport layer materials that can be used in the solar cell. To overcome these challenges, two half‐cells can be processed independently and then diffusion‐bonded through a lamination process. However, current lamination processes for perovskite solar cells suffer from relatively long process times, which can limit throughput when moving toward high‐volume manufacturing. In this study, a custom platform was designed for rapid‐joule heating of perovskite materials and devices. This enabled more than a 99% reduction in lamination time from 26 min to 1 s. Perovskite samples that were laminated in 1 s exhibited comparable values of percent bonded area, interfacial toughness, grain domain size, and X‐ray diffraction spectra to those laminated in greater than 10 min. As a proof‐of‐concept, 18.3% efficient devices were successfully laminated in 5 s. A transient heat transfer model was developed to describe the relationship between the perovskite temperature and the electrical power supplied to the heaters, establishing a baseline for predicting processing conditions in large‐scale manufacturing systems. Ultra‐fast lamination provides a pathway toward scalable roll‐to‐roll or sheet‐to‐sheet manufacturing of PSCs.

heat transfer↗

Interfacial mixing in high-energy-density matter with a multiphysics kinetic model

We have extended here a recently developed multispecies, multitemperature Bhatnagar-Gross-Krook model [Haack et al., J. Stat. Phys. 168, 822 (2017)], to include multiphysics capabilities that enable modeling of a wider range of physical conditions. In terms of geometry, we have extended from the spatially homogeneous setting to one spatial dimension. In terms of the physics, we have included an atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, but for larger length and time scales and for much higher temperatures than can be simulated using molecular dynamics. Relative to molecular dynamics, the kinetic model greatly extends the temperature regime and the spatiotemporal scales over which we are able to model. In our numerical results we observe hydrogen from the ablator material jetting into the fuel during the early stages of the implosion and compare the relative size of various diffusion components (Fickean diffusion, electrodiffusion, and barodiffusion) that drive this process. We also examine kinetic effects, such as anisotropic distributions and velocity separation, in order to determine when this problem can be described with a hydrodynamic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The ABCflux database: Arctic–boreal CO 2 flux observations and ancillary information aggregated to monthly time steps across terrestrial ecosystems

Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO 2 ) fluxes in terrestrial ecosystems across the rapidly warming Arctic–boreal zone (ABZ) have provided valuable information but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic–boreal CO 2 fluxes (ABCflux) that aggregates in situ measurements of terrestrial net ecosystem CO 2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures), and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19 % of the monthly observations), snow diffusion (3 %) and eddy covariance (78 %) techniques. The largest number of observations were collected during the climatological summer (June–August; 32 %), and fewer observations were available for autumn (September–October; 25 %), winter (December–February; 18 %), and spring (March–May; 25 %). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO 2 fluxes and to better estimate the terrestrial ABZ CO 2 budget.

59 BASIC BIOLOGICAL SCIENCES↗

Final Project Report (Oct 2014-Dec 2017): Zeolite Membranes for Krypton/Xenon Separation from Spent Nuclear Fuel Reprocessing Off-Gas

The overall focus of this project is to develop and understand SAPO-34 zeolitic membranes that can separate mixtures of radioisotope krypton-85 and xenon released as off-gases during used nuclear fuel recycling. The primary advantage of separating 85 Kr from Xe is to reduce the volume of radioactive waste for storage. The second advantage is the revenue generated from the sale of high-purity Xe. Zeolite membranes are attractive because of their much lower energy requirements relative to cryogenic distillation, and their high resistance to radiation degradation. We report the detailed study of silicoaluminophosphate zeolite SAPO-34 materials and membranes for this application, due to hypothesized favorable molecular sieving properties. In the 3-year Mission Support project, we developed a novel, high-performance, low-energy intensity, lower-cost zeolite membrane process for Kr/Xe separation during SNF processing; and investigated the underlying molecular adsorption and transport processes in both ‘idealized’ and ‘realistic’ operating conditions to develop reliable synthesis-structure-property relationships for such membranes. Adsorption and diffusion measurements on SAPO-34 crystals indicate their potential for use in Kr-Xe separation membranes, but also highlight competing effects of adsorption and diffusion selectivity. SAPO-34 membranes synthesized on α-alumina substrates via steam-assisted conversion seeding and hydrothermal growth are characterized in detail, with Kr permeances 26 GPU and ideal Kr/Xe selectivities >20 at 298 K after thickness reduction. Post-synthesis cation exchange shows large (>50%) increases in selectivity at ambient or slight sub-ambient conditions. In addition, we confirm that SAPO-34 membrane is stable under radiation exposure and the impact of radiation exposure on membrane performance would not be substantial. We also successfully synthesized hollow-fiber SAPO-34 membranes with the same performance levels as the disk-type and tubular membranes. This important development will allow a very compact and low-cost Kr/Xe separation system. Finally, a detailed process calculation for techno-economic analysis was performed by integrating Maxwell-Stefan model into cross-flow membrane system, in order to estimate the required number of membrane stages and the total cost.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-E-J-61: Electrodynamics and Nano-Scale Fluid Dynamics in Protein Localization of Nuclear Pore Complexes

Purpose: To develop a simulation to catalyze a reevaluation of common assumptions about 3 dimensional diffusive processes and help cell biologists gain a more nuanced, intuitive understanding of the true physical hurdles of protein signaling cascades. Furthermore, to discuss the possibility of intracellular electrodynamics as a critical, unrecognized component of cellular biology and protein dynamics that is necessary for optimal information flow from the cell membrane to the nucleus. Methods: The Unity 3D gaming physics engine was used to build an accurate virtual scale model of the cytoplasm within a few hundred nanometers of the nuclear membrane. A cloud of simulated pERK proteins is controlled by the physics simulation, where diffusion is based on experimentally measured values and the electrodynamics are based on theoretical nano-fluid dynamics. The trajectories of pERK within the cytoplasm and through the 1250 nuclear pores on the nuclear surface is recorded and analyzed. Results: The simulation quickly demonstrates that pERKs moving solely by diffusion will rarely locate and come within capture distance of a nuclear pore. The addition of intracellular electrodynamics between charges on the nuclear pore complexes and on pERKs increases the number of successful translocations by allowing the electro-physical attractive effects to draw in pERKs from the cytoplasm. The effects of changes in intracellular shielding ion concentrations allowed for estimation of the “capture radius” under varying conditions. Conclusion: The simulation allows a shift in perspective that is paramount in attempting to communicate the scale and dynamics of intracellular protein cascade mechanics. This work has allowed researchers to more fully understand the parameters involved in intracellular electrodynamics, such as shielding anion concentration and protein charge. As these effects are still far below the spatial resolution of currently available measurement technology this simulation is invaluable to researchers to visualize and study current hypotheses concerning protein dynamics.

60 APPLIED LIFE SCIENCES↗

Physics based Reliability Models for Supercritical CO 2 Turbomachinery Components

CSP plants using sCO 2 power cycle can potentially achieve high thermal conversion efficiency at low capital cost due to compact turbo-machinery and other components. An sCO 2 expander and improved heat exchanger is expected to provide a major stepping stone for achieving CSP power at $0.06/kW-hr LCOE, energy conversion efficiency > 50%, and total power block cost < $1,200/kW installed. Hybrid gas bearings (HGB) and Dry gas seals (DGS) are two key components that are critical for robust and efficient operation of future sCO 2 power cycles. From a system perspective, HGBs can provide substantial benefits (~3 to 5%) to the modular CSP operation, while use of DGSs limits the leakage losses as well as the windage losses in the generator and can increase overall generator efficiency by ~ 8%. The life limiting mechanisms of HGBs and DGSs in high pressure, high temperature sCO2 environment were not well understood and hence this program aims at developing physics-based performance prediction and life prediction models for these two components. The output of the performance models serves as input to the life models. The performance models were validated using experimental data and using data from literature. DGS and HGB models have shed light on some risks associated with its operation in sCO2 expander and will guide more robust designs tailored towards sCO 2 cycle applications. A DGS was tested in 1MWe sCO2 turbo expander to study its performance and validate the performance models. To understand the effect of high pressures, high temperatures and sCO 2 environment on crack initiation, crack propagation and LCF/HCF life of these components, a novel experimental setup was developed. The LCF test data was used to calibrate and validate the life prediction models. The effect of sCO 2 corrosion on the life of Ni-based super-alloys is quantified using long duration experiments. Advanced microstructure and spectroscopic analyses conducted shed light on some key differences between IN718 and IN617 in terms of oxidation morphology, chemical species diffusion and trapping, the formation of protective corrosion resistant layers and changes in surface properties. A novel crack initiation model was developed which takes into account strain distribution, temperature, hold times and activation energy at the gas-metal interface and predicts the number of cycles to initiate the crack in the Inconel alloy. An advanced 3D fracture mechanics-based crack propagation model was also developed which takes into account the transient loading mission of the component and predict how fast the crack will grow and what will be its final shape and size. High energy X-Ray tomography measurements were conducted to study and validate 3D morphology of crack growth inside the specimen. These two models along with chemical kinetics models obtained from corrosion tests can be combined to predict the total life of the component under given long term operating conditions. A Bayesian Hybrid Modeling (BHM) probabilistic framework was developed to quantify uncertainty in physics-based performance models and in crack initiation and crack propagation models, to calibrate the models and to validate the models with statistical confidence. Advanced coupled physics multi-phase flow models were developed to study condensation, boiling and phase change behavior of CO 2 under various sCO 2 turbomachinery and heat exchanger conditions. A novel experimental setup was also designed to observe phase change of CO 2 during its transition across liquid-gas-supercritical state lines.

14 SOLAR ENERGY↗

Physics based Reliability Models for Supercritical CO 2 Turbomachinery Components

CSP plants using sCO 2 power cycle can potentially achieve high thermal conversion efficiency at low capital cost due to compact turbo-machinery and other components. An sCO 2 expander and improved heat exchanger is expected to provide a major stepping stone for achieving CSP power at 0.06/kW-hr LCOE, energy conversion efficiency > 50 % ,and total power block cost< $1,200/kW installed. Hybrid gas bearings (HGB) and Dry gas seals (DGS) are two key components that are critical for robust and efficient operation of future sCO 2 power cycles. From a system perspective, HGBs can provide substantial benefits (~3 to 5%) to the modular CSP operation, while use of DGSs limits the leakage losses as well as the windage losses in the generator and can increase overall generator efficiency by ~ 8%. The life limiting mechanisms of HGBs and DGSs in high pressure, high temperature sCO 2 environment were not well understood and hence this program aims at developing physics-based performance prediction and life prediction models for these two components. The output of the performance models serves as input to the life models. The performance models were validated using experimental data and using data from literature. DGS and HGB models have shed light on some risks associated with its operation in sCO 2 expander and will guide more robust designs tailored towards sCO 2 cycle applications. A DGS was tested in 1MWe sCO 2 turbo expander to study its performance and validate the performance models. To understand the effect of high pressures, high temperatures and sCO 2 environment on crack initiation, crack propagation and LCF/HCF life of these components, a novel experimental setup was developed. The LCF test data was used to calibrate and validate the life prediction models. The effect of sCO 2 corrosion on the life of Ni-based super-alloys is quantified using long duration experiments. Advanced microstructure and spectroscopic analyses conducted shed light on some key differences between IN718 and IN617 in terms of oxidation morphology, chemical species diffusion and trapping, the formation of protective corrosion resistant layers and changes in surface properties. A novel crack initiation model was developed which takes into account strain distribution, temperature, hold times and activation energy at the gas-metal interface and predicts the number of cycles to initiate the crack in the Inconel alloy. An advanced 3D fracture mechanics-based crack propagation model was also developed which takes into account the transient loading mission of the component and predict how fast the crack will grow and what will be its final shape and size. High energy X-Ray tomography measurements were conducted to study and validate 3D morphology of crack growth inside the specimen. These two models along with chemical kinetics models obtained from corrosion tests can be combined to predict the total life of the component under given long term operating conditions. A Bayesian Hybrid Modeling (BHM) probabilistic framework was developed to quantify uncertainty in physics-based performance models and in crack initiation and crack propagation models, to calibrate the models and to validate the models with statistical confidence. Advanced coupled physics multi-phase flow models were developed to study condensation, boiling and phase change behavior of CO 2 under various sCO 2 turbomachinery and heat exchanger conditions. A novel experimental setup was also designed to observe phase change of CO 2 during its transition across liquid-gas-supercritical state lines.

14 SOLAR ENERGY↗

Planetary Ices and the Linear Mixing Approximation

Here, the validity of the widely used linear mixing approximation (LMA) for the equations of state (EOSs) of planetary ices is investigated at pressure–temperature conditions typical for the interiors of Uranus and Neptune. The basis of this study is ab initio data ranging up to 1000 GPa and 20,000 K, calculated via density functional theory molecular dynamics simulations. In particular, we determine a new EOS for methane and EOS data for the 1:1 binary mixtures of methane, ammonia, and water, as well as their 2:1:4 ternary mixture. Additionally, the self-diffusion coefficients in the ternary mixture are calculated along three different Uranus interior profiles and compared to the values of the pure compounds. We find that deviations of the LMA from the results of the real mixture are generally small; for the thermal EOSs they amount to 4% or less. The diffusion coefficients in the mixture agree with those of the pure compounds within 20% or better. Finally, a new adiabatic model of Uranus with an inner layer of almost pure ices is developed. The model is consistent with the gravity field data and results in a rather cold interior (${T}_{\mathrm{core}}\sim 4000$ K).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Reactor Scale-Up Methodology from Lab-Scale to Pilot-Scale Operations: Numerical Modeling of THFA Dehydration to DHP in Packed-Bed Reactors

This manuscript discusses developing a model-based scale-up methodology for a successful technology transfer of gas-phase catalytic reactors from a lab-scale to a pilot-scale operation. The manuscript demonstrates the methodology for gas-phase dehydration of tetrahydrofurfuryl alcohol (THFA) to dihydropyran (DHP) process over commercial Al 2 O 3 catalysts. A two-dimensional reactor model was developed using COMSOL Multiphysics 6.1 software. The model solves heat and mass transport equations in bed-scale and particle scales simultaneously. This powerful feature enables accurate prediction of the heat and mass transfer limitations in pilot-scale reactors, if any exists. Further, the model uses isothermal lab-scale experimental data to derive and validate the reaction chemistry, flow fields and boundary conditions. The model was then scaled-up to project conversion, selectivity, yield and formation rate of DHP in a pilot-scale reactor. The results highlight the complex nature of chemistry, heat, and mass transfer effects in lab-scale and pilot-scale reactors. The model results inform the possible operational limitations of the pilot-scale reactor and design strategies to improve process efficiency. Although the scale-up approach is explained through the THFA dehydration process, the methodology is applicable to any catalytic packed-bed reactor models for a successful process scale-up.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Identifying multiple synergistic factors on the susceptibility to stress relaxation cracking in variously heat-treated weldments

The 347H austenitic stainless steel has been widely used for pressure vessels and pipeline (PVP) applications due to its excellent creep and corrosion resistance, which fit ideally to the harsh conditions in petrochemical industries, fossil fuel or nuclear power plants, and modern energy storages. However, a failure mode has been commonly observed with cracks emerging at the heat affected zone (HAZ) of weldments during post-weld heat treatment (PWHT) or under intermediate to high temperature service conditions. This phenomenon is termed as Stress Relaxation Cracking (SRC) since the purpose of PWHT is to relieve the welding-induced residual stress fields, or as Stress Age Cracking (SAC) if failure happens during service. A leading literature explanation of this failure suggests that the residual stress relaxation and the precipitation dissolution and/or re-precipitation occur in the same temperature range, which can lead to locally high strains and thus to crack at the grain boundaries. Since in situ spatial measurements of residual stress fields, microstructural evolution, and failure processes are nearly infeasible, this work recourses to a micromechanical finite element framework that models the high temperature failure as the nucleation and growth of grain boundary cavities, whereas various parameters such as thermomechanical loading history and its evolution, the competition of grain-interior dislocation creep and grain-boundary diffusion in failure lifetime, and microstructural heterogeneities (such as the precipitate free zone near grain boundaries) can be quantitatively incorporated. It can be concluded from these microstructure-explicit simulations that an accurate knowledge of residual stress evolution and a carefully calibrated set of material constitutive parameters are the essential prerequisites for lifetime predictions. The understanding of individual governing factors also leads to a mechanistic interpretation of the observed SRC susceptibility C-curves. In conclusion, these results suggest that the criticality of residual stress evolution, but not the precipitation-induced local strains, be the leading factor for SRC.

347H stainless steel weldments↗

Laboratory Verification of Electron-Scale Reconnection Regions Modulated by a Three-Dimensional Instability

We report that during magnetic reconnection in collisionless space plasma, the electron fluid decouples from the magnetic field within narrow current layers, and theoretical models for this process can be distinguished in terms of their predicted current layer widths. From theory, the off-diagonal stress in the electron pressure tensor is related to the thermal noncircular orbit motion of electrons around the magnetic field lines. This stress becomes significant when the width of the reconnecting current layer approaches the small characteristic length scale of the electron motion. To aid in situ spacecraft and numerical investigations of reconnection, the structure of the electron diffusion region is here investigated using the Terrestrial Reconnection Experiment. In agreement with the closely matched kinetic simulations, laboratory observations reveal the presence of electron-scale current layer widths. Although the layers are modulated by a current-driven instability, three-dimensional simulations demonstrate that it is the off-diagonal stress that is responsible for breaking the frozen-in condition of the electron fluid.

79 ASTRONOMY AND ASTROPHYSICS↗