SEARCH · Search NASA

Results for “model results”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Reducing impacts of artificial ponding in modeling salt marshes using a conductivity-formulated subgrid model

The landforms of salt marsh systems are rather complex, with meandering channel networks cutting through low-lying, extensive marsh platforms. Some small-scale topographic details are critical to flooding and draining processes but cannot be reflected in the Digital Elevation Model (DEM) due to limited measurement resolution and vegetation bias. Numerical models need a grid with high enough resolution to capture the small-scale connectivity to correctly model the small-scale flow processes, making the model computationally unaffordable. Here, in this study, we derived a subgrid model, which couples two flow components, one is the regular subgrid model for simulating the subgrid-resolved flows, and the other is a conductivity-formulated porous flow model for modeling small-scale gut flows unresolved at the subgrid level. The model validation was performed using an idealized case with a narrow slot that was unresolved in the coupled flow model. Model results were compared with a slot-resolved regular subgrid model where only surface flows are solved and suggested that, with a proper model parameterization, the coupled flow model is capable of reproducing flooding and draining processes similar to the slot-resolved subgrid model of surface flow. The model was applied to simulating flooding and draining processes in a meso-tidal salt marsh. The model results demonstrate that the coupled subgrid model can be used to reduce the impacts of artificial ponding that usually occur in numerical simulations of salt marshes with numerous narrow gullies and creeks.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Multiscale Modeling of the Mechanical Response of Silicon Carbide Composite Within the Accelerated Fuel Qualification Framework

The accelerated fuel qualification (AFQ) framework has been used for the initial development of multiscale modeling of silicon carbide (SiC) fiber reinforced composite (SiC-SiC). The AFQ framework provides a methodology to leverage physics-informed multiscale modeling along with a reduced set of empirical test data to reduce the time and cost of licensing and qualification of new nuclear fuel systems while maintaining the overall nuclear power plant safety case. SiC-SiC is being proposed for in-core applications, most notably fuel cladding, for current and next-generation nuclear reactors because of its high temperature stability, irradiation tolerance, and ability to withstand many accident conditions. As these composites exhibit multiscale architectures and complex microstructure-based fracture mechanics, it is an appealing use case for the AFQ methodology. While the end goal of this work is a single multiscale model that can be used for predictive in-core performance, current focus is on the individual various length scale models. Four individual models have been initially developed from microscale to engineering system level to capture key physics-based effects across different length scales. These models include a microscale homogenized tow model, a mesoscale fast Fourier transform–based weave model that integrates the homogenized tow model, a mesoscale finite element–based weave model, and a system-level BISON fuel performance model. Results of these models have undergone an initial comparison with separate-effects test data showing a good match to experimental results. By using the AFQ framework during model development, several near-term benefits have been secured including a reduction in development time for the SiC-SiC cladding, more targeted irradiation testing, and a better understanding of uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Modeling Ring Current Proton Fluxes Using Artificial Neural Network and Van Allen Probe Measurements

Abstract Terrestrial ring current dynamics are a critical part of the near‐space environment, in that they directly drive geomagnetic field variations that control particle drifts, and define geomagnetic storms. The present study aims to specify a global and time‐varying distribution of ring current proton using geomagnetic indices and solar wind parameters with their history as input. We train an artificial neural network (ANN) model to reproduce proton fluxes measured by the Radiation Belt Storm Probes Ion Composition Experiment instrument onboard Van Allen Probes. By choosing optimal feature parameters and their history length, the model results show a high correlation and a small error between model specifications and satellite measurements. The modeled results well capture energy‐dependent proton dynamics in association with geomagnetic storms, including inward radial diffusion, acceleration and decay. Our ANN model produces proton fluxes with their corresponding 3D spatiotemporal variations, capturing the latitudinal distribution and local time asymmetry that are consistent with observations and that can further inform theory.

Li, Jinxing↗

DOE OSTI2021

Polarimetric Radar Variables in Snowfall at Ka- and W-Band Frequency Bands: A Comparative Analysis

