Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 2”

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 397 records · Page 22

Chemical timescale effects on detonation convergence

Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.

Adaptive mesh refinement↗

Full-dimensional diabatic potential energy surfaces including dissociation: The {sup 2}E{sup ″} state of NO{sub 3}

A scheme to produce accurate full-dimensional coupled diabatic potential energy surfaces including dissociative regions and suitable for dynamical calculations is proposed. The scheme is successfully applied to model the two-sheeted surface of the {sup 2}E{sup ″} state of the NO{sub 3} radical. An accurate potential energy surface for the NO{sub 3}{sup −} anion ground state is developed as well. Both surfaces are based on high-level ab initio calculations. The model consists of a diabatic potential matrix, which is expanded to higher order in terms of symmetry polynomials of symmetry coordinates. The choice of coordinates is key for the accuracy of the obtained potential energy surfaces and is discussed in detail. A second central aspect is the generation of reference data to fit the expansion coefficients of the model for which a stochastic approach is proposed. A third ingredient is a new and simple scheme to handle problematic regions of the potential energy surfaces, resulting from the massive undersampling by the reference data unavoidable for high-dimensional problems. The final analytical diabatic surfaces are used to compute the lowest vibrational levels of NO{sub 3}{sup −} and the photo-electron detachment spectrum of NO{sub 3}{sup −} leading to the neutral radical in the {sup 2}E{sup ″} state by full dimensional multi-surface wave-packet propagation for NO{sub 3} performed using the Multi-Configuration Time Dependent Hartree method. The achieved agreement of the simulations with available experimental data demonstrates the power of the proposed scheme and the high quality of the obtained potential energy surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Summary of LANL Critical Benchmark Comparison Study and Revisions for Cases Involving HEU, LEU, MIX, and Pu

This report documents results obtained for revisions made to cases involving Highly Enriched Uranium (HEU), Intermediate Enriched Uranium (IEU), a mixture of Pu and Uranium (MIX), as well as Pu cases. A previous summary of revisions for HEU an Pu cases was reported and additional investigations into four cases originally presented therein uncovered further revisions which led to better agreement with other transport codes, those cases are updated in this report. The summary of all cases reported in Reference 2 is updated in this report. In addition, a previous summary of revisions for LEU and MIX was reported, a summary of those revisions in reproduced in this report for a comprehensive summary of changes to benchmarks beginning in fiscal year 2020 to current date. The report focuses on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. Feedback was used to pinpoint review of benchmark input files and to revise them when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, there is an effort tied to this work involving collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository. LANL has a benchmark library of critical experiments from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook modeled for use with MCNP. This collection is now over 1100 benchmarks, referred to as the Whisper-1.1 library because it is used with the sensitivity/uncertainty package, Whisper, which supports nuclear criticality safety validation and is released with MCNP6.2. The collection, originally created several decades ago, is a combination of smaller collections, which has been revised and expanded, by various groups at LANL over the years. The original authors are no longer at the laboratory and little formal documentation of review and revision of these benchmarks exists today. A branch of the benchmark collection was already the subject of a formal review undertaken by the LANL NCS Division and expanded to include XCP Division.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Chemical bonding in americium oxides probed by X-ray spectroscopy

The electronic structure and the chemical state in Am binary oxides and Am-doped UO 2 were studied by means of X-ray absorption spectroscopy at shallow Am core (4d and 5d) edges. In particular, the Am 5f states were probed and the nature of their bonding to the oxygen states was analyzed. The interpretation of the experimental data was supported by the Anderson impurity model (AIM) calculations which took into account the full multiplet structure due to the interaction between 5f electrons as well as the interaction with the core hole. The sensitivity of the branching ratio of the Am 4d 3/2 and 4d 5/2 X-ray absorption lines to the chemical state of Am was shown using Am binary oxides as reference systems. The observed ratio for Am-doped UO 2 suggests that at least at low Am concentrations, americium is in the Am(III) state in the UO 2 lattice. To confirm the validity of the applied AIM approach, the analysis of the Am 4f X-ray photoelectron spectra of AmO 2 and Am 2 O 3 was also performed which revealed a good agreement between experiment and calculations. As a whole, AmO 2 can be classified as the charge-transfer compound with the 5f occupancy (n ƒ ) equal to 5.73 electrons, while Am 2 O 3 is rather a Mott–Hubbard system with n ƒ = 6.05.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal Safety Analyses for the Production of Plutonium-238 at the High Flux Isotope Reactor

