Publications

Cote, Raymond O.; O'Hern, Timothy J.; Torczynski, J.R.

Abstract not provided.

Cote, Raymond O.; O'Hern, Timothy J.; Torczynski, J.R.

Abstract not provided.

Proposed for publication in Measurement Science and Technology.

Torczynski, J.R.; Tortora, Paul R.

Abstract not provided.

Rader, Daniel J.; Trott, Wayne T.; Torczynski, J.R.; Gallis, Michail A.; Castaneda, Jaime N.; Grasser, Thomas W.

A combined experimental/modeling study was conducted to better understand the critical role of gas-surface interactions in rarefied gas flows. An experimental chamber and supporting diagnostics were designed and assembled to allow simultaneous measurements of gas heat flux and inter-plate gas density profiles in an axisymmetric, parallel-plate geometry. Measurements of gas density profiles and heat flux are made under identical conditions, eliminating an important limitation of earlier studies. The use of in situ, electron-beam fluorescence is demonstrated as a means to measure gas density profiles although additional work is required to improve the accuracy of this technique. Heat flux is inferred from temperature-drop measurements using precision thermistors. The system can be operated with a variety of gases (monatomic, diatomic, polyatomic, mixtures) and carefully controlled, well-characterized surfaces of different types (metals, ceramics) and conditions (smooth, rough). The measurements reported here are for 304 stainless steel plates with a standard machined surface coupled with argon, helium, and nitrogen. The resulting heat-flux and gas-density-profile data are analyzed using analytic and computational models to show that a simple Maxwell gas-surface interaction model is adequate to represent all of the observations. Based on this analysis, thermal accommodation coefficients for 304 stainless steel coupled with argon, nitrogen, and helium are determined to be 0.88, 0.80, and 0.38, respectively, with an estimated uncertainty of {+-}0.02.

Journal of Microelectromechanical Systems

Gallis, Michail A.; Torczynski, J.R.

An improved gas-damping model for the out-of-plane motion of a near-substrate microbeam is developed based on the Reynolds equation (RE). A boundary condition for the RE is developed that relates the pressure at the beam edge to the beam motion. The coefficients in this boundary condition are determined from Navier-Stokes slip-jump (NSSJ) simulations for small slip lengths (relative to the gap height) and from direct simulation Monte Carlo (DSMC) molecular gas dynamics simulations for larger slip lengths. This boundary condition significantly improves the accuracy of the RE when the microbeam width is only slightly greater than the gap height between the microbeam and the substrate. The improved RE model is applied to microbeams fabricated using the SUMMiT V process. © 2004 IEEE.

O'Hern, Timothy J.; Torczynski, J.R.; Trujillo, Steven M.

Abstract not provided.

Rader, Daniel J.; Wally, Karl W.; Ellis, C R.; Torczynski, J.R.

Abstract not provided.

O'Hern, Timothy J.; Torczynski, J.R.; Cote, Raymond O.; Castaneda, Jaime N.

