Publications

39 Results
Skip to search filters

Interaction of a planar shock with a dense field of particles in a multiphase shock tube

Beresh, Steven J.; Kearney, S.P.; Trott, Wayne T.; Castaneda, Jaime N.; Pruett, Brian O.; Baer, Melvin B.

A novel multiphase shock tube has been constructed to test the interaction of a planar shock wave with a dense gas-solid field of particles. The particle field is generated by a gravity-fed method that results in a spanwise curtain of 100-micron particles producing a volume fraction of about 15%. Interactions with incident shock Mach numbers of 1.67 and 1.95 are reported. High-speed schlieren imaging is used to reveal the complex wave structure associated with the interaction. After the impingement of the incident shock, transmitted and reflected shocks are observed, which lead to differences in flow properties across the streamwise dimension of the curtain. Tens of microseconds after the onset of the interaction, the particle field begins to propagate downstream, and disperse. The spread of the particle field, as a function of its position, is seen to be nearly identical for both Mach numbers. Immediately downstream of the curtain, the peak pressures associated with the Mach 1.67 and 1.95 interactions are about 35% and 45% greater than tests without particles, respectively. For both Mach numbers tested, the energy and momentum fluxes in the induced flow far downstream are reduced by about 30-40% by the presence of the particle field.

More Details

Surface rheology and interface stability

Mondy, L.A.; Brooks, Carlton F.; Grillet, Anne M.; Moffat, Harry K.; Koehler, Timothy P.; Yaklin, Melissa A.; Cote, Raymond O.; Castaneda, Jaime N.

We have developed a mature laboratory at Sandia to measure interfacial rheology, using a combination of home-built, commercially available, and customized commercial tools. An Interfacial Shear Rheometer (KSV ISR-400) was modified and the software improved to increase sensitivity and reliability. Another shear rheometer, a TA Instruments AR-G2, was equipped with a du Nouey ring, bicone geometry, and a double wall ring. These interfacial attachments were compared to each other and to the ISR. The best results with the AR-G2 were obtained with the du Nouey ring. A Micro-Interfacial Rheometer (MIR) was developed in house to obtain the much higher sensitivity given by a smaller probe. However, it was found to be difficult to apply this technique for highly elastic surfaces. Interfaces also exhibit dilatational rheology when the interface changes area, such as occurs when bubbles grow or shrink. To measure this rheological response we developed a Surface Dilatational Rheometer (SDR), in which changes in surface tension with surface area are measured during the oscillation of the volume of a pendant drop or bubble. All instruments were tested with various surfactant solutions to determine the limitations of each. In addition, foaming capability and foam stability were tested and compared with the rheology data. It was found that there was no clear correlation of surface rheology with foaming/defoaming with different types of surfactants, but, within a family of surfactants, rheology could predict the foam stability. Diffusion of surfactants to the interface and the behavior of polyelectrolytes were two subjects studied with the new equipment. Finally, surface rheological terms were added to a finite element Navier-Stokes solver and preliminary testing of the code completed. Recommendations for improved implementation were given. When completed we plan to use the computations to better interpret the experimental data and account for the effects of the underlying bulk fluid.

More Details

Determination of thermal accommodation coefficients from heat transfer measurements between parallel plates

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

Thermal accommodation coefficients have been derived for a variety of gas-surface combinations using an experimental apparatus developed to measure the pressure dependence of the conductive heat flux between parallel plates at unequal temperature separated by a gas-filled gap. The heat flux is inferred from temperature-difference measurements across the plates in a configuration where the plate temperatures are set with two carefully controlled thermal baths. Temperature-controlled shrouds provide for environmental isolation of the opposing test plates. Since the measured temperature differences in these experiments are very small (typically 0.3 C or less over the entire pressure range), high-precision thermistors are used to acquire the requisite temperature data. High-precision components have also been utilized on the other control and measurement subsystems in this apparatus, including system pressure, gas flow rate, plate alignment, and plate positions. The apparatus also includes the capability for in situ plasma cleaning of the installed test plates. Measured heat-flux results are used in a formula based on Direct Simulation Monte Carlo (DSMC) code calculations to determine the thermal accommodation coefficients. Thermal accommodation coefficients have been determined for three different gases (argon, nitrogen, helium) in contact with various surfaces. Materials include metals and alloys such as aluminum, gold, platinum, and 304 stainless steel. A number of materials important to fabrication of Micro Electro Mechanical Systems (MEMS) devices have also been examined. For most surfaces, coefficient values are near 0.95, 0.85, and 0.45 for argon, nitrogen, and helium, respectively. Only slight differences in accommodation as a function of surface roughness have been seen. Surface contamination appears to have a more significant effect: argon plasma treatment has been observed to reduce thermal accommodation by as much as 0.10 for helium. Mixtures of argon and helium have also been examined, and the results have been compared to DSMC simulations incorporating thermal-accommodation values from single-species experiments.

