Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HEART”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,657 records · Page 92

Consensus Statement on Proton Therapy in Early-Stage and Locally Advanced Non–Small Cell Lung Cancer

Radiation dose escalation has been shown to improve local control and survival in patients with non–small cell lung cancer in some studies, but randomized data have not supported this premise, possibly owing to adverse effects. Because of the physical characteristics of the Bragg peak, proton therapy (PT) delivers minimal exit dose distal to the target volume, resulting in better sparing of normal tissues in comparison to photon-based radiation therapy. This is particularly important for lung cancer given the proximity of the lung, heart, esophagus, major airways, large blood vessels, and spinal cord. However, PT is associated with more uncertainty because of the finite range of the proton beam and motion for thoracic cancers. PT is more costly than traditional photon therapy but may reduce side effects and toxicity-related hospitalization, which has its own associated cost. The cost of PT is decreasing over time because of reduced prices for the building, machine, maintenance, and overhead, as well as newer, shorter treatment programs. PT is improving rapidly as more research is performed particularly with the implementation of 4-dimensional computed tomography–based motion management and intensity modulated PT. Given these controversies, there is much debate in the oncology community about which patients with lung cancer benefit significantly from PT. The Particle Therapy Co-operative Group (PTCOG) Thoracic Subcommittee task group intends to address the issues of PT indications, advantages and limitations, cost-effectiveness, technology improvement, clinical trials, and future research directions. This consensus report can be used to guide clinical practice and indications for PT, insurance approval, and clinical or translational research directions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An Analysis of Plan Robustness for Esophageal Tumors: Comparing Volumetric Modulated Arc Therapy Plans and Spot Scanning Proton Planning

Purpose: Planning studies to compare x-ray and proton techniques and to select the most suitable technique for each patient have been hampered by the nonequivalence of several aspects of treatment planning and delivery. A fair comparison should compare similarly advanced delivery techniques from current clinical practice and also assess the robustness of each technique. The present study therefore compared volumetric modulated arc therapy (VMAT) and single-field optimization (SFO) spot scanning proton therapy plans created using a simultaneous integrated boost (SIB) for dose escalation in midesophageal cancer and analyzed the effect of setup and range uncertainties on these plans. Methods and Materials: For 21 patients, SIB plans with a physical dose prescription of 2 Gy or 2.5 Gy/fraction in 25 fractions to planning target volume (PTV){sub 50Gy} or PTV{sub 62.5Gy} (primary tumor with 0.5 cm margins) were created and evaluated for robustness to random setup errors and proton range errors. Dose–volume metrics were compared for the optimal and uncertainty plans, with P<.05 (Wilcoxon) considered significant. Results: SFO reduced the mean lung dose by 51.4% (range 35.1%-76.1%) and the mean heart dose by 40.9% (range 15.0%-57.4%) compared with VMAT. Proton plan robustness to a 3.5% range error was acceptable. For all patients, the clinical target volume D{sub 98} was 95.0% to 100.4% of the prescribed dose and gross tumor volume (GTV) D{sub 98} was 98.8% to 101%. Setup error robustness was patient anatomy dependent, and the potential minimum dose per fraction was always lower with SFO than with VMAT. The clinical target volume D{sub 98} was lower by 0.6% to 7.8% of the prescribed dose, and the GTV D{sub 98} was lower by 0.3% to 2.2% of the prescribed GTV dose. Conclusions: The SFO plans achieved significant sparing of normal tissue compared with the VMAT plans for midesophageal cancer. The target dose coverage in the SIB proton plans was less robust to random setup errors and might be unacceptable for certain patients. Robust optimization to ensure adequate target coverage of SIB proton plans might be beneficial.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Cerebrovascular Diseases in Childhood Cancer Survivors: Role of the Radiation Dose to Willis Circle Arteries

