Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel methods”

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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Ferroelectricity in high-density H{sub 2}O ice

The origin of longstanding anomalies in experimental studies of the dense solid phases of H{sub 2}O ices VII, VIII, and X is examined using a combination of first-principles theoretical methods. We find that a ferroelectric variant of ice VIII is energetically competitive with the established antiferroelectric form under pressure. The existence of domains of the ferroelectric form within anti-ferroelectric ice can explain previously observed splittings in x-ray diffraction data. The ferroelectric form is stabilized by density and is accompanied by the onset of spontaneous polarization. The presence of local electric fields triggers the preferential parallel orientation of the water molecules in the structure, which could be stabilized in bulk using new high-pressure techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Accelerated fast iterative shrinkage thresholding algorithms for sparsity-regularized cone-beam CT image reconstruction

Purpose: The development of iterative image reconstruction algorithms for cone-beam computed tomography (CBCT) remains an active and important research area. Even with hardware acceleration, the overwhelming majority of the available 3D iterative algorithms that implement nonsmooth regularizers remain computationally burdensome and have not been translated for routine use in time-sensitive applications such as image-guided radiation therapy (IGRT). In this work, two variants of the fast iterative shrinkage thresholding algorithm (FISTA) are proposed and investigated for accelerated iterative image reconstruction in CBCT. Methods: Algorithm acceleration was achieved by replacing the original gradient-descent step in the FISTAs by a subproblem that is solved by use of the ordered subset simultaneous algebraic reconstruction technique (OS-SART). Due to the preconditioning matrix adopted in the OS-SART method, two new weighted proximal problems were introduced and corresponding fast gradient projection-type algorithms were developed for solving them. We also provided efficient numerical implementations of the proposed algorithms that exploit the massive data parallelism of multiple graphics processing units. Results: The improved rates of convergence of the proposed algorithms were quantified in computer-simulation studies and by use of clinical projection data corresponding to an IGRT study. The accelerated FISTAs were shown to possess dramatically improved convergence properties as compared to the standard FISTAs. For example, the number of iterations to achieve a specified reconstruction error could be reduced by an order of magnitude. Volumetric images reconstructed from clinical data were produced in under 4 min. Conclusions: The FISTA achieves a quadratic convergence rate and can therefore potentially reduce the number of iterations required to produce an image of a specified image quality as compared to first-order methods. We have proposed and investigated accelerated FISTAs for use with two nonsmooth penalty functions that will lead to further reductions in image reconstruction times while preserving image quality. Moreover, with the help of a mixed sparsity-regularization, better preservation of soft-tissue structures can be potentially obtained. The algorithms were systematically evaluated by use of computer-simulated and clinical data sets.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Monte Carlo Second- and Third-Order Many-Body Green’s Function Methods with Frequency-Dependent, Nondiagonal Self-Energy

We fully develop the Monte Carlo many-body Green’s function (MC-GF) method with the following enhancements: (1) The truncation order of the perturbation expansion of the Dyson self-energy is raised from the second order (MC-GF2) to the third order (MC-GF3) with the aid of a computerized procedure to enumerate and transform all 84 third-order Goldstone diagrams into Monte Carlo integrable expressions and then into central processing unit (CPU)/graphical processing unit (GPU)-parallel computer codes. (2) An efficient algorithm is proposed that computes all off-diagonal and diagonal elements of the MC-GF2 and MC-GF3 self-energy matrices by common subexpression elimination. (3) The frequency-independent approximation is lifted by introducing a method that computes frequency derivatives of the MC-GF2 and MC-GF3 self-energies up to any arbitrarily high order at nearly no additional computational cost. (4) The imaginary-time integration in the Laplace-transformed expressions of the self-energy is carried out stochastically (instead of using a quadrature in the previous implementations), resulting in a 50- to 200-fold speedup. (5) The efficiency of the redundant-walker convergence acceleration scheme is analyzed numerically, and the guidelines are established to select an optimal number of walkers for maximal efficiency. When such an optimal number is used, the cost per sample is constant of molecular size on either many CPUs or many GPUs. (6) The computational cost to obtain a binding energy within a given statistical uncertainty is observed to increase as (tentatively) O(n 4 ) and O(n 5 ) of molecular size (n is the number of electrons) with and without the diagonal approximation, respectively, as opposed to O(n 5 ) and O(n 6 ) of the corresponding deterministic algorithms. With this method applied to the electron binding energies of C 60 , we show that the third-order corrections to the self-energies are much greater in electron binding energies than in ground-state energies. As a result, they display a sign of oscillatory convergence toward experimental results, not necessarily improving the agreement with increasing perturbation order, justifying MC-GF3 and motivating even higher-order methods.

