Search

The Monte Carlo (MC) method has become an indispensable part of the nuclear radiation research field. Several widely used and well-known MC packages were developed for simulation of radiation transport and interaction with matter. All these MC packages require users to prepare an input script. The input script can become lengthy for complex models. The process of preparing these input scripts is time-consuming and error-prone. In the present work, we have developed an open-source GUI computer program for modelling radiation transport and interaction in multi-segmented slab phantoms using grid-based system for the widely used PHITS MC package. The developed tools would be useful for future users of PHITS MC package and particularly inexperienced users. The present program is distributed under GPL license and all users can freely download, modify and redistribute the program without any restrictions.

Abstract The present work introduced a framework to investigate the effectiveness of proton boron fusion therapy (PBFT) at the cellular level. The framework consisted of a cell array generator program coupled with PHITS Monte Carlo package with a dedicated terminal-based code editor that was developed in this work. The framework enabled users to model large cell arrays with normal, all boron, and random boron filled cytoplasm, to investigate the underlying mechanism of PBFT. It was found that alpha particles and neutrons could be produced in absence of boron mainly because of nuclear reaction induced by proton interaction with 16O, 12C and 14N nuclei. The effectiveness of PBFT is highly dependent on the incident proton energy, source size, cell array size, buffer medium thickness layer, concentration and distribution of boron in the cell array. To quantitatively assess the effectiveness of PBFT, of the total energy deposition by alpha particle for different cases were determined. The number of alpha particle hits in cell cytoplasm and nucleus for normal and 100 ppm boron were determined. The obtained results and the developed tools would be useful for future development of PBFT to objectively determine the effectiveness of this treatment modality.

Proton irradiations are highly sensitive to spatial variations, mainly due to their high linear energy transfer (LET) and densely ionizing nature. In realistic clinical applications, the targets of ionizing radiation are inhomogeneous in terms of geometry and chemical composition (i.e., organs in the human body). One of the main methods for proton range monitoring is to utilize the production of proton induced positron emitting radionuclides; these could be measured precisely with positron emission tomography (PET) systems. One main positron emitting radionuclide that could be used for proton range monitoring and verification was found to be 13N that produces a peak close to the Bragg peak. In the present work, we have employed the Monte Carlo method and Spectral Analysis (SA) technique to investigate the feasibility of utilizing the 13N peak for proton range monitoring and verification in inhomogeneous targets. Two different phantom types, namely, (1) ordinary slab and (2) MIRD anthropomorphic phantoms, were used. We have found that the generated 13N peak in such highly inhomogeneous targets (ordinary slab and human phantom) is close to the actual Bragg peak, when irradiated by incident proton beam. The feasibility of using the SA technique to estimate the distribution of positron emitter was also investigated. The current findings and the developed tools in the present work would be helpful in proton range monitoring and verification in realistic clinical radiation therapy using proton beams.

The proto-diastolic third heart sound (S3) is observed in various hemodynamic conditions in both normal and diseased hearts. We propose a novel, one-degree of freedom mathematical model of mechanical vibrations of heart and blood that generates the third heart sound, implemented in a real-time model of the cardiovascular system (CircAdapt). To examine model functionality, S3 simulations were performed for conditions mimicking the normal heart as well as heart failure with preserved ejection fraction (HFpEF), atrioventricular valve regurgitation (AVR), atrioventricular valve stenosis (AVS) and septal shunts (SS). Simulated S3 showed both qualitative and quantitative agreements with measured S3 in terms of morphology, frequency, and timing. It was shown that ventricular mass, ventricular viscoelastic properties as well as inflow momentum play a key role in the generation of S3. The model indicated that irrespective of cardiac conditions, S3 vibrations are always generated, in both the left and right sides of the heart, albeit at different levels of audibility. S3 intensities increased in HFpEF, AVR and SS, but the changes of acoustic S3 features in AVS were not significant, as compared with the reference simulation. S3 loudness in all simulated conditions was proportional to the level of cardiac output and severity of cardiac conditions. In conclusion, our hemodynamics-driven mathematical model provides a fast and realistic simulation of S3 under various conditions which may be helpful to find new indicators for diagnosis and prognosis of cardiac diseases.

