SEARCH · Engineering Papers

Results for “surface 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,261 records · Page 70

All sources NASA NTRS DOE OSTI

NASA NTRS1994

A Combined Experimental/Computational Study of Flow in Turbine Blade Cooling Passage

Laser velocimetry was utilized to map the velocity field in a serpentine turbine blade cooling passage at Reynolds and Rotation numbers of up to 25.000 and 0.48. These results were used to assess the combined influence of passage curvature and Coriolis force on the secondary velocity field generated. A Navier-Stokes code (NASTAR) was validated against incompressible test data and then used to simulate the effect of buoyancy. The measurements show a net convection from the low pressure surface to high pressure surface. The interaction of the secondary flows induced by the turns and rotation produces swirl at the turns, which persisted beyond 2 hydraulic diameters downstream of the turns. The incompressible flow field predictions agree well with the measured velocities. With radially outward flow, the buoyancy force causes a further increase in velocity on the high pressure surface and a reduction on the low pressure surface. The results were analyzed in relation to the heat transfer measurements of Wagner et al. (1991). Predicted heat transfer is enhanced on the high pressure surfaces and in turns. The incompressible flow simulation underpredicts heat transfer in these locations. Improvements observed in compressible flow simulation indicate that the buoyancy force may be important.

Tse, D. G. N.↗

NASA NTRS1995

Evaluation of numerical schemes for the analysis of sound generation by blade-gust interaction

In this investigation three different numerical algorithms have been utilized to compute the flow about a flat plate in the presence of a transverse gust described by a sinusoidal disturbance. The three schemes include the MacCormack explicit finite difference scheme which is second order accurate in both time and space, the Gottlieb and Turkel modification of MacCormack's scheme which is fourth order accurate in space and second order accurate in time, (referred to as the 2-4 scheme), and a two step scheme developed by Bayliss et al. which has second order temporal accuracy and sixth order spatial accuracy (a 2-6 scheme). The flow field results are obtained with these schemes by using the same code with the only difference being the implementation of the respective solution algorithms. The problem is set up so that the sinusoidal disturbance is imposed at the surface of the flat plate as a surface boundary condition. Thus the problem is treated as scattering problem. The computed results include the time average of the acoustic pressure squared along grid lines five points away from the boundaries; distribution throughout the computational domain is monitored at various times. The numerical results are compared with an exact solution obtained by Atassi, Dusey, and Davis.

Scott, J. N.↗

NASA NTRS1995

In-flight radiometric calibration of AVIRIS in 1994

The AVIRIS sensor must be calibrated at the time it measures spectra from the ER-2 airborne platform in order to achieve research and application objectives that are both quantitative and physically based. However, the operational environment inside the Q-bay of the ER-2 at 20 km altitude differs from that in the AVIRIS laboratory with respect to temperature, pressure, vibration, and high-frequency electromagnetic fields. Experiments at surface calibration targets are used in each flight season to confirm the accuracy of AVIRIS in-flight radiometric calibrations. For these experiments, the MODTRAN radiative transfer code is constrained by using in situ measurements to independently predict the upwelling spectral radiance arriving at AVIRIS for a specific calibration target. AVIRIS calibration is validated in flight by comparing the MODTRAN-predicted radiance to the laboratory-calibrated radiance measured by the AVIRIS sensor for the same time over the calibration target. We present radiometric calibration results for the AVIRIS in-flight calibration experiment held at the beginning of the 1994 flight season.

Green, Robert O.↗

NASA NTRS1995

Imaging spectrometer measurement of water vapor in the 400 to 2500 nm spectral region

