Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic profiling”

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,261 records · Page 70

Suppression of Composition g -modes in Chemically Equilibrating Warm Neutron Stars

We investigate the impact of chemical equilibration and the resulting bulk viscosity on nonradial oscillation modes of warm neutron stars at temperatures up to T ≈ 5 MeV, relevant for protoneutron stars and neutron star postmerger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition g-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition g-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition g-modes. In contrast, the f-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of g-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.

Neutron stars↗

Accounting for forest carbon pool dynamics in product carbon footprints: Challenges and opportunities

Modification and loss of forests due to natural and anthropogenic disturbance contribute an estimated 20% of annual greenhouse gas (GHG) emissions worldwide. Although forest carbon pool modeling rarely suggests a 'carbon neutral' flux profile, the life cycle assessment community and associated product carbon footprint protocols have struggled to account for the GHG emissions associated with forestry, specifically, and land use generally. Principally, this is due to underdeveloped linkages between life cycle inventory (LCI) modeling for wood and forest carbon modeling for a full range of forest types and harvest practices, as well as a lack of transparency in globalized forest supply chains. In this paper, through a comparative study of U.S. and Chinese coated freesheet paper, we develop the initial foundations for a methodology that rescales IPCC methods from the national to the product level, with reference to the approaches in three international product carbon footprint protocols. Due to differences in geographic origin of the wood fiber, the results for two scenarios are highly divergent. This suggests that both wood LCI models and the protocols need further development to capture the range of spatial and temporal dimensions for supply chains (and the associated land use change and modification) for specific product systems. The paper concludes by outlining opportunities to measure and reduce uncertainty in accounting for net emissions of biogenic carbon from forestland, where timber is harvested for consumer products. - Highlights: Black-Right-Pointing-Pointer Typical life cycle assessment practice for consumer products often excludes significant land use change emissions when estimating carbon footprints. Black-Right-Pointing-Pointer The article provides a methodology to rescale IPCC guidelines for product-level carbon footprints. Black-Right-Pointing-Pointer Life cycle inventories and product carbon footprint protocols need more comprehensive land use-related accounting. Black-Right-Pointing-Pointer Interdisciplinary collaboration linking the LCA and forest carbon modeling communities is necessary.

54 ENVIRONMENTAL SCIENCES↗

Transcriptomic Profiles of Zymomonas mobilis 8b to Furfural Acute and Long-Term Stress in Both Glucose and Xylose Conditions

Zymomonas mobilis 8b is an ethanologenic bacterium engineered to utilize both glucose and xylose. The impacts of lignocellulosic hydrolyzate inhibitors on the growth of Zymomonas mobilis 8b have been investigated. However, the molecular responses of these inhibitors have not been completely elucidated yet. In this study, molecular responses to furfural were investigated using transcriptomic approaches of both chip-based microarray and a directional mRNA-Seq. Furfural acute shock time-course experiment with 3 g/L furfural supplemented when cells reached exponential phase and stress response experiment in the presence of 2 g/L furfural from the beginning of fermentation were carried out to study the physiological and transcriptional profiles of short-term and long-term effects of furfural on 8b. Furfural negatively affected 8b growth in terms of final biomass and the fermentation time. Transcriptomic studies indicated that the response of 8b to furfural was dynamic and complex, and differences existed between short-term shock and long-term stress responses. However, the gene function categories were similar with most down-regulated genes related to translation and biosynthesis, while the furfural up-regulated genes were mostly related to general stress responses. Several gene candidates have been identified and genetic studies indicated that expression of ZMO0465 and cysteine synthase operon ZMO0003-0006 driven by its native promoter in a shuttle vector enhanced the furfural tolerance of 8b. In addition, the relationship between microarray and mRNA-Seq was compared with good correlations. The directional mRNA-Seq data not only provided the gene expression profiling, but also can be applied for transcriptional architecture improvement to identify and confirm operons, novel transcripts, hypothetical gene functions, transcriptional start sites, and promoters with different strength.

