Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “M CODES”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Validation of NASA Thermal Ice Protection Computer Codes: The Validation of Antice - Part 3

An experimental program was generated by the Icing Technology Branch at NASA Glenn Research Center to validate two ice protection simulation codes: (1) LEWICE/Thermal for transient electrothermal de-icing and anti-icing simulations, and (2) ANTICE for steady state hot gas and electrothermal anti-icing simulations. An electrothermal ice protection system was designed and constructed integral to a 36 inch chord NACA0012 airfoil. The model was fully instrumented with thermo-couples, RTD'S, and heat flux gages. Tests were conducted at several icing environmental conditions during a two week period at the NASA Glenn Icing Research Tunnel. Experimental results of running-wet and evaporative cases were compared to the ANTICE computer code predictions and are presented in this paper.

Al-Khalil, Kamel M.↗
NASA NTRS1994

A Portable Debugger for Parallel and Distributed Programs

In this paper, we describe the design and implementation of a portable debugger for parallel and distributed programs. The design incorporates a client-server model in order to isolate non-portable debugger code from the user interface. The precise definition of a protocol for client-server interaction permits a high degree of portability of the client user interface. Replication of server components permits the implementation of a debugger for distributed computations. Portability across message passing implementations is achieved with a protocol that dictates the interaction between a message passing library and the debugger. This permits the same debugger to be used both on PVM and MTI programs. The process abstractions used for debugging message-passing programs can be easily adapted to debug HPF programs at the source level. This allows the debugger to present information hidden in tool-generated code in a meaningful manner.

Cheng, Doreen Y.↗
NASA NTRS2004

Transonic Drag Prediction on a DLR-F6 Transport Configuration Using Unstructured Grid Solvers

A second international AIAA Drag Prediction Workshop (DPW-II) was organized and held in Orlando Florida on June 21-22, 2003. The primary purpose was to inves- tigate the code-to-code uncertainty. address the sensitivity of the drag prediction to grid size and quantify the uncertainty in predicting nacelle/pylon drag increments at a transonic cruise condition. This paper presents an in-depth analysis of the DPW-II computational results from three state-of-the-art unstructured grid Navier-Stokes flow solvers exercised on similar families of tetrahedral grids. The flow solvers are USM3D - a tetrahedral cell-centered upwind solver. FUN3D - a tetrahedral node-centered upwind solver, and NSU3D - a general element node-centered central-differenced solver. For the wingbody, the total drag predicted for a constant-lift transonic cruise condition showed a decrease in code-to-code variation with grid refinement as expected. For the same flight condition, the wing/body/nacelle/pylon total drag and the nacelle/pylon drag increment predicted showed an increase in code-to-code variation with grid refinement. Although the range in total drag for the wingbody fine grids was only 5 counts, a code-to-code comparison of surface pressures and surface restricted streamlines indicated that the three solvers were not all converging to the same flow solutions- different shock locations and separation patterns were evident. Similarly, the wing/body/nacelle/pylon solutions did not appear to be converging to the same flow solutions. Overall, grid refinement did not consistently improve the correlation with experimental data for either the wingbody or the wing/body/nacelle pylon configuration. Although the absolute values of total drag predicted by two of the solvers for the medium and fine grids did not compare well with the experiment, the incremental drag predictions were within plus or minus 3 counts of the experimental data. The correlation with experimental incremental drag was not significantly changed by specifying transition. Although the sources of code-to-code variation in force and moment predictions for the three unstructured grid codes have not yet been identified, the current study reinforces the necessity of applying multiple codes to the same application to assess uncertainty.

Lee-Rausch, E. M.↗
NASA NTRS1998

The stability of the RNA bases: implications for the origin of life

