Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heterogeneous computing”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Thermal Performance Analysis of a Geologic Borehole Repository

The Brazilian Nuclear Research Institute (IPEN) proposed a design for the disposal of Disused Sealed Radioactive Sources (DSRS) based on the IAEA Borehole Disposal of Sealed Radioactive Sources (BOSS) design that would allow the entirety of Brazil’s inventory of DSRS to be disposed in a single borehole. The proposed IPEN design allows for 170 waste packages (WPs) containing DSRS (such as Co-60 and Cs-137) to be stacked on top of each other inside the borehole. The primary objective of this work was to evaluate the thermal performance of a conservative approach to the IPEN proposal with the equivalent of two WPs and two different inside configurations using Co-60 as the radioactive heat source. The current WP configuration (heterogeneous) for the IPEN proposal has 60% of the WP volume being occupied by a nuclear radioactive heat source and the remaining 40% as vacant space. The second configuration (homogeneous) considered for this project was a homogeneous case where 100% of the WP volume was occupied by a nuclear radioactive heat source. The computational models for the thermal analyses of the WP configurations with the Co-60 heat source considered three different cooling mechanisms (conduction, radiation, and convection) and the effect of mesh size on the results from the thermal analysis. The results of the analyses yielded maximum temperatures inside the WPs for both of the WP configurations and various mesh sizes. The heterogeneous WP considered the cooling mechanisms of conduction, convection, and radiation. The temperature results from the heterogeneous WP analysis suggest that the model is cooled predominantly by conduction with effect of radiation and natural convection on cooling being negligible. From the thermal analysis comparing the two WP configurations, the results suggest that either WP configuration could be used for the design. The mesh sensitivity results verify the meshes used, and results obtained from the thermal analyses were close to being independent of mesh size. The results from the computational case and analytically-calculated case for the homogeneous WP in benchmarking were almost identical, which indicates that the computational approach used here was successfully verified by the analytical solution.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2011

Equilibrium Speciation of Select Lanthanides in the Presence of Acidic Ligands in Homo- and Heterogeneous Solutions

