Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinated optimization”

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 1,135 records · Page 63

Addressing Human and Organizational Factors in Nuclear Industry Modernization: A Sociotechnically Based Strategic Framework

The modernization of nuclear power plants will require an advanced concept of operations, involving an integrated set of tightly coupled systems in which all stakeholders act in a coordinated manner. For this modernization effort to be enabled, we developed a human and organizational factors approach based on a broad sociotechnical framework. Starting from core human factors principles, we conducted a literature review of the methods and approaches relevant to the modernization problem. These included not only core disciplines such as cognitive systems engineering, systems theoretic accident modeling and processes, human systems integration, resilience engineering, and macroergonomics but also related topics of safety culture and organizational change. From this literature, we developed a conceptual framework centered around the work system with its four interacting components: people, technology, process, and governance. In an effective work system, these four components are jointly optimized according to three systems criteria: efficiency, effectiveness, and safety. System failure may result from excessive emphasis on any one criterion. The actual work of attaining joint optimization in a given work system can be accomplished by utilizing three high-level functions: knowledge elicitation, knowledge representation, and cross-functional integration. Finally, we illustrated the utility of this approach by applying it to practical problems and case studies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Enabling Modular Autonomous Feedback‐Loops in Materials Science through Hierarchical Experimental Laboratory Automation and Orchestration

Abstract Materials acceleration platforms (MAPs) operate on the paradigm of integrating combinatorial synthesis, high‐throughput characterization, automatic analysis, and machine learning. Within a MAP, one or multiple autonomous feedback loops may aim to optimize materials for certain functional properties or to generate new insights. The scope of a given experiment campaign is defined by the range of experiment and analysis actions that are integrated into the experiment framework. Herein, the authors present a method for integrating many actions within a hierarchical experimental laboratory automation and orchestration (HELAO) framework. They demonstrate the capability of orchestrating distributed research instruments that can incorporate data from experiments, simulations, and databases. HELAO interfaces laboratory hardware and software distributed across several computers and operating systems for executing experiments, data analysis, provenance tracking, and autonomous planning. Parallelization is an effective approach for accelerating knowledge generation provided that multiple instruments can be effectively coordinated, which the authors demonstrate with parallel electrochemistry experiments orchestrated by HELAO. Efficient implementation of autonomous research strategies requires device sharing, asynchronous multithreading, and full integration of data management in experimental orchestration, which to the best of the authors’ knowledge, is demonstrated for the first time herein.

36 MATERIALS SCIENCE↗

A Distributed OpenCL Framework using Redundant Computation and Data Replication

Applications written solely in OpenCL or CUDA cannot execute on a cluster as a whole. Most previous approaches that extend these programming models to clusters are based on a common idea: designating a centralized host node and coordinating the other nodes with the host for computation. However, the centralized host node is a serious performance bottleneck when the number of nodes is large. In this paper, we propose a scalable and distributed OpenCL framework called SnuCL-D for large-scale clusters. SnuCL-D's remote device virtualization provides an OpenCL application with an illusion that all compute devices in a cluster are confined in a single node. To reduce the amount of control-message and data communication between nodes, SnuCL-D replicates the OpenCL host program execution and data in each node. We also propose a new OpenCL host API function and a queueing optimization technique that significantly reduce the overhead incurred by the previous centralized approaches. To show the effectiveness of SnuCL-D, we evaluate SnuCL-D with a microbenchmark and eleven benchmark applications on a large-scale CPU cluster and a medium-scale GPU cluster.

Kim, Junghyun↗

Revisiting metastable immiscibility in SiO 2 –Al 2 O 3 : Structure and phase separation of supercooled liquids and glasses

Understanding and controlling liquid–liquid phase separation in aluminosilicates is crucial for optimizing glass properties. However, the metastable nature of aluminosilicates’ phase separation has made it difficult to study experimentally, and uncertainty persists regarding the compositional and temperature extents of the miscibility gap. Here, we present new experimental evidence that suggests a consolute temperature between 1440 and 1590°C and endmember compositions of 7 and 62 mol.% Al 2 O 3 for the phase-separated glasses. Using containerless melt processing, deeply supercooled liquids over the 0–60 mol.% Al 2 O 3 range are probed with in situ small- and wide-angle X-ray scattering, which simultaneously reveals changes in nanoscale density heterogeneity and atomic structure. Correlations between phase separation and atomic coordination environments are compared for liquids and glasses. Pair distribution function analysis shows mean O–(Si + Al) coordination increases with Al 2 O 3 content and decreases with temperature.