ObjectiveA method to estimate absolute left ventricular (LV) pressure and its maximum rate of rise (LV dP/dtmax) from epicardial accelerometer data and machine learning is proposed.MethodsFive acute experiments were performed on pigs. Custom-made accelerometers were sutured epicardially onto the right ventricle, LV, and right atrium. Different pacing configurations and contractility modulations, using isoflurane and dobutamine infusions, were performed to create a wide variety of hemodynamic conditions. Automated beat-by-beat analysis was performed on the acceleration signals to evaluate amplitude, time, and energy-based features. For each sensing location, bootstrap aggregated classification tree ensembles were trained to estimate absolute maximum LV pressure (LVPmax) and LV dP/dtmax using amplitude, time, and energy-based features. After extraction of acceleration and pressure-based features, location specific, bootstrap aggregated classification ensembles were trained to estimate absolute values of LVPmax and its maximum rate of rise (LV dP/dtmax) from acceleration data.ResultsWith a dataset of over 6,000 beats, the algorithm narrowed the selection of 17 predefined features to the most suitable 3 for each sensor location. Validation tests showed the minimal estimation accuracies to be 93% and 86% for LVPmax at estimation intervals of 20 and 10 mmHg, respectively. Models estimating LV dP/dtmax achieved an accuracy of minimal 93 and 87% at estimation intervals of 100 and 200 mmHg/s, respectively. Accuracies were similar for all sensor locations used.ConclusionUnder pre-clinical conditions, the developed estimation method, employing epicardial accelerometers in conjunction with machine learning, can reliably estimate absolute LV pressure and its first derivative.

The interaction of ionizing radiation with matter is a stochastic process and statistical analysis of such a process would be a crucial step in understanding radioactivity. Geiger–Müller (GM) counter is a widely used radiation detector used in nuclear radiation surveying, which produces counts upon exposure to a radioactive source. There are a variety of multi-purpose software that can be used to perform statistical analysis of measured counts from a GM counter. However, statistical analysis is a lengthy, error prone and time-consuming process, which gets more tedious when the number of measurements increases. In the present work, we have developed an open-source and easy-to-use graphical user interface (GUI) computer program named RadStat for statistical analysis of counts measured by a GM counter. RadStat has its own scripting syntaxes and bundled with gnuplot for quick visualization of output results. We believe the present open-source GUI program would be a useful tool for research and teaching of nuclear radiation physics.

The Monte Carlo method is employed in this study to simulate the proton irradiation of a water-gel phantom. Positron-emitting radionuclides such as 11C, 15O, and 13N are scored using the Particle and Heavy Ion Transport Code System Monte Carlo code package. Previously, it was reported that as a result of 16O(p,2p2n)13N nuclear reaction, whose threshold energy is relatively low (5.660 MeV), a 13N peak is formed near the actual Bragg peak. Considering the generated 13N peak, we obtain offset distance values between the 13N peak and the actual Bragg peak for various incident proton energies ranging from 45 to 250 MeV, with an energy interval of 5 MeV. The offset distances fluctuate between 1.0 and 2.0 mm. For example, the offset distances between the 13N peak and the Bragg peak are 2.0, 2.0, and 1.0 mm for incident proton energies of 80, 160, and 240 MeV, respectively. These slight fluctuations for different incident proton energies are due to the relatively stable energy-dependent cross-section data for the 16O(p,2p2n)13N nuclear reaction. Hence, we develop an open-source computer program that performs linear and non-linear interpolations of offset distance data against the incident proton energy, which further reduces the energy interval from 5 to 0.1 MeV. In addition, we perform spectral analysis to reconstruct the 13N Bragg peak, and the results are consistent with those predicted from Monte Carlo computations. Hence, the results are used to generate three-dimensional scatter plots of the 13N radionuclide distribution in the modeled phantom. The obtained results and the developed methodologies will facilitate future investigations into proton range monitoring for therapeutic applications.

