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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-J-225: Histology Study of MR Guided Pulsed Focused Ultrasound On Treatment of Prostate Cancer in Vivo

Purpose: Our previous study demonstrated significant tumor growth delay in the mice treated with pulsed high intensity focused ultrasound (pHIFU). The purpose of this study is to understand the cell killing mechanisms of pHIFU. Methods: Prostate cancer cells (LNCaP), were grown orthotopically in 17 nude mice. Tumor-bearing mice were treated using pHIFU with an acoustic power of 25W, pulse width 100msec and 300 pulses in one sonication under MR guidance. Mutiple sonications were used to cover the whole tumor volume. The temperature (less than 40 degree centigrade in the focal spot) was monitored using MR thermometry. Animals were euthanized at pre-determined time points (n=2) after treatment: 0 hours; 6 hrs; 24 hrs; 48 hrs; 4 days and 7 days. Two tumorbearing mice were used as control. Three tumor-bearing mice were treated with radiation (RT, 2 Gy) using 6 MV photon beams. RT treated mice were euthanized at 0 hr, 6 hrs and 24 hrs. The tumors were processed for immunohistochemical (IHC) staining for PARP (a surrogate of apoptosis). A multispectral imaging analysis system was used to quantify the expression of PARP staining. Cell apoptosis was calculated based on the PARP expression level using the DAB analysis software. Results: Our data showed that PARP related apoptosis peaked at 48 hrs and 7 days in pHIFU treated mice, which is comparable to that for the RT group at 24 hrs. The preliminary results from this study were consistent with our previous study on tumor growth delay using pHIFU. Conclusion: Our results demonstrated that non-thermal pHIFU increased apoptotic tumor cell death through the PARP related pathway. MR guided pHIFU may have a great potential as a safe, noninvasive treatment modality for cancer therapy. This treatment modality may synergize with PARP inhibitors to achieve better therapeutic result.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-80: A Mobile Application for Intra-Operative Electron Radiotherapy Treatment Planning

Purpose: Intraoperative electron radiotherapy (IORT) poses a unique set of challenges for treatment planning. Planning must be performed in a busy operating room environment over a short timeframe often with little advance knowledge of the treatment depth or applicator size. Furthermore, IORT accelerators can have a large number of possible applicators, requiring extensive databooks that must be searched for the appropriate dosimetric parameters. The goal of this work is to develop a software tool to assist in the planning process that is suited to the challenges faced in the IORT environment. Methods: We developed a mobile application using HTML5 and Javascript that can be deployed to tablet devices suitable for use in the operating room. The user selects the desired treatment parameters cone diameter, bevel angle, and energy (a total of 141 datasets) and desired bolus. The application generates an interactive display that allows the user to dynamically select points on the depth-dose curve and to visualize the shape of the corresponding isodose contours. The user can indicate a prescription isodose line or depth. The software performs a monitor unit calculation and generates a PDF report. Results: We present our application, which is now used routinely in our IORT practice. It has been employed successfully in over 23 cases. The interactivity of the isodose distributions was found to be of particular use to physicians who are less-frequent IORT users, as well as for the education of residents and trainees. Conclusion: This software has served as a useful tool in IORT planning, and demonstrates the need for treatment planning tools that are designed for the specialized challenges encountered in IORT. This software is the subject of a license agreement with the IntraOp Medical Corporation. This software is the subject of a license agreement between Massachusetts General Hospital / Partners Healthcare and the IntraOp Medical Corporation. CLW is consulting on software development with the IntraOp Medical Corporation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-BRC-04: Collimator Angle Optimization in Volumetric Modulated Arc Therapy

Purpose: Volumetric modulated arc therapy (VMAT) has revolutionized radiation treatment by decreasing treatment time and monitor units, thus reducing scattered and whole body radiation dose. As the collimator angle changes the apparent leaf gap becomes larger which can impact plan quality, organ at risk (OAR) sparing as well as IMRT QA passing rate which is investigated. Methods: Two sites (prostate and head and neck) that have maximum utilization of VMAT were investigated. Two previously treated VMAT patients were chosen. For each patient 10 plans were created by maintaining constant optimization constraints while varying collimator angles from 0-90 deg at an interval of 10 degrees for the first arc and the appropriate complimentary angle for the second arc. Plans were created with AAA algorithm using 6 MV beam on a Varian IX machine with Millennium 120 MLC. The dose-volume histogram (DVH) for each plan was exported and dosimetric parameters (D98, D95, D50, D2) as well homogeneity index (HI) and conformity index (CI) were computed. Each plan was validated for QA using ArcCheck with gamma index passing criteria of 2%/2 mm and 3%/3 mm. Additionally, normal tissue complication probability (NTCP) for each OAR was computed using Uzan-Nahum software. Results: The CI values for both sites had no impact as target volume coverage in every collimator angle were the same since it was optimized for adequate coverage. The HI which is representative of DVH gradient or dose uniformity in PTV showed a clear trend in both sites. The NTCP for OAR (brain and cochlea) in H&N plan and (bladder and rectum) in prostate plan showed a distinct superiority for collimator angles between 15-30 deg. The gamma passing rates were not correlated with angle. Conclusion: Based on CI, HI, NTCP and gamma passing index, it can be concluded that collimator angles should be maintained within 15–30 deg.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Analysis Code - Data Analysis in 'Leveraging Multiple Statistical Methods for Inverse Prediction in Nuclear Forensics Applications' (LMSMIPNFA) v. 1.0

