• Energy Research
• 2025-2025close
• other engineering and technologies close
• US close

With California's ambitious goal to achieve decarbonization of the electrical grid by the year 2045, significant challenges arise in power system investment planning. Existing modeling methods and software focus on computational efficiency, which is currently achieved by simplifying the associated unit commitment formulation. This may lead to unjustifiable inaccuracies in the cost and constraints of gas-fired generation operations, and may affect both the timing and the extent of investment in new resources, such as renewable energy and energy storage. To address this issue, this paper develops a more detailed and rigorous mixed-integer model, and more importantly, a solution methodology utilizing surrogate level-based Lagrangian relaxation to overcome the combinatorial complexity that results from the enhanced level of model detail. This allows us to optimize a model with approximately 12 million binary and 100 million total variables in under 48 hours. The investment plan is compared with those produced by E3's RESOLVE software, which is currently employed by the California Energy Commission and California Public Utilities Commission. Our model produces an investment plan that differs substantially from that of the existing method and saves California over 12 billion dollars over the investment horizon.

We develop an explicit second order staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge domain boundary values that are significantly more complex than in the case of a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.