There has been a considerable effort over the previous few years to demonstrate and optimize the production of plutonium-238 ( 238 Pu) at the High Flux Isotope Reactor (HFIR). This effort has involved resources from multiple divisions and facilities at the Oak Ridge National Laboratory (ORNL) to demonstrate the fabrication, irradiation, and chemical processing of targets containing neptunium-237 ( 237 Np) dioxide (NpO 2 )/aluminum (Al) cermet pellets. A critical preliminary step to irradiation at the HFIR is to demonstrate the safety of the target under irradiation via documented experiment safety analyses. The steady-state thermal safety analyses of the target are simulated in a finite element model with the COMSOL Multiphysics code that determines, among other crucial parameters, the limiting maximum temperature in the target. Safety analysis efforts for this model discussed in the present report include: (1) initial modeling of single and reduced-length pellet capsules in order to generate an experimental knowledge base that incorporate initial non-linear contact heat transfer and fission gas equations, (2) modeling efforts for prototypical designs of partially loaded and fully loaded targets using limited available knowledge of fabrication and irradiation characteristics, and (3) the most recent and comprehensive modeling effort of a fully coupled thermo-mechanical approach over the entire fully loaded target domain incorporating burn-up dependent irradiation behavior and measured target and pellet properties, hereafter referred to as the production model. These models are used to conservatively determine several important steady-state parameters including target stresses and temperatures, the limiting condition of which is the maximum temperature with respect to the melting point. The single pellet model results provide a basis for the safety of the irradiations, followed by parametric analyses in the initial prototypical designs that were necessary due to the limiting fabrication and irradiation data available. The calculated parameters in the final production target model are the most accurate and comprehensive, while still conservative. Over 210 permutations in irradiation time and position were evaluated, and are supported by the most recent inputs and highest fidelity methodology. The results of these analyses show that the models presented in this report provide a robust and reliable basis for previous, current and future experiment safety analyses. In addition, they reveal an evolving knowledge of the steady-state behavior of the NpO 2 /Al pellets under irradiation for a variety of target encapsulations and potential conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Direct observation of C 3 S particle dissolution using fast nano X-ray computed tomography

