Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large deviations”

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 163 records · Page 9

Microseismic Monitoring at the Farnsworth CO 2 -EOR Field

The Farnsworth Unit in northern Texas is a field site for studying geologic carbon storage during enhanced oil recovery (EOR) using CO 2 . Microseismic monitoring is essential for risk assessment by detecting fluid leakage and fractures. We analyzed borehole microseismic data acquired during CO 2 injection and migration, including data denoising, event detection, event location, magnitude estimation, moment tensor inversion, and stress field inversion. We detected and located two shallow clusters, which occurred during increasing injection pressure. The two shallow clusters were also featured by large b values and tensile cracking moment tensors that are obtained based on a newly developed moment tensor inversion method using single-borehole data. The inverted stress fields at the two clusters showed large deviations from the regional stress field. The results provide evidence for microseismic responses to CO 2 /fluid injection and migration.

02 PETROLEUM↗

Vertical structure of magnetized accretion disks around young stars

We model the vertical structure of the magnetized accretion disks that are subject to viscous and resistive heating and irradiation by the central star. We apply our formalism to the radial structure of the magnetized accretion disks that are threaded by the poloidal magnetic field dragged during the process of star formation, which was developed by Shu and coworkers. We consider disks around low-mass protostars, T Tauri, and FU Orionis stars, as well as two levels of disk magnetization: λ{sub sys}=4 (strongly magnetized disks) and λ{sub sys}=12 (weakly magnetized disks). The rotation rates of strongly magnetized disks have large deviations from Keplerian rotation. In these models, resistive heating dominates the thermal structure for the FU Ori disk, and the T Tauri disk is very thin and cold because it is strongly compressed by magnetic pressure; it may be too thin compared with observations. Instead, in the weakly magnetized disks, rotation velocities are close to Keplerian, and resistive heating is always less than 7% of the viscous heating. In these models, the T Tauri disk has a larger aspect ratio, which is consistent with that inferred from observations. All the disks have spatially extended hot atmospheres where the irradiation flux is absorbed, although most of the mass (∼90%–95%) is in the disk midplane. With the advent of ALMA one expects direct measurements of magnetic fields and their morphology at disk scales. It will then be possible to determine the mass-to-flux ratio of magnetized accretion disks around young stars, an essential parameter for their structure and evolution. Our models contribute to the understanding of the vertical structure and emission of these disks.

79 ASTRONOMY AND ASTROPHYSICS↗

The SAMI Galaxy Survey: galaxy interactions and kinematic anomalies in Abell 119

Galaxy mergers are important events that can determine the fate of a galaxy by changing its morphology, star formation activity and mass growth. Merger systems have commonly been identified from their disturbed morphologies, and we now can employ integral field spectroscopy to detect and analyze the impact of mergers on stellar kinematics as well. We visually classified galaxy morphology using deep images (μ{sub r}=28 mag arcsec{sup −2}) taken by the Blanco 4 m telescope at the Cerro Tololo Inter-American Observatory. In this paper we investigate 63 bright (M{sub r}<−19.3) spectroscopically selected galaxies in Abell 119, of which 53 are early type and 20 show a disturbed morphology by visual inspection. A misalignment between the major axes in the photometric image and the kinematic map is conspicuous in morphologically disturbed galaxies. Our sample is dominated by early-type galaxies, yet it shows a surprisingly tight Tully–Fisher relation except for the morphologically disturbed galaxies which show large deviations. Three out of the eight slow rotators in our sample are morphologically disturbed. The morphologically disturbed galaxies are generally more asymmetric, visually as well as kinematically. Our findings suggest that galaxy interactions, including mergers and perhaps fly-bys, play an important role in determining the orientation and magnitude of a galaxy’s angular momentum.

79 ASTRONOMY AND ASTROPHYSICS↗

KMT-2019-BLG-1715: Planetary Microlensing Event with Three Lens Masses and Two Source Stars