The Airborne Visible-Infrared Imaging Spectrometer (AVIRIS) measures the total upwelling spectral radiance from 400 to 2500 nm sampled at 10 nm intervals. The instrument acquires spectral data at an altitude of 20 km above sea level, as images of 11 by up to 100 km at 17x17 meter spatial sampling. We have developed a nonlinear spectral fitting algorithm coupled with a radiative transfer code to derive the total path water vapor from the spectrum, measured for each spatial element in an AVIRIS image. The algorithm compensates for variation in the surface spectral reflectance and atmospheric aerosols. It uses water vapor absorption bands centered at 940 nm, 1040 nm, and 1380 nm. We analyze data sets with water vapor abundances ranging from 1 to 40 perceptible millimeters. In one data set, the total path water vapor varies from 7 to 21 mm over a distance of less than 10 km. We have analyzed a time series of five images acquired at 12 minute intervals; these show spatially heterogeneous changes of advocated water vapor of 25 percent over 1 hour. The algorithm determines water vapor for images with a range of ground covers, including bare rock and soil, sparse to dense vegetation, snow and ice, open water, and clouds. The precision of the water vapor determination approaches one percent. However, the precision is sensitive to the absolute abundance and the absorption strength of the atmospheric water vapor band analyzed. We have evaluated the accuracy of the algorithm by comparing several surface-based determinations of water vapor at the time of the AVIRIS data acquisition. The agreement between the AVIRIS measured water vapor and the in situ surface radiometer and surface interferometer measured water vapor is 5 to 10 percent.

Green, Robert O.↗

NASA NTRS1998

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations that include transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Using several simplifying assumptions, linearized, stability equations are derived that describe the response of the plate just after bifurcation buckling occurs. These equations are then modified to allow the plate reference surface to be located a distance z(c), from the centroid surface which is convenient for modeling stiffened-plate assemblies. The implementation of the new theory into the VICONOPT buckling and vibration analysis and optimum design program code is described. Either classical plate theory (CPT) or first-order shear-deformation plate theory (SDPT) may be selected in VICONOPT. Comparisons of numerical results for several example problems with different loading states are made. Results from the new curved-plate analysis compare well with closed-form solution results and with results from known example problems in the literature. Finally, a design-optimization study of two different cylindrical shells subject to uniform axial compression is presented.

McGowan, David M.↗

NASA NTRS1997

Physics of Rotation Powered Pulsars and Their Nebulae

This letter is my progress report for the Astrophysics Theory grant. The first goal of the research supported by this grant is to produce a time dependent theory of the unsteady relativistic collisionless shock wave terminating the relativistic wind from a pulsar, and compare the predicted surface brightness fluctuations to Hubble Space Telescope observations of the wisps in the Crab Nebula. The second goal is to model the production of electron-positron pairs over the polar caps of rotation powered pulsars, and use the results to predict the heating of the surface due to particle trapping and bombardment of the atmosphere at the base of the polar field lines. We have succeeded in creating a one dimensional hybrid code, in which the electron-positron pairs incident on the shock structure are modeled as a relativistic, adiabatic, ideal MHD fluid, while the heavy ions are treated as particles using a particle-in-cell algorithm. The electromagnetic fields are evaluated from the currents and charge densities in the pairs and the ions, while the particles and the fluid accelerate in response to the computed self-consistent electromagnetic fields. The results are promising, in that the underlying ion cyclotron instability generates finite amplitude, propagating magnetosonic waves in the pairs, whose wavelengths and frequencies, when translated into physical units, are comparable to the observed running waves of brightness observed by HST near the Crab pulsar. The code is undergoing a number of tests, to assure us that this preliminary correspondence is not an artifact. In the coming year, the observational appearance of the models will be computed and compared to the HST observations of the Crab now in hand, and used to predict the HST results which will be obtained the year after next. WE also developed a one dimensional cascade theory for pair creation over pulsars' polar caps. A linear integral equation describing the synchrotron cascade has been derived and solved by iterative techniques, in the case when a high energy electron moving parallel to a star centered dipole magnetic field initiates the cascade through curvature gamma ray emission.

Arons, Jonathan↗

NASA NTRS2005

A Hexapod Robot to Demonstrate Mesh Walking in a Microgravity Environment