An experimental program was conducted to study a proposed approach for oil reintroduction in the Strategic Petroleum Reserve (SPR). The goal was to assess whether useful oil is rendered unusable through formation of a stable oil-brine emulsion during reintroduction of degassed oil into the brine layer in storage caverns. An earlier report (O'Hern et al., 2003) documented the first stage of the program, in which simulant liquids were used to characterize the buoyant plume that is produced when a jet of crude oil is injected downward into brine. This report documents the final two test series. In the first, the plume hydrodynamics experiments were completed using SPR oil, brine, and sludge. In the second, oil reinjection into brine was run for approximately 6 hours, and sampling of oil, sludge, and brine was performed over the next 3 months so that the long-term effects of oil-sludge mixing could be assessed. For both series, the experiment consisted of a large transparent vessel that is a scale model of the proposed oil-injection process at the SPR. For the plume hydrodynamics experiments, an oil layer was floated on top of a brine layer in the first test series and on top of a sludge layer residing above the brine in the second test series. The oil was injected downward through a tube into the brine at a prescribed depth below the oil-brine or sludge-brine interface. Flow rates were determined by scaling to match the ratio of buoyancy to momentum between the experiment and the SPR. Initially, the momentum of the flow produces a downward jet of oil below the tube end. Subsequently, the oil breaks up into droplets due to shear forces, buoyancy dominates the flow, and a plume of oil droplets rises to the interface. The interface was deflected upward by the impinging oil-brine plume. Videos of this flow were recorded for scaled flow rates that bracket the equivalent pumping rates in an SPR cavern during injection of degassed oil. Image-processing analyses were performed to quantify the penetration depth and width of the oil jet. The measured penetration depths were shallow, as predicted by penetration-depth models, in agreement with the assumption that the flow is buoyancy-dominated, rather than momentum-dominated. The turbulent penetration depth model overpredicted the measured values. Both the oil-brine and oil-sludge-brine systems produced plumes with hydrodynamic characteristics similar to the simulant liquids previously examined, except that the penetration depth was 5-10% longer for the crude oil. An unexpected observation was that centimeter-size oil 'bubbles' (thin oil shells completely filled with brine) were produced in large quantities during oil injection. The mixing experiments also used layers of oil, sludge, and brine from the SPR. Oil was injected at a scaled flow rate corresponding to the nominal SPR oil injection rates. Injection was performed for about 6 hours and was stopped when it was evident that brine was being ingested by the oil withdrawal pump. Sampling probes located throughout the oil, sludge, and brine layers were used to withdraw samples before, during, and after the run. The data show that strong mixing caused the water content in the oil layer to increase sharply during oil injection but that the water content in the oil dropped back to less than 0.5% within 16 hours after injection was terminated. On the other hand, the sediment content in the oil indicated that the sludge and oil appeared to be well mixed. The sediment settled slowly but the oil had not returned to the baseline, as-received, sediment values after approximately 2200 hours (3 months). Ash content analysis indicated that the sediment measured during oil analysis was primarily organic.

Physical Review E - Statistical, Nonlinear, and Soft Matter Physics

Gallis, Michail A.; Torczynski, J.R.; Rader, Daniel J.

The molecular velocity distribution of a gas with heat flow was analyzed using Bird's direct simulation Monte Carlo (DSMC) method. Large numbers of computational molecules represented the gas in DSMC. Chapman-Enskog behavior was obtained for inverse-power-law molecules at continuum nonequilibrium conditions. It was shown that the Sonine-polynomial coefficients differ systematically from their continuum values as the local Knudsen number is increased, at noncontinuum nonequilibrium conditions.

Torczynski, J.R.; Torczynski, J.R.; Gallis, Michail A.

A Micro Electro Mechanical System (MEMS) typically consists of micron-scale parts that move through a gas at atmospheric or reduced pressure. In this situation, the gas-molecule mean free path is comparable to the geometric features of the microsystem, so the gas flow is noncontinuum. When mean-free-path effects cannot be neglected, the Boltzmann equation must be used to describe the gas flow. Solution of the Boltzmann equation is difficult even for the simplest case because of its sevenfold dimensionality (one temporal dimension, three spatial dimensions, and three velocity dimensions) and because of the integral nature of the collision term. The Direct Simulation Monte Carlo (DSMC) method is the method of choice to simulate high-speed noncontinuum flows. However, since DSMC uses computational molecules to represent the gas, the inherent statistical noise must be minimized by sampling large numbers of molecules. Since typical microsystem velocities are low (< 1 m/s) compared to molecular velocities ({approx}400 m/s), the number of molecular samples required to achieve 1% precision can exceed 1010 per cell. The Discrete Velocity Gas (DVG) method, an approach motivated by radiation transport, provides another way to simulate noncontinuum gas flows. Unlike DSMC, the DVG method restricts molecular velocities to have only certain discrete values. The transport of the number density of a velocity state is governed by a discrete Boltzmann equation that has one temporal dimension and three spatial dimensions and a polynomial collision term. Specification and implementation of DVG models are discussed, and DVG models are applied to Couette flow and to Fourier flow. While the DVG results for these benchmark problems are qualitatively correct, the errors in the shear stress and the heat flux can be order-unity even for DVG models with 88 velocity states. It is concluded that the DVG method, as described herein, is not sufficiently accurate to simulate the low-speed gas flows that occur in microsystems.

Torczynski, J.R.; Torczynski, J.R.; O'Hern, Timothy J.; Barney, Jeremy B.; Castaneda, Jaime N.; Cote, Raymond O.

An experimental program is being conducted to study a proposed approach for oil reintroduction in the Strategic Petroleum Reserve (SPR). The goal is to assess whether useful oil is rendered unusable through formation of a stable oil-brine emulsion during reintroduction of degassed oil into the brine layer in storage caverns. This report documents the first stage of the program, in which simulant liquids are used to characterize the buoyant plume that is produced when a jet of crude oil is injected downward from a tube into brine. The experiment consists of a large transparent vessel that is a scale model of the proposed oil injection process at the SPR. An oil layer is floated on top of a brine layer. Silicon oil (Dow Corning 200{reg_sign} Fluid, 5 cSt) is used as the simulant for crude oil to allow visualization of the flow and to avoid flammability and related concerns. Sodium nitrate solution is used as the simulant for brine because it is not corrosive and it can match the density ratio between brine and crude oil. The oil is injected downward through a tube into the brine at a prescribed depth below the oil-brine interface. Flow rates are determined by scaling to match the ratio of buoyancy to momentum between the experiment and the SPR. Initially, the momentum of the flow produces a downward jet of oil below the tube end. Subsequently, the oil breaks up into droplets due to shear forces, buoyancy dominates the flow, and a plume of oil droplets rises to the interface. The interface is deflected upward by the impinging oil-brine plume. Two different diameter injection tubes were used (1/2-inch and 1-inch OD) to vary the scaling. Use of the 1-inch injection tube also assured that turbulent pipe flow was achieved, which was questionable for lower flow rates in the 1/2-inch tube. In addition, a 1/2-inch J-tube was used to direct the buoyant jet upwards rather than downwards to determine whether flow redirection could substantially reduce the oil-plume size and the oil-droplet residence time in the brine. Reductions of these quantities would inhibit emulsion formation by limiting the contact between the oil and the brine. Videos of this flow were recorded for scaled flow rates that bracket the equivalent pumping rates in an SPR cavern. Image-processing analyses were performed to quantify the penetration depth of the oil jet, the width of the jet, and the deflection of the interface. The measured penetration depths are shallow, as predicted by penetration-depth models, in agreement with the assumption that the flow is buoyancy-dominated, rather than momentum-dominated. The turbulent penetration depth model provided a good estimate of the measured values for the 1-inch injection tube but overpredicted the penetration depth for the 1/2-inch injection tube. Adding a virtual origin term would improve the prediction for the 1/2-inch tube for low to nominal injection flow rates but could not capture the rollover seen at high injection flow rates. As expected, the J-tube yielded a much narrower plume because the flow was directed upward, unlike the downward-oriented straight-tube cases where the plume had to reverse direction, leading to a much wider effective plume area. Larger surface deflections were caused by the narrower plume emitted from the J-tube. Although velocity was not measured in these experiments, the video data showed that the J-tube plume was clearly faster than those emitted from the downward-oriented tubes. These results indicate that oil injection tube modifications could inhibit emulsion formation by reducing the amount of contact (both time and area) between the oil and the brine. Future studies will employ crude oil, saturated brine, and interfacial solids (sludge) from actual SPR caverns.

Proposed for publication in Aerosol Science and Technology.

Gallis, Michail A.; Gallis, Michail A.; Rader, Daniel J.; Torczynski, J.R.

A general, approximate expression is described that can be used to predict the thermophoretic force on a free-molecular, motionless, spherical particle suspended in a quiescent gas with a temperature gradient. The thermophoretic force is equal to the product of an order-unity coefficient, the gas-phase translational heat flux, the particle cross-sectional area, and the inverse of the mean molecular speed. Numerical simulations are used to test the accuracy of this expression for monatomic gases, polyatomic gases, and mixtures thereof. Both continuum and noncontinuum conditions are examined; in particular, the effects of low pressure, wall proximity, and high heat flux are investigated. The direct simulation Monte Carlo (DSMC) method is used to calculate the local molecular velocity distribution, and the force-Green's-function method is used to calculate the thermophoretic force. The approximate expression is found to predict the calculated thermophoretic force to within 10% for all cases examined.

Gallis, Michail A.; Gallis, Michail A.; Torczynski, J.R.

Abstract not provided.

Trujillo, Steven M.; Tortora, Paul R.; Ceccio, Steven L.; O'Hern, Timothy J.; Torczynski, J.R.; Trujillo, Steven M.

Abstract not provided.

Liter, Scott L.; Torczynski, J.R.

A novel electrical-impedance tomography (EIT) diagnostic system, including hardware and software, has been developed and used to quantitatively measure material distributions in multiphase flows within electrically-conducting (i.e., industrially relevant or metal) vessels. The EIT system consists of energizing and measuring electronics and seven ring electrodes, which are equally spaced on a thin nonconducting rod that is inserted into the vessel. The vessel wall is grounded and serves as the ground electrode. Voltage-distribution measurements are used to numerically reconstruct the time-averaged impedance distribution within the vessel, from which the material distributions are inferred. Initial proof-of-concept and calibration was completed using a stationary solid-liquid mixture in a steel bench-top standpipe. The EIT system was then deployed in Sandia's pilot-scale slurry bubble-column reactor (SBCR) to measure material distributions of gas-liquid two-phase flows over a range of column pressures and superficial gas flow rates. These two-phase quantitative measurements were validated against an established gamma-densitometry tomography (GDT) diagnostic system, demonstrating agreement to within 0.05 volume fraction for most cases, with a maximum difference of 0.15 volume fraction. Next, the EIT system was combined with the GDT system to measure material distributions of gas-liquid-solid three-phase flows in Sandia's SBCR for two different solids loadings. Accuracy for the three-phase flow measurements is estimated to be within 0.15 volume fraction. The stability of the energizing electronics, the effect of the rod on the surrounding flow field, and the unsteadiness of the liquid temperature all degrade measurement accuracy and need to be explored further. This work demonstrates that EIT may be used to perform quantitative measurements of material distributions in multiphase flows in metal vessels.

Torczynski, J.R.; Gallis, Michail A.; Rader, Daniel J.

An approach is presented to compute the force on a spherical particle in a rarefied flow of a monatomic gas. This approach relies on the development of a Green's function that describes the force on a spherical particle in a delta-function molecular velocity distribution function. The gas-surface interaction model in this development allows incomplete accommodation of energy and tangential momentum. The force from an arbitrary molecular velocity distribution is calculated by computing the moment of the force Green's function in the same way that other macroscopic variables are determined. Since the molecular velocity distribution function is directly determined in the DSMC method, the force Green's function approach can be implemented straightforwardly in DSMC codes. A similar approach yields the heat transfer to a spherical particle in a rarefied gas flow. The force Green's function is demonstrated by application to two problems. First, the drag force on a spherical particle at arbitrary temperature and moving at arbitrary velocity through an equilibrium motionless gas is found analytically and numerically. Second, the thermophoretic force on a motionless particle in a motionless gas with a heat flux is found analytically and numerically. Good agreement is observed in both situations.

Stuecker, John N.; Cesarano, Joseph C.; Shollenberger, K.A.; Roach, R.A.; Torczynski, J.R.; Thomas, Edward V.; Van Ornum, David J.

This report is a summary of the work completed in FY00 for science-based characterization of the processes used to fabricate cermet vias in source feedthrus. In particular, studies were completed to characterize the CND50 cermet slurry, characterize solvent imbibition, and identify critical via filling variables. These three areas of interest are important to several processes pertaining to the production of neutron generator tubes. Rheological characterization of CND50 slurry prepared with 94ND2 and Sandi94 primary powders were also compared. The 94ND2 powder was formerly produced at the GE Pinellas Plant and the Sandi94 is the new replacement powder produced at CeramTec. Processing variables that may effect the via-filling process were also studied and include: the effect of solids loading in the CND50 slurry; the effect of milling time; and the effect of Nuosperse (a slurry ''conditioner''). Imbibition characterization included a combination of experimental, theoretical, and computational strategies to determine solvent migration though complex shapes, specifically vias in the source feedthru component. Critical factors were determined using a controlled set of experiments designed to identify those variables that influence the occurrence of defects within the cermet filled via. These efforts were pursued to increase part production reliability, understand selected fundamental issues that impact the production of slurry-filled parts, and validate the ability of the computational fluid dynamics code, GOMA, to simulate these processes. Suggestions are made for improving the slurry filling of source feedthru vias.

Trujillo, Steven M.; Shollenberger, K.A.; O'Hern, Timothy J.; Grasser, Thomas W.; Torczynski, J.R.

Abstract not provided.

Physics of Fluids

Gallis, Michail A.; Torczynski, J.R.

Abstract not provided.

American Society of Mechanical Engineers, Fluids Engineering Division (Publication) FED

George, D.L.; Shollenberger, K.A.; Torczynski, J.R.

Gamma-densitometry tomography is applied to study the effects of sparger hole geometry, gas flow rate, column pressure, and phase properties on gas volume fraction profiles in bubble columns. Tests are conducted in a column 0.48 m in diameter, using air and mineral oil, superficial gas velocities ranging from 5 to 30 cm s-1, and absolute column pressures from 103 to 517 kPa. Reconstructed gas volume fraction profiles from two sparger geometries are presented. The development length of the gas volume fraction profile is found to increase with gas flow rate and column pressure. Increases in gas flow rate increase the local gas volume fraction preferentially on the column axis, whereas increases in column pressure produce a uniform rise in gas volume fraction across the column. A comparison of results from the two spargers indicates a significant change in development length with the number and size of sparger holes.

Gallis, Michail A.; Torczynski, J.R.

A new scheme to simulate elastic collisions in particle simulation codes is presented. The new scheme aims at simulating the collisions in the highly collisional regime, in which particle simulation techniques typically become computationally expensive.The new scheme is based on the concept of a grid-based collision field. According to this scheme, the particles perform a single collision with the background grid during a time step. The properties of the background field are calculated from the moments of the distribution function accumulated on the grid. The collision operator is based on the Langevin equation. Based on comparisons with other methods, it is found that the Langevin method overestimates the collision frequency for dilute gases.

Gallis, Michail A.; Torczynski, J.R.

Abstract not provided.

International Journal of Multiphase Flow

Shollenberger, K.A.; Torczynski, J.R.; O'Hern, Timothy J.; Shollenberger, K.A.

Experiments are presented in which electrical-impedance tomography (EIT) and gamma-densitometry tomography (GDT) measurements were combined to simultaneously measure the solid, liquid, and gas radial distributions in a vertical three-phase flow. The experimental testbed was a 19.05-cm diameter bubble column in which gas is injected at the bottom and exits out the top while the liquid and solid phases recirculate. The gas phase was air and the liquid phase was deionized water with added electrolytes. Four different particle classes were investigated for the solid phase: 40--100 {micro}m and 120--200 {micro}m glass beads (2.41 g/cm{sup 3}), and 170--260 {micro}m and 200--700 {micro}m polystyrene beads (1.04 g/cm{sup 3}). Superficial gas velocities of 3 to 30 cm/s and solid volume fractions up to 0.30 were examined. For all experimental conditions investigated, the gas distribution showed only a weak dependence on both particle size and density. Average gas volume fraction as a function of superficial gas velocity can be described to within {+-} 0.04 by curve passing through the center of the data. For most cases the solid particle appeared to be radically uniformly dispersed in the liquid.

Torczynski, J.R.; Shollenberger, K.A.; O'Hern, Timothy J.

An electrical-impedance tomography (EIT) system has been developed for quantitative measurements of radial phase distribution profiles in two-phase and three-phase vertical column flows. The EIT system is described along with the computer algorithm used for reconstructing phase volume fraction profiles. EIT measurements were validated by comparison with a gamma-densitometry tomography (GDT) system. The EIT system was used to accurately measure average solid volume fractions up to 0.05 in solid-liquid flows, and radial gas volume fraction profiles in gas-liquid flows with gas volume fractions up to 0.15. In both flows, average phase volume fractions and radial volume fraction profiles from GDT and EIT were in good agreement. A minor modification to the formula used to relate conductivity data to phase volume fractions was found to improve agreement between the methods. GDT and EIT were then applied together to simultaneously measure the solid, liquid, and gas radial distributions within several vertical three-phase flows. For average solid volume fractions up to 0.30, the gas distribution for each gas flow rate was approximately independent of the amount of solids in the column. Measurements made with this EIT system demonstrate that EIT may be used successfully for noninvasive, quantitative measurements of dispersed multiphase flows.