Tricalcium silicate (C 3 S) occupies 50 % to 70 % of ordinary portland cement (OPC) by mass and it is an important component affecting the hydration of OPC [1], [2], [3], [4], [5], [6], [7]. Generally, the hydration of C 3 S is described by two processes: the dissolution of C 3 S particles and the precipitation of hydration products. While it is understood that the dissolution rates of C 3 S vary with time, more precise measurements are needed to understand this process. Many mechanisms have been proposed to explain the time-evolving dissolution rates of C 3 S [2]. The metastable barrier hypothesis suggests that a thin metastable layer of hydrates forms around the C 3 S particle surface and prohibits the access of grains to the aqueous solution [8], [9], [10], [11], [12]. The slow dissolution step hypothesis suggests that the increased ion concentration from the initial reaction delays the C 3 S dissolution [2], [13], [14], [15], [16], [17]. More recent publications suggest that C 3 S may react differently depending on the existence of crystallographic defects [18], [19], [20]. Etch pits are thought to open on the particle surface during the initial reaction and this contributes to the C 3 S dissolution [21], [22]. As hydrates precipitate and cover these highly reactive surfaces, hydration slows down and the induction period starts [18], [23], [24], [25], [26]. Many experiments have been conducted to test the aforementioned mechanisms. Some hydration studies utilize bulk measurements, such as isothermal calorimetry [27], [28], [29], pore solution analysis [30], quasi-elastic neutron scattering [31], and nuclear magnetic resonance spectroscopy [32], [33]. One limitation of these measurements is that they do not provide direct and detailed information on the individual C 3 S particles. Some other studies utilize imaging techniques, such as scanning electron microscopy (SEM) [34], [35], [36] and transmission electron microscopy (TEM) [37]. However, SEM/TEM cannot track the evolution of individual particles throughout hydration [34], [35], [38], [39], [40] and they do not give insights into the microstructure of materials before hydration [34], [35], [39]. This makes it challenging to draw strong conclusions from only SEM or TEM observations. Synchrotron X-ray tomography techniques have been used more broadly in recent years to study cement hydration. They are not only non-destructive but also able to image a sample in full 3D with resolutions that can reach from micron to nanoscale. Nano computed tomography (nCT) is one technique that has been applied to study cement hydration at the nanoscale [26], [41]. A typical nCT can reach a pixel size from 15 to 65 nm, providing enough detail for observing features <1 μm. However, nCT often takes >0.5 h to finish one scan. This makes the application of this technique on continuous scans for in-situ observations challenging. Fast X-ray computerized tomography (fCT) is another technique that has shown success in studying the time-evolving cement microstructures [20], [42], [43], [44], [45], [46], [47], [48], [49]. Due to the high flux of the X-ray beam from the synchrotron ring, fCT allows a scan to be captured within 1 min at a pixel size of 1 μm. This allows a paste sample to be continuously scanned during the hydration process. However, the micron-sized resolutions limit does not provide detailed insights for particles <5 μm [20], [49]. Fortunately, the combination of nCT and fCT has allowed the development of fast nano X-ray computed tomography (fnCT). fnCT can capture a 3D data set in <2 min at a pixel size of 50 nm. This makes this procedure an exciting method to evaluate hydrating pastes. fnCT collects multiple X-ray radiographs at various rotation angles and generates a 3D model of the scanned sample, which is also referred to as a 3D tomography [50], [51]. In one tomography, the X-ray absorptions of different components (e.g., C 3 S and hydrates) differ as functions of density and chemistry [52], [53]. These X-ray absorption contrasts can be used to extract detailed information about the 3D microstructure [26], [54], [55]. In this paper, fnCT is used to collect time-lapse tomographs of hydrating C 3 S paste from 18 min after mixing to 7 h of hydration. The bulk measurements of anhydrous C 3 S, as well as the microstructural changes of individual C 3 S particles, are directly observed, quantified, and discussed. The dissolution behavior of C 3 S particles at various size scales is systematically analyzed and compared. This work aims to find the relationship between the size of C 3 S particle sizes and their dissolution rates. This provides significant insights into the early-age hydration of C 3 S on length and time scales not previously possible. Because of the magnitude of the data and the substantial amount of observations, this work will solely focus on the change in the anhydrous particles. Changes in the hydration products will be reported in future work.

42 ENGINEERING↗

Evidence of medium response to hard probes using correlations of Z bosons with hadrons in heavy ion collisions

The first measurement of pseudorapidity and azimuthal angle distributions relative to the momentum vector of a Z boson for low transverse momentum ($p_\mathrm{T}$) charged hadrons in lead-lead (PbPb) collisions is presented. By studying the hadrons produced in an event with a high-$p_\mathrm{T}$ Z boson (40 $\lt$ $p_\mathrm{T}$ $\lt$ 350 GeV), the analysis probes how the quark-gluon plasma (QGP) medium created in these collisions affects the parton recoiling opposite to the Z boson. Utilizing PbPb data at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_\mathrm{NN}}}$ = 5.02 TeV from 2018 with an integrated luminosity of 1.67 nb$^{-1}$ and proton-proton (pp) data at the same energy from 2017 with 301 pb$^{-1}$, the distributions are examined in bins of charged-hadron $p_\mathrm{T}$. A significant modification of the distributions for charged hadrons in the range 1 $\lt$ $p_\mathrm{T}$ $\lt$ 2 GeV in PbPb collisions is observed when compared to reference measurements from pp collisions. The data provide new information about the correlation between hard and soft particles in heavy ion collisions, which can be used to test predictions of various jet quenching models. The results are consistent with expectations of a hydrodynamic wake created when the QGP is depleted of energy by the parton propagating through it. Based on comparisons of PbPb data with pp references and predictions from theoretical models, this Letter presents the first evidence of medium-recoil and medium-hole effects caused by a hard probe.

CMS↗

