Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-dependent data”

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 559 records · Page 31

Analysis of Crystallization Kinetics

A realistic computer model for polymorphic crystallization (i.e., initial and final phases with identical compositions), which includes time-dependent nucleation and cluster-size-dependent growth rates, is developed and tested by fits to experimental data. Model calculations are used to assess the validity of two of the more common approaches for the analysis of crystallization data. The effects of particle size on transformation kinetics, important for the crystallization of many systems of limited dimension including thin films, fine powders, and nanoparticles, are examined.

Kelton, Kenneth F.↗

Computations of Axisymmetric Flows in Hypersonic Shock Tubes

A time-accurate two-dimensional fluid code is used to compute test times in shock tubes operated at supersonic speeds. Unlike previous studies, this investigation resolves the finer temporal details of the shock-tube flow by making use of modern supercomputers and state-of-the-art computational fluid dynamic solution techniques. The code, besides solving the time-dependent fluid equations, also accounts for the finite rate chemistry in the hypersonic environment. The flowfield solutions are used to estimate relevant shock-tube parameters for laminar flow, such as test times, and to predict density and velocity profiles. Boundary-layer parameters such as bar-delta(sub u), bar-delta(sup *), and bar-tau(sub w), and test time parameters such as bar-tau and particle time of flight t(sub f), are computed and compared with those evaluated by using Mirels' correlations. This article then discusses in detail the effects of flow nonuniformities on particle time-of-flight behind the normal shock and, consequently, on the interpretation of shock-tube data. This article concludes that for accurate interpretation of shock-tube data, a detailed analysis of flowfield parameters, using a computer code such as used in this study, must be performed.

Sharma, Surendra P.↗

Turbulence Modeling for Unsteady Transonic Flows

Conditionally sampled, ensemble-averaged velocity measurements, made with a laser velocimeter, were taken in the flowfield over the rear half of an 18% thick circular arc airfoil at zero incidence tested at M = 0.76 and at a Reynolds number based on chord of 11 x 10(exp 6). Data for one cycle of periodic unsteady flow having a reduced frequency f of 0.49 are analyzed. A series of compression waves, which develop in the early stages of the cycle, strengthen and coalesce into a strong shock wave that moves toward the airfoil leading edge. A thick shear layer forms downstream of the shock wave. The kinetic energy and shear stresses increase dramatically, reach a maximum when dissipation and diffusion of the turbulence exceed production, and then decrease substantially. The response lime of the turbulence to the changes brought about by the shock-wave passage upstream depends on the shock-wave strength and position in the boundary layer. The cycle completes itself when the shock wave passes the midchord, weakens, and the shear layer collapses. Remarkably good comparisons are found with computations that employ the time-dependent Reynolds averaged form of the Navier-Stokes equations using an algebraic eddy viscosity model, developed for steady flows.

Marvin, J. G.↗

Turbulence Modeling for Unsteady Transonic Flows

Conditionally sampled, ensemble-averaged velocity measurements, made with a laser velocimeter, were taken in the flowfield over the rear half of an 18% thick circular arc airfoil at zero incidence tested at M = 0.76 and of a Reynolds number based on chord of 11 x 10(exp 6). Data for one cycle of periodic unsteady flow having a reduced frequency bar-f of 0.49 are analyzed. A series of compression waves, which develop in the early stages of the cycle, strengthen and coalesce into a strong shock wave that moves toward the airfoil leading edge. A thick shear layer forms downstream of the shock wave. The kinetic energy and shear stresses increase dramatically, reach a maximum when dissipation and diffusion of the turbulence exceed production, and then decrease substantially. The response time of the turbulence to the changes brought about by the shock-wave passage upstream depends on the shock-wave strength and position in the boundary layer. The cycle completes itself when the shock wave passes the midchord, weakens, and the shear layer collapses. Remarkably good comparisons are found with computations that employ the time-dependent Reynolds averaged form of the Navier-Stokes equations using an algebraic eddy viscosity model, developed for steady flows.

Marvin, J. G.↗

Global Warming Estimation from MSU: Correction for Drift and Calibration Errors

