SEARCH · Engineering Papers

Results for “F 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Nuclear modules for space electric propulsion

Analysis of interplanetary cargo and piloted missions requires calculations of the performances and masses of subsystems to be integrated in a final design. In a preliminary and scoping stage the designer needs to evaluate options iteratively by using fast computer simulations. The Oak Ridge National Laboratory (ORNL) has been involved in the development of models and calculational procedures for the analysis (neutronic and thermal hydraulic) of power sources for nuclear electric propulsion. The nuclear modules will be integrated into the whole simulation of the nuclear electric propulsion system. The vehicles use either a Brayton direct-conversion cycle, using the heated helium from a NERVA-type reactor, or a potassium Rankine cycle, with the working fluid heated on the secondary side of a heat exchanger and lithium on the primary side coming from a fast reactor. Given a set of input conditions, the codes calculate composition. dimensions, volumes, and masses of the core, reflector, control system, pressure vessel, neutron and gamma shields, as well as the thermal hydraulic conditions of the coolant, clad and fuel. Input conditions are power, core life, pressure and temperature of the coolant at the inlet of the core, either the temperature of the coolant at the outlet of the core or the coolant mass flow and the fluences and integrated doses at the cargo area. Using state-of-the-art neutron cross sections and transport codes, a database was created for the neutronic performance of both reactor designs. The free parameters of the models are the moderator/fuel mass ratio for the NERVA reactor and the enrichment and the pitch of the lattice for the fast reactor. Reactivity and energy balance equations are simultaneously solved to find the reactor design. Thermalhydraulic conditions are calculated by solving the one-dimensional versions of the equations of conservation of mass, energy, and momentum with compressible flow.

Difilippo, F. C.↗

NASA NTRS2000

An Analysis of Water Line Profiles in Star Formation Regions Observed by SWAS

We present spectral line profiles for the 557 GHz 1(sub 1,0) yields 1(sub 0,1) ground-state rotational transition of ortho-H2(16)O for 18 galactic star formation regions observed by SWAS. 2 Water is unambiguously detected in every source. The line profiles exhibit a wide variety of shapes, including single-peaked spectra and self-reversed profiles. We interpret these profiles using a Monte Carlo code to model the radiative transport. The observed variations in the line profiles can be explained by variations in the relative strengths of the bulk flow and small-scale turbulent motions within the clouds. Bulk flow (infall, outflow) must be present in some cloud cores, and in certain cases this bulk flow dominates the turbulent motions.

Ashby, Matthew L. N.↗

NASA NTRS2006

CDUCT-LaRC Status - Shear Layer Refraction and Noise Radiation

A proposed boundary condition accounting for shear layer effects within the Ffowcs Williams-Hawkings radiation module of the CDUCT-LaRC code is investigated. The development and numerical justification of the boundary condition formulation are reviewed. An initial assessment of the effectiveness of the shear layer correction is conducted through comparison with experimental data. Preliminary results indicate that the correction provides physically meaningful modifications of the baseline predicted directivity patterns. Trends of peak directivity steepening and shifting that appeared in predicted patterns were found to follow similar structures in measured data, particularly at higher radiation angles.

Nark, Douglas M.↗

DOE OSTI2025

Supporting Fusion at the Idaho National Laboratory

The Fusion Safety Program at Idaho National Laboratory has supported efforts to develop magnetic and inertial confinement fusion energy devices for commercial energy power development since 1979. In addition to our DOE less than hazard category 3 tritium handling facility STAR, we also have unique capabilities such as an adaptation of the severe accident code MELCOR developed specifically for fusion, and integrated modeling capabilities using the MOOSE Finite Element Code. In this presentation, we cover our unique experimental capabilities for investigating tritium sorption, solubility and permeation including in molten salts. We also show the application of MELCOR-TMAP to predict thermal hydraulics and tritium migration in normal and accident conditions. Finally, we demonstrate high-fidelity modeling capabilities found within the open-source MOOSE finite element framework as well as fuel-cycle modeling using TMAP8.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Free-free matrix-elements for two-photon opacity

