Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

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,279 records · Page 71

Effectiveness of denoising diffusion probabilistic models for fast and high-fidelity whole-event simulation in high-energy heavy-ion experiments

Artificial intelligence (AI) generative models, such as generative adversarial networks (GANs), variational autoencoders, and normalizing flows, have been widely used and studied as efficient alternatives for traditional scientific simulations. However, they have several drawbacks, including training instability and inability to cover the entire data distribution, especially for regions where data are rare. This is particularly challenging for whole-event, full-detector simulations in high-energy heavy-ion experiments, such as sPHENIX at the Relativistic Heavy Ion Collider and Large Hadron Collider experiments, where thousands of particles are produced per event and interact with the detector. This work investigates the effectiveness of denoising diffusion probabilistic models (DDPMs) as an AI-based generative surrogate model for the sPHENIX experiment that includes the heavy-ion event generation and response of the entire calorimeter stack. DDPM performance in sPHENIX simulation data is compared with a popular rival, GANs. Results show that both DDPMs and GANs can reproduce the data distribution where the examples are abundant (low-to-medium calorimeter energies). Nonetheless, DDPMs significantly outperform GANs, especially in high-energy regions where data are rare. Additionally, DDPMs exhibit superior stability compared to GANs. The results are consistent between both central and peripheral centrality heavy-ion collision events. Moreover, DDPMs offer a substantial speedup of approximately a factor of 100 compared to the traditional Geant4 simulation method.

42 ENGINEERING↗

Diffusion in Jammed Particle Packs

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions Φ in the vicinity of the jamming transition at Φ c . Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal “Fickian” diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to Φ c . The time required to recover normal diffusion t * scales as (Φ–Φ c ) –0.5 and the long-time diffusivity D∞~(Φ–Φ c ) 0.5 , or D ∞ ~1/ t *. Here, it is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both t * and D∞ in the limit Φ→Φ c .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

4D tumor centroid tracking using orthogonal 2D dynamic MRI: Implications for radiotherapy planning

Purpose: Current pretreatment, 4D imaging techniques are suboptimal in that they sample breathing motion over a very limited “snapshot” in time. Heretofore, long-duration, 4D motion characterization for radiotherapy planning, margin optimization, and validation have been impractical for safety reasons, requiring invasive markers imaged under x-ray fluoroscopy. To characterize 3D tumor motion and associated variability over durations more consistent with treatments, the authors have developed a practical dynamic MRI (dMRI) technique employing two orthogonal planes acquired in a continuous, interleaved fashion.Methods: 2D balanced steady-state free precession MRI was acquired continuously over 9–14 min at approximately 4 Hz in three healthy volunteers using a commercial 1.5 T system; alternating orthogonal imaging planes (sagittal, coronal, sagittal, etc.) were employed. The 2D in-plane pixel resolution was 2 × 2 mm{sup 2} with a 5 mm slice profile. Simultaneous with image acquisition, the authors monitored a 1D surrogate respiratory signal using a device available with the MRI system. 2D template matching-based anatomic feature registration, or tracking, was performed independently in each orientation. 4D feature tracking at the raw frame rate was derived using spline interpolation.Results: Tracking vascular features in the lung for two volunteers and pancreatic features in one volunteer, the authors have successfully demonstrated this method. Registration error, defined here as the difference between the sagittal and coronal tracking result in the SI direction, ranged from 0.7 to 1.6 mm (1σ) which was less than the acquired image resolution. Although the healthy volunteers were instructed to relax and breathe normally, significantly variable respiration was observed. To demonstrate potential applications of this technique, the authors subsequently explored the intrafraction stability of hypothetical tumoral internal target volumes and 3D spatial probability distribution functions. The surrogate respiratory information allowed the authors to show how this technique can be used to study correlations between internal and external (surrogate) information over these prolonged durations. However, compared against the gold standard of the time stamps in the dMRI frames, the temporal synchronization of the surrogate 1D respiratory information was shown to be likely unreliable.Conclusions: The authors have established viability of a novel and practical pretreatment, 4D tumor centroid tracking method employing a commercially available dynamic MRI sequence. Further developments from the vendor are likely needed to provide a reliably synchronized surrogate 1D respiratory signal, which will likely broaden the utility of this method in the pretreatment radiotherapy planning context.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-T-258: Development of a New Patient Set-Up Monitoring System Using Force Sensing Resistor (FSR) Sensor for the Radiation Therapy

