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,729 records · Page 96

The DaveMLTranslator: An Interface for DAVE-ML Aerodynamic Models

It can take weeks or months to incorporate a new aerodynamic model into a vehicle simulation and validate the performance of the model. The Dynamic Aerospace Vehicle Exchange Markup Language (DAVE-ML) has been proposed as a means to reduce the time required to accomplish this task by defining a standard format for typical components of a flight dynamic model. The purpose of this paper is to describe an object-oriented C++ implementation of a class that interfaces a vehicle subsystem model specified in DAVE-ML and a vehicle simulation. Using the DaveMLTranslator class, aerodynamic or other subsystem models can be automatically imported and verified at run-time, significantly reducing the elapsed time between receipt of a DAVE-ML model and its integration into a simulation environment. The translator performs variable initializations, data table lookups, and mathematical calculations for the aerodynamic build-up, and executes any embedded static check-cases for verification. The implementation is efficient, enabling real-time execution. Simple interface code for the model inputs and outputs is the only requirement to integrate the DaveMLTranslator as a vehicle aerodynamic model. The translator makes use of existing table-lookup utilities from the Langley Standard Real-Time Simulation in C++ (LaSRS++). The design and operation of the translator class is described and comparisons with existing, conventional, C++ aerodynamic models of the same vehicle are given.

Hill, Melissa A.↗

How to detect and correct myocardial creep in myocardial perfusion imaging using Rubidium-82 PET?

Reliability of myocardial blood flow (MBF) quantification in myocardial perfusion imaging (MPI) using PET can majorly be affected by the occurrence of myocardial creep when using pharmacologically induced stress. In this paper, we provide instructions on how to detect and correct for myocardial creep. For example, in each time frame of the PET images the myocardium contour and the observed activity have to be compared to check for misalignments. In addition, we provide an overview of the functionality of commonly used software packages to perform this quality control step as not all software packages currently provide this functionality. Furthermore, important clinical considerations to obtain accurate MBF measurements are given.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Fourth-order vibrational perturbation theory with the Watson Hamiltonian: Report of working equations and preliminary results

A derivation of fourth-order vibrational perturbation theory (VPT4) based on the Watson Hamiltonian in dimensionless rectilinear normal coordinates is presented. Terms that are linear and cubic in the (n k + 1/2), with n k being the zeroth-order harmonic oscillator quantum numbers, appear in fourth order and extend the much simpler second-order vibrational perturbation theory model. The rather involved expressions for the fourth-order terms are derived with Rayleigh-Schrödinger perturbation theory, the process of verifying their correctness is described, and a computer code to generate the VPT4 constants from the potential energy surface derivatives is provided. Here, the paper concludes with numerical examples featuring the H 2 O, Si 2 C, and cyclic-C 3 H 2 molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identifying candidate structured RNAs in CRISPR operons

Noncoding RNAs with secondary structures play important roles in CRISPR-Cas systems. Many of these structures likely remain undiscovered. We used a large-scale comparative genomics approach to predict 156 novel candidate structured RNAs from 36,111 CRISPR-Cas systems. A number of these were found to overlap with coding genes, including palindromic candidates that overlapped with a variety of Cas genes in type I and III systems. Among these 156 candidates, we identified 46 new models of CRISPR direct repeats and 1 tracrRNA. This tracrRNA model occasionally overlapped with predicted cas9 coding regions, emphasizing the importance of expanding our search windows for novel structure RNAs in coding regions. We also demonstrated that the antirepeat sequence in this tracrRNA model can be used to accurately assign thousands of predicted CRISPR arrays to type II-C systems. This study highlights the importance of unbiased identification of candidate structured RNAs across CRISPR-Cas systems.

59 BASIC BIOLOGICAL SCIENCES↗

3D full wave fast wave modeling with realistic HHFW antenna geometry and SOL plasma in NSTX-U