This article incorporates details of our calculations of the observable quantities for the scattering of electrons and positrons from a post transition metal Thallium (Tl), in the energy range of 1 eV ≤ Ei ≤ 1 MeV, using the relativistic Dirac partial wave (phase-shift) analysis employing a complex optical-potential. Absolute differential, integrated elastic and inelastic, transport, total ionization, and total cross sections and a thorough study of the critical minima in the elastic differential cross sections along with the associated angular positions of the maximum polarization points in the Sherman function are provided to study the collision dynamics. The optical potential model incorporates the interactions of the incident electron and/or positron with both the nucleus and the bound electrons of the target atom. In-depth analyses of the spin asymmetry, which are sensitive to phases related interference effect, brought on by the various ingredients of the lepton-atom interaction, are also presented. The performance of the current approach to explain the observations, with the exception of the extremely low energy domain, is shown by a comparison of the previous experimental and theoretical results on this target atom.

New & noteworthyTo the best of our knowledge, this is the first hemodynamic-based heart sound generation model embedded in a complete real-time computational model of the cardiovascular system. Simulated heart sounds are similar to experimental and clinical measurements, both quantitatively and qualitatively. Our model can be used to investigate the relationships between heart sound acoustic features and hemodynamic factors/anatomical parameters.

Abstract A calibration method was proposed in the present work to determine the medium-thickness-dependent proton doses absorbed in cellular components (i.e., cellular cytoplasm and nucleus) in radiobiological experiments. Consideration of the dependency on medium thickness was crucial as the linear energy transfer (LET) of protons could rise to a sharp peak (known as the Bragg peak) towards the end of their ranges. Relationships between the calibration coefficient R vs medium-layer thickness were obtained for incident proton energies of 10, 15, 20, 25, 30 and 35 MeV, and for various medium thicknesses up to 5000 μm, where R was defined as the ratio D A /D E , D A was the absorbed proton dose in cellular components, and D E was the absorbed proton dose in a separate radiation detector. In the present work, D A and D E were determined using the MCNPX (Monte Carlo N-Particle eXtended) code version 2.4.0. For lower incident proton energies (i.e., 10, 15 and 20 MeV), formation of Bragg-peak-like features were noticed in their R-vs-medium-layer-thickness relationships, and large R values of >7 and >6 were obtained for cytoplasm and nucleus of cells, respectively, which highlighted the importance of careful consideration of the medium thickness in radiobiological experiments.

In present work, molecular dynamics (MD) simulation was employed to model the impact of single and double molten aluminum drops onto to an aluminum substrate. The study was related to the pileup process in metal deposition manufacture. Most previous investigations used the continuum scale Volume-of-Fluid (VOF) method to model this phenomenon. Here, the underlying mechanism of pileup process was examined in more details. The effect of impact speed on the height and width of the solidified drops was studied. The penetration depth of molten drops into the solid substrate was also investigated. The present model would be useful for future investigations as the process could now be studied at the atomic scale and the final solidified structure could also be assessed in terms of its defect and crystallinity. Keywords: Molecular dynamics, Pileup, Metal deposition, Atomistic modeling

The Monte Carlo method was employed to simulate realistic treatment situations for photon and proton radiation therapy for a set of Oak Ridge National Laboratory (ORNL) pediatric phantoms for 15, 10, 5 and 1-year olds as well as newborns. Complete radiotherapy situations were simulated using the previously developed NRUrad input code for Monte Carlo N-Particle (MCNP) code package. Each pediatric phantom was irradiated at five different positions, namely, the testes, colon, liver, left lung and brain, and the doses in targeted organs (Dt) were determined using the track length estimate of energy. The dispersed photon and proton doses in non-targeted organs (Dd), namely, the skeleton, skin, brain, spine, left and right lungs were computed. The conversion coefficients (F = Dd/Dt) of the dispersed doses were used to study the dose dispersion in different non-targeted organs for phantoms for 15, 10, 5 and 1-year olds as well as newborns. In general, the F values were larger for younger patients. The F values for non-targeted organs for phantoms for 1-year olds and newborns were significantly larger compared to those for other phantoms. The dispersed doses from proton radiation therapy were also found to be significantly lower than those from conventional photon radiation therapy. For example, the largest F values for the brain were 65.6% and 0.206% of the dose delivered to the left lung (P4) for newborns during photon and proton radiation therapy, respectively. The present results demonstrated that dispersion of photons and generated electrons significantly affected the absorbed doses in non-targeted organs during pediatric photon therapy, and illustrated that proton therapy could in general bring benefits for treatment of pediatric cancer patients.

