Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “map improvement”

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 739 records · Page 41

Revisiting the Supernova Engines in the 3C 397 and W49B Supernova Remnants

The nature of the supernova remnants (SNRs) 3C 397 and W49B has long been a subject of debate, with prior studies offering conflicting interpretations between thermonuclear and core-collapse scenarios. To help settle this debate, we present a systematic, spatially resolved, spectroscopic analysis of both remnants using XMM-Newton. By applying multicomponent thermal models, we derive key physical properties including elemental abundances, ejecta temperatures, ambient densities, and explosion energetics. We compare the inferred metal abundance ratios to a wide range of core-collapse and thermonuclear nucleosynthesis models, including new models whose explosion energies differ from the canonical value of 10 51 erg. We find that the observed Fe/Si and Ca/Si ratios in both SNRs are best matched by certain thermonuclear models. However, no model fully reproduces the complete set of observed abundance patterns. In 3C 397, high Fe enrichment and spatial abundance variations suggest interaction with a dense progenitor environment, and W49B’s composition is overall consistent with a thermonuclear origin; however, both require a low-energy (∼10 50 erg) supernova explosion. We additionally map the Fe Kα line centroid energies and find a spread, with W49B falling within the core-collapse region—highlighting both environmental complexity and the limitations of this diagnostic for supernova classification. Our results highlight the need for caution in relying on any single diagnostic or nucleosynthesis model for supernova typing, underscore the need for improved nucleosynthesis models, and motivate future high-resolution, high-throughput observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Ducted Fuel Injection Provides Consistently Lower Soot Emissions in Sweep to Full-Load Conditions

Earlier studies have proven how ducted fuel injection (DFI) substantially reduces soot for low- and mid-load conditions in heavy-duty engines, without significant adverse effects on other emissions. Nevertheless, no comprehensive DFI study exists showing soot reductions at high- and full-load conditions. This study investigated DFI in a single-cylinder, 1.7-L, optical engine from low- to full-load conditions with a low-net-carbon fuel consisting of 80% renewable diesel and 20% biodiesel. Over the tested load range, DFI reduced engine-out soot by 38.1–63.1% compared to conventional diesel combustion (CDC). This soot reduction occurred without significant detrimental effects on other emission types. Thus, DFI reduced the severity of the soot–NO x tradeoff at all tested conditions. While DFI delivered considerable soot reductions in thepresent study, previous DFI studies at low- and mid-load conditions delivered larger soot reductions (>90%) compared to CDC operation at the same conditions. Therefore, the DFI configuration used here has been deemed nonoptimal (in terms of parameters such as the injector-spray and piston geometries), and several improvements are recommended for future studies with high-load DFI. Further, these improvements include employing better spray-duct alignment, a deeper piston bowl with a smaller injector umbrella angle, and a fuel injector that opens and closes faster. The study also suggests future research to make DFI ready for commercialization, such as metal-engine tests to ensure desirable DFI performance over an engine’s complete speed/load map. Overall, this study supports the continued development and commercialization of DFI to meet upcoming emissions regulations for heavy-duty vehicles. Specifically, multicylinder engine experiments and CFD simulations should be utilized to optimize the performance and clarify the full potential of DFI.

33 ADVANCED PROPULSION SYSTEMS↗

System and method for focal-plane illuminator/detector (FASID) design for improved graded index lenses

The present disclosure relates to a detector system for imaging an optical signal received by a graded index (GRIN) optical element to account for known variations in a graded index distribution of the GRIN optical element. The detector system incorporates a plurality of optical detector elements configured to receive optical rays received by the GRIN optical element at specific locations on a plane of the GRIN optical element. Ray tracing software is configured to receive and map the optical rays to a plurality of additional specific locations on the plane based on the known variations in the graded index distribution of the GRIN optical element. A processor uses algorithms for diagonalization of a linear system matrix to determine information on both an intensity and an angle of the received optical rays at each one of the plurality of specific locations on the plane.

Bude, Jeffrey Devin↗

Data for Genome-Wide Association and Genomic Prediction for Yield and Component Traits of Miscanthus sacchariflorus

