Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo 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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Development and characterization of a three-dimensional radiochromic film stack dosimeter for megavoltage photon beam dosimetry

Purpose: Three-dimensional (3D) dosimeters are particularly useful for verifying the commissioning of treatment planning and delivery systems, especially with the ever-increasing implementation of complex and conformal radiotherapy techniques such as volumetric modulated arc therapy. However, currently available 3D dosimeters require extensive experience to prepare and analyze, and are subject to large measurement uncertainties. This work aims to provide a more readily implementable 3D dosimeter with the development and characterization of a radiochromic film stack dosimeter for megavoltage photon beam dosimetry. Methods: A film stack dosimeter was developed using Gafchromic{sup ®} EBT2 films. The dosimeter consists of 22 films separated by 1 mm-thick spacers. A Virtual Water™ phantom was created that maintains the radial film alignment within a maximum uncertainty of 0.3 mm. The film stack dosimeter was characterized using simulations and measurements of 6 MV fields. The absorbed-dose energy dependence and orientation dependence of the film stack dosimeter were investigated using Monte Carlo simulations. The water equivalence of the dosimeter was determined by comparing percentage-depth-dose (PDD) profiles measured with the film stack dosimeter and simulated using Monte Carlo methods. Film stack dosimeter measurements were verified with thermoluminescent dosimeter (TLD) microcube measurements. The film stack dosimeter was also used to verify the delivery of an intensity-modulated radiation therapy (IMRT) procedure. Results: The absorbed-dose energy response of EBT2 film differs less than 1.5% between the calibration and film stack dosimeter geometries for a 6 MV spectrum. Over a series of beam angles ranging from normal incidence to parallel incidence, the overall variation in the response of the film stack dosimeter is within a range of 2.5%. Relative to the response to a normally incident beam, the film stack dosimeter exhibits a 1% under-response when the beam axis is parallel to the film planes. Measured and simulated PDD profiles agree within a root-mean-square difference of 1.3%. In-field film stack dosimeter and TLD measurements agree within 5%, and measurements in the field penumbra agree within 0.5 mm. Film stack dosimeter and TLD measurements have expanded (k = 2) overall measurement uncertainties of 6.2% and 5.8%, respectively. Film stack dosimeter measurements of an IMRT dose distribution have 98% agreement with the treatment planning system dose calculation, using gamma criteria of 3% and 2 mm. Conclusions: The film stack dosimeter is capable of high-resolution, low-uncertainty 3D dose measurements, and can be readily incorporated into an existing film dosimetry program.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Thermodynamic properties of a hard/soft-magnetic bilayer model

A model for describing the thermodynamic properties of a hard/soft-magnetic bilayer is proposed and thoroughly studied using the Monte Carlo method. Temperature dependences of the heat capacity, total magnetization, magnetizations of the hard- and soft-magnetic layers, total magnetic susceptibility, and susceptibilities of the hard- and soft-magnetic layers have been calculated by this method in the framework of the proposed model. The obtained temperature dependences of the heat capacity and magnetic susceptibility display double maxima that result from the two phase transitions that take place in the system. The influence of system dimensions on the thermodynamic properties of the model has been considered.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

SU-E-T-375: Passive Scattering to Pencil-Beam-Scanning Comparison for Medulloblastoma Proton Therapy: LET Distributions and Radiobiological Implications

Purpose: To compare the linear energy transfer (LET) distributions between passive scattering and pencil beam scanning proton radiation therapy techniques for medulloblastoma patients and study the potential radiobiological implications. Methods: A group of medulloblastoma patients, previously treated with passive scattering (PS) proton craniospinal irradiation followed by prosterior fossa or involved field boost, were selected from the patient database of our institution. Using the beam geometry and planning computed tomography (CT) image sets of the original treatment plans, pencil beam scanning (PBS) treatment plans were generated for the cranial treatment for each patient, with average beam spot size of 8mm (sigma in air at isocenter). 3-dimensional dose and LET distributions were calculated by Monte Carlo methods (TOPAS) both for the original passive scattering and new pencil beam scanning treatment plans. LET volume histograms were calculated for the target and OARs and compared for the two delivery methods. Variable RBE weighted dose distributions and volume histograms were also calculated using a variable dose and LET-based model. Results: Better dose conformity was achieved with PBS planning compared to PS, leading to increased dose coverage for the boost target area and decreased average dose to the structures adjacent to it and critical structures outside the whole brain treatment field. LET values for the target were lower for PBS plans. Elevated LET values for OARs close to the boosted target areas were noticed, due to end of range of proton beams falling inside these structures, resulting in higher RBE weighted dose for these structures compared to the clinical RBE value of 1.1. Conclusion: Transitioning from passive scattering to pencil beam scanning proton radiation treatment can be dosimetrically beneficial for medulloblastoma patients. LET–guided treatment planning could contribute to better decision making for these cases, especially for critical structures at close proximity to the boosted target area.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Atomistic Simulation of the Fission-Fragment-Induced Formation of Defects in a Uranium–Molybdenum Alloy