Algorithms↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Comparison of methods for the measurement of radiation dose distributions in high dose rate (HDR) brachytherapy: Ge-doped optical fiber, EBT3 Gafchromic film, and PRESAGE{sup Registered-Sign} radiochromic plastic

Purpose: Dose distribution measurement in clinical high dose rate (HDR) brachytherapy is challenging, because of the high dose gradients, large dose variations, and small scale, but it is essential to verify accurate treatment planning and treatment equipment performance. The authors compare and evaluate three dosimetry systems for potential use in brachytherapy dose distribution measurement: Ge-doped optical fibers, EBT3 Gafchromic film with multichannel analysis, and the radiochromic material PRESAGE{sup Registered-Sign} with optical-CT readout. Methods: Ge-doped SiO{sub 2} fibers with 6 {mu}m active core and 5.0 mm length were sensitivity-batched and their thermoluminescent properties used via conventional heating and annealing cycles. EBT3 Gafchromic film of 30 {mu}m active thickness was calibrated in three color channels using a nominal 6 MV linear accelerator. A 48-bit transmission scanner and advanced multichannel analysis method were utilized to derive dose measurements. Samples of the solid radiochromic polymer PRESAGE{sup Registered-Sign }, 60 mm diameter and 100 mm height, were analyzed with a parallel beam optical CT scanner. Each dosimetry system was used to measure the dose as a function of radial distance from a Co-60 HDR source, with results compared to Monte Carlo TG-43 model data. Each system was then used to measure the dose distribution along one or more lines through typical clinical dose distributions for cervix brachytherapy, with results compared to treatment planning system (TPS) calculations. Purpose-designed test objects constructed of Solid Water and held within a full-scatter water tank were utilized. Results: All three dosimetry systems reproduced the general shape of the isolated source radial dose function and the TPS dose distribution. However, the dynamic range of EBT3 exceeded those of doped optical fibers and PRESAGE{sup Registered-Sign }, and the latter two suffered from unacceptable noise and artifact. For the experimental conditions used in this study, the useful range from an isolated HDR source was 5-40 mm for fibers, 3-50 mm for EBT3, and 4-21 mm for PRESAGE{sup Registered-Sign }. Fibers demonstrated some over-response at very low dose levels, suffered from volume averaging effects in the dose distribution measurement, and exhibited up to 9% repeatability variation over three repeated measurements. EBT3 demonstrated excellent agreement with Monte Carlo and TPS dose distributions, with up to 3% repeatability over three measurements. PRESAGE{sup Registered-Sign} gave promising results, being the only true 3D dosimeter, but artifacts and noise were apparent. Conclusions: The comparative response of three emerging dosimetry systems for clinical brachytherapy dose distribution measurement has been investigated. Ge-doped optical fibers have excellent spatial resolution for single-direction measurement but are currently too large for complex dose distribution assessment. The use of PRESAGE{sup Registered-Sign} with optical-CT readout gave promising results in the measurement of true 3D dose distributions but further development work is required to reduce noise and improve dynamic range for brachytherapy dose distribution measurements. EBT3 Gafchromic film with multichannel analysis demonstrated accurate and reproducible measurement of dose distributions in HDR brachytherapy. Calibrated dose measurements were possible with agreement within 1.5% of TPS dose calculations. The suitability of EBT3 as a dosimeter for 2D quality control or commissioning work has been demonstrated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Study of Electromagnetic Electron Cyclotron Waves for Kappa Distribution with AC Field in the Magnetosphere of Saturn

Parallel propagating electromagnetic electron cyclotron (EMEC) waves in the extended plasma sheet (~12R{sub S}) and in the outer magnetosphere (~18R{sub S}) of Saturn have been studied. A dispersion relation for parallel propagating relativistic EMEC waves has been applied to the magnetosphere of Saturn, and comparisons have been made with the data of Voyager 1 at these radial distances. The detailed investigations for EMEC waves have been done in the presence of the perpendicular AC electric field, using the kappa distribution function. The relativistic temporal growth rate is calculated by the method of characteristic solution with the data provided by Voyager 1. The effect of the suprathermal electron density, temperature anisotropy, frequency of AC electric field, thermal energy of ions, and relativistic factor on the temporal growth rate of EMEC wave emission has been studied. The simulation results show that the growth of parallel propagating EMEC waves is significantly affected by variations in the temperature anisotropy, electron density, ion thermal energy, and relativistic factor in both the extended plasma sheet and the outer magnetosphere of Saturn. The temperature anisotropy (T⊥/T{sub ║}), ion thermal energy (K{sub B}T{sub ║i}), and electron density (n{sub 0}) have been found to be a major source of free energy for parallel propagating EMEC waves in both regions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Bayesian Spectral Moment Estimation and Uncertainty Quantification

