Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “impact modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High‐Latitude Stratospheric Aerosol Injection to Preserve the Arctic

Abstract Stratospheric aerosol injection (SAI) has been shown in climate models to reduce some impacts of global warming in the Arctic, including the loss of sea ice, permafrost thaw, and reduction of Greenland Ice Sheet (GrIS) mass; SAI at high latitudes could preferentially target these impacts. In this study, we use the Community Earth System Model to simulate two Arctic‐focused SAI strategies, which inject at 60°N latitude each spring with injection rates adjusted to either maintain September Arctic sea ice at 2030 levels (“Arctic Low”) or restore it to 2010 levels (“Arctic High”). Both simulations maintain or restore September sea ice to within 10% of their respective targets, reduce permafrost thaw, and increase GrIS surface mass balance by reducing runoff. Arctic High reduces these impacts more effectively than a globally focused SAI strategy that injects similar quantities of SO 2 at lower latitudes. However, Arctic‐focused SAI is not merely a “reset button” for the Arctic climate, but brings about a novel climate state, including changes to the seasonal cycles of Northern Hemisphere temperature and sea ice and less high‐latitude carbon uptake relative to SSP2‐4.5. Additionally, while Arctic‐focused SAI produces the most cooling near the pole, its effects are not confined to the Arctic, including detectable cooling throughout most of the northern hemisphere for both simulations, increased mid‐latitude sulfur deposition, and a southward shift of the location of the Intertropical Convergence Zone. For these reasons, it would be incorrect to consider Arctic‐focused SAI as “local” geoengineering, even when compared to a globally focused strategy.

Lee, Walker Raymond↗
DOE OSTI2021

Electrification of Pakistan's Transport System: Modeling Electric Vehicle Penetration and Energy Supply Chain Impacts

Initial analysis of the National Electric Vehicle Policy (NEVP) of Pakistan was undertaken by an integrated energy planning team from November 2019 through August 2020. This technical reference document seeks to present the background, inputs and assumptions, methodology, and results of the policy analysis in detail. The intended audience is the technical modeler, analyst, or reviewer who seeks to understand specifically the modeling effort within this initial analytical phase, with the end-goal of interpreting, re-producing, or modifying the simulations or extending the models to conduct follow-on studies subsequently. These studies may examine in greater detail the energy sector, impacts to emissions, charging station infrastructure needs, and overall benefits-cost trade of the policy, among other areas of interest.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Impact of flow-induced beam deflection on beam propagation in ignition scale hohlraums

Experiments examining the amount of specular reflection (or “glint”) within hohlraums containing different gas fill densities have recently been performed. Simulations of these experiments are presented that show using a single flux limiter cannot explain the decrease in glinted power with increasing gas fill density. The hypothesis that flow-induced beam deflection alters laser absorption is presented. A model is proposed that can be implemented into a ray tracing description of the laser commonly used in radiation hydrodynamic codes. It is shown that simulations using this model capture the trend with gas fill density improving agreement with measurements. This formulation is then applied to an ensemble of laser-driven inertial confinement fusion experiments performed at the National Ignition Facility. The proposed model shows little impact on the total x-ray drive on the capsule but a large impact on the resulting implosion symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Physically-based control-oriented modeling for turbocharged stoichiometric spark-ignited engine with cooled EGR and flexible VVT systems

Accurate estimation and prediction of engine gas exchange system and in-cylinder conditions are critical for spark-ignited engine control and diagnostic algorithm development. In this paper, a physically-based, control-oriented model for a 2.8 l turbocharged, variable valve timing (VVT) and low pressure (LP) exhaust gas recirculation (EGR)-utilizing SI engine was developed. The model includes the impact of modulation to any combination of 10 actuators, including the throttle valve, compressor bypass valve, fueling rate, waste-gate, LP EGR valve, number of deactivated cylinders, intake valve open (IVO) timing, intake valve close (IVC) timing, exhaust valve open (EVO) timing and exhaust valve close (EVC) timing. The accuracy of the model in capturing engine dynamics was demonstrated by validating it against high-fidelity engine GT-Power simulation results for various drive cycles, particularly emphasizing elevated loads. In comparison to the open literature, novel contributions of the effort described in this paper includes in-cylinder gas composition modeling and turbine-out pressure estimation.

