Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Software Validation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34Linked to original sources

Dosimetric comparison of Monte Carlo codes (EGS4, MCNP, MCNPX) considering external and internal exposures of the Zubal phantom to electron and photon sources.

This paper aims at comparing dosimetric assessments performed with three Monte Carlo codes: EGS4, MCNP4c2 and MCNPX2.5e, using a realistic voxel phantom, namely the Zubal phantom, in two configurations of exposure. The first one deals with an external irradiation corresponding to the example of a radiological accident. The results are obtained using the EGS4 and the MCNP4c2 codes and expressed in terms of the mean absorbed dose (in Gy per source particle) for brain, lungs, liver and spleen. The second one deals with an internal exposure corresponding to the treatment of a medullary thyroid cancer by 131I-labelled radiopharmaceutical. The results are obtained by EGS4 and MCNPX2.5e and compared in terms of S-values (expressed in mGy per kBq and per hour) for liver, kidney, whole body and thyroid. The results of these two studies are presented and differences between the codes are analysed and discussed.

Computer Simulation↗

Theoretical simulation of residual nuclide products in 208,207,206Pb, NATPb and 209Bi (P,X) reactions at intermediate and high energies.

The independent and cumulative measured yields of residual products in thin lead and bismuth targets irradiated with 0.04-2.6 GeV protons are compared with results by the LAHET, CEM03, LAQGSM03, INCL+ABLA, CASCADE and YIELDX codes, in order to evaluate the predictive power of the codes in this energy region. We found that the predictive power of the tested codes is different but is satisfactory for most of the nuclides in the spallation region, though none of the codes agree well with the data in the whole-mass region of product nuclides and all should be improved further. On the whole, the predictive power of all codes for the data in the fission and fragmentation product regions and, especially, at the borders between spallation and fission and between fission and fragmentation regions is much worse than in the spallation region; therefore, development of better evaporation/fission/fragmentation models is of first priority.

Bismuth↗

Recent enhancements to the MARS15 code.

The MARS code is under continuous development and has recently undergone substantial improvements that further increase its reliability and predictive power in numerous shielding, accelerator, detector and space applications. The major developments and new features of the MARS15 (2004) version described in this paper concern an extended list of elementary particles and arbitrary heavy ions and their interaction cross sections, inclusive and exclusive nuclear event generators, module for modelling particle electromagnetic interactions, enhanced geometry and histogramming options, improved MAD-MARS Beam Line Builder, enhanced graphical user interface and an MPI-based parallelisation of the code.

Algorithms↗

Towards a heavy-ion transport capability in the MARS15 code.

In order to meet the challenges of new accelerator and space projects and further improve modelling of radiation effects in microscopic objects, heavy-ion interaction and transport physics have been recently incorporated into the MARS15 Monte Carlo code. A brief description of new modules is given in comparison with experimental data.

Algorithms↗

Preliminary shielding assessment for the 100 MeV proton linac (KOMAC).

The Proton Engineering Frontier Project is building the Korea Multipurpose Accelerator Complex facilities from 2002 to 2012, which consists of a high-current 100 MeV proton linear accelerator and various beam-lines. This paper provides a preliminary estimate of the shielding required for the 20 mA proton linac and the beam-dump. For an accurate information on secondary neutron production from the guiding magnet and primary heat sink of the beam dump, proton-induced 63Cu and 65Cu cross section data were evaluated and applied to shielding calculations. The required thickness of the concrete was assessed by a simple line-of-sight model for the lateral shielding of the beam-line and the full shielding of the beam dump. Monte Carlo simulations were also performed using the MCNPX code to obtain the source term and attenuation coefficients for the three-dimensional lateral shielding model of the beam-line.

Computer Simulation↗

Radiation transport calculations for 50 MV photon therapy beam using the Monte Carlo code GEANT4.

A new thin transmission target technique for fast dose delivery using narrow scanned photon beams has been developed. High-energy, 50-100 MeV, electron beams of low emittance incident on thin low-Z targets produce narrow and intense high-energy bremsstrahlung beams. However, electrons transmitted through the target are bent from the therapeutic beam by a purging magnet and have to be effectively absorbed in a dedicated electron collector. The electron-photon transport through a treatment head has been studied using the Monte Carlo simulation toolkit Geant4. The Geant4 electromagnetic physics processes have been compared with experimental data of radial dose profiles. The differences between calculated and measured radial dose distributions are approximately 2-10%. Preliminary investigations of the collector design have been carried out in order to minimise secondary electron and photon contamination of the therapeutic beam. The toolkit presented here is promising for further development of narrow photon beam therapy.

Computer Simulation↗

