Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “topological data analysis”

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 577 records · Page 32

Challenges to Computational Aerothermodynamic Simulation and Validation for Planetary Entry Vehicle Analysis

Challenges to computational aerothermodynamic (CA) simulation and validation of hypersonic flow over planetary entry vehicles are discussed. Entry, descent, and landing (EDL) of high mass to Mars is a significant driver of new simulation requirements. These requirements include simulation of large deployable, flexible structures and interactions with reaction control system (RCS) and retro-thruster jets. Simulation of radiation and ablation coupled to the flow solver continues to be a high priority for planetary entry analyses, especially for return to Earth and outer planet missions. Three research areas addressing these challenges are emphasized. The first addresses the need to obtain accurate heating on unstructured tetrahedral grid systems to take advantage of flexibility in grid generation and grid adaptation. A multi-dimensional inviscid flux reconstruction algorithm is defined that is oriented with local flow topology as opposed to grid. The second addresses coupling of radiation and ablation to the hypersonic flow solver--flight- and ground-based data are used to provide limited validation of these multi-physics simulations. The third addresses the challenges of retro-propulsion simulation and the criticality of grid adaptation in this application. The evolution of CA to become a tool for innovation of EDL systems requires a successful resolution of these challenges.

Gnoffo, Peter A.↗

Topological analysis of the CfA redshift survey

We study the topology of large-scale structure in the Center for Astrophysics Redshift Survey, which now includes approximately 12,000 galaxies with limiting magnitude m(sub B) is less than or equal to 15.5. The dense sampling and large volume of this survey allow us to compute the topology on smoothing scales from 6 to 20/h Mpc; we thus examine the topology of structure in both 'nonlinear' and 'linear' regimes. On smoothing scales less than or equal to 10/h Mpc this sample has 3 times the number of resolution elements of samples examined in previous studies. Isodensity surface of the smoothed galaxy density field demonstrate that coherent high-density structures and large voids dominate the galaxy distribution. We compute the genus-threshold density relation for isodensity surfaces of the CfA survey. To quantify phase correlation in these data, we compare the CfA genus with the genus of realizations of Gaussian random fields with the power spectrum measured for the CfA survey. On scales less than or equal to 10/h Mpc the observed genus amplitude is smaller than random phase (96% confidence level). This decrement reflects the degree of phase coherence in the observed galaxy distribution. In other words the genus amplitude on these scales is not good measure of the power spectrum slope. On scales greater than 10/h Mpc, where the galaxy distribution is rougly in the 'linear' regime, the genus ampitude is consistent with the random phase amplitude. The shape of the genus curve reflects the strong coherence in the observed structure; the observed genus curve appears broader than random phase (94% confidence level for smoothing scales less than or equal to 10/h Mpc) because the topolgoy is spongelike over a very large range of density threshold. This departre from random phase consistent with a distribution like a filamentary net of 'walls with holes.' On smoothing scales approaching approximately 20/h Mpc the shape of the CfA genus curve is consistent with random phase. There is very weak evidence for a shift of the genus toward a 'bubble-like' topology. To test cosmological models, we compute the genus for mock CfA surveys drawn from large (L greater than or approximately 400/h Mpc) N-body simulations of three variants of the cold dark matter (CDM) cosmogony. The genus amplitude of the 'standard' CDM model (omega h = 0.5, b = 1.5) differs from the observations (96% confidence level) on smoothing scales is less than or approximately 10/h Mpc. An open CDM model (omega h = 0.2) and a CDM model with nonzero cosmological constant (omega h = 0.24, lambda (sub 0) = 0.6) are consistent with the observed genus amplitude over the full range of smoothing scales. All of these models fail (97% confidence level) to match the broadness of the observed genus curve on smoothing scales is less than or equal to 10/h Mpc.

Vogeley, Michael S.↗

Flight Readiness of Mochii S: Portable Spectroscopic Scanning Electron Microscope Facility on the International Space Station (ISS)