36 MATERIALS SCIENCE↗

Understanding cation-disordered rocksalt oxyfluoride cathodes

Partial fluorine (F) substitution into the oxygen (O) sublattice has been shown to improve cycling stability of cation-disordered Li-excess rocksalt oxide (DRX) cathodes. Detailed understanding on failure mechanisms and key optimization knowledge of fluorinated-DRX (F-DRX), however, are lacking. In the present study, we incorporate different amounts of F into a baseline DRX system, Li 1.2 Ti 0.4 Mn 0.4 O 2 (LTMO3.0, 3.0 denotes nominal Li/Mn ratio), and synthesize two oxyfluoride compounds, Li-rich Li 1.3 Ti 0.3 Mn 0.4 O 1.7 F 0.3 (LTMOF3.25) and Mn-rich Li 1.2 Ti 0.2 Mn 0.6 O 1.8 F 0.2 (LTMOF2.0) with an increased and reduced Li/Mn ratio of 3.25 and 2.0, respectively. Through careful monitoring of chemical and structural evolution, we show that cycling-induced changes are manifested not only by Mn reduction and degradation of its local coordination environment, but also by F enrichment and formation of LiF-type of domains on the surface. A “concerted-densification” based failure mechanism, involving atomic-level changes in both transition-metal cationic sublattice and oxygen/fluorine anionic sublattice, is proposed for the degradation in F-DRX cathode materials. The study reveals that increasing F content accompanied by reduced Li/Mn ratio mitigates the degradation process, offering key design strategies in achieving balanced cathode capacity and stability.

25 ENERGY STORAGE↗

Free-energy calculations for semi-flexible macromolecules: Applications to DNA knotting and looping

We present a method to obtain numerically accurate values of configurational free energies of semiflexible macromolecular systems, based on the technique of thermodynamic integration combined with normal-mode analysis of a reference system subject to harmonic constraints. Compared with previous free-energy calculations that depend on a reference state, our approach introduces two innovations, namely, the use of internal coordinates to constrain the reference states and the ability to freely select these reference states. As a consequence, it is possible to explore systems that undergo substantially larger fluctuations than those considered in previous calculations, including semiflexible biopolymers having arbitrary ratios of contour length L to persistence length P. To validate the method, high accuracy is demonstrated for free energies of prime DNA knots with L/P = 20 and L/P = 40, corresponding to DNA lengths of 3000 and 6000 base pairs, respectively. We then apply the method to study the free-energy landscape for a model of a synaptic nucleoprotein complex containing a pair of looped domains, revealing a bifurcation in the location of optimal synapse (crossover) sites. This transition is relevant to target-site selection by DNA-binding proteins that occupy multiple DNA sites separated by large linear distances along the genome, a problem that arises naturally in gene regulation, DNA recombination, and the action of type-II topoisomerases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

4D beam matrix reconstruction in particle accelerators

Transverse beam parameters in particle accelerators are commonly described using the Twiss parameters, which are experimentally accessible yet inherently limited because they neglect correlations between different transverse coordinates. Such correlations frequently arise from uncompensated cathode magnetic fields or misaligned focusing quadrupoles, affecting beam quality and accelerator performance. To address this limitation, we propose and validate a novel diagnostic method for the complete four-dimensional (4D) transverse beam matrix. Our method involves passing the beam through a beamline comprising both conventional and skew quadrupole magnets, followed by downstream measurements of the resulting two-dimensional (2D) beam profiles. These measurements represent distinct 2D projections of the underlying 4D transverse phase–space distribution. By systematically varying quadrupole strengths, multiple independent projections of the beam phase space are obtained. We reconstruct the original 4D beam matrix from these measured projections using an optimization-based least-square fit, providing fast and robust reconstruction regardless of the specific beamline configuration. Through extensive numerical simulations and realistic particle-tracking studies, we demonstrate the diagnostic’s accuracy, robustness, and capability to achieve reconstruction uncertainties smaller than measurement errors, particularly when employing sufficient numbers of quadrupole scans. This method presents a powerful and flexible approach for comprehensive beam characterization and accelerator tuning.