Benchmarking of MCNP for calculating dose rates at an interim storage facility for nuclear waste.

During the operation of research facilities at Research Centre Jülich, Germany, nuclear waste is stored in drums and other vessels in an interim storage building on-site, which has a concrete shielding at the side walls. Owing to the lack of a well-defined source, measured gamma spectra were unfolded to determine the photon flux on the surface of the containers. The dose rate simulation, including the effects of skyshine, using the Monte Carlo transport code MCNP is compared with the measured dosimetric data at some locations in the vicinity of the interim storage building. The MCNP data for direct radiation confirm the data calculated using a point-kernel method. However, a comparison of the modelled dose rates for direct radiation and skyshine with the measured data demonstrate the need for a more precise definition of the source. Both the measured and the modelled dose rates verified the fact that the legal limits (<1 mSv a(-1)) are met in the area outside the perimeter fence of the storage building to which members of the public have access. Using container surface data (gamma spectra) to define the source may be a useful tool for practical calculations and additionally for benchmarking of computer codes if the discussed critical aspects with respect to the source can be addressed adequately.

Benchmarking↗

Comparison of deterministic and Monte Carlo methods in shielding design.

In shielding calculation, deterministic methods have some advantages and also some disadvantages relative to other kind of codes, such as Monte Carlo. The main advantage is the short computer time needed to find solutions while the disadvantages are related to the often-used build-up factor that is extrapolated from high to low energies or with unknown geometrical conditions, which can lead to significant errors in shielding results. The aim of this work is to investigate how good are some deterministic methods to calculating low-energy shielding, using attenuation coefficients and build-up factor corrections. Commercial software MicroShield 5.05 has been used as the deterministic code while MCNP has been used as the Monte Carlo code. Point and cylindrical sources with slab shield have been defined allowing comparison between the capability of both Monte Carlo and deterministic methods in a day-by-day shielding calculation using sensitivity analysis of significant parameters, such as energy and geometrical conditions.

Algorithms↗

Three-dimensional shielding calculations for the IFMIF neutron source using a coupled Monte Carlo/Deterministic computational scheme.

Shielding calculations for the International Fusion Materials Irradiation Facility (IFMIF) are complicated due to the geometrical complexity of the target system and the large-scale bulk shields around the source target. Three-dimensional shielding calculations were performed by using a newly developed Monte Carlo/Deterministic computational scheme. The neutron-photon fluxes and dose rate distributions in the back wall of the Test Cell and the access/maintenance room are presented and compared with previous shielding calculations. The results demonstrate that this coupled scheme is an useful computational tool for three-dimensional shielding analyses of complex and large nuclear facilities.

Computer Simulation↗

QUADOS intercomparison: a summary of photon and charged particle problems.

QUADOS, a Concerted Action of the European Commission, has promoted an intercomparison aimed at evaluating the use of computational codes for dosimetry in radiation protection and medical physics. This intercomparison was open to all users of radiation transport codes. Eight problems were selected for their relevance to the radiation dosimetry community, five of which involved photon and proton transport. This paper focuses on the analysis of the photon and charged particle problems. The neutron problems were presented in a paper at the NEUDOS9 conference.

Algorithms↗

Comparison of synchrotron radiation calculations between analytical codes (STAC8, photon) and Monte Carlo codes (FLUKA, EGS4).

STAC8 is a significant improvement upon its predecessor PHOTON and is a valuable analytical code for quick and conservative beamline shielding designs for synchrotron radiation (SR) facilities. In order to check the applicability, accuracy and limitations of STAC8, studies were conducted to compare the results of STAC8 with those of PHOTON with calculations using the FLUKA and EGS4 Monte Carlo codes. Doses and spectra for scattered SR in a few beam-target-shield geometries were calculated, with and without photon linear polarisation effects. Areas for expanding the STAC8 capabilities, for example, features of the mirror-reflected lights and double-Compton light calculations and the use of monochromatic light have been identified. Some of these features have been implemented and benchmarked against Monte Carlo calculations. Reasonable agreements were found between the STAC8 and Monte Carlo calculations.

Algorithms↗

MCNPX vs. DORT for SNS shielding design studies.

Radiation transport occurs through the 18 m long access way adjacent to the Spallation Neutron Source accelerator tunnel and the 2.2 m thick massive shielding door which closes the access way. A variety of typical materials for accelerator shielding, such as concrete and steel, were used for construction of the door to study radiation penetration. A comparison was carried out using both Monte Carlo (code MCNPX) and discrete ordinates (code DORT) methods. The beam losses during the accelerator operation are the sources for the radiation calculations. Analyses show that the results from the two methods are in good agreement.

