SEARCH · Engineering Papers

Results for “COUPLINGS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

DOE OSTI2018

GMD Coupling to Power Systems and Disturbance Mitigation

Presentation includes slides on Geomagnetic Disturbance: Ground Fields; Geomagnetic Disturbance: Coupling to Bulk Electric System; Geomagnetic Disturbance: Transformers; GMD Assessment Workflow (TPL-007-1); FERC order 830; Goals; SuperMag (1 min data) Nov. 20-21, 2003 Storm (DST = -422); Spherical Harmonics; Spherical Harmonics Nov. 20-21, 2003 Storm (DST = -422); DST vs HN0,0; Fluctuations vs. DST; Fluctuations; Conclusions and Next Steps; GMD Assessment Workflow (TPL-007-1); EMP E3 Coupling to Texas 2000 Bus Model; E3 Coupling Comparison (total GIC) Varying Ground Zero; E3 Coupling Comparison (total MVAR) Varying Ground Zero; E3 Coupling Comparison (GIC) at Peak Ground Zero; E3 Coupling Comparison (GIC) at Peak Ground Zero; and Conclusion.

58 GEOSCIENCES↗

DOE OSTI2019

Selection of glasses in support of the Sludge Batch 9 variability study for coupled operation with the Salt Waste Processing Facility

Prior to the processing of Sludge Batch 9 (SB9) with Frit 803 at the Defense Waste Processing Facility (DWPF), assessments and experimental glass work were completed to support the qualification of SB9 coupled processing with streams from the Actinide Removal Process (ARP) and Modular Caustic Side Solvent Extraction Unit (MCU). This previous work did not address coupled processing with streams from the Salt Waste Processing Facility (SWPF) because the Product Composition Control System (PCCS) models limited the TiO 2 concentration in glass to 2 weight percent (wt.%). These high activity streams from SWPF include monosodium titanate (MST) and sludge solids from the Sludge Solids Receipt Tank as well as Cs-containing strip effluent. The incorporation of these SWPF streams is expected to reach TiO2 concentrations in glass greater than 2 wt.% based on an MST addition of 0.4 g MST/L of salt solution. If double strike operation is required, then the MST addition is increased to 0.8 g MST/L of salt solution. Therefore, revisions to the models contained in PCCS were completed to allow the evaluation of glasses containing greater than 2 wt.% TiO 2 , but less than 6 wt.%. To support the implementation of the revised PCCS for SB9 – Frit 803 processing at DWPF prior to SWPF startup, the applicability of the new durability models for the SB9 variability study glasses was verified. Since the initial SB9 qualification for ARP-MCU coupled operation, additional Tank 40 analyses have been completed, which served as the basis for an updated projection of the SB9 blend composition in Tank 40. To identify a frit optimized for SB9 coupled operation with SWPF, the Savannah River National Laboratory (SRNL) performed a paper-study by utilizing Measurement Acceptance Region (MAR) assessments with DWPF PCCS. Based on these MAR assessment results, Frit 625 was recommended for processing both single and double MST strike scenarios. Savannah River Remediation (SRR) requested that SRNL perform a variability study with both Frit 803 and Frit 625 for SB9 coupled operation with SWPF. The objective of the variability study is to demonstrate that the Product Consistency Test (PCT) responses of the glasses are acceptable relative to the chemical durability of the Environmental Assessment (EA) reference glass and predictable by the current DWPF PCCS models for durability. This report documents the development of a test matrix that will support the variability study. To provide operational flexibility prior to SWPF startup, the predictability of glasses within the glass region for Frit 625-SB9 coupled operation with ARP-MCU was also confirmed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Coupling CTH to Linear Acoustic Propagation across an Air-Earth Interface