Some concepts in nuclear radiation physics are abstract and intellectually demanding. In the present paper, an "MCHP platform" (MCHP was an acronym for Monte Carlo simulations + Human Phantoms) was proposed to provide assistance to the students through visualization. The platform involved Monte Carlo simulations of interactions between ionizing radiations and the Oak Ridge National Laboratory (ORNL) adult male human phantom. As an example to demonstrate the benefits of the proposed MCHP platform, the present paper investigated the variation of the absorbed photon dose per photon from a 137Cs source in three selected organs, namely, brain, spine and thyroid of an adult male for concrete and lead shields with varying thicknesses. The results were interesting but not readily comprehensible without direct visualization. Graphical visualization snapshots as well as video clips of real time interactions between the photons and the human phantom were presented for the involved cases, and the results were explained with the help of such snapshots and video clips. It is envisaged that, if the platform is found useful and effective by the readers, the readers can also propose examples to be gradually added onto this platform in future, with the ultimate goal of enhancing students' understanding and learning the concepts in an undergraduate nuclear radiation physics course or a related course.

Proton range monitoring and verification is important to enhance the effectiveness of treatment by ensuring that the correct dose is delivered to the correct location. Upon proton irradiation, different positron emitting radioisotopes are produced by the inelastic nuclear interactions of protons with the target elements. Recently, it was reported that the 16O(p,2p2n)13N reaction has a relatively low threshold energy, and it could be potentially used for proton range verification. In the present work, we have proposed an analysis scheme (i.e., algorithm) for the extraction and three-dimensional visualization of positron emitting radioisotopes. The proposed step-by-step analysis scheme was tested using our own experimentally obtained dynamic data from a positron emission mammography (PEM) system (our developed PEMGRAPH system). The experimental irradiation was performed using an azimuthally varying field (AVF) cyclotron with a 80 MeV monoenergetic pencil-like beam. The 3D visualization showed promising results for proton-induced radioisotope distribution. The proposed scheme and developed tools would be useful for the extraction and 3D visualization of positron emitting radioisotopes and in turn for proton range monitoring and verification.

Data security has become indispensable, with a view to keep sensitive information confidential. One important method is through image encryption, upon which features in an image would no longer be visible. The original image with its features could only be restored upon decryption using a set of keys. There are prestigious works in the literature regarding image encryption. However, there is a lack of easy-to-use, GUI-based, user-customizable computer programs for image encryption. In the present work, we developed a GUI-based image encryption and decryption program with server file transfer support, namely, SecureVision. A custom-made random number generator using the equation of an ellipse was developed to randomly shuffle the pixel positions. SecureVision was found to be robust, user-friendly and fast in both encryption and decryption. The program was highly sensitive to the supplied keys, which prevented brute-force attacks. SecureVision provided full user control, where users could modify the program modules to match their desired applications, which was particularly desirable for pedagogical purposes in that interested parties had the freedom to explore the concept of image encryption and decryption. SecureVision is distributed under a GPLv3 license, which would allow everyone to use, modify and distribute the program without any restriction.