Dual-frequency millimeter-wavelength radar observations in snowfall are analyzed in order to evaluate differences in conventional polarimetric radar variables such as differential reflectivity (ZDR) specific differential phase shift (K DP ) and linear depolarization ratio (LDR) at traditional cloud radar frequencies at Ka and W bands (~35 and ~94 GHz, correspondingly). Low radar beam elevation (~5°) measurements were performed at Oliktok Point, Alaska, with a scanning fully polarimetric radar operating in the horizontal–vertical polarization basis. This radar has the same gate spacing and very close beam widths at both frequencies, which largely alleviates uncertainties associated with spatial and temporal data matching. It is shown that observed Ka- and W-band Z DR differences are, on average, less than about 0.5 dB and do not have a pronounced trend as a function of snowfall reflectivity. The observed Z DR differences agree well with modeling results obtained using integration over nonspherical ice particle size distributions. For higher signal-to-noise ratios, KDP data derived from differential phase measurements are approximately scaled as reciprocals of corresponding radar frequencies indicating that the influence of non-Rayleigh scattering effects on this variable is rather limited. This result is also in satisfactory agreement with data obtained by modeling using realistic particle size distributions. Observed Ka- and W-band LDR differences are strongly affected by the radar hardware system polarization “leak” and are generally less than 4 dB. Smaller differences are observed for higher depolarizations, where the polarization “leak” is less pronounced. Realistic assumptions about particle canting and the system polarization isolation lead to modeling results that satisfactorily agree with observational dual-frequency LDR data.

42 ENGINEERING↗

DOE OSTI2023

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Long-term sinking of nuclear waste canisters in salt formations by low-stress creep at high temperature

Abstract Rock salt has a self-sealing capacity, low permeability, and high thermal conductivity, making it a potential host for heat-generating nuclear waste. The feasibility of nuclear waste disposal within salt formations has been investigated mostly for small-sized canisters. Geologic disposal of larger-sized canisters originally designed for spent fuel storage and transportation has lately been examined as a cost-effective alternative. This raises questions about their long-term vertical movement due to their weight and high decay heat. Low-stress creep governs this movement; however, most salt constitutive models do not incorporate it. In this paper, the Norton and the WIPP creep models are compared with the Lux/Wolters/Lerche (LWL) model and a simpler model that combines linear and Norton creep laws (named combined creep model). The LWL and combined creep models consider pressure solution creep, though all incorporate dislocation creep. The models are first applied to creep tests under various stress levels. The LWL and the combined creep models results fit the experimental data well in both high and low stress ranges, whereas the Norton and WIPP models results only fit in higher-stress ranges. The different models are further applied for analyzing long-term canister movement. A sinking rate of $$-4.4\times 10^{-7}$$ - 4.4 × 10 - 7 mm/year was predicted using the Norton and WIPP models versus $$-2.1\times 10^{-2}$$ - 2.1 × 10 - 2 mm/year and $$-3.1\times 10^{-2}$$ - 3.1 × 10 - 2 mm/year using the LWL and the combined creep models, respectively. This comparative study confirms that creep models calibrated exclusively against high-deviatoric stress data might result in an inaccurate estimation of waste packages sinking rate in salt formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Marmot

Marmot is a data formatting and visualization tool for production cost modelling results. It provides an efficient way to view PLEXOS production cost modeling results quickly, while also creating publication ready figures and data tables.

Levie, Daniel↗

DOE OSTI2020

Effect of Particle Spreading Dynamics on Powder Bed Quality in Metal Additive Manufacturing

Powder spreading precedes creation of every new layer in powder bed additive manufacturing (AM). The powder spreading process can lead to powder layer defects such as porosity, poor surface roughness and particle segregation. Therefore, the creation of homogeneous layers is the first task for optimal part printing. Discrete element methods (DEM) powder spreading simulations are typically limited to a single layer and/or small number of particles. Therefore, results from such model configurations may not be generalized to multiple layer processes. In this study, a computationally efficient multi-layer powder spreading DEM simulation model is proposed. The model is calibrated experimentally using static Angle of Repose measurements. The adhesion model parameter, cohesive energy density is related to adhesive surface energy and strain energy release rate parameters. The model results show that interaction between particle and the powder spreading rake leads to noticeable variation in packing density, surface roughness, dynamic angle of repose (AOR), particle size distribution, and particle segregation. Finally, the powder model is experimentally validated using a recoater spreading rig to measure the dynamic AOR at spreading speeds consistent with recoating speeds and layer heights used in AM processes.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Herbaceous Feedstock 2019 State of Technology Report