Purpose: The purpose of this work is to develop a new patient set-up monitoring system using force sensing resistor (FSR) sensors that can confirm pressure of contact surface and evaluate its feasibility. Methods: In this study, we focused on develop the patient set-up monitoring system to compensate for the limitation of existing optical based monitoring system, so the developed system can inform motion in the radiation therapy. The set-up monitoring system was designed consisting of sensor units (FSR sensor), signal conditioning devices (USB cable/interface electronics), a control PC, and a developed analysis software. The sensor unit was made by attaching FSR sensor and dispersing pressure sponge to prevent error which is caused by concentrating specific point. Measured signal from the FSR sensor was sampled to arduino mega 2560 microcontroller, transferred to control PC by using serial communication. The measured data went through normalization process. The normalized data was displayed through the developed graphic user interface (GUI) software. The software was designed to display a single sensor unit intensity (maximum 16 sensors) and display 2D pressure distribution (using 16 sensors) according to the purpose. Results: Changes of pressure value according to motion was confirmed by the developed set-up monitoring system. Very small movement such as little physical change in appearance can be confirmed using a single unit and using 2D pressure distribution. Also, the set-up monitoring system can observe in real time. Conclusion: In this study, we developed the new set-up monitoring system using FSR sensor. Especially, we expect that the new set-up monitoring system is suitable for motion monitoring of blind area that is hard to confirm existing optical system and compensate existing optical based monitoring system. As a further study, an integrated system will be constructed through correlation of existing optical monitoring system. This work was supported by the Industrial R&D program of MOTIE/KEIT. [10048997, Development of the core technology for integrated therapy devices based on real-time MRI guided tumor tracking] and the Mid-career Researcher Program (2014R1A2A1A10050270) through the National Research Foundation of Korea funded by the Ministry of Science, ICT&Future Planning.

60 APPLIED LIFE SCIENCES↗

SU-G-BRC-02: A Novel Multi-Criteria Optimization Approach to Generate Deliverable Intensity-Modulated Radiation Therapy (IMRT) Treatment Plans

Purpose: To present a novel multi-criteria optimization (MCO) solution approach that generates treatment plans with deliverable apertures using column generation. Methods: We demonstrate our method with 10 locally advanced head-and-neck cancer cases retrospectively. In our MCO formulation, we defined an objective function for each structure in the treatment volume. This resulted in 9 objective functions, including 3 distinct objectives for primary target volume, high-risk and low-risk target volumes, 5 objectives for each of the organs-at-risk (OARs) (two parotid glands, spinal cord, brain stem and oral cavity), and one for the non-target non-OAR normal tissue. Conditional value-at-risk (CVaR) constraints were utilized to ensure at least certain fraction of the target volumes receiving the prescription doses. To directly generate deliverable plans, column generation algorithm was embedded within our MCO approach for aperture shape generation. Final dose distributions for all plans were generated using a Monte Carlo kernel-superposition dose calculation. We compared the MCO plans with the clinical plans, which were created by clinicians. Results: At least 95% target coverage was achieved by both MCO plans and clinical plans. However, the average conformity indices of clinical plans and the MCO plans were 1.95 and 1.35, respectively (31% reduction, p<0.01). Compared to the conventional clinical plan, the proposed MCO method achieved average reductions in left parotid mean dose of 5% (p=0.06), right parotid mean dose of 18% (p<0.01), oral cavity mean dose of 21% (p=0.03), spinal cord maximum dose of 20% (p<0.01), brain stem maximum dose of 61% (p<0.01), and normal tissue maximum dose of 5% (p<0.01), respectively. Conclusion: We demonstrated that the proposed MCO method was able to obtain deliverable IMRT treatment plans while achieving significant improvements in dosimetric plan quality.