Zhang, Xu↗
DOE OSTI2024

A multi-model assessment of the early last deglaciation (PMIP4 LDv1): a meltwater perspective

Abstract. The last deglaciation (∼20–11 ka BP) is a period of a major, long-term climate transition from a glacial to interglacial state that features multiple centennial- to decadal-scale abrupt climate variations whose root cause is still not fully understood. To better understand this time period, the Paleoclimate Modelling Intercomparison Project (PMIP) has provided a framework for an internationally coordinated endeavour in simulating the last deglaciation whilst encompassing a broad range of models. Here, we present a multi-model intercomparison of 17 transient simulations of the early part of the last deglaciation (∼20–15 ka BP) from nine different climate models spanning a range of model complexities and uncertain boundary conditions and forcings. The numerous simulations available provide the opportunity to better understand the chain of events and mechanisms of climate changes between 20 and 15 ka BP and our collective ability to simulate them. We conclude that the amount of freshwater forcing and whether it follows the ice sheet reconstruction or induces an inferred Atlantic meridional overturning circulation (AMOC) history, heavily impacts the deglacial climate evolution for each simulation rather than differences in the model physics. The course of the deglaciation is consistent between simulations except when the freshwater forcing is above 0.1 Sv – at least 70 % of the simulations agree that there is warming by 15 ka BP in most places excluding the location of meltwater input. For simulations with freshwater forcings that exceed 0.1 Sv from 18 ka BP, warming is delayed in the North Atlantic and surface air temperature correlations with AMOC strength are much higher. However, we find that the state of the AMOC coming out of the Last Glacial Maximum (LGM) also plays a key role in the AMOC sensitivity to model forcings. In addition, we show that the response of each model to the chosen meltwater scenario depends largely on the sensitivity of the model to the freshwater forcing and other aspects of the experimental design (e.g. CO2 forcing or ice sheet reconstruction). The results provide insight into the ability of our models to simulate the first part of the deglaciation and how choices between uncertain boundary conditions and forcings, with a focus on freshwater fluxes, can impact model outputs. We can use these findings as helpful insight in the design of future simulations of this time period.

Snoll, Brooke (ORCID:0000000230914394)↗
DOE OSTI2020

Performance and Accuracy Implications of Parallel Split Physics-Dynamics Coupling in the Energy Exascale Earth System Atmosphere Model

Simultaneous calculation of atmospheric processes is faster than calculating processes one at a time. This type of parallelism is beneficial or perhaps even necessary to provide good performance on modern supercomputers, which achieve faster performance through increased processor count rather than improved clock speed. The scalability of the Energy Exascale Earth System Model (E3SM) Atmosphere Model (EAM) is limited by the fluid dynamics which scales up to the number of mesh cells in the global mesh. In contrast, the suite of physics parameterizations in EAM is scalable up to the total number of physics columns, which is an order of magnitude greater than the number of mesh cells. A proposed solution to unlocking the greater potential performance from the physics suite is to solve the physics and dynamics in parallel. This work represents a first attempt at parallel splitting of the grid-scale fluid dynamics model and the subgrid-scale physics parameterizations in a global atmosphere model. We will demonstrate that switching to parallel physics-dynamics coupling extends the scalability of the EAM to up to 3 times the previous peak scalability limit and is up to 20% faster than the sequentially split coupling at the highest core counts and the same time step. Decadal simulations of both coupling approaches show very little impact to the model climate. This improved performance does not come without drawbacks, however. Parallel splitting requires a shorter time step and other modifications which largely offset performance gains. A mass fixer is required for conservation. Techniques for mitigating these issues are also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evaluation of the Near-Surface Variables in the HRRR Weather Model Using Observations from the ARM SGP Site