A Comparison of Exact and Model Exchange–Correlation Potentials for Molecules

Accurate exchange–correlation (XC) potentials for three-dimensional systems–via solution of the inverse density functional theory (DFT) problem–are now available to test the quality of DFT approximations. Herein, the exact XC potential for seven molecules–dihydrogen at four different bond-lengths, lithium hydride, water, and ortho-benzyne–are computed from full configuration interaction reference densities. These are compared to model XC potentials from nonlocal (B3LYP, HSE06, SCAN0, and M08-HX) and semilocal/local (SCAN, PBE, and PW92) XC functionals. Whereas for most systems, relative errors in the ground-state densities are O(10 –3 –10 –2 ) , the model XC potentials have much higher errors of O(10 –1 –10 0 ) . Among the model XC functionals, SCAN0 offers the best agreement with the exact XC potential, underlining the significance of satisfying exact conditions as well as including nonlocal effects in XC functionals. This work indicates that tests against the exact XC potential will provide a promising new direction for building more accurate XC functionals for DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-F-T-31: Shape and Isodose Distributions in Co60 HDR Brachytherapy for Different Utero-Vaginal Time Ratios

Purpose: To optimize the dose in bladder and rectum and show the different shapes of the isodose volumes in Co60-HDR brachytherapy, considering different utero and vaginal sources dwell ratio times (TU:TV). Methods: Besides Ir192-HDR, new Co60-HDR sources are being incorporated. We considered different TU:TV times and computed the dosis in bladder, rectum and at the reference points of the Manchester system. Also, we calculated the isodose volume and shape in each case. We used a EZAG-BEBIG Co0.A86 model with TPS HDRplus3.0.4. and LCT42-7, LCT42-2(R,L) applicators. A reference dose RA= 1.00 Gy was given to the A-right point. We considered the TU:TV dwell time ratios 1:0.25, 1:0.33, 1:0.5, 1:1, 1:2, 1:3, and 1:4. Given TU:TV, the stop time at each dwell position is fixed for each applicator. Results: Increasing TU:TV systematically results in a decreasing of the dose in bladder and rectum, e.g. 9% and 27% reduction were found in 1:0.25 with respect to 1:1, while 12% and 34% increase were found in 1:4 with respect to 1:1. Also, the isodose volume parameters height (h), width (w), thickness (t) and volume (hwt) increased from the 1:0.25 case to the 1:4 value: hwt is 25% lower and 31% higher than the 1:1 reference volume in these cases. Also w decreased for higher TU:TV and may compromise the tumoral volume coverage, decreasing 17% in the 1:0.25 case compared to the 1:1 case. The shape of the isodose volume was obtained for the different TU:TV considered. Conclusion: We obtained the shape of isodose volumes for different TU:TV values in gynecological Co60-HDR. We studied the dose reduction in bladder and rectum for different TU:TV ratios. The volume parameters and hwt are strongly dependent on this ratio. This information is useful for a quantitative check of the TPS and as a starting point towards optimization.

60 APPLIED LIFE SCIENCES↗

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model (Wagner and Pruss, 2002) for the thermodynamic properties of water. It further evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations (Cox et al., 1989). The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature and reduced density. Here is the absolute temperature (K), is density (kg/m3), and the subscript refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which is 322 kg/m3). The code solves four basic types of problems, distinguished by the specified inputs: 1, Temperature (K) and density (kg/m3) or reduced density 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, only some of which correspond to valid results. Obtaining valid results depends mainly on appropriate choice of starting values. The default values in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 (Johnson et al., 1992) and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa).

Wolery, ThomasJ.↗

Following O and OH in He/O 2 and He/H 2 O gas mixtures—from the gas phase through the liquid phase to modifications on a biological sample

