Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional Distribution”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Development and Evaluation of Chemistry‐Aerosol‐Climate Model CAM5‐Chem‐MAM7‐MOSAIC: Global Atmospheric Distribution and Radiative Effects of Nitrate Aerosol

Abstract An advanced aerosol treatment, with a focus on semivolatile nitrate formation, is introduced into the Community Atmosphere Model version 5 with interactive chemistry (CAM5‐chem) by coupling the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the 7‐mode Modal Aerosol Module (MAM7). An important feature of MOSAIC is dynamic partitioning of all condensable gases to the different fine and coarse mode aerosols, as governed by mode‐resolved thermodynamics and heterogeneous chemical reactions. Applied in the free‐running mode from 1995 to 2005 with prescribed historical climatological conditions, the model simulates global distributions of sulfate, nitrate, and ammonium in good agreement with observations and previous studies. Inclusion of nitrate resulted in ∼10% higher global average accumulation mode number concentrations, indicating enhanced growth of Aitken mode aerosols from nitrate formation. While the simulated accumulation mode nitrate burdens are high over the anthropogenic source regions, the sea‐salt and dust modes respectively constitute about 74% and 17% of the annual global average nitrate burden. Regional clear‐sky shortwave radiative cooling of up to −5 W m −2 due to nitrate is seen, with a much smaller global average cooling of −0.05 W m −2 . Significant enhancements in regional cloud condensation nuclei (at 0.1% supersaturation) and cloud droplet number concentrations are also attributed to nitrate, causing an additional global average shortwave cooling of −0.8 W m −2 . Taking into consideration of changes in both longwave and shortwave radiation under all‐sky conditions, the net change in the top of the atmosphere radiative fluxes induced by including nitrate aerosol is −0.7 W m −2 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Universal method for the optimization of HDC coating uniformity on non-planar, non-stationary substrates for inertial confinement fusion targets

The thickness uniformity of chemical vapor deposited (CVD) diamond coatings on non-planar, non-stationary substrates depends on both the intrinsic instantaneous coating thickness distribution (ICTD) of the coating conditions used and, if applicable, on the frequency of substrate reorientation. While important for many CVD diamond applications, the relative impact of the ICTD and substrate reorientation on the coating thickness uniformity has not been studied. In this work, we systematically investigate the effect of these factors for microwave-plasma chemical vapor deposition (MPCVD) of diamond (referred to as high density carbon (HDC) in the inertial confinement fusion (ICF) community) coatings on spherical, rolling substrates. This coating technique is used to fabricate capsules for ICF experiments, which require extreme coating uniformity with <0.3 % thickness variation (so-called Mode 1 or M1) to ensure symmetric compression of imploding targets. To extract the otherwise unobservable reorientation timescale (Δt), Monte Carlo simulations were performed using experimental ICTD data as input. This combined approach confirms scaling relationships between the substrate reorientation timescale as well as coating thickness and coating uniformity, as expected from a 3D random walk. Simulations confirm that M1 is Rayleigh-distributed and scales as (Δt) 1/2 , consistent with the randomization of two angles that determine orientation of a sphere. We also demonstrate that, under the conditions studied, Δt is the dominant factor in determining thickness uniformity while the intrinsic ICTD has minimal impact. Finally, experiments show that Δt can be affected by total batch size under constant agitation conditions due to space constraints that limit the capsule reorientation kinetics. In conclusion, this study highlights the utility of a combined experiment-simulation approach as a general methodology for understanding and improving coating uniformity on non-planar, non-stationary substrates.

Capsule↗
DOE OSTI2021

Experimental investigation of the effect of pore size distribution on nano-particle capture efficiency within ceramic particulate filters