The interface between the Earth and the atmosphere forms a strong contrast in material properties. As such, numerical issues can arise when simulating an elastic wavefield across such a boundary when using a numerical simulation scheme. This is exacerbated when two different simulation codes are coupled straddling that interface. In this report we document how we implement the coupling of CTH, a nonlinear shock physics code, to a linearized elastic/acoustic wave propagation algorithm, axiElasti, across the air-earth interface. We first qualitatively verify that this stable coupling between the two algorithms produces expected results with no visible effects of the coupling interface. We then verify the coupling interface quantitatively by checking consistency with results from previous work and with coupled acoustic-elastic seismo-acoustic source inversions in three earth materials.

58 GEOSCIENCES↗

DOE OSTI2020

Coupling of CTF and TRANSFORM using the Functional Mockup Interface

An in-memory coupling between the sub-channel thermal hydraulics code COBRA-TF (CTF), which is included in the Virtual Environment for Reactor Applications (VERA), and the systems code Transient Simulation Framework of Reconfigurable Models (TRANSFORM) was developed in this work. Data exchange is accomplished by using the Functional Mock-Up Interface (FMI), an open standard for coupling dynamic systems models together. The FMI-based coupling necessitated the development of a novel FORTRAN wrapper for communicating with Functional Mock-Up Units. This wrapper will facilitate future FMI-based code couplings with other FORTRAN-based programs. CTF-TRANSFORM is exercised on a simplified Molten Salt Reactor Experiment (MSRE) model in steady-state and transient configurations. The coupled CTF-TRANSFORM model is shown to predict core temperature deltas similar to those available in historical MSRE operational data and is shown to be robust to fast power transients.The coupling paves the way to performing sensitivity and design studies impossible with VERA/CTF alone, in which secondary and tertiary loop operating conditions and design parameters may be perturbed and their impact on the core operating conditions studied.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Coupled Ocean Wave-Atmosphere Models for Offshore Wind Energy Applications

Coupling ocean wave models to mesoscale atmospheric models is necessary to represent the effect of waves on wind turbine hub-height winds. In this report we provide a review of the most widely used ocean waves models and the phased-averaged spectral wave modeling paradigm that they are based on. Methodologies used to couple these wave models to mesoscale atmospheric models are described along with details of existing coupled modeling systems. We summarize impacts on offshore wind resource assessments that have been investigated with such coupled modeling systems to-date. Specifically, coupling in the North and Baltic Seas was shown to have a small negative effect on the offshore wind energy resource at weak-to-moderate wind speeds, with little impact at higher wind speeds. Finally, limitations of these existing coupled modeling systems and impact assessments are discussed, including examples of and the potential for improved parameterizations of the air-sea fluxes, and the use of fine-scaled simulations using phased-resolved wave models.

17 WIND ENERGY↗

DOE OSTI2022

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

An Immersed Finite Element Lagrangian-Eulerian Code-Coupling Framework

This report presents an assessment of immersed Eulerian-Lagrangian code-coupling techniques suitable for use in a broad range of mechanics applications. The coupling algorithm is based on an immersed finite element method that considers the Lagrangian and Eulerian overlap regions in the overall variational formulation. In this report the basic formulation details are presented followed by various aspects of the code-coupling algorithm using OpenIFEM as the Lagrangian/coupling framework. A series of representative test cases that illustrate the code-coupling algorithm are discussed. The current work provides an in-depth investigation into the immersed finite element method for the purposes of providing a rigorous coupling technique that is minimally invasive in the respective Eulerian and Lagrangian codes. A number of extensions to the base immersed finite element method have been examined. These extension include nodal and quadrature-based indicator functions, a Lagrangian volume-fraction calculation in regions of overlap, and the use of penalty constraints between the Lagrangian and Eulerian domains. A unique MPI-based coupling strategy that retains the independent MPI structure of each code has been demonstrated.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Use Cases and Model Development of Thermal Storage Coupling for Advanced Nuclear Reactors