Microwave Sounding Unit (MSU) radiometer observations in Ch 2 (53.74 GHz), made in the nadir direction from sequential, sun-synchronous, polar-orbiting NOAA morning satellites (NOAA 6, 10 and 12 that have approximately 7am/7pm orbital geometry) and. afternoon satellites (NOAA 7, 9, 11 and 14 that have approximately 2am/2pm orbital geometry) are analyzed in this study to derive global temperature trend from 1980 to 1998. In order to remove the discontinuities between the data of the successive satellites and to get a continuous time series, first we have used shortest possible time record of each satellite. In this way we get a preliminary estimate of the global temperature trend of 0.21 K/decade. However, this estimate is affected by systematic time-dependent errors. One such error is the instrument calibration error eo. This error can be inferred whenever there are overlapping measurements made by two satellites over an extended period of time. From the available successive satellite data we have taken the longest possible time record of each satellite to form the time series during the period 1980 to 1998 to this error eo. We find eo can decrease the global temperature trend by approximately 0.07 K/decade. In addition there are systematic time dependent errors ed and ec present in the data that are introduced by the drift in the satellite orbital geometry. ed arises from the diurnal cycle in temperature and ec is the drift related change in the calibration of the MSU. In order to analyze the nature of these drift related errors the multi-satellite Ch 2 data set is partitioned into am and pm subsets to create two independent time series. The error ed can be assessed in the am and pm data of Ch 2 on land and can be eliminated. Observation made in the MSU Ch 1 (50.3 GHz) support this approach. The error ec is obvious only in the difference between the pm and am observations of Ch 2 over the ocean. We have followed two different paths to assess the impact of the error ec on the global temperature trend. In one path the entire error ec is placed in the am data while in the other it is placed in the pm data. Global temperature trend is increased or decreased by approximately 0.03 K/decade depending upon this placement. Taking into account all random errors and systematic errors our analysis of MSU observations leads us to conclude that a conservative estimate of the global warming is 0. 11 (+/-) 0.04 K/decade during 1980 to 1998.

Prabhakara, C.↗

Global Warming Estimation from MSU: Correction for Drift and Calibration Errors

Microwave Sounding Unit (MSU) radiometer observations in Ch 2 (53.74 GHz), made in the nadir direction from sequential, sun-synchronous, polar-orbiting NOAA morning satellites (NOAA 6, 10 and 12 that have about 7am/7pm orbital geometry) and afternoon satellites (NOAA 7, 9, 11 and 14 that have about 2am/2pm orbital geometry) are analyzed in this study to derive global temperature trend from 1980 to 1998. In order to remove the discontinuities between the data of the successive satellites and to get a continuous time series, first we have used shortest possible time record of each satellite. In this way we get a preliminary estimate of the global temperature trend of 0.21 K/decade. However, this estimate is affected by systematic time-dependent errors. One such error is the instrument calibration error. This error can be inferred whenever there are overlapping measurements made by two satellites over an extended period of time. From the available successive satellite data we have taken the longest possible time record of each satellite to form the time series during the period 1980 to 1998 to this error. We find we can decrease the global temperature trend by about 0.07 K/decade. In addition there are systematic time dependent errors present in the data that are introduced by the drift in the satellite orbital geometry arises from the diurnal cycle in temperature which is the drift related change in the calibration of the MSU. In order to analyze the nature of these drift related errors the multi-satellite Ch 2 data set is partitioned into am and pm subsets to create two independent time series. The error can be assessed in the am and pm data of Ch 2 on land and can be eliminated. Observations made in the MSU Ch 1 (50.3 GHz) support this approach. The error is obvious only in the difference between the pm and am observations of Ch 2 over the ocean. We have followed two different paths to assess the impact of the errors on the global temperature trend. In one path the entire error is placed in the am data while in the other it is placed in the pm data. Global temperature trend is increased or decreased by about 0.03 K/decade depending upon this placement. Taking into account all random errors and systematic errors our analysis of MSU observations leads us to conclude that a conservative estimate of the global warming is 0. 11 (+-) 0.04 K/decade during 1980 to 1998.

Prabhakara, C.↗

Thermal Field Imaging Using Ultrasound

