Computational Fluid Dynamics for Phase Change Materials

Last updated

Computational Fluid Dynamics (CFD) modeling and simulation for phase change materials (PCMs) is a technique used 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 shapes are often modeled in CFD software's to study.

Contents

Typically the CFD models 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

Here f represents "other" body forces (per unit volume), such as gravity or centrifugal force. The shear stress term becomes , where 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.[8]

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. [9][10]

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.[11]

Applications

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 [4] 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 [5] 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.

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 [6] 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[7]. 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.

Related Research Articles

Continuum mechanics is a branch of mechanics that deals with the deformation of and transmission of forces through materials modeled as a continuous medium rather than as discrete particles.

<span class="mw-page-title-main">Fluid dynamics</span> Aspects of fluid mechanics involving flow

In physics, physical chemistry and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluids – liquids and gases. It has several subdisciplines, including aerodynamics and hydrodynamics. Fluid dynamics has a wide range of applications, including calculating forces and moments on aircraft, determining the mass flow rate of petroleum through pipelines, predicting weather patterns, understanding nebulae in interstellar space and modelling fission weapon detonation.

The thermal conductivity of a material is a measure of its ability to conduct heat. It is commonly denoted by , , or and is measured in W·m−1·K−1.

Thermal conduction is the diffusion of thermal energy (heat) within one material or between materials in contact. The higher temperature object has molecules with more kinetic energy; collisions between molecules distributes this kinetic energy until an object has the same kinetic energy throughout. Thermal conductivity, frequently represented by k, is a property that relates the rate of heat loss per unit area of a material to its rate of change of temperature. Essentially, it is a value that accounts for any property of the material that could change the way it conducts heat. Heat spontaneously flows along a temperature gradient. For example, heat is conducted from the hotplate of an electric stove to the bottom of a saucepan in contact with it. In the absence of an opposing external driving energy source, within a body or between bodies, temperature differences decay over time, and thermal equilibrium is approached, temperature becoming more uniform.

<span class="mw-page-title-main">Boundary layer</span> Layer of fluid in the immediate vicinity of a bounding surface

In physics and fluid mechanics, a boundary layer is the thin layer of fluid in the immediate vicinity of a bounding surface formed by the fluid flowing along the surface. The fluid's interaction with the wall induces a no-slip boundary condition. The flow velocity then monotonically increases above the surface until it returns to the bulk flow velocity. The thin layer consisting of fluid whose velocity has not yet returned to the bulk flow velocity is called the velocity boundary layer.

Quantum turbulence is the name given to the turbulent flow – the chaotic motion of a fluid at high flow rates – of quantum fluids, such as superfluids. The idea that a form of turbulence might be possible in a superfluid via the quantized vortex lines was first suggested by Richard Feynman. The dynamics of quantum fluids are governed by quantum mechanics, rather than classical physics which govern classical (ordinary) fluids. Some examples of quantum fluids include superfluid helium, Bose–Einstein condensates (BECs), polariton condensates, and nuclear pasta theorized to exist inside neutron stars. Quantum fluids exist at temperatures below the critical temperature at which Bose-Einstein condensation takes place.

<span class="mw-page-title-main">Euler equations (fluid dynamics)</span> Set of quasilinear hyperbolic equations governing adiabatic and inviscid flow

In fluid dynamics, the Euler equations are a set of partial differential equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. In particular, they correspond to the Navier–Stokes equations with zero viscosity and zero thermal conductivity.

<span class="mw-page-title-main">Onsager reciprocal relations</span> Relations between flows and forces, or gradients, in thermodynamic systems

In thermodynamics, the Onsager reciprocal relations express the equality of certain ratios between flows and forces in thermodynamic systems out of equilibrium, but where a notion of local equilibrium exists.

<span class="mw-page-title-main">Smoothed-particle hydrodynamics</span> Method of hydrodynamics simulation

Smoothed-particle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. It was developed by Gingold and Monaghan and Lucy in 1977, initially for astrophysical problems. It has been used in many fields of research, including astrophysics, ballistics, volcanology, and oceanography. It is a meshfree Lagrangian method, and the resolution of the method can easily be adjusted with respect to variables such as density.

<span class="mw-page-title-main">Lattice Boltzmann methods</span> Class of computational fluid dynamics methods

The lattice Boltzmann methods (LBM), originated from the lattice gas automata (LGA) method (Hardy-Pomeau-Pazzis and Frisch-Hasslacher-Pomeau models), is a class of computational fluid dynamics (CFD) methods for fluid simulation. Instead of solving the Navier–Stokes equations directly, a fluid density on a lattice is simulated with streaming and collision (relaxation) processes. The method is versatile as the model fluid can straightforwardly be made to mimic common fluid behaviour like vapour/liquid coexistence, and so fluid systems such as liquid droplets can be simulated. Also, fluids in complex environments such as porous media can be straightforwardly simulated, whereas with complex boundaries other CFD methods can be hard to work with.