60 APPLIED LIFE SCIENCES↗

High Value Opportunities to Advance Automation in Electric Grid Control Rooms

Grid operators and electric utilities are increasingly driven to develop automated processes and introduce decision support tools to reduce cognitive demand on System Operations staff and aid in the operational decision-making process. However, System Operation staff are traditionally hesitant to implement and accept new solutions, automation processes, and tools within the control room. This paper explores high-value opportunities to advance automation in electric grid control rooms for the purpose of improving grid safety, reliability, and resiliency during normal and emergency operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

THE RUNAWAYS AND ISOLATED O-TYPE STAR SPECTROSCOPIC SURVEY OF THE SMC (RIOTS4)

We present the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), a spatially complete survey of uniformly selected field OB stars that covers the entire star-forming body of the Small Magellanic Cloud (SMC). Using the IMACS (Inamori-Magellan Areal Camera and Spectrograph) multislit spectrograph and MIKE (Magellan Inamori Kyocera Echelle) echelle spectrograph on the Magellan telescopes, we obtained spectra of 374 early-type field stars that are at least 28 pc from any other OB candidates. We also obtained spectra of an additional 23 field stars in the SMC bar identified from slightly different photometric criteria. Here, we present the observational catalog of stars in the RIOTS4 survey, including spectral classifications and radial velocities. For three multi-slit fields covering 8% of our sample, we carried out monitoring observations over 9–16 epochs to study binarity, finding a spectroscopic, massive binary frequency of at least ∼60% in this subsample. Classical Oe/Be stars represent a large fraction of RIOTS4 (42%), occurring at much higher frequency than in the Galaxy, consistent with expectation at low metallicity. RIOTS4 confirmed a steep upper initial mass function in the field, apparently caused by the inability of the most massive stars to form in the smallest clusters. Our survey also yields evidence for in situ field OB star formation, and properties of field emission-line star populations, including sgB[e] stars and classical Oe/Be stars. We also discuss the radial velocity distribution and its relation to SMC kinematics and runaway stars. RIOTS4 presents a first quantitative characterization of field OB stars in an external galaxy, including the contributions of sparse, but normal, star formation; runaway stars; and candidate isolated star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

GAMMA–GAMMA ABSORPTION IN THE BROAD LINE REGION RADIATION FIELDS OF GAMMA-RAY BLAZARS

The expected level of γγ absorption in the Broad Line Region (BLR) radiation field of γ -ray loud Flat Spectrum Radio Quasars (FSRQs) is evaluated as a function of the location of the γ -ray emission region. This is done self-consistently with parameters inferred from the shape of the spectral energy distribution (SED) in a single-zone leptonic EC-BLR model scenario. We take into account all geometrical effects both in the calculation of the γγ opacity and the normalization of the BLR radiation energy density. As specific examples, we study the FSRQs 3C279 and PKS 1510-089, keeping the BLR radiation energy density at the location of the emission region fixed at the values inferred from the SED. We confirm previous findings that the optical depth due to γγ absorption in the BLR radiation field exceeds unity for both 3C279 and PKS 1510-089 for locations of the γ -ray emission region inside the inner boundary of the BLR. It decreases monotonically, with distance from the central engine and drops below unity for locations within the BLR. For locations outside the BLR, the BLR radiation energy density required for the production of GeV γ -rays rapidly increases beyond observational constraints, thus making the EC-BLR mechanism implausible. Therefore, in order to avoid significant γγ absorption by the BLR radiation field, the γ -ray emission region must therefore be located near the outer boundary of the BLR.

79 ASTRONOMY AND ASTROPHYSICS↗

The Environment of the Binary Neutron Star Merger GW170817