R code that performs the analysis of a data set presented in the paper ‘Leveraging Multiple Statistical Methods for Inverse Prediction in Nuclear Forensics Applications’ by Lewis, J., Zhang, A., Anderson-Cook, C. It provides functions for doing inverse predictions in this setting using several different statistical methods. The data set is a publicly available data set from a historical Plutonium production experiment.

Lewis, John↗
DOE OSTI2018

Technical Manual and User's Guide for MILDOS, Version 4.1

The MILDOS computer code is used to estimate the radiological impacts of airborne emissions from uranium mining and milling facilities. The code allows users to consider both conventional uranium ore operations and operations associated with in situ recovery (ISR) facilities. The code is used by license applicants and U.S. Nuclear Regulatory Commission staff to perform routine radiological impact and compliance evaluations for various uranium recovery operations. As part of a major code update, the following improvements and additions were made to the code: (a) revised ISR models, (b) support for ores containing thorium-232 (Th-232) and its daughter radionuclides in addition to the currently supported uranium-238 (U-238) and its daughter radionuclides, (c) a revised area source model, (d) the capability to perform sensitivity analysis on specific input parameters, (e) the capability to use current meteorological data provided by the National Climatic Data Center, and (f) an interactive results module.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

MIRC-X: A Highly Sensitive Six-telescope Interferometric Imager at the CHARA Array

Michigan InfraRed Combiner-eXeter (MIRC-X) is a new highly sensitive six-telescope interferometric imager installed at the CHARA Array that provides an angular resolution equivalent of up to a 330 m diameter baseline telescope in J- and H-band wavelengths ((λ/(2B))∼0.6 mas). We upgraded the original Michigan InfraRed Combiner (MIRC) instrument to improve sensitivity and wavelength coverage in two phases. First, a revolutionary sub-electron noise and fast-frame-rate C-RED ONE camera based on an SAPHIRA detector was installed. Second, a new-generation beam combiner was designed and commissioned to (i) maximize sensitivity, (ii) extend the wavelength coverage to J band, and (iii) enable polarization observations. A low-latency and fast-frame-rate control software enables high-efficiency observations and fringe tracking for the forthcoming instruments of the CHARA Array. Since mid-2017, MIRC-X has been offered to the community and has demonstrated best-case H-band sensitivity down to 8.2 correlated magnitude. MIRC-X uses single-mode fibers to coherently combine the light from six telescopes simultaneously with an image-plane combination scheme and delivers a visibility precision better than 1%, and closure phase precision better than 1°. MIRC-X aims at (i) imaging protoplanetary disks, (ii) detecting exoplanets with precise astrometry, and (iii) imaging stellar surfaces and starspots at an unprecedented angular resolution in the near-infrared. In this paper, we present the instrument design, installation, operation, and on-sky results, and demonstrate the imaging capability of MIRC-X on the binary system ι Peg. The purpose of this paper is to provide a solid reference for studies based on MIRC-X data and to inspire future instruments in optical interferometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A 78 day X-ray period detected from NGC 5907 ULX1 by Swift

We report the detection of a 78.1 ± 0.5 day period in the X-ray light curve of the extreme ultraluminous X-ray source NGC 5907 ULX1 (L{sub X,peak}∼5×10{sup 40} erg s{sup −1}), discovered during an extensive monitoring program with Swift. These periodic variations are strong, with the observed flux changing by a factor of ∼3–4 between the peaks and the troughs of the cycle; our simulations suggest that the observed periodicity is detected comfortably in excess of 3σ significance. We discuss possible origins for this X-ray period, but conclude that at the current time we cannot robustly distinguish between orbital and super-orbital variations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A novel tool for improving the data collection process during control room modernization human-system interface testing and evaluation activities - 271