Background and Purpose: The aim of this study was to investigate the role of radiation dose received to the circle of Willis (WC) during radiation therapy (RT) and of potential dose-response modifiers on the risk of stroke after treatment of childhood cancer. Methods: We evaluated the risk factors for stroke in a cohort of 3172 5-year survivors of childhood cancer who were followed up for a median time of 26 years. Radiation doses to the WC and brain structures were estimated for each of the 2202 children who received RT. Results: Fifty-four patients experienced a confirmed stroke; 39 were ischemic. Patients not receiving RT had a stroke risk similar to that of the general population, whereas those who received RT had an 8.5-fold increased risk (95% confidence interval [CI]: 6.3-11.0). The excess of incidence of stroke increased yearly. The dose of radiation to the WC, rather than to other brain structures, was found to be the best predictor of stroke. The relative risk was 15.7 (95% CI: 4.9-50.2) for doses of 40 Gy or more. At 45 years of age, the cumulative stroke incidence was 11.3% (95% CI: 7.1%-17.7%) in patients who received 10 Gy or more to the WC, compared with 1% expected from general population data. Radiation doses received to the heart and neck also increased the risk. Surgery for childhood brain cancer was linked to hemorrhagic strokes in these patients. Conclusion: The WC should be considered as a major organ at risk during RT for childhood brain cancers. The incidence of radiation-induced ischemic stroke strongly increases with long-term follow-up.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Joint Estimation of Cardiac Toxicity and Recurrence Risks After Comprehensive Nodal Photon Versus Proton Therapy for Breast Cancer

Purpose: The study aims to perform joint estimation of the risk of recurrence caused by inadequate radiation dose coverage of lymph node targets and the risk of cardiac toxicity caused by radiation exposure to the heart. Delivered photon plans are compared with realistic proton plans, thereby providing evidence-based estimates of the heterogeneity of treatment effects in consecutive cases for the 2 radiation treatment modalities. Methods and Materials: Forty-one patients referred for postlumpectomy comprehensive nodal photon irradiation for left-sided breast cancer were included. Comparative proton plans were optimized by a spot scanning technique with single-field optimization from 2 en face beams. Cardiotoxicity risk was estimated with the model of Darby et al, and risk of recurrence following a compromise of lymph node coverage was estimated by a linear dose-response model fitted to the recurrence data from the recently published EORTC (European Organisation for Research and Treatment of Cancer) 22922/10925 and NCIC-CTG (National Cancer Institute of Canada Clinical Trials Group) MA.20 randomized controlled trials. Results: Excess absolute risk of cardiac morbidity was small with photon therapy at an attained age of 80 years, with median values of 1.0% (range, 0.2%-2.9%) and 0.5% (range, 0.03%-1.0%) with and without cardiac risk factors, respectively, but even lower with proton therapy (0.13% [range, 0.02%-0.5%] and 0.06% [range, 0.004%-0.3%], respectively). The median estimated excess absolute risk of breast cancer recurrence after 10 years was 0.10% (range, 0.0%-0.9%) with photons and 0.02% (range, 0.0%-0.07%) with protons. The association between age of the patient and benefit from proton therapy was weak, almost non-existing (Spearman rank correlations of −0.15 and −0.30 with and without cardiac risk factors, respectively). Conclusions: Modern photon therapy yields limited risk of cardiac toxicity in most patients, but proton therapy can reduce the predicted risk of cardiac toxicity by up to 2.9% and the risk of breast cancer recurrence by 0.9% in individual patients. Predicted benefit correlates weakly with age. Combined assessment of the risk from cardiac exposure and inadequate target coverage is desirable for rational consideration of competing photon and proton therapy plans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Long-Term Results of a Highly Performing Conformal Electron Therapy Technique for Chest Wall Irradiation After Mastectomy

