Non-Newtonian fluid - Revision history

E6WikiAdmin at 19:31, 22 October 2022

2022-10-22T19:31:04Z

E6WikiAdmin: /* de Vahl Davis example */

2022-10-22T19:30:00Z

‎de Vahl Davis example

E6WikiAdmin: /* de Vahl Davis example */

2022-10-22T19:29:28Z

‎de Vahl Davis example

E6WikiAdmin at 19:28, 22 October 2022

2022-10-22T19:28:19Z

MihaR: Improved grammar

2020-10-21T11:29:38Z

Improved grammar

MihaR: Link to example

2020-10-12T18:32:54Z

Link to example

MihaR: Added tags to figure captions

2020-10-12T18:13:54Z

Added <caption> tags to figure captions

MihaR: Animation display fix

2020-10-08T19:42:27Z

Animation display fix

MihaR: Initial version

2020-10-08T19:37:35Z

Initial version

New page

Click here to return back to [[Fluid Mechanics]]

=Introduction=
Other fluid mechanics examples have assumed that the fluids in question are Newtonian, but that is not always the case when dealing with real world applications. Many significant fluids exhibit non-Newtonian behaviour with examples including blood, foodstuffs, polymer melts... Newtonian fluids have a linear relationship between viscous stress and shear rate but that is no longer true for non-Newtonian fluids.

Some fluids exhibit increasing or decreasing viscosity with increasing shear rate which is the simplest form of non-Newtonian behaviour that we will tackle in the following sections. Some act as solids at small stress but begin flowing when sufficiently disturbed, while others behaviour depends on the duration of movement and applied stress.

=Model formulation=
Viscous forces are a result of velocity differentials within the fluid and are therefore connected with the velocity gradient. We are interested in the symmetric part that excludes the rotational contributions. Shear rate tensor $\dot{u}$ for incompressible fluids can be expressed as
$$
\dot{u}=\frac{1}{2}\left(\boldsymbol{\nabla} \boldsymbol{v}+(\boldsymbol{\nabla} \boldsymbol{v})^{\top}\right),
$$
where $\boldsymbol{v}$ denotes the fluid velocity. We introduce viscosity $\mu$ through a generalized Newton's law of viscosity with viscous stress tensor $\tau$ expressed as
$$
\tau = \mu \dot{u}.
$$
In general $\mu$ is an arbitrary tensor but usually reduces to a scalar value due to isotropic properties of fluids. Non-Newtonian fluids are fluids where viscosity is a function of shear stress and duration of stress. There are many, often quantitative, models of varying complexities used to describe the viscosity of various non-Newtonian fluids.

We limit ourselves to Ostwald-de Waele power law model
$$
\mu = \mu_0 \left\| \dot{u} \right\| ^{n-1} = \mu_0 \left( 2 \dot{u}{:}\dot{u} \right) ^\frac{n-1}{2},
$$

<figure id="fig:devahl_scheme">
[[File:image.png|thumb|Scheme of the de Vahl Davis benchmark test.]]
</figure>

which is one of the simplest non-Newtonian viscosity models. We introduced two parameters used to fit the model to the experimental data for the fluid in question. Viscosity constant $\mu_0$ used as a scaling factor and non-Newtonian power index $n$ that describes the behaviour. Model with $n < 1$ describes shear-thinning (pseudoplastic), $n = 1$ Newtonian and $n > 1$ shear thickening (dilatant) fluids.

We will focus on shear thinning fluids as their behaviour of reduced viscosity at higher velocity gradients leads towards interesting behaviour with decreasing $n$.

=de Vahl Davis example=
We take a closer look at the behaviour of non-Newtonian fluid in a natural convection driven flow within a square cavity. The convection is driven by differentially heated east and west walls with the north and south walls thermally insulated as is schematically shown in <xr id="fig:devahl_scheme"/>.

This is a common benchmark with existing solutions that we can benchmark our meshless implementation against. Convergence and comparison with benchmark article is shown in <xr id="fig:nonNewtonianConvergence"/> where we plot the changes of average Nusselt number with increasing node density. Nusselt number is a dimensionless number that expresses the ratio between convective and conductive heat transfer. We use the average value on the western wall and define it as