Accelerating biomass improvement is a major goal of miscanthus breeding. The development and implementation of genomic-enabled breeding tools, like marker-assisted selection (MAS) and genomic selection, has the potential to improve the efficiency of miscanthus breeding. The present study conducted genome-wide association (GWA) and genomic prediction of biomass yield and 14 yield-components traits in Miscanthus sacchariflorus . We evaluated a diversity panel with 590 accessions of M. sacchariflorus grown across four years in one subtropical and three temperate locations and genotyped with 268,109 single-nucleotide polymorphisms (SNPs). The GWA study identified a total of 835 significant SNPs and 674 candidate genes across all traits and locations. Of the significant SNPs identified, 280 were localized in mapped quantitative trait loci intervals and proximal to SNPs identified for similar traits in previously reported miscanthus studies, providing additional support for the importance of these genomic regions for biomass yield. Our study gave insights into the genetic basis for yield-component traits in M. sacchariflorus that may facilitate marker-assisted breeding for biomass yield. Genomic prediction accuracy for the yield-related traits ranged from 0.15 to 0.52 across all locations and genetic groups. Prediction accuracies within the six genetic groupings of M. sacchariflorus were limited due to low sample sizes. Nevertheless, the Korea/NE China/Russia (N = 237) genetic group had the highest prediction accuracy of all genetic groups (ranging 0.26–0.71), suggesting that with adequate sample sizes, there is strong potential for genomic selection within the genetic groupings of M. sacchariflorus . This study indicated that MAS and genomic prediction will likely be beneficial for conducting population-improvement of M. sacchariflorus .

Biomass Analytics↗

Numerical Investigation of Two-Stage Vapor Compression System with Simultaneous Vapor and Liquid Injection

Two-stage vapor compression systems can be advantageous over single-stage systems by providing improved system performance, lower discharge temperature and reduced throttling losses. These systems employ various intermediate configurations, such as liquid injection and vapor injection. This paper presents a configuration for a R1234ze(E) two-stage air-conditioning system using a turbo compressor with two injection ports, one for vapor injection and the second for liquid injection. The liquid injection is used to cool the motor and electronics. A component-based representation and solution approach was used to simulate the two-stage compression system with simultaneous vapor and liquid injection at steady state. The turbo compressor was represented using a customized performance map. The condenser and the evaporator were modeled using finite-volume approach. A parametric study was conducted to assess the impact of the following three variables on the system performance: vapor injection ratio, condenser air flow rate, and discharge pipe pressure drop. The simulation results show that as the vapor injection ratio increased, the system performance undergoes a tradeoff between an enhanced subcooling effect (and thus enhanced unit refrigeration capacity) and a decreased suction mass flow rate. Maximum COP occurs when the vapor injection ratio was 0.1. The results also show that as the condenser air flow rate increased, both the capacity and power consumption (including fan power) increased monotonically, and COP increased first and then decreased. At 75% load, the COP improvement at the optimum flow rate was marginally less than 0.5%. Lastly, higher discharge pipe pressure drop increased the discharge pressure. It showed very small effect on the overall system performance at the condition selected for the current study.

component-based steady-state simulation↗

The role of the Pt-group dehydrogenation catalyst in alkane metathesis for polyolefin deconstruction

Recent proposed approaches in the depolymerization of waste plastics employ an olefin intermediate to produce alkanes or alkenes using olefin metathesis in tandem chemistry. Here, in this study, we investigated the role of the dehydrogenation catalyst on reaction rate, kinetics, and product distribution in heterogeneous tandem dehydrogenation and olefin metathesis (alkane metathesis) of three different alkane reactants, including polyethylene. We found that many properties to which alkane dehydrogenation rates were sensitive-including metal composition, nanoparticle size, and surface doping of Re species also controlled activity in Tandem D/OM. When comparing Pd, Pt, and Pt 3 Sn 1 , supported Pd in tandem with a Re 2 O 7 olefin metathesis catalyst showed four-fold higher activity (surface area basis) compared to Pt or Pt 3 Sn 1 catalysts on the same support, mainly due to differences in the rate of hydrogenation. Catalyst preparation resulted in metal nanoparticles partially covered by ReO x , as seen from elemental mapping. Co-location of Re 2 O 7 and Pd correlated with increased rates of hydrogenation (i.e., an increase in the rate of alkane formation and simultaneous lowering of the rate of alkene formation), with a reaction order in catalyst study that further supported this conclusion. The Pd and Re 2 O 7 system displayed marked improvement compared to Pt or Pt 3 Sn 1 with Re 2 O 7 , and previous work, in the depolymerization rate of a linear polyethylene feedstock, with over 94 % reduction in polymer molecular weight in 15 h at 190 °C using less catalyst and increased reactant loadings, while keeping solvent to polymer consumption below 2.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plan evaluation for heat resilience: complementary methods to comprehensively assess heat planning in Tempe and Tucson, Arizona