09 BIOMASS FUELS↗

Dynamics of phage-host interactions in Bacteroides fragilis resolved by single-cell transcriptomics

The interactions between lytic phages and their hosts are typically studied in bulk culture, which obscures cell-cell differences in infection susceptibility or expression of protective factors. Here, we use bacterial single-cell RNA sequencing to profile the transcriptomes of ~50,000 cells from cultures of a human pathobiont, Bacteroides fragilis, infected with a lytic bacteriophage. From a single sampling, we quantified the asynchronous progression of phage infection in individual bacterial cells and reconstructed the infection timeline, characterizing both host and phage transcriptomic changes as infection unfolded. Further, we discovered phenotypic subpopulations of bacteria that remained uninfected. Each cell’s vulnerability to phage infection was influenced by expression of multiple genetic loci, most prominently phase-variable capsular polysaccharide (CPS) biosynthesis pathways and an operon predicted to encode fimbrial genes. These findings uncovered genome-wide phase variation and stochasticity that enable bacterial survival and re-growth without acquiring additional mutations. Overall, we establish bacterial single-cell RNA sequencing as a powerful platform for investigating the dynamics of host-phage interactions and revealing the roles of phase variation and stochasticity in bacterial defenses.

Bacteria↗

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

Engineering Redox Flow Battery Electrodes with Spatially Varying Porosity Using Non‐Solvent‐Induced Phase Separation

Redox flow batteries (RFBs) are a promising electrochemical platform for efficiently and reliably delivering electricity to the grid. Within the RFB, porous carbonaceous electrodes facilitate electrochemical reactions and distribute the flowing electrolyte. Tailoring electrode microstructure and surface area can improve RFB performance, lowering costs. Electrodes with spatially varying porosity may increase electrode utilization and provide surface area in reaction‐limited zones; however, the efficacy of such designs remains an open area of research. Herein, a non‐solvent‐induced phase‐separation (NIPS) technique that enables the reproducible synthesis of macrovoid‐free electrodes with well‐defined across‐thickness porosity gradients is described. The monotonically varying porosity profile is quantified and the physical properties and surface chemistries of porosity‐gradient electrodes are compared with macrovoid‐containing electrode, also synthesized by NIPS. Then, the electrochemical and fluid dynamic performance of the porosity‐gradient electrodes is evaluated, exploring the effect of changing the direction of the porosity gradient and benchmarking against the macrovoid‐containing electrode. Lastly, the performance is examined in a vanadium RFB, finding that the porosity‐gradient electrode outperforms the macrovoid electrode, is independent of gradient direction, and performs favorably compared to advanced electrodes in the contemporary literature. It is anticipated that the approach motivates further exploration of microstructurally tailored electrodes in electrochemical systems.

25 ENERGY STORAGE↗

Increases in Cell-Wall Homogalacturonan but Decreases in Xylogalacturonan Accompany Transition from Dormant to Vegetative Stages in Chrysolaena Obovata Rhizophores

The rhizophores of Chrysolaena obovata, a perennial herb native to the Brazilian Cerrado, are underground storage organs that accumulate fructans and support vegetative propagation. Here, we investigated how the content of fructans and the composition of rhizophore cell wall polysaccharides change across developmental stages and contribute to seasonal adaptation. In addition to the accumulation of fructans in vacuoles (~26-55% of the dry mass according to the developmental stage), cell-wall polysaccharide composition changed throughout the growth cycle, reflecting structural adjustments linked to developmental and environmental cues. These cell walls have a type I architecture, typical of eudicots. Although total pectin content remained stable, homogalacturonans increased in proportion during the transition from dormant to vegetative stages, whereas amounts of xylogalacturonans decreased. Xylogalacturonans accumulate within the middle lamella during the dormant phase and are potentially associated with flexibility and water retention. Hemicelluloses consisted mainly of low-substituted arabinoxylans and xyloglucan, consistent with a secondary cell wall profile. Xylans decreased slightly in the transition from the dormant to vegetative stage, while xyloglucan proportions increased. Overall, our findings show that seasonal metabolic changes in C. obovata involve dynamic modifications of rhizophore cell wall polysaccharides, which contribute to organ development and drought adaptation throughout the seasonal cycle.

