Search NASA⌕ Search

SEARCH · Search NASA

Results for “code comparison study”

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 505 records · Page 28

A kinetic study of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas

We present particle-in-cell simulations with Monte Carlo collisions of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas. We study the energy balance in the one-dimensional expansion of a hot-spot by simulating Coulomb collisions, fusion reactions, and bremsstrahlung emission with a Monte Carlo model and inverse bremsstrahlung absorption using a new PIC model. This allows us to self-consistently capture the alpha particle heating and radiative losses in the expanding hot-spot and surrounding cold fuel. After verifying our model in a code-to-code comparison with both kinetic and fluid codes for the case of a deuterium–tritium hot-spot, we simulate the expansion of a proton–boron hot-spot initialized at 200 keV and 1,000 g/cm 3 . Our model predicts that energy radiated by the hot-spot is recaptured by the surrounding high-density opaque fuel reducing the expansion work done by the propagating burn wave. As a result, we find the net fusion energy produced over the course of $20$~ps is twice the initial hot-spot energy independent of whether radiation physics is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aerodynamic analysis of hypersonic waverider aircraft

The purpose of this study is to validate two existing codes used by the Systems Analysis Branch at NASA ARC, and to modify the codes so they can be used to generate and analyze waverider aircraft at on-design and off-design conditions. To generate waverider configurations and perform the on-design analysis, the appropriately named Waverider code is used. The Waverider code is based on the Taylor-Maccoll equations. Validation is accomplished via a comparison with previously published results. The Waverider code is modified to incorporate a fairing to close off the base area of the waverider configuration. This creates a more realistic waverider. The Hypersonic Aircraft Vehicle Optimization Code (HAVOC) is used to perform the off-design analysis of waverider configurations generated by the Waverider code. Various approximate analysis methods are used by HAVOC to predict the aerodynamic characteristics, which are validated via a comparison with experimental results from a hypersonic test model.

Sandlin, Doral R.↗

Reactor Performance and Safety Characteristics of Beryllium-Based Composite Moderators as Replacements for Graphite in mHTGRs

Here, this study evaluates beryllium-based two-phase composite moderators as an alternative to graphite in an evaluation of reactor performance and safety characteristics. Historically, modular high-temperature gas-cooled reactors (mHTGRs) use graphite as a moderator because of its high moderating ratio and reasonable thermal properties; however, graphite has unfavorable properties under irradiation, which can require component replacement and a significant radioactive waste burden. In this assessment, we explore advanced moderators comprised of magnesium oxide (MgO) as the host matrix and beryllium metal and/or beryllium oxide (Be and/or BeO) as the entrained moderating phase. For the reactor performance and thermal-hydraulic safety analysis, the core design model of the General Atomics mHTGR-350 was used to demonstrate the feasibility of a “drop-in” replacement of graphite using the beryllium-based moderators. We employed the neutronics code Serpent to analyze the moderating behavior of the composite moderators with comparisons drawn to graphite. We performed a scoping analysis of accidents for mHTGRs using RELAP to show that these moderators do not present impediments to safety and are expected to stay within temperature limits. Measured thermophysical properties of the composite moderators are used in the thermal-hydraulic assessments. Our analysis reveals that the two-phase composite MgO-matrix beryllium-based moderators are a suitable replacement for graphite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An empirical comparison of software fault tolerance and fault elimination

Reliability is an important concern in the development of software for modern systems. Some researchers have hypothesized that particular fault-handling approaches or techniques are so effective that other approaches or techniques are superfluous. The authors have performed a study that compares two major approaches to the improvement of software, software fault elimination and software fault tolerance, by examination of the fault detection obtained by five techniques: run-time assertions, multi-version voting, functional testing augmented by structural testing, code reading by stepwise abstraction, and static data-flow analysis. This study has focused on characterizing the sets of faults detected by the techniques and on characterizing the relationships between these sets of faults. The results of the study show that none of the techniques studied is necessarily redundant to any combination of the others. Further results reveal strengths and weakness in the fault detection by the techniques studied and suggest directions for future research.

Shimeall, Timothy J.↗

Comparison between computations and experimental data in unsteady three-dimensional transonic aerodynamics, including aeroelastic applications

