Trajectory integrated smoothening of exchange fields for discrete phase simulations

A.S. Simonsen (Corresponding author), N.G. Deen, K. Sørensen, T. Condra, L. Rosendahl

Research output: Contribution to journalArticleAcademicpeer-review

Abstract

Eulerian-Lagrangian simulations often require a two-way coupling, where the discrete phase affects the continuous phase and vice versa. The continuous phase in a CFD-simulation is typically represented on a computational mesh, which consists of a number of discrete cells. As the discrete phase is not directly associated with this mesh, a mapping from the discrete phase to the continuous phase is required in order to couple the two phases. This will distribute the exchange rates from the discrete to the continuous phase, where the exchanged properties can be thermal energy, momentum, species etc. The mapping can be done in several ways, and different methods exist. Some of these include a smoothening step, which ensures that the exchange fields do not contain large spikes between any two neighbouring cells, which typically increases the convergence rate of the CFD simulation. Gaussian smoothening is often applied, where the discrete phase exchange rates are spread to the surrounding cells by a normal distribution with a specified standard deviation. This paper proposes a new weight for distributing the discrete phase exchange rates to the continuous phase cells. This is done by analytically integrating any distribution over the entirety of the particle trajectory, and only requires the probability- and cumulative density functions to be specified. The proposed weight is used in combination with three different smoothening methods: node-based, kernel-based, and cell-based. Common existing weight distributions are also investigated, which are used in combination with the different methods. Each method and weight combination is compared with the analytical ground-truth exchange field, where the distributions are accurately integrated over the cell volumes. The results show that the proposed weight has an average error of ≈ 2 % in the exchange field for the test case used. This is about three times more accurate compared to existing weights, where the exchange field is not integrated along the trajectory of the particle. Alternative distributions were investigated to decrease the computational requirements. It was found that a fourth-order polynomial distribution was about 20 times faster to evaluate, whilst the error induced by switching to this distribution was ≈ 4.5 % compared with the previous ≈ 2 %.

LanguageEnglish
Pages15-23
Number of pages9
JournalComputers & Fluids
Volume186
DOIs
StatePublished - 30 May 2019

Fingerprint

Computational fluid dynamics
Trajectories
Normal distribution
Thermal energy
Probability density function
Momentum
Polynomials

Keywords

  • Gaussian
  • Mapping
  • Smoothening
  • Weighting

Cite this

Simonsen, A.S. ; Deen, N.G. ; Sørensen, K. ; Condra, T. ; Rosendahl, L./ Trajectory integrated smoothening of exchange fields for discrete phase simulations. In: Computers & Fluids. 2019 ; Vol. 186. pp. 15-23
@article{37611fca21dc46129d71ef1b1eeece12,
title = "Trajectory integrated smoothening of exchange fields for discrete phase simulations",
abstract = "Eulerian-Lagrangian simulations often require a two-way coupling, where the discrete phase affects the continuous phase and vice versa. The continuous phase in a CFD-simulation is typically represented on a computational mesh, which consists of a number of discrete cells. As the discrete phase is not directly associated with this mesh, a mapping from the discrete phase to the continuous phase is required in order to couple the two phases. This will distribute the exchange rates from the discrete to the continuous phase, where the exchanged properties can be thermal energy, momentum, species etc. The mapping can be done in several ways, and different methods exist. Some of these include a smoothening step, which ensures that the exchange fields do not contain large spikes between any two neighbouring cells, which typically increases the convergence rate of the CFD simulation. Gaussian smoothening is often applied, where the discrete phase exchange rates are spread to the surrounding cells by a normal distribution with a specified standard deviation. This paper proposes a new weight for distributing the discrete phase exchange rates to the continuous phase cells. This is done by analytically integrating any distribution over the entirety of the particle trajectory, and only requires the probability- and cumulative density functions to be specified. The proposed weight is used in combination with three different smoothening methods: node-based, kernel-based, and cell-based. Common existing weight distributions are also investigated, which are used in combination with the different methods. Each method and weight combination is compared with the analytical ground-truth exchange field, where the distributions are accurately integrated over the cell volumes. The results show that the proposed weight has an average error of ≈ 2 {\%} in the exchange field for the test case used. This is about three times more accurate compared to existing weights, where the exchange field is not integrated along the trajectory of the particle. Alternative distributions were investigated to decrease the computational requirements. It was found that a fourth-order polynomial distribution was about 20 times faster to evaluate, whilst the error induced by switching to this distribution was ≈ 4.5 {\%} compared with the previous ≈ 2 {\%}.",
keywords = "Gaussian, Mapping, Smoothening, Weighting",
author = "A.S. Simonsen and N.G. Deen and K. S{\o}rensen and T. Condra and L. Rosendahl",
year = "2019",
month = "5",
day = "30",
doi = "10.1016/j.compfluid.2019.03.015",
language = "English",
volume = "186",
pages = "15--23",
journal = "Computers & Fluids",
issn = "0045-7930",
publisher = "Elsevier",

}