The molten core material (corium) that can be formed during severe accident scenarios in light water reactors (LWRs) can drop in the water pool and form a “melt jet” (hereafter referred to as “jet”) flow. We presented a numerical model to study the hydrodynamic deformation and fragmentation of a jet in the water pool using the Level Set (LS) method, despite that previous investigations had employed the widely used Volume-Of-Fluid (VOF) method. The effects of jet inlet speed, diameter and surface tension on the jet breakup length were studied and the results were in good agreement with previously performed experiments and with results obtained using the VOF method. The dimensionless jet breakup length was found to generally decrease with increasing jet inlet speed and increase with increasing jet diameter. In conclusion, the present model could successfully simulate the jet breakup in melt coolant interactions, and the results demonstrated a promising use of the LS method. Keywords: Fuel-coolant interactions, Level set method, Breakup length, CFD, FCI

X-ray and γ-ray photons have been widely used for studying radiobiological effects of ionizing radiations. Photons are indirectly ionizing radiations so they need to set in motion electrons (which are a directly ionizing radiation) to perform the ionizations. When the photon dose decreases to below a certain limit, the number of electrons set in motion will become so small that not all cells in an "exposed" cell population can get at least one electron hit. When some cells in a cell population are not hit by a directly ionizing radiation (in other words not irradiated), there will be rescue effect between the irradiated cells and non-irradiated cells, and the resultant radiobiological effect observed for the "exposed" cell population will be different. In the present paper, the mechanisms underlying photon interactions in radiobiological experiments were studied using our developed NRUphoton computer code, which was benchmarked against the MCNP5 code by comparing the photon dose delivered to the cell layer underneath the water medium. The following conclusions were reached: (1) The interaction fractions decreased in the following order: 16O > 12C > 14N > 1H. Bulges in the interaction fractions (versus water medium thickness) were observed, which reflected changes in the energies of the propagating photons due to traversals of different amount of water medium as well as changes in the energy-dependent photon interaction cross-sections. (2) Photoelectric interaction and incoherent scattering dominated for lower-energy (10 keV) and high-energy (100 keV and 1 MeV) incident photons. (3) The fractions of electron ejection from different nuclei were mainly governed by the photoelectric effect cross-sections, and the fractions from the 1s subshell were the largest. (4) The penetration fractions in general decreased with increasing medium thickness, and increased with increasing incident photon energy, the latter being explained by the corresponding reduction in interaction cross-sections. (5) The areas under the angular distribution curves of photons exiting the medium layer and subsequently undergoing interactions within the cell layer became smaller for larger incident photon energies. (6) The number of cells suffering at least one electron hit increased with the administered dose. For larger incident photon energies, the numbers of cells suffering at least one electron hit became smaller, which was attributed to the reduction in the photon interaction cross-section. These results highlighted the importance of the administered dose in radiobiological experiments. In particular, the threshold administered doses at which all cells in the exposed cell array suffered at least one electron hit might provide hints on explaining the intriguing observation that radiation-induced cancers can be statistically detected only above the threshold value of ~100 mSv, and thus on reconciling controversies over the linear no-threshold model.

Monte Carlo method was used to study the characteristics of neutron interactions with cells underneath a water medium layer with varying thickness. The following results were obtained. (1) The fractions of neutron interaction with 1H, 12C, 14N and 16O nuclei in the cell layer were studied. The fraction with 1H increased with increasing medium thickness, while decreased for 12C, 14N and 16O nuclei. The bulges in the interaction fractions with 12C, 14N and 16O nuclei were explained by the resonance spikes in the interaction cross-section data. The interaction fraction decreased in the order: 1H > 16O > 12C > 14N. (2) In general, as the medium thickness increased, the number of "interacting neutrons" which exited the medium and then further interacted with the cell layer increased. (3) The area under the angular distributions for "interacting neutrons" decreased with increasing incident neutron energy. Such results would be useful for deciphering the reasons behind discrepancies among existing results in the literature.