High-temperature origin-of-life theories require that the components of the first genetic material are stable. We therefore have measured the half-lives for the decomposition of the nucleobases. They have been found to be short on the geologic time scale. At 100 degreesC, the growth temperatures of the hyperthermophiles, the half-lives are too short to allow for the adequate accumulation of these compounds (t1/2 for A and G approximately 1 yr; U = 12 yr; C = 19 days). Therefore, unless the origin of life took place extremely rapidly (<100 yr), we conclude that a high-temperature origin of life may be possible, but it cannot involve adenine, uracil, guanine, or cytosine. The rates of hydrolysis at 100 degreesC also suggest that an ocean-boiling asteroid impact would reset the prebiotic clock, requiring prebiotic synthetic processes to begin again. At 0 degreesC, A, U, G, and T appear to be sufficiently stable (t1/2 >/= 10(6) yr) to be involved in a low-temperature origin of life. However, the lack of stability of cytosine at 0 degreesC (t1/2 = 17, 000 yr) raises the possibility that the GC base pair may not have been used in the first genetic material unless life arose quickly (<10(6) yr) after a sterilization event. A two-letter code or an alternative base pair may have been used instead.

Non-NASA Center↗
NASA NTRS2005

Approximate Micromechanics Treatise of Composite Impact

A formalism is described for micromechanic impact of composites. The formalism consists of numerous equations which describe all aspects of impact from impactor and composite conditions to impact contact, damage progression, and penetration or containment. The formalism is based on through-the-thickness displacement increments simulation which makes it convenient to track local damage in terms of microfailure modes and their respective characteristics. A flow chart is provided to cast the formalism (numerous equations) into a computer code for embedment in composite mechanic codes and/or finite element composite structural analysis.

Chamis, Christos C.↗
NASA NTRS2009

Navier-Stokes Simulation of a Heavy Lift Slowed-Rotor Compound Helicopter Configuration

Time accurate numerical simulations were performed using the Reynolds-averaged Navier-Stokes (RANS) flow solver OVERFLOW for a heavy lift, slowed-rotor, compound helicopter configuration, tested at the NASA Langley 14- by 22-Foot Subsonic Tunnel. The primary purpose of these simulations is to provide support for the development of a large field of view Particle Imaging Velocimetry (PIV) flow measurement technique supported by the Subsonic Rotary Wing (SRW) project under the NASA Fundamental Aeronautics program. These simulations provide a better understanding of the rotor and body wake flows and helped to define PIV measurement locations as well as requirements for validation of flow solver codes. The large field PIV system can measure the three-dimensional velocity flow field in a 0.914m by 1.83m plane. PIV measurements were performed upstream and downstream of the vertical tail section and are compared to simulation results. The simulations are also used to better understand the tunnel wall and body/rotor support effects by comparing simulations with and without tunnel floor/ceiling walls and supports. Comparisons are also made to the experimental force and moment data for the body and rotor.

Allan, Brian G.↗
NASA NTRSDate not supplied

Verification Test Suite for Spalart-Allmaras QCR2000 Turbulence Model

The paper presents three benchmark cases for verification of Reynolds-averaged Navier-Stokes solvers. The verification studies focus on a one-equation Spalart-Allmaras model, SA-[neg]-QCR2000, that uses a version of quadratic constitutive relations.The benchmark cases are a two-dimensional subsonic flow around a Joukowski airfoil, a three-dimensional supersonic flow through a square duct, and a three-dimensional flow over a wing-fuselage configuration. The turbulence-model formulation, geometry, flow conditions, grids, and the expected output are described in detail. Reference solutions computed by several established codes are shown

Computational Aerodynamics↗
NASA NTRS2003

Improved models for the evolution of the Coma cluster of galaxies

The analysis by Fitchett & Webster of the observations of the Coma cluster of galaxies has demonstrated that the center of the Coma Cluster consists of two subclusters. Therefore, it is important to construct realistic dynamical models of a galaxy cluster with two mass centers. Our previous N-body models for the Coma Cluster consisted of point masses or particles with simple interaction properties. In the current paper, we employ a more sophisticated N-body code, which includes dynamical friction, mass exchange, and mergers between galaxies.

Dark↗
DOE OSTI2020

Key Technical Issues for Greater-Than-Class-C (GTCC) Waste Disposal - 20164

