Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient flow dynamics”

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

Modeling of Thermal Performance of Multiphase Nuclear Fuel Cell Under Variable Gravity Conditions

A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.

Ding, Z.↗

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗

Constructing Synoptic Maps of Stratospheric Column Ozone from HALOE, SAGE and Balloonsonde Data Using Potential Vorticity Isentropic Coordinate Transformations

In this study we utilize potential vorticity - isentropic (PVI) coordinate transformations as a means of combining ozone data from different sources to construct daily, synthetic three-dimensional ozone fields. This methodology has been used successfully to reconstruct ozone maps in particular regions from aircraft data over the period of the aircraft campaign. We expand this method to create high-resolution daily global maps of profile ozone data, particularly in the lower stratosphere, where high-resolution ozone data are sparse. Ozone climatologies in PVI-space are constructed from satellite-based SAGE II and UARS/HALOE data, both of which-use solar occultation techniques to make high vertical resolution ozone profile measurements, but with low spatial resolution. A climatology from ground-based balloonsonde data is also created. The climatologies are used to establish the relationship between ozone and dynamical variability, which is defined by the potential vorticity (in the form of equivalent latitude) and potential temperature fields. Once a PVI climatology has been created from data taken by one or more instruments, high-resolution daily profile ozone field estimates are constructed based solely on the PVI fields, which are available on a daily basis from NCEP analysis. These profile ozone maps could be used for a variety of applications, including use in conjunction with total ozone maps to create a daily tropospheric ozone product, as input to forecast models, or as a tool for validating independent ozone measurements when correlative data are not available. This technique is limited to regions where the ozone is a long-term tracer and the flow is adiabatic. We evaluate the internal consistency of the technique by transforming the ozone back to physical space and comparing to the original profiles. Biases in the long-term average of the differences are used to identify regions where the technique is consistently introducing errors. Initial results show the technique is useful in the lower stratosphere at most latitudes throughout the year,and in the winter hemisphere in the middle stratosphere. The results are problematic in the summer hemisphere middle stratosphere due to increased ozone photochemistry and weak PV gradients. Alternate techniques in these regions will be discussed. An additional limitation is the quality and resolution of the meteorological data.

Hollandsworth, Stacey M.↗

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E × B drift to cause a nondiffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine f L , the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, f L scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should, in principle, relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear $\hat{s}$ was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of f L with the electron collisionality. Still, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Report on Predictive Analyses of PRD as Function of Anomalies

This report summarizes the work completed in FY-2024 to analyze the power reactivity decrement (PRD) concepts of the ARC-100 core. The PRD has been traditionally defined by the reactivity change from a hot zero power (HZP) to a particular power. Consequently, the PRD accounts for the core reactivity changes due to increase coolant temperature gradient axially and radially across the core, and increased fuel temperature. The coolant temperature gradient leads to sodium and structure density changes, to radial core expansion from assembly flowering and bowing (due to axial and radial temperature gradients within the assemblies), and to control rod driveline thermal expansion. The fuel temperature increase associated with coolant temperature and power increases leads to Doppler effect and axial thermal expansion. In this work, the normal operating Hot Full Power (HFP) state is the only power level of interest, so analyses focus on the PRD calculated from HZP to HFP. The PRD has been used to assess the reactor safety features asymptotically in unprotected accident scenarios, including the loss of heat sink (LOHS), loss of flow (LOF), and transient overpower (TOP) without scram. The PRD concept relies on the “global” reactivity coefficients A, B, and C that are estimated based on “individual” reactivity effects (Doppler, sodium density, etc.). The objectives of this work are: 1) to improve and verify the methodology used to calculate the ABC coefficients used in the PRD, 2) to assess if the PRD can be used to reliably identify abnormal events. This report fulfills the FY-2024 scope of WBS#1.15.8.3 activity, “ANL0120 – Predictive Analyses of PRD as function of deformation”. The PRD concept is described in Sections 2. Additional effort in refining the methodology for core bowing modeling is performed in Sections 3. Then, two verification exercises are proposed in Section 4 to benchmark these coefficients based on direct neutronic-only calculations and on dynamic core transient simulations. Finally, the PRD approach is assessed for the detection of several unexpected events, such as primary flow perturbation or improper fuel loading, in Section 5.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Investigation of Scale-Dependent Groundwater/Surface-water Exchange in Rivers by Gradient Self-Potential Logging: Numerical Modeling and Field Experiments

Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.