Radiotherapy is a common cancer treatment module, where a certain amount of dose will be delivered to the targeted organ. This is achieved usually by photons generated by linear accelerator units. However, radiation scattering within the patient's body and the surrounding environment will lead to dose dispersion to healthy tissues which are not targets of the primary radiation. Determination of the dispersed dose would be important for assessing the risk and biological consequences in different organs or tissues. In the present work, the concept of conversion coefficient (F) of the dispersed dose was developed, in which F = (Dd/Dt), where Dd was the dispersed dose in a non-targeted tissue and Dt is the absorbed dose in the targeted tissue. To quantify Dd and Dt, a comprehensive model was developed using the Monte Carlo N-Particle (MCNP) package to simulate the linear accelerator head, the human phantom, the treatment couch and the radiotherapy treatment room. The present work also demonstrated the feasibility and power of parallel computing through the use of the Message Passing Interface (MPI) version of MCNP5.

Cold atmospheric plasmas (CAPs) are being used in applications related to dentistry. Potential benefits include tooth whitening/bleaching, the sterilization of dental cavities, and root canal disinfection. Generated reactive species, such as hydroxyl (OH) radicals, play a critical role in the effectiveness of CAPs in dentistry. In the present work, the mandibular jaw and teeth were modeled. The propagation of CAP plume in ambient air was dynamically tracked using the level set method. The transport and dispersion OH radicals away from the nozzle and towards the teeth under treatment were also tracked. The distributions of concentration of OH radicals over the teeth were obtained for nozzle to tooth distances of 2 and 4 mm. The discharge of the OH radicals out of the nozzle was found to be asymmetrical. Interestingly, depending on the type of tooth treated, the dispersion of OH radicals out of the nozzle could be altered. The present model and obtained results could be useful for advancements towards a fear-free dentistry using CAPs.

Monte Carlo method has been used to determine the efficiency for proton production and to study the energy and angular distributions of the generated protons. The ENDF library of cross sections is used to simulate the interactions between the neutrons and the atoms in a polyethylene (PE) layer, while the ranges of protons with different energies in PE are determined using the Stopping and Range of Ions in Matter (SRIM) computer code. The efficiency of proton production increases with the PE layer thickness. However the proton escaping from a certain polyethylene volume is highly dependent on the neutron energy and target thickness, except for a very thin PE layer. The energy and angular distributions of protons are also estimated in the present paper, showing that, for the range of energy and thickness considered, the proton flux escaping is dependent on the PE layer thickness, with the presence of an optimal thickness for a fixed primary neutron energy.

Radiochromic films are convenient dosimeters for quality assurance of biological X-ray irradiators. However, their calibration curves (optical density vs. radiation dose) can have relatively short lifespans due to non-catalytic auto-development. In the present work, “corrected” calibration curves were proposed, which were constructed based on the darkness of the unexposed film before non-catalytic development. The “corrected” calibration curves enabled the use of a previous calibration curve which had expired due to non-catalytic development. Keywords: Radiochromic film, EBT3 film, Biological X-ray irradiator, Radiobiology

Commercially available radiochromic films are primarily designed for clinical X-ray dosimetry. These films change color upon exposures to radiation as a result of solid-state polymerization (SSP). Built on a previous model developed for SSP upon exposures to ultraviolet (UV) radiation, a new model was developed in the present work for X-ray-induced coloration in Gafchromic EBT3 films. Monte Carlo simulations using the Monte Carlo N-Particle (MCNP) code were employed to model the transport and interaction of photons and the generated secondary electrons within the film active layer. The films were exposed to continuous-energy photon beams. The dose DE in the external radiation detector (i.e., ionization chamber) was determined and the realistic dose DA in the film active layer was then obtained using the calibration coefficient R (=DA/DE). The finite element method (FEM) was used to solve the classical steady-state Helmholtz equation using the multifrontal massively parallel sparse direct solver (MUMPS). An extensive grid independence test was carried out and the numerical stability of the present model was ensured. The reflected light intensity from the film surface was used to theoretically obtain the net reflective optical density of the film exposed to X-ray. Good agreement was obtained between the experimental and theoretical results of the net reflective optical density of the film. For X-ray doses >~600 cGy, due to the already formed densely cross-linked structure in the active layer of the EBT3 film so further bond formation was less likely, the reflected light intensity from the film surface increased at a relatively lower rate when compared to those for dose values