Abstract Applied cold atmospheric plasma allows for the controlled delivery of reactive oxygen and nitrogen species tailored for specific applications. Through the manipulation of the plasma parameters, feed gases, and careful consideration of the environment surrounding the treatment target, selective chemistries that preferentially influence the target can be produced and delivered. To demonstrate this, the COST reference microscale atmospheric pressure plasma jet is used to study the generation and transport of O and ⋅ OH from the gas phase through the liquid to the biological model target cysteine. Relative and absolute species densities of ⋅ OH and O are measured in the gas phase through laser induced fluorescence (LIF) and two-photon absorption LIF respectively. The transport of these species is followed into the liquid phase by hydrogen peroxide quantification and visualized by a fluorescence assay. Modifications to the model biological sample cysteine exposed to ⋅ OH and H 2 O 2 dominated chemistry (He/H 2 O (0.25%)) and O dominated chemistry (He/O 2 (0.6%)) is measured by FTIR spectroscopy. The origin of these species that modify cysteine is considered through the use of heavy water (H 2 18 O) and mass spectrometry. It is found that the reaction pathways differ significantly for He/O 2 and He/H 2 O. Hydrogen peroxide is formed mainly in the liquid phase in the presence of a substrate for He/O 2 whereas for He/H 2 O it forms in the gas phase. The liquid chemistry resulting from the He/O 2 admixture mainly targets the sulfur moiety of cysteine for oxidation up to irreversible oxidation states, while He/H 2 O treatment leads preferentially to reversible oxidation products. The more O or OH/H 2 O 2 dominated chemistry produced by the two gas admixtures studied offers the possibility to select species for target modification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resource-Efficient Chemistry on Quantum Computers with the Variational Quantum Eigensolver and The Double Unitary Coupled-Cluster approach

Applications of quantum simulation algorithms to obtain electronic energies of molecules on noisy intermediate-scale quantum (NISQ) devices require careful consideration of resources describing the complex electron correlation effects. In modeling second-quantized problems, the biggest challenge confronted is that the number of qubits scales linearly with the size of molecular basis. This poses a significant limitation on the size of the basis sets and the number of correlated electrons included in quantum simulations of chemical processes. To address this issue and to enable more realistic simulations on NISQ computers, we employ the double unitary coupled-cluster (DUCC) method to effectively downfold correlation effects into the reduced-size orbital space, commonly referred to as the active space. Using downfolding techniques, we demonstrate that properly constructed effective Hamiltonians can capture the effect of the whole orbital space in small-size active spaces. Combining the downfolding pre-processing technique with the Variational Quantum Eigensolver, we solve for the ground-state energy of H2 and Li2 in the cc-pVTZ basis using the DUCC-reduced active spaces. We compare these results to full configuration-interaction and high-level coupled-cluster reference calculations.

quantum computing, variational quantum solver, cou↗

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗

Impact of chain microstructure on solution and thin film self-assembly of PCHD-based semi-flexible/flexible diblock copolymers

We examine self-assembly of semi-flexible/flexible block copolymers in a selective solvent iusing a set of diblock copolymers where the chain microstructure of the semi-flexible block is manipulated in order to tune chain stiffness. The reduced conformational space of the semi-flexible block is anticipated to alter the way the chains pack, potentially changing the structure of self-assembled aggregates in comparison to flexible diblock copolymer analogs. Semi-flexible/flexible diblock copolymers comprised of poly(styrene)-block-poly(1,3-cyclohexadiene) (PS-b-PCHD) having systematic changes in chain microstructure, as captured by the ratio of 1,4/1,2-linkages between cyclohexenyl repeat units, and molecular weight of the PCHD blocks were synthesized using anionic polymerization. Moreover, these diblocks were dissolved in tetrahydrofuran (THF), which is a preferential solvent for PS, and the structures formed were examined using laser light scattering and complementary imaging techniques. Results show that PS-b-PCHD copolymers with a chain microstructure of 90% 1,4/10% 1,2 linkages between cyclohexenyl repeat units (referred to simply as 90/10) are able to micellize, forming spherical structures, while diblocks of 70/30 and 50/50 1,4-to-1,2 ratios remain as single chains and ill-defined aggregates, respectively, when dissolved in THF. With inferences drawn from simple structural models, we speculate that this self-assembly behavior arises due to the change in the chain configuration with increasing content of 1,2-links in the backbone. Finally, this renders the chain with higher 1,2 content incapable of swelling in response to solvent and unable to pack into well-defined self-assembled structures.

36 MATERIALS SCIENCE↗

