Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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,225 records · Page 68

Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator

In our previous work Walton and Tran (2023), numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the long time dynamics of the equation given the multiple blow-up time phenomenon. In this second work in the series, we propose numerical treatments for the complete 3-wave kinetic equations, in which the complete, much more complicated collision operators are fully considered based on a novel conservative form of the equation. Here we then derive an implicit finite volume scheme to solve the equation. The new discretization uses an adaptive time-stepping method which allows for the simulations to be carried to very long times. Our computed solutions are compared with previously derived long-time asymptotic estimates for the decay rate of total energy of time-dependent solutions of 3-wave kinetic equations and found to be in excellent agreement.

97 MATHEMATICS AND COMPUTING↗

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗

Conversion of CO 2 from power plant into CaCO 3 nanoparticles

Carbon dioxide (CO 2 ), a main composition of flue gas, represents a significant and largely untapped carbon resource. Herein, mediated by glycine (Gly), we captured and converted CO 2 into CaCO 3 nanoparticles using real flue gas from a power plant, demonstrating for the first time the feasibility of using amino acid to convert CO 2 from power plant flue gasses. The method did not require extraneous energy and CaCO 3 nanoparticles with a size of ∼25 nm were obtained. Moreover, the potential toxicity of CO 2 -converted nanoparticles was investigated. It appeared that both the initial CO 2 loading and the carbamate percentage significantly influence the shape and size of the CaCO 3 particles. Our method was also proven effective for flue gas with varying CO 2 concentrations (4 %, 12 %, and 20 %). By tuning flue gas bubbling time and flow rate to achieve consistent CO 2 loading and carbamate levels, we produced CaCO 3 nanoparticles with similar shapes and sizes across all CO 2 concentrations studied. In addition, our data indicated that although real flue gas contains small amounts of gases like oxygen and CO, they insignificantly influence the shape and size of our nanoparticles but did impact the phase component of CaCO 3 . In conclusion, the toxicity experiments found that CaCO 3 nanoparticles produced from both real flue gas and simulated flue gas exhibited concentration- and time-dependent effects on cell viability.

36 MATERIALS SCIENCE↗

Experimental investigation into coal wettability changes caused by reactions with scCO 2 -H 2 O

Geological CO 2 sequestration (GCS) can help mitigate global warming and enhance methane recovery from coal beds. However, few studies have linked the effects of CO 2 to surface chemistry changes controlling wetting behavior in deep coal beds. Contact angles (CAs) of CO 2 /N 2 -high volatile bituminous coal-water systems were measured under different temperatures and pressures. The surface chemistry and physical structure of coals were characterized to investigate changes in physicochemical properties and their relations with wettability after reactions. For N 2 treatment, the time-dependence of static and dynamic CAs were insignificant, ranging within 4°. For gaseous CO 2 treatment, the static CAs and the average advancing angles increased slightly. With supercritical (sc) CO 2 , both the static and dynamic CAs increased significantly, and θ adv changed to intermediate-wet (92°). Reactions with minerals exposed to scCO 2 resulted in greater surface roughness and heterogeneity, greater contact angle hysteresis and more surface sites occupied by scCO 2 rather than H 2 O. Increases in hydrophobic functional groups and decreases in hydrophilicity were shown by FTIR spectra, reflecting the shedding of polar oxygen-containing functional groups, reduction of hydrogen bonds, and increasing percentage of hydrocarbons. XRD patterns obtained following scCO 2 -treatment showed that crystallite growth and molecular polymerization were higher toward graphite-like. The calculated structural parameters of functional groups and crystallites both showed elevated coal rank. Changes in crystallite structure, notably higher carbon content and decreased negative surface charge, are unfavorable for water-wetting. Finally, this study contributes to understanding surface chemistry changes responsible for decreased wettability during CO 2 -enhanced coal bed methane recovery and GCS in coal reservoirs.

01 COAL, LIGNITE, AND PEAT↗

Spectroscopic and computational investigations of organometallic complexation of group 12 transition metals by methanobactins from Methylocystis sp. SB2