59 BASIC BIOLOGICAL SCIENCES↗

Impact of resistance mutations on inhibitor binding to HIV-1 integrase

Here, HIV-1 integrase (IN) is essential for HIV-1 replication, catalyzing two key reaction steps termed 3' processing and strand transfer. Therefore, IN has become an important target for antiviral drug discovery. However, mutants have emerged, such as E92Q/N155H and G140S/Q148H, which confer resistance to raltegravir (RAL), the first IN strand transfer inhibitor (INSTI) approved by the FDA, and to the recently approved elvitegravir (EVG). To gain insights into the molecular mechanisms of ligand binding and drug resistance, we performed molecular dynamics (MD) simulations of homology models of the HIV-1 IN and four relevant mutants complexed with viral DNA and RAL. The results show that the structure and dynamics of the 140s loop, comprising residues 140 to 149, are strongly influenced by the IN mutations. In the simulation of the G140S/Q148H double mutant, we observe spontaneous dissociation of RAL from the active site, followed by an intrahelical swing-back of the 3' -OH group of nucleotide A17, consistent with the experimental observation that the G140S/Q148H mutant exhibits the highest resistance to RAL compared to other IN mutants. An important hydrogen bond between residues 145 and 148 is present in the wild-type IN but not in the G140S/Q148H mutant, accounting for the structural and dynamical differences of the 140s' loop and ultimately impairing RAL binding in the double mutant. End-point free energy calculations that broadly capture the experimentally known RAL binding profiles elucidate the contributions of the 140s' loop to RAL binding free energies and suggest possible approaches to overcoming drug resistance.

60 APPLIED LIFE SCIENCES↗

Developing a Redox Network for Coastal Saltmarsh Systems in the PFLOTRAN Reaction Model

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Here our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO 2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO 2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO 2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.

54 ENVIRONMENTAL SCIENCES↗

The radial profile and flattening of the Milky Way’s stellar halo to 80 kpc from the SEGUE K-giant survey

We characterize the radial density, metallicity, and flattening profile of the Milky Way’s stellar halo, based on the large sample of spectroscopically confirmed giant stars from SDSS/SEGUE-2, spanning galactocentric radii 10 kpc ≤ r{sub GC} ≤ 80 kpc. After excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream), the sample has 1757 K giants, with a typical metallicity precision of 0.2 dex and a mean distance accuracy of 16%. Compared to blue horizontal branch stars or RR Lyrae variables, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modeling of those data, based on ellipsoidal stellar density models, ν{sub *}(R, z), with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of n = 3.1 ± 0.5 with an effective radius r{sub eff}=15±2 kpc and a flattening of q = 0.7 ± 0.02, or comparably well by an equally flattened broken power law, with radial slopes of α{sub in} = 2.1 ± 0.3 and α{sub out} = 3.8 ± 0.1, with a break radius of r{sub break} = 18 ± 1 kpc; this is largely consistent with earlier work. We find a modest but significant metallicity gradient within the “outer” stellar halo, [Fe/H] decreasing outward. If we allow for a variable flattening q=f(r{sub GC}), we find the distribution of halo giants to be considerably more flattened at small radii, q(10 kpc) = 0.55 ± 0.02, compared to q(>30 kpc) = 0.8 ± 0.03. Remarkably, the data are then very well fit by a single power law with index of 4.2 ± 0.1 on the variable r{sub q}≡√(R{sup 2}+(z/q(r)){sup 2}). In this simple and better-fitting model, there is a break in flattening at ∼20 kpc, instead of a break in the radial density function. While different parameterizations of the radial profile vary in their parameters, their implied density gradient, ∂lnν{sub ∗}/∂lnr, is stable along a direction intermediate between major and minor axis; this gradient is crucial in any dynamical modeling that uses halo stars as tracers.