This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (ALWR) and a high-temperature gas-cooled reactor (HTGR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. Cost functions for the A-LWR were derived from the fully balanced models that were developed based on three different coupling options with three different thermal energy bypass ratios. For the next steps, cost functions for the HTGR will also be derived, and additional nuclear reactors (e.g., a liquid-cooled fast reactor [LFR] or molten-salt reactor [MSR]) will be evaluated for coupling with TES in similar fashion, including the evaluation of their steady-state condition models and cost functions. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state operation with minimum disruption to the nuclear power generation cycle. This report presents the results of steady state models developed using Aspen HYSYS®, wherein the thermal energy bypass for an NPP-TES coupling was varied up to 50%. The various components were sized using Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The current steady-state models and cost functions provide a baseline for additional work focusing on dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools to evaluate the technoeconomic viability and transient operations of TES-coupled A-NPPs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

The Strange Case of Ground-Coupled Airwaves on Seismoacoustic Stations at Local to Near-Regional Scales

Here we investigate the application of ground-coupled airwaves observed by seismoacoustic stations at local to near-regional scales to detect signals of interest and determine back-azimuth information. Ground-coupled airwaves are created from incident pressure waves traveling through the atmosphere that couple to the earth and transmit as a seismic wave with retrograde elliptical motion. Previous studies at sub-local scales (<10 km from a source of interest) found the back-azimuth to the source could be accurately determined from seismoacoustic signals recorded by acoustic and 3-component seismic sensors spatially separated on the order of 10 to 150 m. The potential back-azimuth directions are estimated from the coherent signals between the acoustic and vertical seismic data, via a propagation-induced phase shift of the seismoacoustic signal. A unique solution is then informed by the particle motion of the 3-component seismic station, which was previously found to be less accurate than the seismoacoustic-sensor method. We investigate the applicability of this technique to greater source-receiver distances, from 50-100 km and up to 400 km, which contains pressure waves with tropospheric and stratospheric ray paths, respectively. Specifically, we analyze seismoacoustic sources with ground truth from rocket motor fuel elimination events at the Utah Test and Training Range (UTTR) as well as a 2020 rocket launch in Southern California. From these sources we observe evidence that while coherent signals can be seen from both sources on multiple seismoacoustic station pairs, the determined ground-coupled airwave back-azimuths are more complicated than results at more local scales. Our findings suggest more complex factors including incidence angle, coupling location, subsurface material, and atmospheric propagation effects need to be fully investigated before the ground-coupled airwave back-azimuth determination method can be applied or assessed at these further distances.

58 GEOSCIENCES↗

DOE OSTI2024

Assessing Impacts of Waves on Hub-Height Winds off the U.S. West Coast Using Lidar Buoys and Coupled Modeling Approaches

Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.

17 WIND ENERGY↗

DOE OSTI2025

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB+QRPA calculations with the Gogny interaction

The Zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98Zr. This region includes a benchmarking case and offers insights into nuclear deformation phenomena. To investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially-symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the Quasiparticle Random Phase Approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the K π = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the K π = 0 − and K π = 1 − components. Intrinsic transition densities and associ ated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e. we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Ideas and perspectives: how coupled is the vegetation to the boundary layer?

Understanding the sensitivity of transpiration to stomatal conductance is critical to simulating the water cycle. This sensitivity is a function of the degree of coupling between the vegetation and the atmosphere and is commonly expressed by the decoupling factor. The degree of coupling assumed by models varies considerably and has previously been shown to be a major cause of model disagreement when simulating changes in transpiration in response to elevated CO 2 . The degree of coupling also offers us insight into how different vegetation types control transpiration fluxes, which is fundamental to our understanding of land–atmosphere interactions. To explore this issue, we combined an extensive literature summary from 41 studies with estimates of the decoupling coefficient estimated from FLUXNET data. We found some notable departures from the values previously reported in single-site studies. There was large variability in estimated decoupling coefficients (range 0.05–0.51) for evergreen needleleaf forests. This is a result that was broadly supported by our literature review but contrasts with the early literature which suggests that evergreen needleleaf forests are generally well coupled. Estimates from FLUXNET indicated that evergreen broadleaved forests were the most tightly coupled, differing from our literature review and instead suggesting that it was evergreen needleleaf forests. We also found that the assumption that grasses would be strongly decoupled (due to vegetation stature) was only true for high precipitation sites. These results were robust to assumptions about aerodynamic conductance and, to a lesser extent, energy balance closure. Hence, these data form a benchmarking metric against which to test model assumptions about coupling. Our results identify a clear need to improve the quantification of the processes involved in scaling from the leaf to the whole ecosystem. Progress could be made with targeted measurement campaigns at flux sites and greater site characteristic information across the FLUXNET network.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

Abstract. A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.

58 GEOSCIENCES↗

DOE OSTI2018

The path to CAM6: coupled simulations with CAM5.4 and CAM5.5

This paper documents coupled simulations of two developmental versions of the Community Atmosphere Model (CAM) towards CAM6. The configuration called CAM5.4 introduces new microphysics, aerosol, and ice nucleation changes, among others to CAM. The CAM5.5 configuration represents a more radical departure, as it uses an assumed probability density function (PDF)-based unified cloud parameterization to replace the turbulence, shallow convection, and warm cloud macrophysics in CAM. This assumed PDF method has been widely used in the last decade in atmosphere-only climate simulations but has never been documented in coupled mode. Here, we compare the simulated coupled climates of CAM5.4 and CAM5.5 and compare them to the control coupled simulation produced by CAM5.3. We find that CAM5.5 has lower cloud forcing biases when compared to the control simulations. Improvements are also seen in the simulated amplitude of the Niño-3.4 index, an improved representation of the diurnal cycle of precipitation, subtropical surface wind stresses, and double Intertropical Convergence Zone biases. Degradations are seen in Amazon precipitation as well as slightly colder sea surface temperatures and thinner Arctic sea ice. Simulation of the 20th century results in a credible simulation that ends slightly colder than the control coupled simulation. The authors find this is due to aerosol indirect effects that are slightly stronger in the new version of the model and propose a solution to ameliorate this. Altogether, in these early coupled simulations, CAM5.5 produces a credible climate that is appropriate for science applications and is ready for integration into the National Center for Atmospheric Research's (NCAR's) next-generation climate model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A comparative study of two-way and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US: performance evaluation and impacts of chemistry–meteorology feedbacks on air quality

