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,549 records · Page 86

Characterization of a 1D-imaging high-energy x-ray backlighter driven by the National Ignition Facility Advanced Radiographic Capability laser

Plastic deformation of samples compressed to Mbar pressures at high strain rates at the National Ignition Facility (NIF) forms the basis of ongoing material strength experiments in conditions relevant to meteor impacts, geophysics, armor development, and inertial confinement fusion. Hard x-ray radiography is the primary means of measuring the evolution of these samples, typically employing a slit-collimated high-Z microdot driven by the NIF laser to generate >40 keV x rays [E. Gumbrell et al., Rev. Sci. Instrum. 89, 10G118 (2018) and C. M. Huntington et al., Rev. Sci. Instrum. 89, 10G121 (2018)]. Alternatively, a dysprosium “micro-flag” target driven by the Advanced Radiographic Capability laser (~2 kJ, 10 ps) can deliver significantly higher spatiotemporal resolution [M. P. Hill et al., Rev. Sci. Instrum. 92, 033535 (2021)], especially in high-opacity samples. Here, initial experiments revealed problematic brightness and spectral gradients from this source, but by radiographing a set of diamond-turned, 105 µm-thick Pb test objects and supported by simulations using the 3D Monte Carlo code GEANT4, these geometry-dependent gradients across the field of view are quantified and mitigation strategies are assessed. In addition to significantly enhancing the modulation transfer function compared to the existing system, image stacking from multiple layers of image plate is shown to almost double the signal to noise ratio that will reduce uncertainties in future dynamic strength experiments.

47 OTHER INSTRUMENTATION↗

Adverse Events Involving Radiation Oncology Medical Devices: Comprehensive Analysis of US Food and Drug Administration Data, 1991 to 2015

Purpose: Radiation oncology relies on rapidly evolving technology and highly complex processes. The US Food and Drug Administration collects reports of adverse events related to medical devices. We sought to characterize all events involving radiation oncology devices (RODs) from the US Food and Drug Administration's postmarket surveillance Manufacturer and User Facility Device Experience (MAUDE) database, comparing these with non–radiation oncology devices. Methods and Materials: MAUDE data on RODs from 1991 to 2015 were sorted into 4 product categories (external beam, brachytherapy, planning systems, and simulation systems) and 5 device problem categories (software, mechanical, electrical, user error, and dose delivery impact). Outcomes included whether the device was evaluated by the manufacturer, adverse event type, remedial action, problem code, device age, and time since 510(k) approval. Descriptive statistics were performed with linear regression of time-series data. Results for RODs were compared with those for other devices by the Pearson χ{sup 2} test for categorical data and 2-sample Kolmogorov-Smirnov test for distributions. Results: There were 4234 ROD and 4,985,698 other device adverse event reports. Adverse event reports increased over time, and events involving RODs peaked in 2011. Most ROD reports involved external beam therapy (50.8%), followed by brachytherapy (24.9%) and treatment planning systems (21.6%). The top problem types were software (30.4%), mechanical (20.9%), and user error (20.4%). RODs differed significantly from other devices in each outcome (P<.001). RODs were more likely to be evaluated by the manufacturer after an event (46.9% vs 33.0%) but less likely to be recalled (10.5% vs 37.9%) (P<.001). Device age and time since 510(k) approval were shorter among RODs (P<.001). Conclusions: Compared with other devices, RODs may experience adverse events sooner after manufacture and market approval. Close postmarket surveillance, improved software design, and manufacturer–user training may help mitigate these events.

07 ISOTOPE AND RADIATION SOURCES↗

Direct Estimation of Parameters in ODE Models Using WENDy: Weak-Form Estimation of Nonlinear Dynamics

Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of $$C^{\infty }$$ C ∞ bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at https://github.com/MathBioCU/WENDy .

97 MATHEMATICS AND COMPUTING↗

Chromosome-level genome assembly and transcriptome of the green alga Chromochloris zofingiensis illuminates astaxanthin production