Methanotrophic bacteria catalyze the aerobic oxidation of methane to methanol using Cu-containing enzymes, thereby exerting a modulating influence on the global methane cycle. To facilitate the acquisition of Cu ions, some methanotrophic bacteria secrete small modified peptides known as “methanobactins,” which strongly bind Cu and function as an extracellular Cu recruitment relay, analogous to siderophores and Fe. In addition to Cu, methanobactins form complexes with other late transition metals, including the Group 12 transition metals Zn, Cd, and Hg, although the interplay among solution-phase configurations, metal interactions, and the spectroscopic signatures of methanobactin-metal complexes remains ambiguous. In this study, the complexation of Zn, Cd, and Hg by methanobactin from Methylocystis sp. strain SB2 was studied using a combination of absorbance, fluorescence, extended x-ray absorption fine structure (EXAFS) spectroscopy, and time-dependent density functional theory (TD-DFT) calculations. We report changes in sample absorbance and fluorescence spectral dynamics, which occur on a wide range of experimental timescales and characterize a clear stoichiometric complexation dependence. Mercury L3-edge EXAFS and TD-DFT calculations suggest a linear model for Hg--S coordination, and TD-DFT suggests a tetrahedral model for Zn 2+ and Cd 2+ . We observed an enhancement in the fluorescence of methanobactin upon interaction with transition metals and propose a mechanism of complexation-hindered isomerization drawing inspiration from the wild-type Green Fluorescent Protein active site. Collectively, our results represent the first combined computational and experimental spectroscopy study of methanobactins and shed new light on molecular interactions and dynamics that characterize complexes of methanobactins with Group 12 transition metals.

59 BASIC BIOLOGICAL SCIENCES↗

Revealing the synergistic effects of sequential and simultaneous dual beam irradiations in tungsten via in-situ TEM

Elucidating the synergistic effects of different energetic beams on the radiation response of nuclear materials is critical for developing an improved methodology for their evaluation when exposed to extreme environments. This article describes in-situ sequential (He implantation followed by Kr irradiation and vice-versa) and simultaneous (heavy ion Kr irradiation and He implantation) dual beam irradiations performed on tungsten at 1223 K. Dislocation loop density, average area, and total loop damage as a function of irradiation history and dose/fluence are quantified. The loop Burgers vectors and cavity damage (cavity density, size and total change in volume) are also determined at the final maximum dose for each condition. The loop damage evolution was different in all cases, with the smallest loop damage observed in the simultaneous experiment. Annihilation of Kr generated dislocation loops during He implantation in the Kr + He experiment was an unanticipated observation that may be explained by the dynamic evolution of dislocation loop sink strengths and time-dependent defect fluxes. Dislocation loop raft formation, denuded zones near extended defects, and cavity damage are compared across the different conditions. The phenomena observed and discussed in this work will stimulate further experimental and computational modeling activities leading to improved fundamental understanding of the irradiation response of nuclear materials under reactor-similar environments.

36 MATERIALS SCIENCE↗

A systematic exploration of the hydrolysis products of the uranium trioxide polymorphs and their optical vibrational spectra

Study of the uranium trioxide (UO 3 )-water system is complex with inconclusive results and limited details in the literature. The UO 3 system is home to at least seven structural polymorphs and an amorphous phase. The proposed hydrolysis products of UO 3 are just as numerous, yet investigations of these alteration products are sporadic and generally antiquated, thus requiring systematic investigations. Recent developments in the understanding of UO 3 phase space, aided by improvements in analytical and computational techniques, necessitate more modern investigations into the uranyl hydroxide family and their naturally occurring mineral counterparts. We present findings from a systematic investigation of the products formed via hydrothermal reactions of common UO 3 polymorphs and discuss how the equatorial coordination of the uranyl/uranyl-like ions within the UO 3 precursors leads to differences in the optical vibrational spectra of the resulting hydrolysis products. The hypo-stochiometric nature of α-UO 2 (OH) 2 allows for the formation of multiple unique uranyl sites and a distortion of the unit cell to a lower symmetry. This study provides, for the first time, an analysis of β-UO 2 (OH) 2 using modern techniques and instrumentation (Raman/infrared spectroscopy and powder x-ray diffraction) and lays a foundation for future time-dependent investigations into the structural dependence of the hydrolysis kinetics of the UO 3 phases. In conclusion, given the prevalence of UO 3 at both ends of the nuclear fuel cycle, an understanding of its behavior with water has applications ranging from nuclear forensics to waste management and environmental transport.