The two-way coupled Weather Research and Forecasting and Community Multiscale Air Quality (WRF-CMAQ) model has been developed to more realistically represent the atmosphere by accounting for complex chemistry–meteorology feedbacks. In this study, we present a comparative analysis of two-way (with consideration of both aerosol direct and indirect effects) and offline coupled WRF v3.4 and CMAQ v5.0.2 over the contiguous US. Long-term (5 years from 2008 to 2012) simulations using WRF-CMAQ with both offline and two-way coupling modes are carried out with anthropogenic emissions based on multiple years of the U.S. National Emission Inventory and chemical initial and boundary conditions derived from an advanced Earth system model (i.e., a modified version of the Community Earth System Model/Community Atmospheric Model). The comprehensive model evaluations show that both two-way WRF-CMAQ and WRF-only simulations perform well for major meteorological variables such as temperature at 2 m, relative humidity at 2 m, wind speed at 10 m, precipitation (except for against the National Climatic Data Center data), and shortwave and longwave radiation. Both two-way and offline CMAQ also show good performance for ozone (O3) and fine particulate matter (PM 2.5 ). Due to the consideration of aerosol direct and indirect effects, two-way WRF-CMAQ shows improved performance over offline coupled WRF and CMAQ in terms of spatiotemporal distributions and statistics, especially for radiation, cloud forcing, O 3 , sulfate, nitrate, ammonium, elemental carbon, tropospheric O 3 residual, and column nitrogen dioxide (NO 2 ). For example, the mean biases have been reduced by more than 10 W m –2 for shortwave radiation and cloud radiative forcing and by more than 2 ppb for max 8 h O 3 . However, relatively large biases still exist for cloud predictions, some PM 2.5 species, and PM 10 that warrant follow-up studies to better understand those issues. The impacts of chemistry–meteorological feedbacks are found to play important roles in affecting regional air quality in the US by reducing domain-average concentrations of carbon monoxide (CO), O 3 , nitrogen oxide (NO x ), volatile organic compounds (VOCs), and PM 2.5 by 3.1 % (up to 27.8 %), 4.2 % (up to 16.2 %), 6.6 % (up to 50.9 %), 5.8 % (up to 46.6 %), and 8.6 % (up to 49.1 %), respectively, mainly due to reduced radiation, temperature, and wind speed. The overall performance of the two-way coupled WRF-CMAQ model achieved in this work is generally good or satisfactory and the improved performance for two-way coupled WRF-CMAQ should be considered along with other factors in developing future model applications to inform policy making.