<figure id="fig:nonNewtonianConvergence">
[[File:nonNewtonianConvergence.png|thumb|left|700px|
Convergence of average Nusselt number value with increasing node count and comparison with reference values from reference<ref name="Turan2011">O. Turan, A. Sachdeva, N. Chakraborty in R. J. Poole, Laminar natural convection of power-law Fluids in a square enclosure with differentially heated side walls subjected to constant temperatures, Journal of Non-Newtonian Fluid Mechanics 166, 1049 (2011).</ref>. Convergence is calculated at Ra$=10^4$ and Pr$=10$.]]
</figure>

$$
\mathrm{Nu}=\frac{\Omega_W}{T_{H}-T_{C}}\left|\frac{\partial T}{\partial p_x}\right|_{p_x=0},
$$
where $T$ denotest the demperature, $p_x$ $x$ coordinate and $\Omega_W$ cavity width. Nusselt number is a convenient tool that can be used to reduce the dynamics of a complex multidimensional system into a single scalar value.

We define the problem in terms of Rayleigh and Prandtl dimensionless numbers that need to be slightly modified to remain relevant in non-Newtonian case but specifics are not really relevant for this review. Dimensionless number definitions used in these calculations match reference<ref name="Turan2011">O. Turan, A. Sachdeva, N. Chakraborty in R. J. Poole, Laminar natural convection of power-law Fluids in a square enclosure with differentially heated side walls subjected to constant temperatures, Journal of Non-Newtonian Fluid Mechanics 166, 1049 (2011).</ref> and can be viewed there.

Parameters that lead to oscillatory behaviour in de Vahl Davis cavity can be found in literature<ref name="Kosec2011">G. Kosec in B. Šarler, Solution of a low Prandtl number natural convection benchmark by a local meshless method, International Journal of Numerical Methods for Heat & Fluid Flow 23, 22 (2013).</ref>. Referenced article found oscillatory behaviour in Newtonian fluids at Ra=$5 \cdot 10^4$ and Pr=$0.01$ but we want to start slightly lower to find the transition as the dynamics become wilder with smaller $n$. We chose Ra=$2 \cdot 10^4$ and Pr=$0.01$ as the parameters to be used for all subsequently displayed calculations.

The following figures display or mention dimensionless time that is defined with the same definition as described in reference<ref name="Kosec2011">G. Kosec in B. Šarler, Solution of a low Prandtl number natural convection benchmark by a local meshless method, International Journal of Numerical Methods for Heat & Fluid Flow 23, 22 (2013).</ref>.

<figure id="fig:nusseltEvolutionNN">
[[File:nusseltEvolutionNN.png|thumb|left|700px|
Evolution of average Nusselt number in a differentially heated square cavity for different values of non-Newtonian power index $n$.]]
</figure>

We reuse the Nusselt number to analyse the differences in behaviour for different non-Newtonian power indices $n$ shown in <xr id="fig:nusseltEvolutionNN"/>. The Newtonian case with $n=1$ settles into a stationary case which is expected, as we chose parameters below the previously identified oscillatory solution. Evolutions of Nusselt number become oscillatory as we increase the strenght of non-Newtonian behaviour while the solutions become less and less regular as we reduce the power index $n$ which can be seen in pictures of velocity and temperature distributions shown in <xr id="fig:squareCavityDisplay"/>. Detailed animation of behaviour in the most extreme $n=0.7$ case is shown in <xr id="fig:animation"/>.

Exact implementation of non-Newtonian fluid modeled with Ostwald-de Waele power law can be found in example.
TODO: link to example

<figure id="fig:squareCavityDisplay">
[[File:squareCavity_t10_Ra20000Pr001.png|frame|center|
Velocity and temperature fields within the cavity at $t=10$ for different values of non-Newtonian power index $n$.]]
</figure>