It is often desirable to be able to determine the temperature field in the interiors of opaque fluids forced into convection by externally imposed temperature gradients. To measure the temperature at a point in an opaque fluid in the usual fashion requires insertion of a probe, and to determine the full field therefore requires either the ability to move this probe or the introduction of multiple probes. Neither of these solutions is particularly satisfactory, although they can lead to quite accurate measurements. As an alternative we have investigated the use of ultrasound as a relatively non-intrusive probe of the temperature field in convecting opaque fluids. The temperature dependence of the sound velocity can be sufficiently great to permit a determination of the temperature from timing the traversal of an ultrasound pulse across a chamber. In this paper we will present our results on convecting flows of transparent and opaque fluids. Our experimental cells consist of relatively narrow rectangular cavities made of thermally insulating materials on the sides, and metal top and bottom plates. The ultrasound transducer is powered by a pulser/receiver, the signal output of which goes to a very high speed signal averager. The average of several hundred to several thousand signals is then sent to a computer for storage and analysis. The experimental procedure is to establish a convective flow by imposing a vertical temperature gradient on the chamber, and then to measure, at several regularly spaced locations, the transit time for an ultrasound pulse to traverse the chamber horizontally (parallel to the convecting rolls) and return to the transducer. The transit time is related to the temperature of the fluid through which the sound pulse travels. Knowing the relationship between transit time and temperature (determined in a separate experiment), we can extract the average temperature across the chamber at that location. By changing the location of the transducer it is then possible to find the average temperature at different locations along the chamber, thereby determining the temperature profile along the system. (In the future we will construct an array of transducers. This will give us the capability to determine the temperature profile much more rapidly than at present, an important consideration if time-dependent phenomena are to be studied.) To validate our procedure we introduced encapsulated liquid crystal particles into glycerol. The liquid crystal particles' color varies depending on the temperature of the fluid. A photograph of the fluid through transparent sidewalls therefore gives a picture of the temperature field of the convecting fluid, independent of our ultrasound imaging. A representative result is shown in the Figure 1, which reveals a very satisfying correspondence between the two techniques. Therefore we have a great deal of confidence that the ultrasound imaging approach is indeed measuring the actual temperature profile of the fluid. The technique has also been applied to convecting liquid metal flows, and representative data will be presented from those experiments as well.

Andereck, D.↗

ATIC Flight Data Processing

The first flight of the Advanced Thin Ionization Calorimeter (ATIC) experiment from McMurdo, Antarctica lasted for 16 days, starting in December, 2000. The ATIC instrument consists of a fully active 320-crystal, 960-channel Bismuth Germanate (BGO) calorimeter, 202 scintillator strips in 3 hodoscopes interleaved with a graphite target, and a 4480-pixel silicon matrix charge detector. We have developed an Object Oriented data processing package based on ROOT. In this paper, we will describe the data processing scheme used in handling the accumulated 45 GB of flight data. We will also discuss trigger issues by comparing the measured energy-dependent trigger efficiency with its simulation and calibration issues by considering the time-dependence of housekeeping information, etc.

Ahn, H. S.↗

On the Dynamics of Implicit Linear Multistep Methods

Some new guidelines on the usage of implicit linear multistep methods (LMMs) as time-dependent approaches for obtaining steady-state numerical solutions in computational fluid dynamics (CFD) are explored. The commonly used implicit LMMs in CFD belong to the class of superstable time discretizations. It can be shown that the nonlinear asymptotic behavior in terms of bifurcation diagrams and basins of attractions of these schemes can provide an improved range of initial data and time step over the linearized stability limit.

Yee, H. C.↗

Incorporation of Solar Noble Gases from a Nebula-Derived Atmosphere During Magma Ocean Cooling