The effect of the pore size distribution on size-resolved filtration efficiency was investigated for two ceramic particulate filters using particulate matter (PM) generated by a spark-ignition direct-injection engine fueled with gasoline. The cordierite filter tested had a porosity of 43%, a median pore diameter (d_50) of 12 µm, and a wide pore size distribution with a lognormal standard deviation (s’) of 0.3. The aluminum titanate filter had very similar porosity, d_50, and thickness, but significantly narrower pore size distribution (s’ = 0.1). The testing of two filters under identical experimental conditions enabled the impact of the pore size distribution on filtration performance to be evaluated. Filtration experiments were performed focusing on just the filter wall, starting from a clean filter until the transition to cake filtration (filtration efficiency > 99%). Time-resolved particle size distribution measurements were used to evaluate the progression of filtration performance and estimate trapped mass within the filter. The aluminum titanate filter, with a narrow pore size distribution, exhibited significantly better diffusion capture efficiency. The negative impact of higher flow velocity on diffusion capture efficiency was more pronounced for a narrower pore size distribution. Flow distributions measured using capillary-flow porometry were used to develop a cylindrical pore flow model to understand the impact of the differences in pore size distribution on observed trends in diffusion capture efficiency within a clean filter. The model helps demonstrate how large pores on the high end of a wider pore size distribution can carry a large fraction of the flow, with less efficient capture of small particles by diffusion. The experimental results and data demonstrate that the bubble point diameter and width of the pore size distribution significantly influence diffusion capture efficiencies for filters with very similar median pore diameter, porosity, and thickness.

aerosol filtration, Aftertreatment, ceramic exhaus↗
DOE OSTI2022

Preliminary Thermal-hydraulics and Selected Safety Analysis for Holos-Quad Reactor Design

This report presents thermal hydraulics and safety analyses of the 10 MWe Holos-Quad micro-reactor design developed by HolosGen LLC. These analyses were executed under steady-state and station black-out (SBO) transient conditions using the System Analysis Module (SAM), with a focus on ensuring peak fuel temperatures do not exceed the safety thresholds of 1,250°C for steady state normal operating conditions and 1,600°C for transient conditions. For steady-state temperature predictions, a coupled 1D fluid-to-3D solid heat conduction methodology was used. The steady-state simulations included a model of the central core region where power peaking occurs to predict the maximum steady-state fuel temperature, and a coarse-mesh full-core model to provide the initial condition for the consequent transient SBO decay heat removal simulation. The maximum fuel temperature predicted from the central core region simulations under the normal operating condition was 1,222°C, which is below the design limit 1,250°C. To simulate a SBO transient, a full core model was necessary because of the non-axial symmetric core design and heat loss from the outer core structures to the environment. To reduce the required computational costs, a novel scaled subassembly approach was used to model the full core domain. With the full core steady-state temperature distribution provided as the initial condition, the passive decay heat removal of the core was simulated with the scaled subassembly full core model under SBO transient conditions. This analysis showed that the peak fuel temperatures of the Holos-Quad core remained below their steady-state values during the transient, therefore much lower than the safety limit 1,600°C. This is mainly due to the large thermal inertia of the graphite matrix, relatively low power density, and the large surface-to-volume-ratio of the core. The preliminary analyses presented in this report confirm the inherent safety of the Holos-Quad micro-reactor thermal-hydraulics design, as the peak fuel temperatures under both normal and off-normal operational conditions remain within design and safety limits.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

The Cell Utilized Partitioning Model as a Predictive Tool for Optimizing Counter-Current Chromatography Processes

Counter-current chromatography (CCC) is capable of unique elution modes that isolate analytes using the movement of the stationary phase in addition to moving the mobile phase. These modes include elution-extrusion CCC (EECCC) and dual-mode CCC (DM CCC) that are not possible in traditional solid-liquid chromatography systems. Although EECCC and DM CCC are widely used to recover highly retained components, to our knowledge, optimizing the elution process in these modes with predictive models has not been reported. To address this gap, we developed a predictive model for CCC dubbed the Cell Utilized Partitioning (CUP) model. The CUP model accurately predicts the effluents of multicomponent separations in EECCC and DM CCC modes when compared to experimental data. Furthermore, CUP model simulations were extended to investigate the influence of operating and intrinsic parameters on the yield and productivity, and to compare the separation performances of EECCC and DM CCC in various conditions. The results demonstrate that low distribution constants, usually a KD less than 1, and a selectivity > 1.3, under specific flowrate ranges, increase both productivity and yield. From these results, generalized optimization and scaleup guidelines are proposed that can apply to research settings and to industrial processes to maximize preparative CCC performance.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2023

Bi-Level Linear Programming Model for Automatic Load Shedding: A Distributed Wide-Area Measurement System-based Solution