<figure id="fig:animation">
[[File:squareCavityAnimationNN_n07.mp4|frame|center|
Animation of velocity and temperature fields within the cavity between $t=10$ and $t=11$ for $n=0.7$.]]
</figure>

=References=

<references />

← Older revision		Revision as of 19:31, 22 October 2022
Line 39:		Line 39:
	where $\mu_0$ denotes the viscosity constant, $\boldsymbol{v}$ the fluid velocity and $n$ the non-Newtonian power index. The model parameters $\mu_0$ and $n$ for specific fluids are usually determined from experimental data.		where $\mu_0$ denotes the viscosity constant, $\boldsymbol{v}$ the fluid velocity and $n$ the non-Newtonian power index. The model parameters $\mu_0$ and $n$ for specific fluids are usually determined from experimental data.

−	~~<figure id="fig:devahl_scheme">~~
−	~~[[File:image.png\|thumb\|<caption>Scheme of the de Vahl Davis benchmark test.</caption>]]~~
−	~~</figure>~~

	We will solve the natural convection driven problem in a square cavity that is schematically shown in <xr id="fig:devahl_scheme"/>. The north and south walls are thermally insulated while the east and west walls are kept at constant but different temperatures. Different temperatures lead to differing fluid densities that cause natural convection. The problem is defined with a system of coupled partial differential equations		We will solve the natural convection driven problem in a square cavity that is schematically shown in <xr id="fig:devahl_scheme"/>. The north and south walls are thermally insulated while the east and west walls are kept at constant but different temperatures. Different temperatures lead to differing fluid densities that cause natural convection. The problem is defined with a system of coupled partial differential equations

@@ Line 32: / Line 32: @@
 =de Vahl Davis example=
 We take a closer look at the behaviour of non-Newtonian fluid in a natural convection driven flow within a square cavity. The convection is driven by differentially heated east and west walls with the north and south walls thermally insulated as is schematically shown in <xr id="fig:devahl_scheme"/>.
 This example is similar to the explicit pressure correction described in the [[de Vahl Davis natural convection test]] with the addition of the non-Newtonian viscosity $\mu$ modeled with the Ostwald-de Waele power law