79 ASTRONOMY AND ASTROPHYSICS↗

The Fermi-Pasta-Ulam recurrence and related phenomena for 1D shallow-water waves in a finite basin

Different regimes of the Fermi-Pasta-Ulam (FPU) recurrence are simulated numerically for fully nonlinear 'one-dimensional' potential water waves in a finite-depth flume between two vertical walls. In such systems, the FPU recurrence is closely related to the dynamics of coherent structures approximately corresponding to solitons of the integrable Boussinesq system. A simplest periodic solution of the Boussinesq model, describing a single soliton between the walls, is presented in analytic form in terms of the elliptic Jacobi functions. In the numerical experiments, it is observed that depending on the number of solitons in the flume and their parameters, the FPU recurrence can occur in a simple or complicated manner, or be practically absent. For comparison, the nonlinear dynamics of potential water waves over nonuniform beds is simulated, with initial states taken in the form of several pairs of colliding solitons. With a mild-slope bed profile, a typical phenomenon in the course of evolution is the appearance of relatively high (rogue) waves, while for random, relatively short-correlated bed profiles it is either the appearance of tall waves or the formation of sharp crests at moderate-height waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The FaceBase Consortium: A comprehensive resource for craniofacial researchers

The FaceBase Consortium, funded by the National Institute of Dental and Craniofacial Research, National Institutes of Health, is designed to accelerate understanding of craniofacial developmental biology by generating comprehensive data resources to empower the research community, exploring high-throughput technology, fostering new scientific collaborations among researchers and human/computer interactions, facilitating hypothesis-driven research, and translating science into improved health care to benefit patients. The resources generated by the FaceBase projects include a number of dynamic imaging modalities, genome-wide association studies, software tools for analyzing human facial abnormalities, detailed phenotyping, anatomical and molecular atlases, global and specific gene expression patterns, and transcriptional profiling over the course of embryonic and postnatal development in animal models and humans. The integrated data visualization tools, faceted search infrastructure, and curation provided by the FaceBase Hub offer flexible and intuitive ways to interact with these multidisciplinary data. In parallel, the datasets also offer unique opportunities for new collaborations and training for researchers coming into the field of craniofacial studies. Here we highlight the focus of each spoke project and the integration of datasets contributed by the spokes to facilitate craniofacial research.

60 APPLIED LIFE SCIENCES↗

NRAP-Open-IAM: NRAP Open Source Integrated Assessment Model

