Search NASA⌕ Search

SEARCH · Search NASA

Results for “projected gradient method”

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 217 records · Page 12

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗

Final Report for Project 13-4791: New Mechanistic Models of Creep-Fatigue Crack Growth Interactions for Advanced High Temperature Reactor Components

This project developed and validated a novel, multi-scale, mechanism-based model to quantitatively predict creep-fatigue crack growth and failure for Ni-based Alloy 617 at 800°C. Alloy 617 is a target material for intermediate heat exchangers in Generation IV very high temperature reactor designs, and it is envisioned that this model will aid in the design of safe, long lasting nuclear power plants. The technical effectiveness of the model was shown by demonstrating that experimentally observed crack growth rates can be predicted under both steady state and overload crack growth conditions. Feasibility was considered by incorporating our model into a commercially available finite element method code, ABAQUS, that is commonly used by design engineers. While the focus of the project was specifically on an alloy targeted for Generation IV nuclear reactors, the benefits to the public are expected to be wide reaching. Indeed, creep-fatigue failure is a design consideration for a wide range of high temperature mechanical systems that rely on Ni-based alloys, including industrial gas power turbines, advanced ultra-super critical steam turbines, and aerospace turbine engines. It is envisioned that this new model can be adapted to a wide range of engineering applications.

36 MATERIALS SCIENCE↗

Integral Reactor Containment Condensation Model and Experimental Validation