@@ Line 124: / Line 86: @@
 The code is available as [https://gitlab.com/e62Lab/medusa/blob/dev/examples/nonNewtonian_fluid/nonNewtonian_fluid.cpp nonNewtonian_fluid.cpp] in the Medusa examples.
-Go back to [[Medusa#Examples|Examples]].
+This is a common benchmark with existing solutions that we can compare our meshless implementation against. Convergence and comparison with the benchmark article is shown in <xr id="fig:nonNewtonianConvergence"/> where we plot how the average Nusselt number changes with increasing node density. Nusselt number is a dimensionless number that expresses the ratio between the convective and conductive heat transfer. We use the average value on the western wall and define it as
 =References=
 <references />

@@ Line 30: / Line 30: @@
 We will focus on shear thinning fluids as their behaviour of reduced viscosity at higher velocity gradients leads towards interesting behaviour with decreasing $n$.
-=de Vahl Davis example=
+= non-Newtonian de Vahl Davis case=
 We take a closer look at the behaviour of non-Newtonian fluid in a natural convection driven flow within a square cavity. The convection is driven by differentially heated east and west walls with the north and south walls thermally insulated as is schematically shown in <xr id="fig:devahl_scheme"/>.

@@ Line 33: / Line 33: @@
 We take a closer look at the behaviour of non-Newtonian fluid in a natural convection driven flow within a square cavity. The convection is driven by differentially heated east and west walls with the north and south walls thermally insulated as is schematically shown in <xr id="fig:devahl_scheme"/>.
-This is a common benchmark with existing solutions that we can benchmark our meshless implementation against. Convergence and comparison with benchmark article is shown in <xr id="fig:nonNewtonianConvergence"/> where we plot the changes of average Nusselt number with increasing node density. Nusselt number is a dimensionless number that expresses the ratio between convective and conductive heat transfer. We use the average value on the western wall and define it as
+This is a common benchmark with existing solutions that we can compare our meshless implementation against. Convergence and comparison with the benchmark article is shown in <xr id="fig:nonNewtonianConvergence"/> where we plot how the average Nusselt number changes with increasing node density. Nusselt number is a dimensionless number that expresses the ratio between the convective and conductive heat transfer. We use the average value on the western wall and define it as
 <figure id="fig:nonNewtonianConvergence">
 [[File:nonNewtonianConvergence.png|thumb|left|700px|
-<caption>Convergence</caption> of average Nusselt number value with increasing node count and comparison with reference values from reference<ref name="Turan2011">O. Turan, A. Sachdeva, N. Chakraborty in R. J. Poole, Laminar natural convection of power-law Fluids in a square enclosure with differentially heated side walls subjected to constant temperatures, Journal of Non-Newtonian Fluid Mechanics 166, 1049 (2011).</ref>. Convergence is calculated at Ra$=10^4$ and Pr$=10$.]]
+<caption>Convergence</caption> of the average Nusselt number value with the increasing node count and comparison with reference values from the reference<ref name="Turan2011">O. Turan, A. Sachdeva, N. Chakraborty in R. J. Poole, Laminar natural convection of power-law Fluids in a square enclosure with differentially heated side walls subjected to constant temperatures, Journal of Non-Newtonian Fluid Mechanics 166, 1049 (2011).</ref>. Convergence is calculated at Ra$=10^4$ and Pr$=10$.]]
 </figure>
@@ Line 43: / Line 43: @@
 \mathrm{Nu}=\frac{\Omega_W}{T_{H}-T_{C}}\left|\frac{\partial T}{\partial p_x}\right|_{p_x=0},
 $$
-where $T$ denotest the demperature, $p_x$ $x$ coordinate and $\Omega_W$ cavity width. Nusselt number is a convenient tool that can be used to reduce the dynamics of a complex multidimensional system into a single scalar value.
+where $T$ denotes the temperature, $p_x$ the $x$ coordinate and $\Omega_W$ the cavity width. Nusselt number is a convenient tool that can be used to reduce the dynamics of a complex multidimensional system into a single scalar value.
-We define the problem in terms of Rayleigh and Prandtl dimensionless numbers that need to be slightly modified to remain relevant in non-Newtonian case but specifics are not really relevant for this review. Dimensionless number definitions used in these calculations match reference<ref name="Turan2011">O. Turan, A. Sachdeva, N. Chakraborty in R. J. Poole, Laminar natural convection of power-law Fluids in a square enclosure with differentially heated side walls subjected to constant temperatures, Journal of Non-Newtonian Fluid Mechanics 166, 1049 (2011).</ref> and can be viewed there.
+We define the problem in terms of Rayleigh and Prandtl dimensionless numbers that need to be slightly modified to remain relevant in a non-Newtonian case but their specifics are not really relevant for this review. Dimensionless number definitions used in these calculations match reference<ref name="Turan2011">O. Turan, A. Sachdeva, N. Chakraborty in R. J. Poole, Laminar natural convection of power-law Fluids in a square enclosure with differentially heated side walls subjected to constant temperatures, Journal of Non-Newtonian Fluid Mechanics 166, 1049 (2011).</ref> and can be viewed there.
-Parameters that lead to oscillatory behaviour in de Vahl Davis cavity can be found in literature<ref name="Kosec2011">G. Kosec in B. Šarler, Solution of a low Prandtl number natural convection benchmark by a local meshless method, International Journal of Numerical Methods for Heat & Fluid Flow 23, 22 (2013).</ref>. Referenced article found oscillatory behaviour in Newtonian fluids at Ra=$5 \cdot 10^4$ and Pr=$0.01$ but we want to start slightly lower to find the transition as the dynamics become wilder with smaller $n$. We chose Ra=$2 \cdot 10^4$ and Pr=$0.01$ as the parameters to be used for all subsequently displayed calculations.
+Parameters that lead to an oscillatory behaviour in the de Vahl Davis cavity can be found in literature<ref name="Kosec2011">G. Kosec in B. Šarler, Solution of a low Prandtl number natural convection benchmark by a local meshless method, International Journal of Numerical Methods for Heat & Fluid Flow 23, 22 (2013).</ref>. The referenced article found oscillatory behaviour in Newtonian fluids at Ra=$5 \cdot 10^4$ and Pr=$0.01$, but we want to use a slightly lower Ra value to find the transition as the dynamics become wilder with a smaller $n$. We chose Ra=$2 \cdot 10^4$ and Pr=$0.01$ as the parameters to be used for all subsequently displayed calculations.
-The following figures display or mention dimensionless time that is defined with the same definition as described in reference<ref name="Kosec2011">G. Kosec in B. Šarler, Solution of a low Prandtl number natural convection benchmark by a local meshless method, International Journal of Numerical Methods for Heat & Fluid Flow 23, 22 (2013).</ref>.
+The following figures display or mention dimensionless time that is defined with the same definition as described in the reference<ref name="Kosec2011">G. Kosec in B. Šarler, Solution of a low Prandtl number natural convection benchmark by a local meshless method, International Journal of Numerical Methods for Heat & Fluid Flow 23, 22 (2013).</ref>.
 <figure id="fig:nusseltEvolutionNN">
 [[File:nusseltEvolutionNN.png|thumb|left|700px|
-<caption>Evolution of average Nusselt number in a differentially heated square cavity for different values of non-Newtonian power index $n$.</caption>]]
+<caption>Evolution of the average Nusselt number in a differentially heated square cavity for different values of the non-Newtonian power index $n$.</caption>]]
 </figure>