Abstract Escalating impacts from climate change and urban heat are increasing the urgency for communities to equitably plan for heat resilience. Cities in the desert Southwest are among the hottest and fastest warming in the U.S., placing them on the front lines of heat planning. Urban heat resilience requires an integrated planning approach that coordinates strategies across the network of plans that shape the built environment and risk patterns. To date, few studies have assessed cities’ progress on heat planning. This research is the first to combine two emerging plan evaluation approaches to examine how networks of plans shape urban heat resilience through case studies of Tempe and Tucson, Arizona. The first methodology, Plan Quality Evaluation for Heat Resilience, adapts existing plan quality assessment approaches to heat. We assess whether plans meet 56 criteria across seven principles of high-quality planning and the types of heat strategies included in the plans. The second methodology, the Plan Integration for Resilience Scorecard™ (PIRS™) for Heat, focuses on plan policies that could influence urban heat hazards. We categorize policies by policy tool and heat mitigation strategy and score them based on their heat impact. Scored policies are then mapped to evaluate their spatial distribution and the net effect of the plan network. The resulting PIRS™ for Heat scorecard is compared with heat vulnerability indicators to assess policy alignment with risks. We find that both cities are proactively planning for heat resilience using similar plan and strategy types, however, there are clear and consistent opportunities for improvement. Combining these complementary plan evaluation methods provides a more comprehensive understanding of how plans address heat and a generalizable approach that communities everywhere could use to identify opportunities for improved heat resilience planning.

Environmental Sciences & Ecology↗

A smart alarm for particle accelerator beamline operations

We present the initial results of a proof-of-concept 'smart alarm' for the Continuous Electron Beam Accelerator Facility injector beamline at Jefferson Lab. To minimize machine downtime and improve operational efficiency, an autonomous alarm system able to identify and diagnose unusual machine states is needed. Our approach leverages a trained neural network capable of alerting operators (a) when an anomalous condition exists in the beamline and (b) identifying the element setting that is the root cause. The tool is based on an inverse model that maps beamline readings (diagnostic readbacks) to settings (beamline attributes operators can modify). The model takes as input readings from the machine and computes machine settings which are compared to control setpoints. Instances where predictions differ from setpoints by a user-defined threshold are flagged as anomalous. Given data corresponding to 354 anomalous injector configurations, the model can narrow the root cause of an anomalous condition to three potential candidates with 94.6% accuracy. Furthermore, compared to the current method of identifying anomalous conditions which raises an alarm when machine parameters drift outside their normal tolerances, the data-driven model can identify 83% more anomalous conditions.

43 PARTICLE ACCELERATORS↗

Using automated machine learning for the upscaling of gross primary productivity

Estimating gross primary productivity (GPP) over space and time is fundamental for understanding the response of the terrestrial biosphere to climate change. Eddy covariance flux towers provide in situ estimates of GPP at the ecosystem scale, but their sparse geographical distribution limits larger-scale inference. Machine learning (ML) techniques have been used to address this problem by extrapolating local GPP measurements over space using satellite remote sensing data. However, the accuracy of the regression model can be affected by uncertainties introduced by model selection, parameterization, and choice of explanatory features, among others. Recent advances in automated ML (AutoML) provide a novel automated way to select and synthesize different ML models. In this work, we explore the potential of AutoML by training three major AutoML frameworks on eddy covariance measurements of GPP at 243 globally distributed sites. We compared their ability to predict GPP and its spatial and temporal variability based on different sets of remote sensing explanatory variables. Explanatory variables from only Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data and photosynthetically active radiation explained over 70 % of the monthly variability in GPP, while satellite-derived proxies for canopy structure, photosynthetic activity, environmental stressors, and meteorological variables from reanalysis (ERA5-Land) further improved the frameworks' predictive ability. We found that the AutoML framework Auto-sklearn consistently outperformed other AutoML frameworks as well as a classical random forest regressor in predicting GPP but with small performance differences, reaching an r 2 of up to 0.75. We deployed the best-performing framework to generate global wall-to-wall maps highlighting GPP patterns in good agreement with satellite-derived reference data. This research benchmarks the application of AutoML in GPP estimation and assesses its potential and limitations in quantifying global photosynthetic activity.

54 ENVIRONMENTAL SCIENCES↗

Stochastic reconstruction of fracture network pattern using spatial point processes

