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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

The Far-Infrared Spectrum of Arp 220

ISO/LWS grating observations of the ultraluminous infrared galaxy Arp 220 shows absorption in molecular lines of OH, H 2 0 , CH, NH, and "3, well as in the [0 I] 63 pm line and emission in the [C 111 158 pm line. We have modeled the continuum and the emission/absorption of all observed features by means of a non-local radiative transfer code. The continuum from 25 to 1300 pm is modeled AS A WARM (106 K) NUCLEAR REGION THAT IS OPTICALLY THICK IN THE FAR-INFRARED, attenuated by an extended region (size 2") that is heated mainly through absorption of nuclear infrared radiation. The molecular absorption in the nuclear region is characterized by high excitation due to the high infrared radiation density. The OH column densities are high toward the nucleus and the extended region (about 2 x 10 sup 17 cm sup-2). The H2O column density is also high toward the nucleus (2 - 10 x 1017 cm-2) and lower in the extended region. The column densities in a halo that accounts for the absorption by the lowest lying levels are similar to what are found in the diffuse clouds toward the star forming regions in the Sgr B2 molecular cloud complex near the Galactic Center. Most notable are the high column densities found for NH and NH3 toward the nucleus, with values of about 1.5 x 10supl6 cmsup-2 and about 3 x 10supl6 cmsup-2, respectively, whereas the NH2 column density is lower than about 2 x 10sup15 cmsup-2. A combination of PDRs in the extended region and hot cores with enhanced H20 photodissociation and a possible shock contribution in the nuclei may explain the relative column densities of OH and H20, whereas the nitrogen chemistry may be strongly affected by cosmic ray ionization. The [C II] 158 pm line is well reproduced by our models and its "deficit" relative to the CII/FIR ratio in normal and starburst galaxies is suggested to be mainly a consequence of the dominant non-PDR component of far- infrared radiation, ALTHOUGH OUR MODELS ALONE CANNOT RULE OUT EXTINCTION EFFECTS IN THE NUCLEI.

Gonzalez-Alfonso, Eduardo↗
NASA NTRS1986

Two-phase flows and heat transfer within systems with ambient pressure above the thermodynamic critical pressure

In systems where the design inlet and outlet pressures P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.

Hendricks, R. C.↗
NASA NTRS1986

Two-phase flows and heat transfer within systems with ambient pressure above the thermodynamic critical pressure

In systems where the design inlet and outlet pressure P sub amb are maintained above the thermodynamic critical pressure P sub c, it is often assumed that heat and mass transfer are governed by single-phase relations and that two-phase flows cannot occur. This simple rule of thumb is adequate in many low-power designs but is inadequate for high-performance turbomachines, boilers, and other systems where two-phase regions can exist even though P sub amb P sub c. Heat and mass transfer and rotordynamic-fluid-mechanic restoring forces depend on momentum differences, and those for a two-phase zone can differ significantly from those for a single-phase zone. By using a laminar, variable-property bearing code and a rotating boiler code, pressure and temperature surfaces were determined that illustrate nesting of a two-phase region within a supercritical pressure region. The method of corresponding states is applied to bearings with reasonable rapport.

Hendricks, R. C.↗
NASA NTRS2005

The Far-Infrared Spectrum of Arp 220

