skip to main content


Title: Optimal cooling of an internally heated disc
Motivated by the search for sharp bounds on turbulent heat transfer as well as the design of optimal heat exchangers, we consider incompressible flows that most efficiently cool an internally heated disc. Heat enters via a distributed source, is passively advected and diffused, and exits through the boundary at a fixed temperature. We seek an advecting flow to optimize this exchange. Previous work on energy-constrained cooling with a constant source has conjectured that global optimizers should resemble convection rolls; we prove one-sided bounds on energy-constrained cooling corresponding to, but not resolving, this conjecture. In the case of an enstrophy constraint, our results are more complete: we construct a family of self-similar, tree-like ‘branching flows’ whose cooling we prove is within a logarithm of globally optimal. These results hold for general space- and time-dependent source–sink distributions that add more heat than they remove. Our main technical tool is a non-local Dirichlet-like variational principle for bounding solutions of the inhomogeneous advection–diffusion equation with a divergence-free velocity. This article is part of the theme issue ‘Mathematical problems in physical fluid dynamics (part 1)’.  more » « less
Award ID(s):
2025000
NSF-PAR ID:
10345194
Author(s) / Creator(s):
Date Published:
Journal Name:
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume:
380
Issue:
2225
ISSN:
1364-503X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We consider the problem of optimizing heat transport through an incompress- ible fluid layer. Modeling passive scalar transport by advection-diffusion, we maximize the mean rate of total transport by a divergence-free velocity field. Subject to various boundary conditions and intensity constraints, we prove that the maximal rate of transport scales linearly in the r.m.s. kinetic energy and, up to possible logarithmic corrections, as the one-third power of the mean enstro- phy in the advective regime. This makes rigorous a previous prediction on the near optimality of convection rolls for energy-constrained transport. On the other hand, optimal designs for enstrophy-constrained transport are significantly more difficult to describe: we introduce a “branching” flow design with an unbounded number of degrees of freedom and prove it achieves nearly optimal transport. The main technical tool behind these results is a variational principle for evalu- ating the transport of candidate designs. The principle admits dual formulations for bounding transport from above and below. While the upper bound is closely related to the “background method,” the lower bound reveals a connection be- tween the optimal design problems considered herein and other apparently re- lated model problems from mathematical materials science. These connections serve to motivate designs. 
    more » « less
  2. Abstract This work presents an approach to optimally designing a composite with thermal conductivity enhancers infiltrated with phase change material based on figure of merit (FOM) for thermal management of portable electronic devices. The FOM defines the balance between effective thermal conductivity and energy storage capacity. In this study, thermal conductivity enhancers are in the form of a honeycomb structure. Thermal conductivity enhancers are often used in conjunction with phase change material to enhance the conductivity of the composite medium. Under constrained heat sink volume, the higher volume fraction of thermal conductivity enhancers improves the effective thermal conductivity of the composite, while it reduces the amount of latent heat storage simultaneously. This work arrives at the optimal design of composite for electronic cooling by maximizing the FOM to resolve the stated tradeoff. In this study, the total volume of the composite and the interfacial heat transfer area between the phase change material and thermal conductivity enhancers are constrained for all design points. A benchmarked two-dimensional direct computational fluid dynamics model was employed to investigate the thermal performance of the phase change material and thermal conductivity enhancer composite. Furthermore, assuming conduction-dominated heat transfer in the composite, a simplified effective numerical model that solves the single energy equation with the effective properties of the phase change material and thermal conductivity enhancer has been developed. The effective properties like heat capacity can be obtained by volume averaging; however, effective thermal conductivity (required to calculate FOM) is unknown. The effective thermal conductivity of the composite is obtained by minimizing the error between the transient temperature gradient of direct and simplified model by iteratively varying the effective thermal conductivity. The FOM is maximized to find the optimal volume fraction for the present design. 
    more » « less
  3. UPS (Uninterruptible Power Supply) units and batteries are essential subsystems in data centers or telecom industries to protect equipment from electrical power spikes, surges and power outages. UPS units handle electrical power and dissipate a large amount of heat, and possess a high efficiency. Therefore, cooling units (e.g., CRACs) are needed to manage the thermal reliability of this equipment. On the other hand, battery operating conditions and reliability are closely related to the ambient temperature according to battery manufacturers; reliability increases when the ambient room temperature is around 25ºC. This study analyzed different room configurations and scenarios using the commercial CFD software 6Sigma Room DCXTM. As a first approach, we evaluated the thermal behavior and cooling degradation using standard thermal performance metrics SHI (Supply Heat Index) and RHI (Return Heat Index). These are frequently implemented in data centers to measure the level of mixing between cold and hot air streams. The results from this evaluation showed that standard cooling practices are inefficient, as values for the two metrics differed considerably from industry recommendations. We also considered a metric from the second law of thermodynamics using exergy destruction. This technique allowed us to find the mechanisms that increase entropy generation the most, including viscous shear and air stream mixing. Reducing exergy destruction will result in lessening lost thermodynamic work and thus reduce energy required for cooling. Typically, UPS and batteries are located in different rooms due to the hydrogen generation by the batteries. The integration of both equipment in the same room is a new concept, and this study aims to analyze the thermal performance of the room. Adding controllability showed improvements by reducing the exergy destruction due to viscous dissipation while slightly increasing thermal mixing in the rooms. Ducting the return flows to avoid flow mixing increased pressure drop, but reduced heat transfer between the hot and cold air streams, which in turn, improved the thermal performance. In the study, we determined the optimal configuration and possible strategies to improve cooling while maintaining desirable battery temperatures. 
    more » « less
  4. The most common approach to air cooling of data centers involves the pressurization of the plenum beneath the raised floor and delivery of air flow to racks via perforated floor tiles. This cooling approach is thermodynamically inefficient due in large part to the pressure losses through the tiles. Furthermore, it is difficult to control flow at the aisle and rack level since the flow source is centralized rather than distributed. Distributed cooling systems are more closely coupled to the heat generating racks. In overhead cooling systems, one can distribute flow to distinct aisles by placing the air mover and water cooled heat exchanger directly above an aisle. Two arrangements are possible: (i.) placing the air mover and heat exchanger above the cold aisle and forcing downward flow of cooled air into the cold aisle (Overhead Downward Flow (ODF)), or (ii.) placing the air mover and heat exchanger above the hot aisle and forcing heated air upwards from the hot aisle through the water cooled heat exchanger (Overhead Upward Flow (OUF)). This study focuses on the steady and transient behavior of overhead cooling systems in both ODF and OUF configurations and compares their cooling effectiveness and energy efficiency. The flow and heat transfer inside the servers and heat exchangers are modeled using physics based approaches that result in differential equation based mathematical descriptions. These models are programmed in the MATLAB™ language and embedded within a CFD computational environment (using the commercial code FLUENT™) that computes the steady or instantaneous airflow distribution. The complete computational model is able to simulate the complete flow and thermal field in the airside, the instantaneous temperatures within and pressure drops through the servers, and the instantaneous temperatures within and pressure drops through the overhead cooling system. Instantaneous overall energy consumption (1st Law) and exergy destruction (2nd Law) were used to quantify overall energy efficiency and to identify inefficiencies within the two systems. The server cooling effectiveness, based on an effectiveness-NTU model for the servers, was used to assess the cooling effectiveness of the two overhead cooling approaches 
    more » « less
  5. Abstract Energy exchanges between large-scale ocean currents and mesoscale eddies play an important role in setting the large-scale ocean circulation but are not fully captured in models. To better understand and quantify the ocean energy cycle, we apply along-isopycnal spatial filtering to output from an isopycnal 1/32° primitive equation model with idealized Atlantic and Southern Ocean geometry and topography. We diagnose the energy cycle in two frameworks: 1) a non-thickness-weighted framework, resulting in a Lorenz-like energy cycle, and 2) a thickness-weighted framework, resulting in the Bleck energy cycle. This paper shows that framework 2 is more useful for studying energy pathways when an isopycnal average is used. Next, we investigate the Bleck cycle as a function of filter scale. Baroclinic conversion generates mesoscale eddy kinetic energy over a wide range of scales and peaks near the deformation scale at high latitudes but below the deformation scale at low latitudes. Away from topography, an inverse cascade transfers kinetic energy from the mesoscales to larger scales. The upscale energy transfer peaks near the energy-containing scale at high latitudes but below the deformation scale at low latitudes. Regions downstream of topography are characterized by a downscale kinetic energy transfer, in which mesoscale eddies are generated through barotropic instability. The scale- and flow-dependent energy pathways diagnosed in this paper provide a basis for evaluating and developing scale- and flow-aware mesoscale eddy parameterizations. Significance Statement Blowing winds provide a major energy source for the large-scale ocean circulation. A substantial fraction of this energy is converted to smaller-scale eddies, which swirl through the ocean as sea cyclones. Ocean turbulence causes these eddies to transfer part of their energy back to the large-scale ocean currents. This ocean energy cycle is not fully simulated in numerical models, but it plays an important role in transporting heat, carbon, and nutrients throughout the world’s oceans. The purpose of this study is to quantify the ocean energy cycle by using fine-scale idealized numerical simulations of the Atlantic and Southern Oceans. Our results provide a basis for how to include unrepresented energy exchanges in coarse global climate models. 
    more » « less