Crystal plasticity modeling of irradiation effects on flow stress in pure-iron and iron-copper alloys
Not Available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
Kinetic models based on Fermi's Golden Rule are commonly employed to understand photoinduced electron transfer dynamics at molecule-semiconductor interfaces. Implicit in such second-order perturbative descriptions is the assumption that nuclear relaxation of the photoexcited electron donor is fast compared to electron injection into the semiconductor. This approximation breaks down in systems where electron transfer transitions occur on 100-fs time scale. Here, we present a fourth-order perturbative model that captures the interplay between time-coincident electron transfer and nuclear relaxation processes initiated by light absorption. The model consists of a fairly small number of parameters, which can be derived from standard spectroscopic measurements (e.g., linear absorbance, fluorescence) and/or first-principles electronic structure calculations. Insights provided by the model are illustrated for a two-level donor molecule coupled to both (i) a single acceptor level and (ii) a density of states (DOS) calculated for TiO{sub 2} using a first-principles electronic structure theory. These numerical calculations show that second-order kinetic theories fail to capture basic physical effects when the DOS exhibits narrow maxima near the energy of the molecular excited state. Overall, we conclude that the present fourth-order rate formula constitutes a rigorous and intuitive framework for understanding photoinduced electron transfer dynamics that occur on the 100-fs time scale.
For the High Luminosity Upgrade of the CERN Large Hadron Collider, lower * quadrupole magnets based on advanced Nb3Sn conductors will be installed on each side of the ATLAS and CMS experiment insertion zones. As part of the technological developments needed to achieve the required field gradient of 132.6 T/m within a 150-mm aperture, short length model magnets, named MQXFS, are tested both at the CERN SM18 and FERMI Lab test facilities. The model magnets rely on two types of Nb3Sn conductors (RRP and PIT) and on an innovative bladders and keys design to provide mechanical support against the Lorentz forces. In 2016 and 2017, the powering tests of the first two models MQXFS3 (RRP) and MQXFS5 (PIT) proved that nominal performance (16.5 kA) could be reached with excellent memory of the quench current after thermal cycle. However both magnets showed a slow training behavior with clear observations of voltage disturbances before the quench. Besides MQXFS5 only could reach ultimate current (17.9 kA) whereas erratic behavior was observed on MQXFS3 due to conductor local degradation at the head of one of the coil. In 2018, this limiting coil was changed and the applied azimuthal pre-stress increased. If ultimate current could then be reached, no stable current could be maintained due to identified defect on the outer layer of the new coil. Finally the outcome of the test of the new model MQXFS4, featuring the final RRP conductors that will be used for the series production and variation on the inner layer quench heater designs are here reported in details.
High-speed aerospace engineering applications rely heavily on computational fluid dynamics (CFD) models for design and analysis due to the expense and difficulty of flight tests and experiments. This reliance on CFD models necessitates performing accurate and reliable uncertainty quantification (UQ) of the CFD models. However, it is very computationally expensive to run CFD for hypersonic flows due to the fine grid resolution required to capture the strong shocks and large gradients that are typically present. Furthermore, UQ approaches are “many-query” problems requiring many runs with a wide range of input parameters. One way to enable computationally expensive models to be used in such many-query problems is to employ projection-based reduced-order models (ROMs) in lieu of the (high-fidelity) full-order model. In particular, the least-squares Petrov–Galerkin (LSPG) ROM (equipped with hyper-reduction) has demonstrated the ability to significantly reduce simulation costs while retaining high levels of accuracy on a range of problems including subsonic CFD applications. This allows computationally inexpensive LSPG ROM simulations to replace the full-order model simulations in UQ studies, which makes this many-query task tractable, even for large-scale CFD models. This work presents the first application of LSPG to a hypersonic CFD application. In particular, we present results for LSPG ROMs of the HIFiRE-1 in a three-dimensional, turbulent Mach 7.1 flow, showcasing the ability of the ROM to significantly reduce computational costs while maintaining high levels of accuracy in computed quantities of interest.
An investigation was conducted on an 0.08-scale semispan model of the Chance Vought XF7U-1 airplane in the Langley high-speed 7- by 10-foot tunnel in the Mach number range from 0.40 to 0.97. The results are compared with those obtained with an 0.08-scale sting-mounted complete model tested in the same tunnel and with an 0.026-scale semispan model tested by the wing-flow method. The lift-curve slopes obtained for the 0.08-scale semispan model and the 0.026-scale wing-flow model were in good agreement but both were generally lower than the values obtained for the sting model. The results of an unpublished investigation have shown that tunnel-wall boundary-layer and strut-leakage effects can came the difference noted between the lift-curve slopes of the sting and the semispan data. Fair agreement was obtained among the data of the three models as regard the variation of pitching-moment coefficients with lift coefficient. The agreement between the complete and the semispan models was more favorable with the vertical fine on, because the wall-boundary-layer and strut leakage effects were less severe. In the Mach number range between 0.94 and 0.97, ailavator-control reversal was indicated in the wing-flow data near zero lift; Whereas, these same trends were indicated in the larger scale semispan data at somewhat higher lift coefficients.
Not Available
Brown carbon (BrC) is an absorbing organic aerosol (OA), primarily emitted through biomass burning (BB), which exhibits light absorption unique to both black carbon (BC) and other organic aerosols. Despite many field and laboratory studies seeking to constrain BrC properties, the radiative forcing (RF) of BrC is still highly uncertain. To better understand its climate impact, we introduced BrC to the One-Moment Aerosol (OMA) module of the GISS ModelE Earth system model (ESM). We assessed ModelE sensitivity to primary BrC processed through a novel chemical aging scheme and to secondary BrC formed from biogenic volatile organic compounds (BVOCs). Initial results show that BrC typically contributes a top-of-the-atmosphere (TOA) radiative effect of 0.04 W m −2 . Sensitivity tests indicate that explicitly simulating BrC (separating it from other OAs), including secondary BrC, and simulating chemical bleaching of BrC contribute distinguishable radiative effects and should be accounted for in BrC schemes. This addition of prognostic BrC to ModelE allows greater physical and chemical complexity in OA representation with no apparent trade-off in model performance, as the evaluation of ModelE aerosol optical depth against Aerosol Robotic Network (AERONET) and Moderate Resolution Imaging Spectroradiometer (MODIS) retrieval data, with and without the BrC scheme, reveals similar skill in both cases. Thus, BrC should be explicitly simulated to allow more physically based chemical composition, which is crucial for more detailed OA studies like comparisons to in situ measurement campaigns. We include a summary of best practices for BrC representation within ModelE at the end of this paper.
Use of hydrogen as a fuel by [FeFe]-hydrogenase enzymes in nature requires heterolytic cleavage of the H-H bond into a proton (H+) and hydride (H-), a reaction that is also a critical step in homogeneous catalysts for hydrogenation of C=O and C=N bonds. An understanding of the catalytic oxidation of H2 by hydrogenases provides insights into the design of synthetic catalysts that are sought as cost-effective alternatives to the use of the precious metal platinum in fuel cells. Crystallographic studies on the [FeFe]-hydrogenase enzyme were critical to understanding of its reactivity, but the key H-H cleavage step is not readily observed experimentally in natural hydrogenases. Synthetic biomimics have provided evidence for H2 cleavage leading to hydride transfer to the metal and proton transfer to an amine. Limitations on the precise location of hydrogen atoms by x-ray diffraction can be overcome by use of neutron diffraction, though its use is severely limited by the difficulty of obtaining suitable crystals and by the scarcity of neutron sources. Here we show that an iron complex with a pendant amine in the diphosphine ligand cleaves hydrogen heterolytically under mild conditions, leading to [CpC5F4NFeH(PtBu2NtBu2H)]+BArF4-, [PtBu2NtBu2 = 1,5-di(tert-butyl)-3,7-di(tert-butyl)-1,5-diaza-3,7-diphosphacyclooctane; ArF = 3,5-bis(trifluoromethyl)phenyl]. The Fe-H- - - H-N moiety has a strong dihydrogen bond, with a remarkably short H • • • H distance of 1.489(10) Å between the protic N-Hδ+ and hydridic Fe-Hδ-. The structural data for [CpC5F4NFeH(PtBu2NtBu2H)]+ provide a glimpse of how the H-H bond is oxidized or generated in hydrogenase enzymes, with the pendant amine playing a key role as a proton relay. The iron complex [CpC5F4NFeH(PtBu2NtBu2H)]+BArF4- is an electrocatalyst for oxidation of H2 (1 atm) at 22 °C, so the structural data are obtained on a complex that is a functional model for catalysis by [FeFe]-hydrogenase enzymes. This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences. Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy.
Not Available
Accurate reconstruction of spatial positions of field evaporated ions from detector patterns depends on a correct understanding of the tip shape evolution and evaporation laws of component atoms. Reconstruction artifacts introduced in the heterogeneous tip can be contributed by the assumption on homogeneous evaporation behavior. A level set method based specimen shape evolution model is developed in this study to simulate the evaporation of synthetic layered-structured atom probe tip. The simulation results of the shape evolution by the present level set model agree well with the literature data using other numerical approaches, i.e. finite difference method. The asymmetric evolving shape predicted by the level set model demonstrates that the complex evaporation behavior of heterogeneous tips and interface curvature can potentially lead to the tip reconstruction distortion. Compared with other shape evaporation methods the new method provides considerably higher resolution and computational efficiency with the aid of the intrinsic subgrid scale accuracy.
This Modeling Archive is in support of an NGEE Arctic discussion paper available at doi:10.5194/gmd-9-927-2016. The purpose is to document the simulations to allow verification, reproducibility, and follow-up studies. This dataset contains shell scripts to create the CLM-PFLOTRAN cases, specific input files for PFLOTRAN and CLM, outputs, and python scripts to make the figures using the outputs in the publication. Through these results, we demonstrate that CLM-PFLOTRAN can approximately reproduce CLM results in selected cases for the Arctic, temperate and tropic sites. In addition, the new framework facilitates mechanistic representations of soil biogeochemistry processes in the land surface model. This dataset includes four zipped files and one .pdf file.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
We present a diagnostic analysis of the marine low cloud climatology simulated by two state-of-the-art coupled atmosphere-ocean models: the NCAR Community Earth System Model (CESM) and the NCEP Global Forecasting System (GFS). In both models, the shallow convection and boundary layer turbulence parameterizations have been recently updated: both models now use a mass-flux scheme for the parameterization of shallow convection, and a turbulence parameterization capable of handling Stratocumulus (Sc)-topped Planetary Boundary Layers (PBLs). For shallow convection, both models employ a convective trigger function based on the concept of convective inhibition and both include explicit convective overshooting/penetrative entrainment formulation. For Sc-topped PBL, both models treat explicitly turbulence mixing and cloud-top entrainment driven by cloud-top radiative cooling. Our focus is on the climatological transition from Sc to shallow Cumulus (Cu)-topped PBL in the subtropical eastern oceans. We show that in the CESM the coastal Sc-topped PBLs in the subtropical Eastern Pacific are well-simulated but the climatological transition from Sc to shallow Cu is too abrupt and happens too close to the coast. By contrast, in the GFS coupled simulation the coastal Sc amount and PBL depth are severely underestimated while the transition from Sc to shallow Cu is ³delayed² and offshore Sc cover is too extensive in the subtropical Eastern Pacific. We discuss the possible connections between such differences in the simulations and differences in the parameterizations of shallow convection and boundary layer turbulence in the two models.
Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.
Abstract not provided
In the first year of this contractual effort a hypo-elastic constitutive model was developed and shown to have great potential in modeling the elastic response of parenchyma. This model resides at the macroscopic level of the continuum. In this, the second year of our support, an isotropic dodecahedron is employed as an alveolar model. This is a microscopic model for parenchyma. A hopeful outcome is that the linkage between these two scales of modeling will be a source of insight and inspiration that will aid us in the final year's activity: creating a viscoelastic model for parenchyma.
A low-level model-to-model interface is presented that will enable independent models to be linked into an integrated system of models. The interface is based on a standard set of functions that contain appropriate export and import schemas that enable models to be linked with no changes to the models themselves. These ideas are presented in the context of a specific multiscale material problem that couples atomistic-based molecular dynamics calculations to continuum calculations of fluid ow. These simulations will be used to examine the influence of interactions of the fluid with an adjacent solid on the fluid ow. The interface will also be examined by adding it to an already existing modeling code, Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) and comparing it with our own molecular dynamics code.
A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.