We investigate the gravitational microlensing event KMT-2019-BLG-1715, the light curve of which shows two short-term anomalies from a caustic-crossing binary-lensing light curve: one with a large deviation and the other with a small deviation. We identify five pairs of solutions, in which the anomalies are explained by adding an extra lens or source component in addition to the base binary-lens model. We resolve the degeneracies by applying a method in which the measured flux ratio between the first and second source stars is compared with the flux ratio deduced from the ratio of the source radii. Applying this method leaves a single pair of viable solutions, in both of which the major anomaly is generated by a planetary-mass third body of the lens, and the minor anomaly is generated by a faint second source. A Bayesian analysis indicates that the lens comprises three masses: a planet-mass object with ∼2.6 M {sub J} and binary stars of K and M dwarfs lying in the galactic disk. We point out the possibility that the lens is the blend, and this can be verified by conducting high-resolution follow-up imaging for the resolution of the lens from the source.

79 ASTRONOMY AND ASTROPHYSICS↗

Preconditioning of Interplanetary Space Due to Transient CME Disturbances

Interplanetary space is characteristically structured mainly by high-speed solar wind streams emanating from coronal holes and transient disturbances such as coronal mass ejections (CMEs). While high-speed solar wind streams pose a continuous outflow, CMEs abruptly disrupt the rather steady structure, causing large deviations from the quiet solar wind conditions. For the first time, we give a quantification of the duration of disturbed conditions (preconditioning) for interplanetary space caused by CMEs. To this aim, we investigate the plasma speed component of the solar wind and the impact of in situ detected interplanetary CMEs (ICMEs), compared to different background solar wind models (ESWF, WSA, persistence model) for the time range 2011–2015. We quantify in terms of standard error measures the deviations between modeled background solar wind speed and observed solar wind speed. Using the mean absolute error, we obtain an average deviation for quiet solar activity within a range of 75.1–83.1 km s{sup −1}. Compared to this baseline level, periods within the ICME interval showed an increase of 18%–32% above the expected background, and the period of two days after the ICME displayed an increase of 9%–24%. We obtain a total duration of enhanced deviations over about three and up to six days after the ICME start, which is much longer than the average duration of an ICME disturbance itself (∼1.3 days), concluding that interplanetary space needs ∼2–5 days to recover from the impact of ICMEs. The obtained results have strong implications for studying CME propagation behavior and also for space weather forecasting.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating wind profiles in a numerical weather prediction model with Doppler lidar

We use Doppler lidar wind profiles from six locations around the globe to evaluate the wind profile forecasts in the boundary layer generated by the operational global Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts (ECMWF). The six locations selected cover a variety of surfaces with different characteristics (rural, marine, mountainous urban, and coastal urban). We first validated the Doppler lidar observations at four locations by comparison with co-located radiosonde profiles to ensure that the Doppler lidar observations were of sufficient quality. The two observation types agree well, with the mean absolute error (MAE) in wind speed almost always less than 1 m s –1 . Large deviations in the wind direction were usually only seen for low wind speeds and are due to the wind direction uncertainty increasing rapidly as the wind speed tends to zero. Time–height composites of the wind evaluation with 1 h resolution were generated, and evaluation of the model winds showed that the IFS model performs best over marine and coastal locations, where the mean absolute wind vector error was usually less than 3 m s –1 at all heights within the boundary layer. Larger errors were seen in locations where the surface was more complex, especially in the wind direction. For example, in Granada, which is near a high mountain range, the IFS model failed to capture a commonly occurring mountain breeze, which is highly dependent on the sub-grid-size terrain features that are not resolved by the model. The uncertainty in the wind forecasts increased with forecast lead time, but no increase in the bias was seen. At one location, we conditionally performed the wind evaluation based on the presence or absence of a low-level jet diagnosed from the Doppler lidar observations. The model was able to reproduce the presence of the low-level jet, but the wind speed maximum was about 2 m s –1 lower than observed. This is attributed to the effective vertical resolution of the model being too coarse to create the strong gradients in wind speed observed. Our results show that Doppler lidar is a suitable instrument for evaluating the boundary layer wind profiles in atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Coordination and control – limits in standard representations of multi-reservoir operations in hydrological modeling