This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). In the three years of investigation duration, following the original proposal, the planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). A film flow model was developed for the scaling analysis, and the results suggest that the 1/3 scaled test facility covers large portion of laminar film flow, leading to a lower average heat transfer coefficient comparing to the prototypic value. Although it is conservative in reactor safety analysis, the significant reduction of heat transfer coefficient (50%) could under estimate the prototypic condensation heat transfer rate, resulting in inaccurate prediction of the decay heat removal capability. Further investigation is thus needed to quantify the scaling distortion for safety analysis code validation. Experimental investigations were performed in the existing MASLWR test facility at OrST with minor modifications. A total of 13 containment condensation tests were conducted for pressure ranging from 4 to 21 bar with three different static inventories of non-condensable gas. Condensation and heat transfer rates were evaluated employing several methods, notably from measured temperature gradients in the HTP as well as measured condensate formation rates. A detailed mass and energy accounting was used to assess the various measurement methods and to support simplifying assumptions required for the analysis. Condensation heat fluxes and heat transfer coefficients are calculated and presented as a function of pressure to satisfy the objectives of this investigation. The major conclusions for those tests are summarized below: (1) In the steam blow-down tests, the initial condensation heat transfer process involves the heating-up of the containment heat transfer plate. An inverse heat conduction model was developed to capture the rapid transient transfer characteristics, and the analysis method is applicable to SMR safety analysis. (2) The average condensation heat transfer coefficients for different pressure conditions and non-condensable gas mass fractions were obtained from the integral test facility, through the measurements of the heat conduction rate across the containment heat transfer plate, and from the water condensation rates measurement based on the total energy balance equation. 15 (3) The test results using the measured HTP wall temperatures are considerably lower than popular condensation models would predict mainly due to the side wall conduction effects in the existing MASLWR integral test facility. The data revealed the detailed heat transfer characteristics of the model containment, important to the SMR safety analysis and the validation of associated evaluation model. However this approach, unlike separate effect tests, cannot isolate the condensation heat transfer coefficient over the containment wall, and therefore is not suitable for the assessment of the condensation heat transfer coefficient against system pressure and noncondensable gas mass fraction. (4) The average condensation heat transfer coefficients measured from the water condensation rates through energy balance analysis are appropriate, however, with considerable uncertainties due to the heat loss and temperature distribution on the containment wall. With the consideration of the side wall conduction effects, the results indicate that the measured heat transfer coefficients in the tests is about 20% lower than the prediction of Dehbi’s correlation, mainly due to the side wall conduction effects. The investigation also indicates an increase in the condensation heat transfer coefficient at high containment pressure conditions, but the uncertainties invoked with this method appear to be substantial. (5) Non-condensable gas in the tests has little effects on the condensation heat transfer at high elevation measurement ports. It does affect the bottom measurements near the water level position. The results suggest that the heavier non-condensable gas is accumulated in the lower portion of the containment due to stratification in the narrow containment space. The overall effects of the non-condensable gas on the heat transfer process should thus be negligible for tall containments of narrow condensation spaces in most SMR designs. Therefore, the previous correlations with noncondensable gas effects are not appropriate to those small SMR containments due to the very poor mixing of steam and non-condensable gas. The MELCOR simulation results agree with the experimental data reasonably well. However, it is observed that the MELCOR overpredicts the heat flux for all analyzed tests. The MELCOR predicts that the heat fluxes for CCT’s approximately range from 30 to 45 kW/m2 whereas the experimental data (averaged) ranges from about 25 to 40 kW/m2. This may be due to the limited availability of liquid film models included in MELCOR. Also, it is believed that due to complex test geometry, measured temperature gradients across the heat transfer plate may have been underestimated and thus the heat flux had been underestimated. The MELCOR model predicts a film thickness on the order of 100 microns, which agrees very well with film flow model developed in this study for scaling analysis. However, the expected differences in film thicknesses for near vacuum and near atmospheric test conditions are not significant. Further study on the behavior of condensate film is expected to refine the simulation results. Possible refinements include but are not limited to, the followings: CFD simulation focusing on the liquid film behavior and benchmarking with experimental analyses for simpler geometries. 16 1 INTRODUCTION This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). The experimental results are employed to validate the containment condensation model in reactor containment system safety analysis code for integral SMRs. Such a containment condensation model is important to demonstrate the adequate cooling. In the three years of investigation, following the original proposal, the following planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). The results are applicable to integral Small Modular Reactor (SMR) designs, including NuScale, mPower, Westinghouse SMR, Holtec-160 and other integral reactors with small containments of relatively high pressures under accidental conditions. Testing has been conducted at the OrSU laboratory in the existing MASLWR (Multi-Application Small Light Water Reactor) integral test facility sponsored by the US Department of Energy. Its highpressure stainless steel containment model (~2 MPa) is scaled to the NuScale SMR currently under development at NuScale Power, Inc.. Minor modifications to the model containment have been made to control the non-condensable gas fraction and to utilize the secondary loop stable steam flow for condensation testing. UW-Madison has developed a containment condensation model, which leveraged previous validated containment heat transfer work carried out at UW-Madison, and extended the range of applicability of the model to integral SMR designs that utilize containment vessels of high heat transfer efficiencies. In this final report, the research background and literature survey are presented in Chapter 2 and 3, respectively. The test facility description and modifications are summarized in Chapter 4, and the scaling analysis is introduced in Chapter 5. The tests description, procedures, and data analysis are presented in Chapter 6, while the numerical modeling is presented in Chapter 7, followed by a conclusion section in Chapter 8.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Accuracy of predictions made by machine learned models for biocrude yields obtained from hydrothermal liquefaction of organic wastes

Hydrothermal liquefaction (HTL) has potential for converting abundant wet organic wastes into renewable fuels. Because HTL consists of a complex reaction network, deterministic, physics-based prediction of its biocrude yield is prohibitively difficult. Data-driven methods provide an alternative to the physics-based approach; however, rigorous testing must be performed to ensure the accuracy of predictions made by data-driven methods. To this end, a data set was assembled consisting of 570 data points appearing in the open literature. The data set was divided into training, validation, and test sub-sets and used for evaluating different machine learning regression approaches to predict biocrude yield. Among the tested algorithms, Random Forest and eXtreme Gradient Boosting (XGBoost) predicted biocrude yields in a test set that had not been used for training with the greatest accuracy, with root mean square errors (RMSE) of 8.34 and 8.57, respectively. Further refinement of the Random Forest model reduced its RMSE to 8.07. In comparison, predictions of a series of literature models resulted in RMSE ranging from 9.16 in the most accurate case to 27.6 in the least accurate; most literature models yielded RMSE values > 10. Using biocrude yield predictions from the most accurate Random Forest model and a probabilistic economic analysis found that the model accuracy is sufficient to prioritize allocation of resources based on projected minimum fuel selling price. In our report the models and analysis represent a major advance in the ability to use readily available data to predict biocrude yields on new feedstocks that have not previously been studied.