This dissertation explores lanthanide speciation in liquid solution systems related to separation schemes involving the acidic ligands: bis(2-ethylhexyl) phosphoric acid (HDEHP), lactate, and 8-hydroxyquinoline. Equilibrium speciation of neodymium (Nd 3+ ), sodium (Na+), HDEHP, water, and lactate in the TALSPEAK liquid-liquid extraction system was explored under varied Nd 3+ loading of HDEHP in the organic phase and through extraction from aqueous HCl and lactate media. System speciation was probed through vapor pressure osmometry, visible and Fourier Transform Infrared (FTIR) spectroscopy, 22 Na and 13 C labeled lactate radiotracer distribution measurements, Karl Fischer titrations, and equilibrium pH measurements. Distribution of Nd 3+ , Na + , lactate, and equilibrium pH were modeled using the SXLSQI software to obtain logKNd and logKNa extraction constants under selected conditions. Results showed that high Nd 3+ loading of the HDEHP led to Nd 3+ speciation that departs from the ion exchange mechanism and includes formation of highly aggregated, polynuclear [NdLactate(DEHP) 2 ] x ; (with x > 1). By substituting lanthanum (La 3+ ) for Nd 3+ in this system, NMR scoping experiments using 23 Na, 31 P nuclei and 13 C labeled lactate were performed. Results indicated that this technique is sensitive to changes in system speciation, and that further experiments are warranted. In a homogeneous system representing the TALSPEAK aqueous phase, Lactate protonation behavior at various temperatures was characterized using a combination of potentiometric titration and modeling with the Hyperquad computer program. The temperature dependent deprotonation behavior of lactate showed little change with temperature at 2.0 M NaCl ionic strength. Cloud point extraction is a non-traditional separation technique that starts with a homogeneous phase that becomes heterogeneous by the micellization of surfactants through the increase of temperature. To better understand the behavior of europium (Eu 3+ ) and 8-hydroxyquinoline under cloud point extraction conditions, potentiometric and spectrophotometric titrations coupled with modeling with Hyperquad and SQUAD computer programs were performed to assess europium (Eu 3+ ) and 8-hydroxyquinoline speciation. Experiments in both water and a 1wt% Triton X-114/water mixed solvent were compared to understand the effect of Triton X-114 on the system speciation. Results indicated that increased solvation of 8-hydroxyquinoline by the mixed solvent lead to more stable complexes involving 8-hydroxyquinoline than in water, whereas competition between hydroxide and Triton X-114 for Eu 3+ led to lower stability hydrolysis complexes in the mixed solvent than in water. Lanthanide speciation is challenging due to the trivalent oxidation state that leads to multiple ligand complexes, including some mixed complexes. The complexity of the system demands well-designed and precise experiments that capture the nuances of the chemistry. This work increased the understanding of lanthanide speciation in the explored systems, but more work is required to produce a comprehensive understanding of the speciation involved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Evolution of microscopic heterogeneity and dynamics in choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are an emerging class of non-aqueous solvents that are potentially scalable, easy to prepare and functionalize for many applications ranging from biomass processing to energy storage technologies. Predictive understanding of the fundamental correlations between local structure and macroscopic properties is needed to exploit the large design space and tunability of DESs for specific applications. Here, we employ a range of computational and experimental techniques that span length-scales from molecular to macroscopic and timescales from picoseconds to seconds to study the evolution of structure and dynamics in model DESs, namely Glyceline and Ethaline, starting from the parent compounds. We show that systematic addition of choline chloride leads to microscopic heterogeneities that alter the primary structural relaxation in glycerol and ethylene glycol and result in new dynamic modes that are strongly correlated to the macroscopic properties of the DES formed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Porting DMRG++ Scientific Application to OpenPOWER

With the rapidly changing microprocessor designs and architectural diversity (multi-cores, many-cores, accelerators) for the next generation HPC systems, scientific applications must adapt to the hardware, to exploit the different types of parallelism and resources available in the architecture. To get the benefit of all the in-node hardware threads, it is important to use a single programming model to map and coordinate the available work to the different heterogeneous execution units in the node (e.g., multi-core hardware threads (latency optimized), accelerators (bandwidth optimized), etc.).Our goal is to show that we can manage the node complexity of these systems by using OpenMP for in-node parallelization by exploiting different “programming styles” supported by OpenMP 4.5 to program CPU cores and accelerators. Finding out the suitable programming-style (e.g., SPMD style, multi-level tasks, accelerator programming, nested parallelism, or a combination of these) using the latest features of OpenMP to maximize performance and achieve performance portability across heterogeneous and homogeneous systems is still an open research problem.We developed a mini-application, Kronecker Product (KP), from the original DMRG++ application (sparse matrix algebra) computational motif to experiment with different OpenMP programming styles on an OpenPOWER architecture and present their results in this paper.

Chatterjee, Ronnie↗
DOE OSTI2019

A Real-Time Big Data Control-Theoretical Framework for Cyber-Physical-Human Systems [Book Chapter]

It is known that cyber-physical-human systems naturally arise from interdependent infrastructure systems and smart connected communities. Such applications require ubiquitous information sensing and processing, intelligent machine-to-machine communication for a seamless coordination, as well as intelligent interactions between humans and machines. This chapter presents a control-theoretical framework to model heterogeneous physical dynamic systems, information and communication, as well as cooperative controls and/or distributed optimization of such interconnected systems. It is shown that efficient analytical and computational algorithms can be modularly designed and hierarchically implemented to operate and optimize cyber-physical-human systems, first to quantify individually the input–output relationship of nonlinear dynamic behaviors of every physical subsystem, then to coordinate locally both cyber-physical interactions of neighboring agents as well as physical-human interactions, and finally to dynamically model and optimize the overall networked system. The hierarchical structure makes the overall optimization and control problem scalable and solvable. Moreover, the three levels integrate individual designs and optimization, distributed cooperative optimization, and decision-making through real-time, data-driven, model-based learning and control. Specifically, one of the contributions of the chapter is to demonstrate how the combination of dissipativity theory and cooperative control serves as a natural framework and promising tools to analyze, optimize, and control such large-scale system. Application to digital grid is investigated as an illustrative example.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Catalytic consequences of hydrogen addition events and solvent-adsorbate interactions during guaiacol-H 2 reactions at the H 2 O-Ru(0 0 0 1) interface