The structure of PX{sub 3} (X = Cl, Br, I) molecular liquids from X-ray diffraction, molecular dynamics simulations, and reverse Monte Carlo modeling

Synchrotron X-ray diffraction measurements have been conducted on liquid phosphorus trichloride, tribromide, and triiodide. Molecular Dynamics simulations for these molecular liquids were performed with a dual purpose: (1) to establish whether existing intermolecular potential functions can provide a picture that is consistent with diffraction data and (2) to generate reliable starting configurations for subsequent Reverse Monte Carlo modelling. Structural models (i.e., sets of coordinates of thousands of atoms) that were fully consistent with experimental diffraction information, within errors, have been prepared by means of the Reverse Monte Carlo method. Comparison with reference systems, generated by hard sphere-like Monte Carlo simulations, was also carried out to demonstrate the extent to which simple space filling effects determine the structure of the liquids (and thus, also estimating the information content of measured data). Total scattering structure factors, partial radial distribution functions and orientational correlations as a function of distances between the molecular centres have been calculated from the models. In general, more or less antiparallel arrangements of the primary molecular axes that are found to be the most favourable orientation of two neighbouring molecules. In liquid PBr{sub 3} electrostatic interactions seem to play a more important role in determining intermolecular correlations than in the other two liquids; molecular arrangements in both PCl{sub 3} and PI{sub 3} are largely driven by steric effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CORPSE model with litter decomposition parameters derived from the LIDET dataset

This is a version of the CORPSE model (Carbon, Organisms, Rhizosphere and Protection in the Soil Environment, Sulman et al. 2014) that uses litter decomposition parameters derived from a modified Monte Carlo simulation using the LIDET litter decomposition dataset (Long-term Intersite Decomposition Experiment Team, Harmon 2013). The code also includes the Baseline parameters, and the eight other best parameter sets identified in a modified Monte Carlo simulation. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of files:The folder "Input Files" contains one folder for each LIDET site with data necessary to run the model. Note that "(site)" in the filenames below indicates where the LIDET site code appears (see Table 1 for site codes). Data streams include: CORPSE_full_spinup_litter.csv, CORPSE_full_spinup_rhizo.csv, CORPSE_full_spinup_bulk.csv, litterbag_init_100g_6spp.csv: initial C and N (kg C or N/m2) pool values for each soil layer, the litterbag_init_100_6spp.csv file is for the litterbag layer and is the same file for all sites. All initial C and N files have the same columns (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 soilT (site) DOY274start.csv: Average daily soil temperature (oC) interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). soilT (site) DOY274start.csv: Average daily soil volumetric water content (VWC) scalar interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). litter production.csv: Average daily litter production values for each site, data sources listed in Table S3 of related publication. litter (site) CN.csv: C:N ratio for each species from LIDET dataset (Table 2, Harmon 2013). (site).csv: Table indicating number of observations for each species decomposed at each site. Instructions: Save the model code ("CORPSE_LIDET.R") and "Input Files" folder in the same folder. Also make a folder for the model output (e.g., "results_Baseline") in the same folder. Set the working directory (setwd) in the model code to the folder with the files saved in step #1. Select the parameter set to use for the litter and litterbag compartments, comment out all other parameter sets. Run code. Output will be saved in the folder made in step 1. Output destination can be changed as necessary in code section called "Running the model." Table 1 LIDET sites and site codes used in model files. Site Code - Site AND - H.J. Andrews Experimental Forest BNZ - Bonanza Creek Experimental Forest BSF - Blodgett Research Forest CDR - Cedar Creek Natural History Area CPR - Central Plains Experimental Range HBR - Hubbard Brook Experimental Forest HFR - Harvard Forest JUN - Juneau KBS - Kellogg Biological Station KNZ - Konza Prairie Research Natural Area NWT - Niwot Ridge/Green Lakes Valley OLY - Olympic National Park OLY Conifer forest SEV - Sevilleta National Wildlife Refuge SMR - Santa Margarita Ecological Reserve UFL - University of Florida VCR - Virginia Coast Reserve Table 2 LIDET species and species codes used in model files (6 common species). Species - Species Code Sugar maple (Acer saccharum) - ACSA Drypetes (Drypetes glauca) - DRGL Red pine (Pinus resinosa) - PIRE Chestnut oak (Quercus prinus) - QUPR Western redcedar (Thuja plicata) - THPL Wheat (Triticum aestivum) - TRAE References:Bonan, G. B., Hartman, M. D., Parton, W. J., & Wieder, W. R. (2013). Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4). Global Change Biology, 19(3), 957-974. https://doi.org/https://doi.org/10.1111/gcb.12031 Harmon, M. (2013). LTER Intersite Fine Litter Decomposition Experiment (LIDET), 1990 to 2002. Long-Term Ecological Research. Forest Science Data Bank, Corvallis, OR. [Data set]. Accessed http://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023. https://doi.org/10.6073/pasta/f35f56bea52d78b6a1ecf1952b4889c5. Sulman, B. N., Phillips, R. P., Oishi, A. C., Shevliakova, E., & Pacala, S. W. (2014). Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2. Nature Climate Change, 4, 1099 - 1102. https://doi.org/10.1038/nclimate2436

