Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surrogates”

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,369 records · Page 76

Combining synchrotron X-ray diffraction, mechanistic modeling and machine learning for in situ subsurface temperature quantification during laser melting

Laser melting, such as that encountered during additive manufacturing, produces extreme gradients of temperature in both space and time, which in turn influence microstructural development in the material. Qualification and model validation of the process itself and the resulting material necessitate the ability to characterize these temperature fields. However, well established means to directly probe the material temperature below the surface of an alloy while it is being processed are limited. To address this gap in characterization capabilities, a novel means is presented to extract subsurface temperature-distribution metrics, with uncertainty, from in situ synchrotron X-ray diffraction measurements to provide quantitative temperature evolution data during laser melting. Temperature-distribution metrics are determined using Gaussian process regression supervised machine-learning surrogate models trained with a combination of mechanistic modeling (heat transfer and fluid flow) and X-ray diffraction simulation. The trained surrogate model uncertainties are found to range from 5 to 15% depending on the metric and current temperature. The surrogate models are then applied to experimental data to extract temperature metrics from an Inconel 625 nickel superalloy wall specimen during laser melting. The maximum temperatures of the solid phase in the diffraction volume through melting and cooling are found to reach the solidus temperature as expected, with the mean and minimum temperatures found to be several hundred degrees less. The extracted temperature metrics near melting are determined to be more accurate because of the lower relative levels of mechanical elastic strains. However, uncertainties for temperature metrics during cooling are increased due to the effects of thermomechanical stress.

36 MATERIALS SCIENCE↗

A General Framework for Bounding Approximate Dynamic Programming Schemes

For years, there has been interest in approximation methods for solving dynamic programming problems, because of the inherent complexity in computing optimal solutions characterized by Bellman’s principle of optimality. A wide range of approximate dynamic programming (ADP) methods now exists. It is of great interest to guarantee that the performance of an ADP scheme be at least some known fraction, say ß , of optimal. This letter introduces a general approach to bounding the performance of ADP methods, in this sense, in the stochastic setting. The approach is based on new results for bounding greedy solutions in string optimization problems, where one has to choose a string (ordered set) of actions to maximize an objective function. This bounding technique is inspired by submodularity theory, but submodularity is not required for establishing bounds. Instead, the bounding is based on quantifying certain notions of curvature of string functions; the smaller the curvatures the better the bound. The key insight is that any ADP scheme is a greedy scheme for some surrogate string objective function that coincides in its optimal solution and value with those of the original optimal control problem. The ADP scheme then yields to the bounding technique mentioned above, and the curvatures of the surrogate objective determine the value ß of the bound. The surrogate objective and its curvatures depend on the specific ADP.

discrete event systems↗

A Data-Driven Global Sensitivity Analysis Framework for Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.

14 SOLAR ENERGY↗

AN AUTOMATED MACHINE LEARNING-GENETIC ALGORITHM FRAMEWORK WITH ACTIVE LEARNING FOR DESIGN OPTIMIZATION

The use of machine learning (ML)-based surrogate models is a promising technique to significantly accelerate simulation-driven design optimization of internal combustion (IC) engines, due to the high computational cost of running computational fluid dynamics (CFD) simulations. However, training the ML models requires hyperparameter selection, which is often done using trial-and-error and domain expertise. Another challenge is that the data required to train these models are often unknown a priori. In this work, we present an automated hyperparameter selection technique coupled with an active learning approach to address these challenges. The technique presented in this study involves the use of a Bayesian approach to optimize the hyperparameters of the base learners that make up a super learner model. In addition to performing hyperparameter optimization (HPO), an active learning approach is employed, where the process of data generation using simulations, ML training, and surrogate optimization is performed repeatedly to refine the solution in the vicinity of the predicted optimum. The proposed approach is applied to the optimization of a compression ignition engine with control parameters relating to fuel injection, in-cylinder flow, and thermodynamic conditions. It is demonstrated that by automatically selecting the best values of the hyperparameters, a 1.6% improvement in merit value is obtained, compared to an improvement of 1.0% with default hyperparameters. Overall, the framework introduced in this study reduces the need for technical expertise in training ML models for optimization while also reducing the number of simulations needed for performing surrogate-based design optimization.

