No abstract

Spray cooling with phase change offers the means to achieve the highest rates of heat transfer from microelectronic components and other high energy density devices. The extreme complexity of the flow created by the impact of millions of droplets per second per cm 2 creates a need for a heat transfer design model which incorporates enough physical detail to yield accurate predictions while being sufficiently simplified to allow its use in routine design computations. The spray cooling group at West Virginia University is pursuing a coordinated program of laboratory experiments and computational simulations to develop a Monte Carlo-based spray cooling model that will satisfy this requirement. This paper reports progress in both the laboratory and computational fluid dynamics (CFD) phases of this study, including comparisons between single droplet impact results at We values of 161 and 633 to validate the CFD simulations and experimental measurements. A key goal of the present work is to determine the thickness of the thin liquid layer remaining in the crater formed when a liquid droplet impacts a surface covered by a preexisting liquid film. The volume of liquid in this thin liquid film in the droplet impact crater is one factor that will influence the onset of critical heat flux (CHF). Based on the present data the preexisting static film thickness is observed to have more of an effect than the initial drop parameters on the minimum liquid volume under the crater, but the crater lifetime depends on the initial drop conditions. The comparisons between experiments and simulations for these two cases show promise, but more refinement in experimental and computational technique are needed to achieve more consistent determinations of the volume of liquid under a crater. NomenclatureD = Drop diameter D U = upper crown diameter D V = Initial drop volume Fr = Froude number = H = Preexisting static liquid pool thickness H crown = Height of droplet impact crown h o = Centerline liquid film thickness h avg = Average liquid film thickness Oh = Ohnesorge number = Re = Reynolds number = 2 r = radius, measured from center of droplet impact R B = Radius of outside bottom of crown t = Time We = Weber number = V = Impact velocity Vol = Volume liquid under the impact crater Y = splashing-deposition criteria ρ = Density μ = Viscosity σ = Surface tension

Spray cooling is a topic of current interest for its ability to uniformly remove high levels of waste heat for densely packed microelectronics. A Monte-Carlo (MC) spray cooling simulation model is under development that is based on empirical data to be a cost effective design tool that will predict accurate heat fluxes based on nozzle conditions and heater geometry. This work reports spray and single drop experiments with the goal of computing the volume beneath a drop impact cavity (sub-cavity volume) created by a single impinging droplet on an initial liquid layer. A relevant test plan for the single drop experiments in terms of We and Re numbers was created through utilization of Phase Doppler Anemometry to characterize a water spray generated by a nozzle of interest for varying flow conditions. Liquid thickness profiles of the sub-cavity formed by a single impinging drop onto a range of initial liquid film thicknesses were measured versus time via a non-contact optical thickness sensor. Time dependent sub-cavity volumes were computed by integrating these sub-cavity liquid film thickness profiles measured radially outward from the impact centerline. It is found that higher We and lower h 0 * result in a more radially uniform sub-cavity surface contour versus time, except for regions near the outer bottom cavity radius, where the liquid film was thinner. The sub-cavity volume was found to be nearly constant for a majority of the cavity lifetime and increased with We and h 0 * . These results will be incorporated in future work into the MC model to improve its predictive capability. Nomenclatured = Drop diameter Greek Letters D = Arithmetic mean droplet diameter ρ = Density D 32 = Sauter mean droplet diameter η = Index of refraction Fr = Froude number = μ = Dynamic viscosity h = Liquid film thickness σ = Surface tension R = Radial location τ = Dimensionless time (t•V axial /d) Re = Reynolds number = Superscripts T = Temperature * = Dimensionless parameter t = Time Subscripts V = Arithmetic mean velocity 0 = Initial condition Vol = Sub-cavity volume axial = Axial velocity component We = Weber number = c = Cavity z = Axial standoff distance from the nozzle tip r = Radial velocity component s = Spray

No abstract

A methodology is developed for predicting the performance of near-surface CO(2) leak detection systems at geologic sequestration sites. The methodology integrates site characterization and modeling to predict the statistical properties of natural CO(2) fluxes, the transport of CO(2) from potential subsurface leakage points, and the detection of CO(2) surface fluxes by the monitoring network. The probability of leak detection is computed as the probability that the leakage signal is sufficient to increase the total flux beyond a statistically determined threshold. The methodology is illustrated for a highly idealized site monitored with CO(2) accumulation chamber measurements taken on a uniform grid. The TOUGH2 code is used to predict the spatial profile of surface CO(2) fluxes resulting from different leakage rates and different soil permeabilities. A response surface is fit to the TOUGH2 results to allow interpolation across a continuous range of values of permeability and leakage rate. The spatial distribution of leakage probability is assumed uniform in this application. Nonlinear, nonmonotonic relationships of network performance to soil permeability and network density are evident. In general, dense networks (with ∼10-20 m between monitors) are required to ensure a moderate to high probability of leak detection.