Computer Simulation↗

Development of a connection-method calculation utility for the DORT-TORT code.

A code package has been developed to utilise connection-method calculation for both 2-D and 3-D SN codes, DORT-TORT, extending further the capabilities of DORSET/TORSET. It also includes support utilities to make input preparation easier. Potential uses of the method are, for example, to allow the coupling of one problem to the end of another by transforming output fluxes of the former into an external source term for the latter. For evaluating the code's applicability in a more specific way, the coupling calculation has been performed using a two-leg duct-streaming problem with a 14 MeV neutron source. The RZ-RZ (DORT) result shows good overall agreement with the reference MCNP4B result. However, it is also found that the result is influenced by the dimension of the calculation region and that some transformation error has occurred. For evaluating the 3-D transformation, full preparation of first collision source is required for TORT.

Algorithms↗

Radiation safety design for the J-PARC Project.

The High-Intensity Proton Accelerator Project, named J-PARC, is in progress, with the aim of enabling studies on the latest basic science and the advancement of nuclear technology. In the project, a high-energy proton accelerator complex with the world's highest instantaneous intensity is under construction. In order to establish a reasonable shielding design, both simplified and detailed design methods were used in the shielding design of J-PARC. This paper reviews the present status of the radiation safety design study for J-PARC.

Computer Simulation↗

Skyshine analysis using energy and angular dependent dose-contribution fluxes obtained from air-over-ground adjoint calculation.

A generalised and convenient skyshine dose analysis method has been developed based on forward-adjoint folding technique. In the method, the air penetration data were prepared by performing an adjoint DOT3.5 calculation with cylindrical air-over-ground geometry having an adjoint point source (importance of unit flux to dose rate at detection point) in the centre. The accuracy of the present method was certified by comparing with DOT3.5 forward calculation. The adjoint flux data can be used as generalised radiation skyshine data for all sorts of nuclear facilities. Moreover, the present method supplies plenty of energy-angular dependent contribution flux data, which will be useful for detailed shielding design of facilities.

Air↗

Therapeutic step and shoot proton beam spot-scanning with a multi-leaf collimator: a Monte Carlo study.

In step and shoot spot-scanning, a small-diameter proton beam is magnetically swept and varied in energy in order to cover the tumour. Initial estimates of the beam size indicate that additional collimating hardware will be needed for lower energy proton beams in order to achieve a clinically acceptable lateral dose falloff at the edge of the proton beam. In this report, we present dosimetric data from Monte Carlo simulations with a model of a simple multileaf collimator which indicate that such a device may be used to improve the lateral dose falloff. The dosimetric quantities relevant to the clinical usefulness of the device are studied, including lateral penumbra, leaf transmission and scalloping effect. Multileaf collimation is compared with a differential spot-weighting technique of sharpening the lateral dose falloff.

Body Burden↗

A comparison of skyshine computational methods.

A variety of methods employing radiation transport and point-kernel codes have been used to model two skyshine problems. The first problem is a 1 MeV point source of photons on the surface of the earth inside a 2 m tall and 1 m radius silo having black walls. The skyshine radiation downfield from the point source was estimated with and without a 30-cm-thick concrete lid on the silo. The second benchmark problem is to estimate the skyshine radiation downfield from 12 cylindrical canisters emplaced in a low-level radioactive waste trench. The canisters are filled with ion-exchange resin with a representative radionuclide loading, largely 60Co, 134Cs and 137Cs. The solution methods include use of the MCNP code to solve the problem by directly employing variance reduction techniques, the single-scatter point kernel code GGG-GP, the QADMOD-GP point kernel code, the COHORT Monte Carlo code, the NAC International version of the SKYSHINE-III code, the KSU hybrid method and the associated KSU skyshine codes.

Air↗

Role of contact and noncontact mapping in the curative ablation of tachyarrhythmias.

Assessment of the timing of electrical activation recorded by multiple electrodes positioned in various locations within the heart has been the conventional method for mapping cardiac arrhythmias. This technique requires fluoroscopy for catheter manipulation, which in addition to being harmful (ionizing radiation), is inadequate for visualizing the complex three-dimensional cardiac anatomy and lacks reproducibility regarding localization of sites of interest. Because of these limitations, several new mapping systems that can function in a complimentary role to the conventional mapping technique, or can be used independently, have been developed. These new mapping strategies have unique advantages. They overcome the limitations of fluoroscopy by creating accurate three-dimensional intracardiac maps. The ability to localize and accurately display intracardiac catheter positioning and ablation lesion sites facilitate increasingly complex catheter ablation procedures.

Body Surface Potential Mapping↗