Owoyele, Opeoluwa↗

Evaluating performance of high efficiency mist eliminators

Processing liquid wastes frequently generates off gas streams with high humidity and liquid aerosols. Droplet laden air streams can be produced from tank mixing or sparging and processes such as reforming or evaporative volume reduction. Unfortunately these wet air streams represent a genuine threat to HEPA filters. High efficiency mist eliminators (HEME) are one option for removal of liquid aerosols with high dissolved or suspended solids content. HEMEs have been used extensively in industrial applications, however they have not seen widespread use in the nuclear industry. Filtering efficiency data along with loading curves are not readily available for these units and data that exist are not easily translated to operational parameters in liquid waste treatment plants. A specialized test stand has been developed to evaluate the performance of HEME elements under use conditions of a US DOE facility. HEME elements were tested at three volumetric flow rates using aerosols produced from an iron-rich waste surrogate. The challenge aerosol included submicron particles produced from Laskin nozzles and super micron particles produced from a hollow cone spray nozzle. Test conditions included ambient temperature and relative humidities greater than 95%. Data collected during testing HEME elements from three different manufacturers included volumetric flow rate, differential temperature across the filter housing, downstream relative humidity, and differential pressure (dP) across the filter element. Filter challenge was discontinued at three intermediate dPs and the filter to allow determining filter efficiency using dioctyl phthalate and then with dry surrogate aerosols. Filtering efficiencies of the clean HEME, the clean HEME loaded with water, and the HEME at maximum dP were also collected using the two test aerosols. Results of the testing included differential pressure vs. time loading curves for the nine elements tested along with the mass of moisture and solid material on each element at final dP. Plots of overall filtering efficiencies for DOP (spherical aerosol) and dry surrogate (aspherical aerosols) at specified dPs were computed for each filter. Filtering efficiencies were determined as a function of particle size. Curves were also generated showing the most penetrating particle size as a function of dP. A preliminary set of tests was conducted to evaluate spray location, duration, pressure, and wash volume for in-place cleaning the interior surface (reducing dP) of the HEME element. A variety of nozzle designs were evaluated and test results demonstrated the potential to overload the HEME (saturate filter medium) resulting in very high dPs and extensive drain times. At least one combination of spray nozzle design, spray location on the surface of the element, and spray time/pressure was successful in achieving extension of operational life. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

TU-F-17A-03: A 4D Lung Phantom for Coupled Registration/Segmentation Evaluation

Purpose: Coupling the processes of segmentation and registration (regmentation) is a recent development that allows improved efficiency and accuracy for both steps and may improve the clinical feasibility of online adaptive radiotherapy. Presented is a multimodality animal tissue model designed specifically to provide a ground truth to simultaneously evaluate segmentation and registration errors during respiratory motion. Methods: Tumor surrogates were constructed from vacuum sealed hydrated natural sea sponges with catheters used for the injection of PET radiotracer. These contained two compartments allowing for two concentrations of radiotracer mimicking both tumor and background signals. The lungs were inflated to different volumes using an air pump and flow valve and scanned using PET/CT and MRI. Anatomical landmarks were used to evaluate the registration accuracy using an automated bifurcation tracking pipeline for reproducibility. The bifurcation tracking accuracy was assessed using virtual deformations of 2.6 cm, 5.2 cm and 7.8 cm of a CT scan of a corresponding human thorax. Bifurcations were detected in the deformed dataset and compared to known deformation coordinates for 76 points. Results: The bifurcation tracking accuracy was found to have a mean error of −0.94, 0.79 and −0.57 voxels in the left-right, anterior-posterior and inferior-superior axes using a 1×1×5 mm3 resolution after the CT volume was deformed 7.8 cm. The tumor surrogates provided a segmentation ground truth after being registered to the phantom image. Conclusion: A swine lung model in conjunction with vacuum sealed sponges and a bifurcation tracking algorithm is presented that is MRI, PET and CT compatible and anatomically and kinetically realistic. Corresponding software for tracking anatomical landmarks within the phantom shows sub-voxel accuracy. Vacuum sealed sponges provide realistic tumor surrogate with a known boundary. A ground truth with minimal uncertainty is thus realized that can be used for comparing the performance of registration and segmentation algorithms.