The ISS (International Space Station) currently lacks the capability to image and chemically analyze nano-to-micron scale particles from numerous engineering systems. To identify these particles, we must wait for a re-entry vehicle to return them from low earth orbit for ground-based SEM (Scanning Electron Microscope) / EDS (Energy Dispersive X-Ray Spectroscopy) analysis. This may take months, potentially delaying the affected system. Having an EDS-equipped SEM (Mochii S) aboard the ISS will accelerate response time thereby enhancing crew and vehicle safety by rapid and accurate identification of microscopic threats, especially in time-critical situations.The Mochii S payload will be stationed in the Japanese Experiment Module (JEM) powered by 120 VAC (Volts Alternating Current) inverter and connected to station Ethernet and WiFi (Fig. 1). To date the Mochii S payload has undergone testing for command and data handling, power quality, flight vibration, and radiation testing at Johnson Space Center (JSC). Mochii's high-RPM (Revolutions Per Minute) rotating vacuum pumps and high voltage systems have been reviewed to meet safety standards by JSC (Johnson Space Center) Engineering. Topology of the system in the JEM module has been baselined by ISS Safety and JAXA (Japan Space Exploration Agency). Digital controls to and from ISS over Joint Station LAN (Local Area Network) uplink have been simulated and the latencies and data rates have been found to be sufficient for successful operation of the payload from ground.Transporting sensitive electron optical instruments aboard a rocket that sustains 7G acceleration for 8 minutes and then operating it the unique microgravity (micro-g) environment is no trivial matter. To meet strict safety requirements and increase robustness for mission success, over 500 unique verifications must be completed before the payload is certified for spaceflight. Two of which will be discussed in detail are: vibroacoustic testing and magnetic susceptibility shielding and validation.

Mochii↗

Data analysis techniques, differential cross sections, and spin density matrix elements for the reaction γp → Φp

High-statistics measurements of differential cross sections and spin density matrix elements for the reaction γ p → Φp have been made using the CLAS detector at Jefferson Lab. We cover center-of-mass energies (√s) from 1.97 to 2.84 GeV, with an extensive coverage in the Φ production angle. The high statistics of the data sample made it necessary to carefully account for the interplay between the Φ natural lineshape and effects of the detector resolution, that are found to be comparable in magnitude. We study both the charged- (Φ → K⁺K⁻) and neutral- (Φ → K 0 S K 0 L ) $K\bar{K}$ decay modes of the Φ. Further, for the charged mode, we differentiate between the cases where the final K⁻ track is directly detected or its momentum reconstructed as the total missing momentum in the event. The two charged-mode topologies and the neutral-mode have different resolutions and are calibrated against each other. Extensive usage is made of kinematic fitting to improve the reconstructed Φ mass resolution. Our final results are reported in 10- and mostly 30-MeV-wide √s bins for the charged- and the neutral-mode, respectively. Possible effects from K⁺Λ* channels with p$K\bar{K}$ final-states are discussed. These present results constitute the most precise and extensive Φ photoproduction measurements to date and in conjunction with the ω photoproduction results recently published by CLAS, will greatly improve our understanding of low energy vector meson photoproduction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analysis and Mitigation of Cascading Failures Using a Stochastic Interaction Graph with Eigen-analysis

In studies on complex network systems using graph theory, eigen-analysis is typically performed on an undirected graph model of the network. However, when analyzing cascading failures in a power system, the interactions among failures suggest the need for a directed graph beyond the topology of the power system to model directions of failure propagation. To accurately quantify failure interactions for effective mitigation strategies, this paper proposes a stochastic interaction graph model and associated eigen-analysis. Different types of modes on failure propagations are defined and characterized by the eigenvalues of a stochastic interaction matrix, whose absolute values are unity, zero, or in between. Finding and interpreting these modes helps identify the probable patterns of failure propagation, either local or widespread, and the participating components based on eigenvectors. Then, by lowering the failure probabilities of critical components highly participating in a mode of widespread failures, cascading can be mitigated. Here, the validity of the proposed stochastic interaction graph model, eigen-analysis and the resulting mitigation strategies is demonstrated using simulated cascading failure data on an NPCC 140-bus system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

High Energy Yield Distributed Architecture for Large Commercial and Utility-Scale PV Systems

