Search

Mimicking vascular systems in living beings, designers have realized microvascular composites to achieve thermal regulation and other functionalities, such as electromagnetic modulation, sensing, and healing. Such material systems avail circulating fluids through embedded vasculatures to accomplish the mentioned functionalities that benefit various aerospace, military, and civilian applications. Although heat transfer is a mature field, control of thermal characteristics in synthetic microvascular systems via circulating fluids is new, and a theoretical underpinning is lacking. What will benefit designers are predictive mathematical models and an in-depth qualitative understanding of vascular-based active cooling/heating. So, the central focus of this paper is to address the remarked knowledge gap. \emph{First}, we present a reduced-order model with broad applicability, allowing the inlet temperature to differ from the ambient temperature. \emph{Second}, we apply mathematical analysis tools to this reduced-order model and reveal many heat transfer properties of fluid-sequestered vascular systems. We derive point-wise properties (minimum, maximum, and comparison principles) and global properties (e.g., bounds on performance metrics such as the mean surface temperature and thermal efficiency). These newfound results deepen our understanding of active cooling/heating and propel the perfecting of thermal regulation systems.

We derive a numerical method for Darcy flow, hence also for Poisson's equation in mixed (first order) form, based on discrete exterior calculus (DEC). Exterior calculus is a generalization of vector calculus to smooth manifolds and DEC is one of its discretizations on simplicial complexes such as triangle and tetrahedral meshes. DEC is a coordinate invariant discretization, in that it does not depend on the embedding of the simplices or the whole mesh. We start by rewriting the governing equations of Darcy flow using the language of exterior calculus. This yields a formulation in terms of flux differential form and pressure. The numerical method is then derived by using the framework provided by DEC for discretizing differential forms and operators that act on forms. We also develop a discretization for spatially dependent Hodge star that varies with the permeability of the medium. This also allows us to address discontinuous permeability. The matrix representation for our discrete non-homogeneous Hodge star is diagonal, with positive diagonal entries. The resulting linear system of equations for flux and pressure are saddle type, with a diagonal matrix as the top left block. The performance of the proposed numerical method is illustrated on many standard test problems. These include patch tests in two and three dimensions, comparison with analytically known solution in two dimensions, layered medium with alternating permeability values, and a test with a change in permeability along the flow direction. We also show numerical evidence of convergence of the flux and the pressure. A convergence experiment is also included for Darcy flow on a surface. A short introduction to the relevant parts of smooth and discrete exterior calculus is included in this paper. We also include a discussion of the boundary condition in terms of exterior calculus., Comment: Added numerical experiment for flow on a surface. Other small changes in meshing related comments