43 PARTICLE ACCELERATORS↗

THE COORDINATED RADIO AND INFRARED SURVEY FOR HIGH-MASS STAR FORMATION. II. SOURCE CATALOG

The CORNISH project is the highest resolution radio continuum survey of the Galactic plane to date. It is the 5 GHz radio continuum part of a series of multi-wavelength surveys that focus on the northern GLIMPSE region (10 Degree-Sign < l < 65 Degree-Sign ), observed by the Spitzer satellite in the mid-infrared. Observations with the Very Large Array in B and BnA configurations have yielded a 1.''5 resolution Stokes I map with a root mean square noise level better than 0.4 mJy beam{sup -1}. Here we describe the data-processing methods and data characteristics, and present a new, uniform catalog of compact radio emission. This includes an implementation of automatic deconvolution that provides much more reliable imaging than standard CLEANing. A rigorous investigation of the noise characteristics and reliability of source detection has been carried out. We show that the survey is optimized to detect emission on size scales up to 14'' and for unresolved sources the catalog is more than 90% complete at a flux density of 3.9 mJy. We have detected 3062 sources above a 7{sigma} detection limit and present their ensemble properties. The catalog is highly reliable away from regions containing poorly sampled extended emission, which comprise less than 2% of the survey area. Imaging problems have been mitigated by down-weighting the shortest spacings and potential artifacts flagged via a rigorous manual inspection with reference to the Spitzer infrared data. We present images of the most common source types found: H II regions, planetary nebulae, and radio galaxies. The CORNISH data and catalog are available online at http://cornish.leeds.ac.uk.

79 ASTRONOMY AND ASTROPHYSICS↗

Coordinated platooning with multiple speeds

In a platoon, vehicles travel one after another with small intervehicle distances; trailing vehicles in a platoon save fuel because they experience less aerodynamic drag. This work presents a coordinated platooning model with multiple speed options that integrates scheduling, routing, speed selection, and platoon formation/dissolution in a mixed-integer linear program that minimizes the total fuel consumed by a set of vehicles while traveling between their respective origins and destinations. The performance of this model is numerically tested on a grid network and the Chicago-area highway network. We find that the fuel-savings factor of a multivehicle system significantly depends on the time each vehicle is allowed to stay in the network; this time affects vehicles’ available speed choices, possible routes, and the amount of time for coordinating platoon formation. For problem instances with a large number of vehicles, we propose and test a heuristic decomposed approach that applies a clustering algorithm to partition the set of vehicles and then routes each group separately. When the set of vehicles is large and the available computational time is small, the decomposed approach finds significantly better solutions than does the full model.

33 ADVANCED PROPULSION SYSTEMS↗

High-throughput screening for discovery of benchtop separations systems for selected rare earth elements

Abstract Rare earth (RE) elements (scandium, yttrium, and the lanthanides) are critical for their role in sustainable energy technologies. Problems with their supply chain have motivated research to improve separations methods to recycle these elements from end of life technology. Toward this goal, we report the synthesis and characterization of the ligand tris[(1-hydroxy-2-oxo-1,2-dihydropyridine-3-carboxamido)ethyl]amine, H 3 1·TFA (TFA = trifluoroacetic acid), and complexes 1·RE (RE = La, Nd, Dy). A high-throughput experimentation (HTE) screen was developed to quantitatively determine the precipitation of 1·RE as a function of pH as well as equivalents of H 3 1·TFA. This method rapidly determines optimal conditions for the separation of RE mixtures, while minimizing materials consumption. The HTE-predicted conditions are used to achieve the lab-scale separation of Nd/Dy ( SF Nd/Dy = 213 ± 34) and La/Nd ( SF La/Nd = 16.2 ± 0.2) mixtures in acidic aqueous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-E-J-148: Tools for Development of 4D Proton CT