Load shedding is currently implemented as a two-step based approach. In the first step, manual load shedding is taken place, were system operators, using estimates, inform distribution utilities of predicted stressful conditions. Information provided include the potential use of energy reserves, as well as load shedding amount. In a second step, automatic load shedding is done. The latter is realized using protection relays. While considering frequency variation, pre-defined values of load to be shed and correspondent number of stages for such to be realized are transformed into relay settings. Under-frequency protection relays use only local measurements towards decision making, thus operate in a decentralized architecture. Decision making is done in milliseconds plus breaker time. While this approach has provided much system reliability, considering the new smart grid paradigm, where system dynamics are much faster due to increasing renewable resources penetration, in some operating conditions it will generate sub-optimal solutions, such as islanding. Phasor measurement units provide a source of information which can be useful for this problem. Centralized architecture-based solutions for automatic load shedding, as present in the state-of-the-art, require though total processing times which are not acceptable for real-life implementation. In this work, considering the above, a bi-level linear programming model is presented. The model is implemented considering a distributed architecture while leveraging phasor measurement units data. The upper-level model estimates the current system state. Results of this model are embedded in a lower-level model, which decision variables are the location and load value to be shed. Easy-to-implement model, built-on the classic weighted least squares solution, highlight potential aspects towards real-life applications.

Bretas, Arturo Suman↗
DOE OSTI2025

Measurement of three-dimensional inclusive muon-neutrino charged-current cross sections on argon with the MicroBooNE detector

We report the measurement of the triple-differential cross section d 3 σ/dE vis d cos(θ μ )dP μ for inclusive muon-neutrino charged-current scattering on argon. This measurement utilizes data from 6.4 x 10 20 protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located along the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8 GeV. The mapping from reconstructed kinematics to truth quantities is validated within uncertainties by comparing the distribution of reconstructed hadronic energy in data to that of the model prediction in different muon scattering angle bins after applying a conditional constraint from the muon momentum distribution in data. The success of this validation provides confidence that the energy transfer in the MicroBooNE detector is well-modeled within simulation uncertainties, enabling a reliable unfolding to a triple-differential cross section defined at the nominal neutrino flux over muon momentum, muon scattering angle, and visible neutrino energy. This validation not only supports accurate cross-section extraction, but also establishes a critical foundation for tuning interaction models used in future neutrino oscillation measurements. The unfolded measurement covers an extensive phase space, providing a wealth of information useful for future liquid argon time projection chamber experiments measuring neutrino oscillations. Comparisons against a number of commonly used model predictions are included and their performance in different parts of the available phase-space is discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Numerical framework for integrated additive manufacturing-compression molding (AM-CM) of thermoplastic composites

Additive manufacturing-compression molding (AM-CM) has emerged as a transformative technology in advanced composite manufacturing. Additive manufacturing (AM) offers high design flexibility and the ability to produce complex geometries with precisely aligned fibers in the preferred orientation. Compression molding (CM) enhances composite materials by providing excellent dimensional stability, reduced porosity, high production rates, and a smooth surface finish. Despite these advantages, extensive integrated analysis is required to optimize processing conditions for improved fiber orientation distribution (FOD) and porosity control. Here, this study develops a comprehensive numerical model to simulate the AM-CM manufacturing process. The model isolates the effects of both the AM and CM phases while also capturing their integration. Additionally, it accounts for heat transfer, temperature-dependent viscosity, and fiber orientation in the extruded fiber-filled polymer, accurately representing material behavior during processing. This approach enables the analysis of interactions between deposited beads of complex strand shapes and their interface regions after full compression. Moreover, the model predicts key parameters such as polymer flowability, fiber orientation, and temperature evolution in AM-CM parts. By optimizing processing conditions, it facilitates a controlled and predictable microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparison of Probabilistic Forecasts for Predictive Voltage Control