The U.S. Department of Energy (DOE) promotes the production of advanced liquid transportation fuels from lignocellulosic biomass by funding fundamental and applied research that advances the State of Technology (SOT). As part of its involvement with this mission, Idaho National Laboratory (INL) completes an annual SOT report for biomass feedstock logistics. This report summarizes supply system impacts of Bioenergy Technologies Office (BETO)-funded research and development efforts at INL and elsewhere (such as the High-Tonnage Feedstock Logistics projects (Webb et al. 2013a, Webb et al. 2013b, Webb et al. 2013c, Webb and Sokhansanj 2014, Sokhansanj et al. 2014)) that lead to improvements in feedstock supply systems. These include improvements to and observed performance of innovative harvest and collection methods, storage technologies, transportation and handling approaches, and advanced preprocessing technologies. Biomass quality and variability, and the interface between feedstock quality and conversion performance are key drivers in addition to delivered feedstock cost. In this report, we estimate the benefits of R&D improvements to individual supply system unit operations and present the status of feedstock logistics technology development for converting biomass into biofuels. These analyses are supported by experimental data where possible and help to align the SOT relative to the cost goals defined in the Multi-Year Program Plan. The 2019 Herbaceous SOT incorporates several technology changes in feedstock preprocessing and introduces opportunities from the integrated landscape management (ILM) strategy and increased grower participation to reduce biomass access costs, while maintaining or improving grower profitability. During FY18 uneven flow from the horizontal bale grinder was identified as a significant issue limiting preprocessing system throughput. Based on FSL-funded research at INL, the 2019 Herbaceous SOT replaces the horizontal bale grinder used in the first stage size reduction with a bale processor. The improved uniformity of biomass flow entering the PDU eliminated slugging flow from the first stage size reduction and improved the throughput of downstream operations. In order to achieve moisture reduction through frictional heating during grinding (which allowed elimination of the costly rotary drum dryer in previous SOTs), the second stage grinder was changed from a rotary shear, which does not remove moisture, back to a hammer mill. Finally, the 2019 Herbaceous SOT introduces modified three-pass and two-pass corn stover supply curves derived from the BT16 resource assessment, based on FY19 modeling results (WBS 4.2.1.20) quantifying economic benefits of ILM in the supply area, together with modeling results (WBS 1.2.1.5) identifying ILM strategies to increase grower participation. The 2019 Herbaceous SOT report documents the current modeled cost of an herbaceous feedstock supply system from harvest to the pretreatment reactor throat for hydrocarbon fuel production via biochemical conversion, based on equipment and processes now available or potentially available in the near term. The modeled cost also considers both the required quality and the availability of the biomass resources. The 2019 Herbaceous SOT predicts a modeled delivered feedstock cost of $81.37 /dry ton (2016$); this is a $2.30/dry ton (2016$) decrease from the 2018 Herbaceous SOT. Technology improvements that contributed to this modeled cost reduction include reduced cost for the new preprocessing design and quantification of the opportunities of the integrated landscape management (ILM) strategy and an increased grower participation rate to reduce the grower payment portion of biomass access costs, while maintaining or improving grower profitability. A greenhouse gas emissions (GHG) assessment was completed by Argonne National Laboratory using the 2019 Greenhouse Gases, Regulated Emissions, and Energy use in Transportation model, estimating an increase of 14.89 kg CO2e/ton from the 2018 SOT (69.27 kg CO2e/ton in 2018 to 84.16 kg CO2e/ton in 2019). The increase of energy consumption during preprocessing along with higher transportation distance to access low cost biomass from further distance contributed to the increase of GHG emissions in the 2019 Herbaceous SOT. The reason for the increased transportation distances was the cost tradeoff of going farther from the biorefinery to access the cheaper ILM-derived counties (the cheaper price outweighed the cost of increased supply radius).

09 BIOMASS FUELS↗

DOE OSTI2020

Clustering and α-capture reaction rate from ab initio 20 Ne