Purpose: To develop tools for performing 4D proton computed tomography (CT). Methods: A suitable patient with a tumor in the right lower lobe was selected from a set of 4D CT scans. The volumetric CT images formed the basis for calculating the parameters of a breathing model that allows reconstruction of a static reference CT and CT images in each breathing phase. The images were imported into the TOPAS Monte Carlo simulation platform for simulating an experimental proton CT scan with 45 projections spaced by 4 degree intervals. Each projection acquired data for 2 seconds followed by a gantry rotation for 2 seconds without acquisition. The scan covered 180 degrees with individual protons passing through a 9-cm slab of the patient’s lung covering the moving tumor. An initial proton energy sufficient for penetrating the patient from all directions was determined. Performing the proton CT simulation, TOPAS provided output of the proton energy and coordinates registered in two planes before and after the patient, respectively. The set of projection data was then used with an iterative reconstruction algorithm to generate a volumetric proton CT image set of the static reference image and the image obtained under breathing motion, respectively. Results: An initial proton energy of 230 MeV was found to be sufficient, while for an initial energy of 200 MeV a substantial number of protons did not penetrate the patient. The reconstruction of the static reference image set provided sufficient detail for treatment planning. Conclusion: We have developed tools to perform studies of proton CT in the presence of lung motion based on the TOPAS simulation toolkit. This will allow to optimize 4D reconstruction algorithms by synchronizing the acquired proton CT data with a breathing signal and utilizing a breathing model obtained prior to the proton CT scan. This research has been supported by the National Institute Of Biomedical Imaging And Bioengineering of the National Institutes of Health under Award Number R01EB013118.

60 APPLIED LIFE SCIENCES↗

A dynamic protein interactome drives energy conservation and electron flux in Thermococcus kodakarensis

ABSTRACT Life is supported by energy gains fueled by catabolism of a wide range of substrates, each reliant on the selective partitioning of electrons through redox ( red uction and ox idation) reactions. Electron flux through tunable and regulated protein interactions provides dynamic routes for energy conservation, but how electron flux is regulated in vivo , particularly for archaeal metabolisms that support rapid growth at the thermodynamic limits of life, is poorly understood. Identification of bona fide in vivo protein assemblies and how such assemblies dictate the totality of electron flux is critical to our understanding of the regulation imposed on metabolism, energy production, and energy conservation. Here, 25 key proteins in central metabolic redox pathways in the model, genetically accessible, hyperthermophilic archaeon Thermococcus kodakarensis , were purified to reveal an extensive, dynamic, and tightly interconnected network of protein interactions that responds to environmental cues (such as the availability of various reductive sinks) to direct electron flux to maximize energetic gains. Interactions connecting disparate functions suggest many catabolic and anabolic activities occur in spatial proximity in vivo , and while protein complexes have been historically defined under optimal conditions, many of these complexes appear to maintain alternative partnerships in changing conditions. The totality of the results obtained redefines our understanding of in vivo assemblies driving ancient metabolic strategies supporting the growth of modern Archaea. IMPORTANCE Given the potential for rational genetic manipulations of biofuel- and biotech-promising archaea to yield transformative results for major markets, it is a priority to define how the metabolisms of such species are controlled, at least in part, by in vivo protein assemblies, and from such, define routes of energy flux that can be most efficiently altered toward biofuel or biotechnological gains. Proteinaceous electron carriers (PECs, such as ferredoxins) offer the potential for specific protein–protein interactions to coordinate selective reductive flow. Employing the model, genetically accessible, hyperthermophilic archaeon, Thermococcus kodakarensis , we establish the metabolic protein interactome of 25 key redox proteins, revealing that each redox active protein has a dynamic partnership profile, suggesting catabolic and anabolic activities may occur in concert and in temporal and spatial proximity in vivo . These results reveal critical importance in evaluating the newly identified partnerships and their role and utility in providing regulated redox flux in T. kodakarensis .

Williams, Sere A. (ORCID:0000000235509590)↗

Optimizing Control for Efficient Load Shifting with Thermal Energy Storage in Existing HVAC Systems