The JPL Micro-Robot Explorer (MRE) Spiderbot is a robot that takes advantage of its small size to perform precision tasks suitable for space applications. The Spiderbot is a legged robot that can traverse harsh terrain otherwise inaccessible to wheeled robots. A team of Spiderbots can network and can exhibit collaborative efforts to SUCCeSSfUlly complete a set of tasks. The Spiderbot is designed and developed to demonstrate hexapods that can walk on flat surfaces, crawl on meshes, and assemble simple structures. The robot has six legs consisting of two spring-compliant joints and a gripping actuator. A hard-coded set of gaits allows the robot to move smoothly in a zero-gravity environment along the mesh. The primary objective of this project is to create a Spiderbot that traverses a flexible, deployable mesh, for use in space repair. Verification of this task will take place aboard a zero-gravity test flight. The secondary objective of this project is to adapt feedback from the joints to allow the robot to test each arm for a successful grip of the mesh. The end result of this research lends itself to a fault-tolerant robot suitable for a wide variety of space applications.

Foor, David C.↗

NASA NTRS2012

As-Flown Shuttle Micro-Meteoroid Orbital Debris (MMOD) Assessment History

The final flight of the Space Shuttle program was completed on July 21st of 2011. The reusable nature of the Orbiter vehicles has provided NASA with a unique opportunity to inspect and sample spacecraft surfaces exposed to the MMOD environment in low Earth orbit. This paper will provide details of MMOD damages found on Orbiter surfaces after recent flights, as well as putting these damages in context by providing historical comparisons of recent damage to previous flights. The Bumper threat assessment computer code is used by NASA to determine spacecraft MMOD penetration risk. The tool was used before each mission to estimate pre-flight MMOD risk for the Shuttle Program Office using planned vehicles attitudes and exposure times. This paper summarizes the efforts of the authors to assess the expected number MMOD impacts of selected particles sizes of 52 shuttle missions ranging from STS-50 (1992) through STS-133 (2011) using as-flown attitude data. The missions encompass the Shuttle/Mir precursors to the International Space Station (ISS) as well as the ISS assembly flights. The paper includes a comparison of observed MMOD impact damage to the predicted number of impacts from Bumper code assessments. This data can be used to validate engineering models of the orbital debris and micrometeoroid environment, such NASA s ORDEM debris model and MEM meteoroid model.

Hyde, James L.↗

NASA NTRS2011

Effects of Non-Equilibrium Chemistry and Darcy-Forchheimer Flow of Pyrolysis Gas for a Charring Ablator

The Fully Implicit Ablation and Thermal Response code, FIAT, simulates pyrolysis and ablation of thermal protection materials and systems. The governing equations, which include energy conservation, a three-component decomposition model, and a surface energy balance, are solved with a moving grid. This work describes new modeling capabilities that are added to a special version of FIAT. These capabilities include a time-dependent pyrolysis gas flow momentum equation with Darcy-Forchheimer terms and pyrolysis gas species conservation equations with finite-rate homogeneous chemical reactions. The total energy conservation equation is also enhanced for consistency with these new additions. Parametric studies are performed using this enhanced version of FIAT. Two groups of analyses of Phenolic Impregnated Carbon Ablator (PICA) are presented. In the first group, an Orion flight environment for a proposed Lunar-return trajectory is considered. In the second group, various test conditions for arcjet models are examined. The central focus of these parametric studies is to understand the effect of pyrolysis gas momentum transfer on PICA material in-depth thermal responses with finite-rate, equilibrium, or frozen homogeneous gas chemistry. Results are presented, discussed, and compared with those predicted by the baseline PICA/FIAT ablation and thermal response model developed by the Orion Thermal Protection System Advanced Development Project.

Ablation↗

NASA NTRSDate not supplied

Recent Improvements to the LAURA and HARA Codes

This paper describes recent improvements to the LAURA and HARA codes. LAURA is a CFD code for aerothermodynamics, and HARA evaluates the shock-layer radiation that provides the radiative source term for the flowfield energy equations and radiative heating to a surface. The next release of LAURA and HARA includes a variety of new capabilities. These new capabilities include an automated uncertainty quantification workflow for radiative heat transfer, options for specifying surface roughness and turbulent transition location in the algebraic turbulence models, and improved grid and solution interpolation techniques. Additionally, the computational efficiency of both LAURA and HARA have been improved. Optimization of the MPI communication routines in LAURA are shown to improve the parallel efficiency of the primary flow when running with multiple processes per block, and recent optimization of HARA leverage graphics processing unit (GPU) acceleration in the radiation calculations. Using GPU acceleration of HARA is shown to decrease the cost of the radiation line-of-sight calculation by approximately one order of magnitude for a 10.5 km/s Earth entry simulation.