The U.S. Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) Program is developing a scientific basis to extend the existing U.S. nuclear power plant (NPP) operating life beyond the current 60-year licensing period and to ensure NPP long-term reliability, productivity, safety, and security. Under the Advanced Instrumentation, Information, and Control (II and C) Systems Technologies pathway, the LWRS Program conducts targeted research and development to address long-term aging and obsolescence of existing instrumentation and control technologies, as well as to implement and test new digital II and C technology that enables broad innovation and business improvement in the NPP operating model. To this end, the role of human factors engineering (HFE) is important as part of ensuring that the control room modernization (CRM) activities account for the operator's cognitive and physiological characteristics. Part of the HFE process for CRM entails conducting tests and evaluations (T and Es) that will inform the design of new human-system interface (HSI) technologies. Qualitative data that captures operator feedback and important human-system interactions are an important part of T and E. However, the analyzing this qualitative data in order to provide actionable design recommendations can be time-consuming and error prone. This paper presents an early version of a tool, which uses a combination of natural language processing and supervised machine learning to support the analysis of qualitative data. The objectives of this tool is to improve the process of (1) sifting through freeform text for improved design issue identification, (2) maintaining consistency in theme coding conventions, and (3) readily providing actionable design guidance that can be traceable to state-of-the-art HFE guidelines such as NUREG-0700. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

MCNP ® Code Version 6.3.1 Theory & User Manual

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.1, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

High Speed Tribology of the Head-Disk Interface by Computer Simulation Final Report CRADA No. TC-297-92-C

The purpose of this project is to develop and validate a molecular dynamics model of nanoscale tribological processes occurring inside the head-disk interface of magnetic recording devices. Within the three years of the project the partners will: write a new molecular dynamics computer simulation code for growth, sliding and impact inside the head-disk interface; identify accurate models for chemical bonding within coatings used to protect magnetic surfaces, especially amorphous carbon thin films; simulate the growth of amorphous carbon thin films and analyse the resulting chemical bonding; simulate friction and wear processes relevant to the head-disk interface; and validate the model through comparison with AFM and nanoindenter studies of friction and wear of amorphous carbon thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The ocean model for E3SM global applications: Omega version 0.1.0 – a new high-performance computing code for exascale architectures

This paper introduces Omega, the Ocean Model for E3SM Global Applications. Omega is a new ocean model designed to run efficiently on high performance computing (HPC) platforms, including exascale heterogeneous architectures with accelerators, such as Graphics Processing Units (GPUs). Omega is written in C and uses the Kokkos performance portability library. These were chosen because they are well-supported and will help future-proof Omega for upcoming HPC architectures. Omega will eventually replace the Model for Prediction Across Scales-Ocean (MPAS-Ocean) in the US Department of Energy's (DOE's) Energy Exascale Earth System Model (E3SM). Omega runs on unstructured horizontal meshes with variable-resolution capability and implements the same horizontal discretization as MPAS-Ocean. This work documents the design and performance of Omega Version 0.1.0 (Omega-V0), which solves the shallow water equations with passive tracers and is the first step towards the full primitive equation ocean model. On Central Processing Units (CPUs), Omega-V0 is 1.4 times faster than MPAS-Ocean with the same configuration. Omega-V0 is more efficient on GPUs than CPUs on a per-watt basis – by a factor of 5.3 on Frontier and 3.6 on Aurora, two of the world's fastest exascale computers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Measurement of deuterium density profiles in the H-mode steep gradient region using charge exchange recombination spectroscopy on DIII-D

Recent completion of a thirty two channel main-ion (deuterium) charge exchange recombination spectroscopy (CER) diagnostic on the DIII-D tokamak enables detailed comparisons between impurity and main-ion temperature, density, and toroidal rotation. In a H-mode DIII-D discharge, these new measurement capabilities are used to provide the deuterium density profile, demonstrate the importance of profile alignment between Thomson scattering and CER diagnostics, and aid in determining the electron temperature at the separatrix. Sixteen sightlines cover the core of the plasma and another sixteen are densely packed towards the plasma edge, providing high resolution measurements across the pedestal and steep gradient region in H-mode plasmas. Extracting useful physical quantities such as deuterium density is challenging due to multiple photoemission processes. Finally, these challenges are overcome using a detailed fitting model and by forward modeling the photoemission using the FIDASIM code, which implements a comprehensive collisional radiative model. Published by AIP Publishing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

The influence of Hall physics on power-flow along a coaxial transmission line