60 APPLIED LIFE SCIENCES↗

SU-E-J-155: Automatic Quantitative Decision Making Metric for 4DCT Image Quality

Purpose: To develop a quantitative decision making metric for automatically detecting irregular breathing using a large patient population that received phase-sorted 4DCT. Methods: This study employed two patient cohorts. Cohort#1 contained 256 patients who received a phasesorted 4DCT. Cohort#2 contained 86 patients who received three weekly phase-sorted 4DCT scans. A previously published technique used a single abdominal surrogate to calculate the ratio of extreme inhalation tidal volume to normal inhalation tidal volume, referred to as the κ metric. Since a single surrogate is standard for phase-sorted 4DCT in radiation oncology clinical practice, tidal volume was not quantified. Without tidal volume, the absolute κ metric could not be determined, so a relative κ (κrel) metric was defined based on the measured surrogate amplitude instead of tidal volume. Receiver operator characteristic (ROC) curves were used to quantitatively determine the optimal cutoff value (jk) and efficiency cutoff value (τk) of κrel to automatically identify irregular breathing that would reduce the image quality of phase-sorted 4DCT. Discriminatory accuracy (area under the ROC curve) of κrel was calculated by a trapezoidal numeric integration technique. Results: The discriminatory accuracy of ?rel was found to be 0.746. The key values of jk and tk were calculated to be 1.45 and 1.72 respectively. For values of ?rel such that jk≤κrel≤τk, the decision to reacquire the 4DCT would be at the discretion of the physician. This accounted for only 11.9% of the patients in this study. The magnitude of κrel held consistent over 3 weeks for 73% of the patients in cohort#3. Conclusion: The decision making metric, ?rel, was shown to be an accurate classifier of irregular breathing patients in a large patient population. This work provided an automatic quantitative decision making metric to quickly and accurately assess the extent to which irregular breathing is occurring during phase-sorted 4DCT.

60 APPLIED LIFE SCIENCES↗

SU-E-J-163: A Biomechanical Lung Model for Respiratory Motion Study

Purpose: This work presents a biomechanical model to investigate the complex respiratory motion for the lung tumor tracking in radiosurgery by computer simulation. Methods: The models include networked massspring-dampers to describe the tumor motion, different types of surrogate signals, and the force generated by the diaphragm. Each mass-springdamper has the same mechanical structure and each model can have different numbers of mass-spring-dampers. Both linear and nonlinear stiffness parameters were considered, and the damping ratio was tuned in a range so that the tumor motion was over-damped (no natural tumor oscillation occurs without force from the diaphragm). The simulation was run by using ODE45 (ordinary differential equations by Runge-Kutta method) in MATLAB, and all time courses of motions and inputs (force) were generated and compared. Results: The curvature of the motion time courses around their peaks was sensitive to the damping ratio. Therefore, the damping ratio can be determined based on the clinical data of a high sampling rate. The peak values of different signals and the time the peaks occurred were compared, and it was found that the diaphragm force had a time lead over the tumor motion, and the lead time (0.1–0.4 seconds) depended on the distance between the tumor and the diaphragm. Conclusion: We reported a model based analysis approach for the spatial and temporal relation between the motion of the lung tumor and the surrogate signals. Due to the phase lead of the diaphragm in comparing with the lung tumor motion, the measurement of diaphragm motion (or its electromyography signal) can be used as a beam gating signal in radiosurgery, and it can also be an additional surrogate signal for better tumor motion tracking. The research is funded by the American Cancer Society (ACS) grant. The grant name is: Frameless SRS Based on Robotic Head Motion Cancellation. The grant number is: RSG-13-313-01-CCE.

