Search NASA⌕ Search

SEARCH · Search NASA

Results for “C Code”

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,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Computer-controlled finishing via dynamically constraint position-velocity-time scheduler

In a Computer Numerical Controlled (CNC) finishing process, the target material removal from an optical surface is guided by the convolution between the influence function of a machine tool and its dwell time at certain points over the surface. To reduce dynamics stressing and increase machining efficiency, the dwell time must be converted to varying velocities, which are the actual inputs to the machine tool controller. Conventionally, the conversion assumed constant acceleration and relied on linear motion interpolation, which caused discontinuities in velocities. This unsmooth motion affects the material removal distribution, and, thus, the accuracy of the finished surface shape. Many modern CNC machines support the smoother, cubic-polynomial interpolated Position-Velocity-Time (PVT) motion mode; however, the conventional scheduler may fail to provide suitable velocities for the PVT. Here in this study answers this challenge by proposing a novel PVT-based velocity scheduler that achieves smooth motion while considering CNC dynamic limits. Firstly, the principle of the PVT is explained, and the PVT-based velocity scheduler is formulated. Secondly, a quadratic programming is used to optimize the velocities by imposing the CNC dynamic constraints and the C 1 continuities (zeroth and first derivatives are continuous) simultaneously. Thirdly, the smoothness and accuracy of the scheduled velocities are studied on different kinds of tool paths via simulation. Finally, a sub-0.3 nm level surface finishing experiment using ion beam figuring is demonstrated to verify the feasibility of the proposed method. The PVT-based scheduler and simulator code is open-sourced.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Expanding the genetic toolkit: adenine and cytosine base editors for gene disruption in Aspergillus niger

Despite revolutionizing fungal genetic engineering, conventional CRISPR/Cas9-mediated knockouts rely on DNA double-strand breaks (DSBs), which can cause unwanted insertions and deletions, chromosomal abnormalities, and cytotoxicity. Base editors such as adenine base editors (ABEs), which convert A‧T to G‧C, and cytosine base editors (CBEs), which convert C‧G to T‧A, offer a safer alternative by enabling predictable, target-specific single-nucleotide changes without introducing DSBs. To overcome the limitations of traditional genome editing in filamentous fungi, we developed efficient base-editing systems in Aspergillus niger . For the first time, we constructed an ABE in A. niger , achieving up to 80% editing efficiency and inducing predictable A-to-G mutations at the intended intron sites, disrupting gene function through mRNA mis-splicing. We also developed a highly efficient CBE system, capable of introducing premature stop codons with 50–100% efficiency. To broaden the editing scope, we implemented a Cas9-NG variant recognizing a relaxed PAM sequence requiring only a single guanine (G), enabling editing at start codons and splice sites. Leveraging this expanded scope, we established gene disruption approaches by targeting start codons via ABE-mediated A-to-G conversions (ATG-to-GTG and ATG-to-ACG) and CBE-mediated C-to-T conversion (ATG-to-ATA). Additionally, our base-editing systems enable multiplex gRNA delivery and marker-free editing of multiple genes. Collectively, the scope-expanding strategies increase the number of genes targetable for disruption by base-editing in A. niger by 26.3% and enable near-complete coverage of 96% of the coding genes. Overall, this work demonstrates the potential of ABE and CBE systems as versatile, efficient, and safer alternatives to DSBs-based gene disruption in filamentous fungi.

Aspergillus↗
DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

FY17 Status Report on the Initial Development of a Constitutive Model for Grade 91 Steel

Grade 91 is a candidate structural material for high temperature advanced reactor applications. Existing ASME Section III, Subsection HB, Subpart B simplified design rules based on elastic analysis are setup as conservative screening tools with the intent to supplement these screening rules with full inelastic analysis when required. The Code provides general guidelines for suitable inelastic models but does not provide constitutive model implementations. This report describes the development of an inelastic constitutive model for Gr. 91 steel aimed at fulfilling the ASME Code requirements and being included into a new Section III Code appendix, HBB-Z. A large database of over 300 experiments on Gr. 91 was collected and converted to a standard XML form. Five families of Gr. 91 material models were identified in the literature. Of these five, two are potentially suitable for use in the ASME code. These two models were implemented and evaluated against the experimental database. Both models have deficiencies so the report develops a framework for developing and calibrating an improved model. This required creating a new modeling method for representing changes in material rate sensitivity across the full ASME allowable temperature range for Gr. 91 structural components: room temperature to 650° C. On top of this framework for rate sensitivity the report describes calibrating a model for work hardening and softening in the material using genetic algorithm optimization. Future work will focus on improving this trial model by including tension/compression asymmetry observed in experiments and necessary to capture material ratcheting under zero mean stress and by improving the optimization and analysis framework.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase equilibria of advanced technology uranium silicide-based nuclear fuel