-We reuse the Nusselt number to analyse the differences in behaviour for different non-Newtonian power indices $n$ shown in <xr id="fig:nusseltEvolutionNN"/>. The Newtonian case with $n=1$ settles into a stationary case which is expected, as we chose parameters below the previously identified oscillatory solution. Evolutions of Nusselt number become oscillatory as we increase the strenght of non-Newtonian behaviour while the solutions become less and less regular as we reduce the power index $n$ which can be seen in pictures of velocity and temperature distributions shown in <xr id="fig:squareCavityDisplay"/>. Detailed animation of behaviour in the most extreme $n=0.7$ case is shown in <xr id="fig:animation"/>.
+We reuse the Nusselt number to analyse the differences in behaviour for different non-Newtonian power indices $n$ shown in <xr id="fig:nusseltEvolutionNN"/>. The Newtonian case with $n=1$ settles into a stationary case which is expected, as we chose parameters below the previously identified oscillatory solution. The time evolutions of the average Nusselt number become oscillatory as we increase the strenght of non-Newtonian behaviour, while the solutions become less and less regular as we reduce the power index $n$, which can be seen in pictures of velocity and temperature distributions shown in <xr id="fig:squareCavityDisplay"/>. Detailed animation of behaviour in the most extreme $n=0.7$ case is shown in <xr id="fig:animation"/>.
-Exact implementation of non-Newtonian fluid modeled with Ostwald-de Waele power law can be found in [[Non-Newtonian_fluid_example|example]].
+Exact implementation of the non-Newtonian fluid modeled with Ostwald-de Waele power law can be found in the [[Non-Newtonian_fluid_example|example]].
 <figure id="fig:squareCavityDisplay">
 [[File:squareCavity_t10_Ra20000Pr001.png|frame|center|
-<caption>Velocity and temperature fields within the cavity at $t=10$ for different values of non-Newtonian power index $n$.</caption>]]
+<caption>Velocity and temperature fields within the cavity at $t=10$ for different values of the non-Newtonian power index $n$.</caption>]]
 </figure>

@@ Line 67: / Line 67: @@
 <figure id="fig:animation">
-[[File:squareCavityAnimationNN_n07.mp4|frame|center|
+[[File:squareCavityAnimationNN_n07.mp4|thumb|center|upright=2|
 Animation of velocity and temperature fields within the cavity between $t=10$ and $t=11$ for $n=0.7$.]]
 </figure>