This paper reports the significant advancement of our ability to model and to understand how RF waves interact with the SOL plasma, by developing for the full torus 3D scrape-off layer (SOL) plasma simulation together with the antenna and core plasma in the NSTX-U device in HHFW frequency regime for a cold plasma model. The present paper extends the previous modeling (Bertelli et al 2020 AIP Conf. Proc. 2254 030001) to a full 3D device geometry including realistic antenna geometry in order to capture a more realistic experimental situation including the fully 3D wave field effects and the antenna plasma interaction in the SOL plasma, and, at the same time, resulting core wave propagation. Here, the central tool of the present work is the Petra-M code, which is a newly developed state-of-the-art generic electromagnetic simulation tool for modeling RF wave propagation based on MFEM (http://mfem.org), an open source scalable C++ finite element method library.

3D RF modeling↗

Twisted bilayer graphene. IV. Exact insulator ground states and phase diagram

Here, we derive the exact insulator ground states of the projected Hamiltonian of magic-angle twisted bilayer graphene (TBG) flat bands with Coulomb interactions in various limits, and study the perturbations away from these limits. We define the (first) chiral limit where the AA stacking hopping is zero, and a flat limit with exactly flat bands. In the chiral-flat limit, the TBG Hamiltonian has a U(4) × U(4) symmetry, and we find that the exact ground states at integer filling – 4 ≤ $\textit{ν}$ ≤ 4 relative to charge neutrality are Chern insulators of Chern numbers $ν_C$ = 4 – |$\textit{ν}$ |, 2 – |$\textit{ν}$ |, $\cdots$, |$\textit{ν}$ | – 4, all of which are degenerate. This confirms recent experiments where Chern insulators are found to be competitive low-energy states of TBG. When the chiral-flat limit is reduced to the nonchiral-flat limit which has a U(4) symmetry, we find $\textit{ν}$ = 0 , ± 2 has exact ground states of Chern number 0, while $\textit{ν}$ = ±1, ± 3 has perturbative ground states of Chern number $ν_C$ = ± 1, which are U(4) ferromagnetic. In the chiral-nonflat limit with a different U(4) symmetry, different Chern number states are degenerate up to second-order perturbations. In the realistic nonchiral-nonflat case, we find that the perturbative insulator states with Chern number $ν_C$ = 0 (0 < |$ν_C$| < 4 – |$\textit{ν}$ |) at integer fillings $\textit{ν}$ are fully (partially) intervalley coherent, while the insulator states with Chern number |$ν_C$| = 4 – |$\textit{ν}$ | are valley polarized. However, for 0 < |$ν_C$| ≤ 4 – |$\textit{ν}$ |, the fully intervalley coherent states are highly competitive (0.005 meV/electron higher). At nonzero magnetic field |$\textit{B}$ | > 0, a first-order phase transition for $\textit{ν}$ = ± 1, ± 2 from Chern number $ν_C$ = sgn ($\textit{νB}$)(2 – |$\textit{ν}$ |) to $ν_C$ = sgn ($\textit{νB}$)(4 – |$\textit{ν}$ |) is expected, which agrees with recent experimental observations. Lastly, the TBG Hamiltonian reduces into an extended Hubbard model in the stabilizer code limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Accelerating DCA++ (Dynamical Cluster Approximation) Scientific Application on the Summit supercomputer

Optimizing scientific applications on today’s accelerator-based high performance computing systems can be challenging, especially when multiple GPUs and CPUs with heterogeneous memories and persistent non-volatile memories are present. An example is Summit, an accelerator-based system at the Oak Ridge Leadership Computing Facility (OLCF) that is rated as the world’s fastest supercomputer to-date. New strategies are thus needed to expose the parallelism in legacy applications, while being amenable to efficient mapping to the underlying architecture.In this paper we discuss our experiences and strategies to port a scientific application, DCA++, to Summit. DCA++ is a highperformance research application that solves quantum manybody problems with a cutting edge quantum cluster algorithm, the dynamical cluster approximation.Our strategies aim to synergize the strengths of the different programming models in the code. These include: (a) streamlining the interactions between the CPU threads and the GPUs, (b) implementing computing kernels on the GPUs and decreasing CPU-GPU memory transfers, (c) allowing asynchronous GPU communications, and (d) increasing compute intensity by combining linear algebraic operations.Full-scale production runs using all 4600 Summit nodes attained a peak performance of 73.5 PFLOPS with a mixed precision implementation.We observed a perfect strong and weak scaling for the quantum Monte Carlo solver in DCA++, while encountering about 2× input/output (I/O) and MPI communication overhead on the time-to-solution for the full machine run. Our hardware agnostic optimizations are designed to alleviate the communication and I/O challenges observed, while improving the compute intensity and obtaining optimal performance on a complex, hybrid architecture like Summit.

Balduzzi, Giovanni↗

Utah FORGE: Composite 3D Seismic Velocity Model

This is a composite 3D seismic velocity that was constructed from compiled information from several local studies regarding seismic velocities and structural information. This seismic velocity model is provided in NonLinLoc format (slow_len), which is readily usable in NonLinLoc software. Other model formats and versions of the model can be produced using the Python script provided with this data set. Details on how the model was created and prior velocity and structural information was used is provided in the accompanying documentation.

15 GEOTHERMAL ENERGY↗

Comparative study of decay heat calculations with FLUKA and MCNP/CINDER2008

In designing a high-power tungsten target, decay heat driven temperature rise in the spallation volume is a safety concern during maintenance and in loss of coolant accidents. As tungsten hydrates and becomes volatile in steam at above 800 ◦C, it is important to keep the target temperature below this threshold when active cooling is unavailable. Decay heat in a target is calculated with particle transport simulation codes combined with transmutation codes. The calculated decay heat usually differs depending on the nuclear cross sections and the decay particle transport models built in the code architecture. In this paper, we calculated decay heat of a water-cooled tungsten target using popular particle transport codes, FLUKA and MCNP6® paired with CINDER2008 and compared the results. The target-moderator-reflector (TMR) system is modeled with a water-cooled solid tungsten target, water premoderators, liquid hydrogen cold moderators and beryllium reflectors. Water-cooled stainless-steel shielding is modeled around the TMR system. The tungsten volume is clad with a thin layer of erosion/corrosion resistant material. This study provides information about the uncertainty range in decay heat prediction of high-power spallation targets for hazard analysis.

Lee, Yong Joong↗

Report on FY 2020 creep, fatigue and creep fatigue testing of Alloy 709 base metal at ORNL

The testing activities and research in support of ASME Code qualification of Alloy 709, an advanced austenitic steel, are being carried out at Oak Ridge National Laboratory (ORNL), Argonne National Laboratory, and Idaho National Laboratory. This report summarizes the status and results of FY 2020 planned testing at ORNL. Uniaxial tensile tests on the electroslag remelt (ESR) solution-annealed plate (heat number 58776-3RBC) with additional heat treatment were performed for the baseline mechanical properties evaluation for ASME code qualification of Alloy 709. The tensile properties were found to meet ASME SA-213 specifications and were comparable to the Nippon Steel NF709 data generated for an ASME Section I Code Case for seamless tubing. ORNL was tasked to carry out a subset of the Code Case testing for creep rupture. Creep rupture data from 18 tests were generated on a solution-annealed ESR plate with a solution annealing temperature of 1100°C. There are total of 19 intermediate and long-term creep rupture tests of Alloy 709 ongoing at ORNL. A preliminary fatigue design curve at 760°C was developed for Alloy 709, and the results show that the fatigue design curve of Alloy 709 is comparable to that of Alloy 800H at 760°C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Case Study: Leveraging GenAI to Build AI-based Surrogates and Regressors for Modeling Radio Frequency Heating in Fusion Energy Science

This work presents a detailed case study on using Generative AI (GenAI) to develop AI surrogates for simulation models in fusion energy research. The scope includes the methodology, implementation, and results of using GenAI to assist in model development and optimization, comparing these results with previous manually developed models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CSAPR2 CMAC 2.0 level c1 data.

Raw data from ARM precipitation radars must be corrected for atmospheric phenomena and instrument characteristics (e.g., attenuation, clutter) to retrieve precipitation properties. The Corrected Moments in Antenna Coordinates Version 2 (CMAC2) value-added product (VAP) is a set of algorithms and code that does such corrections, and it also retrieves precipitation quantities from the radar measurements. Similar to the X-SAPR radars at the ARM main site, C-SAPR2 radar data also needs corrections and improvements. CMAC 2.0 has been updated to work with the ARM C-SAPR 2 radar. Data and fields that have been processed to: correct for velocity aliasing, unfold and generate a cross-polarimetric phase difference that is monotonically increasing, removing impulses caused by non-uniform beam filling and phase shift on backscatter, recalculate specific differential phase using a 20-point Sobel filter on the aforementioned phase, correct for liquid path attenuation using the polarimetric signals, estimate rainfall rates at the gate using the specific attenuation. In addition, CMAC writes the data out into a community-standard format netCDF File using the CF/Radial conventions. The data are therefore compatible with new and existing National Center for Atmospheric Research (NCAR) tools such as RadXConvert for converting to a variety of popular file formats.

54 ENVIRONMENTAL SCIENCES↗

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-2

Jet-A spray combustion has been evaluated in gas turbine combustion with the use of propane chemical kinetics as the first approximation for the chemical reactions. Here, the numerical solutions are obtained by using the KIVA-2 computer code. The KIVA-2 code is the most developed of the available multidimensional combustion computer programs for application of the in-cylinder combustion dynamics of internal combustion engines. The released version of KIVA-2 assumes that 12 chemical species are present; the code uses an Arrhenius kinetic-controlled combustion model governed by a four-step global chemical reaction and six equilibrium reactions. Researchers efforts involve the addition of Jet-A thermophysical properties and the implementation of detailed reaction mechanisms for propane oxidation. Three different detailed reaction mechanism models are considered. The first model consists of 131 reactions and 45 species. This is considered as the full mechanism which is developed through the study of chemical kinetics of propane combustion in an enclosed chamber. The full mechanism is evaluated by comparing calculated ignition delay times with available shock tube data. However, these detailed reactions occupy too much computer memory and CPU time for the computation. Therefore, it only serves as a benchmark case by which to evaluate other simplified models. Two possible simplified models were tested in the existing computer code KIVA-2 for the same conditions as used with the full mechanism. One model is obtained through a sensitivity analysis using LSENS, the general kinetics and sensitivity analysis program code of D. A. Bittker and K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. The other model is based on the work published by C. K. Westbrook and F. L. Dryer.

Nguyen, H. L.↗

Mach number validation of a new zonal CFD method (ZAP2D) for airfoil simulations

A closed-loop overlapped velocity coupling procedure has been utilized to combine a two-dimensional potential-flow panel code and a Navier-Stokes code. The fully coupled two-zone code (ZAP2D) has been used to compute the flow past a NACA 0012 airfoil at Mach numbers ranging from 0.3 to 0.84 near the two-dimensional airfoil C(lmax) point for a Reynolds number of 3 million. For these cases, the grid domain size can be reduced to 3 chord lengths with less than 3-percent loss in accuracy for freestream Mach numbers through 0.8. Earlier validation work with ZAP2D has demonstrated a reduction in the required Navier-Stokes computation time by a factor of 4 for subsonic Mach numbers. For this more challenging condition of high lift and Mach number, the saving in CPU time is reduced to a factor of 2.

Strash, Daniel J.↗

Fighter-like vortical flow experiments/computations

The objective of this research project is to determine the ability of Euler and Navier-Stokes codes to predict vortex/shock-dominated flow that is representative of modern fighter aircraft. The motivation for this project is fourfold: (1) current fighter aircraft are capable of operating beyond C(sub L(sub MAX)); (2) high angle-of-attack vortex/shock-dominated flows are not well understood; (3) current design methods are of the trial-and-error type; and (4) current data bases are inadequate for Computational Fluid Dynamics (CFD) validation. All data and results are presented in viewgraph format.

Erickson, Larry L.↗

The first ultraviolet spectropolarimetric study of NGC 1068

We present the first ultraviolet spectropolarimetric observations of the Type 2 Seyfert galaxy NGC 1068 obtained by the Wisconsin Ultraviolet Photo Polarimeter Experiment (WUPPE). In this communication we discuss the continuum polarization in the ultraviolet. Shortward of 2500 A the polarization becomes constant at 16 percent at a position angle of 112 deg. This is in general agreement with the 1983 prediction of Miller and Antonucci (1983) and their thick torus model, in which a hidden Type 1 Seyfert nucleus is seen in light scattered by electrons outside the obscured nuclei.

Code, Arthur D.↗

COSMIC monthly progress report

Activities of the Computer Software Management and Information Center (COSMIC) are summarized for the month of August, 1993. Tables showing the current inventory of programs available from COSMIC are presented and program processing and evaluation activities are discussed. Ten articles were prepared for publication in the NASA Tech Brief Journal. These articles (included in this report) describe the following software items: (1) MOM3D - A Method of Moments Code for Electromagnetic Scattering (UNIX Version); (2) EM-Animate - Computer Program for Displaying and Animating the Steady-State Time-Harmonic Electromagnetic Near Field and Surface-Current Solutions; (3) MOM3D - A Method of Moments Code for Electromagnetic Scattering (IBM PC Version); (4) M414 - MIL-STD-414 Variable Sampling Procedures Computer Program; (5) MEDOF - Minimum Euclidean Distance Optimal Filter; (6) CLIPS 6.0 - C Language Integrated Production System, Version 6.0 (Macintosh Version); (7) CLIPS 6.0 - C Language Integrated Production System, Version 6.0 (IBM PC Version); (8) CLIPS 6.0 - C Language Integrated Production System, Version 6.0 (UNIX Version); (9) CLIPS 6.0 - C Language Integrated Production System, Version 6.0 (DEC VAX VMS Version); and (10) TFSSRA - Thick Frequency Selective Surface with Rectangular Apertures. Activities in the areas of marketing, customer service, benefits identification, maintenance and support, and dissemination are also described along with a budget summary.

Source record↗

Characterization of Metal Matrix Composites

Experimental methods were developed, adapted, and applied to the characterization of a metal matrix composite system, namely, silicon carbide/aluminim (SCS-2/6061 Al), and its constituents. The silicon carbide fiber was characterized by determining its modulus, strength, and coefficient of thermal expansion. The aluminum matrix was characterized thermomechanically up to 399 C (750 F) at two strain rates. The unidirectional SiC/Al composite was characterized mechanically under longitudinal, transverse, and in-plane shear loading up to 399 C (750 F). Isothermal and non-isothermal creep behavior was also measured. The applicability of a proposed set of multifactor thermoviscoplastic nonlinear constitutive relations and a computer code was investigated. Agreement between predictions and experimental results was shown in a few cases. The elastoplastic thermomechanical behavior of the composite was also described by a number of new analytical models developed or adapted for the material system studied. These models include the rule of mixtures, composite cylinder model with various thermoelastoplastic analyses and a model based on average field theory. In most cases satisfactory agreement was demonstrated between analytical predictions and experimental results for the cases of stress-strain behavior and thermal deformation behavior at different temperatures. In addition, some models yielded detailed three-dimensional stress distributions in the constituents within the composite.

Daniel, I. M.↗