Abstract. Major multi-reservoir cascades represent a primary mechanism for dealing with hydrologic variability and extremes within institutionally complex river basins worldwide. These coordinated management processes fundamentally reshape water balance dynamics. Yet, multi-reservoir coordination processes have been largely ignored in the increasingly sophisticated representations of reservoir operations within large-scale hydrological models. The aim of this paper is twofold, namely (i) to provide evidence that the common modeling practice of parameterizing each reservoir in a cascade independently from the others is a significant approximation and (ii) to demonstrate potential unintended consequences of this independence approximation when simulating the dynamics of hydrological extremes in complex reservoir cascades. We explore these questions using the Water Balance Model, which features detailed representations of the human infrastructure coupled to the natural processes that shape water balance dynamics. It is applied to the Upper Snake River basin in the western US and its heavily regulated multi-reservoir cascade. We employ a time-varying sensitivity analysis that utilizes the method of Morris factor screening to explicitly track how the dominant release rule parameters evolve both along the cascade and in time according to seasonal high- and low-flow events. This enables us to address aim (i) by demonstrating how the progressive and cumulative dominance of upstream releases significantly dampens the ability of downstream reservoir rules' parameters to influence flow conditions. We address aim (ii) by comparing simulation results with observed reservoir operations during critical low-flow and high-flow events in the basin. Our time-varying parameter sensitivity analysis with the method of Morris clarifies how independent single-reservoir parameterizations and their tacit assumption of independence leads to reservoir release behaviors that generate artificial water shortages and flooding, whereas the observed coordinated cascade operations avoided these outcomes for the same events. To further explore the role of (non-)coordination in the large deviations from the observed operations, we use an offline multi-reservoir water balance model in which adding basic coordination mechanisms drawn from the observed emergency operations is sufficient to correct the deficiencies of the independently parameterized reservoir rules from the hydrological model. These results demonstrate the importance of understanding the state–space context in which reservoir releases occur and where operational coordination plays a crucial role in avoiding or mitigating water-related extremes. Understanding how major infrastructure is coordinated and controlled in major river basins is essential for properly assessing future flood and drought hazards in a changing world.

54 ENVIRONMENTAL SCIENCES↗

Linear Theory of Electron-Plasma Waves at Arbitrary Collisionality

The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite- Laguerre decomposition of the velocity dependence of the electron distribution function. The damping rate, frequency, and eigenmode spectrum of electron-plasma waves are found as functions of the collision frequency and wavelength. A comparison is made between the collisionless Landau damping limit, the Lenard-Bernstein and Dougherty collision operators, and the electron-ion collision operator, finding large deviations in the damping rates and eigenmode spectra. A purely damped entropy mode, characteristic of a plasma where pitch-angle scattering effects are dominant with respect to collisionless effects, is shown to emerge numerically, and its dispersion relation is analytically derived. It is shown that such a mode is absent when simplified collision operators are used, and that like-particle collisions strongly influence the damping rate of the entropy mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiation Tolerance of Silicon Carbide Joints at LWR Relevant Temperature

The development of reliable joints in reactor service conditions is a critical requirement for SiC based accident tolerant fuels and core structures for light water reactor (LWR) applications. The effects of neutron irradiation on the mechanical properties is an important consideration under normal operation. Though various SiC joints retained the strength following neutron irradiation at ∼500 and ∼800 deg. C to ∼3 dpa, the irradiation tolerance of the joints at the LWR relevant temperature is still an unknown. This study presents the strengths of the candidate SiC joints following irradiation at ∼300 deg. C. The bonding strength was investigated using torsional shear testing which was developed for miniature irradiation specimens in previous work. Torsional shear strengths of various SiC joints were evaluated following irradiation at 270-310 deg. C to 8.7 dpa. Strength degradation by the irradiation likely occurred in the Mo diffusion bond. No obvious strength degradation was found in the MAX phase and hybrid polymer/CVD joints. The irradiation caused large deviation of the strength of Ti diffusion bond, TEP slurry and tape joints. The cause of this deviation will be investigated based on microstructural analysis in future work. (authors)