Photovoltaic energy has come a long way from its early days when small cells were fabricated in the labs in the early 50s. While it took a long time for PV modules to mature into what we now see in the market and in rooftops across the world, the power electronics has been lagging behind & was always the weakest point in the PV plant, being responsible for the majority of failures & robbing the PV plants of their full potential by losing precious Watts in the process of converting power from DC to AC. Since the early days of PV power, it has been recognized that the modules need to operate at their maximum power point in order to extract as much power as is available from the PV modules at a given time of day and year. Numerous maximum power point devices were developed to do just this. While effective at extracting as much power as they can, soon after there was a realization that in a PV installation where a central inverter is used, many losses occur, hence energy yield losses. This has led to thinking about alternative architectures where the maximum power point tracking (MPPT) is performed much closer to the modules rather than at a central inverter. Several architectures were proposed ranging from a microinverter, microconverter, string inverter, string optimizer, and to a multi-string inverter or optimizer. While these architectures promise to increase the energy yield by 5% to 25% (depending where MPPT is located), & recovering losses due to shading, soiling, mismatch … there were no real studies and implementations in the field to corroborate these numbers. The program that GE Global Research proposed & received funding for is aimed at a type of distributed architectures, namely the dc architectures where only the MPPT function is distributed across the PV plant. Further, the project is concerned with choosing the right architecture for large scale commercial & utility scale plants but not residential installations. Findings from Task 1- Define Requirements & Develop Architecture- showed that for large scale commercial & utility PV plants, locating the MPPT at the string level with a string size of >3kW is the most cost effective solution & that the increase in energy yield of 6% to 8% is more than enough to make up for the increased cost due to additional power electronics. Our study also shows that in addition to an increase in energy yield, having a dc/dc converter at the string level brings additional benefits such as increased array visibility, allowing for increased monitoring and diagnostics; high and constant dc bus voltage which allows for reducing losses further and optimizing the inverter; and ability to detect ground and arc faults. Task 2 focused on designing a suitable dc/dc converter that meets the requirements defined in Task 1: Efficiency, Cost, & Reliability. For this task, we surveyed and simulated all candidate topologies known by the practitioners in the field. We also generated a novel topology, simulated it, & compared it to existing topologies. Our novel topology being a partial power processing converter came out on top due to its simplicity (fewer parts), it high efficiency, its reliability, and its lower cost. In the remaining tasks of the project, we built the hardware, tested it first in the lab, then in the field and once we worked out all the bugs and issues and were confident our system was robust, we installed it on the rooftop and started collecting the data by mid-Q3 2012. We shared our work with the scientific community in the form of conference & Journal papers, hence contributing to the general understanding of this technology area. The methods, techniques, and hardware developed under this project have all been tested in the field over long periods of time, and under many conditions and have been shown to be robust and reliable. The converters built and developed under this program have been in operation in the field since mid-June 2012, and no hardware failures have been recorded. While some software issues have been seen, none have been showstoppers. The cost analysis of the hardware showed that it meets the cost target set in the requirement document developed in the trade-off study. Improving the performance of solar installations for commercial and utility scale plants contributes to reducing their cost of ownership, hence spreading them to a larger section of the society and making them affordable. Our project contributed to the understanding of the distributed architectures, showed that PV plants’ cost of electricity can be reduced through additional power electronics which allows for harvesting more energy from the plant. When the project started in 2009, the cost of PV energy was still high, hence there was a real need to increase the energy yield, hence reduce the levelized cost of electricity (LCOE). While costs have come down significantly, the need for increasing the energy yield is still relevant, hence the project results, hardware, and techniques developed are still pertinent today. By taking this project from concept to rooftop demonstration, we showed that the technology is mature enough, robust enough, and has the potential to be used by many more customers, hence benefiting the public in general through the spread of PV energy, which ultimately contributes to reducing our dependency on fossil fuel, improving the environment, and creating more jobs.

14 SOLAR ENERGY↗

Nonlinear force-free extrapolation of the coronal magnetic field based on the magnetohydrodynamic relaxation method

We develop a nonlinear force-free field (NLFFF) extrapolation code based on the magnetohydrodynamic (MHD) relaxation method. We extend the classical MHD relaxation method in two important ways. First, we introduce an algorithm initially proposed by Dedner et al. to effectively clean the numerical errors associated with ∇ · B . Second, the multigrid type method is implemented in our NLFFF to perform direct analysis of the high-resolution magnetogram data. As a result of these two implementations, we successfully extrapolated the high resolution force-free field introduced by Low and Lou with better accuracy in a drastically shorter time. We also applied our extrapolation method to the MHD solution obtained from the flux-emergence simulation by Magara. We found that NLFFF extrapolation may be less effective for reproducing areas higher than a half-domain, where some magnetic loops are found in a state of continuous upward expansion. However, an inverse S-shaped structure consisting of the sheared and twisted loops formed in the lower region can be captured well through our NLFFF extrapolation method. We further discuss how well these sheared and twisted fields are reconstructed by estimating the magnetic topology and twist quantitatively.