Purpose: To evaluate locoregional control and survival after mastectomy, as well as toxicity, in patients irradiated by a previously described postmastectomy highly conformal electron beam radiation therapy technique (PMERT). Methods and Materials: We included all women irradiated by postmastectomy electron beam radiation therapy for nonmetastatic breast cancer between 2007 and 2011 in our department. Acute and late toxicities were retrospectively assessed using Common Terminology Criteria for Adverse Events version 3.0 criteria. Results: Among the 796 women included, 10.1% were triple-negative, 18.8% HER2-positive, and 24.6% received neoadjuvant chemotherapy (CT). Multifocal lesions were observed in 51.3% of women, and 64.6% had at least 1 involved lymph node (LN). Internal mammary chain, supraclavicular, infraclavicular, and axillary LNs were treated in 85.6%, 88.3%, 77.9%, and 14.9% of cases, respectively. With a median follow-up of 64 months (range, 6-102 months), 5-year locoregional recurrence–free survival and overall survival were 90% (95% confidence interval 88.1%-92.4%) and 90.9% (95% confidence interval 88.9%-93%), respectively. Early skin toxicity was scored as grade 1 in 58.5% of patients, grade 2 in 35.9%, and grade 3 in 4.5%. Concomitant CT was associated with increased grade 3 toxicity (P<.001). At long-term follow-up, 29.8% of patients presented temporary or permanent hyperpigmentation or telangiectasia or fibrosis (grade 1: 23.6%; grade 2: 5.2%; grade 3: 1%), with higher rates among smokers (P=.06); 274 patients (34.4%) underwent breast reconstruction. Only 24 patients (3%) had early esophagitis of grade 1. Only 3 patients developed ischemic heart disease: all had been treated by anthracycline-based CT with or without trastuzumab, all had been irradiated to the left chest wall and LN, and all presented numerous cardiovascular risk factors (2-4 factors). Conclusions: This study demonstrated the good efficacy of this technique in terms of locoregional control and survival, and good short-term and long-term safety. Longer follow-up is required to analyze chronic cardiac events.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Individualized Selection of Beam Angles and Treatment Isocenter in Tangential Breast Intensity Modulated Radiation Therapy

Purpose and Objective: Propose a novel method for individualized selection of beam angles and treatment isocenter in tangential breast intensity modulated radiation therapy (IMRT). Methods and Materials: For each patient, beam and isocenter selection starts with the fully automatic generation of a large database of IMRT plans (up to 847 in this study); each of these plans belongs to a unique combination of isocenter position, lateral beam angle, and medial beam angle. The imposed hard planning constraint on patient maximum dose may result in plans with unacceptable target dose delivery. Such plans are excluded from further analyses. Owing to differences in beam setup, database plans differ in mean doses to organs at risk (OARs). These mean doses are used to construct 2-dimensional graphs, showing relationships between: (1) contralateral breast dose and ipsilateral lung dose; and (2) contralateral breast dose and heart dose (analyzed only for left-sided). The graphs can be used for selection of the isocenter and beam angles with the optimal, patient-specific tradeoffs between the mean OAR doses. For 30 previously treated patients (15 left-sided and 15 right-sided tumors), graphs were generated considering only the clinically applied isocenter with 121 tangential beam angle pairs. For 20 of the 30 patients, 6 alternative isocenters were also investigated. Results: Computation time for automatic generation of 121 IMRT plans took on average 30 minutes. The generated graphs demonstrated large variations in tradeoffs between conflicting OAR objectives, depending on beam angles and patient anatomy. For patients with isocenter optimization, 847 IMRT plans were considered. Adding isocenter position optimization next to beam angle optimization had a small impact on the final plan quality. Conclusion: A method is proposed for individualized selection of beam angles in tangential breast IMRT. This may be especially important for patients with cardiac risk factors or an enhanced risk for the development of contralateral breast cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Intrinsic Bulk Quantum Oscillations in a Bulk Unconventional Insulator SmB 6