Engineering↗

A new model of resonance in the winter stratosphere

It is generally accepted that the planetary waves observed in the winter stratosphere are primarily a response to dynamical forcing from the troposphere. Nevertheless, the mechanism by which wave amplitudes sometimes become larger remains uncertain. It is possible that anomalously large waves in the stratosphere might simply be the result of anomalously large tropospheric forcing. However, it has also been suggested that they are a response to a stratospheric-tropospheric cavity being in a near-resonant configuration. It has been suggested that nonlinear self-tuning effects could play an important role in the behavior of such a cavity. Self-tuning may occur when a system starts to one side of resonance, such that the mean-state change induced by growing waves brings the system closer to resonance. A new model of the stratospheric cavity is introduced and is then used to re-examine the possibility of wave growth in the real atmosphere and in atmospheric models due to self-tuning effects. The new model is based on the picture of the winter-time stratosphere which has been revealed by the observations of Ertel's potential vorticity, Q. Isentropic maps of Q show two rather distinct regions, the first containing the circumpolar vortex, where gradients of Q are large and Rossby waves may propagate easily. Surrounding this is a second, low-latitude region where the gradients are generally weak and where, because the Eulerian-mean flow is comparable with their phase speed, Rossby waves must be continually breaking. As the waves are observed to grow the relative sizes of these two regions change in time. This leads to the interesting possibility that self-tuning, mainly due to irreversible changes in the size of the polar vortex, is taking place.

Haynes, P. H.↗

Measurements of Local Heat Transfer and Pressure on Six 2-Inch-Diameter Blunt Bodies at a Mach Number of 4.95 and at Reynolds Numbers Per Foot up to 81 x 10(exp 6)

Measurements of the local heat transfer and pressure distribution have been made on six 2-inch-diameter, blunt, axially symmetric bodies in the Langley gas dynamics laboratory at a Mach number of 4.95 and at Reynolds numbers per foot up to 81 x 10(exp 6). During the investigation laminar flow was observed over a hemispherical-nosed body having a surface finish from 10 to 20 microinches at the highest test Reynolds number per foot (for this configuration) of 77.4 x 10(exp 6). Though it was repeatedly possible to measure completely laminar flow at this Reynolds number for the hemisphere, it was not possible to observe completely laminar flow on the flat-nosed body for similar conditions. The significance of this phenomenon is obscured by the observation that the effects of particle impacts on the surface in causing roughness were more pronounced on the flat-nosed body. For engineering purposes, a method developed by M. Richard Dennison while employed by Lockheed Aircraft Corporation appears to be a reasonable procedure for estimating turbulent heat transfer provided transition occurs at a forward location on the body. For rearward-transition locations, the method is much poorer for the hemispherical nose than for the flat nose. The pressures measured on the hemisphere agreed very well with those of the modified Newtonian theory, whereas the pressures on all other bodies, except on the flat-nosed body, were bracketed by modified Newtonian theory both with and without centrifugal forces. For the hemisphere, the stagnation-point velocity gradient agreed very well with Newtonian theory. The stagnation-point velocity gradient for the flat- nosed model was 0.31 of the value for the hemispherical-nosed model. If a Newtonian type of flow is assumed, the ratio 0.31 will be independent of Much number and real-gas effects.

Cooper, Morton↗

Data Base for CFD Validation