79 ASTRONOMY AND ASTROPHYSICS↗

HLPW-4: Wall-Modeled LES and Lattice-Boltzmann Technology Focus Group Workshop Summary

A summary of the nine submissions to the Wall-Modeled LES and Lattice-Boltzmann(WMLESLB) Technical Focus Group (TFG) at the 4th High lift Prediction Workshop is provided. The focus of this TFG was to assess the current capabilities of WMLES and LB methods on a complex high-lift configuration across a wide range of angles-of-attack. Analysis of the submitted data suggests that>250𝑀spatial degrees of freedom are needed to accurately predict pitching moments at high angles-of-attack due to large pressure gradients present on the outboard slat and main element for 𝛼 >17◦(corrected for free-air). While some scatter is reported in pitching moment coefficient at the low-angles of attack (𝛼 <11◦) - excellent agreement is observed between submissions near the 𝐶𝐿,max state. Objective superiority of WMLES methods over RANS can be seen in terms of lack of excess outboard separation; a majority of the WMLES and LB submissions predict wedge-shaped separation patterns consistent with the experimental oil flow. The in-tunnel simulations show excellent agreement with the experiment in terms of a) integrated loads, b) surface flow-topology, and c) mechanism for the onset of inboard stall. Further evidence is provided to demonstrate both qualitative and quantitative superiority of the WMLES submissions over RANS.

TTT↗

STUDY OF THE THREE-DIMENSIONAL CORONAL MAGNETIC FIELD OF ACTIVE REGION 11117 AROUND THE TIME OF A CONFINED FLARE USING A DATA-DRIVEN CESE-MHD MODEL

We apply a data-driven magnetohydrodynamics (MHD) model to investigate the three-dimensional (3D) magnetic field of NOAA active region (AR) 11117 around the time of a C-class confined flare that occurred on 2010 October 25. The MHD model, based on the spacetime conservation-element and solution-element scheme, is designed to focus on the magnetic field evolution and to consider a simplified solar atomsphere with finite plasma {beta}. Magnetic vector-field data derived from the observations at the photosphere is inputted directly to constrain the model. Assuming that the dynamic evolution of the coronal magnetic field can be approximated by successive equilibria, we solve a time sequence of MHD equilibria based on a set of vector magnetograms for AR 11117 taken by the Helioseismic and Magnetic Imager on board the Solar Dynamic Observatory around the time of the flare. The model qualitatively reproduces the basic structures of the 3D magnetic field, as supported by the visual similarity between the field lines and the coronal loops observed by the Atmospheric Imaging Assembly, which shows that the coronal field can indeed be well characterized by the MHD equilibrium in most cases. The magnetic configuration changes very little during the studied time interval of 2 hr. A topological analysis reveals that the small flare is correlated with a bald patch (BP, where the magnetic field is tangent to the photosphere), suggesting that the energy release of the flare can be understood by magnetic reconnection associated with the BP separatrices. The total magnetic flux and energy keep increasing slightly in spite of the flare, while the computed magnetic free energy drops during the flare by {approx}10{sup 30} erg, which seems to be adequate in providing the energy budget of a minor C-class confined flare.

79 ASTRONOMY AND ASTROPHYSICS↗

DISTRIBUTION OF PLASMOIDS IN POST-CORONAL MASS EJECTION CURRENT SHEETS

Recently, the fragmentation of a current sheet in the high-Lundquist-number regime caused by the plasmoid instability has been proposed as a possible mechanism for fast reconnection. In this work, we investigate this scenario by comparing the distribution of plasmoids obtained from Large Angle and Spectrometric Coronagraph (LASCO) observational data of a coronal mass ejection event with a resistive magnetohydrodynamic simulation of a similar event. The LASCO/C2 data are analyzed using visual inspection, whereas the numerical data are analyzed using both visual inspection and a more precise topological method. Contrasting the observational data with numerical data analyzed with both methods, we identify a major limitation of the visual inspection method, due to the difficulty in resolving smaller plasmoids. This result raises questions about reports of log-normal distributions of plasmoids and other coherent features in the recent literature. Based on nonlinear scaling relations of the plasmoid instability, we infer a lower bound on the current sheet width, assuming the underlying mechanism of current sheet broadening is resistive diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗

WISP: Watching grid Infrastructure Stealthily through Proxies (Final Technical Report)

The complex interdependencies of cyber systems (sensors and communications), physical grids and associated electricity market operations make protecting electric power grids a significant challenge. The energy sector is constantly under new, targeted, advanced and dangerous cyber-attacks that have the potential to result in the loss of human life. These threats are further exacerbated by our need to modernize the grid. One focus of cyber security research in smart grids is the securing of the SCADA system through advanced intrusion detection systems (IDS) and bad data detection algorithms in state estimation. These methods either require full knowledge of the system topology and parameters or fail to understand the physical behaviors under attack. WISP (Watching grid Infrastructure Stealthily through Proxies) is designed to provide additional protection to the power grid using only publicly available data. In particular, WISP exploits the spatio-temporal nature of the real time locational marginal prices (LMPs), in conjunction with other information such as bids, weather, outages and load data to analyze anomalous power pricing behaviors and then correlate those observations to localize regions of interest and identify potential cyber events. WISP is non-intrusive as the tool is deployed as a service in the Cloud or on premise and provides reliable information to system operators for enhanced situational awareness, without impeding energy delivery functions. The WISP technology comprises three modules: the data-driven anomaly detection core, the vulnerability and risk analysis and the root cause analysis. The data-driven anomaly detection core performs the tasks of feature selection, anomaly detection and attack region localization. The vulnerability and risk analysis module provides system level information of the vulnerable variables and times, assisting the operators in selecting monitoring and protection nodes. The root cause analysis module takes the detection results and identifies potential operational conditions that contribute to the detected anomalies. In Phase I, we have demonstrated the feasibility and effectiveness of WISP. We developed a realistic electricity market simulator capable of generating normal and attack market data under various operational conditions. We developed a series of cyber-attack detection and analysis algorithms and evaluated them under multiple data sources. Finally, we integrated all modules into an end-to-end software, providing functions for data management, data analytics and visualization. Specifically, we have achieved: (i) real-time data acceptance from external utility interfaces with >99% acceptance rate; (ii) high performance anomaly detection algorithms with >98% detection accuracy and <0.1% false alarm rate; and (iii) ultra-low computing delay <50 milliseconds. Additionally, our team developed algorithms to identify the vulnerable variables in electricity market operations and root cause analysis functions to identify major contributors to the price spikes. These ancillary modules are necessary when deploying WISP in real world industry environment. In Phase II, we have demonstrated the effectiveness of WISP software on realistic largescale power systems. We performed red team testing for the Phase I WISP software and identified software vulnerabilities and implemented corresponding mitigation solutions. We adapted the electricity market simulator for the Texas synthetic 2000-bus system and generated datasets for the false data injection attacks. We created database and visualization interfaces for the Texas system and the ISO New England system. We performed software optimization in terms of operation efficiency, computing speed and detection accuracy. Finally, we tested the software on the Texas system and the ISO New England system and evaluated the detection performance. Overall, we achieved above 89% detection rate, below 3% false alarm rate and below 37 seconds of end-to-end detection delay.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Detection of Spectral Evolution in the Bursts Emitted During the 2008-2009 Active Episode of SGR J1550 - 5418

In early October 2008, the Soft Gamma Repeater SGRJ1550 - 5418 (1E1547.0 - 5408, AXJ155052 - 5418, PSR J1550 - 5418) became active, emitting a series of bursts which triggered the Fermi Gamma-ray Burst Monitor (GBM) after which a second especially intense activity period commenced in 2009 January and a third, less active period was detected in 2009 March-April. Here we analyze the GBM data of all the bursts from the first and last active episodes. We performed temporal and spectral analysis for all events and found that their temporal characteristics are very similar to the ones of other SGR bursts, as well the ones reported for the bursts of the main episode (average burst durations 170ms). In addition, we used our sample of bursts to quantify the systematic uncertainties of the GBM location algorithm for soft gamma-ray transients to less than or equal to 8 degrees. Our spectral analysis indicates significant spectral evolution between the first and last set of events. Although the 2008 October events are best fit with a single blackbody function, for the 2009 bursts an Optically Thin Thermal Bremsstrahlung (OTTB) is clearly preferred. We attribute this evolution to changes in the magnetic field topology of the source, possibly due to effects following the very energetic main bursting episode.