The finding of bulk quantum oscillations in the Kondo insulator SmB 6 proved a considerable surprise. Subsequent measurements of bulk quantum oscillations in other correlated insulators including YbB 12 lent support to our discovery of a class of bulk unconventional insulators that host bulk quantum oscillations. Here we perform a series of experiments to examine evidence for the intrinsic character of bulk quantum oscillations in floating zone-grown single crystals of SmB 6 that have been the subject of our quantum oscillation studies. We present results of thermodynamic, transport, and composition analysis experiments on pristine floating zone-grown single crystals of SmB 6 and compare quantum oscillations with metallic LaB 6 and elemental aluminum. These results establish the intrinsic origin of quantum oscillations from the insulating bulk of floating zone-grown SmB 6 . The similarity of the Fermi surface in insulating SmB 6 with the conduction-electron Fermi surface in metallic hexaborides is at the heart of a theoretical mystery.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Determining the impacts of venoarterial extracorporeal membrane oxygenation on cerebral oxygenation using a one-dimensional blood flow simulator

Venoarterial extracorporeal membrane oxygenation (VA-ECMO) is a mechanical system that provides rapid and short-term support for patients with cardiac failure. In many patients, pulmonary function is also impaired, resulting in poorly-oxygenated cardiac outflow competing against well-oxygenated VA-ECMO outflow, a condition known as North-South syndrome. North-South syndrome is a primary concern because of its potential to cause cerebral hypoxia, which has a critical influence on neurological complications often seen in this patient population. In order to reduce ischemic neurological complications, it is important to understand how clinical decisions regarding VA-ECMO parameters influence blood oxygenation. In our work, we studied the impacts of flow rate and cannulation site on oxygenation using a one-dimensional (1D) model to simulate blood flow. Our model was initially tested by comparing blood flow results to those observed from experimental work in VA-ECMO patients. The 1D model was combined with a two-phase flow model to simulate oxygenation. Additionally, the influence of various other clinician-tunable parameters on oxygenation in the common carotid arteries (CCAs) were tested, including, blood viscosity, cannula position within the insertion artery, heart rate, and systemic vascular resistance (SVR), as well as geometrical changes such as arterial radius and length. Our results indicated that blood oxygenation to the brain strongly depended on the cannula insertion site and the VA-ECMO flow rate with a weaker but potentially significant dependence on arterial radius. During femoral cannulation, VA-ECMO flow rates greater than ~4.9 L/min were needed to perfuse the CCAs. However, axillary and central cannulation began to perfuse the CCAs at significantly lower flow (~1 L/min). These results may help explain the incidence of cerebral hypoxia in this patient population and the common need to change cannulation strategies during treatment to address this clinical problem. While this work describes patient-averaged results, determining these relationships between VA-ECMO parameters and cerebral hypoxia is an important step towards future work to develop patient-specific models that clinicians can use to improve outcomes.

1D blood flow↗

DFT and microkinetic comparison of Pt, Pd and Rh-catalyzed ammonia oxidation

Ammonia oxidation is the heart of the Ostwald process and is important in emissions control. Catalytic behaviors are a function of conditions and are observed to vary across the platinum group metals (PGMs) Pt, Pd, Rh. Here, we combine density functional theory computations and microkinetic modeling to rationalize these dependencies. We compute reactions over model (2 1 1) and (1 1 1) surfaces of PGMs. Binding energies are similar on Pd and Pt and generally greater on Rh, while activation energies vary across all metals. Rates on (2 1 1) surfaces are greater than (1 1 1) surfaces. Additionally, the stepped Pt is most active and stepped Rh most selective to N 2 at ammonia slip conditions, while at Ostwald process conditions, stepped Pd is most active and stepped Pt most selective to NO. Degree of rate and selectivity control analysis provides insights into the reactions limiting performance of PGMs. Both activation barriers and surface coverages influence rates and selectivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A numerical approach for simulating fluid structure interaction of flexible thin shells undergoing arbitrarily large deformations in complex domains