Catalytic reactions of biomass-derived phenolics and H 2 occur on transition metal surfaces via competitive C–O cleavage and ring saturation pathways, with both requiring multiple hydrogen addition events before forming their respective rate limiting transition states. These events are markedly affected by solvent chemical identity, with polar protic solvents ionizing hydrogen adatoms (H*) to interfacial protons (H + ) and opening up new catalytic routes. Here, we establish the reaction coordinate space for guaiacol-H 2 reactions on Ru(0 0 0 1) using density functional theory and describe the atomic-scale effect of a polar protic solvent, H 2 O. Coupled H + and H* attack leads to quasi-equilibrated enol and keto intermediates as the precursors for C–O cleavage and ring saturation, respectively. For C–O cleavage, H 2 O solvent enables a lower energy pathway via concomitant transfer of the hydroxyl H + to the methoxy oxygen during C–OCH 3 cleavage, forming a charge separated [Ru(s)–(C 6 H 5 O – )…(H + )…OCH 3 ] transition state and reducing the barrier by up to 0.8 eV as compared to unassisted C–OCH 3 cleavage. For ring saturation, H* attack onto an unsaturated meta carbon is rate limiting with no direct solvent participation, suggesting that protic polar solvents selectively promote the C–O cleavage pathway. Taken together, we show that activating guaiacol for either C–O bond cleavage or ring saturation product formation depends on the reactive hydrogen identity (H* or H + ), enol/keto isomerization equilibrium, and accessibility of the proton assisted Car–OCH 3 cleavage transition state. Here, all such factors are tunable via changes to the solvent or metal identity.

09 BIOMASS FUELS↗
DOE OSTI2021

A Holistic Approach for Elucidating Local Structure, Dynamics, and Speciation in Molten Salts with High Structural Disorder

To examine ion solvation, exchange, and speciation for minority components in molten salts (MS) typically found as corrosion products, we propose a multimodal approach combining extended X-ray absorption fine structure (EXAFS) spectroscopy, optical spectroscopy, ab initio molecular dynamics (AIMD) simulations, and rate theory of ion exchange. Going beyond conventional EXAFS analysis, our method can accurately quantify populations of different coordination states of ions with highly disordered coordination environments via linear combination fitting of the EXAFS spectra of these coordination states computed from AIMD to the experimental EXAFS spectrum. In a case study of dilute Ni(II) dissolved in the ZnCl 2 +KCl melts, our method reveals heterogeneous distributions of coordination states of Ni(II) that are sensitive to variations in temperature and melt composition. These results are fully explained by the difference in the chloride exchange dynamics at varied temperatures and melt compositions. Finally, this insight will enable a better understanding and control of ion solubility and transport in MS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Evaluating microstructure evolution in an SOFC electrode using digital volume correlation