Comparisons were made of computed and experimental data in three-dimensional unsteady transonic aerodynamics, including aeroelastic applications. The computer code LTRAN3, which is based on small-disturbance aerodynamic theory, was used to obtain the aerodynamic data. A procedure based on the U-g method was developed to compute flutter boundaries by using the unsteady aerodynamic coefficients obtained from LTRAN3. The experimental data were obtained from available NASA publications. All the studies were conducted for thin, unswept, rectangular wings with circular-arc cross sections. Numerical and experimental steady and unsteady aerodynamic data were compared for a wing with an aspect ratio of 3 and a thickness ratio of 5% at Mach numbers of 0.7 and 0.9. Flutter data were compared for a wing with an aspect ratio of 5. Two thickness ratios, 6% at Mach numbers of 0.715, 0.851, and 0.913, and 4% at Mach number of 0.904, were considered. Based on the unsteady aerodynamic data obtained from LTRAN3, flutter boundaries were computed; they were compared with those obtained from experiments and the code NASTRAN, which uses linear aerodynamics.

Guruswamy, P.↗

Experimental and computational studies of unsteady viscous flows

This paper describes progress in validating incompressible Navier-Stokes codes using image-processed flow-field data. Solutions of the Navier-Stokes equations using a predictor/corrector scheme are compared directly with experimental data. A uniform asymptotic expansion is used as an initial condition in order to correctly capture the initial growth of the vortex sheet. With the aid of an experimental study in a tow tank using the particle streak visualization method, measured and computed velocity vectors are compared directly. Comparison of the velocity data with computations shows little effect of initial conditions and varying grid geometry on the predictive capabilities.

Davis, S.↗

Evaluation of a CFD Method for Aerodynamic Database Development using the Hyper-X Stack Configuration

A computational fluid dynamic (CFD) study is performed on the Hyper-X (X-43A) Launch Vehicle stack configuration in support of the aerodynamic database generation in the transonic to hypersonic flow regime. The main aim of the study is the evaluation of a CFD method that can be used to support aerodynamic database development for similar future configurations. The CFD method uses the NASA Langley Research Center developed TetrUSS software, which is based on tetrahedral, unstructured grids. The Navier-Stokes computational method is first evaluated against a set of wind tunnel test data to gain confidence in the code s application to hypersonic Mach number flows. The evaluation includes comparison of the longitudinal stability derivatives on the complete stack configuration (which includes the X-43A/Hyper-X Research Vehicle, the launch vehicle and an adapter connecting the two), detailed surface pressure distributions at selected locations on the stack body and component (rudder, elevons) forces and moments. The CFD method is further used to predict the stack aerodynamic performance at flow conditions where no experimental data is available as well as for component loads for mechanical design and aero-elastic analyses. An excellent match between the computed and the test data over a range of flow conditions provides a computational tool that may be used for future similar hypersonic configurations with confidence.

Parikh, Paresh↗

High-Fidelity Simulations of a Quadrotor Vehicle for Urban Air Mobility

High-fidelity computational fluid dynamics simulations have been carried out to analyze NASA’s quadrotor air taxi concept for urban air mobility. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. The flow solver has been loosely coupled with the rotorcraft comprehensive analysis code. The vehicle simulated is a six-passenger quadrotor for air taxi operations. The quadrotor has been simulated in hover and cruise, with flight conditions representative of an air taxi mission. Initially, the quadcopter vehicle is analyzed using low-fidelity tools. Then, a complete analysis using high-fidelity simulations of the quadcopter is performed in conjunction with collective pitch control. Airloads, wake geometry, and performance are studied and compared. A comparison between isolated rotors and complete vehicle simulations is presented. NASA’s quadrotor is one of NASA’s concept vehicles designed in support of aircraft development for vertical take-off and landing air taxi operations.

ARMD↗

Quantitative Studies of the Optical and UV Spectra of Galactic Early B Supergiants