Abstract The performance of version 4 of the NOAA High-Resolution Rapid Refresh (HRRR) numerical weather prediction model for near-surface variables, including wind, humidity, temperature, surface latent and sensible fluxes, and longwave and shortwave radiative fluxes, is examined over the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) region. The study evaluated the model’s bias and bias-corrected mean absolute error relative to the observations on different time scales. Forecasts of near-surface geophysical variables at five SGP sites (HRRR at 3-km scale) were found to agree well with observations, but some consistent observation–forecast differences also occurred. Sensible and latent heat fluxes are the most challenging variables to be reproduced. The diurnal cycle is the main temporal scale affecting observation–forecast differences of the near-surface variables, and almost all of the variables showed different biases throughout the diurnal cycle. Results show that the overestimation of downward shortwave and the underestimation of downward longwave radiative flux are the two major biases found in this study. The timing and magnitude of downward longwave flux, wind speed, and sensible and latent heat fluxes are also different with contributions from model representations, data assimilation limitations, and differences in scales between HRRR and SGP sites. The positive bias in downward shortwave and negative bias in longwave radiation suggests that the model is underestimating cloud fraction in the study domain. The study concludes by showing a brief comparison with version 3 of the HRRR and shows that version 4 has better performance in almost all near-surface variables. Significance Statement A correct representation of the near-surface variables is important for numerical weather prediction models. This study investigates the capability of the latest NOAA High-Resolution Rapid Refresh (HRRRv4) model in simulating the near-surface variables by comparing against the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) in situ observations. Among others, we find that the surface heat fluxes, such as sensible and latent heat fluxes, are the most difficult variables to be reproduced. This study also shows that the diurnal cycle has the dominant impact on the model’s performance, which means the majority of the outputted near-surface variables have the strong diurnal cycle in their bias errors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Perspective on the Impact of Group Delay Dispersion in Future Terahertz Wireless Systems

This article discusses the challenges and opportunities of managing group delay dispersion (GDD), and its relation to the performance standards of future sixth-generation (6G) wireless communication systems utilizing terahertz frequency waves. The unique susceptibilities of 6G systems to GDD are described, along with a quantitative description of the sources of GDD, including multipath, rough surface scattering, intelligent reflecting surfaces, and propagation through the atmosphere. An experimental case-study is presented that confirms previous models quantifying the impact of atmospheric GDD. Several GDD manipulation strategies are presented, illustrating their hindered effectiveness in the 6G context. Conversely, some benefits of leveraging GDD to enhance 6G systems, such as improved security and simplified hardware, are also discussed. Finally, a perspective on using photonic GDD control devices is provided, revealing quantitative benefits that may unburden existing equalization schemes. Here, the article argues that GDD will uniquely and significantly impact some 6G systems, but that its careful consideration along with new mitigation strategies, including photonic devices, will help optimize system performance. The conclusion provides a perspective to guide future research in this area.

Strecker, Karl↗
DOE OSTI2022

Tracking the environmental impacts of ecological engineering on coastal wetlands with numerical modeling and remote sensing

Coastal wetlands are the most valuable ecosystems on the earth but facing severe degradation and losses owing to climate change and anthropogenic activities. Many ecological engineering projects (EEP) have been conducted to mitigate the degradation of coastal wetlands. However, the geomorphological impacts of EEP on coastal wetlands have not been well documented. In this study, a method employed a process-based hydrodynamic model and remote sensing (RS) was developed to evaluate the impacts of EEP on the geomorphological change of a prototype Ramsar site. Results demonstrated that RS has great potential in improving the quality of bathymetry data for the numerical model with a decrease of RMSE from 0.52 m to 0.3 m. It also showed good capacity in trend detection of geomorphological change spatially. Results showed the Chongming Dongtan wetland experienced erosion with an annual rate of -0.035 m/yr from 2013 to 2016 after the implementation of EEP. The deposition rate changed significantly in the area within 200 m of the EEP. It is found that the EEP modified the composition of vegetation, sediment transportation, as well as substrate stability, affecting the geomorphological change of coastal wetlands. It is suggested that the EEP with moderate anthropogenic disturbance is a direct and effective way to recover the coastal habitats for waterbirds. However, the modification of the coastal wetland ecosystem by EEP will lead to the potential vulnerability to global climate change. Therefore, how to mitigate the advantages and disadvantages of the EEP is needed to be further studied to find a more sustainable way for coastal management.