Microalgae have potential to help meet energy and food demands without exacerbating environmental problems. There is interest in the unicellular green alga Chromochloris zofingiensis, because it produces lipids for biofuels and a highly valuable carotenoid nutraceutical, astaxanthin. Here, to advance understanding of its biology and facilitate commercial development, we present a C. zofingiensis chromosome-level nuclear genome, organelle genomes, and transcriptome from diverse growth conditions. The assembly, derived from a combination of short- and long-read sequencing in conjunction with optical mapping, revealed a compact genome of ~58 Mbp distributed over 19 chromosomes containing 15,274 predicted protein-coding genes. The genome has uniform gene density over chromosomes, low repetitive sequence content (~6%), and a high fraction of protein-coding sequence (~39%) with relatively long coding exons and few coding introns. Functional annotation of gene models identified orthologous families for the majority (~73%) of genes. Synteny analysis uncovered localized but scrambled blocks of genes in putative orthologous relationships with other green algae. Two genes encoding beta-ketolase (BKT), the key enzyme synthesizing astaxanthin, were found in the genome, and both were up-regulated by high light. Isolation and molecular analysis of astaxanthin-deficient mutants showed that BKT1 is required for the production of astaxanthin. Moreover, the transcriptome under high light exposure revealed candidate genes that could be involved in critical yet missing steps of astaxanthin biosynthesis, including ABC transporters, cytochrome P450 enzymes, and an acyltransferase. Finally, the high-quality genome and transcriptome provide insight into the green algal lineage and carotenoid production.

60 APPLIED LIFE SCIENCES↗

The Sensitivity of GPz Estimates of Photo-z Posterior PDFs to Realistically Complex Training Set Imperfections

The accurate estimation of photometric redshifts is crucial to many upcoming galaxy surveys, for example, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Almost all Rubin extragalactic and cosmological science requires accurate and precise calculation of photometric redshifts; many diverse approaches to this problem are currently in the process of being developed, validated, and tested. In this work, we use the photometric redshift code GPz to examine two realistically complex training set imperfections scenarios for machine learning based photometric redshift calculation: (i) where the spectroscopic training set has a very different distribution in color–magnitude space to the test set, and (ii) where the effect of emission line confusion causes a fraction of the training spectroscopic sample to not have the true redshift. By evaluating the sensitivity of GPz to a range of increasingly severe imperfections, with a range of metrics (both of photo-z point estimates as well as posterior probability distribution functions, PDFs), we quantify the degree to which predictions get worse with higher degrees of degradation. In particular, we find that there is a substantial drop-off in photo-z quality when line-confusion goes above ~1%, and sample incompleteness below a redshift of 1.5, for an experimental setup using data from the Buzzard Flock synthetic sky catalogs.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-F-SPS-03: Direct Measurement of Organ Doses Resulting From Head and Cervical Spine Trauma CT Protocols

Purpose: This retrospective study analyzes the exposure history of emergency department (ED) patients undergoing head and cervical spine trauma computed tomography (CT) studies. This study investigated dose levels received by trauma patients and addressed any potential concerns regarding radiation dose issues. Methods: Under proper IRB approval, a cohort of 300 trauma cases of head and cervical spine trauma CT scans received in the ED was studied. The radiological image viewing software of the hospital was used to view patient images and image data. The following parameters were extracted: the imaging history of patients, the reported dose metrics from the scanner including the volumetric CT Dose Index (CTDIvol) and Dose Length Product (DLP). A postmortem subject was scanned using the same scan techniques utilized in a standard clinical head and cervical spine trauma CT protocol with 120 kVp and 280 mAs. The CTDIvol was recorded for the subject and the organ doses were measured using optically stimulated luminescent (OSL) dosimeters. Typical organ doses to the brain, thyroid, lens, salivary glands, and skin, based on the cadaver studies, were then calculated and reported for the cohort. Results: The CTDIvol reported by the CT scanner was 25.5 mGy for the postmortem subject. The average CTDIvol from the patient cohort was 34.1 mGy. From these metrics, typical average organ doses in mGy were found to be: Brain (44.57), Thyroid (33.40), Lens (82.45), Salivary Glands (61.29), Skin (47.50). The imaging history of the cohort showed that on average trauma patients received 26.1 scans over a lifetime. Conclusion: The average number of scans received on average by trauma ED patients shows that radiation doses in trauma patients may be a concern. Available dose tracking software would be helpful to track doses in trauma ED patients, highlighting the importance of minimizing unnecessary scans and keeping doses ALARA.