We present Hubble Space Telescope (HST) and Chandra imaging, combined with Very Large Telescope MUSE integral field spectroscopy of the counterpart and host galaxy of the first binary neutron star merger detected via gravitational-wave emission by LIGO and Virgo, GW170817. The host galaxy, NGC 4993, is an S0 galaxy at z = 0.009783. There is evidence for large, face-on spiral shells in continuum imaging, and edge-on spiral features visible in nebular emission lines. This suggests that NGC 4993 has undergone a relatively recent (≲1 Gyr) “dry” merger. This merger may provide the fuel for a weak active nucleus seen in Chandra imaging. At the location of the counterpart, HST imaging implies there is no globular or young stellar cluster, with a limit of a few thousand solar masses for any young system. The population in the vicinity is predominantly old with ≲1% of any light arising from a population with ages <500 Myr. Both the host galaxy properties and those of the transient location are consistent with the distributions seen for short-duration gamma-ray bursts, although the source position lies well within the effective radius (r{sub e}∼3 kpc), providing an r {sub e}-normalized offset that is closer than ∼90% of short GRBs. For the long delay time implied by the stellar population, this suggests that the kick velocity was significantly less than the galaxy escape velocity. We do not see any narrow host galaxy interstellar medium features within the counterpart spectrum, implying low extinction, and that the binary may lie in front of the bulk of the host galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimation and uncertainty quantification for the output from quantum simulators

The problem of estimating certain distributions over {0, 1} d is considered here. The distribution represents a quantum system of d qubits, where there are non-trivial dependencies between the qubits. A maximum entropy approach is adopted to reconstruct the distribution from exact moments or observed empirical moments. The Robbins Monro algorithm is used to solve the intractable maximum entropy problem, by constructing an unbiased estimator of the un-normalized target with a sequential Monte Carlo sampler at each iteration. In the case of empirical moments, this coincides with a maximum likelihood estimator. A Bayesian formulation is also considered in order to quantify uncertainty a posteriori. Several approaches are proposed in order to tackle this challenging problem, according to recently developed methodologies. In particular, unbiased estimators of the gradient of the log posterior are constructed and used within a provably convergent Langevin-based Markov chain Monte Carlo method. The methods are illustrated on classically simulated output from quantum simulators.

97 MATHEMATICS AND COMPUTING↗

SU-F-T-396: Impact of Shoulder Deformation for Head and Neck VMAT

Purpose: For head and neck VMAT (HN-VMAT), variations of position and deformation of patient’s shoulders is a concern to affect inaccuracy of dose distribution. It has been reported that the setup error of the shoulders was variable from 5 mm – 1 cm. The beams of the HN-VMAT pass through the shoulders. We assessed the impact of shoulder deformation to dose distribution for HN-VMAT. Methods: One HN-VMAT plan was generated using a patient’s CT. The patient’s CT was deformed using ImSimQA (Oncology Systems Limited, Shrewsbury, Shropshire, UK) to generate several patterns of the shoulders’ deformations when the right and left humeral heads were shifted with 3, 6, and 15 mm in the superior and inferior directions (SI), 3, 5, and 15 mm in the anterior and posterior directions (AP), and 5 and 15 mm in the right or left direction (LR). DVH comparison was performed in the different deformation patterns. The dosimetric parameters of D95% for CTV70Gy, CTV60Gy and CTV54Gy and dmax for Spinal cord were also measured. Gamma index evaluation (Criteria: 3%/2mm) was performed to exhibit clinically tolerable area in the comparison. Results: DVH comparison shows similar for all structures. As the comparison for the dosimetric parameters, the variations of D95% in the LR and AP were within 1%. There were larger variations in the SI than those in the other directions, however were within 1.5%. In gamma index evaluation, the small spots with higher gamma index values were appeared when the shift was 6 mm, however the pass ratio was 99.13%. Conclusion: HN-VMAT should be robust for shoulder deformation and geometric accuracy within 6 mm from patient’s setup and image-guided radiotherapy may be clinically acceptable for target dose coverage or normal tissue dose sparing.

60 APPLIED LIFE SCIENCES↗

Revealing Decision Conservativeness Through Inverse Distributionally Robust Optimization