Note: This is the last version (a2.6.1) of NRAP-Open-IAM released during NRAP Phase II in 2022. The latest version of NRAP-Open-IAM is available here: https://edx.netl.doe.gov/dataset/phase-iii-nrap-open-iam NRAP-Open-IAM is an open-source software product that enables quantification of containment effectiveness and leakage risk at storage sites in the context of system uncertainties and variability. NRAP-Open-IAM represents the next-generation in a line of systems-based computational models developed for quantitative geological carbon storage (GCS) risk assessment. The model comprises a set of reduced-order and analytical models of various components of the GCS system, potential leakage pathways, receptors of concern including impact to groundwater resources and the atmosphere, a framework to support stochastic simulation, time stepping, uncertainty quantification, other analytical functionality for scenario and risk-performance evaluation, and a basic graphical user interface to support scenario development, data input simulation definition, and basic post-processing and results display. As the NRAP Open-IAM functionality continues to evolve, we continue to add to its capability to develop quantitative, probabilistic, and time-dependent profiles of the evolution of risk at a GCS site and evaluate the influence of uncertain parameters on uncertainty in predicted risk. It can be used to quantify the dynamics of reservoir saturation plume and pressure-affected area, for evaluation of the area of potential groundwater impact (i.e., Area of Review) and monitoring requirements to support cost and regulatory analysis, and for consideration of different post-injection site care and closure scenarios. This submission contains the current version of NRAP-Open-IAM available for evaluation and testing. To use the NRAP-Open-IAM, download the source code (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/4c24a3da-3b40-4ffe-9892-c807ae9f8760) then open the NRAP-Open-IAM user's guide (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669) to read more about the tool. Installation instructions for Windows, Mac, and Linux can be found in the "installers" folder of the extracted NRAP-Open-IAM folder and describe setup of environment (e.g., Python libraries) needed for proper work of the tool. Test of installation can be done by running "python openiam_setup_tests.py" in the "setup" folder. The installation test also runs a test suite to see if the NRAP-Open-IAM has been installed correctly. To run the test suite separately, run "python iam_test.py" in the "test" folder. User's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669 Developer's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/3bc6ee7d-609d-4eb6-80ba-fa6130ee0313 Reservoir simulation data used in some examples distributed with NRAP-Open-IAM: - Kimberlina: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/eb62cece-61b2-4037-9b6d-32407dde2ab8 - Kimberlina (compartmentalized): https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/366f9530-3b32-4b84-affe-ab2df1d9a8b5 - FutureGen 2.0: https://edx.netl.doe.gov/dataset/futuregen-2-0-1008-simulation-reservoir-lookup-table NRAP-Open-IAM GitLab repository: https://gitlab.com/NRAP/OpenIAM Related publications: - Bacon, D., Yonkofski, C., Brown, C., Demirkanli, D. and Whiting, J., 2019. Risk-based post injection site care and monitoring for commercial-scale carbon storage: Reevaluation of the FutureGen 2.0 site using NRAP-Open-IAM and DREAM. International Journal of Greenhouse Gas Control 90: 102784. - Bacon, D. Demirkanli, D., and White, S., 2020. Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site. International Journal of Greenhouse Gas Control 102: 103153. - Harp, D., Oldenburg, C., and Pawar, R., 2019. A metric for evaluating conformance robustness during geologic CO2 sequestration operations. International Journal of Greenhouse Gas Control 85: 100-108. - Lackey, G., Vasylkivska, V., Huerta, N., King, S., and Dilmore, R., 2019. Managing well leakage risks at a geologic carbon storage site with many wells, International Journal of Greenhouse Gas Control, 88 :182-194. - Vasylkivska, V., Dilmore, R., Lackey, G., Zhang, Y., King, S., Bacon, D., Chen, B., Mansoor, K., and Harp, D., 2021. NRAP-Open-IAM: A flexible open-source integrated assessment model for geologic carbon storage risk assessment and management, Environmental Modelling & Software, 143: 105114. Presentations: - Chen, B., Harp, D., and Pawar, R., A data assimilation approach (ES-MDA) coupling with NRAP-Open-IAM for quantifying uncertainty reduction in geological CO2 sequestration. AGUFM 2019: T44A-02. - Chen, B., and Harp, D., Improving risk analysis precision for geologic CO2 sequestration by quantifying the uncertainty reduction before and after acquiring monitoring data. 14th Greenhouse Gas Control Technologies Conference, Melbourne, Australia, 2018, pp. 21-26. - Harp, D., National Risk Assessment Partnership Task 2: Containment Assurance. No. LA-UR-19-28654, Los Alamos National Laboratory (LANL), Los Alamos, NM (United States), 2019. - Vasylkivska, V., King, S., Bacon, D., Harp, D., Chen, B., Mansoor, K., Onishi, T., Yang, Y., Zhang, Y., and Keating, E., NRAP-Open-IAM: An open-source integrated assessment model, poster, Mastering the Subsurface Through Technology Innovation, Partnerships and Collaboration: Carbon Storage and Oil and Natural Gas Technologies Review Meeting, Pittsburgh, PA, August 13-16, 2018. - Vasylkivska, V., Lackey, G., King, S., Wentworth, A., Huerta, N., Creason, C., DiGiulio, J., Yang, Y., and Dilmore, R., Long-term risk analysis of a geologic CO2 storage project during the post-injection period, SIAM Conference on Computational Science and Engineering, Spokane, WA, February 25-March 1, 2019. - Vasylkivska, V., Overview of the NRAP-Open-IAM tool for carbon storage (beta release), 2019 Annual NRAP Tool Users Meeting, Pittsburgh, PA, August 27, 2019. - Vasylkivska, V., Bacon, D., Chen, B., Dilmore, R., Harp, D., King, S., Lackey, G., Lindner, E., Liu, G., Mansoor, K. and Zhang, Y., NRAP-Open-IAM: A new, open-source code for integrated assessment of geologic carbon storage containment effectiveness and leakage risk, poster, American Geophysical Union Fall Meeting 2020 (virtual meeting), December 2020. - Vasylkivska, V., NRAP open-source integrated assessment model and relevant application, oral presentation, NRAP workshop "NRAP Tools for Geologic Carbon Storage Risk-Based Decision Making" held in conjunction with Groundwater Protection Council (GWPC) 2021 Annual Forum (virtual meeting), Salt Lake City, UT, September 2021. - Vasylkivska, V., NRAP-Open-IAM: open-source integrated assessment model, digital poster/demonstration, software demonstration session, 2022 Carbon Management Project Review Meeting, August 16, 2022