Degradation mechanisms within solid oxide fuel cells (SOFC) during thermal cycling limit operational start-up times and cell lifetime, and must therefore be better understood and mitigated. This work explores such mechanisms using digital volume correlation (DVC) techniques applied to lab-based X-ray tomograms where the microstructural evolution is evaluated during the operational cycling of a Ni–YSZ/YSZ cell. To emulate reduced start-up times, five tomograms were collected over four operational thermal cycles to 750 °C at various ramp-rates: 3, 10, 20 and 30 °C min –1 . Two key features are observed in the 3D strain profiles. Firstly, during low ramp-rates the material produces microscopic channels towards cluster-points. This is thought to be caused by the ceramic skeleton inhibiting the effects of sintering. Secondly, previously unseen macroscopic ‘waves’ developed after high-rate cycling, consisting of linear regions of compression and tension throughout the sample. These wave features decay away from a heterogeneous defect which is thought to be responsible for the non-uniform strain profile. This work demonstrates the first use of sub-micron DVC computations applied to an SOFC exposed to operationally relevant temperatures. Lastly, these findings will assist in the development of new electrode materials from their fabrication to operation, ultimately supporting commercial viability of SOFCs.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Nanoconfinement and mass transport in metal–organic frameworks

The ubiquity of metal–organic frameworks in recent scientific literature underscores their highly versatile nature. MOFs have been developed for use in a wide array of applications, including: sensors, catalysis, separations, drug delivery, and electrochemical processes. Often overlooked in the discussion of MOF-based materials is the mass transport of guest molecules within the pores and channels. Given the wide distribution of pore sizes, linker functionalization, and crystal sizes, molecular diffusion within MOFs can be highly dependent on the MOF-guest system. In this paper, we discuss the major factors that govern the mass transport of molecules through MOFs at both the intracrystalline and intercrystalline scale; provide an overview of the experimental and computational methods used to measure guest diffusivity within MOFs; and highlight the relevance of mass transfer in the applications of MOFs in electrochemical systems, separations, and heterogeneous catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Transformed Fourier and Fick equations for the control of heat and mass diffusion

We review recent advances in the control of diffusion processes in thermodynamics and life sciences through geometric transforms in the Fourier and Fick equations, which govern heat and mass diffusion, respectively. We propose to further encompass transport properties in the transformed equations, whereby the temperature is governed by a three-dimensional, time-dependent, anisotropic heterogeneous convection-diffusion equation, which is a parabolic partial differential equation combining the diffusion equation and the advection equation. We perform two dimensional finite element computations for cloaks, concentrators and rotators of a complex shape in the transient regime. We precise that in contrast to invisibility cloaks for waves, the temperature (or mass concentration) inside a diffusion cloak crucially depends upon time, its distance from the source, and the diffusivity of the invisibility region. However, heat (or mass) diffusion outside cloaks, concentrators and rotators is unaffected by their presence, whatever their shape or position. Finally, we propose simplified designs of layered cylindrical and spherical diffusion cloaks that might foster experimental efforts in thermal and biochemical metamaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Homogenization and Upscaling for Diffusion, Heat Conduction, and Wave Propagation in Heterogeneous Materials

We present a general homogenization method for diffusion, heat conduction, and wave propagation in a periodic heterogeneous material with piecewise constants. The method is relevant to the frequently encountered upscaling issues for heterogeneous materials. The dispersion relation for each problem is first expressed in the general form where the frequency (or wavenumber k) is expanded in terms of the wavenumber k (or frequency ). A general homogenization model can be directly obtained with any given dispersion relation. Next step we study the unit cell of the heterogeneous material and derive the exact dispersion relation. The final homogenized equations include both leading order terms (effective properties) and high order contributions that represent the effect of the microscopic heterogeneity on the macroscopic behavior. That effect can be lumped into a single dimensionless heterogeneity parameter ß, which is bounded between and has a universal expression for all three problems. Numerical examples validate the proposed method and demonstrate a significant computational saving.

Xu, Zhijie↗
DOE OSTI2025

Preliminary Study on Fine-Grained Power and Energy Measurements on Grace Hopper GH200 with Open-Source Performance Tools