Recent opacity measurements have inspired a close study of the two-photon contributions to the opacity of hot plasmas. The absorption and emission of radiation is controlled by dipole matrix-elements of electrons in an atom or ion. This paper describes two independent methods to calculate matrix-elements needed for the two-photon opacity and tests the results by the f-sum rule. The usual f-sum rule is extended to a matrix f-sum that offers a rigorous test for bound-bound, bound-free and free-free transitions. An additional higher-order sum-rule for the two-photon transition amplitudes is described. In this work, we obtain a simple parametric representation of a key plasma density effect on the matrix-elements. The perturbation theory calculation of two-photon cross-sections is compared to an independent method based on the solution of the time-dependent Schroedinger equation for an atom or ion in a high-frequency electromagnetic field. This is described as a high frequency Stark effect or AC Stark effect. Two-photon cross sections calculated with the AC Stark code agree with perturbation theory to within about 5%. In addition to this cross check, the AC Stark code is well suited to evaluating important questions such as the variation of two-photon opacity for different elements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Long-duration direct drive hydrodynamics experiments on the National Ignition Facility: Platform development and numerical modeling with CHIC

Here, we report on a novel planar direct-drive platform for hydrodynamics experiments on the National Ignition Facility (NIF). Its commissioning has been performed as part of the NIF Discovery Science Program. This platform enables the use of a 30 ns drive at an average intensity of 200 TW/cm 2 , creating a planar shock and ablation front over a 2 mm radius. To benchmark the performance of this design, the planarity of both the shock and ablation fronts has been measured between 26 ns and 28 ns after the start of the laser drive in a 3 mm-thick CH foil. The platform was then used to measure late-time Rayleigh-Taylor instability (RTI) growth at the ablation front for a 2D-rippled 300 μm-thick CH foil. Simultaneously, a numerical platform has been developed with the CHIC radiation hydrodynamics code at the CELIA laboratory. The CHIC numerical platform allows, for the first time, a complete simulation of the experiments over 30 ns to be performed. Large-scale simulations recover the trajectory and the 2D RTI growth measurements. They are additionally compared with half-mode simulations performed with identical parameters. We demonstrate that both numerical techniques fit with analytical modeling of RTI growth and discuss plans for future campaigns.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗

NASA NTRS1990

Thermal interaction between an impinging hot jet and a conducting solid surface

Powered-lift aircraft may produce severe high-temperature environments which are potentially damaging to a landing surface or the aircraft. The interaction betweean the high temperature flow field and a nonadiabatic landing surface is analyzed with a coupled computational fluid dynamics/solid thermal conduction computer code, HOTJET. The HOTJET code couples time-accurate, implicit, factored solution schemes for the governing fluid dynamics equations (Reynolds-averaged Navier-Stokes equations) to the unsteady thermal conduction equation, which governs heat flux within a solid. HOTJET is validated against exact solutions to the thermal conduction and Navier-Stokes equations. First-of-a-kind results are included which show the impact of surface material properties on the fluid physics and the coupled fluid/material thermal fields.

Abeloff, P. A.↗

NASA NTRS1991

Performance of popular turbulence models for attached and separated adverse pressure gradient flows

The performance of four popular eddy-viscosity turbulence models under adverse pressure gradient conditions is investigated. The Baldwin-Lomax, the Johnson-King, the Baldwin-Barth, and the Wilcox-omega models have been implemented into the INS code, which solves the incompressible Reynolds-averaged Navier-Stokes equations. Results are shown for the well known Samuel-Joubert flow and two new flowfields, recently reported by D. M. Driver (1991). The two new flowfields pose a stronger test of the models than the Samuel-Joubert flow, because of the more severe retardation of the boundary layer, including separation in one case. A detailed comparison of the numerical results and the experimental data is shown.

Menter, F. R.↗

NASA NTRS2001

NAS Grid Benchmarks: A Tool for Grid Space Exploration