42 ENGINEERING↗

Pathways to Carbon Liberation: A Systems Approach to Understanding Carbon Transformations and Losses from Thawing Permafrost

Our objective in this project was to discover how microbial communities mediate the fate of carbon in thawing permafrost under climate change. We proposed a systems approach integrating (a) molecular microbial and viral ecology, (b) molecular organic chemistry and stable and radiocarbon isotopes, and (c) state-of-the-art modeling, along a chronosequence of permafrost thaw in subarctic Sweden. The fate of carbon (C) in thawing permafrost is an outstanding challenge of modern biogeochemistry and climate change. Permafrost C pools are large (~1700 PgC), and C dynamics of thawing permafrost complex: old C decomposes as it is liberated from thawing permafrost as CO 2 or CH 4 , with a significant fraction cycling through lake sediments, even as new C accumulates due to thaw-initiated ecological succession. Our work is allowing better prediction of the net effect of these processes. Microbes mediate C dynamics in thawing permafrost, but a mechanistic understanding of how to scale microbial population dynamics, genomic potential, and expression to ecosystem-scale processes has been missing. A key question was: What is the interplay of microbial communities and organic matter chemical structure in the decomposition/preservation of organic C across a thaw gradient? And intriguingly, what (if any) is the role of phage (viruses that infect prokaryotic cells) in mediating these processes? Viruses appear to play a large role in driving oceanic function, but these phenomena are virtually unstudied in terrestrial systems. This endeavor linking microbial and viral dynamics, organic geochemistry and trace gas production will improve models of C cycling in permafrost systems, and clarify the fate of C under future climates. Technical Approach – We conducted a study at Stordalen Mire, Sweden, along a permafrost thaw chronosequence encompassing permafrost palsas and their initial collapse stage, thawing bog sites (dominated by Sphagnum moss), fully-thawed and inundated fen sites (dominated by Eriophorum spp or Carex spp.), and lakes. Our project used cutting edge technologies in both biogeochemistry and molecular microbial ecology to advance systems biology research on microbial carbon cycling through: (a) systems-level mapping of chemical states and ages of organic matter (via FT-ICR-MS and 14 C analysis) along thaw gradients to associated microbial communities, biochemical potentials, and activities (via meta-genomics, -transcriptomics, –proteomics, and viral genomics), and to CO 2 and CH 4 fluxes; (b) experimental incubations using Quantum Dot Probing, to test key hypotheses, arising from (a), about particular microbially-driven biochemical degradation pathways; (c) bioinformatics designed to simultaneously enhance the DOE Knowledgebase, and (d) integrated ecosystem C-cycle modeling testable by soil organic chemical and microbial data. Results from this project spanned 4 categories: (i) integrative across disciplines, and disciplinarily-focused in (ii) microbiology, (iii) biogeochemistry, and (iv) modeling. In integrated understanding of the system, we answered the questions posed above, mapping resident microbes to carbon transformations to reveal a carbon-processing shift with permafrost thaw, and discovering abundant, diverse and habitat-specific viruses that appear to be impacting carbon cycling, while also relating microbial communities to modeled processes. In microbiology-specific papers, we dug into viromes, microbial community networks, and informatic methods, while also introducing a new permafrost-associated phylum and its ecogenomic context. Relevance to DOE FOA - With a focus on discovering how microbial communities mediate the fate of old and new carbon in permafrost systems with climate change, our work directly responds to the call for “–omics driven basic research on the contribution of... microbial communities to C cycling processes in terrestrial ecosystems.” We bridge all three of the focal areas of the call, including (i) “Systems biology studies”, and development of –omics approaches (ii) “to investigate microbial community functional processes”, and (iii) “for imaging and analysis of microbially-mediated carbon cycling.”

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermally Insulating Transparent Barrier (THINNER) coatings on single pane windows