geomorphological change, coastal wetlands, Ecologi↗
DOE OSTI2020

A conceptual model for how to design for building envelope characteristics. Impact of thermal comfort intervals and thermal mass on commercial buildings in U.S. climates

The paper presents a simplified conceptual model for energy demand calculations based on building envelope characteristics, thermal mass and local climate. It is based on a network model and lumped analysis of the dynamic process. Characteristic parameters for the buildings are suggested; Driving temperature (DT), Driving temperature difference, (DTD), External Load Temperature (ELT), and Thermal Load Resistance (TLR). The Building Envelope Performance ( BEP 0 ), based on a controlled constant indoor temperature is introduced. Solution techniques using stable explicit forward differences based on analytical solutions are derived. The conceptual model has been used for mapping the Driving temperature difference and introduced two performance factors and . The first factor represents the effect of thermal comfort interval and thermal mass on the energy demand. The latter represents the ratio between cooling and heating energy demand. These three parameters and factors have been visualized on U.S. maps and enable a possibility to communicate the demand of energy, and cooling and the coupling to building characteristics, in a concise way.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Effect of macroscopic surface defects on dynamic damage: An experimental and numerical study

This study examines the impact of macroscopic surface defects on the dynamic ductile damage behavior of polycrystalline metals using plate-impact experiments. Defects of various shapes (flat, round, and point) were manufactured on the free or impact surfaces of annealed copper specimens. The experiments were diagnosed with photon Doppler velocimetry measurements and soft recovery techniques. The experimental results revealed that defect shape and location significantly affect velocity–time profiles and void distribution. In order to understand the dynamics of shock propagation and corresponding ductile damage evolution within the target specimens, the impact experiments were modeled with a calibrated Tepla model, a dynamic ductile damage model for polycrystalline metals [Nguyen et al., Int. J. Solids Struct. 329, 113833 (2026)]. Overall, our resulting simulated velocity showed good agreement with measured velocity, and our simulated porosity distributions qualitatively matched experimental data. Based on our simulation results, defects on the free surface were found to distort rarefaction waves and therefore the corresponding spall planes. On the other hand, defects on the impact surface generated a delayed shock when the flyer plate and the defected target area were in contact, leading to distortion of the spall plane. The distortion of the spall plane resulted in a non-uniform distribution of voids within the defected specimens. Lastly, we examined the role of local work hardening due to the defect manufacturing process on the velocity and porosity distribution, using Tepla simulations with a simplified representation of local hardening. This investigation highlights the importance of defect geometry, location, and local hardening associated with defect generation in dynamic ductile damage processes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Comparing computational times for simulations when using PBPK model template and stand-alone implementations of PBPK models

Introduction We previously developed a PBPK model template that consists of a single model “superstructure” with equations and logic found in many physiologically based pharmacokinetic (PBPK) models. Using the template, one can implement PBPK models with different combinations of structures and features. Methods To identify factors that influence computational time required for PBPK model simulations, we conducted timing experiments using various implementations of PBPK models for dichloromethane and chloroform, including template and stand-alone implementations, and simulating four different exposure scenarios. For each experiment, we measured the required computational time and evaluated the impacts of including various model features (e.g., number of output variables calculated) and incorporating various design choices (e.g., different methods for estimating blood concentrations). Results We observed that model implementations that treat body weight and dependent quantities as constant (fixed) parameters can result in a 30% time savings compared with options that treat body weight and dependent quantities as time-varying. We also observed that decreasing the number of state variables by 36% in our PBPK model template led to a decrease of 20–35% in computational time. Other factors, such as the number of output variables, the method for implementing conditional statements, and the method for estimating blood concentrations, did not have large impacts on simulation time. In general, simulations with PBPK model template implementations of models required more time than simulations with stand-alone implementations, but the flexibility and (human) time savings in preparing and reviewing a model implemented using the PBPK model template may justify the increases in computational time requirements. Conclusion Our findings concerning how PBPK model design and implementation decisions impact computational speed can benefit anyone seeking to develop, improve, or apply a PBPK model, with or without the PBPK model template.