This project developed the integration of a direct-contact heat exchanger (DCHX) based thermal energy storage (TES) system with a chiller–air handling unit (AHU) plant to evaluate its potential for reducing building energy costs. Detailed physical models of the TES unit, building envelope, and HVAC components were developed alongside simplified control-oriented models to support both high-fidelity simulation and real-time optimization. Two control strategies were implemented and compared: a rule-based control (RBC) aligned with utility time-of-use (ToU) rates, and a model predictive control (MPC) framework leveraging forecasts of building load, weather, and internal gains.Simulation results show that the RBC strategy reduced daily electricity costs by around 30% by shifting cooling production from on-peak to off-peak hours. In contrast, the MPC strategy achieved significantly greater performance, reducing daily operating costs by up to 44% and peak-hour costs by more than 60%. Both strategies maintained indoor thermal comfort within acceptable limits, with MPC further improving load distribution and reducing equipment cycling.The outcomes confirm that TES integration, particularly when coordinated with advanced predictive control, can provide substantial cost savings and on-peak demand reduction. . These findings directly support the U.S. Department of Energy’s goals for grid-interactive efficient buildings and demonstrate the potential of TES-enabled HVAC systems for scalable deployment across the commercial building.

99 GENERAL AND MISCELLANEOUS↗

Transfer and Transition: Interagency Coordination for Managing Public Lands at UMTRCA Title II Sites in Wyoming – 16614

By the end of fiscal year 2025, the U.S. Department of Energy (DOE) Office of Legacy Management (LM) is anticipating adding 17 sites remediated under Title II of the Uranium Mill Tailings Radiation Control Act (UMTRCA) to the current inventory of 90 sites that it manages. Among the new sites are ones where federal public lands occur within the proposed long-term care boundary, the boundary determined by the Nuclear Regulatory Commission and LM as necessary to maintain site protectiveness for the entombed uranium mill tailings and residual groundwater contamination. For these sites, public land withdrawals for land and minerals will need to be established. LM’s primary mission at UMTRCA sites is to protect the public and the environment from exposure to contamination at the sites. For the sites with public lands or federally controlled minerals that will be transferring to LM, the Office will apply to the Department of the Interior (DOI) Bureau of Land Management (BLM) for new, public land and mineral withdrawals. At most current LM UMTRCA sites that involved public lands and minerals, DOI granted DOE “full administrative jurisdiction” and permanent withdrawals. Hence, these withdrawals are, permanently, no longer subject to public land, mining, and mineral-leasing laws and regulations. LM is coordinating with DOI/BLM in Wyoming to permanently withdraw full and partial jurisdiction at future UMTRCA Title II sites in that state. This approach would allow LM to fully administer surface lands and minerals, where necessary, and DOI and LM to administer surface lands and leasable minerals where it would not jeopardize sites’ radiological safety and long-term public and environmental protection. This “shared-jurisdiction approach” will meet LM’s strategic goal of protecting human health and the environment but also allow BLM to fulfill their mission to “manage and conserve the lands under the mandate of multiple-use and sustained yield.” In addition, LM could also fulfill the fourth goal of its Strategic Plan, to optimize land use and assets.

54 ENVIRONMENTAL SCIENCES↗

Transfer and Transition: Interagency Coordination for Managing Public Lands at UMTRCA Title II Sites in Wyoming - 16614