The increasing adoption of tightly integrated, heterogeneous architectures, combined with the slowdown of Moore’s law, has made application power and energy-driven optimizations critical to efficiently use high-performance computing systems. This paper introduces a newly developed open-source toolkit that seamlessly integrates the Linux real-time hardware monitoring program hwmon with the Performance Application Programming Interface and the Score-P performance measurement system, thereby enabling fine-grained power and energy measurements for high-performance computing applications. Our primary target platform is the Wombat test bed, which is a system based on the NVIDIA GH200 superchip. The toolkit can capture transient power peaks with high temporal resolution (50 ms) and, thanks to Score-P integration, can map power metrics to specific code regions, thereby providing actionable information on power-intensive operations and inefficiencies. The toolkit also provides a holistic view of both the power and the energy consumption of the entire GH200 superchip by covering all major components: the Grace CPU, the Hopper GPU, and the I/O subsystem. Experiments that use Locally Self-consistent Multiple Scattering, which is an application for first-principles calculations of materials developed at Oak Ridge National Laboratory, have demonstrated the tool’s ability to identify transient power spikes and uncover opportunities for energy-aware optimizations. Additionally, we introduce a Python-based utility for converting Open Trace Format 2 traces to Parquet format, thus enabling advanced data analysis for numerical integration methods applied to power data for accurate energy profiling.