This paper explores predictive cooperative voltage control in distribution systems with highly variable sources such as photovoltaics (PV). The goal is to maintain the voltage profile within the limits despite the fluctuations due to sudden changes in solar power generation. The predictive voltage control method relies on probabilistic solar power and load forecasts to select the optimum Voltage Regulator (VR) taps appropriately. VR taps are selected to minimize the risk of voltage violation. A modified version of the IEEE 123 system is used as the case study. A 100% penetration of solar power is assumed for the distribution system with profiles for solar generation and loads added to the system. Three different probabilistic forecast models (Quantile Regression (QR), Gaussian distribution and volatility forecasting using Generalized Autoregressive Conditional Heteroskedasticity (GARCH)) are explored in this study. The results for the VR taps and Voltage Deviation Index (VDI) are compared to find the most effective forecast model.

Panamtash, Hossein↗
DOE OSTI2025

Hopping mediated transport between finite pools of redox proteins

Transport reactions in biology involve the flow of particles—electrons, ions, or molecules—between reservoirs. Here, we explore how electron transport between finite reservoirs depends on the nature of the reservoirs, including their size, occupancy, and interactions. We compare the transport kinetics produced by narrowband and wideband infinite reservoir models (described earlier) with a finite narrowband reservoir model. The transport between finite reservoirs is found to depend on both the initial charge distribution and the number of carriers present. Whether or not a steady-state transport regime is accessed prior to reaching the equilibrium charge distribution depends on these initial conditions.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Alpine permafrost could account for a quarter of thawed carbon based on Plio-Pleistocene paleoclimate analogue

Estimates of the permafrost-climate feedback vary in magnitude and sign, partly because permafrost carbon stability in warmer-than-present conditions is not well constrained. Here we use a Plio-Pleistocene lacustrine reconstruction of mean annual air temperature (MAAT) from the Tibetan Plateau, the largest alpine permafrost region on the Earth, to constrain past and future changes in permafrost carbon storage. Clumped isotope-temperatures (Δ 47 -T) indicate warmer MAAT (~1.2 °C) prior to 2.7 Ma, and support a permafrost-free environment on the northern Tibetan Plateau in a warmer-than-present climate. Δ 47 -T indicate ~8.1 °C cooling from 2.7 Ma, coincident with Northern Hemisphere glacial intensification. Combined with climate models and global permafrost distribution, these results indicate, under conditions similar to mid-Pliocene Warm period (3.3–3.0 Ma), ~60% of alpine permafrost containing ~85 petagrams of carbon may be vulnerable to thawing compared to ~20% of circumarctic permafrost. This estimate highlights ~25% of permafrost carbon and the permafrost-climate feedback could originate in alpine areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Investigation of Nanoparticle Degradation in Hydrogen Fuel Cell Systems through Automated Electron Microscopy

Proton exchange membrane fuel cells (PEMFC) are promising devices for the deployment of hydrogen-powered heavy-duty vehicles, providing a higher efficiency for similar driving range and fueling time than the existing ones. However, PEMFCs still encounter durability challenges mainly due to catalyst degradation in the cathode. Mitigating these performance losses requires a better understanding of the degradation mechanisms under heavy-duty accelerated stress tests (ASTs) [1]. Scanning transmission electron microscopy (STEM) combined with energy dispersive X-ray spectroscopy (EDS) are key tools for the analysis of Pt and PtCo nanoparticle size, spatial distribution and composition [2]. Electron tomography is also used to determine the rate and type of degradation of catalyst nanoparticles as a function of their position on the carbon support. In this work, automated data acquisition software, paired with a custom Python code, have been used to study the effect of different accelerated stress tests (ASTs) on nanoparticle coarsening [2]. Figure 1 shows high-angle annular dark-field (HAADF)-STEM images and EDS maps comparing the cathodes of membrane electrode assemblies (MEAs) following an electrocatalyst AST performed under H2/N2 with that of the heavy-duty AST performed under H2/air. We will discuss how AST conditions affect considerably the spatial distribution of the nanoparticles across the electrode between the membrane and microporous layer. Although the median particle size increased more in the MEA aged under the heavy-duty AST, as determined using a high-throughput image analysis, the quantitative EDS measurements demonstrate that the electrocatalyst AST resulted in more Pt and Co dissolution from the cathode, which is another important indicator of electrocatalyst degradation. We will further present the impact of the relative humidity (% RH) on the degradation mechanisms demonstrated using the same approach. Electron tomography has been used to distinguish the Pt nanoparticles residing on the carbon support surface (exterior) from those within the pore structure (interior) in order to determine the relative stability of interior and exterior nanoparticles. As shown in Figure 2, we will compare the Pt catalyst particle size at the beginning of test (BOT) and end of test (EOT), and discuss the importance of automating the electron tomography workflow, i.e. acquisition, reconstruction, and visualization, to increase sampling and determine the standard deviation of these measurement. The outlook for utilizing low-dose cryo-tomography for limiting damage to the catalyst, support, and especially proton-conducting ionomer will also be discussed [3].