36 MATERIALS SCIENCE↗

How the “main condition” of phase stability can explain the effect of the velocity deviation of secondary electrons in DC-biased single-sided multipactors

In this work, a “main condition” for phase stability has been employed to investigate the effects of the velocity deviation of the electrons in DC-biased single-sided multipactors (MPs). In a previous study [M. Mostajeran, Phys. Plasmas 21, 053108 (2014)], a stability equation was derived, where the secondary electron was assumed to have zero initial velocity and the phase deviation from the resonant phase was considered. In this work, both deviations in phase and velocity from the resonant condition are taken into account, assuming nonzero initial velocity for the secondary electrons. Using the main condition for stability, it is shown that MP discharge can rise in situations, where large velocity deviations from initial velocity and large phase deviations from resonant phase exist. This is contrary to what can be predicted on the basis of the “simple stability condition.” This result is further confirmed by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rapid Crossing of a Resonance During Tune Modulation by Synchrotron Oscillations

It had become a standard practice to constrain the particle’s tune footprint while designing the ring lattice so that the particle tunes fit between harmful resonances that limit ring dynamic aperture (DA). This approach, known as “tune confinement”, puts tight limits on the magnitude of the tune shifts with amplitude and with momentum. The latter requires labor-intensive optimization of the off-energy dynamic aperture and control over the corresponding tune footprint for the large energy deviations and is a prerequisite to achieving reasonable lifetime if it is dominated by the Touschek scattering. As nonlinearities of the modern ring lattices are much enhanced as compared with the previous generation of synchrotron light sources, it is becoming more and more difficult to keep the off-momentum tune footprint confined. While one could reach reasonable on-energy DA suitable for injection into the machine, the off-energy DA optimization is much more complex and needs to satisfy more constraints so to reach large energy acceptance and secure decent lifetime at the high ring current. In certain cases when the lattice solution with the confined off-energy tune footprint cannot be found, one may ask whether crossing of a resonance stopband from off-momentum particle will lead to loss of this particle. Intuitively, if the particle crosses the stopband quickly, one may expect that the betatron oscillation amplitude will not increase substantially thereby keeping the particle within the machine acceptance. At the same time the tight tolerances of the modern lattice elements as designed and produced afford much narrower resonance stopbands being compared with the machines built decades ago. Under these circumstances the question above may be asked as what are the conditions for operating a ring with a large chromatic tune footprint, which crosses resonances, while not observing particle losses. In this paper we try to answer the questions above by studying the model of a particle with large initial momentum deviation (as, for example, driven by the Touschek scattering event), which crosses stopband of a half-integer resonance while executing synchrotron oscillations in the highly chromatic ring lattice.

36 MATERIALS SCIENCE↗

Impact of buildings on surface solar radiation over urban Beijing

The rugged surface of an urban area due to varying buildings can interact with solar beams and affect both the magnitude and spatiotemporal distribution of surface solar fluxes. Here we systematically examine the impact of buildings on downward surface solar fluxes over urban Beijing by using a 3-D radiation parameterization that accounts for 3-D building structures vs. the conventional plane-parallel scheme. We find that the resulting downward surface solar flux deviations between the 3-D and the plane-parallel schemes are generally ±1–10 W m -2 at 800 m grid resolution and within ±1 W m -2 at 4 km resolution. Pairs of positive–negative flux deviations on different sides of buildings are resolved at 800 m resolution, while they offset each other at 4 km resolution. Flux deviations from the unobstructed horizontal surface at 4 km resolution are positive around noon but negative in the early morning and late afternoon. The corresponding deviations at 800 m resolution, in contrast, show diurnal variations that are strongly dependent on the location of the grids relative to the buildings. Both the magnitude and spatiotemporal variations of flux deviations are largely dominated by the direct flux. Furthermore, we find that flux deviations can potentially be an order of magnitude larger by using a finer grid resolution. Atmospheric aerosols can reduce the magnitude of downward surface solar flux deviations by 10–65 %, while the surface albedo generally has a rather moderate impact on flux deviations. The results imply that the effect of buildings on downward surface solar fluxes may not be critically significant in mesoscale atmospheric models with a grid resolution of 4 km or coarser. However, the effect can play a crucial role in meso-urban atmospheric models as well as microscale urban dispersion models with resolutions of 1 m to 1 km.