We introduce a new framework for studying clustering and for calculating α partial widths using ab initio wave functions. We demonstrate the formalism for 20 Ne, by calculating the overlap between the 16 O + α cluster configuration and states in 20 Ne computed in the ab initio symmetry-adapted no-core shell model. We present spectroscopic amplitudes and spectroscopic factors, and compare those to no-core symplectic shell-model results in larger model spaces, to gain insight into the underlying physics that drives α clustering. Specifically, we report on the α partial width of the lowest 1 – resonance in 20 Ne, which is found to be in good agreement with experiment.We also present first no-core shell-model estimates for asymptotic normalization coefficients for the ground state, as well as for the first excited 4 + state in 20 Ne that lies in a close proximity to the α + 16 O threshold. This outcome highlights the importance of correlations for developing cluster structures and for describing α widths. The widths can then be used to calculate α-capture reaction rates for narrow resonances of interest to astrophysics. Finally, we explore the reaction rate for the α-capture reaction 16 O(α, γ ) 20 Ne at astrophysically relevant temperatures and determine its impact on simulated x-ray burst abundances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

From Data to Insights: A Covariate Analysis of the IARPA BRIAR Dataset for Multimodal Biometric Recognition Algorithms at Altitude and Range

This paper examines covariate effects on fused whole body biometrics performance in the IARPA BRIAR dataset, specifically focusing on UAV platforms, elevated positions, and distances up to 1000 meters. The dataset includes outdoor videos compared with indoor images and controlled gait recordings. Normalized raw fusion scores relate directly to predicted false accept rates (FAR), offering an intuitive means for interpreting model results. A linear model is developed to predict biometric algorithm scores, analyzing their performance to identify the most influential covariates on accuracy at altitude and range. Weather factors like temperature, wind speed, solar loading, and turbulence are also investigated in this analysis. The study found that resolution and camera distance best predicted accuracy and findings can guide future research and development efforts in long-range/elevated/UAV biometrics and support the creation of more reliable and robust systems for national security and other critical domains.

Bolme, David↗

DOE OSTI2020

Analytic Thermal Model of an Optical Fiber Based Gamma Thermometer and its Application in a University Research Reactor

This paper describes and validates, by comparison with numerical modeling results, an analytical model of thermal transport in an optical fiber based gamma thermometer (OFBGT) that is appropriate for use in university research reactors. The maximum temperature difference between the thermal mass and the outer sheath ( Δ T ) for the OFBGT design that we have considered is approximately 50 ∘ C , for the OFBGT in the Central Irradiation Facility of the Ohio State University Research Reactor (OSURR) with the reactor operating at full power (450 kW). The maximum value of Δ T that is predicted by the analytic model for the OFBGT design is smaller by approximately 1.1 °C than the value of Δ T which is predicted by the numerical model, for the same OFBGT design, but including all the details of the design. We have used the analytical model of thermal transport in an OFBGT to determine a normalized Modulation Transfer Function M T F ′ ( k ) ) for the OFBGT. We conclude that MTF ′ ( k OSU ) > 0.99 , where k OSU is the spatial frequency for the axial dependence of the reactor power distribution in the

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Experimental Setup and Learning-Based AI Model for Developing Accurate PV Inverter Models [Slides]

The integration of power electronics-based interfaces presents challenges due to the absence of detailed models and the high computational complexity. Generic models used in system studies lack accuracy in capturing converter dynamics. This paper proposes a data-driven approach developed from experimental setup data. This approach enhances accuracy in photovoltaic inverter modeling. We used two types of PV inverters in the experiment. The recorded experimental data undergo processing through a machine learning model. Results from the model trained through machine learning is also presented.

14 SOLAR ENERGY↗

DOE OSTI2023

Continental Scale Hydrostratigraphy: Basin‐Scale Testing of Alternative Data‐Driven Approaches