60 APPLIED LIFE SCIENCES↗

SU-E-J-186: Using 4DCT-Based Motion Modeling to Predict Motion and Duty Cycle On Successive Days of Gated Radiotherapy

Purpose: To determine if 4DCT-based motion modeling and external surrogate motion measured during treatment simulation can enhance prediction of residual tumor motion and duty cycle during treatment delivery. Methods: This experiment was conducted using simultaneously recorded tumor and external surrogate motion acquired over multiple fractions of lung cancer radiotherapy. These breathing traces were combined with the XCAT phantom to simulate CT images. Data from the first day was used to estimate the residual tumor motion and duty cycle both directly from the 4DCT (the current clinical standard), and from external-surrogate based motion modeling. The accuracy of these estimated residual tumor motions and duty cycles are evaluated by comparing to the measured internal/external motions from other treatment days. Results: All calculations were done for 25% and 50% duty cycles. The results indicated that duty cycle derived from 4DCT information alone is not enough to accurately predict duty cycles during treatment. Residual tumor motion was determined from the recorded data and compared with the estimated residual tumor motion from 4DCT. Relative differences in residual tumor motion varied from −30% to 55%, suggesting that more information is required to properly predict residual tumor motion. Compared to estimations made from 4DCT, in three out of four patients examined, the 30 seconds of motion modeling data was able to predict the duty cycle with better accuracy than 4DCT. No improvement was observed in prediction of residual tumor motion for this dataset. Conclusion: Motion modeling during simulation has the potential to enhance 4DCT and provide more information about target motion, duty cycles, and delivered dose. Based on these four patients, 30 seconds of motion modeling data produced improve duty cycle estimations but showed no measurable improvement in residual tumor motion prediction. More patient data is needed to verify this Result. I would like to acknowledge funding from MRA, VARIAN Medical Systems, Inc.

60 APPLIED LIFE SCIENCES↗

Liver 4DMRI: A retrospective image-based sorting method

Purpose: Four-dimensional magnetic resonance imaging (4DMRI) is an emerging technique in radiotherapy treatment planning for organ motion quantification. In this paper, the authors present a novel 4DMRI retrospective image-based sorting method, providing reduced motion artifacts than using a standard monodimensional external respiratory surrogate. Methods: Serial interleaved 2D multislice MRI data were acquired from 24 liver cases (6 volunteers + 18 patients) to test the proposed 4DMRI sorting. Image similarity based on mutual information was applied to automatically identify a stable reference phase and sort the image sequence retrospectively, without the use of additional image or surrogate data to describe breathing motion. Results: The image-based 4DMRI provided a smoother liver profile than that obtained from standard resorting based on an external surrogate. Reduced motion artifacts were observed in image-based 4DMRI datasets with a fitting error of the liver profile measuring 1.2 ± 0.9 mm (median ± interquartile range) vs 2.1 ± 1.7 mm of the standard method. Conclusions: The authors present a novel methodology to derive a patient-specific 4DMRI model to describe organ motion due to breathing, with improved image quality in 4D reconstruction.

60 APPLIED LIFE SCIENCES↗

Technical Report: TG-142 compliant and comprehensive quality assurance tests for respiratory gating