LAURA HARA CFD 5.6↗

NASA NTRS2015

Superboom Caustic Analysis and Measurement Program (SCAMP) Final Report

The objectives of the Superboom Caustic Analysis and Measurement (SCAMP) Program were to develop and validate, via flight-test measurements, analytical models for sonic boom signatures in and around focal zones as they are expected to occur during commercial aircraft transition from subsonic to supersonic flight, and to apply these models to focus boom prediction of low-boom aircraft designs. The SCAMP program has successfully investigated sonic boom focusing both analytically and experimentally, while gathering a comprehensive empirical flight test and acoustic dataset, and developing a suite of focused sonic boom prediction tools. An experimental flight and acoustic measurement test was designed during the initial year of the SCAMP program, with execution of the SCAMP flight test occurring in May 2011. The current SCAMP team, led by Wyle, includes partners from the Boeing Company, Pennsylvania State University, Gulfstream Aerospace, Eagle Aeronautics, and Central Washington University. Numerous collaborators have also participated by supporting the experiment with human and equipment resources at their own expense. The experiment involved precision flight of a McDonnell Douglas (now Boeing) F-18B executing different maneuvers that created focused sonic booms. The maneuvers were designed to center on the flight regime expected for commercial supersonic aircraft transonic transition, and also span a range of caustic curvatures in order to provide a variety of conditions for code validations. The SCAMP experiment was designed to capture concurrent F-18B on-board flight instrumentation data, high-fidelity ground-based and airborne acoustic data, and surface and upper air meteorological data. Close coordination with NASA Dryden resulted in the development of new experimental instrumentation and techniques to facilitate the SCAMP flight-test execution, including the development of an F-18B Mach rate cockpit display, TG-14 powered glider in-flight sonic boom measurement instrumentation and "Where's the Focus?" (WTF) software for near-real time way-point computation accounting for local atmospherics. In May 2011, 13 F-18B flights were conducted during 5 flying days over a 2 week period. A densely populated 10,000 ft-long ground acoustic array with 125-ft microphone spacing was designed to capture pre-, focus, and post-focus regions. The ground-based acoustic array was placed in a nominally east-west orientation in the remote Cuddeback lakebed region, north of Edwards AFB. This area was carefully selected to avoid placing focused booms on populated areas or solar power facilities. For the SCAMP measurement campaign, approvals were obtained to temporarily extend the Black Mountain supersonic corridor northward by three miles. The SCAMP flight tests successfully captured 70 boom events, with 61 focus passes, and 9 calibration passes. Seventeen of the focus passes and three of the calibration passes were laterally offset; with the others being centerline flights. Airborne incoming sonic boom wave measurements were measured by the TG-14 for 10 of the F-18B flight passes including one maximum focus signature, several N-u combinations, several overlapped N-u signatures, and several evanescent waves. During the 27-month program, the SCAMP team developed a suite of integrated computer codes with sonic boom focusing predictive capabilities: PCBoom, Lossy Nonlinear Tricomi Equation Method (LNTE) and the Nonlinear Progressive wave Equation (NPE) method. PCBoom propagates the rays through the atmosphere and, in addition to legacy focus signature prediction based on the Gill-Seebass method, provides input source characteristics and propagation parameters to LNTE and NPE. LNTE, a Tricomi solver that incorporates atmospheric losses, computes the focus signature at the focus, and computes the focus signature in the vicinity of the focal zone, including the evanescent and post-focus zones. LNTE signature auralization from low-boom vehicle designs has been demonstrated in the NASA Langley Interior Effects Room (IER). The NPE has also been validated for use in prediction of focused ground boom signatures in sonic boom focal zones. The NPE formulation has the capability to incorporate atmospheric turbulence in the predictions. This has been applied to sonic boom propagation in the past. Prediction of turbulence effects on focal zone signatures was not, however, explored during the SCAMP program.

Page, Juliet↗

NASA NTRS1975

The NATA code; theory and analysis. Volume 2: User's manual