We present an approach for benchmarking services provided by computational Grids. It is based on the NAS Parallel Benchmarks (NPB) and is called NAS Grid Benchmark (NGB) in this paper. We present NGB as a data flow graph encapsulating an instance of an NPB code in each graph node, which communicates with other nodes by sending/receiving initialization data. These nodes may be mapped to the same or different Grid machines. Like NPB, NGB will specify several different classes (problem sizes). NGB also specifies the generic Grid services sufficient for running the bench-mark. The implementor has the freedom to choose any specific Grid environment. However, we describe a reference implementation in Java, and present some scenarios for using NGB.

Frumkin, Michael↗

NASA NTRS2004

How Consistent are Recent Variations in the Tropical Energy and Water Cycle Resolved by Satellite Measurements?

One notable aspect of Earth's climate is that although the planet appears to be very close to radiative balance at top-of-atmosphere (TOA), the atmosphere itself and underlying surface are not. Profound exchanges of energy between the atmosphere and oceans, land and cryosphere occur over a range of time scales. Recent evidence from broadband satellite measurements suggests that even these TOA fluxes contain some detectable variations. Our ability to measure and reconstruct radiative fluxes at the surface and at the top of atmosphere is improving rapidly. One question is 'How consistent, physically, are these diverse remotely-sensed data sets'? The answer is of crucial importance to understanding climate processes, improving physical models, and improving remote sensing algorithms. In this work we will evaluate two recently released estimates of radiative fluxes, focusing primarily on surface estimates. The International Satellite Cloud Climatology Project 'FD' radiative flux profiles are available from mid-1983 to near present and have been constructed by driving the radiative transfer physics from the Goddard Institute for Space Studies (GISS) global model with ISCCP clouds and TOVS (TIROS Operational Vertical Sounder)thermodynamic profiles. Full and clear sky SW and LW fluxes are produced. A similar product from the NASA/GEWEX Surface Radiation Budget Project using different radiative flux codes and thermodynamics from the NASA/Goddard Earth Observing System (GEOS-1) assimilation model makes a similar calculation of surface fluxes. However this data set currently extends only through 1995. We also employ precipitation measurements from the Global Precipitation Climatology Project (GPCP) and the Tropical Rainfall Measuring Mission (TRMM). Finally, ocean evaporation estimates from the Special Sensor Microwave Imager (SSM/I) are considered as well as derived evaporation from the NCAR/NCEP Reanalysis. Additional information is included in the original extended abstract.

Robertson, F. R.↗

DOE OSTI2021

Study of the Effect of Reflections on High-Power, 110-GHz Pulsed Gyrotron Operation

The effect of reflection is studied experimentally and theoretically on a high-power 110-GHz gyrotron operating in the TE 22,6 mode in 3 μs pulses at 96 kV, 40 A. The experimental setup allows variation of the reflected power from 0 to 33% over a range of gyrotron operating conditions. The phase of the reflection is varied by translating the reflector along the axis. Operating at a higher efficiency point, at 4.40 T with 940kW of output power, reflected power exceeding 11% causes a switch from operation in the TE 22,6 to simultaneous operation in the TE 22,6 and TE 21,6 modes with a large decrease of the total gyrotron output power. This switching effect is in good agreement with simulations using the MAGY code. Operating at a more stable point, 4.44 T with 580kW of output power, when the reflection is increased, the output power remains in the TE 22,6 mode but it decreases monotonically with increasing reflection, dropping to 200kW at 33% reflection. Furthermore, at a reflection above 22%, a power modulation at 25 to 30MHz is observed, independent of the phase of the reflected wave. Such a modulated signal may be useful in spectroscopic and other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Analysis and Simulation of One-Dimensional Transport Models for Lithium Symmetric Cells