The presence of solar noble gases in the deep interior of the Earth is inferred from the Ne isotopic compositions of MORB (Mid-ocean Ridge Basalts) and OIB (Oceanic Island Basalt); Ar data may also consistent with a solar component in the deep mantle. Models of the transport and distribution of noble gases in the earth's mantle allow for the presence of solar Ar/Ne and Xe/Ne ratios and permit the calculation of lower mantle noble gas concentrations. These mantle data and models also indicate that the Earth suffered early (0.7 to 2 x 10(exp 8) yr) and large (greater than 99 percent) losses of noble gases from the interior, a result previously concluded for atmospheric Xe. We have pursued the suggestion that solar noble gases were incorporated in the forming Earth from a massive, nebula-derived atmosphere which promoted large-scale melting, so that gases from this atmosphere dissolved in the magma ocean and were mixed downward. Models of a primitive atmosphere captured from the solar nebula and supported by accretion luminosity indicate that pressures at the Earth's surface were adequate (and largely more than the required 100 Atm) to dissolve sufficient gases. We have calculated the coupled evolution of the magma ocean and the overlying atmosphere under conditions corresponding to the cessation (or severe attenuation) of the sustaining accretion luminosity, prior to the complete removal of the solar nebula. Such a condition was likely to obtain, for instance, when most of the unaccumulated mass resided in large bodies which were only sporadically accreted. The luminosity supporting the atmosphere is then that provided by the cooling Earth, consideration of which sets a lower limit to the time required to solidify the mantle and terminate the incorporation of atmospheric gases within it. In our initial calculations, we have fixed the nebula temperature at To = 300K, a value likely to be appropriate for nebular temperatures at lAU in the early planet-building epoch. We treated the background (nebula) pressure as an adjustable, time-dependent parameter. Additional information is contained within the original extended abstract.

Woolum, D. S.↗

Structured Course Objects in a Digital Library

We are developing an Undergraduate Digital Library Framework (UDLF) that will support creation/archiving of courses and reuse of existing course material to evolve courses. UDLF supports the publication of course materials for later instantiation for a specific offering and allows the addition of time-dependent and student-specific information and structures. Instructors and, depending on permissions, students can access the general course materials or the materials for a specific offering. We are building a reference implementation based on NCSTRL+, a digital library derived from NCSTRL. Digital objects in NCSTRL+ are called buckets, self-contained entities that carry their own methods for access and display. Current bucket implementations have a two level structure of packages and elements. This is not a rich enough structure for course objects in UDLF. Typically, courses can only be modeled as a multilevel hierarchy and among different courses, both the syntax and semantics of terms may vary. Therefore, we need a mechanism to define, within a particular library, course models, their constituent objects, and the associated semantics in a flexible, extensible way. In this paper, we describe our approach to define and implement these multilayered course objects. We use XML technology to emulate complex data structures within the NCSTRL+ buckets. We have developed authoring and browsing tools to manipulate these course objects. In our current implementation a user downloading an XML based course bucket also downloads the XML-aware tools: an applet that enables the user to edit or browse the bucket. We claim that XML provides an effective means to represent multi-level structure of a course bucket.

Maly, K.↗

Issues on Reproducibility/Reliability of Magnetic NDE Methods

One of the critical elements related to the practicality of any NDE technique is its reproducibility under nominally the same inspection conditions. The results of certain test methodologies, however, are not always repeatable and understanding the origin of the irreproducibility is often as critical as obtaining reproducible results. One example is the characterization of residual stress in structural ferromagnets using the magnetoacoustic (MAC) method. Although it has not been widely publicized, the test results of this method are known to be time-dependent. Two distinct types of time dependencies have been observed during testing. The first type has a clearly definable relaxation time, while no such trend has been observed for the second. The purpose of the present study is to systematically investigate the time dependence of the second type, to find out the range and, if possible, the origin of the variation in the test results. For this, MAC curves were obtained under various stress levels and the tests were repeated over time. Particular attention was given to whether noise in the measuring device or a change in the laboratory environment could have been a contributing factor. The steel samples used for the study were cut from C- and U-class railroad wheels. The MAC behavior of these samples was reported previously. Each steel sample was first machined to be a cylindrical rod of 3.175 cm (1.25 in) in diameter and 26.67 cm (10.5 in) in length. The center portion of the samples were further machined to form a pair of flat and parallel surfaces for the launch and reflection of the ultrasonic pulses. Data acquisition involved two major elements; magnetic and acoustic measurements. Throughout the experiment the net magnetic induction, B, was measured by integrating the induction pickup coil output using an integrating fluxmeter. The acoustic measurements employed the phase-locked technique which will be described in the following.

Namkung, M.↗

Cameleon calcium indicator reports cytoplasmic calcium dynamics in Arabidopsis guard cells