More Details

Development of a multiphase shock tube for energetic materials characterization

Baer, Melvin B.; Cooper, Marcia A.; Castaneda, Jaime N.; Beresh, Steven J.; Pruett, Brian O.; Kearney, S.P.; Trott, Wayne T.

A novel multiphase shock tube to study particle dynamics in gas-solid flows has been constructed and tested. Currently, there is a gap in data for flows having particle volume fractions between the dusty and granular regimes. The primary purpose of this new facility is to fill that gap by providing high quality data of shock-particle interactions in flows having dense gas particle volume fractions. Towards this end, the facility aims to drive a shock into a spatially isotropic field, or curtain, of particles. Through bench-top experimentation, a method emerged for achieving this challenging task that involves the use of a gravity-fed contoured particle seeder. The seeding method is capable of producing fields of spatially isotropic particles having volume fractions of about 1 to 35%. The use of the seeder in combination with the shock tube allows for the testing of the impingement of a planar shock on a dense field of particles. The first experiments in the multiphase shock tube have been conducted and the facility is now operational.

More Details

Investigation of interfacial rheology & foam stability

Grillet, Anne M.; Koehler, Timothy P.; Reichert, Matthew D.; Yaklin, Melissa A.; Brooks, Carlton F.; Mondy, L.A.; Castaneda, Jaime N.; Cote, Raymond O.

The rheology at gas-liquid interfaces strongly influences the stability and dynamics of foams and emulsions. Several experimental techniques are employed to characterize the rheology at liquid-gas interfaces with an emphasis on the non-Newtonian behavior of surfactant-laden interfaces. The focus is to relate the interfacial rheology to the foamability and foam stability of various aqueous systems. An interfacial stress rheometer (ISR) is used to measure the steady and dynamic rheology by applying an external magnetic field to actuate a magnetic needle suspended at the interface. Results are compared with those from a double wall ring attachment to a rotational rheometer (TA Instruments AR-G2). Micro-interfacial rheology (MIR) is also performed using optical tweezers to manipulate suspended microparticle probes at the interface to investigate the steady and dynamic rheology. Additionally, a surface dilatational rheometer (SDR) is used to periodically oscillate the volume of a pendant drop or buoyant bubble. Applying the Young-Laplace equation to the drop shape, a time-dependent surface tension can be calculated and used to determine the effective dilatational viscosity of an interface. Using the ISR, double wall ring, SDR, and MIR, a wide range of sensitivity in surface forces (fN to nN) can be explored as each experimental method has different sensitivities. Measurements will be compared to foam stability.

More Details

Diagnostic development for determining the joint temperature/soot statistics in hydrocarbon-fueled pool fires : LDRD final report

Frederickson, Kraig; Grasser, Thomas W.; Castaneda, Jaime N.; Hewson, John C.; Luketa, Anay L.

A joint temperature/soot laser-based optical diagnostic was developed for the determination of the joint temperature/soot probability density function (PDF) for hydrocarbon-fueled meter-scale turbulent pool fires. This Laboratory Directed Research and Development (LDRD) effort was in support of the Advanced Simulation and Computing (ASC) program which seeks to produce computational models for the simulation of fire environments for risk assessment and analysis. The development of this laser-based optical diagnostic is motivated by the need for highly-resolved spatio-temporal information for which traditional diagnostic probes, such as thermocouples, are ill-suited. The in-flame gas temperature is determined from the shape of the nitrogen Coherent Anti-Stokes Raman Scattering (CARS) signature and the soot volume fraction is extracted from the intensity of the Laser-Induced Incandescence (LII) image of the CARS probed region. The current state of the diagnostic will be discussed including the uncertainty and physical limits of the measurements as well as the future applications of this probe.

More Details