Fracture spatial patterns can strongly affect fluid flow in the subsurface. Proximity and distribution of fractures control reservoir flow behavior over various length scales. In many studies, however, simplified geometrical patterns are generated for fractures which may lead to unrealistic subsurface models. Here we introduce a new method for characterization and modeling of fracture spatial patterns based on outcrop observations. We use Ripley's K-function to characterize the arrangement of fracture barycenters and intersection points over various length scales. In addition, we use semivariograms to quantify spatial correlation in fracture intensity maps. Using this information, we develop a stochastic algorithm that generates two-dimensional fracture network realizations with spatial properties similar to those of a real fracture network measured in the field. Numerical simulation models indicate that the generated fracture realizations exhibit similar flow behaviors as that of the original fracture network. Furthermore, such modeling tools expand and improve our capability in building representative fracture models and in quantification of uncertainty in naturally and hydraulically fractured reservoirs.

58 GEOSCIENCES↗

Metall: A persistent memory allocator for data-centric analytics

Data analytics applications transform raw input data into analytics-specific data structures before performing analytics. Unfortunately, such data ingestion steps are often more expensive than analytics. In addition, various types of NVRAM devices are already used in many HPC systems today. Such devices will be useful for storing and reusing data structures beyond a single process life cycle. We developed Metall, a persistent memory allocator built on top of the memory-mapped file mechanism. Metall enables applications to transparently allocate custom C++ data structures into various types of persistent memories. Metall incorporates a concise and high-performance memory management algorithm inspired by Supermalloc and the rich C++ interface developed by Boost.Interprocess library. On a dynamic graph construction workload, Metall achieved up to 11.7x and 48.3x performance improvements over Boost.Interprocess and memkind (PMEM kind), respectively. We also demonstrate Metall’s high adaptability by integrating Metall into a graph processing framework, GraphBLAS Template Library. Here this study’s outcomes indicate that Metall will be a strong tool for accelerating future large-scale data analytics by allowing applications to leverage persistent memory efficiently.

97 MATHEMATICS AND COMPUTING↗

Relation between sampling, sensitivity and precision in strain mapping using the Geometric Phase Analysis method in Scanning Transmission Electron Microscopy

The sensitivity and the precision of the Geometric Phase Analysis (GPA) method for strain characterization is a topic widely discussed in the literature and is usually difficult to quantify. Indeed, the GPA precision is intricately linked to the resolution of the strain maps defined when masking the periodic reflections in Fourier space. In this study an additional parameter, sampling, is proposed to be analyzed regarding the precision of GPA by developing the concept of a phase noise in the GPA equations. Both experimentally and theoretically, the following article demonstrates how the precision, and the sensitivity of the GPA method is improved when using a larger pixel spacing to record an electron micrograph in Scanning Transmission Electron Microscopy (STEM). In conclusion, the counterintuitive concept of increasing the field of view to improve the GPA precision results in an extension of the application of strain characterization methods in STEM towards low deformation levels.

47 OTHER INSTRUMENTATION↗

On-the-fly closed-loop materials discovery via Bayesian active learning

Active learning—the field of machine learning (ML) dedicated to optimal experiment design—has played a part in science as far back as the 18th century when Laplace used it to guide his discovery of celestial mechanics. In this work, we focus a closed-loop, active learning-driven autonomous system on another major challenge, the discovery of advanced materials against the exceedingly complex synthesis-processes-structure-property landscape. We demonstrate an autonomous materials discovery methodology for functional inorganic compounds which allow scientists to fail smarter, learn faster, and spend less resources in their studies, while simultaneously improving trust in scientific results and machine learning tools. This robot science enables science-over-the-network, reducing the economic impact of scientists being physically separated from their labs. The real-time closed-loop, autonomous system for materials exploration and optimization (CAMEO) is implemented at the synchrotron beamline to accelerate the interconnected tasks of phase mapping and property optimization, with each cycle taking seconds to minutes. We also demonstrate an embodiment of human-machine interaction, where human-in-the-loop is called to play a contributing role within each cycle. This work has resulted in the discovery of a novel epitaxial nanocomposite phase-change memory material.

36 MATERIALS SCIENCE↗

Joint Focal Mechanism Inversion Using Downhole and Surface Monitoring at the Decatur, Illinois, CO2 Injection Site