Conventional silica aerogel monoliths can provide remarkable thermal insulation but the presence of large pores (> 30 nm) tends to scatter visible light and render the material opaque or translucent instead of transparent. In addition, they are prone to cracking during synthesis and handling which makes them difficult to integrate in products and in particular in window solutions. This project developed two new solgel synthesis methods using silica precursors or preformed silica nanoparticles and ambient drying to produce mesoporous organo-silica monoliths. The monoliths were (i) thermally insulating, (ii) optically transparent, (iii) flexible, and (iv) hydrophobic. They feature porosity ranging from 50% to 90% with narrow pore size distribution with pore less than 20 nm resulting in excellent optical clarity (haze < 2%) and very low thermal conductivity (< 30 mW/mK). Interestingly, not only porosity but also pore size and mass fractal dimension were found to affect the thermal conductivity of the mesoporous silica. The superior transparency of the monoliths was shown to be attributed to dependent scattering among silica nanoparticles. Flexibility was achieved through trimethylchlorosilane surface modification. Furthermore, process scale-up and integration of the ambigel monoliths into window solutions using optically clear adhesives were also demonstrated in 6”x6” double-pane windows. The aerogel and the window solution were shown to be durable to accelerated aging under UV, moisture, and/or temperature gradient and thermal cycling. Overall, the window solution achieved the technical performance and the cost target of $10/sqft set for the SHIELD program. However, scaling to industry relevant scale (> 10’x10’) remains a challenge due to requirements on the drying process and fume hood size and to the propensity of large aerogel slab to crack during drying.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Linking international technical specifications for acoustic characterization of marine energy converter sounds with environmental compliance criteria

As new ocean energy technologies emerge and are deployed for testing and operations, sound emissions are a potential concern for environmental effects to marine life. Consistent acoustic measurement and data analysis methods can help promote comparisons of technologies and transferability between project sites. In 2022, acoustic emissions from a prototype scale wave energy converter (WEC) were characterized for a range of environmental conditions and power generation states in the coastal waters off southern California using a set of international technical specifications. Results from the international technical specification analyses were applied to United States regulatory threshold criteria for acoustic impacts to marine mammals and examined in the context of European underwater noise monitoring guidelines. Weighted 24 hour cumulative sound exposure levels SEL24h calculated from the highest power generation state WEC sound pressure levels were often more than 20 dB below threshold criteria for temporary threshold shifts in five relevant marine mammal hearing groups. Following European Union recommendations for analyses and reporting, WEC sound characterization in third octave bands centered at 63 Hz and 125 Hz show clear spatial decay of WEC-generated noise, with more pronounced attenuation at 63 Hz, and a less marked but still detectable gradient at 125 Hz, collectively suggesting a relatively confined acoustic footprint under the observed conditions. The value of the international technical specification approach is highlighted by the isolation of WEC sounds from the surrounding soundscape. This allows for a robust characterization of acoustic emissions through a range of device power generation and sea states. Furthermore, in threshold-based regulatory contexts like the U.S., this facilitates direct evaluation of source contributions, while in broader monitoring frameworks used in the E.U. it provides a reproducible foundation for assessing the contribution of emerging ocean energy technologies to the underwater acoustic environment.

Haxel, Joseph H. (ORCID:0000000273864761)↗

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗

Thermal / structural analysis of the HB 650 thermal shield

Fermilab’s PIP-II project’s superconducting linear accelerator will drive the next generation of particle accelerators through a revolution in beam intensity. Key to beam intensification are the high efficiencies of niobium-tin superconducting radiofrequency (SRF) cavities operating at cryogenic temperatures near 5 Kelvin. A multifaceted approach is employed to achieve and maintain the extreme temperature. Vacuum provides the first barrier to thermal convection. Physical thermal intercept zones further isolate the exterior shell vacuum vessel at 300 Kelvin from the supercooled beamline. The first thermal intercept lies just inside the exterior vacuum vessel forming the 40 Kelvin zone. This 40K thermal shield is the focus of this investigation. In preparation for operation, the thermal shield is cooled from 300K to 40K in a slow process over 2-1/2 days. The factors limiting the cooling rate are a high thermal gradient which produces thermal strain and high mechanical stress. D ecreasing the time required to reach operational temperatures is the desired research goal. Physical testing data and finite element method (FEM) computer models from the Single Spoke Resonator – 1 (SSR1) and Linac Coherent Light Source (LCLS-II) cryomodules will be evaluated to determine needed design changes. The thermal shield design for the High Beta 650 MHz (HB650) will be modified and analyzed using FEM computer models to determine maximum cooling rate with a target of 50 thermal cycles of life. Maximum allowable stress will be determined through low cycle fatigue calculations.