High speed (GHZ), ultra-high pressure (GPA) sensor array fabricated in integrated CMOS+MEMS process

Proceedings of the IEEE International Conference on Micro Electro Mechanical Systems (MEMS)

Okandan, Murat O.; Olsson, Roy H.; Baker, Michael; Resnick, Paul J.; Hill, Thomas A.; Lackey, Chad; Pearson, Sean P.; Castaneda, Jaime N.; Trott, Wayne T.; Jones, David A.

In order to observe and quantify pressure levels generated during testing of energetic materials, a sensor array with high temporal resolution (∼1 ns) and extremely high pressure range (> 1 GPa) is needed. We have developed such a sensor array which utilizes a novel integrated high performance CMOS+MEMS process. ©2009 IEEE.

More Details

Experiments for foam model development and validation

Mondy, L.A.; Gorby, Allen D.; Cote, Raymond O.; Castaneda, Jaime N.; Thompson, Kyle R.; Rao, Rekha R.; Moffat, Harry K.; Kraynik, Andrew M.; Russick, Edward M.; Adolf, Douglas B.; Grillet, Anne M.; Brotherton, Christopher M.; Bourdon, Christopher B.

A series of experiments has been performed to allow observation of the foaming process and the collection of temperature, rise rate, and microstructural data. Microfocus video is used in conjunction with particle image velocimetry (PIV) to elucidate the boundary condition at the wall. Rheology, reaction kinetics and density measurements complement the flow visualization. X-ray computed tomography (CT) is used to examine the cured foams to determine density gradients. These data provide input to a continuum level finite element model of the blowing process.

More Details

Wetting and free surface flow modeling for potting and encapsulation

Adolf, Douglas B.; Castaneda, Jaime N.; Kraynik, Andrew M.; Noble, David R.; Sun, Amy C.; Cote, Raymond O.; Grillet, Anne M.; Notz, Patrick N.; Brooks, Carlton F.; Givler, R.C.; Hopkins, Matthew M.; Mondy, L.A.; Rao, Rekha R.

As part of an effort to reduce costs and improve quality control in encapsulation and potting processes the Technology Initiative Project ''Defect Free Manufacturing and Assembly'' has completed a computational modeling study of flows representative of those seen in these processes. Flow solutions are obtained using a coupled, finite-element-based, numerical method based on the GOMA/ARIA suite of Sandia flow solvers. The evolution of the free surface is solved with an advanced level set algorithm. This approach incorporates novel methods for representing surface tension and wetting forces that affect the evolution of the free surface. In addition, two commercially available codes, ProCAST and MOLDFLOW, are also used on geometries representing encapsulation processes at the Kansas City Plant. Visual observations of the flow in several geometries are recorded in the laboratory and compared to the models. Wetting properties for the materials in these experiments are measured using a unique flowthrough goniometer.

More Details

Experimental measurements of thermal accommodation coefficients for microscale gas-phase heat transfer

Collection of Technical Papers - 39th AIAA Thermophysics Conference

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

An experimental apparatus is described that measures gas-surface thermal accommodation coefficients from the pressure dependence of the conductive heat flux between parallel plates separated by a gas-filled gap. Heat flux between the plates is inferred from measurements of temperature drop between the plate surface and an adjacent temperature-controlled water bath. Thermal accommodation coefficients are determined from the pressure dependence of the heat flux at a fixed plate separation. The apparatus is designed to conduct tests with a variety of gases in contact with interchangeable, well-characterized surfaces of various materials (e.g., metals, ceramics, semiconductors) with various surface finishes (e.g., smooth, rough). Experiments are reported for three gases (argon, nitrogen, and helium) in contact with pairs of 304 stainless steel plates prepared with one of two finishes: lathe-machined or mirror-polished. For argon and nitrogen, the measured accommodation coefficients for machined and polished plates are near unity and independent of finish to within experimental uncertainty. For helium, the accommodation coefficients are much lower and show a slight variation with surface roughness. Two different methods are used to determine the accommodation coefficient from experimental data: the Sherman-Lees formula and the GTR formula. These approaches yield values of 0.87 and 0.94 for argon, 0.80 and 0.86 for nitrogen, 0.36 and 0.38 for helium with the machined finish, and 0.40 and 0.42 for helium with the polished finish, respectively, with an uncertainty of ±0.02. The GTR values for argon and nitrogen are generally in better agreement with the results of other investigators than the Sherman-Lees values are, and both helium results are in reasonable agreement with values in the literature.