Bernstein, Amanda S.↗
DOE OSTI2023

Supporting Nuclear Energy Research with MLOps

Artificial intelligence (AI) and machine learning (ML) have been shown to be increasingly helpful tools in a growing number of use-cases relevant to scientific research, despite significant software-related obstacles. There exist large technical costs to setting up, using, and maintaining AI/ML models in production. This often prevents researchers from utilizing these models in their work. The growing field of machine learning operations (MLOps) aims to automate much of the AI/ML life cycle while increasing access to these models. This paper presents our initial work in creating a nuclear energy MLOps platform for use by researchers at Idaho National Laboratory (INL) and aims to reduce the barriers of using AI/ML in scientific research. We hope to promote the integration of the latest AI/ML technologies into researchers' workflows and create more opportunity for scientific innovation. In this paper we discuss how our MLOps efforts aim to increase usage and the impact of AI/ML models created by engineers. We also present several use-cases that are currently integrated. Finally, we evaluate the maturity of our project as well as our plans for future functionality.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

SUHMO: an adaptive mesh refinement SUbglacial Hydrology MOdel v1.0

Abstract. Water flowing under ice sheets and glaciers can have a strong influence on ice dynamics, particularly through pressure changes, suggesting that a comprehensive ice sheet model should include the effect of basal hydrology. Modeling subglacial hydrology remains a challenge, however, mainly due to the range of spatial and temporal scales involved – from subglacial channels to vast subglacial lakes. Additionally, networks of subglacial drainage channels dynamically evolve over time. To address some of these challenges, we have developed an adaptive mesh refinement (AMR) model based on the Chombo software framework. We extend the model proposed by Sommers et al. (2018) with a small but significant change to accommodate the transition from unresolved to resolved flow features. We handle the strong nonlinearities present in the equations by resorting to an efficient nonlinear full approximation scheme multigrid (FAS-MG) algorithm. We outline the details of the algorithm and present convergence analysis results demonstrating its good performance. Additionally, we present results validating our approach, using test cases from the Subglacial Hydrology Model Intercomparison Project (SHMIP) (de Fleurian et al., 2018). We finish by presenting a more complex, 100 km-by-100 km synthetic test case with peaks and valleys that we use to investigate the effective pressure distribution as the number of AMR levels increases. These preliminary results suggest that a minimum spatial resolution is needed to properly capture channel features, but additional work is required to precisely quantify this and its impact on accurately modeling the coupled ice sheet–hydrology system. The efficiency of our approach, relying on localized refinement, is also demonstrated. Future work will include coupling the SUbglacial Hydrology MOdel (SUHMO) with the BISICLES AMR ice sheet model (Cornford et al., 2013), both built on the same numerical framework.

58 GEOSCIENCES↗
DOE OSTI2022

Synthetic Data Generation Using Machine Learning

Robust machine learning techniques for image analysis require a substantial amount of data to yield confident results. In the nuclear domain, data scarcity is a substantial challenge because there are so few facilities worldwide. This research focuses on being able alleviate the data scarcity problem by generating synthetic data to bridge the gap between large and small datasets. This work achieves that goal using a Generative Adversarial Network (GAN) architectural approach, by training a model on real-world data and expands that small dataset through synthetic data amendments. Model performance is impacted by the size of the real-world dataset and the number of training epochs utilized. This means that 1) It is important to develop your GAN to be optimized with the specific data type, and 2) approaches taken when training the GAN should be specialized to encompass important aspects of the dataset that it is generating. By taking a step to improve dataset sizes in this way, the gap between models trained by parties with significant amount of data and those without access to large data, closes, allowing for robust analyses of satellite imagery for nuclear domain applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Power System Operational Impacts of Electric Vehicle Dynamic Wireless Charging