ISO/LWS grating observations of the ultraluminous infrared galaxy Arp 220 shows absorption in molecular lines of OH, H(sub 2)O, CH, NH, and NH(sub 3), as well as in the [O I] 63 micron line and emission in the [C II] 158 micron line. We have modeled the continuum and the emission/absorption of all observed features by means of a non-local radiative transfer code. The continuum from 25 to 1300 microns is modeled as a warm (106 K) nuclear region that is optically thick in the far-infrared, attenuated by an extended region (size 2") that is heated mainly through absorption of nuclear infrared radiation. The molecular absorption in the nuclear region is characterized by high excitation due to the high infrared radiation density. The OH column densities are high toward the nucleus (2 - 6 x 10(exp 17) cm(exp -2)) and the extended region (approximately 2 x 10(exp 17) cm(exp -2)). The H(sub 2)O column density is also high toward the nucleus (2 - 10 x 10(exp 17) cm(exp -2)) and lower in the extended region. The column densities in a halo that accounts for the absorption by the lowest lying levels are similar to what are found in the diffuse clouds toward the star forming regions in the Sgr B2 molecular cloud complex near the Galactic Center. Most notable are the high column densities found for NH and NH(sub 3) toward the nucleus, with values of approximately 1.5 x 10(exp 16) cm(exp -2) and approximately 3 x 10(exp 16) cm(exp -2), respectively, whereas the NH(sub 2) column density is lower than approximately 2 x 10(exp 15) cm(exp -2). A combination of PDRs in the extended region and hot cores with enhanced H(sub 2)O photodissociation and a possible shock contribution in the nuclei may explain the relative column densities of OH and H(sub 2)O, whereas the nitrogen chemistry may be strongly affected by cosmic ray ionization. The [C II] 158 micron line is well reproduced by our models and its deficit relative to the CII/FIR ratio in normal and starburst galaxies is suggested to be mainly a consequence of the dominant non-PDR component of far-infrared radiation, although our models alone cannot rule out extinction effects in the nuclei.

Gonzalez-Alfonso, Eduardo↗
DOE OSTI2024

Towards reverse mode automatic differentiation of Kokkos-based codes

Derivative computation is a key component of optimization, sensitivity analysis, uncertainty quantification, and the solving of nonlinear problems. Automatic differentiation (AD) is a powerful technique for evaluating such derivatives, and in recent years, has been integrated into programming environments such as Jax, PyTorch, and TensorFlow to support derivative computations needed for training of machine learning models, facilitating wide-spread use of these technologies. The C++ language has become the de facto standard for scientific computing due to numerous factors, yet language complexity has made the wide-spread adoption of AD technologies for C++ difficult, hampering the incorporation of powerful differentiable programming approaches into C++ scientific simulations. This is exacerbated by the increasing emergence of architectures, such as GPUs, with limited memory capabilities and requiring massive thread-level concurrency. C++ AD tools must effectively use these environments to bring novel scientific simulations to next-generation DOE experimental and observational facilities. In this project, we investigated source transformation-based automatic differentiation using LLVM compiler infrastructure to automatically generate portable and efficient gradient computations of Kokkos-based code. We have demonstrated that our proposed strategy is feasible by investigating the usage of a prototype LLVM-based source transformation tool to generate gradients of simple functions made of sequences of simple Kokkos parallel regions. Speedups of up to 500x compared to Sacado were observed on NVIDIA V100 GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Comparative analysis of energy deposition modes available in Serpent 2 within the framework of the supercritical water reactor - Fuel qualification test reactor physics benchmark

A joint European Canadian Chinese development of a supercritical water-cooled small modular reactor (SCW-SMR) technology is in progress since September 2020 in the framework of a Horizon 2020 project called ECC-SMART. As a main purpose of the project, proper estimates of energy deposition and its spatial distribution are prerequisites for the accurate analysis of safety related parameters of the SCW-SMR concept under development. A supercritical water reactor fuel computational benchmark model, provided by Canadian Nuclear Laboratories, was applied for detailed comparison of different energy deposition calculation options available in the Serpent 2 Monte Carlo code. The effect of energy deposition options on the normalization of the results as well as on the spatial distribution of the energy deposition are discussed. Consistent energy deposition calculation methods are presented between three Monte Carlo codes, viz., Serpent 2, MCNP6 and OpenMC. Although resource-intensive, the use of the coupled neutron-photon transport mode of Serpent 2 is recommended for accurate spatial and quantitative characterization of energy deposition in the SCW-SMR fuel assemblies, accounting for both neutron and photon heating of all the materials. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗
DOE OSTI2022

Understanding and controlling capsule symmetry in near vacuum hohlraums at the National Ignition Facility

The near vacuum hohlraum platform is an inertial confinement fusion design at the National Ignition Facility (NIF) that uses the lowest practical density of helium gas of 30 μg/cc to fill the hohlraum, which is ten times lower than now used routinely. This has several advantages, such as high laser coupling; however, the inability to understand and simulate the symmetry of the imploded capsule has limited the use of this platform. This work presents the first simulations that are able to accurately capture the highly prolate implosion seen experimentally without unphysical, ad hoc model changes. While previous investigations attributed this asymmetry to multi-species interpenetration in the hohlraum, we find that this alone has little effect on symmetry. Instead, it is the presence of crossed-beam energy transfer (CBET), occurring with no applied wavelength shift between the laser beams, that increases the laser power to the inner cones and causes a more prolate implosion. The effect of CBET is increased in the simulation model when the hohlraum laser entrance hole hardware is included. Using this understanding, CBET is exploited by shifting the inner-beam wavelength by –0.75 Å (at 1ω) with respect to the outer-beams. This transfers laser power to the outer-beams in contrast to positive wavelength shifts as done routinely on NIF and produces a round capsule implosion in our simulations. Furthermore, this work shows the possibility of the near vacuum hohlraum as a viable experimental platform.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The accuracy of multi-group models for nonlocal electron transport in magnetized plasmas

In the extreme conditions of inertial confinement fusion experiments, heat flow plays a vital role, but local diffusive models frequently break down and overestimate the heat flow. The situation becomes more complicated again in the significant magnetic fields generated during laser–plasma interactions or in magnetized fusion schemes. Accurate non-local and magnetized heat flow computations can be carried out using Vlasov–Fokker–Planck (VFP) simulations, but these are computationally expensive. There is, therefore, significant interest in using faster multi-group models to accurately calculate the non-local heat flow in magnetized plasmas. We benchmark two such multi-group models for calculating the heat flow, M1 and hybrid-AWBS-BGK, against diffusive models and full VFP simulations, before applying the models to realistic example test cases, both magnetized and unmagnetized. We find that the multi-group models generally perform very well for moderate non-localities up to kλmfp∼0.01, but the computational cost increases dramatically. hybrid-AWBS-BGK performs more effectively than M1 at high non-localities, up to kλmfp∼1, due to its adaptive solver and robust P1 closure, but tends to fail in very strong magnetic fields. Both codes are much faster than VFP simulations but are still slow in steep temperature gradients.

Arran, C. (ORCID:0000000286448118)↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: A. identifying options for the control parameters*

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter-erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding. This part A of the present study reports new results for largely the same ITER cases, confirming the strong correlations reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne [Ne/s], and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In the companion report, part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction 2 point modelling, Rev2PM. The input to the reduced model is a value of the variable pair (T et @q ⟂,pk ,Inj Ne ) and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , n e,Xpt (the electron density at the poloidal location of the X-point), etc. In part B the reduced model is quantitatively characterized using one half of the code cases; it is then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half of the cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Aeroelastic code comparison using the IEA 22MW reference turbine

Reference wind turbine designs and the associated aeroelastic models are widely used in both research and industry. Reference models representing future concepts are of particular interest. Current state of the art aeroelastic tools are relied upon to design the next generation of large wind turbines. However, modelling assumptions may be invalidated by upcoming very large turbines, and different aeroelastic tools may give inconsistent results. A 22MW turbine model has been defined as part of International Energy Agency (IEA) Wind Task 55 on Reference Wind Turbines and Farms to represent future turbines to be deployed in the 2030s. In this study, an aeroelastic model of this turbine has been created in four tools; Bladed, HAWC2, OpenFAST, and QBlade. Code comparisons are presented for steady state operation, linear stability analysis, and time domain power production simulations in steady and turbulent wind. Generally, the codes show a good agreement, but with some differences present in the linear stability analysis, periodic azimuthal variation, and time domain simulations. The models are a good basis for further study with the IEA 22MW turbine, and further code comparison exercises.

17 WIND ENERGY↗
DOE OSTI2023

Fundamental cosmology from ANDES precision spectroscopy

ABSTRACT Fundamental cosmology observations, such as the detection of the redshift drift and tests of the universality of physical laws, are key science and design drivers of the ArmazoNes high Dispersion Echelle Spectrograph (ANDES), an Extremely Large Telescope instrument. While separate forecasts for each of them have been reported, we have developed Fisher matrix based forecast tools combining both of these observables. We demonstrate the synergies between the two ANDES data sets, quantifying the improvements in cosmology and fundamental physics parameter constraints for two separate theoretical paradigms. We publicly release this forecast code, which is one of the tools for the optimization of the ANDES observing strategy.

Marques, C. M. J.↗
DOE OSTI2024

An ascomycete H4 variant with an unknown function

Histone variants leading to altered nucleosome structure, dynamics and DNA accessibility occur frequently, albeit rarely for H4. We carried out a comprehensive in silico scrutiny of fungal genomes, which revealed the presence of a novel H4 variant (H4E) in the ascomycetes, throughout the Pezizomycotina, in basal species of the Taphrinomycotina and also in the Glomeromycota. The coding cognate genes show a specific intron/exon organization, different from H4 canonical genes. H4Es diverge from canonical H4s mainly in the N- and C-terminal extensions, showing marked differences in the distribution and number of Lys and Arg residues, which may result in novel post-translational modifications. In Aspergillus nidulans (Pezizomycotina, Eurotiomycetes) the H4E variant protein level is low in mycelia. However, the encoding gene is well expressed at 37°C under nitrogen starvation. H4E localizes to the nucleus and interacts with H3, but its absence or overexpression does not result in any detectable phenotype. Deletion of only one of the of the two canonical H4 genes results in a strikingly impaired growth phenotype, which indicates that H4E cannot replace this canonical histone. Thus, an H4 variant is present throughout a whole subphylum of the ascomycetes, but with hitherto no experimentally detectable function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The landscape of Chlamydomonas histone H3 lysine 4 methylation reveals both constant features and dynamic changes during the diurnal cycle

Chromatin modifications are epigenetic regulatory features with major roles in various cellular events, yet they remain understudied in algae. Here, we interrogated the genome-wide distribution pattern of mono- and trimethylated histone H3 lysine 4 (H3K4) using chromatin-immunoprecipitation followed by deep-sequencing (ChIP-seq) during key phases of the Chlamydomonas cell cycle: early G 1 phase, Zeitgeber Time 1 (ZT1), when cells initiate biomass accumulation, S/M phase (ZT13) when cells are replicating DNA and undergoing mitosis, and late G 0 phase (ZT23) when they are quiescent. Tri-methylated H3K4 was predominantly enriched at transcription start sites of the majority of protein coding genes (85%). The likelihood of a gene being marked by H3K4me3 correlated with it being transcribed at some point during the life cycle but not necessarily by continuous active transcription, as exemplified by early zygotic genes, which may remain transcriptionally dormant for thousands of generations between sexual cycles. The exceptions to this rule were around 120 loci, some of which encode non-poly-adenylated transcripts, such as small nuclear RNAs and replication-dependent histones that had H3K4me3 peaks only when they were being transcribed. Mono-methylated H3K4 was the default state for the vast majority of histones that were bound outside of transcription start sites and terminator regions of genes. A small fraction of the genome that was depleted of any H3 lysine 4 methylation was enriched for DNA cytosine methylation and the genes within these DNA methylation islands were poorly expressed. Besides marking protein coding genes, H3K4me3 ChIP-seq data served also as a annotation tool for validation of hundreds of long non-coding RNA genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

The Right Time at the Right Places

In my undergraduate studies at the Massachusetts Institute of Technology (MIT), the short-lived chemical physics major allowed me to evade a number of courses required for chemistry majors. Thus, it was possible to take many physics courses and most of the advanced PhD-level courses in physical chemistry. I also took the introductory electrical engineering course in computer programming. The latter allowed me to write lots of computer code as a part of my (passing, but largely unsuccessful) senior thesis directed graciously by Professor Walter Thorson. As recommended by MIT Professor John C. Slater, I moved to Stanford University with Professor Frank Harris as my PhD supervisor. Frank was the perfect advisor for me, providing very close direction during the first year, and then allowing me to develop more and more independently. Within a few days of my twenty-fifth birthday, I became an assistant professor of chemistry at the University of California, Berkeley. Eighteen years later, I moved to the University of Georgia as director of a new research institute.

autobiography↗
DOE OSTI2019

Inelastic analysis procedure based on the Grade 91 unified viscoplastic constitutive model for ASME implementation

This report describes a design by inelastic analysis procedure for Grade 91 to be proposed for incorporation in the ASME Boiler and Pressure Vessel Code. The package of proposed additions includes guidance on inelastic modeling and a reference constitutive model for Grade 91 to be used with the existing Section III, Division 5, Subsection HB, Subpart B rules for design by inelastic analysis. The current Code specifies acceptance criteria for design by inelastic analysis but does not provide an acceptable constitutive model for Grade 91 or any of the other Class A materials. The proposed Nonmandatory Appendix would provide such a reference model, greatly easing the use of design by inelastic analysis. As the inelastic methods produce significantly less overconservative designs than the design by elastic analysis rules this proposed Code change could reduce the design and fabrication cost of future high temperature reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Development of a multiaxial deformation measure and creep-fatigue damage summation for multiple load cycle types in support of an improved creep-fatigue design method

This report describes methods for extending the elastic perfectly plastic (EPP) simplified model test (SMT) approach to cover multiaxial load conditions and methods to combine EPP+SMT use fractions for multiple loading cycles. The completion of this work, in conjunction with companion reports on constructing and validating SMT design curves, marks the development of a complete methodology for assessing creep-fatigue damage in high temperature nuclear reactor structural components. This EPP+SMT method aims to better account for elastic follow up effects, simplify the design process, and reduce overconservatism when compared to current simplified bounding approaches in the ASME Boiler and Pressure Vessel Code Section III, Division 5, Subsection HB, Subpart B covering high temperature nuclear reactors.

36 MATERIALS SCIENCE↗