Amichi, Lynda↗
DOE OSTI2024

Alteration of phycobilisome excitation energy transfer properties in response to attenuations in peripheral electron flow

In Synechocystis sp. PCC 6803 ( S . 6803), two types of phycobilisome (PBS) complexes, CpcG-PBS and CpcL-PBS, function to harvest light energy for photosynthetic reaction centers (RCs), photosystem I (PSI) and photosystem II (PSII). The compositional differences between these two forms of PBS and their specificity for RCs have led to suggestions that they may differ in function. To address this question, we examined how PBS-RC interactions, and the transfer of excitation energy from PBS to RCs, might be adjusted under conditions where electron demand and photon availability are modulated. The CpcG-PBS, CpcL-PBS, and RC complexes were isolated from a S . 6803 strain defective in expression of flavodiiron 1 (oxygen reduction reaction 1, ORR1) grown under varied light regimes. The energy transfer preference from CpcL-PBS to either PSI or PSII was investigated by in vitro crosslinking and 77 K fluorescence emission spectroscopy to assess energy transfer efficiency under photoexcitation. While the results demonstrate that the transfer of excitation energy from CpcL-PBS favors PSI over PSII in WT strains as previously shown, the preference of CpcL-PBS switches from PSI to PSII in ORR1 strains. Surprisingly, this change in preference was reproduced when ORR1 CpcL-PBS was crosslinked with WT RCs, or when WT CpcL-PBS was cross-crosslinked with ORR1 RCs, indicating there are physical modifications to both PBS and RCs that mediate the preference switch. In contrast, the analysis with ORR1 CpcG-PBS shows similar preferences to WT. Additionally, PBS populations in ORR1 shifted to a greater proportion of CpcL-PBS relative to CpcG-PBS. These results demonstrate that under conditions where electron utilization changes, there is a tuning of the excitation energy allocation from CpcL-PBS to RCs to manage the energy distribution for photosynthesis under dynamic flux conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Li-ion Battery Material phase prediction through Hierarchical Curriculum Learning

Li-ion Batteries (LIB), one of the most efficient energy storage devices, are widely adopted in many industrial applications. Imaging data of these battery electrodes obtained from X-ray tomography can explain the distribution of material constituents and allow reconstructions to study electron transport pathways. Therefore, it can eventually help quantify various associated properties of electrodes (e.g., volume-specific surface area, porosity) which determine the performance of batteries. However, these images often suffer from low image contrast between multiple material constituents , making it difficult for humans to distinguish and characterize these constituents through visualization. A minor error in detecting distributions among the material constituents can lead to a high error in the calculated parameters of material properties.We present a novel hierarchical curriculum learning framework to address the complex task of estimating material constituent distribution in battery electrodes. To provide spatially smooth prediction, our framework comprises three modules: (i) an uncertainty-aware model trained to yield inferences conditioned upon global knowledge of material distribution, (ii) a technique to capture relatively more fine-grained (local) distributional signals, (iii) an aggregator to appropriately fuse the local and global effects towards obtaining the final distribution.

Tabassum, Anika↗
DOE OSTI2026

Marine Aerosol to Refinery Emissions: Transport and Evolution of CCN in the Houston Metropolitan Area and Their Impact on Cloud Formation