More Details

Evaluation of aluminum participation in the development of reactive waves in shock compressed HMX

Proceedings of the 13th International Detonation Symposium, IDS 2006

Pahl, Robert J.; Trott, Wayne T.; Castaneda, Jaime N.; Marley, Stephen K.; Snedigar, Shane

Nanometric aluminum (123nm, spherical) was mixed with two different sieve-cut sizes of HMX (106-150 μm and 212-300 μm), and a series of gas gun tests were conducted to compare reactive wave development in pure HMX to that of aluminized HMX. In the absence of added metal, 4-mm-thick, low-density (68% of theoretical maximum density) pressings of the 106-150 μm HMX respond to modest shock loading by developing distinctive reactive waves that exhibit both temporal and meso-scale spatial fluctuations. Similar pressings of Al/HMX containing 10% aluminum (by mass) show an initial suppression of the usual wave growth seen in HMX samples. The suppression is then followed by an induction period where it is hypothesized that a phase change in the aluminum may occur. Data from VISAR, line-ORVIS, and 2-color pyrometry are given and discussed, and numerical modeling of inert sucrose is used to aid the explanation of the resulting data.

More Details

Modeling injection molding of net-shape active ceramic components

Rao, Rekha R.; Brooks, Carlton F.; Cote, Raymond O.; Castaneda, Jaime N.; Mondy, L.A.; Noble, David R.; Hopkins, Matthew M.; Notz, Patrick N.; Halbleib, Laura L.; Yang, Pin Y.; Burns, George B.; Grillet, Anne M.

To reduce costs and hazardous wastes associated with the production of lead-based active ceramic components, an injection molding process is being investigated to replace the current machining process. Here, lead zirconate titanate (PZT) ceramic particles are suspended in a thermoplastic resin and are injected into a mold and allowed to cool. The part is then bisque fired and sintered to complete the densification process. To help design this new process we use a finite element model to describe the injection molding of the ceramic paste. Flow solutions are obtained using a coupled, finite-element based, Newton-Raphson numerical method based on the GOMA/ARIA suite of Sandia flow solvers. The evolution of the free surface is solved with an advanced level set algorithm. This approach incorporates novel methods for representing surface tension and wetting forces that affect the evolution of the free surface. Thermal, rheological, and wetting properties of the PZT paste are measured for use as input to the model. The viscosity of the PZT is highly dependent both on temperature and shear rate. One challenge in modeling the injection process is coming up with appropriate constitutive equations that capture relevant phenomenology without being too computationally complex. For this reason we model the material as a Carreau fluid and a WLF temperature dependence. Two-dimensional (2D) modeling is performed to explore the effects of the shear in isothermal conditions. Results indicate that very low viscosity regions exist near walls and that these results look similar in terms of meniscus shape and fill times to a simple Newtonian constitutive equation at the shear-thinned viscosity for the paste. These results allow us to pick a representative viscosity to use in fully three-dimensional (3D) simulation, which because of numerical complexities are restricted to using a Newtonian constitutive equation. Further 2D modeling at nonisothermal conditions shows that the choice of representative Newtonian viscosity is dependent on the amount of heating of the initially room temperature mold. An early 3D transient model shows that the initial design of the distributor is sub-optimal. However, these simulations take several months to run on 4 processors of an HP workstation using a preconditioner/solver combination of ILUT/GMRES with fill factors of 3 and PSPG stabilization. Therefore, several modifications to the distributor geometry and orientations of the vents and molds have been investigated using much faster 3D steady-state simulations. The pressure distribution for these steady-state calculations is examined for three different distributor designs to see if this can indicate which geometry has the superior design. The second modification, with a longer distributor, is shown to have flatter, more monotonic isobars perpendicular to the flow direction indicating a better filling process. The effects of the distributor modifications, as well as effects of the mold orientation, have also been examined with laboratory experiments in which the flow of a viscous Newtonian oil entering transparent molds is recorded visually. Here, the flow front is flatter and voids are reduced for the second geometry compared to the original geometry. A horizontal orientation, as opposed to the planned vertical orientation, results in fewer voids. Recently, the Navier-Stokes equations have been stabilized with the Dohrman-Bochev PSPP stabilization method, allowing us to calculate transient 3D simulations with computational times on the order of days instead of months. Validation simulations are performed and compared to the experiments. Many of the trends of the experiments are captured by the level set modeling, though quantitative agreement is lacking mainly due to the high value of the gas phase viscosity necessary for numerical stability, though physically unrealistic. More correct trends are predicted for the vertical model than the horizontal model, which is serendipitous as the actual mold is held in a vertical geometry. The full, transient mold filling calculations indicate that the flow front is flatter and voids may be reduced for the second geometry compared to the original geometry. The validated model is used to predict mold filling for the actual process with the material properties for the PZT paste, the original distributor geometry, and the mold in a vertical orientation. This calculation shows that voids may be trapped at the four corners of the mold opposite the distributor.