Extended-MHD simulations of a coaxial transmission line are performed in axisymmetric cylindrical geometry, in particular, in examining the influence of Hall physics on a plasma layer initialized against the anode versus the cathode, for which an MHD model is insensitive. The results indicate that Hall physics is required in order to model an electron E × B drift current in the electrode plasma, which is parallel to the anode current and opposite the cathode current. This results in confinement of the electrode plasma when initialized against the cathode and expansion of the plasma layer when initialized against the anode. The expansion in the anode-initialized case results in filaments of plasma bridging the gap, causing substantial power-flow losses. These results represent the first fluid simulations of power-flow, to our knowledge, that, by including Hall physics, recover fundamental aspects of anode and cathode dynamics predicted by kinetic theory while simulating over a dynamic range (nine orders of magnitude density variation from solid-density electrodes down to low-density electrode plasma) which is prohibitive for Particle-In-Cell (PIC) codes. In conclusion, this work demonstrates the need for further development of extended-MHD and two-fluid modeling of power-flow dynamics, which, possibly through hybridization with a PIC code, will eventually culminate in a code with reliable predictive capability for power-flow coupling and energy losses in pulsed-power systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Linear response of a Hall magnetic drift wave for verification of Hall MHD algorithms

Numerical implementations of Hall magnetohydrodynamics (Hall MHD) can be challenging due to the nonlinear multidimensional nature of the Hall term. Here, a model problem is presented that couples the hydrodynamic motion of the plasma to Hall MHD evolution of the magnetic field. The Hall MHD equations are linearized about unperturbed solutions in both cylindrical and Cartesian coordinates in two dimensions. The magnetic field is assumed to lie in the ignorable direction, and the linear response about the unperturbed solution is considered. The resulting ordinary differential equation is used to numerically compute the eigenfunctions and eigenfrequencies of the mode. The resulting eigenfunctions do not make the local wave approximation but are instead global solutions that depend on the spatial dependence of the unperturbed Alfvén speed. Hall MHD simulations are then performed in the Ares multiphysics code and shown to agree with the predicted phase velocities of the wave, and the simulated solutions are shown to numerically converge to the semianalytic modes. By varying the background density of the plasma (and correspondingly, the ion inertial length), the importance of Hall physics can be varied. This allows the test problem to transition from the classical MHD limit to the extreme Hall MHD limit. Furthermore, this problem is a useful tool for the verification of Hall MHD routines implemented in various codes, and the robustness of a routine can be tested in regimes in which Hall physics is dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulation of divertor heat flux width on EAST by BOUT++ transport code

The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Clusterdv: a simple density-based clustering method that is robust, general and automatic

Abstract Motivation How to partition a dataset into a set of distinct clusters is a ubiquitous and challenging problem. The fact that data vary widely in features such as cluster shape, cluster number, density distribution, background noise, outliers and degree of overlap, makes it difficult to find a single algorithm that can be broadly applied. One recent method, clusterdp, based on search of density peaks, can be applied successfully to cluster many kinds of data, but it is not fully automatic, and fails on some simple data distributions. Results We propose an alternative approach, clusterdv, which estimates density dips between points, and allows robust determination of cluster number and distribution across a wide range of data, without any manual parameter adjustment. We show that this method is able to solve a range of synthetic and experimental datasets, where the underlying structure is known, and identifies consistent and meaningful clusters in new behavioral data. Availability and implementation The clusterdv is implemented in Matlab. Its source code, together with example datasets are available on: https://github.com/jcbmarques/clusterdv. Supplementary information Supplementary data are available at Bioinformatics online.

Marques, João C.↗
DOE OSTI2019

Search for α decay of Te 104 with a novel recoil-decay scintillation detector

A search for superallowed α decay of N = Z nuclei Te 104 and Xe 108 was carried out using a novel recoil-decay scintillator detector at the tandem accelerator facility at the Japan Atomic Energy Agency (JAEA). In this study, inorganic crystal scintillation material YAP:Ce (yttrium aluminum perovskite) coupled to a position-sensitive photomultiplier tube (PSPMT) was implemented for the first time in a radioactive decay experiment. Residues from the fusion-evaporation reaction Ni 58 + Fe 54 → Xe * 112 were separated by the JAEA Recoil Mass Separator (RMS) and implanted into the YAP:Ce crystal. α decays of neutron-deficient tellurium isotopes were identified and proton emission of I 109 was observed. The α -decay chain Xe 109 → Te 105 → Sn 101 was recorded with a time interval of 960 ns between two α pulses. Position localization in the crystal for decays and ions in the energy range from hundreds of keV to 60 MeV was achieved with an accuracy of 0.67 mm, proving that this detector is capable of making temporal and spatial correlations for fast decay events. No conclusive evidence was found for the decay chain Xe 108 → Te 104 → Sn 100 within 3 days of experiment. However, two events were observed with properties consistent with the reported observation at the Fragment Mass Analyzer (FMA), but with a separation between signals of less than 4 ns. The cross section limit of 130 pb was obtained for production of two events of Xe 108 , about an order of magnitude below the expectation based on earlier cross section measurements and the HIVAP fusion-evaporation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