The main objective of this paper is to analyze transport models for lithium symmetric cells and arrive at an efficient model and code to simulate the same. Two one-dimensional models are considered. The first model uses the dilute solution Nernst–Planck equations in conjunction with the electroneutrality assumption (EN-NP model). For binary electrolytes, an analytical solution for electrolyte concentration is derived and compared with numerical solutions by an approximate finite volume method. The second approach relaxes the electroneutrality assumption by way of Poisson’s equation for the electrostatic potential (PNP model). The computational difficulties of the PNP model are tackled using the approximate finite volume method, with demonstrated convergence characteristics even with bulk dimensions several orders of magnitude larger than the characteristic double layer size. A robust code is developed for the PNP model. The computationally efficient transport models can facilitate simulations, physical understanding, and analysis. This is illustrated by a case study in which these models are coupled with modified kinetic models, which are then parameterized with experimental voltage response data using a systems-level approach. The estimated parameters provide further insight into the electrochemical phenomena underpinning voltage transitions in symmetric cells.

25 ENERGY STORAGE↗

DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1977

L-band tone-code-data transponder calibration

The objectives of this program were to identify and quantify factors which affect the performance of the L-band tone-code-data ranging transponders. Specific objectives included the following: (1) assemble the L-band ranging transponder, previously deployed in Hawaii for the tracking of the ATS-5 satellite, at the GE Radio-Optical Observatory; (2) configure the observatory to conduct calibration exercises with the transponder; and (3) conduct sufficient calibration experiments to demonstrate factors which degrade transponder accuracy, precision, and reliability, to quantify these factors where possible, and to verify long term transponder stability under controlled conditions.

Brisken, A. F.↗

NASA NTRS1991

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form were iterated in time until a steady solution was reached. Evidence was given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at the domain boundaries was proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects were accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. The flow in a developing S-duct was then solved in the laminar regime in a Reynolds number (Re) of 790 and in the turbulent regime at Re equals 40,000 by using the Baldwin-Lomax model. The Stanitz elbow was then solved by using an invicid version of the same code at M sub inlet equals 0.4. Grid dependence and convergence rate were investigated, showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re equals 2.5 times 10(exp 6) was solved with the Baldwin-Lomax and the q-omega models. Both approaches show satisfactory agreement with experiments, although the q-omega model was slightly more accurate.

Michelassi, V.↗

NASA NTRS1991

Hypersonic waverider analysis - A comparison of numerical and experimental results

In this study two waverider configurations are analyzed numerically and experimentally. The geometries are designed using an inverse design and optimization code with the optimization procedure driven by maximized L/D and limited by geometric constraints. Euler and Navier-Stokes numerical simulations are performed with a partially flux-split finite-difference code developed at NASA Ames, F3D. Experimental data is obtained at NASA Langley in the Unitary Plan Wind Tunnel (UPWT) and the 20-Inch Mach 6 wind tunnel for flow visualization, force/moment, and surface pressure measurements at both on- and off-design conditions. Euler simulations are compared with analytic results obtained from conical flow theory, and Navier-Stokes simulations are compared with experimental results both qualitatively and quantitatively for on- and off-design conditions.

Jones, K. D.↗

NASA NTRS1994

Operating in the space plasma environment: A spacecraft charging study of the Solar X-ray Imager

This study presents the results of a spacecraft charging effects protection study conducted on the Solar X-ray Imager (SXI). The SXI is being developed by NASA Marshall Space Flight Center for NOAA's Space Environment Laboratory, and will be used to aid in forecasting energetic particle events and geomagnetic storms. Images will provide information on the intensity and location of solar flares, coronal mass ejections, and high speed solar streams. The SXI will be flown on a next-generation GOES sometime in the mid to late 1990's. Charging due to the encounter with a worst-case magnetic substorm environment is modeled using the NASCAP/GEO computer code. Charging levels of exterior surfaces and the floating potential of the spacecraft relative to plasma are determined as a function of spacecraft design, operational configuration, and orbital conditions. Areas where large surface voltage gradients exist on or near the SXI are identified as possible arc-discharge sites. Results of the charging analysis are then used to develop design recommendations that will limit the effects of spacecraft charging on the SXI operation.

Herr, Joel L.↗