In fluid dynamics, turbulence kinetic energy (TKE) is the mean kinetic energy per unit mass associated with eddies in turbulent flow. Physically, the turbulence kinetic energy is characterized by measured root-mean-square (RMS) velocity fluctuations. In the Reynolds-averaged Navier Stokes equations, the turbulence kinetic energy can be calculated based on the closure method, i.e. a turbulence model.

The derivation of the Navier–Stokes equations as well as their application and formulation for different families of fluids, is an important exercise in fluid dynamics with applications in mechanical engineering, physics, chemistry, heat transfer, and electrical engineering. A proof explaining the properties and bounds of the equations, such as Navier–Stokes existence and smoothness, is one of the important unsolved problems in mathematics.

In fluid dynamics, Airy wave theory gives a linearised description of the propagation of gravity waves on the surface of a homogeneous fluid layer. The theory assumes that the fluid layer has a uniform mean depth, and that the fluid flow is inviscid, incompressible and irrotational. This theory was first published, in correct form, by George Biddell Airy in the 19th century.

Relativistic heat conduction refers to the modelling of heat conduction in a way compatible with special relativity. In special relativity, the usual heat equation for non-relativistic heat conduction must be modified, as it leads to faster-than-light signal propagation. Relativistic heat conduction, therefore, encompasses a set of models for heat propagation in continuous media that are consistent with relativistic causality, namely the principle that an effect must be within the light-cone associated to its cause. Any reasonable relativistic model for heat conduction must also be stable, in the sense that differences in temperature propagate both slower than light and are damped over time.

<span class="mw-page-title-main">Enthalpy of fusion</span> Enthalpy change when a substance melts

In thermodynamics, the enthalpy of fusion of a substance, also known as (latent) heat of fusion, is the change in its enthalpy resulting from providing energy, typically heat, to a specific quantity of the substance to change its state from a solid to a liquid, at constant pressure.

In engineering, physics, and chemistry, the study of transport phenomena concerns the exchange of mass, energy, charge, momentum and angular momentum between observed and studied systems. While it draws from fields as diverse as continuum mechanics and thermodynamics, it places a heavy emphasis on the commonalities between the topics covered. Mass, momentum, and heat transport all share a very similar mathematical framework, and the parallels between them are exploited in the study of transport phenomena to draw deep mathematical connections that often provide very useful tools in the analysis of one field that are directly derived from the others.

The multiphase particle-in-cell method (MP-PIC) is a numerical method for modeling particle-fluid and particle-particle interactions in a computational fluid dynamics (CFD) calculation. The MP-PIC method achieves greater stability than its particle-in-cell predecessor by simultaneously treating the solid particles as computational particles and as a continuum. In the MP-PIC approach, the particle properties are mapped from the Lagrangian coordinates to an Eulerian grid through the use of interpolation functions. After evaluation of the continuum derivative terms, the particle properties are mapped back to the individual particles. This method has proven to be stable in dense particle flows, computationally efficient, and physically accurate. This has allowed the MP-PIC method to be used as particle-flow solver for the simulation of industrial-scale chemical processes involving particle-fluid flows.

In fluid mechanics, non-dimensionalization of the Navier–Stokes equations is the conversion of the Navier–Stokes equation to a nondimensional form. This technique can ease the analysis of the problem at hand, and reduce the number of free parameters. Small or large sizes of certain dimensionless parameters indicate the importance of certain terms in the equations for the studied flow. This may provide possibilities to neglect terms in the considered flow. Further, non-dimensionalized Navier–Stokes equations can be beneficial if one is posed with similar physical situations – that is problems where the only changes are those of the basic dimensions of the system.

Heat transfer physics describes the kinetics of energy storage, transport, and energy transformation by principal energy carriers: phonons, electrons, fluid particles, and photons. Heat is thermal energy stored in temperature-dependent motion of particles including electrons, atomic nuclei, individual atoms, and molecules. Heat is transferred to and from matter by the principal energy carriers. The state of energy stored within matter, or transported by the carriers, is described by a combination of classical and quantum statistical mechanics. The energy is different made (converted) among various carriers. The heat transfer processes are governed by the rates at which various related physical phenomena occur, such as the rate of particle collisions in classical mechanics. These various states and kinetics determine the heat transfer, i.e., the net rate of energy storage or transport. Governing these process from the atomic level to macroscale are the laws of thermodynamics, including conservation of energy.

CFD stands for computational fluid dynamics. As per this technique, the governing differential equations of a flow system or thermal system are known in the form of Navier–Stokes equations, thermal energy equation and species equation with an appropriate equation of state. In the past few years, CFD has been playing an increasingly important role in building design, following its continuing development for over a quarter of a century. The information provided by CFD can be used to analyse the impact of building exhausts to the environment, to predict smoke and fire risks in buildings, to quantify indoor environment quality, and to design natural ventilation systems.

References