Purpose: To develop and establish a comprehensive gating commissioning and quality assurance procedure in compliance with TG-142. Methods: Eight Varian TrueBeam Linacs were used for this study. Gating commissioning included an end-to-end test and baseline establishment. The end-to-end test was performed using a CIRS dynamic thoracic phantom with a moving cylinder inside the lung, which was used for carrying both optically simulated luminescence detectors (OSLDs) and Gafchromic EBT2 films while the target is moving, for a point dose check and 2D profile check. In addition, baselines were established for beam-on temporal delay and calibration of the surrogate, for both megavoltage (MV) and kilovoltage (kV) beams. A motion simulation device (MotionSim) was used to provide periodic motion on a platform, in synchronizing with a surrogate motion. The overall accuracy and uncertainties were analyzed and compared. Results: The OSLD readings were within 5% compared to the planned dose (within measurement uncertainty) for both phase and amplitude gated deliveries. Film results showed less than 3% agreement to the predicted dose with a standard sinusoid motion. The gate-on temporal accuracy was averaged at 139 ± 10 ms for MV beams and 92 ± 11 ms for kV beams. The temporal delay of the surrogate motion depends on the motion speed and was averaged at 54.6 ± 3.1 ms for slow, 24.9 ± 2.9 ms for intermediate, and 23.0 ± 20.1 ms for fast speed. Conclusions: A comprehensive gating commissioning procedure was introduced for verifying the output accuracy and establishing the temporal accuracy baselines with respiratory gating. The baselines are needed for routine quality assurance tests, as suggested by TG-142.

60 APPLIED LIFE SCIENCES↗

SU-F-T-91: Development of Real Time Abdominal Compression Force (ACF) Monitoring System

Purpose: Hard-plate based abdominal compression is known to be effective, but no explicit method exists to quantify abdominal compression force (ACF) and maintain the proper ACF through the whole procedure. In addition, even with compression, it is necessary to do 4D CT to manage residual motion but, 4D CT is often not possible due to reduced surrogating sensitivity. In this study, we developed and evaluated a system that both monitors ACF in real time and provides surrogating signal even under compression. The system can also provide visual-biofeedback. Methods: The system developed consists of a compression plate, an ACF monitoring unit and a visual-biofeedback device. The ACF monitoring unit contains a thin air balloon in the size of compression plate and a gas pressure sensor. The unit is attached to the bottom of the plate thus, placed between the plate and the patient when compression is applied, and detects compression pressure. For reliability test, 3 volunteers were directed to take several different breathing patterns and the ACF variation was compared with the respiratory flow and external respiratory signal to assure that the system provides corresponding behavior. In addition, guiding waveform were generated based on free breathing, and then applied for evaluating the effectiveness of visual-biofeedback. Results: We could monitor ACF variation in real time and confirmed that the data was correlated with both respiratory flow data and external respiratory signal. Even under abdominal compression, in addition, it was possible to make the subjects successfully follow the guide patterns using the visual biofeedback system. Conclusion: The developed real time ACF monitoring system was found to be functional as intended and consistent. With the capability of both providing real time surrogating signal under compression and enabling visual-biofeedback, it is considered that the system would improve the quality of respiratory motion management in radiation therapy. This research was supported by the Mid-career Researcher Program through NRF funded by the Ministry of Science, ICT & Future Planning of Korea (NRF-2014R1A2A1A10050270) and by the Radiation Technology R&D program through the National Research Foundation of Korea funded by the Ministry of Science, ICT & Future Planning (No. 2013M2A2A7038291)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-F-T-113: Inherent Functional Dependence of Spinal Cord Doses of Variable Irradiated Volumes in Spine SBRT