Trajectory integrated smoothening of exchange fields for discrete phase simulations. / Simonsen, A.S. (Corresponding author); Deen, N.G.; Sørensen, K.; Condra, T.; Rosendahl, L.

In: Computers & Fluids, Vol. 186, 30.05.2019, p. 15-23.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Trajectory integrated smoothening of exchange fields for discrete phase simulations

AU - Simonsen,A.S.

AU - Deen,N.G.

AU - Sørensen,K.

AU - Condra,T.

AU - Rosendahl,L.

PY - 2019/5/30

Y1 - 2019/5/30

N2 - Eulerian-Lagrangian simulations often require a two-way coupling, where the discrete phase affects the continuous phase and vice versa. The continuous phase in a CFD-simulation is typically represented on a computational mesh, which consists of a number of discrete cells. As the discrete phase is not directly associated with this mesh, a mapping from the discrete phase to the continuous phase is required in order to couple the two phases. This will distribute the exchange rates from the discrete to the continuous phase, where the exchanged properties can be thermal energy, momentum, species etc. The mapping can be done in several ways, and different methods exist. Some of these include a smoothening step, which ensures that the exchange fields do not contain large spikes between any two neighbouring cells, which typically increases the convergence rate of the CFD simulation. Gaussian smoothening is often applied, where the discrete phase exchange rates are spread to the surrounding cells by a normal distribution with a specified standard deviation. This paper proposes a new weight for distributing the discrete phase exchange rates to the continuous phase cells. This is done by analytically integrating any distribution over the entirety of the particle trajectory, and only requires the probability- and cumulative density functions to be specified. The proposed weight is used in combination with three different smoothening methods: node-based, kernel-based, and cell-based. Common existing weight distributions are also investigated, which are used in combination with the different methods. Each method and weight combination is compared with the analytical ground-truth exchange field, where the distributions are accurately integrated over the cell volumes. The results show that the proposed weight has an average error of ≈ 2 % in the exchange field for the test case used. This is about three times more accurate compared to existing weights, where the exchange field is not integrated along the trajectory of the particle. Alternative distributions were investigated to decrease the computational requirements. It was found that a fourth-order polynomial distribution was about 20 times faster to evaluate, whilst the error induced by switching to this distribution was ≈ 4.5 % compared with the previous ≈ 2 %.

AB - Eulerian-Lagrangian simulations often require a two-way coupling, where the discrete phase affects the continuous phase and vice versa. The continuous phase in a CFD-simulation is typically represented on a computational mesh, which consists of a number of discrete cells. As the discrete phase is not directly associated with this mesh, a mapping from the discrete phase to the continuous phase is required in order to couple the two phases. This will distribute the exchange rates from the discrete to the continuous phase, where the exchanged properties can be thermal energy, momentum, species etc. The mapping can be done in several ways, and different methods exist. Some of these include a smoothening step, which ensures that the exchange fields do not contain large spikes between any two neighbouring cells, which typically increases the convergence rate of the CFD simulation. Gaussian smoothening is often applied, where the discrete phase exchange rates are spread to the surrounding cells by a normal distribution with a specified standard deviation. This paper proposes a new weight for distributing the discrete phase exchange rates to the continuous phase cells. This is done by analytically integrating any distribution over the entirety of the particle trajectory, and only requires the probability- and cumulative density functions to be specified. The proposed weight is used in combination with three different smoothening methods: node-based, kernel-based, and cell-based. Common existing weight distributions are also investigated, which are used in combination with the different methods. Each method and weight combination is compared with the analytical ground-truth exchange field, where the distributions are accurately integrated over the cell volumes. The results show that the proposed weight has an average error of ≈ 2 % in the exchange field for the test case used. This is about three times more accurate compared to existing weights, where the exchange field is not integrated along the trajectory of the particle. Alternative distributions were investigated to decrease the computational requirements. It was found that a fourth-order polynomial distribution was about 20 times faster to evaluate, whilst the error induced by switching to this distribution was ≈ 4.5 % compared with the previous ≈ 2 %.

KW - Gaussian

KW - Mapping

KW - Smoothening

KW - Weighting

UR - http://www.scopus.com/inward/record.url?scp=85063901181&partnerID=8YFLogxK

U2 - 10.1016/j.compfluid.2019.03.015

DO - 10.1016/j.compfluid.2019.03.015

M3 - Article

VL - 186

SP - 15

EP - 23

JO - Computers & Fluids

T2 - Computers & Fluids

JF - Computers & Fluids

SN - 0045-7930

ER -

Simonsen AS, Deen NG, Sørensen K, Condra T, Rosendahl L. Trajectory integrated smoothening of exchange fields for discrete phase simulations. Computers & Fluids. 2019 May 30;186:15-23. Available from, DOI: 10.1016/j.compfluid.2019.03.015