This paper introduces Inverse Distributionally Robust Optimization (I-DRO) as a method to infer the conservativeness level of a decision-maker, represented by the size of a Wasserstein metric-based ambiguity set, from the optimal decisions made using Forward Distributionally Robust Optimization (F-DRO). By leveraging the Karush-Kuhn-Tucker (KKT) conditions of the convex F-DRO model, we formulate I-DRO as a bi-linear program, which can be solved using off-the-shelf optimization solvers. Additionally, this formulation exhibits several advantageous properties. We demonstrate that I-DRO not only guarantees the existence and uniqueness of an optimal solution but also establishes the necessary and sufficient conditions for this optimal solution to accurately match the actual conservativeness level in F-DRO. Furthermore, we identify three extreme scenarios that may impact I-DRO effectiveness. Our case study applies F-DRO for power system scheduling under uncertainty and employs I-DRO to recover the conservativeness level of system operators. Numerical experiments based on an IEEE 5-bus system and a realistic NYISO 11-zone system demonstrate I-DRO performance in both normal and extreme scenarios. An extended version of this paper with additional analyses is available at li2024revealing.

distributionally robust optimization↗

FORMATION AND EVOLUTION OF BLUE STRAGGLERS IN 47 TUCANAE

Blue stragglers (BSS) are stars whose position in the color–magnitude diagram (CMD) places them above the main sequence (MS) turn-off (TO) point of a star cluster. Using data from the core of 47 Tuc in the ultraviolet (UV), we have identified various stellar populations in the CMD, and used their radial distributions to study the evolution and origin of BSS, and obtain a dynamical estimate of the mass of BSS systems. When we separate the BSS into two samples by their magnitude, we find that the bright BSS show a much more centrally concentrated radial distribution and thus higher mass estimate (over twice the TO mass for these BSS systems), suggesting an origin involving triple or multiple stellar systems. In contrast, the faint BSS are less concentrated, with a radial distribution similar to the MS binaries, pointing to the MS binaries as the likely progenitors of these BSS. Putting our data together with available photometric data in the visible and using MESA evolutionary models, we calculate the expected number of stars in each evolutionary stage for the normal evolution of stars and the number of stars coming from the evolution of BSS. The results indicate that BSS have a post-MS evolution comparable to that of a normal star of the same mass and a MS BSS lifetime of about 200–300 Myr. We also find that the excess population of asymptotic giant branch stars in 47 Tuc is due to evolved BSS.

79 ASTRONOMY AND ASTROPHYSICS↗

Sci—Thur AM: YIS - 04: Gold Nanoparticle Enhanced Arc Radiotherapy: A Monte Carlo Feasibility Study

Introduction: The use of gold nanoparticles (GNPs) in radiotherapy has shown promise for therapeutic enhancement. In this study, we explore the feasibility of enhancing radiotherapy with GNPs in an arc-therapy context. We use Monte Carlo simulations to quantify the macroscopic dose-enhancement ratio (DER) and tumour to normal tissue ratio (TNTR) as functions of photon energy over various tumour and body geometries. Methods: GNP-enhanced arc radiotherapy (GEART) was simulated using the PENELOPE Monte Carlo code and penEasy main program. We simulated 360° arc-therapy with monoenergetic photon energies 50 – 1000 keV and several clinical spectra used to treat a spherical tumour containing uniformly distributed GNPs in a cylindrical tissue phantom. Various geometries were used to simulate different tumour sizes and depths. Voxel dose was used to calculate DERs and TNTRs. Inhomogeneity effects were examined through skull dose in brain tumour treatment simulations. Results: Below 100 keV, DERs greater than 2.0 were observed. Compared to 6 MV, tumour dose at low energies was more conformai, with lower normal tissue dose and higher TNTRs. Both the DER and TNTR increased with increasing cylinder radius and decreasing tumour radius. The inclusion of bone showed excellent tumour conformality at low energies, though with an increase in skull dose (40% of tumour dose with 100 keV compared to 25% with 6 MV). Conclusions: Even in the presence of inhomogeneities, our results show promise for the treatment of deep-seated tumours with low-energy GEART, with greater tumour dose conformality and lower normal tissue dose than 6 MV.

