A generalized thermal/dynamic wall function: Part 1
In a previous post i wrote about an extension of the Reichardt law of the wall to pressure gradient effects. That was derived by assuming the Reichardt profile as a solution for the case without pressure gradient and using integration by parts. In particular, given the the Reichardt function ( is the Von Karman constant):
with:
the pressure gradient sensitized formula for the velocity is given by:
where is the nondimensional pressure gradient.
However, while useful for its purpose, this formulation turns out to have at least two drawbacks, both of which depending, somehow, on the fact that Reichardt didn't actually solved the underlying problem exactly. The first one is that the base profile (i.e., even without the pressure gradient effect) does not even reproduce the correct profile derivatives at the wall. Hence, as a main consequence, the formula might not be suitable as toy model for certain theoretical purposes. The second drawback, which is also another consequence of the same fact, is that it can be used as thermal wall function only for , as higher values will tend to amplify the part of the profile which is uncorrect.
To overcome such issues, i follow here the same steps as before but with a different formulation for , using one which can be integrated exactly. We start from the temperature equation, which is more general:
where (constant pressure specific heat), (dynamic viscosity), (Prandtl number), (turbulent Prandtl number) and (defined in the following) are all constant. The fact that this equation correctly represents a large set of velocity/temperature viscous boundary conditions should be self-evident. For example, gives the classical linear relationship for the temperature (which becomes the one for the velocity when ). Also, with leads to the classical wall-function formulation (for certain forms of ). According to the specific equation form, can have different meanings. For example, in the case of the velocity, it can include the wall-parallel pressure gradient (note that y is assumed to be the wall-normal coordinate, x the wall parallel one), the buoyancy term, etc. In general, it can include any of the terms that have been neglected in the LHS for the given equation. The only assumption here is that they are assumed constant along y (as will be shown in the following, it turns out that this can be actually relaxed).
The first step to proceed is integrating the equation above once. This produces (using ):
where (the underscore meaning that it is evaluated at the wall, i.e., for ). Multiplying by and rearranging, we obtain:
where:
The same equation, interpreted for the velocity, would read:
where . This confirms that the temperature solution can also be used for the velocity by a proper interpretation of the terms. Going back to the temperature equation, we then end up with the following equation for the nondimensional temperature derivative:
Then, by a further integration:
Note that i temporarily suspended the hypothesis that is constant. Now, let us also assume that we know the first integral, that is:
is a known function. Then, the solution becomes:
And after integration by parts:
For the sake of conciseness, let us also assume:
so that the full solution is:
While we have just moved the original problem to that of solving the two integrals for and , those will actually turn out to be easier to solve. Note also that the parenthesis in is a polynomial of the same order of (assuming that it is approximated with a certain degree polynomial). Particularly relevant is the case with constant (as is typically assumed). However, note that ,for the that we are going to derive below, a linear approximation is also integrable, as well as any other polynomial of higher degree. Obviously, things rapidly become cumbersome, with the degree of the polynomial, so for the present purposes it is sufficient to go back to the case constant. In this case is simply the integral of times . Note that Popovac and Hanjalic have shown (for the velocity) that this is actually relevant for several practical cases. Nonetheless, is worth remembering that this is not a costriction of the present approach (probably it wasn't neither when was based on the Reichardt profile).
So, at this point we are left with determining . To this end, we need a model for . Note that this is not different from the approach used by Van Driest to determine his velocity profile. It's just that the mixing length hypothesis he used led to a non integrable function. The same also happened to Reichardt, who then started modificating his function until something integrable came out.
To overcome this limit, here we abandon those formulations and use instead the one from Musker:
with being, in its original formulation, a constant (see below). It is worth mentioning that this functional form respects both the behavior for and the behavior for .
Note also that the Spalart-Allmaras model has a simlar behavior near the wall (see Allmaras et al., ICCFD 2012):
which then fails to recover the correct behavior for (However, Spalart and Allmaras explicitly mention in their original 1994 work that thay didn't notice any meaningful difference deriving from using this specific form instead of a correct one, say, the mixing length). This just to say that the Musker form is not only meaningful, but also representative of the near wall behavior of the Spalart-Allmaras model.
To solve for with the Musker turbulent viscosity, we first consider the case , as this already has a known solution, provided in its most general form by Monkewitz et al. (POF 20, 105102, 2008):
where:
The constant is found by calibration with the turbulent velocity profile. Monkewitz et al. use . In contrast, i calibrated the constant using the mixing length model with a constant , resulting in .
Now, the whole point of the procedure is: how to extend such formulation to the case ? See part 2
with:
the pressure gradient sensitized formula for the velocity is given by:
where is the nondimensional pressure gradient.
However, while useful for its purpose, this formulation turns out to have at least two drawbacks, both of which depending, somehow, on the fact that Reichardt didn't actually solved the underlying problem exactly. The first one is that the base profile (i.e., even without the pressure gradient effect) does not even reproduce the correct profile derivatives at the wall. Hence, as a main consequence, the formula might not be suitable as toy model for certain theoretical purposes. The second drawback, which is also another consequence of the same fact, is that it can be used as thermal wall function only for , as higher values will tend to amplify the part of the profile which is uncorrect.
To overcome such issues, i follow here the same steps as before but with a different formulation for , using one which can be integrated exactly. We start from the temperature equation, which is more general:
where (constant pressure specific heat), (dynamic viscosity), (Prandtl number), (turbulent Prandtl number) and (defined in the following) are all constant. The fact that this equation correctly represents a large set of velocity/temperature viscous boundary conditions should be self-evident. For example, gives the classical linear relationship for the temperature (which becomes the one for the velocity when ). Also, with leads to the classical wall-function formulation (for certain forms of ). According to the specific equation form, can have different meanings. For example, in the case of the velocity, it can include the wall-parallel pressure gradient (note that y is assumed to be the wall-normal coordinate, x the wall parallel one), the buoyancy term, etc. In general, it can include any of the terms that have been neglected in the LHS for the given equation. The only assumption here is that they are assumed constant along y (as will be shown in the following, it turns out that this can be actually relaxed).
The first step to proceed is integrating the equation above once. This produces (using ):
where (the underscore meaning that it is evaluated at the wall, i.e., for ). Multiplying by and rearranging, we obtain:
where:
The same equation, interpreted for the velocity, would read:
where . This confirms that the temperature solution can also be used for the velocity by a proper interpretation of the terms. Going back to the temperature equation, we then end up with the following equation for the nondimensional temperature derivative:
Then, by a further integration:
Note that i temporarily suspended the hypothesis that is constant. Now, let us also assume that we know the first integral, that is:
is a known function. Then, the solution becomes:
And after integration by parts:
For the sake of conciseness, let us also assume:
so that the full solution is:
While we have just moved the original problem to that of solving the two integrals for and , those will actually turn out to be easier to solve. Note also that the parenthesis in is a polynomial of the same order of (assuming that it is approximated with a certain degree polynomial). Particularly relevant is the case with constant (as is typically assumed). However, note that ,for the that we are going to derive below, a linear approximation is also integrable, as well as any other polynomial of higher degree. Obviously, things rapidly become cumbersome, with the degree of the polynomial, so for the present purposes it is sufficient to go back to the case constant. In this case is simply the integral of times . Note that Popovac and Hanjalic have shown (for the velocity) that this is actually relevant for several practical cases. Nonetheless, is worth remembering that this is not a costriction of the present approach (probably it wasn't neither when was based on the Reichardt profile).
So, at this point we are left with determining . To this end, we need a model for . Note that this is not different from the approach used by Van Driest to determine his velocity profile. It's just that the mixing length hypothesis he used led to a non integrable function. The same also happened to Reichardt, who then started modificating his function until something integrable came out.
To overcome this limit, here we abandon those formulations and use instead the one from Musker:
with being, in its original formulation, a constant (see below). It is worth mentioning that this functional form respects both the behavior for and the behavior for .
Note also that the Spalart-Allmaras model has a simlar behavior near the wall (see Allmaras et al., ICCFD 2012):
which then fails to recover the correct behavior for (However, Spalart and Allmaras explicitly mention in their original 1994 work that thay didn't notice any meaningful difference deriving from using this specific form instead of a correct one, say, the mixing length). This just to say that the Musker form is not only meaningful, but also representative of the near wall behavior of the Spalart-Allmaras model.
To solve for with the Musker turbulent viscosity, we first consider the case , as this already has a known solution, provided in its most general form by Monkewitz et al. (POF 20, 105102, 2008):
where:
The constant is found by calibration with the turbulent velocity profile. Monkewitz et al. use . In contrast, i calibrated the constant using the mixing length model with a constant , resulting in .
Now, the whole point of the procedure is: how to extend such formulation to the case ? See part 2
Total Comments 0