The phases in uranium-silicide binary system were evaluated in regards to their stabilities, phase boundaries, crystal structures, and phase transitions. The results from this study were used in combination with a well assessed literature to optimize the U-Si phase diagram using the CALPHAD method. A thermodynamic database was developed, which could be used to guide nuclear fuel fabrication, could be incorporated into other nuclear fuel thermodynamic databases, or could be used to generate data required by fuel performance codes to model fuel behavior in normal or off-normal reactor operations. The U 3 Si 2 and U 3 Si 5 phases were modeled using the Compound Energy Formalism model with 3 sublattices to account for the variation in composition. The crystal structure used for the USi phase was the tetragonal with an I4/mmm space. Above 450°C, the U 3 Si 5 phase was modeled. The composition of the USi 2 phase was adjusted to USi 1.84 . The calculated invariant reactions and the enthalpy of formation for the stoichiometric phases were in agreement with experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Vadose Zone Model for B-3A/B Pond for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort; carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP) simulator officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the B-3A/B Ponds model. The scope of this ECF is to document the development and results of the B-3 A/B Ponds vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Vadose Zone Model for SALDS for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP) is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the State-Approved Land Disposal Site (SALDS) model. The scope of this ECF is to document the development and results of the SALDS vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Dose-Dependent Cortical Thinning After Partial Brain Irradiation in High-Grade Glioma

Purpose: Radiation-induced cognitive deficits may be mediated by tissue damage to cortical regions. Volumetric changes in cortex can be reliably measured using high-resolution magnetic resonance imaging (MRI). We used these methods to study the association between radiation therapy (RT) dose and change in cortical thickness in high-grade glioma (HGG) patients. Methods and Materials: We performed a voxel-wise analysis of MRI from 15 HGG patients who underwent fractionated partial brain RT. Three-dimensional MRI was acquired pre- and 1 year post RT. Cortex was parceled with well-validated segmentation software. Surgical cavities were censored. Each cortical voxel was assigned a change in cortical thickness between time points, RT dose value, and neuroanatomic label by lobe. Effects of dose, neuroanatomic location, age, and chemotherapy on cortical thickness were tested using linear mixed effects (LME) modeling. Results: Cortical atrophy was seen after 1 year post RT with greater effects at higher doses. Estimates from LME modeling showed that cortical thickness decreased by −0.0033 mm (P<.001) for every 1-Gy increase in RT dose. Temporal and limbic cortex exhibited the largest changes in cortical thickness per Gy compared to that in other regions (P<.001). Age and chemotherapy were not significantly associated with change in cortical thickness. Conclusions: We found dose-dependent thinning of the cerebral cortex, with varying neuroanatomical regional sensitivity, 1 year after fractionated partial brain RT. The magnitude of thinning parallels 1-year atrophy rates seen in neurodegenerative diseases and may contribute to cognitive decline following high-dose RT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Pushered single shell implosions for mix and radiation trapping studies using high-Z layers on National Ignition Facility

Pushered Single Shells (PSSs) are another approach to Inertial Confinement Fusion implosions that employ high-Z materials in the innermost capsule layer (pusher) as a means to enhance radiation trapping and lower core ignition requirements. However, adding high-Z materials can also increase losses due to mix, provide extra tamping, and make the capsule emission opaque to x-ray diagnostics. The first PSS implosions performed on the National Ignition Facility use plastic ablators with a germanium (Ge) dopant as a high-Z surrogate in the pusher to isolate the effects of high-Z mix and radiation trapping without changing tamping. Using a 2-shock laser pulse, the PSS implosions are designed and symmetrized to reach 3.7 keV core temperatures. A low concentration (2.8%) Ge dopant is added to the innermost layer, and the resulting effects on mix and x-ray opacity are observed. The method of separated reactants is used to infer information about mixing between the deuterated plastic pusher and the capsule fill gas (25% tritium) from the resulting nuclear DT reactions. Radiation transport is studied via capsule emission x-ray spectroscopy and imaging. Both nuclear and x-ray data corroborate the hypothesis that the addition of Ge strongly affects the mix region through radiation losses but has a minimal effect on the core and the warm, unmixed regions. Simulations using diffusive and turbulent mix models agree qualitatively with data, but quantitative agreement may require hybrid mix models that can model the transitional regime between turbulence and diffusion. Simulations matching the observables show increased core radiation trapping when Ge is added.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2015