von Kienlin, Andreas↗

A novel invariant mass method to isolate resonance backgrounds from the chiral magnetic effect

The Chiral Magnetic Effect (CME) refers to charge separation along a strong magnetic field, due to topological charge fluctuations in QCD. Charge correlation (Δγ) signals consistent with CME have been first observed almost a decade ago. It has also been known since then that the Δγ is contaminated by a major background from resonance decays coupled with elliptic flow. In this contribution, we propose differential Δγ measurements as function of the pair invariant mass (m inv ). The Δγ in the high minv region is essentially free of resonance backgrounds. In the low m inv region, the Δγ backgrounds show resonance peaks. The CME signal, presumably smooth in m inv , may thus be extracted from a two-component model fit. We demonstrate the feasibility and effectiveness of this novel method by using the AMPT and toy-model Monte-Carlo simulations. We also discuss an application of the method in data analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Knowledge-guided graph machine learning for spatially distributed prediction of daily discharge and nitrogen export dynamics

Spatially distributed prediction of streamflow and nitrogen export dynamics is essential for precision management of agricultural watersheds. While temporal deep learning models such as Long Short-Term Memory (LSTM) have shown strong performance at basin scales, their ability to generalize spatially is limited by insufficient representation of spatial dependencies and flow paths, particularly under data-scarce conditions. To address this gap, we propose HydroGraphNet, a knowledge-guided graph machine learning framework that integrates process-based knowledge and explicit spatial learning into temporal modeling. This framework incorporates directed graph topology to encode watershed connectivity and upstream inflows, with mass balance constraints to improve physical consistency. To enhance generalization in sparsely monitored regions, HydroGraphNet is pretrained on synthetic data generated by the SWAT+ (Soil and Water Assessment Tool Plus) model. We evaluated HydroGraphNet in the Upper Sangamon River Basin (44 HUC-12 subwatersheds, 2001–2020) against two LSTM baselines: a lumped basin-level model and a distributed variant. When benchmarked on SWAT+ simulations in pretraining, HydroGraphNet improved test NSEs by 8.9% (discharge) and 13.7% (NO₃–N load) in temporal extrapolation, and by 27.1% and 34.7% in spatial extrapolation, relative to the Lumped LSTM baseline. After fine-tuning with USGS monitoring data, the model achieved mean test NSE (KGE) scores of 0.768 (0.861) for discharge and 0.626 (0.664) for NO₃–N load, substantially outperforming baselines. Attribution analysis further highlighted the importance of upstream inflow representation and graph-based spatial learning in capturing cross-subwatershed dependencies. The model also reproduced seasonal hydrological and biogeochemical patterns consistent with known processes, demonstrating its robustness and process fidelity for spatially distributed prediction. Altogether, HydroGraphNet advances the integration of physical knowledge and spatially explicit learning in hydrological modeling, offering a generalizable framework for distributed modeling to support spatially targeted water quality management in data-scarce watersheds.

54 ENVIRONMENTAL SCIENCES↗

Understanding Complex Magnetic Spin Textures with Simulation-Assisted Lorentz Transmission Electron Microscopy

There is an increased interest in topologically nontrivial magnetic spin textures such as skyrmions and chiral domain-wall solitons, both from a point of fundamental physics understanding as well as potential technological interest in low-power memory applications. In order to control their behavior, it is necessary to understand their complex spin texture at the nanoscale. Lorentz transmission electron microscopy (LTEM) is a suitable technique for studying these systems due to its high spatial resolution and capability to simultaneously characterize magnetic texture and microstructure. In this work, we present the application of PyLorentz, an open-source software suite that we have developed, for quantitative image analysis of Neel-type skyrmions in thin-film heterostructures. PyLorentz enhances LTEM capabilities by enabling reconstruction of magnetic induction maps from experimental images, as well as simulating LTEM images using micromagnetic simulation data. We demonstrate this for simulated Neel skyrmions as well as experimental data from [Pt/Co/W] multilayer heterostructures. Finally, we also show how simulation-assisted LTEM analysis is crucial for understanding these complex magnetic spin textures, in which the reconstructed magnetic induction map (seen in the LTEM images) differs significantly from the magnetization configuration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental search for the chiral magnetic effect in relativistic heavy-ion collisions: A perspective