Abstract Integrated hydrological modeling is an effective method for understanding interactions between parts of the hydrologic cycle, quantifying water resources, and furthering knowledge of hydrologic processes. However, these models are dependent on robust and accurate datasets that physically represent spatial characteristics as model inputs. This study evaluates multiple data‐driven approaches for estimating hydraulic conductivity and subsurface properties at the continental‐scale, constructed from existing subsurface dataset components. Each subsurface configuration represents upper (unconfined) hydrogeology, lower (confined) hydrogeology, and the presence of a vertical flow barrier. Configurations are tested in two large‐scale U.S. watersheds using an integrated model. Model results are compared to observed streamflow and steady state water table depth (WTD). We provide model results for a range of configurations and show that both WTD and surface water partitioning are important indicators of performance. We also show that geology data source, total subsurface depth, anisotropy, and inclusion of a vertical flow barrier are the most important considerations for subsurface configurations. While a range of configurations proved viable, we provide a recommended Selected National Configuration 1 km resolution subsurface dataset for use in distributed large‐and continental‐scale hydrologic modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Correction to “Frank–Kasper Phases of Diblock Copolymer Melts Studied with the DPD Model: SCF Results”

In a recent study, (1) we performed the self-consistent field (SCF) calculations of a dissipative particle dynamics (DPD) model of conformationally asymmetric diblock copolymer A-B melts to study the stability of various Frank–Kasper (FK) phases formed by such A-B melts. While all the equations in ref (1) are correct, we did not calculate the dimensionless (mean-field) Helmholtz free energies per chain βf c = βu c,χ + βu c,κ – s c /k B and the dimensionless stresses dβf c /dθ of ordered phases accurately. For example, the dimensionless (nonbonded) internal energy per chain due to the system compressibility βu c,κ = (N/2κ)∑ q βû 0 (|q|)(ϕ̂ A (q) + ϕ̂ B (q) – δ q,0 )(ϕ̂ A (−q) + ϕ̂ B (−q) – δ– q,0 ) should be calculated in the reciprocal space as given (note that βû 0 (|q|) is obtained analytically, and that since ϕ A (r) is real 𝜙̂ A (−𝐪)= $\overline{𝜙̂ A}$(𝐪), i.e., the complex conjugate of ϕ̂ A (q)); this corresponds to numerically evaluating βu c,κ = (N/2κ)∫dr(ϕ A (r) + ϕ B (r) – 1)ψ(r) with ψ(r) ≡ (∫dr′βu 0 (|r–r′|)(ϕ A (r′) + ϕ B (r′) – 1)/V in the real space via the (composite) trapezoidal rule. While the trapezoidal rule is less accurate than the Romberg integration (2) in general, the opposite occurs when the integration domain is periodic; this was pointed out by Matsen in ref (3) and explained in detail in refs (4and5). With our ignorance of this exceptional convergence property of the trapezoidal rule for periodic functions, we calculated all convolutions having the form of ψ(r) in the reciprocal space via the fast Fourier transforms (6) but those having the form of βu c,κ in the real space via Romberg integration; the latter leads to inaccurate results and unfortunately incorrect phase diagrams reported in ref (1), which are corrected here.

He, Juntong↗

DOE OSTI2021

How structural differences influence cross-model consistency: An electric sector case study

Multiple models are often employed to describe a range of possible outcomes for one or more scenarios, yielding insights into causal relationships and their uncertainties. Electric sector capacity expansion scenarios are a common topic of such efforts due to the economic influence of the electric sector, but model results typically span a broad solution space despite efforts to harmonize input assumptions, making decision implications difficult to discern. This study investigates the relationship between input harmonization and cross-model scenario consistency under disparate electric sector scenarios. We compare cross-model consistency between two state-of-the-art electric sector capacity expansion models (GCAM-USA and ReEDS) for six electric sector scenarios comprising alternate assumptions about fossil fuel resource availability, technology innovation, and long-term economy-wide transitions under four harmonization configurations varying model representations of electricity demand, fuel prices, renewable resources, and capacity retirements. These comparisons reveal that cross-model consistency can vary across scenarios under a given harmonization configuration, suggesting that harmonization efforts must often be scenario-specific if comparable cross-model consistency is desired. Model structural differences can hinder consistency, and the impact of these differences can depend on the scenario. Ultimately, thorough harmonization can reveal insights into cross-model consistency, which can be used to tighten uncertainty bounds and improve the decision-making implications of multi-model activities.

Cohen, Stuart↗

DOE OSTI2021

A data-driven peridynamic continuum model for upscaling molecular dynamics

Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials’ behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.

homogenization↗