Uranium trioxide↗

A Kinetic Model of the Long-Term Corrosion of Glass-Ceramic Materials

Multiphase waste forms show promise for increased waste loading and for the ability to dispose of contaminated solid and particulate waste through direct densification. However, achieving predictive capability for long-term durability of multiphase waste forms, and thus assessing their possible deployment, requires expanding the current, limited knowledge base. Here, we describe the development of a corrosion model of a two-phase waste form consisting of crystals of known volume fraction embedded in a glass matrix. This model accounts for the dissolution of both the crystalline and glass phases as well as the hydration of the glass phase through an ion exchange reaction. Because of the large difference in solubility between the two phases, the reactive surface of the crystalline phase is a function of the extent of dissolution of the glass phase in this model. Model parameterization was performed using corrosion data, such as from single-pass flow-through tests, for the individual phases. The parameterized corrosion model was evaluated against static dissolution test data for a glass-ceramic multiphase waste form. This evaluation demonstrated the model’s ability to reproduce the time-dependent release of key tracers of glass and crystalline phase dissolution. Hence, the development of a kinetic model provides a pathway for long-term durability predictions and thus the use of multiphase waste forms in nuclear cleanup missions.

Kerisit, Sebastien N.↗

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗

Spatially-resolved lithiation dynamics from operando X-ray diffraction and electrochemical modeling of lithium-ion cells

In this study, energy dispersive X-ray diffraction is used to profile the time evolution of ordered Li x C 6 phases in solid electrodes of lithium-ion cells charged at rates between 0.2 and 4.7C (where 1C corresponds to full discharge in 1 h). The methods for quantifying lithium concentration in these phases from the acquired diffraction patterns are described. Compact expressions for time-dependent concentration gradients in the solid electrodes using orthogonal polynomial expansions are presented. Experimentally, these gradients persisted in lithiated graphite electrodes even after the cells rested at open-circuit for over 9 h. A multiphase electrochemical model of graphite intercalation captured many of the observed behaviors, including the progression of phase transitions and the persistent gradients at zero current. However, the magnitude of concentration gradients in both the oxide cathode and graphite anode is underestimated by the model, even at moderate currents.

25 ENERGY STORAGE↗

A pseudo-kinetic model to simulate phase changes in gas hydrate bearing sediments

Herein, modeling of the phase transitions anticipated in gas hydrate bearing sediments (GHBS) is critical for a proper understanding of time-dependent changes in states and volumes (e.g. the production of methane from this type of soils). We propose a new pseudo-kinetic approach to simulate the typical phase changes anticipated in GHBS, using published experimental results involving gas hydrate dissociation that are the basis of a widely used kinetic model. The proposed pseudo-kinetic model is formulated in the pressure-temperature (P-T) plane and assumes a rate of gas hydrate dissociation (or formation) proportional to the distance between the current state and the phase boundary. The model consists of only one parameter and is simple to implement in numerical simulators. A similar concept is used to model ice formation/thawing phenomena, but based on the ice/liquid-water phase boundary. We implemented the pseudo-kinetic model in a fully coupled thermo-hydro-chemo-mechanical (THCM) finite element code and validated it against experimental results performed on the dissociation of synthetic gas hydrate. We also evaluated the pseudo-kinetic model using synthetic cases covering several scenarios associated with gas hydrate formation/dissociation and ice formation/thawing. The model successfully reproduced the gas production test from a natural GHBS core from Korea (scaled gas venting experiment over 14 h), and also the formation of gas hydrate and ice in permafrost in Alaska (over 2 × 10 6 years). -The analyses show the versatility of the proposed pseudo-kinetic approach by applying it to model the different types of phase transitions typically encounter in GHBS. The simple formulation, easy implementation in numerical simulator, and reduced number of parameters (only one per phase change) make this model an attractive option for simulating phase transformations in problems involving GHBS.