54 ENVIRONMENTAL SCIENCES↗

Adaptive Gain-based Stable Power Smoothing of a DFIG

In a power system that has a high wind penetration, the output power fluctuation of a large-scale wind turbine generator (WTG) caused by the varying wind speed increases the maximum frequency deviation, which is an important metric to assess the quality of electricity, because of the reduced system inertia. This paper proposes a stable power-smoothing scheme of a doubly-fed induction generator (DFIG) that can suppress the maximum frequency deviation, particularly for a power system with a high wind penetration. To do this, the proposed scheme employs an additional control loop relying on the system frequency deviation that operates in combination with the maximum power point tracking control loop. To improve the power-smoothing capability while guaranteeing the stable operation of a DFIG, the gain of the additional loop is modified with the rotor speed and frequency deviation. The gain is set to be high if the rotor speed and/or frequency deviation is large. Here, the simulation results based on the IEEE 14-bus system demonstrate that the proposed scheme significantly lessens the output power fluctuation of a WTG under various scenarios by modifying the gain with the rotor speed and frequency deviation, and thereby it can regulate the frequency deviation within a narrow range.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Finite-time correlations boost large voltage angle fluctuations in electric power grids

Abstract Decarbonization in the energy sector has been accompanied by an increased penetration of new renewable energy sources in electric power systems. Such sources differ from traditional productions in that, first, they induce larger, undispatchable fluctuations in power generation and second, they lack inertia. Recent measurements have indeed reported long, non-Gaussian tails in the distribution of local voltage frequency data. Large frequency deviations may induce grid instabilities, leading in worst-case scenarios to cascading failures and large-scale blackouts. In this article, we investigate how correlated noise disturbances, characterized by the cumulants of their distribution, propagate through meshed, high-voltage power grids. For a single source of fluctuations, we show that long noise correlation times boost non-Gaussian voltage angle fluctuations so that they propagate similarly to Gaussian fluctuations over the entire network. However, they vanish faster, over short distances if the noise fluctuates rapidly. We furthermore demonstrate that a Berry–Esseen theorem leads to the vanishing of non-Gaussianities as the number of uncorrelated noise sources increases. Our predictions are corroborated by numerical simulations on realistic models of power grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SU-E-T-98: An Analysis of TG-51 Electron Beam Calibration Correction Factor Uncertainty

Purpose: To analyze the uncertainty of the TG-51 electron beam calibration correction factors for farmer type ion chambers currently used by institutions visited by IROC Houston. Methods: TG-51 calibration data were collected from 181 institutions visited by IROC Houston physicists for 1174 and 197 distinct electron beams from modern Varian and Elekta accelerators, respectively. Data collected and analyzed included ion chamber make and model, nominal energy, N{sub D,w}, I{sub 50}, R{sub 50}, k’R{sub 50}, d{sub ref}, P{sub gr} and pdd(d{sub ref}). k’R{sub 50} data for parallel plate chambers were excluded from the analysis. Results: Unlike photon beams, electron nominal energy is a poor indicator of the actual energy as evidenced by the range of R{sub 50} values for each electron beam energy (6–22MeV). The large range in R{sub 50} values resulted k’R{sub 50} values with a small standard deviation but large range between maximum value used and minimum value (0.001–0.029) used for a specific Varian nominal energy. Varian data showed more variability in k’R{sub 50} values than the Elekta data (0.001–0.014). Using the observed range of R{sub 50} values, the maximum spread in k’R{sub 50} values was determined by IROC Houston and compared to the spread of k’R{sub 50} values used in the community. For Elekta linacs the spreads were equivalent, but for Varian energies of 6 to 16MeV, the community spread was 2 to 6 times larger. Community P{sub gr} values had a much larger range of values for 6 and 9 MeV values than predicted. The range in Varian pdd(d{sub ref} ) used by the community for low energies was large, (1.4–4.9 percent), when it should have been very close to unity. Exradin, PTW Roos and PTW farmer chambers N{sub D,w} values showed the largest spread, ≥11 percent. Conclusion: While the vast majority of electron beam calibration correction factors used are accurate, there is a surprising spread in some of the values used.