58 GEOSCIENCES↗

DOE OSTI2017

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.

58 GEOSCIENCES↗

DOE OSTI2016

Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model

We address the inverse problem of inferring the basal geothermal heat flux from surface velocity observations using a steady-state thermomechanically coupled nonlinear Stokes ice flow model. This is a challenging inverse problem since the map from basal heat flux to surface velocity observables is indirect: the heat flux is a boundary condition for the thermal advection–diffusion equation, which couples to the nonlinear Stokes ice flow equations; together they determine the surface ice flow velocity. This multiphysics inverse problem is formulated as a nonlinear least-squares optimization problem with a cost functional that includes the data misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to render the problem well posed. We derive adjoint-based gradient and Hessian expressions for the resulting partial differential equation (PDE)-constrained optimization problem and propose an inexact Newton method for its solution. As a consequence of the Petrov–Galerkin discretization of the energy equation, we show that discretization and differentiation do not commute; that is, the order in which we discretize the cost functional and differentiate it affects the correctness of the gradient. Using two- and three-dimensional model problems, we study the prospects for and limitations of the inference of the geothermal heat flux field from surface velocity observations. The results show that the reconstruction improves as the noise level in the observations decreases and that short-wavelength variations in the geothermal heat flux are difficult to recover. We analyze the ill-posedness of the inverse problem as a function of the number of observations by examining the spectrum of the Hessian of the cost functional. Motivated by the popularity of operator-split or staggered solvers for forward multiphysics problems – i.e., those that drop two-way coupling terms to yield a one-way coupled forward Jacobian – we study the effect on the inversion of a one-way coupling of the adjoint energy and Stokes equations. Here, we show that taking such a one-way coupled approach for the adjoint equations can lead to an incorrect gradient and premature termination of optimization iterations. This is due to loss of a descent direction stemming from inconsistency of the gradient with the contours of the cost functional. Nevertheless, one may still obtain a reasonable approximate inverse solution particularly if important features of the reconstructed solution emerge early in optimization iterations, before the premature termination.

15 GEOTHERMAL ENERGY↗

NASA NTRS1976

Stability of elastic bending and torsion of uniform cantilever rotor blades in hover with variable structural coupling

The stability of elastic flap bending, lead-lag bending, and torsion of uniform, untwisted, cantilever rotor blades without chordwise offsets between the elastic, mass, tension, and areodynamic center axes is investigated for the hovering flight condition. The equations of motion are obtained by simplifying the general, nonlinear, partial differential equations of motion of an elastic rotating cantilever blade. The equations are adapted for a linearized stability analysis in the hovering flight condition by prescribing aerodynamic forces, applying Galerkin's method, and linearizing the resulting ordinary differential equations about the equilibrium operating condition. The aerodynamic forces are obtained from strip theory based on a quasi-steady approximation of two-dimensional unsteady airfoil theory. Six coupled mode shapes, calculated from free vibration about the equilibrium operating condition, are used in the linearized stability analysis. The study emphasizes the effects of two types of structural coupling that strongly influence the stability of hingeless rotor blades. The first structural coupling is the linear coupling between flap and lead-lag bending of the rotor blade. The second structural coupling is a nonlinear coupling between flap bending, lead-lag bending, and torsion deflections. Results are obtained for a wide variety of hingeless rotor configurations and operating conditions in order to provide a reasonably complete picture of hingeless rotor blade stability characteristics.

Hodges, D. H., Roberta.↗