|
[Sponsors] |
How to modify viscosity field within a specific region? |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14 ![]() |
Suppose a moving solid body in a two-fluid (e.g. Gas + Liquid), the viscosity of the mixture of two-fluid in the fluid region is:
Code:
mu_fluid = alpha_G* mu_G + alpha_L* mu_L; Code:
alpha_G + alpha_L = 1.0; Code:
mu = alpha_G* mu_G + alpha_L* mu_L + alpha_S* mu_S; rho = alpha_G* rho_G + alpha_L* rho_L + alpha_S* rho_S; Code:
alpha_G + alpha_L + alpha_S =1 Since volume fraction field (alpha_S) can be explicitly calculated, before solving the momentum equation and alpha equation (for gas and liquid), we can model these three-phase flow with a modified InterFoam. The problem remaining is that: how do we effectively modify the fluid properties inside the solid region, if the cells (ID) inside the solid body are known? The first way comes to my mind is rather straight-forward: Code:
for (auto& cID : cellsInSolid) { // fluid-solid boundary cells + inner solid cells mu[cID] = alpha_S* mu_S + alpha_G* mu_G + alpha_L* mu_L; rho[cID] = alpha_S* rho_S + alpha_G* rho_G + alpha_L* rho_L; } Thanks in advance. ![]() |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Robert
Join Date: May 2015
Location: Bremen, GER
Posts: 292
Rep Power: 12 ![]() |
You can try using a fvOptions source. Take a look at these examples.
__________________
If you liked my answer to your question, please consider leaving a "Like" in return ![]() |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
Senior Member
David Long
Join Date: May 2012
Location: Germany
Posts: 104
Rep Power: 14 ![]() |
Quote:
Can you explain a bit more in detail? |
||
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
p_rgh initial residual no change with different settings | manuc | OpenFOAM Running, Solving & CFD | 3 | June 26, 2018 15:53 |
Allocate a specific region to a specific processor in a parallel run | Saideep | OpenFOAM Pre-Processing | 4 | October 5, 2017 06:04 |
Moving mesh | Niklas Wikstrom (Wikstrom) | OpenFOAM Running, Solving & CFD | 122 | June 15, 2014 06:20 |
User field function for specific heat deifnition | crevoise | STAR-CCM+ | 2 | January 28, 2013 04:25 |
How to modify the viscosity model | mpml | OpenFOAM Running, Solving & CFD | 4 | October 13, 2010 07:44 |