View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0001261 | OpenFOAM | Bug | public | 2014-04-05 11:21 | 2014-04-08 18:25 |
| Reporter | Assigned To | will | |||
| Priority | normal | Severity | major | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Summary | 0001261: Wrong solid viscosity in JohnsonJacksonParticleSlip boundary condition | ||||
| Description | The JohnsonJacksonParticleSlip boundary condition uses the solid bulk viscosity defined in the file "thermophysicalProperties.particles" instead of using the effective solid viscosity calculated by the kinetic theory model. | ||||
| Tags | No tags attached. | ||||
|
|
Yes, I can confirm this also. The problem is in lines 200-203 in JohnsonJacksonParticleSlipFvPatchVectorField.C. const scalarField nu ( phased.nu()->boundaryField()[patch().index()] ); This uses thermo-nu and not kinetic theory nu. |
|
|
Thanks for the report. I've fixed this in commit eb086ae45503be15714cf2b491df0edfcaa91c92 by looking up the turbulent viscosity calculated by the kinetic theory turbulence model. This includes both the shear and frictional components. It isn't terribly clear in my reference*, whether it should be the combined viscosity, or just the shear value used in the boundary condition. If you think it should be the latter, please let us know. * http://www.sciencedirect.com/science/article/pii/S0009250908003904 |