SIMULATION OF ASTRONOMICAL IMAGES FROM OPTICAL SURVEY TELESCOPES USING A COMPREHENSIVE PHOTON MONTE CARLO APPROACH

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons s{sup −1}, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to (1) construct kilometer scale phase screens necessary for wide-field telescopes, (2) reproduce atmospheric point-spread function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, (3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and (4) recover system effective area predicted from analytic photometry integrals. This new code, the Photon Simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

TILDA-V: A full-differential code for proton tracking in biological matter

Understanding the radiation-induced effects at the cellular level is of prime importance for predicting the fate of irradiated biological organisms. Thus, whether it is in radiobiology to identify the DNA critical lesions or in medicine to adapt the radio-therapeutic protocols, an accurate knowledge of the numerous interactions induced by charged particles in living matter is required. Monte-Carlo track-structure simulations represent the most suitable and powerful tools, in particular for modelling the full slowing-down of the ionizing particles in biological matter. Furthermore more of the existing codes are based on semi-empirical cross sections as well as the use of water as surrogate of the biological matter.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2023

Snekmer: a scalable pipeline for protein sequence fingerprinting based on amino acid recoding

Abstract Motivation The vast expansion of sequence data generated from single organisms and microbiomes has precipitated the need for faster and more sensitive methods to assess evolutionary and functional relationships between proteins. Representing proteins as sets of short peptide sequences (kmers) has been used for rapid, accurate classification of proteins into functional categories; however, this approach employs an exact-match methodology and thus may be limited in terms of sensitivity and coverage. We have previously used similarity groupings, based on the chemical properties of amino acids, to form reduced character sets and recode proteins. This amino acid recoding (AAR) approach simplifies the construction of protein representations in the form of kmer vectors, which can link sequences with distant sequence similarity and provide accurate classification of problematic protein families. Results Here, we describe Snekmer, a software tool for recoding proteins into AAR kmer vectors and performing either (i) construction of supervised classification models trained on input protein families or (ii) clustering for de novo determination of protein families. We provide examples of the operation of the tool against a set of nitrogen cycling families originally collected using both standard hidden Markov models and a larger set of proteins from Uniprot and demonstrate that our method accurately differentiates these sequences in both operation modes. Availability and implementation Snekmer is written in Python using Snakemake. Code and data used in this article, along with tutorial notebooks, are available at http://github.com/PNNL-CompBio/Snekmer under an open-source BSD-3 license. Supplementary information Supplementary data are available at Bioinformatics Advances online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Unsupervised Segmentation and Clustering Workflow for Efficient Processing of 4D-STEM and 5D-STEM Data

Four-dimensional scanning transmission electron microscopy (4D-STEM) enables mapping of diffraction information with nanometer-scale spatial resolution, offering detailed insight into local structure, orientation, and strain. However, as data dimensionality and sampling density increase, particularly for in situ scanning diffraction experiments (5D-STEM), robust segmentation of structurally consistent behavior across sequential measurements becomes essential for efficient and physically meaningful analysis. Here, we introduce a clustering framework that identifies crystallographically distinct domains from 4D-STEM datasets. By using local diffraction-pattern similarity as a metric, the method extracts closed contours delineating spatially contiguous regions. This approach produces cluster-averaged diffraction patterns that improve signal quality while reducing data volume by orders of magnitude, enabling rapid and accurate orientation, phase, and strain mapping. We demonstrate the applicability of this approach to in situ liquid-cell 4D-STEM data of gold nanoparticle growth. Our method provides a scalable and generalizable route for spatially coherent segmentation, data compression, and quantitative structure–strain mapping across diverse 4D-STEM modalities. The full analysis code and example workflows are publicly available to support reproducibility and reuse.

4D-STEM↗
DOE OSTI2024

Matter power spectra in modified gravity: a comparative study of approximations and N -body simulations

ABSTRACT Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade’s largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon, and K-mouflage. In particular, we compare the predictions from N-body simulations solving the full scalar field equation, two N-body codes with approximate time integration schemes, a parametrized modified N-body implementation, and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in 2 per cent agreement at $z\le 1$ and $k\le 1~h\,{\rm Mpc}^{-1}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.

Bose, B. (ORCID:0000000319658614)↗
DOE OSTI2015

SU-E-T-761: TOMOMC, A Monte Carlo-Based Planning VerificationTool for Helical Tomotherapy

Purpose: Present a new Monte Carlo code (TOMOMC) to calculate 3D dose distributions for patients undergoing helical tomotherapy treatments. TOMOMC performs CT-based dose calculations using the actual dynamic variables of the machine (couch motion, gantry rotation, and MLC sequences). Methods: TOMOMC is based on the GEPTS (Gama Electron and Positron Transport System) general-purpose Monte Carlo system (Chibani and Li, Med. Phys. 29, 2002, 835). First, beam models for the Hi-Art Tomotherpy machine were developed for the different beam widths (1, 2.5 and 5 cm). The beam model accounts for the exact geometry and composition of the different components of the linac head (target, primary collimator, jaws and MLCs). The beams models were benchmarked by comparing calculated Pdds and lateral/transversal dose profiles with ionization chamber measurements in water. See figures 1–3. The MLC model was tuned in such a way that tongue and groove effect, inter-leaf and intra-leaf transmission are modeled correctly. See figure 4. Results: By simulating the exact patient anatomy and the actual treatment delivery conditions (couch motion, gantry rotation and MLC sinogram), TOMOMC is able to calculate the 3D patient dose distribution which is in principal more accurate than the one from the treatment planning system (TPS) since it relies on the Monte Carlo method (gold standard). Dose volume parameters based on the Monte Carlo dose distribution can also be compared to those produced by the TPS. Attached figures show isodose lines for a H&N patient calculated by TOMOMC (transverse and sagittal views). Analysis of differences between TOMOMC and TPS is an ongoing work for different anatomic sites. Conclusion: A new Monte Carlo code (TOMOMC) was developed for Tomotherapy patient-specific QA. The next step in this project is implementing GPU computing to speed up Monte Carlo simulation and make Monte Carlo-based treatment verification a practical solution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-72: A Clinical Usable Integrated Contouring Quality Evaluation Software for Radiotherapy

Purpose: To introduce the Auto Contour Evaluation (ACE) software, which is the clinical usable, user friendly, efficient and all-in-one toolbox for automatically identify common contouring errors in radiotherapy treatment planning using supervised machine learning techniques. Methods: ACE is developed with C# using Microsoft .Net framework and Windows Presentation Foundation (WPF) for elegant GUI design and smooth GUI transition animations through the integration of graphics engines and high dots per inch (DPI) settings on modern high resolution monitors. The industrial standard software design pattern, Model-View-ViewModel (MVVM) pattern, is chosen to be the major architecture of ACE for neat coding structure, deep modularization, easy maintainability and seamless communication with other clinical software. ACE consists of 1) a patient data importing module integrated with clinical patient database server, 2) a 2D DICOM image and RT structure simultaneously displaying module, 3) a 3D RT structure visualization module using Visualization Toolkit or VTK library and 4) a contour evaluation module using supervised pattern recognition algorithms to detect contouring errors and display detection results. ACE relies on supervised learning algorithms to handle all image processing and data processing jobs. Implementations of related algorithms are powered by Accord.Net scientific computing library for better efficiency and effectiveness. Results: ACE can take patient’s CT images and RT structures from commercial treatment planning software via direct user input or from patients’ database. All functionalities including 2D and 3D image visualization and RT contours error detection have been demonstrated with real clinical patient cases. Conclusion: ACE implements supervised learning algorithms and combines image processing and graphical visualization modules for RT contours verification. ACE has great potential for automated radiotherapy contouring quality verification. Structured with MVVM pattern, it is highly maintainable and extensible, and support smooth connections with other clinical software tools.

60 APPLIED LIFE SCIENCES↗

Refine your search to explore more results.