The formation of defects in a uranium–molybdenum alloy induced by xenon and zirconium ion irradiation with energies typical of fission fragments has been simulated within the two-temperature atomistic model and Monte Carlo method. This has allowed estimating the total number of primary radiation defects formed at the “thermal spike” stage and in collision cascades for different temperatures. The calculated threshold stopping powers of ions responsible for the formation of defects at the thermal spike stage are in good agreement with the available theoretical and experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Effect of Exchange Interaction Constants on the Magnetization Reversal in a Hard/Soft Magnetic Bilayer Model

A model is developed to study a hard/soft magnetic bilayer in an external magnetic field by the Monte Carlo method. The influence of the constants of the exchange intralayer and interlayer interactions of the hard magnetic and soft magnetic layers on the magnetization reversal in the magnetic bilayer is investigated. The formation of a spin spring is shown to be mainly determined by the exchange interaction in the soft magnetic layer. Phase diagrams are plotted for the magnetic bilayer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

The Optimal Wavelengths for Light Absorption Spectroscopy Measurements Based on Genetic Algorithm–Particle Swarm Optimization

To select the optimal wavelengths in the light extinction spectroscopy measurement, genetic algorithm–particle swarm optimization (GAPSO) based on genetic algorithm (GA) and particle swarm optimization (PSO) is adopted. The change of the optimal wavelength positions in different feature size parameters and distribution parameters is evaluated. Moreover, the Monte Carlo method based on random probability is used to identify the number of optimal wavelengths, and good inversion effects of the particle size distribution are obtained. The method proved to have the advantage of resisting noise. In order to verify the feasibility of the algorithm, spectra with bands ranging from 200 to 1000 nm are computed. Based on this, the measured data of standard particles are used to verify the algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Minibeam radiation therapy for the management of osteosarcomas: A Monte Carlo study

Purpose: Minibeam radiation therapy (MBRT) exploits the well-established tissue-sparing effect provided by the combination of submillimetric field sizes and a spatial fractionation of the dose. The aim of this work is to evaluate the feasibility and potential therapeutic gain of MBRT, in comparison with conventional radiotherapy, for osteosarcoma treatments. Methods: Monte Carlo simulations (PENELOPE/PENEASY code) were used as a method to study the dose distributions resulting from MBRT irradiations of a rat femur and a realistic human femur phantoms. As a figure of merit, peak and valley doses and peak-to-valley dose ratios (PVDR) were assessed. Conversion of absorbed dose to normalized total dose (NTD) was performed in the human case. Several field sizes and irradiation geometries were evaluated. Results: It is feasible to deliver a uniform dose distribution in the target while the healthy tissue benefits from a spatial fractionation of the dose. Very high PVDR values (⩾20) were achieved in the entrance beam path in the rat case. PVDR values ranged from 2 to 9 in the human phantom. NTD{sub 2.0} of 87 Gy might be reached in the tumor in the human femur while the healthy tissues might receive valley NTD{sub 2.0} lower than 20 Gy. The doses in the tumor and healthy tissues might be significantly higher and lower than the ones commonly delivered used in conventional radiotherapy. Conclusions: The obtained dose distributions indicate that a gain in normal tissue sparing might be expected. This would allow the use of higher (and potentially curative) doses in the tumor. Biological experiments are warranted.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Laser pulse driven control of charge and spin order in the two-dimensional Kondo lattice

Fast and dynamical control of quantum phases is highly desired in the application of quantum phenomena in devices. On the experimental side, ultrafast laser pulses provide an ideal platform to induce femtosecond dynamics in a variety of materials. Here we show that a laser pulse driven heavy fermion system can be tuned to dynamically evolve into a phase, which is not present in equilibrium. Using the state-of-the-art time-dependent variational Monte Carlo method, we perform numerical simulations of realistic laser pulses applied to the Kondo lattice model. By tracking spin and charge degrees of freedom we identify a dynamical phase transition, within the partially Kondo screened phase, from a charge ordered into a metallic state, while preserving the spin order of the system. Furthermore, we propose using high-harmonic generation and nonequilibrium optical conductivity to identify the dynamical phase transition in an ultrafast experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2017

First Simultaneous Extraction of Spin-Dependent Parton Distributions and Fragmentation Functions from a Global QCD Analysis

