Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JOINT”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Large-deviation joint statistics of the finite-time Lyapunov spectrum in isotropic turbulence

One of the hallmarks of turbulent flows is the chaotic behavior of fluid particle paths with exponentially growing separation among them while their distance does not exceed the viscous range. The maximal (positive) Lyapunov exponent represents the average strength of the exponential growth rate, while fluctuations in the rate of growth are characterized by the finite-time Lyapunov exponents (FTLEs). In the last decade or so, the notion of Lagrangian coherent structures (which are often computed using FTLEs) has gained attention as a tool for visualizing coherent trajectory patterns in a flow and distinguishing regions of the flow with different mixing properties. A quantitative statistical characterization of FTLEs can be accomplished using the statistical theory of large deviations, based on the so-called Cramér function. To obtain the Cramér function from data, we use both the method based on measuring moments and measuring histograms and introduce a finite-size correction to the histogram-based method. We generalize the existing univariate formalism to the joint distributions of the two FTLEs needed to fully specify the Lyapunov spectrum in 3D flows. The joint Cramér function of turbulence is measured from two direct numerical simulation datasets of isotropic turbulence. Results are compared with joint statistics of FTLEs computed using only the symmetric part of the velocity gradient tensor, as well as with joint statistics of instantaneous strain-rate eigenvalues. When using only the strain contribution of the velocity gradient, the maximal FTLE nearly doubles in magnitude, highlighting the role of rotation in de-correlating the fluid deformations along particle paths. We also extend the large-deviation theory to study the statistics of the ratio of FTLEs. The most likely ratio of the FTLEs λ{sub 1} : λ{sub 2} : λ{sub 3} is shown to be about 4:1:−5, compared to about 8:3:−11 when using only the strain-rate tensor for calculating fluid volume deformations. The results serve to characterize the fundamental statistical and geometric structure of turbulence at small scales including cumulative, time integrated effects. These are important for deformable particles such as droplets and polymers advected by turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Prospects for joint observations of gravitational waves and gamma rays from merging neutron star binaries

The detection of the events GW150914 and GW151226, both consistent with the merger of a binary black hole system (BBH), opened the era of gravitational wave (GW) astronomy. Besides BBHs, the most promising GW sources are the coalescences of binary systems formed by two neutron stars or a neutron star and a black hole. These mergers are thought to be connected with short Gamma Ray Bursts (GRBs), therefore combined observations of GW and electromagnetic (EM) signals could definitively probe this association. We present a detailed study on the expectations for joint GW and high-energy EM observations of coalescences of binary systems of neutron stars with Advanced Virgo and LIGO and with the Fermi gamma-ray telescope. To this scope, we designed a dedicated Montecarlo simulation pipeline for the multimessenger emission and detection by GW and gamma-ray instruments, considering the evolution of the GW detector sensitivities. We show that the expected rate of joint detection is low during the Advanced Virgo and Advanced LIGO 2016–2017 run; however, as the interferometers approach their final design sensitivities, the rate will increase by ∼ a factor of ten. Future joint observations will help to constrain the association between short GRBs and binary systems and to solve the puzzle of the progenitors of GWs. Comparison of the joint detection rate with the ones predicted in this paper will help to constrain the geometry of the GRB jet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Development of thermal joints for conduction-cooling applications

Abstract In support of the development of conduction-cooled superconducting radiofrequency (SRF) niobium cavities for use in continuous-wave linear accelerators, an experimental study of thermal contact resistance was performed on bolted joints using high-purity niobium, aluminium, and copper along with Apiezon N grease and indium foil as interfacial materials. The geometry of each joint investigated aims at replicating the design adopted in conduction-cooled SRF cavities currently under development at Jefferson Lab. The materials’ thermal conductivity and the joints’ thermal resistance were measured in the temperature range of 3.5 - 10 K. The results showed that the low thermal contact resistance of 2-4 K·cm 2 /W at 4.3 K achieved with Al-In-Nb and Al-Apiezon N-Cu joints should allow the operation of conduction-cooled SRF cavities at a targeted accelerating gradient of ∼10 MV/m.

Lewis, J. [Old Dominion University]↗
DOE OSTI2022

Joint Planning of Electricity Transmission and Hydrogen Transportation Networks