07 ISOTOPE AND RADIATION SOURCES↗

Testing quantum gravity using pulsed optomechanical systems

An interesting idea, dating back to Feynman [Report from Chapel Hill Conference, edited by C. M. DeWitt and D. Rickles (1957)], argues that quantum mechanics may break down for large masses if one entertains the possibility that gravity can be “classical,” thereby leading to predictions different from conventional low-energy quantum gravity. Despite the technical difficulty in testing such deviations, a large number of experimental proposals have been put forward due to the high level of fundamental interest. Here, we consider the Schrödinger-Newton (SN) theory and the correlated worldline (CWL) theory, and show that they can be distinguished from conventional quantum mechanics, as well as each other, by performing pulsed optomechanics experiments. For CWL specifically we develop a framework resembling the commonly used “Heisenberg-picture” treatment of coupled oscillators, allowing one to perform simple calculations for such systems without delving into the deeper path-integral formalism. We find that discriminating between the theories will be very difficult until experimental control over low frequency quantum optomechanical systems is pushed much further. Furthermore, the predicted departures of SN and CWL from quantum mechanics occur at the same scale, so both alternative models could in principle be probed by a single experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗

SU-E-T-401: Evaluation of TG-43 Dose Calculation Accuracy for SAVI-Based Accelerated Partial Breast Irradiation (APBI) Via Monte Carlo Simulations

Purpose: The current standard TG-43 dose calculation method for SAVI-based Accelerated Partial Breast Irradiation (APBI) assumes an ideal geometry of infinite homogeneous water. However, in SAVI treatments, the air cavity inside the device and the short source-to-skin distance raise concerns about the dose accuracy of the TG-43 method. This study is to evaluate TG-43 dose calculation accuracy in SAVI treatments using Monte Carlo (MC) simulations. Methods: We recalculated the dose distributions of 15 APBI patients treated with SAVI devices, including five cases with a size of 6–1, five with 8−1 and five with 10−1, using our in-house developed fast MC dose package for HDR brachytherapy (gBMC). A phase-space file was used to model the Ir-192 HDR source. For each case, the patient CT was converted into a voxelized phantom and the dwell positions and times were extracted from treatment plans for MC dose calculations. Clinically relevant dosimetric parameters of the recalculated dose were compared to those computed via the TG-43 approach. Results: A systematic overestimation of doses was found for the 15 cases in TG-43 results, with D90, V150, and V200 for PTV-eval 2.8±1.8%, 2.0±2.2%, and 1.8±3.5% higher than MC results. TG-43 also overestimated the dose to skin with the maximum dose 4.4±8.4% higher on average. The relatively large standard deviation seen in the difference of maximum skin dose is partially ascribed to the statistical uncertainty of MC simulations when computing the maximum dose. It took gBMC ∼1 minute to compute dose for a SAVI plan. Conclusion: The high efficiency of our gBMC package facilitated the studies with a relatively large number of cases. An overestimation of TG-43 doses was found when using this MC package to recompute doses in SAVI cases. Clinical utilization of TG-43 dose calculation method in this scenario should be aware of this fact.

60 APPLIED LIFE SCIENCES↗