AoR↗

Argonne Bubble Experiment Thermal Model Development III

This report describes the continuation of the work reported in “Argonne Bubble Experiment Thermal Model Development” and “Argonne Bubble Experiment Thermal Model Development II”. The experiment was performed at Argonne National Laboratory (ANL) in 2014. A rastered 35 MeV electron beam deposited power in a solution of uranyl sulfate, generating heat and radiolytic gas bubbles. Irradiations were performed at beam power levels between 6 and 15 kW. Solution temperatures were measured by thermocouples, and gas bubble behavior was recorded. The previous report2 described the Monte-Carlo N-Particle (MCNP) calculations and Computational Fluid Dynamics (CFD) analysis performed on the as-built solution vessel geometry. The CFD simulations in the current analysis were performed using Ansys Fluent, Ver. 17.2. The same power profiles determined from MCNP calculations in earlier work were used for the 12 and 15 kW simulations. The primary goal of the current work is to calculate the temperature profiles for the 12 and 15 kW cases using reasonable estimates for the gas generation rate, based on images of the bubbles recorded during the irradiations. Temperature profiles resulting from the CFD calculations are compared to experimental measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tools Assessing Performance (FY2022 Quarter 4 Report)

We were finally successful in acquiring TTU’s 200-meter tower data for the duration of the Isolated Building Wake Experiment (IBWE). This data will provide an independent measurement of the undisturbed wind profile including heights far above those that are typical of meteorological sensor networks. These measurements will characterize the mean and turbulence throughout the atmospheric surface layer, which will be particularly useful for validation of highfidelity computational fluid dynamics models. This dataset can be shared freely among the TAP team but not outside of it as TTU sells this data to the wind power industry. The journal article overview has been submitted to LANL’s publication review and will be submitted to the Journal of Wind Engineering and Industrial Aerodynamics once the review is complete. We are nearing completion on another two articles. The first discusses analyses of the evolution of the measurements of mean wind speed along the centerline of the building as a function of atmospheric stability. Rather than using individual time steps, all of the wind speed measurements are normalized using the upwind measurement at 2 building heights (H) and filtered by wind direction and atmospheric stability. The second discusses the development of the diffusive wake model and uses the averaged profiles from the first article to validate the behavior of the model. The portion of the section below on the diffusive wake model discussing comparisons with the averaged profiles, will be included in the validation section of this paper.

17 WIND ENERGY↗

Dispersion mechanisms of a tidal river junction in the Sacramento–San Joaquin Delta, California