Cytoplasmic free calcium ([Ca2+]cyt) acts as a stimulus-induced second messenger in plant cells and multiple signal transduction pathways regulate [Ca2+]cyt in stomatal guard cells. Measuring [Ca2+]cyt in guard cells has previously required loading of calcium-sensitive dyes using invasive and technically difficult micro-injection techniques. To circumvent these problems, we have constitutively expressed the pH-independent, green fluorescent protein-based calcium indicator yellow cameleon 2.1 in Arabidopsis thaliana (Miyawaki et al. 1999; Proc. Natl. Acad. Sci. USA 96, 2135-2140). This yellow cameleon calcium indicator was expressed in guard cells and accumulated predominantly in the cytoplasm. Fluorescence ratio imaging of yellow cameleon 2.1 allowed time-dependent measurements of [Ca2+]cyt in Arabidopsis guard cells. Application of extracellular calcium or the hormone abscisic acid (ABA) induced repetitive [Ca2+]cyt transients in guard cells. [Ca2+]cyt changes could be semi-quantitatively determined following correction of the calibration procedure for chloroplast autofluorescence. Extracellular calcium induced repetitive [Ca2+]cyt transients with peak values of up to approximately 1.5 microM, whereas ABA-induced [Ca2+]cyt transients had peak values up to approximately 0.6 microM. These values are similar to stimulus-induced [Ca2+]cyt changes previously reported in plant cells using ratiometric dyes or aequorin. In some guard cells perfused with low extracellular KCl concentrations, spontaneous calcium transients were observed. As yellow cameleon 2.1 was expressed in all guard cells, [Ca2+]cyt was measured independently in the two guard cells of single stomates for the first time. ABA-induced, calcium-induced or spontaneous [Ca2+]cyt increases were not necessarily synchronized in the two guard cells. Overall, these data demonstrate that that GFP-based cameleon calcium indicators are suitable to measure [Ca2+]cyt changes in guard cells and enable the pattern of [Ca2+]cyt dynamics to be measured with a high level of reproducibility in Arabidopsis cells. This technical advance in combination with cell biological and molecular genetic approaches will become an invaluable tool in the dissection of plant cell signal transduction pathways.

Non-NASA Center↗

Altering the cellular mechanical force balance results in integrated changes in cell, cytoskeletal and nuclear shape

Studies were carried out with capillary endothelial cells cultured on fibronectin (FN)-coated dishes in order to analyze the mechanism of cell and nuclear shape control by extracellular matrix (ECM). To examine the role of the cytoskeleton in shape determination independent of changes in transmembrane osmotic pressure, membranes of adherent cells were permeabilized with saponin (25 micrograms/ml) using a buffer that maintains the functional integrity of contractile microfilaments. Real-time videomicroscopic studies revealed that addition of 250 microM ATP resulted in time-dependent retraction and rounding of permeabilized cells and nuclei in a manner similar to that observed in intact living cells following detachment using trypsin-EDTA. Computerized image analysis confirmed that permeabilized cells remained essentially rigid in the absence of ATP and that retraction was stimulated in a dose-dependent manner as the concentration of ATP was raised from 10 to 250 microM. Maximal rounding occurred by 30 min with projected cell and nuclear areas being reduced by 69 and 41%, respectively. ATP-induced rounding was also accompanied by a redistribution of microfilaments resulting in formation of a dense net of F-actin surrounding retracted nuclei. Importantly, ATP-stimulated changes in cell, cytoskeletal, and nuclear form were prevented in permeabilized cells using a synthetic myosin peptide (IRICRKG) that has been previously shown to inhibit actomyosin filament sliding in muscle. In contrast, both the rate and extent of cell and nuclear rounding were increased in permeabilized cells exposed to ATP when the soluble FN peptide, GRGDSP, was used to dislodge immobilized FN from cell surface integrin receptors.(ABSTRACT TRUNCATED AT 250 WORDS).

Non-NASA Center↗

CARES/Life Ceramics Durability Evaluation Software Enhanced for Cyclic Fatigue