Juice, Stephanie↗

Diffusion Dominant Solute Transport Modelling In Deep Repository Under The Effect of Emplacement Media Degradation - 13285

Deep geologic disposal of high activity and long-lived radioactive waste is being actively considered and pursued in many countries, where low permeability geological formations are used to provide long term waste contaminant with minimum impact to the environment and risk to the biosphere. A multi-barrier approach that makes use of both engineered and natural barriers (i.e. geological formations) is often used to further enhance the containment performance of the repository. As the deep repository system subjects to a variety of thermo-hydro-chemo-mechanical (THCM) effects over its long 'operational' lifespan (e.g. 0.1 to 1.0 million years, the integrity of the barrier system will decrease over time (e.g. fracturing in rock or clay)). This is broadly referred as media degradation in the present study. This modelling study examines the effects of media degradation on diffusion dominant solute transport in fractured media that are typical of deep geological environment. In particular, reactive solute transport through fractured media is studied using a 2-D model, that considers advection and diffusion, to explore the coupled effects of kinetic and equilibrium chemical processes, while the effects of degradation is studied using a pore network model that considers the media diffusivity and network changes. Model results are presented to demonstrate the use of a 3D pore-network model, using a novel architecture, to calculate macroscopic properties of the medium such as diffusivity, subject to pore space changes as the media degrade. Results from a reactive transport model of a representative geological waste disposal package are also presented to demonstrate the effect of media property change on the solute migration behaviour, illustrating the complex interplay between kinetic biogeochemical processes and diffusion dominant transport. The initial modelling results demonstrate the feasibility of a coupled modelling approach (using pore-network model and reactive transport model) to examine the long term behaviour of deep geological repositories with media property change under complex geochemical conditions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling the effectiveness of remediation efforts in contaminated urban areas: An EMRAS II Urban Areas Working Group exercise - 15631

The Urban Areas Working Group was part of the IAEA's EMRAS II (Environmental Modelling for Radiation Safety) Programme during 2009-2011. The Working Group's goal was to test and improve the capabilities of models used in assessment of radioactive contamination in urban settings. The 'contaminant transport and countermeasures' exercise was focused on modeling the effectiveness of potential countermeasures or remediation efforts with respect to reduction of human exposures and doses. The exercise was based on a hypothetical situation involving an initial concentration of Co-60 or Pu-239 in air in a specified urban area. Participants were given detailed input information, including descriptions of test locations in two types of urban areas, (1) a business area with buildings and asphalt, and (2) a park area near an apartment town. Several kinds of initial weather conditions (dry, light rain, and heavy rain) were explored. Modeling endpoints included contamination densities, dose rates (from each contributing surface and total), annual and cumulative doses (internal and external) for specified reference individuals, and effectiveness of countermeasures and remediation efforts in terms of dose rate reduction and dose reduction. External doses were emphasized for Co-60 and internal doses for Pu-239. Approaches and model predictions from five models are compared, and explanations for similarities and differences among model predictions are offered. Understanding the effects of different approaches, assumptions, and parameter values is essential to a comparison of results from different models and to understanding both the usefulness and the limitations of any given model. (authors)

54 ENVIRONMENTAL SCIENCES↗