58 GEOSCIENCES↗

Investigation of competitive sorption and plasticization of hyperaged CANAL ladder polymers for acid gas purification

Identifying membrane materials that have exceptional separation performance and stability to complex CO 2 -containing mixtures is a pressing topic in separation science. In this work, the membrane separation performance for a recently discovered class of contorted polymers synthesized via catalytic arene-norbornene annulation (CANAL) polymerization is presented. These CANAL polymers achieve high CO 2 /CH 4 selectivity of 68 after physical aging (up to ∼1 year), which significantly augments the size-sieving capabilities of the membranes. Binary CO 2 /CH 4 and ternary H 2 S/CO 2 /CH 4 testing result in a 41 % and 50 % enhancement in selectivities, respectively, for hyperaged (∼1 year) contorted CANAL polymers, highlighting their size-sieving capabilities. The remarkably high CO 2 /CH 4 mixed-gas and combined acid gas (CAG, (CO 2 +H 2 S)/CH 4 ) selectivities of 88 and 95, respectively, for CANAL-Me-S 5 F in particular surpass both the 2018 CO 2 /CH 4 mixed-gas and CAG upper bounds. Performance stability was also investigated for high concentrations of CO 2 , revealing a reduction in CO 2 /CH 4 mixed-gas selectivity without compromising CO 2 permeability, suggesting strong sorption of CO 2 but minimal plasticization effects for short testing periods. In addition, time-dependent plasticization shows negligible effects on CO 2 /CH 4 mixed-gas selectivity despite an increase in CO 2 permeability when exposed to high concentrations of plasticizing CO 2 over an extended period of 170 h. This study provides valuable insights into hyperaged CANAL polymers and their performance in practical industrial processes.

03 NATURAL GAS↗

Incorporating dynamic recrystallization into a crystal plasticity model for high-temperature deformation of Ti-6Al-4V

During hot deformation of (α + β) titanium alloys, the simultaneous action of strain and temperature in the (α + β) regime facilitates dynamic recovery, dynamic recrystallization (DRX), and phase transformations via non-equilibrium paths. DRX is manifested in the form of fine recrystallized α or β grains. Here, in the present study, a two-phase crystal plasticity finite element framework (CP-DRX) has been developed which incorporates DRX kinetics into the crystal plasticity (CP) model to predict the flow characteristics of Ti-6Al-4V alloys during thermo-mechanical processing. The CP slip system parameters, as well as elastic properties from both α and β phases of Ti, are calibrated for different strain rate conditions. An EBSD-informed two-phase microstructure representation has been utilized in the CP-DRX framework to explore the time-dependent evolution of DRX microstructure and crystal orientation for different strain rate conditions. The proposed CP-DRX can capture the evolution of crystal orientation and plastic flow stress–strain response of polycrystalline Ti-6Al-4V during the deformation process. Furthermore, the proposed model is able to capture the softening behavior, observed in average stress–strain response from experiments performed using a Gleeble thermomechanical simulator and predict the recrystallization texture.

(α+ β)-Titanium (Ti) Alloys↗

Investigation on gate oxide reliability under gate bias screening for commercial SiC planar and trench MOSFETs

This paper evaluates the impact of the fast high gate-voltage screening technique on gate oxide reliability of commercial 1.2 kV 4H-SiC power metal-oxide-semiconductor field-effect transistors (MOSFETs) with planar and trench gate structures. The measurements are conducted on packaged devices with the intent that the results will be applicable to wafer-level screening. The threshold voltage (V th ) of SiC MOSFETs is measured before and after various screening treatments and recovery process. In addition, constant-voltage time-dependent dielectric breakdown (TDDB) measurements are performed on SiC MOSFETs to obtain the intrinsic lifetime of gate oxide. The objective of this study is to determine the optimal screening conditions and improve screening efficiency without degradation of gate oxide reliability, such as V th shift and reduced oxide intrinsic lifetime. Due to the differences in the gate oxidation process and structural design of SiC planar and trench MOSFETs, the two types of devices exhibit different oxide reliability in the screening process. The recommended screening conditions obtained in this paper reveal that SiC trench MOSFETs can accept higher screening voltages compared to SiC planar MOSFETs. Hence, it can be concluded that more efficient screening techniques can be adopted for SiC MOSFETs with thicker gate oxide to meet the requirements of industrial applications.