The abundance and uneven distribution of renewable energy might cause congestion and curtailment in electric power systems. Transmission expansions can potentially alleviate transmission congestion and reduce renewable energy curtailment. On the other hand, with the substantial cost reduction of electrolyzer technology and the continuing rise of hydrogen demand, converting surplus renewable energy to hydrogen provides synergistic benefits to power and hydrogen systems, but consequently requires joint planning of the two systems. To this end, we propose a joint planning approach for power transmission and hydrogen transportation networks to coordinately optimize the investment and operation of power and hydrogen infrastructure. In our proposed model, detailed truck routing, pipeline, and hydrogen storage are formulated to quantify the flexibility of hydrogen transportation system. A robust joint planning approach is also proposed to address various uncertainties from renewable energy, electric load, and hydrogen demand. Furthermore, our numerical simulations show that the proposed joint planning model can save the total system cost, reduce renewable energy curtailment, and increase the utilization level of transmission lines.

08 HYDROGEN↗
DOE OSTI2019

Joint Chance Constraints in AC Optimal Power Flow: Improving Bounds through Learning

This paper considers distribution systems with a high penetration of distributed, renewable generation and addresses the problem of incorporating the associated uncertainty into the optimal operation of these networks. Joint chance constraints, which satisfy multiple constraints simultaneously with a prescribed probability, are one way to incorporate uncertainty across sets of constraints, leading to a chance-constrained optimal power flow problem. Departing from the computationally-heavy scenario-based approaches or approximations that transform the joint constraint into conservative deterministic constraints, this paper develops a scalable, data-driven approach which learns operational trends in a power network, eliminates zero-probability events (e.g., inactive constraints), and accurately and efficiently approximates bounds on the joint chance constraint iteratively. In particular, the proposed framework improves upon the classic methods based on the union bound (or Boole's inequality) by generating a much less conservative set of single chance constraints that also guarantees the satisfaction of the original joint constraint. The proposed framework is evaluated numerically using the IEEE 37-node test feeder, focusing on the problem of voltage regulation in distribution grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

A Model for Joint Probabilistic Forecast of Solar Photovoltaic Power and Outdoor Temperature

In this paper, a stochastic model is proposed for a joint statistical description of solar photovoltaic (PV) power and outdoor temperature. The underlying correlation emerges from solar irradiance that is responsible in part for both the variability in solar PV power and temperature. The proposed model can be used to capture the uncertainty in solar PV power and its correlation with the electric power consumption of thermostatically controlled loads. First, a model for solar PV power that explicitly incorporates the stochasticity due to clouds via a regime-switching process between the three classes of sunny , overcast and partly cloudy is proposed. Then, the relationship between temperature and solar power is postulated using a second-order Volterra model. This joint modeling is leveraged to develop a joint probabilistic forecasting method for solar PV power and temperature. Real-world datasets that include solar PV power and temperature measurements in California are analyzed and the effectiveness of the joint model in providing probabilistic forecasts is verified. The proposed forecasting methodology outperforms several reference methods thus portraying that the underlying correlation between temperature and solar PV power is well defined and only requires a simple lower-complexity sampling space.

14 SOLAR ENERGY↗
DOE OSTI2022

The use of subchondral bone topography to approximate the center of rotation of the elbow joint in dogs

Abstract Objectives The aim of this study was to compare the approximate center of rotation in normal and diseased elbows in dogs. Study design Cross‐sectional study Sample population Computer tomography scans of nine dogs with unilateral fragmented medial coronoid process (FMCP). Methods A board certified radiologist confirmed that each dog had unilateral FMCP, and a normal contralateral elbow. Digital 3D models of all elbow joints were uploaded into a surgical planning software package. Four axes approximating the center of rotation (COR) of elbow joints were generated using five geometric shapes based on subchondral topography of the humeral condyle radius and ulna. Images showing the locations where axes exited the medial and lateral cortex of the humeral condyle were captured and imported into a second software package, for measurement of distances between exit points and the origin of a system of axes. Results In normal joints 20/27 (74%) axes exited the medial cortex, and 25/27 (93%) axes exited the lateral cortex cranial and distal to the medial and lateral epicondyles, respectively. In diseased joints 22/27 (81%) axes exited medial cortex and 19/27 (70%) axes exited the lateral cortex, caudal and distal to the medial and lateral epicondyles, respectively. Conclusion Based on CT‐ derived geometry, the COR of elbow affected with FMCP was generally more caudal than normal. Clinical significance External landmarks approximating the location of the elbow COR are provided, and while not validated, may assist in planning, creation, and assessment of procedures for FMCP.

Berger, Chen↗
DOE OSTI2018

Quadri-joint inversion: Method and application to the Big Sky 9C 3D data set in northern Montana

The Kevin Dome [Formula: see text] storage project, located in northern Montana, attempted to characterize the Duperow Formation as a potential long-term storage zone for injected [Formula: see text]. A multicomponent (9C) seismic survey was acquired for the Big Sky Carbon Sequestration Partnership over a portion of the Kevin Dome using P- and S-wave sources. Prestack migrated PP, PS, SH, and SV data sets were generated. We then applied several stratigraphic inversion workflows using one or several kinds of seismic wavefield at the same time resulting in joint inversions of each data set. The aim of our study is to demonstrate the benefits of doing quadri-joint inversion of PP-, PS-, SH-, and SV-wavefields for the recovery of the elastic earth parameters, especially the S-wave impedance and density. These are crucial parameters because they can help determine lithology and porefill in the reservoir characterization workflow. Because the inversion workflow always uses the original seismic data recorded in its own time domain, it is necessary to compute registration laws between PP-PS-, PP-SH-, and PP-SV-wavefields using a time shift computation procedure (warping) based on inverted S-wave impedances from inversion of a single wavefield. This generated a significant improvement over methods that rely on attempting to match trace waveforms that may have a different phase, frequency content, and polarity. Finally, we wanted to investigate the reliability of the quadri-joint inversion results in the Bakken/Banff Formations, which have less lateral geologic variation than the underlying Duperow target. This interval shares many of the geophysical characterization challenges common to shale reservoirs in other North American basins. We computed geomechanical parameters, such as Poisson’s ratio and Young’s modulus, which are a proxy for brittleness. Comparison of these results with independent laboratory measurements in the Bakken interval demonstrates the superiority of the quadri-joint inversion method to the traditional inversion using P-wave data only.

Geochemistry & Geophysics↗
DOE OSTI2013

Joint Analysis of Surface Deformation Data for Reservoir Management (CRADA Final Report)

This was a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Chevron North America Exploration and Production Company, a division of Chevron U.S.A. Inc. (Chevron), to develop an operational tool to carry out joint inversion of global positioning system (GPS), Interferometric Synthetic Aperture Radar (InSAR) and tilt, temperature, and electrical resistivity data to monitor the temporal and spatial evolution of fluid volume and temperature in a thermal diatomite oilfield. The original goal of the research and development carried out under the CRADA was to conduct joint inversion of GPS, InSAR and tilt data sets from Chevron’s McKittrick field. Temperature and electrical resistivity data were added to the joint inversion scheme at later stages of the project. The overall objective of the joint inversion was to provide improved mapping over time of changes in the spatial distributions of reservoir fluid volume and temperature within the reservoir in response to injection and production. This CRADA was originally designated as a three (3) year project, and was to be carried out in four phases. Each task was numbered according to the project phase it fell under. The four phases are identified below. Phases 3 and 4 were developed based on the results obtained in Phase 2, and the detailed work plan and refined budget were added by amendment.

Foxall, William↗
DOE OSTI2017

Computing a Non-trivial Lower Bound on the Joint Entropy between Two Images

In this report, a non-trivial lower bound on the joint entropy of two non-identical images is developed, which is greater than the individual entropies of the images. The lower bound is the least joint entropy possible among all pairs of images that have the same histograms as those of the given images. New algorithms are presented to compute the joint entropy lower bound with a computation time proportional to S log S where S is the number of histogram bins of the images. This is faster than the traditional methods of computing the exact joint entropy with a computation time that is quadratic in S .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Corrosion Control in Carbon Fiber Reinforced Plastic (CFRP) Composite-Aluminum Closure Panel Hem Joints

This project targeted the implementation of CFRP/Al closure panels as a drop-in replacement for all-aluminum closures used in automobile production today. The project focused on the development of corrosion fundamentals involving CFRP materials and the use of those materials in mixed-material joints. Numerous studies were completed to benchmark current materials and test mixed metal joints. Good correlation was found between outdoor exposures and laboratory cyclic corrosion testing. The insights from the corrosion testing fundamentals and benchmarking phases informed the development of low-cure adhesive and electrocoats and conductive primers for the CFRP. Adhesive and electrocoat materials were successfully developed meeting low-cure targets for bake (150C/10min). These materials were scaled-up and used to produce 5 prototype liftgates to further evaluate CFRP/Al joints in a full-scale part. The lower cure materials enabled mitigation of CTE mismatch and reduced the required bake temperature, resulting in reduced part distortion. Corrosion testing of the liftgates revealed unexpected increases in corrosion relative to the standard bake materials. While rigorous evaluation of the parts was not possible due to numerous complexities associated with prototype testing, the results provided significant insights into the material behavior of CFRP/Al joints. Key learnings from the project included a better understanding of the system complexity (interplay of substrate/primer conductivity, production of full-scale prototypes), advances in coating/adhesive design for low-cure applications involving CFRP, and advancement of corrosion test protocols for CFRP.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonlinear Characterization of a Joint Exhibiting a Reduction in Damping at High Energy

Even with the advent of additive manufacturing, the vast majority of complex structures are comprised of individual components held together with bolted joints. However, bolted joints present a challenge for mechanical design as they are a source of nonlinearity and increase the uncertainty in the overall behavior of the system in a dynamic environment. While many advances have been made in the ability to accurately model and test bolted joints, it is still an open area of research. Modes of vibration that exercise bolted joints typically exhibit nonlinear behavior where, with increased excitation level, the natural frequency decreases (i.e. softens) and the damping increases. However, the system under study for this work has an axial mode which does not follow this trend; it does soften as expected, but, after an initial increase, the apparent damping decreases with excitation amplitude. At the highest excitation level, the frequency of the mode decreases to that of a nearby bending mode and the response is amplified nearly 500% above that at lower levels. It is unclear whether the decrease in damping is due to the coupling of the two modes or if it is a characteristic of the axial mode. Therefore, the objective of this project is to investigate the coupling between the axial and bending modes and the dynamics leading to the decrease in damping.

42 ENGINEERING↗
DOE OSTI2022

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2021

Characterization and High Temperature Electrical Properties of Brazed Joints of La0.8Sr0.2CrO3 with Nickel and Nickel Alloys

This study investigates the high-temperature electrical properties and microstructural characteristics of brazed joints combining La₀.₈Sr₀.₂CrO₃ (LSC) with nickel and nickel-based alloys. Using nickel-based filler metals (NB LM and NB 130), the joints were tested in controlled atmospheres (95% N₂ and 5% H₂). Electrical resistivity was analyzed across a range of temperatures (25–727°C), and contact resistance was evaluated using a composite model. Results indicate that joint resistance decreases with increasing temperature, with optimal brazing achieved at 1160°C. Diffusion processes involving Si, Ni, O, and Sr were identified, with silicon oxide formation contributing to joint bonding.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2019

Determining the Tolerances in Fitting for Joints with Interference

A formula for the theoretical tolerance in fitting is derived by analysis of the components in the model of parametric failure for joints with interference. The formula takes account of the life of the joint and the specified probability of fault-free operation. It is found that selection of the fit without a margin of precision does not ensure the necessary reliability of the joint, while correct selection of the margin of precision extends the fault-free operation of the joint.

42 ENGINEERING↗
DOE OSTI2018

Capacity Optimization of a Centralized Charging Station in Joint Operation with a Wind Farm

In the context of large-scale wind power integration and rapid development of electric vehicles (EVs), a joint operation pattern was proposed to use a centralized charging station (CCS) to address high uncertainties incurred by wind power integration. Additionally, this would directly remove the significant indeterminacy of wind power. Because the CCS is adjacent to a wind power gathering station, it could work jointly with wind farms to operate in power system as an independent enterprise. By combining the actual operational characteristics of the wind farm and the CCS, a multidimensional operating index evaluation system was created for the joint system. Based on the wind farm’s known capacity, a multi-target capacity planning model was established for CCS to maximize the probability of realizing diverse indices of the system based on dependent-chance goal programming. In this model, planning and dispatching are combined to improve the feasibility of results in the operational stage. This approach takes the effects of randomness into account for wind power and battery swapping. The model was solved through combining Monte Carlo simulation and a genetic algorithm (GA) based on segmented encoding. Finally, as shown in the simulation results, the proposed model could comprehensively include factors such as initial investment, wind power price, battery life, etc., to optimize CCS capacity and to ultimately improve the operating indices of the joint system.

17 WIND ENERGY↗
DOE OSTI2024

Convex Relaxations of Maximal Load Delivery for Multi-Contingency Analysis of Joint Electric Power and Natural Gas Transmission Networks

Recent increases in gas-fired power generation have engendered increased interdependencies between natural gas and power transmission systems. These interdependencies have amplified existing vulnerabilities in gas and power grids, where disruptions can require the curtailment of load in one or both systems. Although typically operated independently, coordination of these systems during severe disruptions can allow for targeted delivery to lifeline services, including gas delivery for residential heating and power delivery for critical facilities. To address the challenge of estimating maximum joint network capacities under such disruptions, we consider the task of determining feasible steady-state operating points for severely damaged systems while ensuring the maximal delivery of gas and power loads simultaneously, represented mathematically as the nonconvex joint Maximal Load Delivery (MLD) problem. To increase its tractability, we present a mixed-integer convex relaxation of the MLD problem. Then, to demonstrate the relaxation’s effectiveness in determining bounds on network capacities, exact and relaxed MLD formulations are compared across various multi-contingency scenarios on nine joint networks ranging in size from 25 to 1191 nodes. The relaxation-based methodology is observed to accurately and efficiently estimate the impacts of severe joint network disruptions, often converging to the relaxed MLD problem’s globally optimal solution within ten seconds.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