Purpose: Spinal cord tolerance for SBRT has been recommended for the maximum point dose level or at irradiated volumes such as 0.35 mL or 10% of contoured volumes. In this study, we investigated an inherent functional relationship that associates these dose surrogates for irradiated spinal cord volumes of up to 3.0 mL. Methods: A hidden variable termed as Effective Dose Radius (EDR) was formulated based on a dose fall-off model to correlate dose at irradiated spinal cord volumes ranging from 0 mL (point maximum) to 3.0 mL. A cohort of 15 spine SBRT cases was randomly selected to derive an EDR-parameterized formula. The mean prescription dose for the studied cases was 21.0±8.0 Gy (range, 10–40Gy) delivered in 3±1 fractions with target volumes of 39.1 ± 70.6 mL. Linear regression and variance analysis were performed for the fitting parameters of variable EDR values. Results: No direct correlation was found between the dose at maximum point and doses at variable spinal cord volumes. For example, Pearson R{sup 2} = 0.643 and R{sup 2}= 0.491 were obtained when correlating the point maximum dose with the spinal cord dose at 1 mL and 3 mL, respectively. However, near perfect correlation (R{sup 2} ≥0.99) was obtained when corresponding parameterized EDRs. Specifically, Pearson R{sup 2}= 0.996 and R{sup 2} = 0.990 were obtained when correlating EDR (maximum point dose) with EDR (dose at 1 mL) and EDR(dose at 3 mL), respectively. As a result, high confidence level look-up tables were established to correlate spinal cord doses at the maximum point to any finite irradiated volumes. Conclusion: An inherent functional relationship was demonstrated for spine SBRT. Such a relationship unifies dose surrogates at variable cord volumes and proves that a single dose surrogate (e.g. point maximum dose) is mathematically sufficient in constraining the overall spinal cord dose tolerance for SBRT.

60 APPLIED LIFE SCIENCES↗

SU-G-JeP1-08: Dual Modality Verification for Respiratory Gating Using New Real- Time Tumor Tracking Radiotherapy System

Purpose: The respirato ry gating system combined the TrueBeam and a new real-time tumor-tracking radiotherapy system (RTRT) was installed. The RTRT system consists of two x-ray tubes and color image intensifiers. Using fluoroscopic images, the fiducial marker which was implanted near the tumor was tracked and was used as the internal surrogate for respiratory gating. The purposes of this study was to develop the verification technique of the respiratory gating with the new RTRT using cine electronic portal image device images (EPIDs) of TrueBeam and log files of the RTRT. Methods: A patient who underwent respiratory gated SBRT of the lung using the RTRT were enrolled in this study. For a patient, the log files of three-dimensional coordinate of fiducial marker used as an internal surrogate were acquired using the RTRT. Simultaneously, the cine EPIDs were acquired during respiratory gated radiotherapy. The data acquisition was performed for one field at five sessions during the course of SBRT. The residual motion errors were calculated using the log files (E{sub log}). The fiducial marker used as an internal surrogate into the cine EPIDs was automatically extracted by in-house software based on the template-matching algorithm. The differences between the the marker positions of cine EPIDs and digitally reconstructed radiograph were calculated (E{sub EPID}). Results: Marker detection on EPID using in-house software was influenced by low image contrast. For one field during the course of SBRT, the respiratory gating using the RTRT showed the mean ± S.D. of 95{sup th} percentile E{sub EPID} were 1.3 ± 0.3 mm,1.1 ± 0.5 mm,and those of E{sub log} were 1.5 ± 0.2 mm, 1.1 ± 0.2 mm in LR and SI directions, respectively. Conclusion: We have developed the verification method of respiratory gating combined TrueBeam and new real-time tumor-tracking radiotherapy system using EPIDs and log files.

60 APPLIED LIFE SCIENCES↗

MO-FG-BRA-06: Electromagnetic Beacon Insertion in Lung Cancer Patients and Resultant Surrogacy Errors for Dynamic MLC Tumour Tracking