In branching channel networks, such as in the Sacramento–San Joaquin River Delta, junction flow dynamics contribute to dispersion of ecologically important entities such as fish, pollutants, nutrients, salt, sediment, and phytoplankton. Flow transport through a junction largely arises from velocity phasing in the form of divergent flow between junction channels for a portion of the tidal cycle. Field observations in the Georgiana Slough junction, which is composed of the North and South Mokelumne rivers, Georgiana Slough, and the Mokelumne River, show that flow phasing differences between these rivers arise from operational, riverine, and tidal forcing. A combination of Acoustic Doppler Current Profile (ADCP) boat transecting and moored ADCPs over a spring–neap tidal cycle (May to June 2012) monitored the variability of spatial and temporal velocity, respectively. Two complementary drifter studies enabled assessment of local transport through the junction to identify small-scale intrajunction dynamics. We supplemented field results with numerical simulations using the SUNTANS model to demonstrate the importance of phasing offsets for junction transport and dispersion. Different phasing of inflows to the junction resulted in scalar patchiness that is characteristic of MacVean and Stacey’s (2011) advective tidal trapping. Furthermore, we observed small-scale junction flow features including a recirculation zone and shear layer, which play an important role in intra-junction mixing over time scales shorter than the tidal cycle (i.e., super-tidal time scales). Thus, the study period spanned open- and closed-gate operations at the Delta Cross Channel. Synthesis of field observations and modeling efforts suggest that management operations related to the Delta Cross Channel can strongly affect transport in the Delta by modifying the relative contributions of tidal and riverine flows, thereby changing the junction flow phasing.

42 ENGINEERING↗

Uptake of N 2 O 5 by aqueous aerosol unveiled using chemically accurate many-body potentials

The reactive uptake of N 2 O 5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N 2 O 5 in water. The free energy profile highlights that N 2 O 5 is selectively adsorbed to the liquid–vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction–diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction–diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N 2 O 5 in more complex solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collaborative Research: The Role of Grain Boundary Structure and Chemistry in Materials Failure

Throughout the one year duration of this award significant progress was made in understanding the effect of grain boundaries in polycrystalline materials. The starting point of this work is that each grain boundary can be characterized by an interface energy and respective interface parameter, which depend on the boundary chemical composition and structure. As this is a new point of view for treating grain boundaries, an interdisciplinary approach was followed to prove the existence of such interface terms. Throughout the project, hence, theory, experiments and simulations were used. Molecular dynamics and dislocation density-function dynamics were performed to understand dislocation-grain boundary interactions in pure Al crystalline materials, as potentials for Al are well established. The formation of dislocation pile-ups was observed and the strain profiles were in agreement with interfacial gradient plasticity, which is the theoretical framework that allows the introduction of mechanically induced interface energies. After verifying numerically the existence of mechanically induced interface energies, nanoindentation experiments were performed near grain boundaries of Fe-Si crystals that contained three grain boundaries. By indenting in close proximity to the boundary these experiments allowed for measuring the critical stress at which deformation was transmitted across the boundary to the adjacent grains, either through dislocation transmission or absorption and subsequent emission. This critical stress, also termed grain boundary yield stress, was in accordance with the value predicted by interfacial gradient plasticity, further verifying that mechanically induced interface energies can help define the material response. In the Fe-Si crystals the Si content was kept constant at 3wt%, but three specimen were considered each of which was characterized by a different misorientation at the grain boundary. By relating the experiments with the analytical gradient plasticity expressions it was possible to obtain first estimates for the interface parameter as ~200N/m. This value was not proportional to the misorientation angle, but depended also on the grain boundary plane. The continuation of this project is performing more detailed experiments on additional Fe-Si and Fe-C systems (both nanoindentation and in-situ transmission electron microscopy) to further understand the physical origin of such interface parameters. Furthermore, molecular dynamics nanoindentation simulations are also performed on similar systems to better understand the dislocation grain boundary interactions.

36 MATERIALS SCIENCE↗