A fusion of global-scale, river-centric datasets for rapid analysis and modeling frameworks: Veins of the Earth
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Objectives/Scope: Wave resource characterization is a critical step for wave energy converter deployment in the coastal ocean and relies on long-term, high-resolution wave datasets. This study presents a detailed modeling study of the wave resource along the U.S. West Coast (Washington, Oregon, and California), a coastal region that was identified with high wave energy potential in earlier studies. Methods, Procedures, Process: The wave hindcast covers a 32-year period from 1979 to 2010 and is based on a multi-resolution, unstructured-grid SWAN model framework. Model configuration closely follows and meets the requirements recommended by the International Electrotechnical Commission Technical Specification (IEC TS) for wave energy resource assessment and characterization (Class 2 - feasibility study). The model domain covers the entire U.S. Exclusive Economic Zone (EEZ) in the West Coast and has a spatial resolution varying from ~300 m in the nearshore region (20 km from the shoreline) to ~2500 m within the EEZ and ~5000 m at the open boundary, which extends beyond the EEZ. The model was forced by hourly 2-D wave spectra produced by a two-way nested WaveWatch III model, which covers the global ocean domain and the broader U.S. West Coast region domain with spatial resolutions of 0.5 degree and 10 arc-minutes, respectively. Both wave models are forced by hourly, 0.5-degree wind forcing obtained from NCEP’s Climate Forecast System Reanalysis (CFSR) product. Results, Observations, Conclusions: The standard model output for the SWAN model includes 3-hourly output for the six IEC wave resource parameters (e.g., omnidirectional wave power) at each grid point and hourly 2-D spectra at more than 50 NDBC buoys. Extensive model validation was achieved by comparing the six model-predicted IEC parameters with those derived from field observations at representative NDBC buoys. The error statistics indicated the model’s satisfactory performance. Further analyses were conducted to systematically evaluate the temporal and spatial distributions of wave energy potential and wave climate along the U.S. West Coast. Results suggest that Washington and Oregon coasts have similar nearshore wave resource, which is significantly higher than resources in Southern California. Strong seasonal variations are also observed, e.g., high wave energy tends to occur in the winter months. In summary, this study produced the first high-resolution, comprehensive dataset on wave energy distribution along the U.S. West Coast. Novel/Additive Information: The results are being used by the National Renewable Energy Laboratory to update the MHK Atlas, which was originally derived from NOAA’s 4-arc-minute WaveWatch III model output. In addition, the monthly averaged wave energy climatology dataset can be readily shared to support a variety of research and application efforts within the EEZ of the U.S. West Coast.
Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20𝑛 0 . The uncertainties of the 𝜒EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.
Explore the source record for details and available documents.
Energy modeling frameworks and scenario analysis help us explore the potential impact of our actions and uncertainties in future energy systems. Despite their importance, there is no systematic procedure for evaluating the scenario development process. In a literature review, we identify two core elements of the scenario development process: internal consistency and diversity which are oftentimes missing from scenarios. Here, to address this gap, we create the Internal consistency and Diversity Scenario Development (IDSD) comparative framework which aims to assess the feasibility and diversity of scenarios for a given energy model. With this framework, we review commonly used energy models and demonstrate our framework on their scenarios. The IDSD comparative framework can serve several purposes absent from previous scenario development work by aiding energy modelers and report writers in crafting high-quality scenarios. First, the IDSD is a reflective tool which can improve the quality of the scenario development process, enabling a comparative assessment of energy models and scenarios. Second, the IDSD can provide guidance to modeling frameworks with existing scenarios and those still in development; this feedback will enable modelers to improve the development and the communication of the limitations of their scenarios. Third, this study has highlighted areas for improvement in the scenario development of some commonly used energy model frameworks. Finally, there is a complete lack of explanation regarding the stakeholder selection process. Addressing these identified items could increase opportunities for advanced energy technology uptake and improve our options for achieving a more resilient energy system.
Anionic molecular models for nonhydrolyzed and partially hydrolyzed aluminum and gallium framework sites on silica, M[OSi(OtBu) 3 ] 4 - and HOM[OSi(OtBu) 3 ] 3 - (where M=Al or Ga), were synthesized from anionic chlorides Li{M[OSi(OtBu) 3 ] 3 Cl} in salt metathesis reactions. Sequestration of lithium cations with [12]crown-4 afforded charge-separated ion pairs composed of monomeric anions M[OSi(OtBu) 3 ] 4 - with outer-sphere [([12]crown-4) 2 Li] + cations, and hydroxides {HOM[OSi(OtBu) 3 ] 3 } with pendant [([12]crown-4)Li] + cations. These molecular models were characterized by single-crystal X-ray diffraction, vibrational spectroscopy, mass spectrometry and NMR spectroscopy. Upon treatment of monomeric [([12]crown-4)Li]{HOM[OSi(OtBu) 3 ] 3 } complexes with benzyl alcohol, benzyloxide complexes were formed, modeling a possible pathway for the formation of active sites for Meerwin-Ponndorf-Verley (MPV) transfer hydrogenations with Al/Ga-doped silica catalysts.
We propose a vector alternative to the typical raster based population modeling framework. When compared with rasters, vectors are more precise, have the ability to hold more information, and are more conducive to areal constructs such as building and parcel outlines. While rasters have traditionally provided computational efficiency, much of this efficiency is reduced at finer resolutions and computational resources are more plentiful today. Herein we describe the approach and implementation methodology. We also describe the output data stack for the United States and provide examples and applications.
With support from the Department of Energy’s Water Power Technologies Office, Pacific Northwest National Laboratory (PNNL) has developed new tools that incorporate cutting-edge climate and hydrological science capabilities to assess the potential long-term impacts of future climate conditions on unregulated streamflow and water temperature within watershed-river-reservoir systems. The objectives of this project were achieved by enhancing key hydrologic and hydrodynamic models and transferring them to a high-performance computing environment to provide a high-spatiotemporal resolution, multi-scale modeling framework. The new modeling framework has the potential to quantify risks related climate change impacts on runoff, unregulated streamflow, and water temperature. Initial development and demonstration of the modeling framework was conducted under historical and future climate conditions in the Columbia River Basin in the Pacific Northwest and the Connecticut River Basin in New England.
The James Webb Space Telescope (JWST) and Atacama Large Millimeter/submillimeter Array have detected emission lines from the ionized interstellar medium (ISM) in some of the first galaxies at z ≳ 6. These measurements present an opportunity to better understand galaxy assembly histories and may allow important tests of state-of-the-art galaxy formation simulations. It is challenging, however, to model these lines in their proper cosmological context. In order to meet this challenge, we introduce a novel subgrid line emission modeling framework. The framework uses the high-z zoom-in simulation suite from the Feedback in Realistic Environments (FIRE) collaboration. The line emission signals from H II regions within each simulated FIRE galaxy are modeled using the semianalytic HIIL INES code. A machine learning approach is then used to determine the conditional probability distribution for the line luminosity to stellar-mass ratio from the H II regions around each simulated stellar particle. This conditional probability distribution can then be applied to predict the line luminosities around stellar particles in lower-resolution, yet larger volume cosmological simulations. As an example, we apply this approach to the IllustrisTNG simulations at z = 6. The resulting predictions for the [O II ], [O III ], and Balmer line luminosities as a function of star formation rate agree well with current observations. Our predictions differ, however, from related works in the literature, which lack detailed subgrid ISM models. This highlights the importance of our multiscale simulation modeling framework. Finally, we provide forecasts for future line luminosity function measurements from the JWST and quantify the cosmic variance in such surveys.
The development of battery electric vehicles (BEVs) is accelerating due to their environmental advantages over gasoline and diesel-powered vehicles, including a decrease in air pollution and an increase in energy efficiency. The deployment of charging infrastructure will need to increase to keep pace with demand, especially for large commercial vehicles for which few public chargers currently exist. In this paper, a new flexible framework is proposed for optimizing the placement of charging stations for BEVs, within which different physical models and optimization techniques may be used. Furthermore, a set of metrics is suggested to help enforce complex constraints and facilitate direct comparison between different optimization techniques. Unlike many existing charger placement techniques, the proposed method directly considers the historical driving patterns on a vehicle-by-vehicle basis, using transparent models to assess impacts of candidate charger placements, thus improving the explainability of the results. In the developed framework, modeled BEVs are first generated along the road network to mimic historical traffic data and are simulated traveling along a given route according to a simplified vehicle model. During the simulation, the charger placement problem is initially relaxed to allow vehicles to charge at any node along the road network, and vehicle states are tracked to assess areas of high charging demand. Charging stations are then placed based on the results of the relaxed simulation, and suggested placements are evaluated via road network simulation with fixed charger locations. This proposed framework is applied to a sample problem of placing charging stations along five major highway corridors for Class 8 over-the-road electric trucks. A novel mixed integer programming (MIP) formulation is proposed to optimize charger placements based upon the expected charging demand. Constraints were imposed on the final placement results to limit expected wait times at each station and ensure a minimum threshold of trucking routes are viable for BEVs. The results demonstrate the flexibility and potential effectiveness of the developed model-based framework for scalable charger station deployment.
Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. Here in this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.
State-dependent models can be used to represent the system recovery process as a series of stochastic transitions from lower to higher functional states. However, the applications of these models have been limited in scope and there is a lack of a generalized recovery modeling framework. A generalized framework would permit a robust forecasting of systems and system-of-systems recovery under multiple hazards, and more broadly, would contribute to community disaster preparedness. This paper develops a generalized post hazard-event recovery modeling framework based on state-dependent Markov-type processes. We then apply the proposed framework to solve a spectrum of problems that range from hind-casting single-system recovery following a single hazard event to forecasting post-event trajectories under multiple hazards and modeling the recovery of a system-of-systems. First, Markov chains are used to hind-cast the observed recovery for a portfolio of buildings affected by the 2014 South Napa, California, earthquake. Next, Markov processes are used to formulate a parametric post hazard-event recovery model, which can be updated using Bayesian statistics when relevant datasets become available. Semi-Markov processes are then used to develop a more general model of single hazard recovery, which accounts for the intensity of the loading and level of damage caused by the event. Semi-Markov processes with non-renewal features are then used to account for multihazard interactions in a post-event recovery model, and applied to a case study that involves a community in Charleston, South Carolina. Lastly, Markov-type processes are combined with Bayesian networks to model the recovery of residential, commercial, educational, and industrial buildings (system-of-systems) following a hazard event. Overall, these applications demonstrate the versatility of the Markov framework towards handling recovery problems with varying levels of complexity.
To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.
Abstract The constitutive relations of the Richardson‐Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free‐form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson‐Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network‐based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.
Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.