43 PARTICLE ACCELERATORS↗

Numerical Analysis of High-Order Modes in SRF Resonators for Particle Accelerators

Over the past decades, superconducting technology has rapidly evolved towards high accelerating gradients and low surface resistance, making it possible to operate particle accelerators with high average beam currents and large duty factors. However, RF losses due to coherent excitation of the HOM become the limiting factor for these regimes. Unlike the cavity operating mode, which is tuned separately, the HOM parameters can significantly vary from one cavity to another due to finite mechanical tolerances during the manufacturing process. Thus, it is of utmost importance to know the HOM parameter spread in advance in order to predict unexpected cryogenic losses, overheating of beam line components and maintain stable beam dynamics. In this paper, we present a method for generating cavity geometry with an arbitrary spread of mechanical imperfections and numerically evaluating HOM statistics. Knowing the spread of HOM parameters, we calculated the probability of resonant HOM losses in SRF accelerating cavities used in CW beam current machines such as the PIP-II and LCLS-II linacs, as well as for the SRF crab-cavity for the ILC project. Finally, we present experimental results of HOM spectra measurements in hundreds of 1.3 GHz cavities installed in LCLS-II cryomodules. Studying the effects of HOM excitation results in specifications of the SRF cavity and cryomodule and can significantly impact the efficiency and reliability of the machine operation.

Lunin, Andrei [Fermilab] (ORCID:0000000290096792)↗

Recent Advances in Assessing Environmental Effects of Marine Renewable Energy Around the World

Marine renewable energy (MRE) is increasingly of interest to coastal nations as a source of renewable energy that can support climate change mitigation goals as well as provide secure locally-produced energy for coastal and island communities. MRE extracts power from tidal streams, waves, ocean currents, run of rivers, and gradients in the ocean, with specialized devices developed and tested for each energy resource. Alongside development of MRE technologies and systems, first in Europe and then in North America, Australia, Asia, and other regions, it has been universally recognized that there is also a need to examine potential effects on marine animals, habitats, ecosystem processes, local communities and other sea users, to ensure that the MRE industry can be developed in a responsible and sustainable manner. This paper looks at the status of assessment and monitoring for potential environmental effects associated with MRE projects around the world. Over 80 projects were identified worldwide as having been tested, demonstrated, or commercially deployed with associated environmental monitoring. Five of the projects that represent tidal stream, wave, and run of river projects are examined in more detail to determine the types of data and information collected for those projects, the outputs of the monitoring campaigns, and the actions taken as a result of the data collection and analysis. Recommendations are provided for standardization of the monitoring approaches, instruments, and analysis methods at MRE project sites worldwide.

marine renewable energy, environmental effects, ri↗

TU-H-206-07: Assessment of Geometric Distortion in EPI with a SPAMM Tagged Acquisition

Purpose: Echo planar imaging (EPI) is known to exhibit gross geometric distortion caused by multiple factors, including B0 inhomgeneity and transient eddy currents. However, diffusion weighted (DW) EPI has become indispensable for diagnosis and therapy assessment. We propose a methodology for quantifying distortion in EPI sequences that does not require the use of dedicated spatial accuracy phantoms, enabling flexibility in phantom design for QA of distortion effects in EPI protocols. Methods: The proposed methodology utilizes a saturation technique known as Spatial Modulation of Magnetization (SPAMM) that tags the imaging subject with saturated grid lines. Originally intended for tracking cardiac motion, these grids are applied to assess differences between diffusion weighting directions and b-values, or against a more geometrically robust sequence such as fast spin echo (FSE). The saturation preparation sequence consists of binomially weighted (e.g. 1-3-3-1) pulses interleaved with gradient blips along the frequency encode direction, followed by the same sequence with gradient blips in the phase encode direction. Three phantoms were assessed with these sequences: a spherical head-sized phantom, a large shimming phantom, and a modified PET ACR phantom that included compartments of water, air, oil, and Teflon. Each phantom was acquired with three sequences using parameters from a clinically appropriate protocol (22 cm head or 46 cm abdomen): a conventional DW-EPI sequence (3 DW directions), and both the DW-EPI and FSE sequences with tagging. Differences in grid locations were visualized with minimum intensity projection between images, and measured using intersecting locations on the grids. Results: Grid lines were clearly visualized on tagged images and enabled quantification of distortions. Maximum eddy current induced errors of 10.8 to 14.8 mm were observed in areas away from isocenter with DW gradients applied in various directions. Conclusion: SPAMM tagging provides a promising mechanism for assessing spatial distribution of distortion in EPI sequences.