We present a Bayesian spectral fitting method developed for spectroscopic data analysis, particularly (but not solely) in the context of fusion energy research. The presented techniques are particularly valuable to estimating moments and corresponding uncertainties whenever the spectra result from line-integrated measurements in nonuniform plasmas, for which the approximation of atomic line shapes being ideal Gaussians gives poor estimates. We decompose multiple, potentially overlapping spectral lines into a sum of Gauss-Hermite polynomials, whose properties allow efficient truncation and uncertainty quantification, often with only three free parameters per atomic emission line. Tests with both synthetic and experimental data demonstrate the effectiveness and robustness where more standard nonlinear fitting routines may experience difficulties. A parallelized version of our implementation is publicly released under an open source license.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Parallel Radioisotope Collection and Analysis in Response to the Fukushima Release

Two independent air samplers were operated at Pacific Northwest National Laboratory in parallel during the collection of samples from the Fukushima reactor releases. One system is an automated aerosol collection and analysis unit, while the other was a manual sampler of higher daily air volume. The samples collected each day showed excellent correlation in radionuclide activity, although some variations were seen. These variations illustrate the reproducibility of the air sample radionuclide measurements made by the Comprehensive Nuclear-Test-Ban Treaty International Monitoring System (IMS) and show a simple way to acquire useful parallel samples for scientific purposes. In particular, a party wishing to have a “copy” of a sample acquired by the verification regime of the treaty could employ this method and have results similar to the IMS station at low cost and even higher sensitivity.

Atmospheric Aerosol Sampling↗
DOE OSTI2019

MPI jobs within MPI jobs: A practical way of enabling task-level fault-tolerance in HPC workflows

While the use of workflows for HPC is growing, MPI interoperability remains a challenge for workflow management systems. The MPI standard and/or its implementations provide a number of ways to build multiple-programs-multiple-data (MPMD) applications. These methods present limitations related to fault tolerance, and are not easy to use. In this paper, we advocate for a novel MPI_Comm_launch function acting as the parallel counterpart of a system(3) call. MPI_Comm_launch allows a child MPI application to be launched inside the resources originally held by processes of a parent MPI application. Two important aspects of MPI_Comm_launch is that it pauses the calling process, and runs the child processes on the parent’s CPU cores, but in an isolated manner with respect to memory. This function makes it easier to build MPMD applications with well-decoupled subtasks. We show how this feature can provide better flexibility and better fault tolerance in ensemble simulations and HPC workflows. Here, we report results showing 2x throughput improvement for application workflows with faults, and scaling results for challenging workloads up to 256 nodes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2016

Self-Assembly and Critical Aggregation Concentration Measurements of ABA Triblock Copolymers with Varying B Block Types: Model Development, Prediction, and Validation

The dissipative particle dynamics (DPD) simulation technique is a coarse-grained (CG) molecular dynamics-based approach that can effectively capture the hydrodynamics of complex systems while retaining essential information about the structural properties of the molecular species. An advantageous feature of DPD is that it utilizes soft repulsive interactions between the beads, which are CG representation of groups of atoms or molecules. In this study, we used the DPD simulation technique to study the aggregation characteristics of ABA triblock copolymers in aqueous medium. Pluronic polymers (PEG-PPO-PEG) were modeled as two segments of hydrophilic beads and one segment of hydrophobic beads. Tyrosine-derived PEG 5K -b-oligo(desaminotyrosyl tyrosine octyl ester-suberate)-b-PEG 5K (PEG5K-oligo(DTO-SA)-PEG 5K ) block copolymers possess alternate rigid and flexible components along the hydrophobic oligo(DTO-SA) chain, and were modeled as two segments of hydrophilic beads and one segment of hydrophobic, alternate soft and hard beads. The formation, structure, and morphology of the initial aggregation of the polymer molecules in aqueous medium were investigated by following the aggregation dynamics. The dimensions of the aggregates predicted by the computational approach were in good agreement with corresponding results from experiments, for the Pluronic and PEG 5K -oligo(DTO-SA)-PEG 5K block copolymers. In addition, DPD simulations were utilized to determine the critical aggregation concentration (CAC), which was compared with corresponding results from an experimental approach. For Pluronic polymers F68, F88, F108, and F127, the computational results agreed well with experimental measurements of the CAC measurements. For PEG 5K -b-oligo(DTO-SA)-b-PEG 5K block polymers, the complexity in polymer structure made it difficult to directly determine their CAC values via the CG scheme. Therefore, we determined CAC values of a series of triblock copolymers with 3–8 DTO-SA units using DPD simulations, and used these results to predict the CAC values of triblock copolymers with higher molecular weights by extrapolation. In parallel, a PEG 5K -b-oligo(DTO-SA)-b-PEG 5K block copolymer was synthesized, and the CAC value was determined experimentally using the pyrene method. The experimental CAC value agreed well with the CAC value predicted by simulation. Furthermore, these results validate our CG models, and demonstrate an avenue to simulate and predict aggregation characteristics of ABA amphiphilic triblock copolymers with complex structures.

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