We present a new numerical methodology for simulating fluid–structure interaction (FSI) problems involving thin flexible bodies in an incompressible fluid. The FSI algorithm uses the Dirichlet–Neumann partitioning technique. The curvilinear immersed boundary method (CURVIB) is coupled with a rotation-free finite element (FE) model for thin shells enabling the efficient simulation of FSI problems with arbitrarily large deformation. Turbulent flow problems are handled using large-eddy simulation with the dynamic Smagorinsky model in conjunction with a wall model to reconstruct boundary conditions near immersed boundaries. The CURVIB and FE solvers are coupled together on the flexible solid–fluid interfaces where the structural nodal positions, displacements, velocities and loads are calculated and exchanged between the two solvers. Loose and strong coupling FSI schemes are employed enhanced by the Aitken acceleration technique to ensure robust coupling and fast convergence especially for low mass ratio problems. The coupled CURVIB-FE-FSI method is validated by applying it to simulate two FSI problems involving thin flexible structures: 1) vortex-induced vibrations of a cantilever mounted in the wake of a square cylinder at different mass ratios and at low Reynolds number; and 2) the more challenging high Reynolds number problem involving the oscillation of an inverted elastic flag. For both cases the computed results are in excellent agreement with previous numerical simulations and/or experiential measurements. Grid convergence tests/studies are carried out for both the cantilever and inverted flag problems, which show that the CURVIB-FE-FSI method provides their convergence. Finally, the capability of the new methodology in simulations of complex cardiovascular flows is demonstrated by applying it to simulate the FSI of a tri-leaflet, prosthetic heart valve in an anatomic aorta and under physiologic pulsatile conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accurate derivative evaluation for any Grad–Shafranov solver

We present a numerical scheme that can be combined with any fixed boundary finite element based Poisson or Grad–Shafranov solver to compute the first and second partial derivatives of the solution to these equations with the same order of convergence as the solution itself. At the heart of our scheme is an efficient and accurate computation of the Dirichlet to Neumann map through the evaluation of a singular volume integral and the solution to a Fredholm integral equation of the second kind. Our numerical method is particularly useful for magnetic confinement fusion simulations, since it allows the evaluation of quantities such as the magnetic field, the parallel current density and the magnetic curvature with much higher accuracy than has been previously feasible on the affordable coarse grids that are usually implemented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adaptive spectral solution method for the Landau and Lenard-Balescu equations

In this paper, we present an adaptive spectral method for solving the Landau/Fokker-Planck equation for electron-ion systems. The heart of the algorithm is an expansion in Laguerre polynomials, which has several advantages, including automatic conservation of both energy and particles without the need for any special discretization or time-stepping schemes. One drawback of such an expansion is the O(N 3 ) memory requirement, where N is the number of polynomials used. This can impose an inconvenient limit in cases of practical interest, such as when two particle species have widely separated temperatures. The algorithm we describe here addresses this problem by periodically re-projecting the solution onto a judicious choice of new basis functions that are still Laguerre polynomials but have arguments adapted to the current physical conditions. This results in a reduction in the number of polynomials needed, at the expense of increased solution time. Because the equations are solved with little difficulty, this added time is not of much concern compared to the savings in memory. To demonstrate the algorithm, we solve several relaxation problems that could not be computed with the spectral method without re-projection. Another major advantage of this method is that it can be used for collision operators more complicated than that of the Landau equation, and we demonstrate this here by using it to solve the non-degenerate quantum Lenard-Balescu (QLB) equation for a hydrogen plasma. We conclude with some comparisons of temperature relaxation problems solved with the latter equation and the Landau equation with a Coulomb logarithm inspired by the properties of the QLB operator. We find that with this choice of Coulomb logarithm, there is little difference between using the two equations for these particular systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗

A review on advances in oxidative coupling of methane (OCM) for industrial use and prospects of CO 2 –H 2 O splitting integration