60 APPLIED LIFE SCIENCES↗

Radiopure ZnMoO{sub 4} scintillating bolometers for the LUMINEU double-beta experiment

The results of R&D of radiopure zinc molybdate (ZnMoO{sub 4}) based scintillating bolometers for the LUMINEU (Luminescent Underground Molybdenum Investigation for NEUtrino mass and nature) double-beta decay experiment are presented. A dedicated two-stage molybdenum purification technique (sublimation in vacuum and recrystallization from aqueous solutions) and an advanced directional solidification method (the low-thermal-gradient Czochralski technique) were utilized to produce high optical quality large mass (∼1 kg) ZnMoO{sub 4} crystal boules and first {sup 100}Mo (99.5%) enriched Zn{sup 100}MoO{sub 4} crystal scintillator (mass of ∼0.2 kg). Scintillating bolometers based on ZnMoO{sub 4} (≈ 0.33 kg) and Zn{sup 100}MoO{sub 4} (≈ 0.06 kg) scintillation elements and high purity Ge wafers were tested in the EDELWEISS set-up at the Modane Underground Laboratory (France). Long term low temperature tests demonstrate excellent detectors’ performance and effectiveness of the purification and solidification procedures for the achievement of high radiopurity of the material, in particular with a bulk activity of {sup 228}Th and {sup 226}Ra below 4 µBq/kg. The adopted protocol was used to produce for the first time a large volume Zn{sup 100}MoO{sub 4} crystal scintillator (mass of ∼1.4 kg, {sup 100}Mo enrichment is 99.5%) to search for neutrinoless double-beta decay of {sup 100}Mo in the framework of the LUMINEU project.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗

Quantum Mechanical Simulations of Complex Nanostructures for Photovoltaic Applications

A quantitative understanding of the electronic excitations in nanostructures, especially complex nanostructures, is crucial for making new-generation photovoltaic (PV) cells based on nanotechnology, which have high efficiency and low cost. Yet current quantum mechanical simulation methods are either computationally too expensive or not accurate and reliable enough, hindering the rational design of the nanoscale PV cells. The PI seeks to develop new methodologies to overcome the challenges in this very difficult and long-lasting problem, pushing the field forward so that electronic excitations can be accurately predicted for systems involving thousands of atoms. The primary objective of this project is to develop new approaches for electronic excitation calculations that are more accurate than traditional density functional theory (DFT) and are applicable to systems larger than what current beyond-DFT methods can treat. In this proposal, the PI will first address the excited-state problem within the DFT framework to obtain quasiparticle energies from both Kohn-Sham (KS) eigenvalues and orbitals; and the electron-hole binding energy will be computed based on screened Coulomb interaction of corresponding DFT orbitals. The accuracy of these approaches will be examined against many-body methods of GW/BSE and quantum Monte Carlo (QMC). The PI will also work on improving the accuracy and efficiency of the GW/BSE and QMC methods in electronic excitation computations by using better KS orbitals obtained from orbital-dependent DFT as inputs. Then an extended QMC database of ground- and excited-state properties will be developed, and this will be spot checked and supplemented with data from GW/BSE calculations. The investigation will subsequently focus on the development of an improved exchange-correlation (XC) density functional beyond the current generalized gradient approximation (GGA) level of parameterization, with parameters fitted to the QMC database. This will allow the ground-state properties of focus systems to be more precisely predicted using DFT. These new developments will then be applied to investigate a chosen set of complex nanostructures that have great potential for opening new routes in designing materials with improved transport, electronic, and optical properties for PV and other optoelectronic usages: (1) Hybrid interfaces between materials with distinct electronic and optical properties, such as organic molecules (conjugated polymers, e.g. P3HT) and inorganic semiconducting materials (Si and ZnO). Complicated interface structures, including interface bonding configurations, compositional and geometrical blending patterns, interfacial defects, and various sizes and shapes of inorganic nanomaterials, will be considered for the purpose of understanding the working mechanisms of present organic/nano PV systems and designing optimum interface structures for fast charge separation and injection. (2) Complex-structured semiconducting nanomaterials that could induce charge separation without pn- or hetero-junctions. The new methodology will allow the PI to investigate the performance of realistic semiconducting nanomaterials of internal (impurities, defects, etc.) and external (uneven surface, mechanical twisting and bending, surface chemistry, etc.) complexities on optical absorption and charge transport against charge trapping and recombination. Of particular interest is whether such structural complexity in a single material could even be beneficial for PV usage, for example, charge separation through morphology control. Successful completion of the proposed DFT methodology would have a far-reaching impact on our ability to study and exploit the nature of electronic excitations in complex materials, advancing the design of next-generation electronic and optoelectronic devices in all facets of renewable energy conversion and storage, including photovoltaics, thermoelectricity, photochemistry, etc.