60 APPLIED LIFE SCIENCES↗

Smart Methane Emission Detection System Development (Final Report)

Working with the Department of Energy's National Energy Technology Laboratory, Southwest Research Institute® (SwRI®) developed a system to identify methane leaks reliably, accurately, and autonomously at critical midstream sections of the natural gas distribution network in real-time for the purpose of mitigating methane emissions using Optical Gas Imaging (OGI) cameras. SwRI's Smart Leak Detection – Methane (SLED/M) adds a high degree of automation to the process of methane leak detection to minimize sources of human error, minimize response time to a leak event, and maximize midstream visibility. Furthermore, SwRI has been working towards integrating Quantitative OGI (QOGI) capabilities into this existing technology. By leveraging Deep Learning, SwRI now has the capability to estimate fugitive emission leak rates quickly and reliably, which allows operators to detect emissions, quantify leak rate, prioritize repairs, and validate the repairs in a single instrument. The next generation QOGI technology leverages the same cameras used in Leak Detection and Repair (LDAR) programs, with improvements in safety and speed for traditional quantification-based repairs, ultimately leading to less overhead cost for the operators. The goals for this research were to develop two types of models with the following goals: Run in real-time on the edge (≥ 12 Hz), Classification: Achieve less than 5% false positive detection, Classification: Achieve ≥ 95% methane plume detection rate, Regression: achieve ≤ 10 standard cubic feet per hour (scfh) prediction > 70% of the time. In order to achieve these results, multiple infrared (IR) and other sensors were investigated in tandem with the midwave IR (MWIR) OGI to provide additional information to train the underlying models. Information on atmospheric conditions including humidity, temperature, pressure, and solar radiation was provided by a weather station. Several machine learning and deep learning architectures and methods, including looking at quantized classification networks and regressions networks, were explored. As further data was collected, curated, and labeled, it allowed for more refined regressive networks to be adequately trained, leading to better insight into the true flow rates being observed. An important valuable deliverable of this research effort was the development of an advanced network which underwent multiple iterations capable of giving a continuous output. The current network has a predicted mean average percentage error (MAPE) of 12.3% just outside our target goal of 10.00%, but an accuracy of 97.78% at ±50 scfh, well within the overall goal for the Department of Energy (DOE) program. Upon closer inspection, it was observed that more than 10% of datapoints contributing to the MAPE predictions were the result of low flow rate predictions and are beyond the sensitivity of instrument measurement as a result of normal operational variation and noise.

03 NATURAL GAS↗

SU-E-J-273: Simulation of the Radiation Response of Hypoxic Tumors

Purpose: In radiotherapy, it is important to predict the response of tumour to irradiation prior to the treatment. Mathematical modelling of tumour control probability (TCP) based on the dose distribution, medical imaging and other biological information may help to improve this prediction and to optimize the treatment plan. The aim of this work is to develop an image based 3D multiscale radiobiological model, which describes the growth and the response to radiotherapy of hypoxic tumors. Methods: The computer model is based on voxels, containing tumour, normal (including capillary) and dead cells. Killing of tumour cells due to irradiation is calculated by the Linear Quadratic Model (extended for hypoxia), and the proliferation and resorption of cells are modelled by exponential laws. The initial shape of the tumours is taken from CT images and the initial vascular and cell density information from PET and/or MR images. Including the fractionation regime and the physical dose distribution of the radiation treatment, the model simulates the spatial-temporal evolution of the tumor. Additionally, the dose distribution may be biologically optimized. Results: The model describes the appearance of hypoxia during tumour growth and the reoxygenation processes during radiotherapy. Among other parameters, the TCP is calculated for different dose distributions. The results are in accordance with published results. Conclusion: The simulation model may contribute to the understanding of the influence of biological parameters on tumor response during treatment, and specifically on TCP. It may be used to implement dose-painting approaches. Experimental and clinical validation is needed. This study is supported by a grant from the Ministry of Education of Chile, Programa Mece Educacion Superior (2)

60 APPLIED LIFE SCIENCES↗