We undertake an optical and ultraviolet spectroscopic analysis of a sample of 20 Galactic B0-B5 supergiants of luminosity classes Ia, Ib, Iab, and II. Fundamental stellar parameters are obtained from optical diagnostics and a critical comparison of the model predictions to observed UV spectral features is made. Methods. Fundamental parameters (e.g., T(sub eff), log L(sub *), mass-loss rates and CNO abundances) are derived for individual stars using CMFGEN, a nLTE, line-blanketed model atmosphere code. The impact of these newly derived parameters on the Galactic B supergiant Ten scale, mass discrepancy, and wind-momentum luminosity relation is examined. Results. The B supergiant temperature scale derived here shows a reduction of about 1000-3000 K compared to previous results using unblanketed codes. Mass-loss rate estimates are in good agreement with predicted theoretical values, and all of the 20 BO-B5 supergiants analysed show evidence of CNO processing. A mass discrepancy still exists between spectroscopic and evolutionary masses, with the largest discrepancy occuring at log (L/(solar)L approx. 5.4. The observed WLR values calculated for B0-B0.7 supergiants are higher than predicted values, whereas the reverse is true for B1-B5 supergiants. This means that the discrepancy between observed and theoretical values cannot be resolved by adopting clumped (i.e., lower) mass-loss rates as for O stars. The most surprising result is that, although CMFGEN succeeds in reproducing the optical stellar spectrum accurately, it fails to precisely reproduce key UV diagnostics, such as the N v and C IV P Cygni profiles. This problem arises because the models are not ionised enough and fail to reproduce the full extent of the observed absorption trough of the P Cygni profiles. Conclusions. Newly-derived fundamental parameters for early B supergiants are in good agreement with similar work in the field. The most significant discovery, however, is the failure of CMFGEN to predict the correct ionisation fraction for some ions. Such findings add further support to revising the current standard model of massive star winds, as our understanding of these winds is incomplete without a precise knowledge of the ionisation structure and distribution of clumping in the wind. Key words. techniques: spectroscopic - stars: mass-loss - stars: supergiants - stars: abundances - stars: atmospheres - stars: fundamental parameters

Searle, S. C.↗

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data for: Patel et al. Carbon flux estimates are sensitive to data source: A comparison of field and lab temperature sensitivity data

This dataset contains data and code used for the paper "Carbon flux estimates are sensitive to data source: A comparison of field and lab temperature sensitivity data" [DOI COMING SOON]A large literature exists on mechanisms driving soil production of the greenhouse gases CO2 and CH4. Measurements of these gases’ fluxes are often performed using closed-chamber incubations in the laboratory or in situ, i.e., in the field. Although it is common knowledge that measurements obtained through field studies vs. laboratory incubations can diverge because of the vastly different conditions of these environments, few studies have systematically examined these patterns. It is crucial to understand the magnitude and reasons for any differences, as these data are used to parametrize and benchmark ecosystem- to global-scale models, which are then susceptible to the biases of the source data. Here, we specifically examine how greenhouse gas measurements may be influenced by whether the measurement/incubation was conducted in the field vs. laboratory, focusing on CO2 and CH4 measurements. We use Q10 of greenhouse gas flux (temperature sensitivity) for our analyses, because of the ubiquity of this metric in biological and Earth system sciences and its importance to many modeling frameworks. We predicted that laboratory measurements would be less variable, but also less representative of true field conditions. However, there was greater variability in the Q10 values calculated from lab-based measurements of CO2 fluxes, because lab experiments explore extremes rarely seen in situ, and reflect the physical and chemical disturbances occurring during sampling, transport, and incubation. Overall, respiration Q10 values were significantly greater in laboratory incubations (mean = 4.19) than field measurements (mean = 3.05), with strong influences of incubation temperature and climate region/biome. However, this was in part because field measurements typically represent total respiration (Rs), whereas lab incubations typically represent heterotrophic respiration (Rh), making direct comparisons difficult to interpret. Focusing only on Rh-derived Q10, these values showed almost identical distributions across laboratory (n = 1110) and field (n = 581) experiments, providing strong support for using the former as an experimental proxy for the latter, although we caution that geographic biases in the extant data make this conclusion tentative. Due to a smaller sample size of CH4 Q10 data, we were unable to perform a comparable robust analysis, but we expect similar interactions with soil temperature, moisture, and environmental/climatic variables. Our results here suggest the need for more concerted efforts to document and standardize these data, including sample and site metadata. This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/kaizadp/field_lab_q10 (DOI: 10.5281/zenodo.7106554).

54 ENVIRONMENTAL SCIENCES↗

Reduced cost alternatives to premise wiring using ATM and microcellular technologies

The cost of premises wiring keeps increasing due to personnel moves, new equipment, capacity upgrades etc. It would be desirable to have a wireless interface from the workstations to the fixed network, so as to minimize the wiring changes needed. New technologies such as microcellular personal communication systems are promising to bring down the cost of wireless communication. Another promising technology is Code Division Multiple Access (CDMA), which could dramatically increase the bandwidth available for wireless connections. In addition, Asynchronous Transfer Mode (ATM) technology is emerging as a technique for integrated management of voice, data, and video traffic on a single network. The focus of this investigation will be to assess the future utility of these new technologies for reducing the premise wiring cost at KSC. One of the issues to be studied is the cost comparison of 'old' versus 'new,' especially as time and technology progress. An additional issue for closer study is a feasible time-line for progress in technological capability.

Gejji, Raghvendra R.↗

A Study of Flow Separation in Transonic Flow Using Inviscid and Viscous Computational Fluid Dynamics (CFD) Schemes

A comparison of flow separation in transonic flows is made using various computational schemes which solve the Euler and the Navier-Stokes equations of fluid mechanics. The flows examined are computed using several simple two-dimensional configurations including a backward facing step and a bump in a channel. Comparison of the results obtained using shock fitting and flux vector splitting methods are presented and the results obtained using the Euler codes are compared to results on the same configurations using a code which solves the Navier-Stokes equations.

Rhodes, J. A.↗

Nonnuclear Experimental Capabilities to Support Design, Development, and Demonstration of Microreactors

This work provides a summary of selected experimental capabilities being developed to support nonnuclear testing and demonstration of technology in support of microreactors under the U.S. Department of Energy’s (DOE’s) Microreactor Program. Major capabilities include the Single Primary Heat Extraction and Removal Emulator (SPHERE) and the Microreactor Agile Non-nuclear Experimental Test Bed (MAGNET). The SPHERE facility allows for controlled testing of the steady-state and transient heat rejection capabilities of a single heat pipe using electrical heaters that simulate nuclear heating. The facility is capable of monitoring axial temperature profiles along the heat pipe and surrounding test articles during startup, steady-state operation, and transients. Instrumentation includes noncontact infrared thermal imaging, surface thermocouples, spatially distributed fiber optic temperature and strain sensors, electrical power meters, and a water-cooled, gas-gap calorimeter for quantifying heat rejection from the heat pipe. The facility can be operated under both vacuum and inert-gas conditions. The MAGNET facility is a large-scale, 250-kW electrically heated microreactor test bed to enable nonnuclear experimental evaluation of a variety of microreactor concepts. It can be supplied to electrically heat a scaled section of a microreactor and further test the capabilities of heat rejection systems. The initial MAGNET experiments will support technology maturation and reduce uncertainty and risk associated with the design, operation, and deployment of monolithic heat pipe–based reactors. However, this test bed can broadly be applied to multiple microreactor concepts to evaluate a wide range of thermal-hydraulic and structural phenomena such as interface coupling with power conversion units and other collocated systems. MAGNET can evaluate integral thermomechanical effects during electrical heating of an array of heat pipes in a larger test article. Examples of initial testing will include thermal stresses in the monolith and the impact of debonding of a heat pipe from the core block and how that failure could impact surrounding heat pipes, i.e., understanding the potential for cascading failure. This work also discusses some modeling capabilities that can support experiment design, analysis, and interpretation, including the heat pipe code Sockeye and a comparison of thermal-structural simulations performed using ABAQUS and STAR-CCM+.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SERPENT 2 and TRIPOLI-4 neutronics benchmarking of the Jules Horowitz reactor

Jules Horowitz Reactor (JHR) is a future center for nuclear material testing in Europe, currently under construction at CEA Cadarache in southern France. It has been modelled extensively with CEA's HORUS3D/N neutron calculation package, that includes a Monte-Carlo neutron transport code TRIPOLI-4. In order to support the future computational analyses on the reactor core physics, VTT's Monte Carlo neutronics code SERPENT 2 has been used to create an alternative model of the JHR, taking advantage of the pre-existing Constructive Solid Geometry (CSG) structures on which the TRIPOLI-4 model was built. This work describes a benchmark, in which two JHR core configurations were analyzed using both of these models. The first configuration is a fresh fuel core at the beginning of life, in which aluminum based dummy-devices are installed in the reactor test positions. The second configuration represents the end-of-cycle equilibrium core with a heterogeneous burnup distribution, in which CALIPSO devices containing cladding material samples are installed at the in-core experimental positions, while the reflector positions are occupied by DEN-REP devices containing UO{sub 2} fuel. Neutron flux spectra and heat deposition are the parameters of interest in this study, due to their high importance in creating the desired experimental irradiation environment. The impact of delayed gamma energy deposition will be studied in detail, since the two models account for it with different methods. This work aims to establish a SERPENT 2 model of JHR, understanding its differences and agreements with TRIPOLI-4, and discussing its potential future applications. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Lebedev Scheme for Ultrasound Simulation in Composites

The growing use of composite materials for aerospace applications has resulted in a need for quantitative nondestructive evaluation (NDE) methods appropriate for characterizing damage in composite components. NDE simulation tools, such as ultrasound models, can aid in enabling optimized inspection methods and establishing confidence in inspection capabilities. In this paper, a mathematical approach using the Lebedev Finite Difference (LFD) method is presented for ultrasonic wave simulation in composites. Boundary condition equations for implementing stress-free boundaries (necessary for simulation of NDE scenarios) are also presented. Quantitative comparisons between LFD guided wave ultrasound simulation results, experimental guided wave data, and dispersion curves are described. Additionally, stability tests are performed to establish the LFD code behavior in the presence of stress-free boundaries and low-symmetry anisotropy. Results show that LFD is an appropriate approach for simulating ultrasound in anisotropic composite materials and that the method is stable in the presence of low-symmetry anisotropy and stress-free boundaries. Studies presented in this paper include guided wave simulation in hexagonal, monoclinic, triclinic and layered composite laminates.

Guided wave↗

Comparison of Geant4-DNA and RITRACKS/RITCARD: Microdosimetry, Nanodosimetry and DNA Damage Predictions for Low to High LET Ions

Low linear energy transfer (LET) ionizing radiations induce homogeneously distributed DNA damage at the cellular scale, while high LET ionizing radiations are well known to create complex biological damage including clustered double strand breaks (DSB). Clustered DSB are critical events that are difficult to repair, leading to increased biological outcomes compared to low LET radiation exposure. Recent studies have shown that DNA compaction, which can be actively modified by the cell to allow DNA repair, influences the yield of DSB [1]. The goal of this work is to determine the impact of DNA compaction on chromosome aberration formation, a biomarker of radiation-induced cancer risk. The RITRACKS track structure code [2], together with the RITCARD code [3], allow one to simulate the transport of high LET ions, calculate DNA damage and repair, and chromosome aberration formation, but does not include detailed geometrical models of DNA. The Geant4-DNA toolkit provides such detailed models with different levels of DNA compaction [1,4], but does not model chromosome aberration formation. In this work, these highly detailed DNA geometries are used to compute DNA breaks with Geant4-DNA and combine the results with the RITCARD tool to compute chromosome aberrations. As a first step, we compared the transport models used in RITRACKS and Geant4-DNA and newly released heavy ion model transport (G4DNARuddIonisationExtendedModel Geant4 v11.2) to ensure consistency between the physical stage of the two radiation transport codes. We compared energy deposition at the micrometric and nanometric scales for different ions (H 250 MeV, H 150 Me, He 250 MeV/n, C 290 MeV/n, O 350 MeV/n, O 55 MeV/n, Si 170MeV/n, Ti 300 MeV/n, Ti 300 MeV/n, Fe 600 MeV/n, Fe 450 MeV/n and Fe 300 MeV/n) with LET ranging from ~0.4 keV/μm up to 235 keV/μm. Excellent agreement is found for both microdosimetric and nanodosimetric spectra for all ion types and energies. The calculation of DNA breaks and their complexity with the Geant4-DNA is an ongoing work, and comparison with RITCARD will be presented based on result availability at the time of the conference. [1] Tang et al. (2019). Med. Phys., 46(3), 1501-1511. [2] Plante et al. (2008). NJP, 10(12), 125020. [3] Plante et al. (2019). Rad. Res., 192(3), 282-298. [4] Incerti et al. (2018). Med. Phys., 45(8), e722-e739.

Floriane Poignant↗

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↗