In this paper, we perform the first global QCD analysis of polarized inclusive and semi-inclusive deep-inelastic scattering and single-inclusive $e^+e^-$ annihilation data, simultaneously fitting the parton distribution and fragmentation functions using the iterative Monte Carlo method. Without imposing SU(3) symmetry relations, we find the strange polarization to be very small, consistent with zero for both inclusive and semi-inclusive data, which provides a resolution to the strange quark polarization puzzle. Finally, the combined analysis also allows the direct extraction from data of the isovector and octet axial charges, and is consistent with a small SU(2) flavor asymmetry in the polarized sea.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

SU-F-T-610: Comparison of Output Factors for Small Radiation Fields Used in SBRT Treatment

Purpose: In order to fundamentally understand our previous dose verification results between measurements and calculations from treatment planning system (TPS) for SBRT plans for different sized targets, the goal of the present work was to compare output factors for small fields measured using EDR2 films with TPS and Monet Carlo (MC) simulations. Methods: 6MV beam was delivered to EDR2 films for each of the following field sizes; 1×1 cm{sup 2}, 1.5×1.5 cm{sup 2}, 2×2 cm{sup 2}, 3×3 cm{sup 2}, 4×4 cm{sup 2}, 5×5 cm{sup 2} and 10×10 cm{sup 2}. The films were developed in a film processer, then scanned with a Vidar VXR-16 scanner and analyzed using RIT113 version 6.1. A standard calibration curve was obtained with the 6MV beam and was used to get absolute dose for measured field sizes. Similar plans for all fields sizes mentioned above were generated using Eclipse with the Analytical Anisotropic Algorithm. Similarly, MC simulations were carried out using the MCSIM, an in-house MC code for different field sizes. Output factors normalized to 10×10 cm{sup 2} reference field were calculated for different field sizes in all the three cases and compared. Results: For field sizes ranging from 1×1 cm{sup 2} to 2×2 cm{sup 2}, the differences in output factors between measurements (films), TPS and MC simulations were within 0.22%. For field sizes ranging from 3×3cm{sup 2} to 5×5cm{sup 2}, differences in output factors were within 0.10%. Conclusion: No clinically significant difference was obtained in output factors for different field sizes acquired from films, TPS and MC simulations. Our results showed that the output factors are predicted accurately from TPS when compared to the actual measurements and superior dose calculation Monte Carlo method. This study would help us in understanding our previously obtained dose verification results for small fields used in the SBRT treatment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Residual entropy of ice III from Monte Carlo simulation

We calculated the residual entropy of ice III as a function of the occupation probabilities of hydrogen positions α and β assuming equal energies of all configurations. To do this, a discrete ice model with Bjerrum defect energy penalty and harmonic terms to constrain the occupation probabilities was simulated by the Metropolis Monte Carlo method for a range of temperatures and sizes followed by thermodynamic integration and extrapolation to N = ∞. Similarly as for other ices, the residual entropies are slightly higher than the mean-field (no-loop) approximation. However, the corrections caused by fluctuation of energies of ice samples calculated using molecular models of water are too large for accurate determination of the chemical potential and phase equilibria.

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

Poster - 16: Time-resolved diode dosimetry for in vivo proton therapy range verification: calibration through numerical modeling

Purpose: A method to refine the implementation of an in vivo, adaptive proton therapy range verification methodology was investigated. Simulation experiments and in-phantom measurements were compared to validate the calibration procedure of a time-resolved diode dosimetry technique. Methods: A silicon diode array system has been developed and experimentally tested in phantom for passively scattered proton beam range verification by correlating properties of the detector signal to the water equivalent path length (WEPL). The implementation of this system requires a set of calibration measurements to establish a beam-specific diode response to WEPL fit for the selected ‘scout’ beam in a solid water phantom. This process is both tedious, as it necessitates a separate set of measurements for every ‘scout’ beam that may be appropriate to the clinical case, as well as inconvenient due to limited access to the clinical beamline. The diode response to WEPL relationship for a given ‘scout’ beam may be determined within a simulation environment, facilitating the applicability of this dosimetry technique. Measurements for three ‘scout’ beams were compared against simulated detector response with Monte Carlo methods using the Tool for Particle Simulation (TOPAS). Results: Detector response in water equivalent plastic was successfully validated against simulation for spread out Bragg peaks of range 10 cm, 15 cm, and 21 cm (168 MeV, 177 MeV, and 210 MeV) with adjusted R{sup 2} of 0.998. Conclusion: Feasibility has been shown for performing calibration of detector response for a given ‘scout’ beam through simulation for the time resolved diode dosimetry technique.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Hydrodynamic fluctuations near a Hopf bifurcation: Stochastic onset of vortex shedding behind a circular cylinder