The discovery of shale gas reserves has encouraged the development of direct methods for methane conversion into valuable chemicals, offering an alternative to indirect approaches that involve an energy-intensive and intermittent syngas production step, leading to high CO 2 emissions. Amongst the direct methods, the oxidative coupling of methane (OCM) is a potential pathway to reduce CO 2 emissions and can produce commodity chemicals such as ethylene, a chemical regarded as central to the petrochemical industry. Even though OCM has been studied for over four decades, the technology still has not found commercial application. Amongst the challenges regarding industrial deployment of OCM, the most significant one is the requirement of a high ethylene yield of 30 % which is currently reported to be around 20 %. Moreover, the highly exothermic nature of the process and controlling the carbon selectivity over oxides of carbon (COx) is the heart of the problem. Numerous researchers have presented promising results in terms of catalysts, reactor designs and feeding strategies for OCM. However, due to lack of inclusiveness in the results, none of the combination of catalysts, reactors and system optimizations has been able to bring about its industrial viability. The current paper presents an extensive review of the noteworthy attempts to achieve industrial targets for OCM. Moreover, a comprehensive criteria is presented which highlights the desired end state for the industrial deployment of OCM technology. Furthermore, the criteria is based on literature survey and a comparison with industrially deployed ethylene production plants i.e., naphtha or ethane steam cracker plants. Finally, a novel integration technology is presented which includes a combination of OCM and CO 2 -H 2 O splitting in a chemical looping reactor design to enable efficient energy utilization and minimal heat losses to the environment.

CO2 Splitting↗

Snowmass2021 theory frontier white paper: Astrophysical and cosmological probes of dark matter

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring thermal state in mixed immersive virtual environments

Combining immersive virtual environment (IVE) with a controlled environment is a potential solution for analyzing human thermal experience during building design. Existing studies in this field have not adequately analyzed scenarios involving stabilized comfortable and uncomfortable temperature conditions using both thermal state votes and physiological responses, or the influence of the seasons. By combining IVE with a climate chamber, called mixed IVE (MIVE) in this study, experiments were conducted to test the hypothesis that participants' virtual experience did not significantly alter their thermal experience compared to their in-situ experience. Response variables were the control temperature distribution, the thermal state vote (at temperature steps 18.3 °C, 23.8 °C, and 29.4 °C), and physiological responses (heart rate and skin temperature). The results show that the first two response variables were not significantly different between the MIVE and in-situ settings (except for one case). Due to the heat development of the head mounted display device, the mean forehead skin temperature in the MIVE experiments was significantly higher than that in the in-situ experiments in most cases. Furthermore, such difference in skin temperature did not seem to affect general thermal state votes. In addition, significant skin temperature differences at some locations were also observed between the MIVE and in-situ settings.

42 ENGINEERING↗

A variance deconvolution estimator for efficient uncertainty quantification in Monte Carlo radiation transport applications

Monte Carlo simulations are at the heart of many high-fidelity simulations and analyses for radiation transport systems. As is the case with any complex computational model, it is important to propagate sources of input uncertainty and characterize how they affect model output. Unfortunately, uncertainty quantification (UQ) is made difficult by the stochastic variability that Monte Carlo transport solvers introduce. The standard method to avoid corrupting the UQ statistics with the transport solver noise is to increase the number of particle histories, resulting in very high computational costs. In this contribution, we propose and analyze a sampling estimator based on the law of total variance to compute UQ variance even in the presence of residual noise from Monte Carlo transport calculations. Here, we rigorously derive the statistical properties of the new variance estimator, compare its performance to that of the standard method, and demonstrate its use on neutral particle transport model problems involving both attenuation and scattering physics. We illustrate, both analytically and numerically, the estimator’s statistical performance as a function of available computational budget and the distribution of that budget between UQ samples and particle histories. We show analytically and corroborate numerically that the new estimator is unbiased, unlike the standard approach, and is more accurate and precise than the standard estimator for the same computational budget.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