Computational Fluid Dynamics for Phase Change Materials
Computational Fluid Dynamics (CFD) modeling and simulation for phase change materials (PCMs) is a technique to analyze the performance and behavior of PCMs. The CFD models have been successful in studying and analyzing the air quality, natural ventilation and stratified ventilation, air flow initiated by buoyancy forces and temperature space for the systems integrated with PCMs. Simple shapes like flat plates, cylinders or annular tubes, fins, macro- and micro-encapsulations with containers of different shape are often modeled in CFD software's to study. Typically the CFD models generally include Reynold’s Averaged  Navier-Stokes equation (RANS) modeling and Large Eddy Simulation (LES). Conservation equations of mass, momentum and energy (Navier – Stokes) are linearised, discretised, and applied to finite volumes to obtain a detailed solution for field distributions of air pressure, velocity and temperature for both indoor spaces integrated with PCMs. Governing Equations. Mass Equation. where Energy Equation. where Navier Stokes equation. formula_3 Here f represents "other" body forces (per unit volume), such as gravity or centrifugal force. The shear stress term formula_4 becomes formula_5, where formula_6 is the vector Laplacian. Boussinesq eddy-viscosity approximation. where Assumptions. commonly used assumptions are Phase Change  Model. Two main thermal characteristics of phase change are the enthalpy-temperature relationship and temperature hysteresis. PCMs tend to have varying enthalpy temperature relationships due to the fact that they are blends of different materials, but pure PCMs have a more localized relationship, which can be approximated by single values for the enthalpy and phase change temperature. Hysteresis is the phenomenon which causes the PCM to melt and freezes in different temperature ranges and with different enthalpies, which results in a different temperature-enthalpy curve for melting and freezing. Hysteresis is related to the chemical and kinetic properties of the  material. The commonly used enthalpy-porosity model in commercial CFD codes assumes, a linear enthalpy-temperature relationship and ignores hysteresis. The alternate is to use enthalpy-porosity method. When used to simulate PCM sails and a PCM plate-fin unit it produce reasonable temperature prediction in global space temperature terms. However there are inaccuracies in transient simulations where time dependent PCM and local wall and air temperatures are of interest. This is over come by use of source terms that considers hysteresis and varying enthalpy-temperature relationship.  CFD-DEM model are also used sometimes. Phase motion of discrete solids  or particles is obtained by the Discrete Element Method (DEM) which applies Newton's laws of motion to every particle and the flow of continuum fluid is described by the local averaged Navier–Stokes equations that can be solved by the traditional Computational Fluid Dynamics (CFD).CFDEMcoupling (DCS Computing GmbH) is one such open source toolbox for CFD-DEM coupling. Process. The Governing equations are discretized using an explicit Finite Volume Method. The velocity-pressure coupling is resolved by adopting a Fractional Step Method. The adoption of the enthalpy method allows working with a fixed grid instead of an interface tracking  method. The momentum source term intended to model the presence of solid is only needed in the control volumes that contain solid and liquid, not in the pure solid containing volumes. The final form of the source term coefficient(S)depends on the approximation adopted for the behavior of the flow in the “mushy zone” (where mixed solid and liquid states are present). However, in the case of constant phase change temperature, the solid-liquid interface should be of infinitesimal width (although it cannot be thinner than one control volume width in our simulations); therefore, the formulation used for the source  term is not very important in a physical sense, as long as it manages to bring the velocity to zero in mostly solid control volumes and to vanish if the volume contains pure liquid. The various CFD codes[1-3] has been employed for the modeling and simulation of the PCM system to understand the heat transfer mechanism, solidification and melting process, distribution of temperature profile and prediction of the air flow. Various commercial packages have been coupled with the CFD analysis to appreciate the feasibility of evaluating the behavior of PCM integrated system. The simulated heat transfer behavior of the  PCM in Mobilized Thermal Energy Storage, during the charging process can be successfully conducted by CFD modeling  The Volume-Of-Fluid method is employed to solve for the temperature distribution in the multiphase, 2-dimensional pressure-based model. It accounts for the heat transfer mechanism, melting time, and the influence of the structure in charging process using Fluent 12.1. The governing equations employed are mass conversion and continuity equations. Integral, quasi-1D calculations have been reported  mainly for conduction-dominated problem using CFD simulation. It was reported that out of three geometries (cubic, cylindrical and spherical), the spherical capsule  will have the maximum heat for the heat transfer fluid. Also it is concluded that salt hydrates based PCMs are the better choice over organic PCMs. Applications. The systems are developed in such manner that phase change materials are in the shell portion of the module and passage for the flow of air through the tubes. Conjugate steady state CFD heat transfer analysis has been carried out  to analyze the flow and temperature variation of heat transfer fluid in the system. It paves the way for selection and assessment of the geometrical and flow parameters, PCM solidification characteristics for  the given boundary conditions The comparative analysis, to further enhance the effectiveness of shell and tube PCMs has also been accomplished via CFD analysis. Various CFD models with different configuration such as pins embedded on a tube with heat transfer fluid (HTF) flowing in it, with PCM surrounding the tube, fins embedded instead of pins and different configurations of fins on the tube are analyzed, by employing ANSYS code.