The chiral magnetic effect (CME) refers to generation of the electric current along a magnetic field in a chirally imbalanced system of quarks. The latter is predicted by quantum chromodynamics to arise from quark interaction with nontrivial topological fluctuations of the vacuum gluonic field. The CME has been actively searched for in relativistic heavy-ion collisions, where such gluonic field fluctuations and a strong magnetic field are believed to be present. The CME-sensitive observables are unfortunately subject to a possibly large non-CME background, and firm conclusions on a CME observation have not yet been reached. In this perspective, we review the experimental status and progress in the CME search, from the initial measurements more than a decade ago to the dedicated program of isobar collisions in 2018 and the release of the isobar blind analysis result in 2022 to intriguing hints of a possible CME signal in Au + Au collisions, and discuss future prospects of a potential CME discovery in the anticipated high-statistic Au + Au collision data at the Relativistic Heavy-Ion Collider by 2025. We hope such a perspective will help sharpening our focus on the fundamental physics of the CME and steer its experimental search.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Topology-Aware Performance Optimization and Modeling of Adaptive Mesh Refinement Codes for Exascale

Here, we introduce a topology-aware performance optimization and modeling workflow for AMR simulation that includes two new modeling tools, ProgrAMR and Mota Mapper, which interface with the BoxLib AMR framework and the SSTmacro network simulator. ProgrAMR allows us to generate and model the execution of task dependency graphs from high-level specifications of AMR-based applications, which we demonstrate by analyzing two example AMR-based multigrid solvers with varying degrees of asynchrony. Mota Mapper generates multiobjective, network topology-aware box mappings, which we apply to optimize the data layout for the example multigrid solvers. While the sensitivity of these solvers to layout and execution strategy appears to be modest for balanced scenarios, the impact of better mapping algorithms can be significant when performance is highly constrained by network hop latency. Furthermore, we show that network latency in the multigrid bottom solve is the main contributing factor preventing good scaling on exascale-class machines.

97 MATHEMATICS AND COMPUTING↗

TH-CD-207A-03: A Surface Deformation Driven Respiratory Model for Organ Motion Tracking in Lung Cancer Radiotherapy

Purpose: To propose and validate a novel real-time surface-mesh-based internal organ-external surface motion and deformation tracking method for lung cancer radiotherapy. Methods: Deformation vector fields (DVFs) which characterizes the internal and external motion are obtained by registering the internal organ and tumor contours and external surface meshes to a reference phase in the 4D CT images using a recent developed local topology preserved non-rigid point matching algorithm (TOP). A composite matrix is constructed by combing the estimated internal and external DVFs. Principle component analysis (PCA) is then applied on the composite matrix to extract principal motion characteristics and finally yield the respiratory motion model parameters which correlates the internal and external motion and deformation. The accuracy of the respiratory motion model is evaluated using a 4D NURBS-based cardiac-torso (NCAT) synthetic phantom and three lung cancer cases. The center of mass (COM) difference is used to measure the tumor motion tracking accuracy, and the Dice’s coefficient (DC), percent error (PE) and Housdourf’s distance (HD) are used to measure the agreement between the predicted and ground truth tumor shape. Results: The mean COM is 0.84±0.49mm and 0.50±0.47mm for the phantom and patient data respectively. The mean DC, PE and HD are 0.93±0.01, 0.13±0.03 and 1.24±0.34 voxels for the phantom, and 0.91±0.04, 0.17±0.07 and 3.93±2.12 voxels for the three lung cancer patients, respectively. Conclusions: We have proposed and validate a real-time surface-mesh-based organ motion and deformation tracking method with an internal-external motion modeling. The preliminary results conducted on a synthetic 4D NCAT phantom and 4D CT images from three lung cancer cases show that the proposed method is reliable and accurate in tracking both the tumor motion trajectory and deformation, which can serve as a potential tool for real-time organ motion and deformation monitoring in lung cancer radiotherapy. This work is supported in part by grant from VARIAN MEDICAL SYSTEMS INC, the National Natural Science Foundation of China (no 81428019 and no 81301940), the Guangdong Natural Science Foundation (2015A030313302)and the 2015 Pearl River S&T Nova Program of Guangzhou (201506010096).

60 APPLIED LIFE SCIENCES↗