The electrification of the transportation sector poses an opportunity for reducing greenhouse gas (GHG) emissions from passenger vehicles. Electric vehicle (EV) charging through dynamic wireless power transfer (DWPT), known as roadway electrification, could shift EV demand profiles to better coincide with renewable electricity generation. However, this would be a very large new load and few studies evaluate the regional impacts of DWPT charging in a power transmission system. This paper defines methods that address dataset generation for passenger vehicle trips and models to evaluate regional impacts for this emerging technology. Household vehicle miles traveled (VMT) data form localized EV demand profiles through discrete-event simulation. This data serves as exogenous inputs for a Production Cost Model (PCM) of a synthetic transmission system based on the Electric Reliability Council of Texas's (ERCOT) network. EV charging methods are compared for both a 2018 baseline generation mixture and a high-renewable generation case incorporating 20 GW of installed solar photovoltaic (PV) capacity. The PCM employs unit commitment and economic dispatch (UC&ED) models to compare financial, environmental, and grid reliability impacts from EV charging across passenger EV adoption levels. In-transit charging could reduce grid operational costs by as much as 1.49%, with up to $13.7B saved in annual vehicle operational costs for consumers compared to gas-powered vehicles. Health impacts analysis from power plant and vehicle tailpipe emissions from this study show net health benefits increase by 40% for in-transit charging coupled with high renewable generation. Renewable resources provide an avenue for cost-effective in-transit charging with reduced emissions. The combination of dataset generation and open-source power system modeling establish a foundation for the holistic evaluation of regional DWPT impacts.

dynamic wireless power transfer↗
DOE OSTI2023

Emissions and Energy Impacts of the Inflation Reduction Act

If goals set under the Paris Agreement are met, the world may hold warming well below 2 degrees C (1); however, parties are not on track to deliver these commitments (2), increasing focus on policy implementation to close the gap between ambition and action. Recently, the US government passed its most prominent piece of climate legislation to date - the Inflation Reduction Act of 2022 (IRA) - designed to invest in a wide range of programs that, among other provisions, incentivize clean energy and carbon management, encourage electrification and efficiency measures, reduce methane emissions, promote domestic supply chains, and address environmental justice concerns (3). IRA's scope and complexity make modeling important to understand impacts on emissions and energy systems. We leverage results from nine independent, state-of-the-art models to examine potential implications of key IRA provisions, showing economy-wide emissions reductions between 43 and 48% below 2005 levels by 2035.

electricity↗
DOE OSTI2021

Dynamic Linking of Upstream Energy and Freight Demands for Bio and Fossil Energy Pathways in the Global Change Analysis Model

Comprehensive study of the environmental impacts associated with demand for an energy resource or carrier in any one sector requires a full consideration of the direct and indirect impacts on the rest of the regional and global energy system. It is important to consider the energy that is consumed in producing primary energy resources and energy consumed in transporting these energy resources internationally and domestically in order to produce a complete picture of these impacts. Biofuels are especially complex since they have feedbacks not just to the energy system but also to regional and global crop markets. Different modeling and analysis strategies have been applied to consider these impacts, which occur “upstream” of final energy consumption, with some success. Traditional life cycle models allow for rich technological detail in linking upstream impacts of energy demand with downstream final consumption; however, they cannot by themselves represent dynamic economic feedbacks across multiple sectors and regions over time. Computable general equilibrium (CGE) modeling does account for economic feedback among all sectors and regions; however, CGE modeling does not lend itself to physical representation of technological detail. Partial equilibrium (PE) modeling is a heterogeneous category describing economic models that focus on a subset of the economy, and they differ in coverage and ease in incorporating more sectors and economic links. In this study, we present a strategy for dynamically including the direct and indirect impacts of energy demand with physical and technology detail by explicitly adding these upstream energy and transportation links to GCAM, a PE model of global energy, land use, and emissions. We incorporate the following inter-sectoral linkages: energy inputs to crop production, energy inputs to fossil resource production, and freight transport requirements of energy and agricultural commodities. We assess the implications of explicitly including these links by measuring the global impacts of increased corn ethanol demand in the United States with and without these links included. Although the net global impact of the upstream links on energy and emissions are relatively modest in the scenarios we studied, the inclusion of these links illustrates interesting trade-offs in energy and transportation demand among fossil fuel and agriculture sectors. These sectoral interactions suggest that this level of modeling detail could be important in evaluating future analytical questions.

Sampedro Martinez de Estivariz, Jon↗