Here, we investigate hydrodynamic fluctuations in the flow past a circular cylinder near the critical Reynolds number Re c for the onset of vortex shedding. Starting from the fluctuating Navier-Stokes equations, we perform a perturbation expansion around Re c to derive analytical expressions for the statistics of the fluctuating lift force. Molecular-level simulations using the direct simulation Monte Carlo method support the theoretical predictions of the lift power spectrum and amplitude distribution. Notably, we have been able to collect sufficient statistics at distances Re ⁡/ Re c – 1 = O ⁡(10 –3 ) from the instability that confirm the appearance of non-Gaussian fluctuations, and we observe that they are associated with intermittent vortex shedding. These results emphasize how unavoidable thermal-noise-induced fluctuations become dramatically amplified in the vicinity of oscillatory flow instabilities and that their onset is fundamentally stochastic.

42 ENGINEERING↗
DOE OSTI2015

A signed particle formulation of non-relativistic quantum mechanics

A formulation of non-relativistic quantum mechanics in terms of Newtonian particles is presented in the shape of a set of three postulates. In this new theory, quantum systems are described by ensembles of signed particles which behave as field-less classical objects which carry a negative or positive sign and interact with an external potential by means of creation and annihilation events only. This approach is shown to be a generalization of the signed particle Wigner Monte Carlo method which reconstructs the time-dependent Wigner quasi-distribution function of a system and, therefore, the corresponding Schrödinger time-dependent wave-function. Its classical limit is discussed and a physical interpretation, based on experimental evidences coming from quantum tomography, is suggested. Moreover, in order to show the advantages brought by this novel formulation, a straightforward extension to relativistic effects is discussed. To conclude, quantum tunnelling numerical experiments are performed to show the validity of the suggested approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗
DOE OSTI2014

SU-C-BRE-05: PTV Margin Determination Based On Tumor Radiobiological Characteristics and Geometric Uncertainties Derived From Daily Cone- Beam CT Images

Purpose: To determine required PTV margins for ≤1% loss in mean population TCP using systematic (Σ) and random (σ) errors calculated from daily cone-beam CT (CBCT) images of head and neck patients. Methods: Daily CBCT images were acquired for 50 head and neck patients. The CBCT image sets acquired at each fraction were registered with planning CT to obtain positional errors for each patient for each fraction. Systematic and random errors were calculated from data collected for 50 patients as described in IPEM On Target report. CTV delineation uncertainty of 2mm is added quadratically to systematic error. Assuming a spherical target volume, the dose in each voxel of target volume is summed for each fraction in the treatment by shifting the dose grid to calculate mean population TCP inclusive of geometric uncertainties using a Monte Carlo method. These simulations were repeated for the set of Σ and σ in each axis for different PTV margins and drop in TCP for each margin are obtained. In order to study the effect of dose-response curve on PTV margins, two different σα of 0.048 Gy-1 and 0.218 Gy-1 representing steep and shallow dose-response curves are studied. Σ were 2.5, 2.5, 2.1 mm and σ were 0.3, 0.3 0.2 mm respectively in x, y and z axis respectively. Results: PTV margins based on tumor radiobiological characteristics are 4.8, 4.8 and 4 mm in x, y and z axis assuming 25 treatment fractions for σα 0.048 Gy-1 (steep) and 4.2,4.2 and 2.2 for σα of 0.218 Gy-1 (shallow). While the TCP-based margins did not differ much in x and y axis, it is considerably smaller in z axis for shallow DRC. Conclusion: TCP based margins are substantially smaller than physical dose-based margin recipes. This study also demonstrates the importance of considering tumor radiobiological characteristics while deriving margins.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Phase transitions in the antiferromagnetic Ising model on a body-centered cubic lattice with interactions between next-to-nearest neighbors

Phase transitions in the antiferromagnetic Ising model on a body-centered cubic lattice are studied on the basis of the replica algorithm by the Monte Carlo method and histogram analysis taking into account the interaction of next-to-nearest neighbors. The phase diagram of the dependence of the critical temperature on the intensity of interaction of the next-to-nearest neighbors is constructed. It is found that a second-order phase transition is realized in this model in the investigated interval of the intensities of interaction of next-to-nearest neighbors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Incorporating Experimental Information in the Total Monte Carlo Methodology Using File Weights

Some criticism has been directed towards the Total Monte Carlo method because experimental information has not been taken into account in a statistically well-founded manner. In this work, a Bayesian calibration method is implemented by assigning weights to the random nuclear data files and the method is illustratively applied to a few applications. In some considered cases, the estimated nuclear data uncertainties are significantly reduced and the central values are significantly shifted. The study suggests that the method can be applied both to estimate uncertainties in a more justified way and in the search for better central values. Some improvements are however necessary; for example, the treatment of outliers and cross-experimental correlations should be more rigorous and random files that are intended to be prior files should be generated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