View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0004207 | OpenFOAM | Bug | public | 2025-01-30 15:20 | 2025-01-30 15:37 |
| Reporter | Tobermory | Assigned To | henry | ||
| Priority | low | Severity | minor | Reproducibility | always |
| Status | closed | Resolution | no change required | ||
| Platform | GNU/Linux | OS | Other | OS Version | (please specify) |
| Product Version | 12 | ||||
| Summary | 0004207: fixedShearStress boundary sets the kinematic shear stress | ||||
| Description | I am not sure whether this is truly a "bug", or just a different convention. It is confusing though and I don't see why it was done the way it was (unless the class was developed originally for an incompressible solver). The fixedShearStressFvPatchVectorField class is described in its header as: "Set a constant shear stress as tau0 = -nuEff dU/dn" and indeed the class constructor has a default value with dimensions of velocity squared: tau0_(dict.lookupOrDefault<vector>("tau", sqr(dimVelocity), Zero)) My problem is that this is a *kinematic* shear stress, and it's more normal to use dynamic stress, i.e. rho times the above, especially in a compressible solver. You might counter "well just enter the kinematic value", which is fine if you know the density, i.e. pressure & T, a-priori, but sometimes you might not ... | ||||
| Steps To Reproduce | Read the code. | ||||
| Tags | No tags attached. | ||||