Radiochromic films change color upon exposures to radiation doses as a result of solid-state polymerization (SSP). Commercially available radiochromic films are primarily designed for, and have become widely used in, clinical X-ray dosimetry. However, many intriguing properties of radiochromic films are not yet fully understood. The present work aimed at developing a theoretical model at both atomic and macroscopic scales to provide a platform for future works to understand these intriguing properties. Despite the fact that radiochromic films were primarily designed for clinical X-ray dosimetry, dose-response curves for the Gafchromic EBT3 film obtained for ultraviolet (UV) radiation were employed to develop our model in order to avoid complications of ionization, non-uniform energy deposition, as well as dispersed doses caused by secondary electrons set in motion by the indirectly ionizing X-ray photons, which might introduce added uncertainties to the model and overshadow the basic SSP processes. The active layer in the EBT3 film consisted of diacetylene (DA) pentacosa-10,12-diynoate monomers, which were modelled using molecular dynamics (MD). The degrees of SSP in the atomic scale upon different UV exposures were obtained to determine the absorption coefficients of the active layer, which were then input into the finite element method (FEM). The classical steady-state Helmholtz equation was engaged to model the reflection from the active layer using the FEM technique. The multifrontal massively parallel sparse direct solver (MUMPS) was employed to solve the present numerical problem. Very good agreement between experimentally and theoretically obtained coloration in terms of net reflective optical density was achieved for different UV exposures. In particular, for UV exposures larger than ~40 J/cm2, the reflected light intensity decreased at a lower rate when compared to other UV exposure values, which was explained by the densely cross-linked structure under near-complete polymerization, and thus the lower efficiency for further bond formation between DA monomer strands.

Cold plasmas were proposed for treatment of leukemia. In the present work, conceptual designs of mixing chambers that increased the contact between the two fluids (plasma and blood) through addition of obstacles within rectangular-block-shaped chambers were proposed and the dynamic mixing between the plasma and blood were studied using the level set method coupled with heat transfer. Enhancement of mixing between blood and plasma in the presence of obstacles was demonstrated. Continuous tracking of fluid mixing with determination of temperature distributions was enabled by the present model, which would be a useful tool for future development of cold plasma devices for treatment of blood-related diseases such as leukemia.

Previously, we have developed DynamicMC for modeling relative movement of Oak Ridge National Laboratory phantom in a radiation field for the Monte Carlo N-Particle package (Health Physics. 2023,124(4):301–309). Using this software, three-dimensional dose distributions in a phantom irradiated by a certain mono-energetic (Mono E) source can be deduced through its graphical user interface. In this study, we extended DynamicMC to be used in combination with the Particle and Heavy Ion Transport code System (PHITS) by providing it with a higher flexibility for dynamic movement for an anthropomorphic phantom. For this purpose, we implemented four new functions into the software, which are (1) to generate not only Mono E sources but also those having an energy spectrum of an arbitrary radioisotope (2) to calculate the absorbed doses for several radiologically important organs (3) to automatically average the calculated absorbed doses along the path of the phantom and (4) to generate user-defined slab shielding materials. The first and third items utilize the PHITS-specific modalities named radioisotope-source and sumtally functions, respectively. The computational cost and complexity can be dramatically reduced with these features. We anticipate that the present work and the developed open-source tools will be in the interest of nuclear radiation physics community for research and teaching purposes. [ABSTRACT FROM AUTHOR]

The Monte Carlo (MC) method is a powerful tool for modeling nuclear radiation interaction with matter. A variety of MC software packages has been developed, especially for applications in radiation therapy. Most widely used MC packages require users to write their own input scripts for their systems, which can be a time consuming and error prone process and requires extensive user experience. In the present work, we have developed a graphical user interface (GUI) bundled with a custom-made 3D OpenGL visualizer for PHITS MC package. The current version focuses on modeling proton induced positron emitting radioisotopes, which in turn can be used for verification of proton ranges in proton therapy. The developed GUI program does not require extensive user experience. The present open-source program is distributed under GPLv3 license that allows users to freely download, modify, recompile and redistribute the program. [ABSTRACT FROM AUTHOR]