36 MATERIALS SCIENCE↗

Direct observation of key aluminum hydroxide prenucleation oligomers for gibbsite nucleation and crystallization in sodium aluminate solution by liquid ToF-SIMS

The mechanism of gibbsite (aluminum hydroxide) crystallization from highly alkaline solutions such as Bayer liquors remains poorly understood, where aluminum (Al) transforms from largely tetrahedrally coordinated aluminate monomers in sodium aluminate solutions into a network of octahedra in gibbsite crystals. A variety of traditional analytical approaches applied to this system do not readily reveal the presence of higher-order oligomeric intermediates. To overcome this limitation, we employed in-situ liquid Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) to examine the Al species present in concentrated sodium aluminate solutions favoring crystallization of gibbsite or sodium aluminate. A complex mixture of Al oligomers was found. By comparing the change in the relative concentration of Al oligomers with +1 and -1 charge, we were able to identify three major Al oligomer candidates, iso-tetramers, iso-pentamers, and cyclic-hexamers, for the nucleation and crystallization of gibbsite. The concentrations of iso-tetramers and iso-pentamers significantly surpass those of cyclic-hexamers. Time-dependent in-situ Raman spectroscopy analysis indicated that the appearance of gibbsite coincided with the peak concentration of these oligomers. The Density-functional theory (DFT) calculation suggests that the formation of iso-oligomers is more favorable than that of cyclic-hexamers. The combined results suggest that iso-tetramers and iso-pentamers play the most substantial role in the nucleation and growth of gibbsite in the sodium aluminate solutions. Our findings also suggest that the oligomers that promote gibbsite crystallization are more stable in dilute sodium aluminate solutions, making these solutions particularly suitable for efficient gibbsite crystallization. In conclusion, our study fills a major knowledge gap in understanding Al speciation that leads to the nucleation and crystallization of gibbsite in concentrated sodium aluminate solutions.

Bayer liquor↗

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO 3

This work herein demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO 3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO 3 .

36 MATERIALS SCIENCE↗

Development of a Reference Database for Beta-Delayed Neutron Emission

Beta-delayed neutron emission is important for nuclear structure and astrophysics as well as for reactor applications. Significant advances in nuclear experimental techniques in the past two decades have led to a wealth of new measurements that remain to be incorporated in the databases. In this work, we report on a coordinated effort to compile and evaluate all the available β-delayed neutron emission data. The different measurement techniques have been assessed and the data have been compared with semi-microscopic and microscopic-macroscopic models. The new microscopic database has been tested against aggregate total delayed neutron yields, time-dependent group parameters in 6-and 8-group re-presentation, and aggregate delayed neutron spectra. New recommendations of macroscopic delayed-neutron data for fissile materials of interest to applications are also presented. The new Reference Database for Beta-Delayed Neutron Emission Data is available online.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

237 Np Fission Spectrum Cumulative Fission Product Yield Measurement Using Godiva IV Critical Assembly

Precise integral measurement of fast neutron-induced fission product yields for various actinides is of high interest for applied nuclear science. The goal of this effort is to improve uncertainties in fission product yield values of 237 Np. Fission was induced in a NpO 2 (NO 3 ) target using the Godiva IV critical assembly in burst mode. The irradiated sample was transferred to a high-resolution γ-ray detection setup within 50 minutes. γ-ray list mode data was collected from 50 minutes to 1 week after the irradiation. γ-ray spectroscopy was performed to analyze the time dependent γ-ray yields using an automated peak search algorithm to identify isotopes by their decay γ-ray energy and half-life. Finally, the initial activity for each isotope identified was used to calculate their fission product yield.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