More Details

Shock-induced explosive chemistry in a deterministic sample configuration

Trott, Wayne T.; Baer, Melvin B.; Castaneda, Jaime N.; Tappan, Alexander S.

Explosive initiation and energy release have been studied in two sample geometries designed to minimize stochastic behavior in shock-loading experiments. These sample concepts include a design with explosive material occupying the hole locations of a close-packed bed of inert spheres and a design that utilizes infiltration of a liquid explosive into a well-defined inert matrix. Wave profiles transmitted by these samples in gas-gun impact experiments have been characterized by both velocity interferometry diagnostics and three-dimensional numerical simulations. Highly organized wave structures associated with the characteristic length scales of the deterministic samples have been observed. Initiation and reaction growth in an inert matrix filled with sensitized nitromethane (a homogeneous explosive material) result in wave profiles similar to those observed with heterogeneous explosives. Comparison of experimental and numerical results indicates that energetic material studies in deterministic sample geometries can provide an important new tool for validation of models of energy release in numerical simulations of explosive initiation and performance.

More Details

Measurements of thermal accommodation coefficients

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

A previously-developed experimental facility has been used to determine gas-surface thermal accommodation coefficients from the pressure dependence of the heat flux between parallel plates of similar material but different surface finish. Heat flux between the plates is inferred from measurements of temperature drop between the plate surface and an adjacent temperature-controlled water bath. Thermal accommodation measurements were determined from the pressure dependence of the heat flux for a fixed plate separation. Measurements of argon and nitrogen in contact with standard machined (lathed) or polished 304 stainless steel plates are indistinguishable within experimental uncertainty. Thus, the accommodation coefficient of 304 stainless steel with nitrogen and argon is estimated to be 0.80 {+-} 0.02 and 0.87 {+-} 0.02, respectively, independent of the surface roughness within the range likely to be encountered in engineering practice. Measurements of the accommodation of helium showed a slight variation with 304 stainless steel surface roughness: 0.36 {+-} 0.02 for a standard machine finish and 0.40 {+-} 0.02 for a polished finish. Planned tests with carbon-nanotube-coated plates will be performed when 304 stainless-steel blanks have been successfully coated.

More Details

Mold-filling experiments for validation of modeling encapsulation. Part 1, "wine glass" mold

Mondy, L.A.; Grillet, Anne M.; Cote, Raymond O.; Castaneda, Jaime N.

The C6 project 'Encapsulation Processes' has been designed to obtain experimental measurements for discovery of phenomena critical to improving these processes, as well as data required in the verification and validation plan (Rao et al. 2001) for model validation of flow in progressively complex geometries. We have observed and recorded the flow of clear, Newtonian liquids and opaque, rheologically complex suspensions in two mold geometries. The first geometry is a simple wineglass geometry in a cylinder and is reported here in Part 1. The results in a more realistic encapsulation geometry are reported in Part 2.

More Details

Microscale rarefied gas dynamics and surface interactions for EUVL and MEMS applications

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.

More Details

Investigation of oil injection into brine for the Strategic Petroleum Reserve : hydrodynamics and mixing experiments with SPR liquids

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.

More Details

Investigation of oil injection into brine for the strategic petroleum reserve : hydrodynamics experiments with simulant liquids

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.

More Details

Microscale Shock Wave Physics Using Photonic Driver Techniques

Setchell, Robert E.; Trott, Wayne T.; Castaneda, Jaime N.; Farnsworth, Archie V.; Berry, Dante M.