The flow behind the shock wave formed around objects which fly at hypervelocity behaves differently from that of a perfect gas. Molecules become vibrationally excited, dissociated, and ionized. The hot gas may emit or absorb radiation. When the atoms produced by dissociation reach the wall surface, chemical reactions, including recombination, may occur. The thermochemical phenomena of vibration, dissociation, ionization, surface chemical reaction, and radiation are referred to commonly as high-temperature real-gas phenomena. The phenomena cause changes in the dynamic behavior of the flow and the surface pressure and heat transfer distribution around the object. The character of a real gas is described by the internal degrees of freedom and state of constituent molecules; nitrogen and oxygen for air. The internal energy states, rotation, vibration and electronic, of the molecules are excited and, in the limit, the molecular bonds are exceeded and the gas dissociated into atomic and, possibly, ionic constituents. The process of energy transfer causing excitation, dissociation and recombination is a rate process controlled by particle collisions. Binary, two-body, collisions are sufficient to cause internal excitation, dissociation and ionization while three-body collisions are required to recombine the particles into molecular constituents. If the rates of energy transfer are fast with respect to the local fluid dynamic time scale the gas is in, or nearly in, equilibrium. If the energy transfer rates are very slow the gas can be described as frozen. In all other instances, wherein any of the energy exchange rates are comparable to the local fluid time scale, the gas will be thermally or chemically reacting and out of equilibrium. Real gas thermochemical nonequilibrium processes are important in the determination of aerodynamic heating; both convective (including wall catalytic effects) and radiative heating. To illustrate this we consider the hypervelocity flow over a bluff body typical of an atmospheric entry vehicle or an aerospace transfer vehicle (ASTV). The qualitative aspects of a hypersonic flow field over a bluff body are discussed in two parts, forebody and afterbody, with attention to which particular physical effects must be included in an analysis. This will indicate what type of numerical modeling will be adequate in each region of the flow. A bluff forebody flow field is dominated by the presence of the strong bow shock wave and the consequent heating, and chemical reaction of the gas. At high altitude hypersonic flight conditions the thermal excitation and chemical reaction of the gas occur slowly enough that a significant portion of the flow field is in a state of thermochemical nonequilibrium. A second important effect is the presence of the thick boundary layer along the forebody surface. In this region there are large thermal and chemical species gradients due to the interaction of the gas with the wall. Also at high altitudes the shock wave and the boundary layer may become so thick that they merge; in this case the entire shock layer is dominated by viscous effects.

Deiwert, George S.↗

Diffusivity Measurements Made Instant and Easy

A compact common path interferometer (CPI) system has been developed to measure the diffusivity of liquid pairs. The CPI is an optical technique that can be used to measure changes in the gradient of the refraction index of transparent materials. It uses a shearing interferometer that shares the same optical path from a laser light source to the final imaging plane. The molecular diffusion coefficient of liquids can be determined from the physical relations between changes in the optical path length and liquid phase properties. When the data obtained by using the CPI have been compared with similar results from other techniques, the instrument has been demonstrated to be far superior to other instruments for measuring the diffusivity of miscible liquids while staying very compact and robust (ref. 1). Because of its compactness and ease of use, the CPI has been adopted for use in studies of interface dynamics as well as other diffusion-controlled process applications (ref. 2). This progress will permit experiments in microgravity that can quantitatively answer basic science questions about mass and thermal diffusion and their effect in transport processes. This instrument is a spinoff of a diagnostic development for microgravity fluid physics experiments at the NASA Glenn Research Center that has used optics and electronics existing in the fluid physics laboratory for feasibility studies. Optical diagnostic techniques have become an integral part of many areas of measurement applications in industrial and research laboratories. Many types of interferometers and their phase-shifted versions have been used as instruments for measuring optical wave fronts for lens testing and combustion and fluid flow diagnostics. One of these, the point diffraction interferometer, is considered to be robust (see, for example, ref. 3) because it has a common-path design. The point diffraction interferometer is difficult to align and has a limited measurement range for liquid-phase applications. Interferometry and schlieren techniques have been widely used for many years for gas-flow-visualization. On the other hand, the CPI, which is basically a Wollaston prism Polarizer in combination with an analyser (crystalline quartz retardation plate) is used instead of a point diffraction interferometer unit. The advantage of using the CPI over other optical techniques is that it can make quantitative measurements in liquids with the large index of refraction variations that often occur in interface dynamics studies. This can be simply accomplished by using different prism polarizers for each particular experimental condition.