By the end of fiscal year 2025, the U.S. Department of Energy (DOE) Office of Legacy Management (LM) is anticipating adding 17 sites remediated under Title II of the Uranium Mill Tailings Radiation Control Act (UMTRCA) to the current inventory of 90 sites that it manages. Among the new sites are ones where federal public lands occur within the proposed long-term care boundary, the boundary determined by the Nuclear Regulatory Commission and LM as necessary to maintain site protectiveness for the entombed uranium mill tailings and residual groundwater contamination. For these sites, public land withdrawals for land and minerals will need to be established. LM's primary mission at UMTRCA sites is to protect the public and the environment from exposure to contamination at the sites. For the sites with public lands or federally controlled minerals that will be transferring to LM, the Office will apply to the Department of the Interior (DOI) Bureau of Land Management (BLM) for new, public land and mineral withdrawals. At most current LM UMTRCA sites that involved public lands and minerals, DOI granted DOE 'full administrative jurisdiction' and permanent withdrawals. Hence, these withdrawals are, permanently, no longer subject to public land, mining, and mineral-leasing laws and regulations. LM is coordinating with DOI/BLM in Wyoming to permanently withdraw full and partial jurisdiction at future UMTRCA Title II sites in that state. This approach would allow LM to fully administer surface lands and minerals, where necessary, and DOI and LM to administer surface lands and leasable minerals where it would not jeopardize sites' radiological safety and long-term public and environmental protection. This 'shared-jurisdiction approach' will meet LM's strategic goal of protecting human health and the environment but also allow BLM to fulfill their mission to 'manage and conserve the lands under the mandate of multiple-use and sustained yield'. In addition, LM could also fulfill the fourth goal of its Strategic Plan, to optimize land use and assets. The portions of the sites where LM will ask for full jurisdiction are those that contain the disposal cell and any adjacent lands and minerals required to maintain and protect it. It is anticipated that the shared jurisdiction portion of the withdrawal would largely coincide with land where known or projected contaminated groundwater from the past ore processing exists or will migrate over time. Examples of BLM-authorized activities on the shared jurisdiction portions of the withdrawals may include livestock grazing, wildlife habitat improvement projects, recreation, and rights of ways for roads, pipelines, utilities, and communication sites. In addition, BLM could authorize subsurface uses such as oil and gas, coal, or geothermal development or groundwater pumping from deeper aquifers if the proposal would not disturb groundwater contamination. For areas of the withdrawal where the agencies will share jurisdictional authority, DOE LM and BLM will develop an interagency agreement (IA) to define acceptable uses and associated roles and responsibilities. To ensure maximum protection of human health and the environment, each site's IA will include a provision allowing DOE LM the right of first refusal against any action it believes might affect radiological safety at the site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fluorinated Covalent Organic Frameworks: A Novel Pathway to Enhance Hydrogen Sorption and Control Isosteric Heats of Adsorption; HyMARC Seed Project Final Report

Organic framework materials, as a unique sub-class of carbon-based sorbents, have gained increasing attention for promising attributes toward gas storage. Calculations of the hydrogen capacity of optimized frameworks show potential to achieve greater than 60 g/L storage of hydrogen, placing metal-organic frameworks (MOFs) and COFs near the top of the class of porous materials. However, at the ensemble level, poor stacking creates a quasi-amorphous material with low structural integrity and low effective surface area. The chemical versatility of COFs allows for additional methods for producing long-range order that specifically target the inter-layer interactions in COFs but that leave key pore-accessible linker sites open. These modifications (either during synthesis or post-synthetic) enable a host of strategies to both improve crystalline order for more stable and higher surface, as well as add metals with open coordination for enhanced H 2 binding enthalpy.

08 HYDROGEN↗

Optimizing Control for Efficient Load Shifting with Thermal Energy Storage in Existing HVAC Systems

This project developed the integration of a direct-contact heat exchanger (DCHX) based thermal energy storage (TES) system with a chiller–air handling unit (AHU) plant to evaluate its potential for reducing building energy costs. Detailed physical models of the TES unit, building envelope, and HVAC components were developed alongside simplified control-oriented models to support both high-fidelity simulation and real-time optimization. Two control strategies were implemented and compared: a rule-based control (RBC) aligned with utility time-of-use (ToU) rates, and a model predictive control (MPC) framework leveraging forecasts of building load, weather, and internal gains. Simulation results show that the RBC strategy reduced daily electricity costs by around 30% by shifting cooling production from on-peak to off-peak hours. In contrast, the MPC strategy achieved significantly greater performance, reducing daily operating costs by up to 44% and peak-hour costs by more than 60%. Both strategies maintained indoor thermal comfort within acceptable limits, with MPC further improving load distribution and reducing equipment cycling. The outcomes confirm that TES integration, particularly when coordinated with advanced predictive control, can provide substantial cost savings and on-peak demand reduction. . These findings directly support the U.S. Department of Energy’s goals for grid-interactive efficient buildings and demonstrate the potential of TES-enabled HVAC systems for scalable deployment across the commercial building.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Change-point Detection and Image Segmentation for Time Series of Astrophysical Images

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multiband image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all of the change points (in time domain), as well as all of the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real data sets show that the proposed method enjoys very promising empirical properties. Applications to two real data sets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it.

79 ASTRONOMY AND ASTROPHYSICS↗