This report summarizes a multiyear effort to establish a new capability for determining dynamic material properties. By utilizing a significant reduction in experimental length and time scales, this new capability addresses both the high per-experiment costs of current methods and the inability of these methods to characterize materials having very small dimensions. Possible applications include bulk-processed materials with minimal dimensions, very scarce or hazardous materials, and materials that can only be made with microscale dimensions. Based on earlier work to develop laser-based techniques for detonating explosives, the current study examined the laser acceleration, or photonic driving, of small metal discs (''flyers'') that can generate controlled, planar shockwaves in test materials upon impact. Sub-nanosecond interferometric diagnostics were developed previously to examine the motion and impact of laser-driven flyers. To address a broad range of materials and stress states, photonic driving levels must be scaled up considerably from the levels used in earlier studies. Higher driving levels, however, increase concerns over laser-induced damage in optics and excessive heating of laser-accelerated materials. Sufficiently high levels require custom beam-shaping optics to ensure planar acceleration of flyers. The present study involved the development and evaluation of photonic driving systems at two driving levels, numerical simulations of flyer acceleration and impact using the CTH hydrodynamics code, design and fabrication of launch assemblies, improvements in diagnostic instrumentation, and validation experiments on both bulk and thin-film materials having well-established shock properties. The primary conclusion is that photonic driving techniques are viable additions to the methods currently used to obtain dynamic material properties. Improvements in launch conditions and diagnostics can certainly be made, but the main challenge to future applications will be the successful design and fabrication of test assemblies for materials of interest.

More Details

Preliminary Investigation of the Thermal Decomposition of Ablefoam and EF-AR20 Foam (Ablefoam Replacement)

Ulibarri, Tamara A.; Erickson, Kenneth L.; Wiemann, Dora K.; Erickson, Kenneth L.; Castaneda, Jaime N.; Borek, Theodore T.; Renlund, Anita M.; Miller, Jill C.

Preliminary thermal decomposition experiments with Ablefoam and EF-AR20 foam (Ablefoam replacement) were done to determine the important chemical and associated physical phenomena that should be investigated to develop the foam decomposition chemistry sub-models that are required in numerical simulations of the fire-induced response of foam-filled engineered systems for nuclear safety applications. Although the two epoxy foams are physically and chemically similar, the thermal decomposition of each foam involves different chemical mechanisms, and the associated physical behavior of the foams, particularly ''foaming'' and ''liquefaction,'' have significant implications for modeling. A simplified decomposition chemistry sub-model is suggested that, subject to certain caveats, may be appropriate for ''scoping-type'' calculations.

More Details

Dispersive Velocity Measurements in Heterogeneous Materials

Trott, Wayne T.; Castaneda, Jaime N.; Baer, Melvin B.; Chhabildas, Lalit C.; Knudson, Marcus D.; Davis, Jean-Paul D.; Asay, James R.

In order to provide real-time data for validation of three dimensional numerical simulations of heterogeneous materials subjected to impact loading, an optically recording velocity interferometer system (ORVIS) has been adapted to a line-imaging instrument capable of generating precise mesoscopic scale measurements of spatially resolved velocity variations during dynamic deformation. Combining independently variable target magnification and interferometer fringe spacing, this instrument can probe a velocity field along line segments up to 15 mm in length. In high magnification operation, spatial resolution better than 10 {micro}m can be achieved. For events appropriate to short recording times, streak camera recording can provide temporal resolution better than 0.2 ns. A robust method for extracting spatially resolved velocity-time profiles from streak camera image data has been developed and incorporated into a computer program that utilizes a standard VISAR analysis platform. The use of line-imaging ORVIS to obtain measurements of the mesoscopic scale dynamic response of shocked samples has been demonstrated on several different classes of heterogeneous materials. Studies have focused on pressed, granular sugar as a simulant material for the widely used explosive HMX. For low-density (65% theoretical maximum density) pressings of sugar, material response has been investigated as a function of both impact velocity and changes in particle size distribution. The experimental results provide a consistent picture of the dispersive nature of the wave transmitted through these samples and reveal both transverse and longitudinal wave structures on mesoscopic scales. This observed behavior is consistent with the highly structured mesoscopic response predicted by 3-D simulations. Preliminary line-imaging ORVIS measurements on HMX as well as other heterogeneous materials such as foam and glass-reinforced polyester are also discussed.

More Details
39 Results
39 Results