The NATA code is a computer program for calculating quasi-one-dimensional gas flow in axisymmetric nozzles and rectangular channels, primarily to describe conditions in electric archeated wind tunnels. The program provides solutions based on frozen chemistry, chemical equilibrium, and nonequilibrium flow with finite reaction rates. The shear and heat flux on the nozzle wall are calculated and boundary layer displacement effects on the inviscid flow are taken into account. The program contains compiled-in thermochemical, chemical kinetic and transport cross section data for high-temperature air, CO2-N2-Ar mixtures, helium, and argon. It calculates stagnation conditions on axisymmetric or two-dimensional models and conditions on the flat surface of a blunt wedge. Included in the report are: definitions of the inputs and outputs; precoded data on gas models, reactions, thermodynamic and transport properties of species, and nozzle geometries; explanations of diagnostic outputs and code abort conditions; test problems; and a user's manual for an auxiliary program (NOZFIT) used to set up analytical curvefits to nozzle profiles.

Bade, W. L.↗

NASA NTRS1992

Improvements to the missile aerodynamic prediction code DEMON3

The computer program DEMON3 was developed for the aerodynamic analysis of nonconventional supersonic configurations comprising a body with noncircular cross section and up to two wing or fin sections. Within a wing or fin section, the lifting surfaces may be cruciform, triform, planar, or low profile layouts; the planforms of the lifting surfaces allow for breaks in sweep. The body and fin sections are modeled by triplet and constant u-velocity panels, respectively, accounting for mutual body-fin interference. Fin thickness effects are included for the use of supersonic planar source panels. One of the unique features of DEMON3 is the modeling of high angle of attack vortical effects associated with the lifting surfaces and the body. In addition, shock expansion and Newtonian pressure calculation methods can be optionally engaged. These two dimensional nonlinear methods are augmented by aerodynamic interference determined from the linear panel methods. Depending on the geometric details of the body, the DEMON3 program can be used to analyze nonconventional configurations at angles of attack up to 25 degrees for Mach numbers from 1.1 to 6. Calculative results and comparisons with experimental data demonstrate the capabilities of DEMON3. Limitations and deficiencies are listed.

Dillenius, Marnix F. E.↗

NASA NTRS1994

User's manual for the NASA Lewis ice accretion/heat transfer prediction code with electrothermal deicer input

A version of LEWICE has been developed that incorporates a recently developed electrothermal deicer code, developed at the University of Toledo by William B. Wright. This was accomplished, in essence, by replacing a subroutine in LEWICE, called EBAL, which balanced the energies at the ice surface, with a subroutine called UTICE. UTICE performs this same energy balance, as well as handles all the time-timperature transients below the ice surface, for all of the layers of a composite blade as well as the ice layer itself. This new addition is set up in such a fashion that a user may specify any number of heaters, any heater chordwise length, and any heater gap desired. The heaters may be fired in unison, or they may be cycled with periods independent of each other. The heater intensity may also be varied. In addition, the user may specify any number of layers and thicknesses depthwise into the blade. Thus, the new addition has maximum flexibility in modeling virtually any electrothermal deicer installed into any airfoil. It should be noted that the model simulates both shedding and runback. With the runback capability, it can simulate the anti-icing mode of heater performance, as well as detect icing downstream of the heaters due to runback in unprotected portions of the airfoil. This version of LEWICE can be run in three modes. In mode 1, no conduction heat transfer is modeled (which would be equivalent to the original version of LEWICE). In mode 2, all heat transfer is considered due to conduction but no heaters are firing. In mode 3, conduction heat transfer where the heaters are engaged is modeled, with subsequent ice shedding. When run in the first mode, there is virtually identical agreement with the original version of LEWICE in the prediction of accreted ice shapes. The code may be run in the second mode to determine the effects of conduction on the ice accretion process.

Masiulaniec, Konstanty C.↗

NASA NTRS1992

TIGER: Turbomachinery interactive grid generation