Purpose: To assess endo-bronchial electromagnetic beacon insertion and to quantify the geometric accuracy of using beacons as a surrogate for tumour motion in real-time multileaf collimator (MLC) tracking of lung tumours. Methods: The LIGHT SABR trial is a world-first clinical trial in which the MLC leaves move with lung tumours in real time on a standard linear accelerator. Tracking is performed based on implanted electromagnetic beacons (CalypsoTM, Varian Medical Systems, USA) as a surrogate for tumour motion. Five patients have been treated and have each had three beacons implanted endo-bronchially under fluoroscopic guidance. The centre of mass (C.O.M) has been used to adapt the MLC in real-time. The geometric error in using the beacon C.O.M as a surrogate for tumour motion was measured by measuring the tumour and beacon C.O.M in all phases of the respiratory cycle of a 4DCT. The surrogacy error was defined as the difference in beacon and tumour C.O.M relative to the reference phase (maximum exhale). Results: All five patients have had three beacons successfully implanted with no migration between simulation and end of treatment. Beacon placement relative to tumour C.O.M varied from 14 to 74 mm and in one patient spanned two lobes. Surrogacy error was measured in each patient on the simulation 4DCT and ranged from 0 to 3 mm. Surrogacy error as measured on 4DCT was subject to artefacts in mid-ventilation phases. Surrogacy error was a function of breathing phase and was typically larger at maximum inhale. Conclusion: Beacon placement and thus surrogacy error is a major component of geometric uncertainty in MLC tracking of lung tumours. Surrogacy error must be measured on each patient and incorporated into margin calculation. Reduction of surrogacy error is limited by airway anatomy, however should be taken into consideration when performing beacon insertion and planning. This research is funded by Varian Medical Systems via a collaborative research agreement.

60 APPLIED LIFE SCIENCES↗

Matlab Code for Multiscale Topology Optimization

This Matlab script implements a simple topology optimization algorithm that minimizes the compliance of a structure within a 3D prismatic design domain. It uses a simple linear elastic finite element solver to solve the equilibrium equations, and the optimality criterion (00) method to update the design iteratively. This functionality is identical to, or the simplest possible extension from 2D to 3D of, a 99-line Matlab code released in 2001 in a journal article by Sigmund. What is new in this code is our implementation of simple surrogate models for several truss micro-architectures, which converts this code to a multiscale design tool, whereas the original code is single-scale. The development of these surrogate models and their implementation is described in a journal article we intend to submit. We plan to include the full source code in an Appendix in this paper. The code’s primary function is as a pedagogical tool rather than as a finished software product. Sigmund’s 2001 paper describes the implementation line-by-line. We are including this updated version of the code to allow the reader to understand exactly what changes must be made to a “standard” topology optimization code to implement the surrogate models and thus multiscale design optimization.

Watts, Seth↗

Neural Phase Simulation

The Neural Phase Simulation (NPS) is a package of codes for simulating microstructure evolution and accelerated molecular dynamics with deep neural-networks based surrogate models. NPS is designed to offer quantitatively accurate and computationally efficient simulation capabilities by leveraging modern machine-learning techniques. The primary intended use cases of NPS are training neural network surrogate models, though performing simulations on a single node is also supported. The NPS surrogate models can be trained from ground truth simulation methods, which are supposed to be accurate but expensive, such as molecular dynamics, phase field methods, kinetic Monte Carlo and discrete dislocation dynamics.

Zhou, Fei↗

Fusion RF Modeling Machine Learning (FusionML_RF) v1.0

FusionML_RF consists of multiple codes and trained machine learning (ML) models that perform low-cost output modeling from the Genray-CQL3D. Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. For example, completing a single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time. On the other hand, these ML models achieve ~ms of inference time with high accuracy across the input parameter space. This software collection consists of multiple components. (1) codes that use ML methods and precomputed Genray-CQL3D simulation output to build regression models that enable approximate computations of Genray-CLQ3D outputs from arbitrary but physically meaningful input parameters (surrogate modeling); (2) three trained models created by the team, using a database of 16,000+ GENRAY/CQL3D simulations, to study the performance of ML models for surrogate modeling; (3) codes that load the trained models and simulation data, and then compute mean squared error between the models' predictions and the ground truth of simulation output data. This collection is being made available in conjunction with a scientific publication about the work to promote reusability and provide an artifact of the scientific work.

Bai, Zhe↗