Hernandez Mendoza, Oscar [ORNL] (ORCID:00000002538↗
DOE OSTI2026

Identification of mechanisms driving heterogeneous void growth in ductile aluminum

Void growth plays a central role in ductile fracture, yet the specific mechanisms that control this remain obscure. Classical models, such as those proposed by Rice and Tracey in 1969, are able to capture average rates of void growth, but cannot capture the heterogeneity of individual void growth. Building on recent work, the present study employs laboratory-based diffraction contrast tomography and in-situ x-ray computed tomography to investigate the effect of grain structure and other microstructural factors on void growth in an Al-2219 alloy. Crystal plasticity finite element (CP-FE) modeling is used alongside experimental data to evaluate the contributions of local mechanical states, grain orientation, grain size, and neighboring microstructural features. No strong linear relationships are found with any of the considered descriptors and void growth rate. Potential complex nonlinear relationships are explored with the use of a random forest regression model, which identifies initial void volume, void aspect ratio, local normal stress state, local shear stress state, and local equivalent plastic strain (EQPS) as features that most improve void growth rate predictions. The combination of these analyses suggests that these features should be prioritized to improve models of void growth.

Diffraction contrast tomography (DCT)↗
DOE OSTI2023

Lithium Metal Battery Quality Control via Transformer–CNN Segmentation

Lithium metal battery (LMB) has the potential to be the next-generation battery system because of its high theoretical energy density. However, defects known as dendrites are formed by heterogeneous lithium (Li) plating, which hinders the development and utilization of LMBs. Non-destructive techniques to observe the dendrite morphology often use X-ray computed tomography (XCT) to provide cross-sectional views. To retrieve three-dimensional structures inside a battery, image segmentation becomes essential to quantitatively analyze XCT images. This work proposes a new semantic segmentation approach using a transformer-based neural network called TransforCNN that is capable of segmenting out dendrites from XCT data. In addition, we compare the performance of the proposed TransforCNN with three other algorithms, U-Net, Y-Net, and E-Net, consisting of an ensemble network model for XCT analysis. Our results show the advantages of using TransforCNN when evaluating over-segmentation metrics, such as mean intersection over union (mIoU) and mean Dice similarity coefficient (mDSC), as well as through several qualitatively comparative visualizations.

Quenum, Jerome (ORCID:0000000271265853)↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2021

ExaWorks: Workflows for Exascale

Exascale computers will offer transformative capabilities to combine data-driven and learning-based approaches with traditional simulation applications to accelerate scientific discovery and insight. These software combinations and integrations, however, are difficult to achieve due to challenges of coordination and deployment of heterogeneous software components on diverse and massive platforms. We present the ExaWorks project, which can address many of these challenges: ExaWorks is leading a co-design process to create a workflow Software Development Toolkit (SDK) consisting of a wide range of workflow management tools that can be composed and interoperate through common interfaces. We describe the initial set of tools and interfaces supported by the SDK, efforts to make them easier to apply to complex science challenges, and examples of their application to exemplar cases. Furthermore, we discuss how our project is working with the workflows community, large computing facilities as well as HPC platform vendors to sustainably address the requirements of workflows at the exascale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Speckle contrast diffuse correlation tomography of complex turbid medium flow

Purpose: Developed herein is a three-dimensional (3D) flow contrast imaging system leveraging advancements in the extension of laser speckle contrast imaging theories to deep tissues along with our recently developed finite-element diffuse correlation tomography (DCT) reconstruction scheme. This technique, termed speckle contrast diffuse correlation tomography (scDCT), enables incorporation of complex optical property heterogeneities and sample boundaries. When combined with a reflectance-based design, this system facilitates a rapid segue into flow contrast imaging of larger, in vivo applications such as humans. Methods: A highly sensitive CCD camera was integrated into a reflectance-based optical system. Four long-coherence laser source positions were coupled to an optical switch for sequencing of tomographic data acquisition providing multiple projections through the sample. This system was investigated through incorporation of liquid and solid tissue-like phantoms exhibiting optical properties and flow characteristics typical of human tissues. Computer simulations were also performed for comparisons. A uniquely encountered smear correction algorithm was employed to correct point-source illumination contributions during image capture with the frame-transfer CCD and reflectance setup. Results: Measurements with scDCT on a homogeneous liquid phantom showed that speckle contrast-based deep flow indices were within 12% of those from standard DCT. Inclusion of a solid phantom submerged below the liquid phantom surface allowed for heterogeneity detection and validation. The heterogeneity was identified successfully by reconstructed 3D flow contrast tomography with scDCT. The heterogeneity center and dimensions and averaged relative flow (within 3%) and localization were in agreement with actuality and computer simulations, respectively. Conclusions: A custom cost-effective CCD-based reflectance 3D flow imaging system demonstrated rapid acquisition of dense boundary data and, with further studies, a high potential for translatability to real tissues with arbitrary boundaries. A requisite correction was also found for measurements in the fashion of scDCT to recover accurate speckle contrast of deep tissues.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2010

Deformable image registration of heterogeneous human lung incorporating the bronchial tree

Purpose: To investigate the effect of the bronchial tree on the accuracy of biomechanical-based deformable image registration of human lungs. Methods: Three dimensional finite element models have been developed using four dimensional computed tomography image data of ten lung cancer patients. Each model is built of a body, left and right lungs, tumor, and bronchial trees. Triangular shell elements are used for the bronchial trees while tetrahedral elements are used for other components. Hyperelastic material properties based on experimental investigation on human lungs are used for the lung parenchyma. Different material properties are assigned for the bronchial tree using five values for the modulus of elasticity of 0.01, 0.12, 0.5, 10, and 18 MPa. Lungs are modeled to slide inside chest cavities by applying frictionless contact surfaces between each lung and corresponding chest cavity. The accuracy of the models is examined using an average of 40 bronchial bifurcation points identified on inhale and exhale images. Relative accuracy is evaluated by comparing the displacement of all nodes within the lungs as well as the dosimetric difference at the exhale position predicted by the model. Results: There is no significant effect of bronchial tree on the model accuracy based on the bifurcation points analysis. However, on the local level, using an average of 38 000 nodes, there is a maximum difference of 8.5 mm in the deformation of the bronchial trees, as the modulus of elasticity of the bronchial trees increases from 0.01 to 18 MPa; however, more than 96% of nodes are within a 2.5 mm difference in each direction. The average dose difference at the predicted exhale position is less than 35 cGy between the models. Conclusions: The bronchial tree has little effect on the global deformation and the accuracy of deformable image registration of lungs. Hence, the homogenous model is a reasonable assumption. Since there are some local deformation differences between nodes as the material properties of the bronchial tree change that may affect the accuracy of dosimetric results, heterogeneity may be required for a smaller scale modeling of lungs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