Rashidnia, Nasser↗

Hydrological and Thermal Dynamics of a Supra‐Permafrost Subterranean Estuary

Subterranean estuaries (STEs), where groundwater interacts with seawater, influence surface and subsurface coastal ecology and biogeochemistry. In Arctic-STEs overlying permafrost, groundwater flow and heat transport determine the fate of organic matter. Yet, direct observations of groundwater flow and heat and solute transport processes in Arctic STEs remain limited. This study characterized groundwater flow paths and fluxes and heat transport within an Arctic-STE along Alaska's Beaufort Sea coast during thawing, summer, and freeze-up. Intertidal seabed temperature-depth profiles collected along a 10-m transect captured the active groundwater flow period, from thaw and flow onset in mid-June to freeze-up in late-September. During this period, aquifer geometry evolved non-uniformly due to spatially varying thaw rates across the STE (mean (m) thaw depths–beach: 0.25 to 0.55–0.6 m on 20 June, 25 July, 1 October; seabed: 0.6–0.9 m from 25 July to 1 October). Groundwater and surface water levels, salinity, and subsurface temperature profiles measured over tidal time scales were interpreted alongside groundwater flow-heat transport numerical simulations. Fresh groundwater discharge was sporadic during thawing (m: 0.32 m 3 /day/m), abundant in summer (m: 0.45 m 3 /day/m), and was largely absent during freeze-up. During freeze-up, groundwater flow was driven exclusively by seawater recirculation via tidal pumping (from thawing to summer to freeze-up: 0.00025–0.15–0.5 m3/day/m) and convection. Heat advection dominated near aquatic interfaces (shaping intertidal ice), and conduction controlled vertical temperature gradients in low-flow and unsaturated sediments. These findings will help predict how prolonged summers will alter Arctic-STE cryo-hydrology and biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗

EVALUATION OF OPERATION TEMPERATURES UNDER NATURAL CONVECTION HELIUM FLOW IN A CONFINED CAVITY

The Spallation Neutron Source (SNS) is a high-power accelerator-based pulsed neutron source led by Oak Ridge National Laboratory (ORNL) to achieve high fluxes of neutrons for scientific experiments. Active and passive cooling of the systems and parts forming the SNS have been considered to warrant the safe operation of the facility. The diverse cooling systems make use of conjugated heat transfer mechanisms to provide a stable operation temperature for all components in the machine. Thermal power deposited into stainless-steel piping lines due to particle radiation may reach values of up to 1.2 W/cc in the regions located closer to the center of the lower IRP. These energy deposition levels, in not actively cooled components, such as the transfer-line-outer-vacuum-layer may increase the temperature of the component beyond design requirement limitations. The evaluation of the operation temperatures for the former components relies in the assumption that a low-pressure helium atmosphere provides enough heat removal capacity based on natural convection phenomena. In this work, the evaluation of steady state temperatures in components such as the CMS transfer lines has been evaluated using computational fluid dynamics (CFD), analytical correlations and experimental measurements. The companion experiments were conducted in a closed helium system at pressures varying from 1.1 to 1.5 bar. A copper rod was affixed horizontally between viewing windows and heated at constant power, and measurements were made of both the rod temperature and ambient temperature via a system fiberoptic distributed temperature sensors and RTDs. It was found that the measured heat transfer coefficients agree well with the predictions of Churchill and Chu correlations across the range of cases considered. Additionally, the ambient helium volume above the rod was imaged via background oriented schlieren (BOS), and these data was used to determine the line-averaged density gradients in this region. These gradients were compared to simulation data to validate the predictions of natural convection simulations.

