Numerical investigation of porosity effect on a PCM's thermal performance in a porous rectangular channel via thermal lattice Boltzmann method

Mabrouk, Riheb; Naji, Hassane; Dhahri, Hacen; Hammouda, Seifeddine Ben; Younsi, Zohir

doi:10.1016/j.icheatmasstransfer.2020.104992

Cited by 16 publications

(10 citation statements)

References 46 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…These two findings are in agreement with the symbolic and numerical solutions, and they are resumed by Eqs. ( 70) and (71). These solutions are valid in the heterogeneous series subject to the plug flow u ¼ u n1n , but also in the stratified parabolic solutions with u ¼ u s1s , because MðyÞ is piece-wise constant and Eq.…”

Section: The Advective-diffusive Flux Fli=fmrmentioning

confidence: 95%

Mass-balance and locality versus accuracy with the new boundary and interface-conjugate approaches in advection-diffusion lattice Boltzmann method

Ginzburg

Silva

2021

Physics of Fluids

View full text Add to dashboard Cite

We introduce two new approaches, called A-LSOB and N-MR, for boundary and interface-conjugate conditions on flat or curved surface shapes in the advection-diffusion lattice Boltzmann method (LBM). The Local Second-Order, single-node A-LSOB enhances the existing Dirichlet and Neumann normal boundary treatments with respect to locality, accuracy, and P eclet parametrization. The normal-multi-reflection (N-MR) improves the directional flux schemes via a local release of their nonphysical tangential constraints. The A-LSOB and N-MR restore all first-and second-order derivatives from the nodal non-equilibrium solution, and they are conditioned to be exact on a piece-wise parabolic profile in a uniform arbitrary-oriented tangential velocity field. Additionally, the most compact and accurate single-node parabolic schemes for diffusion and flow in grid-inclined pipes are introduced. In simulations, the global mass-conservation solvability condition of the steady-state, two-relaxation-time (S-TRT) formulation is adjusted with either (i) a uniform mass-source or (ii) a corrective surface-flux. We conclude that (i) the surface-flux counterbalance is more accurate than the bulk one, (ii) the A-LSOB Dirichlet schemes are more accurate than the directional ones in the high P eclet regime, (iii) the directional Neumann advective-diffusive flux scheme shows the best conservation properties and then the best performance both in the tangential no-slip and interface-perpendicular flow, and (iv) the directional non-equilibrium diffusive flux extrapolation is the least conserving and accurate. The error P eclet dependency, Neumann invariance over an additive constant, and truncation isotropy guide this analysis. Our methodology extends from the d2q9 isotropic S-TRT to 3D anisotropic matrix collisions, Robin boundary condition, and the transient LBM.

show abstract

Section: The Advective-diffusive Flux Fli=fmrmentioning

confidence: 95%

Mass-balance and locality versus accuracy with the new boundary and interface-conjugate approaches in advection-diffusion lattice Boltzmann method

Ginzburg

Silva

2021

Physics of Fluids

View full text Add to dashboard Cite

show abstract

“…As the entropy generation reflected system irreversibility causing by heat transfer irreversibility (HTI) and fluid friction irreversibility (FFI), the local and average entropy generation rate under the LTNE assumption can be formulated as follows [16,17]:…”

Section: Entropy Generation Ratementioning

confidence: 99%

“…D and τ being the strain and stress tensors, respectively. The average entropy generation Ns ave per unit of volume can be calculated as [17]:…”

Section: Entropy Generation Ratementioning

confidence: 99%

“…Despite this weakness, the SRT approach remains the most popular and has shown its ability to handle phase change in porous media due to its simplicity, numerical precision and computational efficiency while mitigating the numerical diffusion [11,19] across phase interface. Thereby, the LBEs are expressed as [11,19]: (17) g s,i (x + e i δt, t + δt) − g s,i (x, t) streaming = − g s,i (x, t) − g eq s,i (x, t) /τ T,s collision term + (1 + δt∂ t /2)δtSr i,s source term (18) Computation 2022, 10, 3 6 of 14 where τ v (= 3ν + 0.5) and τ T,f;s are the dimensionless collision time for velocity and for temperatures [17,19]:…”

Section: Lb Equationsmentioning

confidence: 99%

“…δt, c (= δx/δt = 1; δx = δt) and α e,f,s being the lattice time step, the streaming speed and the effective diffusivity, respectively. Under the D2Q9 model (Figure 2), the local equilibrium distribution functions f eq i (x, t) and g eq i,f;s can be expressed as follows [4,11,17]:…”

Section: Lb Equationsmentioning

confidence: 99%

See 2 more Smart Citations

On Numerical Modeling of Thermal Performance Enhancementof a Heat Thermal Energy Storage System Using a Phase Change Material and a Porous Foam

et al. 2022

Self Cite

View full text Add to dashboard Cite

In this investigation, a comprehensive numerical analysis of the flow involved in an open-ended straight channel fully filled with a porous metal foam saturated and a phase change material (paraffin) has been performed using a single relaxation time lattice Boltzmann method (SRT-LBM) at the representative elementary volume (REV) scale. The enthalpy-based approach with three density functions has been employed to cope with the governing equations under the local thermal non-equilibrium (LTNE) condition. The in-house code has been validated through a comparison with a previous case in literature. The pore per inch density (10≤PPI≤60) and porosity (0.7≤ε≤0.9) effects of the metal structure were analyzed during melting/solidifying phenomena at two Reynolds numbers (Re = 200 and 400). The relevant findings are discussed for the LTNE intensity and the entropy generation rate (Ns). Through the simulations, the LTNE hypothesis turned out to be secure and valid. In addition, it is maximum for small PPI value (=10) whatever the parameters deemed. On the other hand, high porosity (=0.9) is advised to reduce the system’s irreversibility. However, at a moderate Re (=200), a small PPI (=10) would be appropriate to mitigate the system irreversibility during the charging case, while a large value (PPI = 60) might be advised for the discharging case. In this context, it can be stated that during the melting period, low porosity (=0.7) with low PPI (=10) improves thermal performance, reduces the system irreversibility and speeds up the melting rate, while for high porosity (=0.9), a moderate PPI (=30) should be used during the melting process to achieve an optimal system.

show abstract

Lattice Boltzmann method for heat transfer in phase change materials: a review

Kumar

Panda

Ghodke

et al. 2023

J Therm Anal Calorim

View full text Add to dashboard Cite

Numerical investigation of porosity effect on a PCM's thermal performance in a porous rectangular channel via thermal lattice Boltzmann method

Cited by 16 publications

References 46 publications

Mass-balance and locality versus accuracy with the new boundary and interface-conjugate approaches in advection-diffusion lattice Boltzmann method

Mass-balance and locality versus accuracy with the new boundary and interface-conjugate approaches in advection-diffusion lattice Boltzmann method

On Numerical Modeling of Thermal Performance Enhancementof a Heat Thermal Energy Storage System Using a Phase Change Material and a Porous Foam

Lattice Boltzmann method for heat transfer in phase change materials: a review

Contact Info

Product

Resources

About