60 APPLIED LIFE SCIENCES↗

Multispecies Hypersonic Flow Simulator

MHYSA is a conservative finite-difference code to solve the 1D, 2D, and 3D multispecies Navier-Stokes equations for hypersonic flow applications. It supports Cartesian grids only and has the option of including immersed boundaries. - Solves 1D Euler and 2D and 3D Navier-Stokes equations. - Solves the PDEs over Cartesian grids. - Written entirely in C and uses the MPICH library. It also uses OpenMP threads but this is a work-inprogress. - Can be compiled with PETSc (http://www.mcs.anl.gov/petsc/), if available, where it can use PETSc's time integration module TS (http://www.mcs.anl.gov/petsc/petsc-current/src/ts/). - For 3-dimensional simulations, the immersed boundaries can be used to solve over non-Cartesian geometries. MHYSA has been developed to be scalable, and apart from the usual functionalities to solve a system of PDEs on distributed architectures, it provides scalable file I/O functions.

Ghosh, Debojyoti↗

SPALLMAP

Boiler tubes in steam power plants experience exfoliation of oxide grown on the inner side of the tubes. In extreme cases, the exfoliation cause significant tube blockages that lead to forced power plant outages. It is thus desired to predict through modeling the propensity of exfoliation events in order to inform power plant operators of possible tube blockages. SpallMap solves for the stress-strain equations in an axisymmetric geometry, tracking the stress/strain evolution during boiler operation including outages at one-location along a boiler tube and compares it with scale damage criteria represented by Armitt diagram. The SPALLMAP code contains modules developed for oxide growth, stress analysis, and classical fracture mechanics correlations by taking into account the following phenomena and features, (a) Non-uniform thermal expansion coefficient of oxides and metal substrates, (b) Plant operation schedule with periodic alternate full-load and partial-load regimes, (c) axisymmetric formulation for cylindrical tubes, (d) Multiple oxide layers, (e) oxide-growth induced stresses, and (f) damage criteria from classical fracture mechanics. The computer program is written in FORTRAN90. Its modular structure was sought for allowing the best flexibility in updating the program by implementing new constitutive equations due to availability of new material property data and/or new physical phenomena.

Sabau, Adrian↗

SISSO (Sure-independence-screening sparsifying-operator regressor)

Implementation of the SISSO regression algorithm in MATLAB. The SISSO regression algorithm iteratively selects model features from candidates, converging even when the number of possible features is much greater than the number of available data points. Includes a test script that validates the algorithm by replicating the results using data and procedure from https://analytics-toolkit.nomad-coe.eu/hub/user-redirect/notebooks/tutorials/compressed_sensing.ipynb. Python code from the SISSO regressor in 'sisso.py', from the link above, was used as the basis for developing the MATLAB implementation. The SISSO regression algorithm is detailed by the original authors in R. Ouyang, S. Curtarolo, E. Ahmetcik et al., Phys. Rev. Mater. 2, 083802 (2018), R. Ouyang, E. Ahmetcik, C. Carbogno, M. Scheffler, and L. M. Ghiringhelli, J. Phys.: Mater. 2, 024002 (2019).

Gasper, Paul↗

Gene targets for engineering osmotolerance in Caldicellulosiruptor bescii

Abstract Background Caldicellulosiruptor bescii , a promising biocatalyst being developed for use in consolidated bioprocessing of lignocellulosic materials to ethanol, grows poorly and has reduced conversion at elevated medium osmolarities. Increasing tolerance to elevated fermentation osmolarities is desired to enable performance necessary of a consolidated bioprocessing (CBP) biocatalyst. Results Two strains of C. bescii showing growth phenotypes in elevated osmolarity conditions were identified. The first strain, ORCB001, carried a deletion of the FapR fatty acid biosynthesis and malonyl-CoA metabolism repressor and had a severe growth defect when grown in high-osmolarity conditions—introduced as the addition of either ethanol, NaCl, glycerol, or glucose to growth media. The second strain, ORCB002, displayed a growth rate over three times higher than its genetic parent when grown in high-osmolarity medium. Unexpectedly, a genetic complement ORCB002 exhibited improved growth, failing to revert the observed phenotype, and suggesting that mutations other than the deleted transcription factor (the fruR/cra gene) are responsible for the growth phenotype observed in ORCB002. Genome resequencing identified several other genomic alterations (three deleted regions, three substitution mutations, one silent mutation, and one frameshift mutation), which may be responsible for the observed increase in osmolarity tolerance in the fruR / cra -deficient strain, including a substitution mutation in dnaK , a gene previously implicated in osmoresistance in bacteria. Differential expression analysis and transcription factor binding site inference indicates that FapR negatively regulates malonyl-CoA and fatty acid biosynthesis, as it does in many other bacteria. FruR/Cra regulates neighboring fructose metabolism genes, as well as other genes in global manner. Conclusions Two systems able to effect tolerance to elevated osmolarities in C. bescii are identified. The first is fatty acid biosynthesis. The other is likely the result of one or more unintended, secondary mutations present in another transcription factor deletion strain. Though the locus/loci and mechanism(s) responsible remain unknown, candidate mutations are identified, including a mutation in the dnaK chaperone coding sequence. These results illustrate both the promise of targeted regulatory manipulation for osmotolerance (in the case of fapR ) and the challenges (in the case of fruR/cra ).

59 BASIC BIOLOGICAL SCIENCES↗

Singe Library Review: XCAP Review Report

Singe is a computational library for simulating thermonuclear reactions written in C++14. It was inspired by the earlier Singe library, developed by Joann Campbell. From the beginning of the Singe project, the goal has been to create a library that provides the necessary functionality, runs efficiently, and is flexible enough to connect to a variety of multi-physics host codes. Such a library allows subject-matter experts to focus on development within the library, while library staff can support host-code staff in integrating Singe into their multi-physics codes. A common library used across multiple multi-physics codes allows those subject-matter experts to roll out new features to all host codes simultaneously, and eases the process of code-to-code comparison between codes that are leveraging the same physics in different ways. The Singe effort was established by Joann Campbell, who provided significant guidance on the goals, requirements, and philosophy of the library. It has been developed by Brendan K. Krueger, CJ Solomon, Jesse Canfield, and Terry Adams, with computer science support from Tyler Reddy and Karen Tsai. Discussions with multiple potential stakeholders have been critical to the planning of such a cross-cutting project, including representatives from multiple host codes and Brian Haines who helped define usage and needs for the inertial confinement fusion (ICF) program. This document is based on the XCAP review process. It is an updated version of the initial document distributed to the reviewers, memo XCP-2:21-005. We want to acknowledge the time and effort of the reviewers who participated in this process: Philipp Edelmann, Alan Harrison, Christopher Malone, Edward Norris, Ryosuke Park, Galen Shipman, CJ Solomon, and any others who may have contributed to the feedback that was provided to us.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmark Exercise for the Control Rod Swelling Evaluation

The VTR core has six reactivity control assemblies and three safety assemblies. The control assemblies or primary control rods are adjusted during the normal operation to balance the core reactivity and to control the reactor power. A typical control assembly radial layout is presented in Figure 1. The figure shows the swelled absorber (B 4 C) rod. Initially, helium gas fills the gap between the pin and the cladding before irradiation swelling takes place. For VTR, HT9 steel was selected as the cladding and duct material. The main neutron absorbing material used in the VTR is B 4 C. When residing in the core, the neutronics, thermophysical, and mechanical properties of the materials used in a control assembly will degrade due to accumulated neutron damage. Material degradation limits how long a control assembly can reside in the core. Many phenomena affect the control assembly lifetime, such as the loss of reactivity worth due to B 4 C depletion, the mechanical interaction of the absorber rod and the cladding due to B 4 C swelling, the helium gas buildup in the pin due to B-10 capture, etc. B 4 C swelling, which causes closure of the gap between the absorber rod and the cladding, is usually considered as the main limiting factor from past experience. An initial study was conducted at PNNL to evaluate the irradiation behavior of a VTR control assembly. The evaluation was performed using the CNRD2 code that was initially developed for the FFTF. The study also included an assessment of the VTR control assembly and focused on a 61-pin control assembly design, which is different from that used (37-pin design) in the core design study. The study conducted by PNNL was reviewed independently by ANL. A Python script referred to as the Control Assembly Evaluation Script (CAES) was developed for the independent review and additional assessment of 37-pin control assembly design. The script has focused on the assessment of the absorber rod swelling for its importance in determining the control assembly lifetime. CAES uses geometry, neutronics, materials data as input to predict the swelling of the absorber rod during its residence in the reactor core. The results from CAES showed some non-negligible differences against the PNNL results. Some of the differences can be attributed to the different interpretation of the control rod assembly dimensions. To resolve this issue, a benchmark exercise was proposed. The benchmark specification was developed by PNNL. The benchmark exercise was performed independently at PNNL and ANL using different codes/scripts (CRND2 and CAES). This memo documents the results calculated using the different codes. However, this report is limited to presenting the results obtained. Further investigation of the cause of the observed difference will be performed as part of future activities, pending continuation of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of dispensing hardware for safe fueling of heavy duty vehicles

The development of safe dispensing equipment for the fueling of heavy duty (HD) vehicles is critical to the expansion of this newly and quickly expanding market. This paper discusses the development of a HD dispenser and nozzles assembly (nozzle, hose, breakaway) for these new, larger vehicles where flow rates are more than double compared to light duty (LD) vehicles. This equipment must operate at nominal pressures of 700 bar, -40o C gas temperature, and average flow rate of 5-10 kg/min at a high throughput commercial hydrogen fueling station without leaking hydrogen. The project surveyed HD vehicle manufacturers, station developers, and component suppliers to determine the basic specifications of the dispensing equipment and nozzle assembly. The team also examined existing codes and standards to determine necessary changes to accommodate HD components. From this information, the team developed a set of specifications which will be used to design the dispensing equipment. In order to meet these goals, the team performed computational fluid dynamic, pressure modelling, and temperature analysis in order to determine the necessary parameters to meet existing safety standards modified for HD fueling. The team also considered user, operational, and maintenance requirements, such as freeze lock which has been an issue which prevents the removal of the nozzle from LD vehicles. The team also performed a failure mode and effects analysis (FMEA) to identify the possible failures in the design. The dispenser and nozzle assembly will be tested separately, and then installed on an innovative, HD fueling station which will use a HD vehicle simulator to test the entire system.

08 HYDROGEN↗

Status of the CERBERUS Evaluation for the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook

Modeling & Simulation (M&S) tools are used to analyze advanced reactor designs and the safety of current nuclear operations. As computers continue to improve, we are able to enhance resolution in our calculations. Therefore, the limitations of simulation capability are in the quality of data that is being used, including our ability to quantify the uncertainty and sensitivity of that data. In order to model systems of interest with increasing accuracy, the industry must improve key nuclear data measurements. The International Criticality Safety Benchmark Evaluation Project (ICSBEP) compiles and evaluates experiment data in a handbook that can be used by criticality safety engineers and others to validate computer codes and cross section libraries at nuclear facilities. Both critical and subcritical experiments are included in the handbook. These experiments, along with differential measurements, can help improve the quality of nuclear data. Concerns regarding the accuracy of Cu nuclear data have been published. The large values and trend of C-E for the Zeus intermediate energy benchmark, being one of the primary examples. Furthermore, very few experiments have been designed to be sensitive to Cu (as shown in Figure 1), so an integral, critical experiment is needed to help resolve these differences.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-C-201-06: Small Field Correction Factors for the MicroDiamond Detector in the Gamma Knife-Model C Derived Using Monte Carlo Methods

Purpose: To determine small field correction factors for PTW’s microDiamond detector in Elekta’s Gamma Knife Model-C unit. These factors allow the microDiamond to be used in QA measurements of output factors in the Gamma Knife Model-C; additionally, the results also contribute to the discussion on the water equivalence of the relatively-new microDiamond detector and its overall effectiveness in small field applications. Methods: The small field correction factors were calculated as k correction factors according to the Alfonso formalism. An MC model of the Gamma Knife and microDiamond was built with the EGSnrc code system, using BEAMnrc and DOSRZnrc user codes. Validation of the model was accomplished by simulating field output factors and measurement ratios for an available ABS plastic phantom and then comparing simulated results to film measurements, detector measurements, and treatment planning system (TPS) data. Once validated, the final k factors were determined by applying the model to a more waterlike solid water phantom. Results: During validation, all MC methods agreed with experiment within the stated uncertainties: MC determined field output factors agreed within 0.6% of the TPS and 1.4% of film; and MC simulated measurement ratios matched physically measured ratios within 1%. The final k correction factors for the PTW microDiamond in the solid water phantom approached unity to within 0.4%±1.7% for all the helmet sizes except the 4 mm; the 4 mm helmet size over-responded by 3.2%±1.7%, resulting in a k factor of 0.969. Conclusion: Similar to what has been found in the Gamma Knife Perfexion, the PTW microDiamond requires little to no corrections except for the smallest 4 mm field. The over-response can be corrected via the Alfonso formalism using the correction factors determined in this work. Using the MC calculated correction factors, the PTW microDiamond detector is an effective dosimeter in all available helmet sizes. The authors would like to thank PTW (Friedberg, Germany) for providing the PTW microDiamond detector for this research.

07 ISOTOPE AND RADIATION SOURCES↗

Multiscale modeling of fission gas behavior in U 3 Si 2 under LWR conditions

In this work, we present a model of fission gas behavior in U 3 Si 2 under light water reactor (LWR) conditions for application in engineering fuel performance codes. The model includes components for intra-granular and inter-granular behavior of fission gases. The intra-granular component is based on cluster dynamics and computes the evolution of intra-granular fission gas bubbles and swelling coupled to gas diffusion to grain boundaries. The inter-granular component describes the evolution of grain-boundary fission gas bubbles coupled to fission gas release. Given the lack of experimental data for U 3 Si 2 under LWR conditions, the model is informed with parameters calculated via atomistic simulations. In particular, we derive fission gas diffusivities through density functional theory calculations, and the re-solution rate of fission gas atoms from intra-granular bubbles through binary collision approximation calculations. Further, the developed model is applied to the simulation of an experiment for U 3 Si 2 irradiated under LWR conditions available from the literature. Results point out a credible representation of fission gas swelling and release in U 3 Si 2 . Finally, we perform a sensitivity analysis for the various model parameters. Based on the sensitivity analysis, indications are derived that can help in addressing future research on the characterization of the physical parameters relative to fission gas behavior in U 3 Si 2 . The developed model is intended to provide a suitable infrastructure for the engineering scale calculation of fission gas behavior in U 3 Si 2 that exploits a multiscale approach to fill the experimental data gap and can be progressively improved as new lower-length scale calculations and validation data become available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗

SPEARS: A Database-Invariant Spectral modeling API

The Spectral Physics Environment for Advanced Remote Sensing (SPEARS) application programming interface (API) is a Python-based, line-by-line, local thermal equilibrium (LTE) spectral modeling code which is optimized for simultaneously synthesizing optical spectra from any combination of fundamental spectroscopic databases. In this article, we contribute two novel spectral modeling techniques to the scientific literature. First we describe how SPEARS integrates a physics-based collisional model for calculating pressure broadening in the absence of available broadening coefficients. With this collisional model implementation, a generalized approach to fundamental spectroscopic databases can be achieved across multiple databases. We also detail our adaptive grid mesh algorithm developed to make the code scalable for simulating large spectral bandwidths at high spectral fidelity using intuitive grid parameters. Here, we present comparisons to other modeling tools, experiments, and provide a discussion on the SPEARS user interface.

47 OTHER INSTRUMENTATION↗