CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Surface shear force in OpenFOAM-v1706

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   August 20, 2017, 23:41
Post Surface shear force in OpenFOAM-v1706
New Member
Vishwesh Ravi Shrimali
Join Date: Aug 2017
Posts: 1
Rep Power: 0
vishwesh is on a distinguished road

I am working on a rivulet panel case and am trying to see the effect of surface shear force on a falling waterfilm due to crossflow. I have prepared the test case and had updated the reactingParcelFilmFoam solver according to my test-case.
Now, in OpenFOAM 2.1.1, there is a separate kinematic force - surfaceShearForce which can be enabled or disabled using the surfaceFilmProperties file in constant folder. In OpenFOAM versions 3.0.x and after (I tried both OpenFOAM 3.0.x and v1706), the surfaceShearForce is absent and instead of that, a member function has been included in the laminar turbulence model (src/regionModels/surfaceFilmModels/submodels/kinematic/filmTurbulenceModel/laminar/laminar.C)

 tmp<fvVectorMatrix> laminar::Su(volVectorField& U) const
  125 {
  126     // local reference to film model
  127     const kinematicSingleLayer& film =
  128         static_cast<const kinematicSingleLayer&>(filmModel_);
  130     // local references to film fields
  131     const volScalarField& mu =;
  132     const volVectorField& Uw = film.Uw();
  133     const volScalarField& delta =;
  134     const volVectorField& Up = film.UPrimary();
  135     const volScalarField& rhop = film.rhoPrimary();
  137     // employ simple coeff-based model
  138     volScalarField Cs("Cs", Cf_*rhop*mag(Up - U));
  139     volScalarField Cw("Cw", mu/(0.3333*(delta + film.deltaSmall())));
  140     Cw.min(5000.0);
  142     return
  143     (
  144        - fvm::Sp(Cs, U) + Cs*Up // surface contribution
  145        - fvm::Sp(Cw, U) + Cw*Uw // wall contribution
  146     );
  147 }
(Taken from: OpenFOAM v1706 source code)

If this is compared with the earlier surfaceShearForce code in OpenFOAM 2.1.1, it can be seen that the equations and the procedure is more or less same, except a few changes here and there (like in the case of Cw.min).

When I am running the case in OpenFOAM 2.1.1, I am getting a deviation in the film but when I am running the case in OpenFOAM 3.0.x and OpenFOAM -v1706, there is no such deviation.
I have also modified the laminar code to find out if the function is being called or not and it is being called.

Can anyone help me out in figuring why this difference is there?
vishwesh is offline   Reply With Quote


openfoam, rivulentpanel, surfaceshearforce

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Map of the OpenFOAM Forum - Understanding where to post your questions! wyldckat OpenFOAM 9 March 30, 2017 05:19
impact force (shear stress) by waterjet impingement herozhu Main CFD Forum 0 May 7, 2014 20:39
Access the surface tension force at a patch michielm OpenFOAM Programming & Development 1 December 13, 2013 04:50
[ICEM] Problems with coedge curves and surfaces tommymoose ANSYS Meshing & Geometry 0 August 5, 2011 16:02
gravitational force for free surface flow Jongtae Kim Main CFD Forum 1 July 2, 2000 11:57

All times are GMT -4. The time now is 23:21.