The Experiment of Sea Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) campaign aimed to untangle the impacts of contrasting aerosol sources on cloud microphysical properties within deep convective cells using airborne observations with the National Research Council Canada (NRC) Convair CV‐580 research aircraft. Horizontal and vertical gradients of aerosol number size distributions, total aerosol concentrations and cloud condensation nuclei (CCN) spectra were measured in the lower troposphere to quantify the impact of different aerosol sources on aerosol‐cloud interactions. This study focuses on a research flight dedicated to characterizing the aerosol and CCN properties in the Houston Metropolitan region and identifies five main categories of aerosols based on characteristics of aerosol number size distributions, their CCN properties and meteorological conditions. These categories encompassed more than two orders of magnitude differences in aerosol and CCN concentrations, yet their hygroscopic properties remained similar. Aerosol number size distributions and effective hygroscopicity parameters are used to generate continuous CCN spectra to represent the major aerosol types. The different CCN spectra are then incorporated into a 1‐D aerosol‐cloud parcel model using a large range of updrafts selected in the range of those observed during the ESCAPE campaign to assess the impact of the major aerosol sources in Houston on deep convective cloud microphysical properties.

aerosol-cloud interactions↗
DOE OSTI2025

Intrinsic Kinetics of Polyethylene Terephthalate Pyrolysis via Micropyrolysis and Multivariate Chromatographic Analysis

This study provides an in-depth investigation of the primary decomposition of polyethylene terephthalate (PET) via pyrolysis, employing an experimental-analytic workflow that integrates design of experiments (DoE), micropyrolysis coupled with comprehensive two-dimensional gas chromatography (GC×GC), and multivariate data analysis to verify intrinsic kinetic conditions and elucidate evolving product distributions for mapping key reaction pathways. Peaks that could not be identified using commercial spectral libraries were assigned using Mass Frontier simulations, enabling the identification of divinyl terephthalate, ethyl vinyl terephthalate, and 2-(benzoyloxy)ethyl vinyl terephthalate. A polar×polar (non-orthogonal) column set tailored for the detection of carboxylic acids enhanced the quantification of benzoic acid, 4-vinylbenzoic acid, 4-ethylbenzoic acid, and methylbenzoic acid by up to 6-fold relative to an orthogonal column combination (non-polar×mid-polar). Moreover, pyrolysis variables were systematically evaluated using a Box- Behnken design (BBD), encompassing pyrolysis temperature (500−600 °C), sample weight (50−150 μg), and carrier gas flow rate (100−300 mL min −1 ). Among these, pyrolysis temperature was the only statistically significant factor influencing product yields, ranging from 58.78 to 84.26 wt %. In contrast, neither the sample weight nor the carrier gas flow rate had a significant effect on product yields within the evaluated experimental space. At 600 °C, the major pyrolysis products were benzoic acid (up to 20.20 ± 1.46 wt %) and CO 2 (up to 21.28 ± 1.46 wt %), which can be produced through decarboxylation reactions. These findings underscore the critical importance of selecting appropriate analytical columns for the accurate quantification of heteroatomcontaining products such as carboxylic acids, which may otherwise be underestimated or undetected due to their reactivity with the stationary phase of non-polar and mid-polar columns, as well as other GC components. They also highlight the importance of selecting pyrolysis conditions for investigating the primary decomposition of PET under an isothermal kinetically limited regime.

aromatic compounds↗
DOE OSTI2025

Networked Microgrid Topology Reconfiguration to Promote Fairness in Proactive Load Shedding

Increasing occurrences of natural disasters and grid emergency events consistently challenge the safe and reliable operations of power systems. During such emergency situations, system operators may proactively shed load to mitigate risks. However, uncoordinated implementation of load shedding may disrupt electricity supply and even lead to cascading failures. Meanwhile, it is crucial to address potential biases affecting different customers when executing load shedding. This paper addresses the dynamic topology reconfiguration problem for networked microgrids with distributed energy resources under emergency conditions. Specifically, we propose a novel rolling-horizon optimization model that integrates fairness-aware constraints into the networked microgrid topology reconfiguration. Unlike existing approaches that focus solely on efficiency or apply fairness considerations in static settings, our method explicitly incorporates temporal fairness constraints to restrict repeated or excessive load curtailment for load blocks. Moreover, the fairness-aware constraints are specifically developed for the context of dynamic networked microgrid topology reconfiguration, and are designed to be convex or amenable to linear reformulations, which offers a more tractable alternative to traditional models with non-convex formulations. Numerical studies on a modified IEEE 13-bus system and a larger-sized SMART-DS networked microgrid system demonstrate the performance of the proposed algorithm towards more fairness-aware networked microgrid topology reconfiguration decision-making.

24 POWER TRANSMISSION AND DISTRIBUTION↗