Greater-than-Class C (GTCC) waste is low-level waste (LLW) that exceeds the Class C concentrations tabulated in Title 10, Code of Federal Regulations (CFR), Part 61. Disposal of GTCC waste in the near surface (i.e. upper 30 m of the earth's surface) is prohibited in the United States. Though GTCC waste disposal is generally prohibited, the Commission can approve disposal according to 10 CFR 61.55(a)(2)(iv) on a case-by-case basis. The waste classification tables, developed in the early 1980's, were based on model projections of dose to an inadvertent intruder in agricultural, construction, and discovery exposure scenarios. The assumptions and parameters used were documented in a series of public reports. Though different disposal facility designs and site conditions were considered, the waste classification tables were based on a disposal facility design located in a specific environment. The tables were based on shallow (i.e., top few meters) trench facility designs and did not consider deeper facilities. To determine the suitability of GTCC waste disposal in the near-surface, site-specific analyses must account for differences between GTCC waste and Class A, B, and C LLW (hereafter, traditional LLW). GTCC waste can have concentrations of radionuclides that are much higher than traditional LLW. Because of these higher concentrations, processes that are typically not significant for traditional LLW may be significant with respect to disposal of GTCC waste. These processes include, but are not limited to, heat generation, criticality, and radiolysis. The form of the waste as well as the barriers to release of the waste (e.g. waste package) could be substantially different than they are for traditional LLW. These barriers need to be considered when assessing the impacts of accidents during receipt and placement and in evaluating long-term performance. Some GTCC radioactivity is either embedded in stainless steel or contained in stainless steel barriers. Stainless steel can have very low corrosion rates under a variety of environmental conditions. Finally, and possibly most importantly, GTCC waste would likely need to be disposed deeper than traditional LLW to reduce the probability of disturbance. The waste classification table values of traditional LLW are based on an inadvertent intruder excavating into the waste and bringing some of the material to the land surface. If waste is deeper than approximately 5 m, the excavation scenario becomes very unlikely. Therefore, other intruder scenarios, such as drilling exposure scenarios, need to be evaluated. This paper summarizes key technical issues for the disposal of GTCC waste. Our previous study shows that certain GTCC waste may be suitable for near surface disposal whereas others may not. This study may support ongoing technical analyses assessing potential disposal of GTCC in a near surface disposal facility. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Electric Readiness in Residential Energy Code: Technical Brief

This technical brief provides requirements for electric readiness that could be incorporated into model residential energy codes. It provides background on the basis and benefits of the provisions, and model code language that can be plugged into the IECC or adapted into other energy codes. This brief has been updated to reflect changes in the 2024 IECC.

2021 IECC↗
NASA NTRS2018

HST/GHRS Observations of Cool, Low-Gravity Stars. VI. Mass-Loss Rates and Wind Parameters for M Giants

The photon-scattering winds of M giants absorb parts of the chromospheric emission lines and produce self-reversed spectral features in high-resolution Hubble Space Telescope (HST)/GHRS spectra. These spectra provide an opportunity to assess fundamental parameters of the wind, including flow and turbulent velocities, the optical depth of the wind above the region of photon creation, and the star’s mass-loss rate. This paper is the last paper in the series “GHRS Observations of Cool, Low-Gravity Stars”; the last several have compared empirical measurements of spectral emission lines with models of the winds and mass loss of K giants and supergiants. We have used the Sobolev with Exact Integration radiative transfer code, along with simple models of the outer atmosphere and wind, to determine and compare the wind characteristics of the two M-giant stars, gamma Cru (M3.5III) and mu Gem (M3IIIab), with previously derived values for low-gravity K-stars. The analysis specifies the wind parameters and calculates line profiles for the Mg II resonance lines, in addition to a range of unblended Fe II lines. Our line sample covers a large range of wind opacities and, therefore, probes a range of heights in the atmosphere. Our results show that mu Gem has a slower and more turbulent wind than gamma Cru. Also, mu Gem has a weaker chromosphere, in terms of surface flux, with respect to gamma Cru. This suggests that mu Gem is more evolved than gamma Cru. Comparing the two M giants in this work with previously studied K-giant and supergiant stars (alpha Tau, gamma Dra, and lambda Vel) reveals that the M giants have slower winds than the earlier giants, but exhibit higher mass-loss rates. Our results are interpreted in the context of the winds being driven by Alfvén waves.

chromospheric emission lines↗
DOE OSTI2020

A discrete integral transform for rapid spectral synthesis

Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- & Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed. A Python implementation of the method is compared against state-of-the-art spectral code RADIS, and is since adopted as RADIS's default synthesis method. The synthesis of a benchmark CO2 spectrum consisting of 1.8 M spectral lines and 200k spectral points took only 3.1 s using the proposed method (1011 lines × spectral points/s), a factor ~300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. Finally, an experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2~3 orders performance increase, achieving up to 5 ∙ 10 14 lines × spectral points/s.

42 ENGINEERING↗
DOE OSTI2022

Zero power experiment for molten salt reactor - code verification and validation request

Molten Salt Reactors (MSRs) are one of the six perspective Generation IV nuclear technologies that should fulfil sustainability, economics, safety, and proliferation resistance goals. In recent years these concepts have attracted increased interest in research as well as in industry. In the context of developing new, innovative technology, a zero-power experiment would be an ideal option for improving the understanding of MSR systems under low-risk and low-cost conditions. The given research was focused on cross-verification of the nuclear codes for MSR modelling and simulations as required for the development of the first design of a potential facility and a draft of an experimental programme. The analysis included the evaluation of lattice and Monte Carlo codes. The study showed that all considered nuclear tools were in good agreement with each other. Although this does not mean the results will match the real experimental data, which is essential for the regulatory approval of a novel technology. Currently, experimental data is almost absent for MSR technology. Hence, a zero power experiment is needed to assure we can rely on modelling and simulation studies for further development of MSR technology and deliver real-world experiments to support the licensing of future MSRs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Electrostatic Relativistic Fluid Models of Electron Emission in a Warm Diode

A semi-analytic fluid model has been developed for characterizing relativistic electron emission across a warm diode gap. Here we demonstrate the use of this model in (i) verifying multi-fluid codes in modeling compressible relativistic electron flows (the EMPIRE-Fluid code is used as an example; see also Ref. 1), (ii) elucidating key physics mechanisms characterizing the influence of compressibility and relativistic injection speed of the electron flow, and (iii) characterizing the regimes over which a fluid model recovers physically reasonable solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Variations in lithium vapor cave performance predictions due to radial transport and recycling assumptions

The lithium vapor cave is a detached divertor design that uses a single private flux region baffle to contain a dense cloud of lithium vapor to dissipate heat flux. Plasma flows are created via fuel gas puffing in order to minimize lithium contamination of the main plasma. Significant modeling using the 2D edge code SOLPS-ITER has already been performed, predicting that sufficient target heat flux reductions (q$^{max}_{Target}$ MW m −2 ) with ( $n$ Li /$n$ e ) LCFS < 0.05 is possible in a case with 90 MW m −2 unmitigated heat flux in NSTX-U. Low heat flux and low upstream concentration was found with a variety of combinations of divertor geometries, target recycling coefficients, upstream plasma parameters, lithium evaporation locations and deuterium fueling locations, with variations in performance found for each design choice. However, the most universal uncertainty of SOLPS-ITER simulations has until now remained unaddressed systematically, namely the cross-field anomalous particle and heat diffusivities. This article aims to bound the uncertainty in lithium concentration prediction as a result of the assumed deuterium cross-field transport and recycling. For the simulations presented here, a factor of 2.1 increase in the upstream lithium density prediction across a factor of four decrease to the assumed deuterium radial particle diffusivity is found. This result is compared across different assumed deuterium recycling coefficients, known to reduce with lithium injection. Upstream lithium density is found to vary by a factor of 2.4 across a feasible range of recycling coefficients at similar $n$$^{OMP,sep}_{e}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inverse bremsstrahlung absorption rate for super-Gaussian electron distribution functions including plasma screening

Here we provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)], and we find the absorption rate under these conditions is increased by as much as ≈ 30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z = 6 and Te = 400 eV the Langdon prediction for the absorption is in error by a factor of ≈ 2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to ≈ 30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order “M” and the screening length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