A three dimensional, interactive grid generation code, TIGER, is being developed for analysis of flows around ducted or unducted propellers. TIGER is a customized grid generator that combines new technology with methods from general grid generation codes. The code generates multiple block, structured grids around multiple blade rows with a hub and shroud for either C grid or H grid topologies. The code is intended for use with a Euler/Navier-Stokes solver also being developed, but is general enough for use with other flow solvers. TIGER features a silicon graphics interactive graphics environment that displays a pop-up window, graphics window, and text window. The geometry is read as a discrete set of points with options for several industrial standard formats and NASA standard formats. Various splines are available for defining the surface geometries. Grid generation is done either interactively or through a batch mode operation using history files from a previously generated grid. The batch mode operation can be done either with a graphical display of the interactive session or with no graphics so that the code can be run on another computer system. Run time can be significantly reduced by running on a Cray-YMP.

Soni, Bharat K.↗

NASA NTRS2001

Apparatus and process for freeform fabrication of composite reinforcement preforms

A solid freeform fabrication process and apparatus for making a three-dimensional reinforcement shape. The process comprises the steps of (1) operating a multiple-channel material deposition device for dispensing a liquid adhesive composition and selected reinforcement materials at predetermined proportions onto a work surface; (2) during the material deposition process, moving the deposition device and the work surface relative to each other in an X-Y plane defined by first and second directions and in a Z direction orthogonal to the X-Y plane so that the materials are deposited to form a first layer of the shape; (3) repeating these steps to deposit multiple layers for forming a three-dimensional preform shape; and (4) periodically hardening the adhesive to rigidize individual layers of the preform. These steps are preferably executed under the control of a computer system by taking additional steps of (5) creating a geometry of the shape on the computer with the geometry including a plurality of segments defining the preform shape and each segment being preferably coded with a reinforcement composition defining a specific proportion of different reinforcement materials; (6) generating programmed signals corresponding to each of the segments in a predetermined sequence; and (7) moving the deposition device and the work surface relative to each other in response to these programmed signals. Preferably, the system is also operated to generate a support structure for any un-supported feature of the 3-D preform shape.

Yang, Junsheng↗

NASA NTRS2013

Translating MAPGEN to ASPEN for MER

This software translates MAPGEN (Europa and APGEN) domains to ASPEN, and the resulting domain can be used to perform planning for the Mars Exploration Rover (MER). In other words, this is a conversion of two distinct planning languages (both declarative and procedural) to a third (declarative) planning language in order to solve the problem of faithful translation from mixed-domain representations into the ASPEN Modeling Language. The MAPGEN planning system is an example of a hybrid procedural/declarative system where the advantages of each are leveraged to produce an effective planner/scheduler for MER tactical planning. The adaptation of the planning system (ASPEN) was investigated, and, with some translation, much of the procedural knowledge encoding is amenable to declarative knowledge encoding. The approach was to compose translators from the core languages used for adapting MAGPEN, which consists of Europa and APGEN. Europa is a constraint- based planner/scheduler where domains are encoded using a declarative model. APGEN is also constraint-based, in that it tracks constraints on resources and states and other variables. Domains are encoded in both constraints and code snippets that execute according to a forward sweep through the plan. Europa and APGEN communicate to each other using proxy activities in APGEN that represent constraints and/or tokens in Europa. The composition of a translator from Europa to ASPEN was fairly straightforward, as ASPEN is also a declarative planning system, and the specific uses of Europa for the MER domain matched ASPEN s native encoding fairly closely. On the other hand, translating from APGEN to ASPEN was considerably more involved. On the surface, the types of activities and resources one encodes in APGEN appear to match oneto- one to the activities, state variables, and resources in ASPEN. But, when looking into the definitions of how resources are profiled and activities are expanded, one sees code snippets that access various information available during planning for the moment in time being planned to decide at the time what the appropriate profile or expansion is. APGEN is actually a forward (in time) sweeping discrete event simulator, where the model is composed of code snippets that are artfully interleaved by the engine to produce a plan/schedule. To solve this problem, representative code is simulated as a declarative series of task expansions. Predominantly, three types of procedural models were translated: loops, if statements, and code blocks. Loops and if statements were handled using controlled task expansion, and code blocks were handled using constraint networks that maintained the generation of results based on what the order of execution would be for a procedural representation. One advantage with respect to performance for MAPGEN is the use of APGEN s GUI. This GUI is written in C++ and Motif, and performs very well for large plans.

Rabideau, Gregg R.↗

NASA NTRS2013

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