14 SOLAR ENERGY↗

IRNet: A General Purpose Deep Risidual Regression Framework for Materials Discovery

Materials discovery is crucial for making scientific advances in many domains. Collections of data from experiments and first-principle computations have spurred interest in applying machine learning methods to create predictive models capable of mapping from composition and crystal structures to materials properties. Generally, these are regression problems with the input being a 1D vector composed of numerical attributes representing the material composition and/or crystal structure. While neural networks consisting of fully connected layers have been applied to such problems, their performance often suffers from the vanishing gradient problem when network depth is increased. Hence, predictive modeling for such tasks has been mainly limited to traditional machine learning techniques such as Random Forest. In this paper, we study and propose design principles for building deep regression networks composed of fully connected layers with numerical vectors as input. We introduce a novel deep regression network with individual residual learning, IRNet, that places shortcut connections after each layer so that each layer learns the residual mapping between its output and input. We use the problem of learning properties of inorganic materials from numerical attributes derived from material composition and/or crystal structure to compare IRNet’s performance against that of other machine learning techniques. Using multiple datasets from the Open Quantum Materials Database (OQMD) and Materials Project for training and evaluation, we show that IRNet provides significantly better prediction performance than the state-of-the-art machine learning approaches currently used by domain scientists. We also show that IRNet’s use of individual residual learning leads to better convergence during the training phase than when shortcut connections are between multi-layer stacks while maintaining the same number of parameters.

deep learning↗

Comparisons Between Array Derived Dynamic Strain Rate (ADDS) and Fiber‐Optic Distributed Acoustic Sensing (DAS) Strain Rate

Abstract Distributed acoustic sensing (DAS) strain rate and particle velocity can be compared through approximate scaling with medium velocity. We instead performed a direct comparison between array derived dynamic strain (ADDS) rate and DAS strain rate for six frequency bands. The PoroTomo project at Brady's Hot Springs, Nevada, deployed a 240‐geophone 3C array co‐located with fiber‐optic DAS system and 8.7 km of buried cable. We selected subsets of the geophone array to create four smaller arrays and computed ADDS. The horizontal components of the ADDS were rotated into the direction of the fiber‐optic cable and then compared with the observed DAS strain rates. From three example regional earthquakes of local magnitudes 2.9, 4.1, and 4.3, the ADDS are found to be coherent with DAS for frequencies ≤1 Hz. For frequencies >1‐Hz, this correlation decays quickly. Small differences between linear and areal dynamic strains at 1‐Hz suggest poor signal‐to‐noise or localized strain that is perturbed by shallow heterogeneities compare to the average strain propagating across the geophone array. The implication is that around 1‐Hz, straight fiber DAS is measuring axial strain along the fiber and can provide good approximations to translational particle motions. However, above 1‐Hz, DAS becomes more sensitive to shallow velocity gradients that can be beneficial for geophysical imaging yet becomes a limitation for traditional seismic analysis methods depending on absolute amplitude and phase from translational particle motions.

58 GEOSCIENCES↗