Dominguez-Ontiveros, Elvis [ORNL] (ORCID:000000018↗

An "Inefficient Fin" Non-Dimensional Parameter to Measure Gas Temperatures Efficiently

A gas containment vessel that is not in thermal equilibrium with the bulk gas can affect its temperature measurement. The physical nature of many gas dynamics experiments often makes the accurate measurement of temperature a challenge. The environment itself typically requires that the thermocouple be sheathed, both to protect the wires and hot junction of the instrument from their environment, and to provide a smooth, rigid surface for pressure sealing of the enclosure. However, that enclosure may also be much colder than the gas to be sensed, or vice-versa. Either way, the effect of such gradients is to potentially skew the temperature measurements themselves, since heat may then be conducted by the instrument. Thermocouple designers traditionally address this problem by insulating the sheath from the thermocouple leads and hot junction as much as possible. The thermocouple leads are typically packed in a ceramic powder inside the sheath, protecting them somewhat from temperature gradients along the sheath, but there is no effective mechanism to shield the sheath from the enclosure body itself. Standard practice dictates that thermocouples be used in installations that do not present large thermal gradients along the probe. If this conduction dominates heat transfer near the tip of the probe, then temperature measurements may be expected to be skewed. While the same problem may be experienced in the measurement of temperature at various points within a solid in a gradient, it tends to be aggravated in the measurements of gas temperature, since heat transfer dependent on convection is often less efficient than conduction along the thermocouple. The proposed solution is an inefficient fin thermocouple probe. Conventional wisdom suggests that in many experiments where gas flows through an enclosure (e.g., flow in pipe, manifold, nozzle, etc.), the thermocouple be introduced flush to the surface, so as not to interfere with the flow. In practice, however, many such experiments take place where the flow is already turbulent, so that a protruding thermocouple probe has a negligible effect on the flow characteristics. The key question then becomes just how far into the flow should a thermocouple protrude in order to properly sense the gas temperature at that point. Modeling the thermocouple as an "inefficient fin" directly addresses this question. The appropriate assumptions in this case are: one-dimensional conduction along the fin; steady-state, constant, and homogeneous thermal conductivity; negligible radiation; and a uniform, constant heat transfer coefficient over the probe surface. It is noted that the nature of the ceramic-filled probe makes the key assumption of homogeneous thermal conductivity that much more conservative.

Lemieux, Patrick↗

Evapotranspiration Partitioning of Eucalyptus benthamii and Pinus taeda During Early Stand Development

Increasing demand for bioenergy from intensively managed woody crops raises concerns of increased evapotranspiration and potential decreases in water yield. Pinus taeda (loblolly pine) is currently the most cultivated species in the southeastern USA, the country’s wood basket. However, Eucalyptus species could achieve greater productivity but with unknown ramifications for water budgets. To address the knowledge gap, we determined annual water budgets of loblolly pine and Eucalyptus benthamii (eucalypt) from growing years 3 through 5 in a replicated (n = 3) two-factor design comparing species and groundwater depth. Paired plots were established across a depth-to-groundwater gradient from shallow (~ 2 m) to deep (~ 8 m). Hydrologic budgets were constructed by measuring precipitation, interception, soil evaporation, and transpiration. Eucalypt evapotranspiration and above-ground biomass production for growing years 3 through 5 were on average 25 and 14% greater than pine, respectively; however, evapotranspiration did not differ across groundwater depths. At the end of growing year 5, eucalypt had higher transpiration and evapotranspiration per unit area than pine. Soil evaporation was substantial in young plantations (nearly 500 mm) in growing year 3 but declined as the canopy closed. Partitioning of evapotranspiration components in developing bioenergy plantations was dynamic due to canopy development driven decreases in soil evaporation and increases in transpiration and interception; total evapotranspiration was less variable from year to year. Water use efficiency (WUE, kg biomass/m 3 H 2 O) per unit evapotranspiration was similar between species, but WUE per unit transpiration was higher in pine. As a result, considering total evapotranspiration in young plantations can affect WUE interpretations.

59 BASIC BIOLOGICAL SCIENCES↗

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