Near-field radiative heat transfer between two parallel SiO{sub 2} plates with and without microcavities

Near-to-far-field radiative heat transfer between two macroscopic SiO{sub 2} plates—with and without microcavities—was observed using a highly precise and accurate optical gap-measurement method. The experiments, conducted near 300 K, measured heat transfer as a function of gap separation from 1.0 μm to 50 μm and also as a function of temperature differences between 4.1 and 19.5 K. The gap-dependent heat flux was in excellent agreement with theoretical predictions. Furthermore, the effects of microcavities on the plate surfaces were clearly observed and significant enhancement of near-field radiative heat transfer was confirmed between gold-coated microcavities with narrow vacuum separation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Performance evaluation of OpenFOAM on many-core architectures

In this article application of Open Source Field Operation and Manipulation (OpenFOAM) C++ libraries for solving engineering problems on many-core architectures is presented. Objective of this article is to present scalability of OpenFOAM on parallel platforms solving real engineering problems of fluid dynamics. Scalability test of OpenFOAM is performed using various hardware and different implementation of standard PCG and PBiCG Krylov iterative methods. Speed up of various implementations of linear solvers using GPU and MIC accelerators are presented in this paper. Numerical experiments of 3D lid-driven cavity flow for several cases with various number of cells are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

LibERI—A portable and performant multi-GPU accelerated library for electron repulsion integrals via OpenMP offloading and standard language parallelism

A portable and performant graphics processing unit (GPU)-accelerated library for electron repulsion integral (ERI) evaluation, named LibERI, has been developed and implemented via directive-based (e.g., OpenMP and OpenACC) and standard language parallelism (e.g., Fortran DO CONCURRENT). Offloaded ERIs consist of integrals over low and high contraction s, p, and d functions using the rotated-axis and Rys quadrature methods. GPU codes are factorized based on previous developments with two layers of integral screening and quartet presorting. In this work, the density screening is moved to the GPU to enhance the computational efficacy for large molecular systems. Here, the L-shells in the Pople basis set are also separated into pure S and P shells to increase the ERI homogeneity and reduce atomic operations and the memory footprint. LibERI is compatible with any quantum chemistry drivers supporting the MolSSI Driver Interface. Benchmark calculations of LibERI interfaced with the GAMESS software package were carried out on various GPU architectures and molecular systems. The results show that the LibERI performance is comparable to other state-of-the-art GPU-accelerated codes (e.g., TeraChem and GMSHPC) and, in some cases, outperforms conventionally developed ERI CUDA kernels (e.g., QUICK) while fully maintaining portability.

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

Large Scale GPU Accelerated PPMLR-MHD Simulations for Space Weather Forecast

PPMLR-MHD is a new magnetohydrodynamics (MHD) model used to simulate the interactions of the solar wind with the magnetosphere, which has been proved to be the key element of the space weather cause-and-effect chain process from the Sun to Earth. Compared to existing MHD methods, PPMLR-MHD achieves the advantage of high order spatial accuracy and low numerical dissipation. However, the accuracy comes at a cost. On one hand, this method requires more intensive computation. On the other hand, more boundary data is subject to be transferred during the process of simulation. In this work, we present a parallel hybrid solution of the PPMLR-MHD model implemented using the computing capabilities of both CPUs and GPUs. We demonstrate that our optimized implementation alleviates the data transfer overhead by using GPU Direct technology and can scale up to 151 processes and achieve significant performance gains by distributing the workload among the CPUs and GPUs on Titan at Oak Ridge National Laboratory. The performance results show that our implementation is fast enough to carry out highly accurate MHD simulations in real time.

Computer Science↗
DOE OSTI2017