The CARES/Life computer program predicts the probability of a monolithic ceramic component's failure as a function of time in service. The program has many features and options for materials evaluation and component design. It couples commercial finite element programs--which resolve a component's temperature and stress distribution--to reliability evaluation and fracture mechanics routines for modeling strength-limiting defects. The capability, flexibility, and uniqueness of CARES/Life have attracted many users representing a broad range of interests and has resulted in numerous awards for technological achievements and technology transfer. Recent work with CARES/Life was directed at enhancing the program s capabilities with regards to cyclic fatigue. Only in the last few years have ceramics been recognized to be susceptible to enhanced degradation from cyclic loading. To account for cyclic loads, researchers at the NASA Lewis Research Center developed a crack growth model that combines the Power Law (time-dependent) and the Walker Law (cycle-dependent) crack growth models. This combined model has the characteristics of Power Law behavior (decreased damage) at high R ratios (minimum load/maximum load) and of Walker law behavior (increased damage) at low R ratios. In addition, a parameter estimation methodology for constant-amplitude, steady-state cyclic fatigue experiments was developed using nonlinear least squares and a modified Levenberg-Marquardt algorithm. This methodology is used to give best estimates of parameter values from cyclic fatigue specimen rupture data (usually tensile or flexure bar specimens) for a relatively small number of specimens. Methodology to account for runout data (unfailed specimens over the duration of the experiment) was also included.

Nemeth, Noel N.↗

Constituent-Based Life Models Being Developed for SiC/SiC Composites

For the successful utilization of ceramic matrix composites (CMC) as hot-section components in advanced aeropropulsion engines, the CMC constituents will need to be tailored and optimized to meet all the critical property demands of each component. Under the High-Speed Research (HSR) and Advanced High-Temperature Engine Materials Technology (HITEMP) Programs, the NASA Glenn Research Center at Lewis Field initiated research to develop mechanistic models for key CMC thermostructural properties. These models would describe the effects of different constituent factors (composition, geometry, and volume fraction) and of potential application conditions (stress, time, temperature, and environment) on these properties. Particular focus was placed on both analytical and numerical modeling of state-of-the-art SiC/SiC composites where the primary load-bearing constituents are stoichiometric SiC fibers in a complex multiphase SiC matrix produced by chemical vapor infiltration and melt infiltration. Recent studies have resulted in computer-based numerical models for the elastic modulus, thermal expansion, and thermal conductivity properties of the SiC/SiC system. Additional studies have generated analytical and empirical models for the time dependence of composite rupture strength at temperatures above 2200 F (1200 C), where CMC's have an important thermostructural advantage over current nickel-based superalloys. These life models utilize thermal activation theory and fiber stress-rupture results measured at Glenn to generate Larson-Miller (L M) plots of fiber rupture strength versus q, a single time- and temperature-dependent parameter. Assuming a worse case in which the SiC matrix is cracked, rupture is then controlled by the time-dependent fracture characteristics of the fiber bundles bridging the matrix cracks. With this as the controlling mechanism, one can then use simple composite theory and the fiber L M plots to predict CMC rupture strength versus the q parameter. The dashed line shows the predicted rupture strength of a SiC/SiC composite that is reinforced by a state-of-the-art stoichiometric SiC fiber. For the q parameter, time is in hours and temperature in degrees kelvin. To generate these predictions, a two-dimensional 0 /90 composite with approximately 16 percent fiber in the applied stress direction and an air test environment were assumed. As such, it is possible to compare the model predictions against limited stress-rupture data for this CMC as shown by the data points. The good agreement confirms the rupture model at least for the selected CMC and test conditions. Thus for this particular SiC/SiC composite, one can estimate a 1000-hr rupture strength of approximately 12 ksi at 2400 F (T = 1588 K and q = 39 700 K). At lower CMC application stresses, the SiC matrix is typically uncracked, so both the fiber and matrix constituents share the composite load. In this case, CMC rupture is controlled by the constituent with the longest rupture time based on the creep rate of the composite. Measured Monkman-Grant plots of rupture time versus creep rate for the two SiC constituents have been used to develop CMC life models for this important application condition. NASA and DOD are currently using this information to establish application and material goals for more advanced CMC's that can be used at even higher temperatures.

DiCarlo, James A.↗

Monitoring Land Use and Degradation Using Satellite and Airborne Data

Data on the Manix Basin Area of the Mojave Desert have been collected by AVIRIS, AIRSAR, and Landsat TM. The purpose of this study is to use a fusion of a time series of satellite data with airborne data to provide a context for the airborne data. The satellite data time series will additionally help to validate the observation and analysis of time-dependent processes observed in the single AVIRIS image of fields abandoned for different periods of time.

land↗

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