ABSTRACT The three-year CO2 injection period at the Illinois Basin - Decatur Project site (Decatur, Illinois, United States) produced a number of microseismic events distributed in very distinct spatiotemporal clusters with different orientations. Further characterization of the microseismicity encompasses the determination of the event source mechanisms. Initially, the microseismic monitoring network consisted solely of borehole sensors, but has been extended with surface sensors, thereby significantly improving the data coverage over the focal sphere. This article focuses on 23 events from the northernmost microseismic cluster (about 2 km from the injection point) and takes advantage of both, surface and downhole, recordings. The resulting strike-slip east–west-oriented focal planes are all consistent with the east–west orientation of the cluster in map view. The injection-related increase of pore pressure is far below the formation fracture pressure; however, small stress-field changes associated with the pore-pressure increase may reach as far as to the investigated cluster location. Monte Carlo modeling of the slip reactivation potential within this cluster showed that the observed maximum stress-field orientation of N068° is the optimum orientation for fault reactivation of the east–west-oriented cluster. Our results suggest that the east–west orientation of the investigated cluster is the main reason for its activation, even though the cluster is about 2 km away from the low-pressure injection point.

Geochemistry & Geophysics↗

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved Mobilome Delineation in Fragmented Genomes

The mobilome of a microbe, i.e., its set of mobile elements, has major effects on its ecology, and is important to delineate properly in each genome. This becomes more challenging for incomplete genomes, and even more so for metagenome-assembled genomes (MAGs), where misbinning of scaffolds and other losses can occur. Genomic islands (GIs), which integrate into the host chromosome, are a major component of the mobilome. Our GI-detection software TIGER, unique in its precise mapping of GI termini, was applied to 74,561 genomes from 2,473 microbial species, each species containing at least one MAG and one isolate genome. A species-normalized deficit of ∼1.6 GIs/genome was measured for MAGs relative to isolates. To test whether this undercount was due to the higher fragmentation of MAG genomes, TIGER was updated to enable detection of split GIs whose termini are on separate scaffolds or that wrap around the origin of a circular replicon. This doubled GI yields, and the new split GIs matched the quality of single-scaffold GIs, except that highly fragmented GIs may lack central portions. Cross-scaffold search is an important upgrade to GI detection as fragmented genomes increasingly dominate public databases. TIGER2 better captures MAG microdiversity, recovering niche-defining GIs and supporting microbiome research aims such as virus-host linking and ecological assessment.

54 ENVIRONMENTAL SCIENCES↗

Drell-Yan q T resummation of fiducial power corrections at N 3 LL

We consider Drell-Yan production pp → V*X → LX at small q T << Q, where q T and Q are the total transverse momentum and invariant mass of the leptonic final state L. Experimental measurements require fiducial cuts on L, which in general introduce enhanced, linear power corrections in q T /Q. We show that they can be unambiguously predicted from factorization, and resummed to the same order as the leading-power contribution. For the fiducial q T spectrum, they constitute the complete linear power corrections. We thus obtain predictions for the fiducial q T spectrum to N 3 LL and next-to-leading-power in q T /Q. Matching to full NNLO ($α$$^{2}_{s}$), we find that the linear power corrections are indeed the dominant ones, and once included by factorization, the remaining fixed-order corrections become almost negligible below q T ≲ 40 GeV. We also discuss the implications for more complicated observables, and provide predictions for the fiducial Φ* spectrum at N 3 LL+NNLO. We find excellent agreement with ATLAS and CMS measurements of q T and Φ*. We also consider the $p$$^{ℓ}_{T}$ spectrum. We show that it develops leptonic power corrections in q T /(Q – 2$p$$^{ℓ}_{T}$), which diverge near the Jacobian peak $p$$^{ℓ}_{T}$ ~ Q/2 and must be kept to all powers to obtain a meaningful result there. Doing so, we obtain for the first time an analytically resummed result for the $p$$^{ℓ}_{T}$ spectrum around the Jacobian peak at N 3 LL+NNLO. Our method is based on performing a complete tensor decomposition for hadronic and leptonic tensors. We show that in practice this is equivalent to often-used recoil prescriptions, for which our results now provide rigorous, formal justification. Our tensor decomposition yields nine Lorentz-scalar hadronic structure functions, which for Z/γ* → ℓℓ or W → ℓν directly map onto the commonly used angular coefficients, but also holds for arbitrary leptonic final states. In particular, for suitably defined Born-projected leptons it still yields a LO-like angular decomposition even when including QED final-state radiation. Finally, we also discuss the application to q T subtractions. Including the unambiguously predicted fiducial power corrections significantly improves their performance, and in particular makes them applicable near kinematic edges where they otherwise break down due to large leptonic power corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantifying the unknown impact of segmentation uncertainty on image-based simulations

Abstract Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.

42 ENGINEERING↗