A Unit Testing Framework for Scientific Legacy Code

Large-scale scientific applications play important roles in supporting research. However, it is often very expensive and time-consuming to make changes to, maintain and evolve the scientific code due to its complexity and poor programming skills of researchers. Therefore, in order to visualize scientific code architecture to optimize software design, understand undocumented source code, and analyze software flow and functionality, we first introduce a unit testing framework (UTF). Then, because such infrastructure's performance is very crucial in practical use since the scientific legacy applications simulate instances in a long period of time, we improve the UTF by applying Message Passing based Parallelization and parallel I/O operations. Furthermore, due to the scientific code has enormous state data and the I/O capacity on the server is limited, we apply in situ data analysis method to encounter fewer resource limitations, and adopt signal processing to greatly reduce data transfer. Last, we demonstrated the correctness and high-efficiency of our framework for legacy Earth model on Titan supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

WE-G-17A-08: Electron Gun Operation for in Line MRI-Linac Configurations: An Assessment of Beam Fidelity and Recovery Techniques for Different SIDs and Magnetic Field Strengths

Purpose: To test the functionality of medical electron guns within the fringe field of a purpose built superconducting MRI magnet, and to test different recovery techniques for a variety of imaging field strengths and SIDs. Methods: Three different electron guns were simulated using Finite Element Modelling; a standard diode gun, a standard triode gun, and a novel diode gun designed to operate within parallel magnetic fields. The approximate working regime of each gun was established by assessing exit current in constant magnetic fields of varying strength and defining ‘working’ as less than 10% change in injection current. Next, the 1.0T MRI magnet was simulated within Comsol Multiphysics. The coil currents in this model were also scaled to produce field strengths of .5, 1, 1.5 and 3T. Various magnetic shield configurations were simulated, varying the SID from 800 to 1300mm. The average magnetic field within the gun region was assessed together with the distortion in the imaging volume - greater than 150uT distortion was considered unacceptable. Results: The conventional guns functioned in fields of less than 7.5mT. Conversely, the redesigned diode required fields greater than .1T to function correctly. Magnetic shielding was feasible for SIDS of greater than 1000mm for field strengths of .5T and 1T, and 1100mm for 1.5 and 3.0T. Beyond these limits shielding resulted in unacceptable MRI distortion. In contrast, the redesigned diode could perform acceptably for SIDs of less than 812, 896, 931, and 974mm for imaging strengths of 0.5, 1.0, 1.5, 3.0T. Conclusions: For in-line MRIlinac configurations where the electron gun is operating in low field regions, shielding is a straight forward option. However, as magnetic field strength increases and the SID is reduced, shielding results in too great a distortion in the MRI and redesigning the electron optics is the preferable solution. The authors would like to acknowledge funding from the National Health and Research Council (AUS), National Institute of Health (NIH), and Cancer Institute NSW.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Ion recombination correction in carbon ion beams

Purpose: In this work, ion recombination is studied as a function of energy and depth in carbon ion beams. Methods: Measurements were performed in three different passively scattered carbon ion beams with energies of 62 MeV/n, 135 MeV/n, and 290 MeV/n using various types of plane-parallel ionization chambers. Experimental results were compared with two analytical models for initial recombination. One model is generally used for photon beams and the other model, developed by Jaffé, takes into account the ionization density along the ion track. An investigation was carried out to ascertain the effect on the ion recombination correction with varying ionization chamber orientation with respect to the direction of the ion tracks. The variation of the ion recombination correction factors as a function of depth was studied for a Markus ionization chamber in the 62 MeV/n nonmodulated carbon ion beam. This variation can be related to the depth distribution of linear energy transfer. Results: Results show that the theory for photon beams is not applicable to carbon ion beams. On the other hand, by optimizing the value of the ionization density and the initial mean-square radius, good agreement is found between Jaffé’s theory and the experimental results. As predicted by Jaffé’s theory, the results confirm that ion recombination corrections strongly decrease with an increasing angle between the ion tracks and the electric field lines. For the Markus ionization chamber, the variation of the ion recombination correction factor with depth was modeled adequately by a sigmoid function, which is approximately constant in the plateau and strongly increasing in the Bragg peak region to values of up to 1.06. Except in the distal edge region, all experimental results are accurately described by Jaffé’s theory. Conclusions: Experimental results confirm that ion recombination in the investigated carbon ion beams is dominated by initial recombination. Ion recombination corrections are found to be significant and cannot be neglected for reference dosimetry and for the determination of depth dose curves in carbon ion beams.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