|
[Sponsors] |
contactAngle BC in contactAngleCavity tutorial is not working as expected |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 21, 2023, 07:36 |
contactAngle BC in contactAngleCavity tutorial is not working as expected
|
#1 |
New Member
Iason Tsiapkinis
Join Date: Jul 2019
Posts: 8
Rep Power: 7 |
Hello everybory,
this is a crosspost from https://develop.openfoam.com/Develop.../-/issues/2690 Maybe someone can help me here or point me in the right direction? I am trying to figure out, why the contactAngleCavity tutorial is not working as intended. Summary In the tutorial case contactAngleCavity, the free surface of a cavity is simulated with the `interfaceTrackingFvMesh` library. In the simulation a contact angle of 70° is set at the walls via the file `0.orig/contactAngle`. The contact angle after running the tutorial, however, is ~80°. Steps to reproduce
Code:
pointHistoryLeft1 { type pointHistory; refHistoryPoint (0 0.01 0); fileName leftPoint.txt; } pointHistoryLeft2 { type pointHistory; refHistoryPoint (0.00025 0.01 0); fileName leftPointRef.txt; } pointHistoryRight1 { type pointHistory; refHistoryPoint (0.01 0.01 0); fileName rightPoint.txt; } pointHistoryRight2 { type pointHistory; refHistoryPoint (0.00975 0.01 0); fileName rightPointRef.txt; } Example case I attached the tutorial case with slight modifications in the functionObjects, together with the post-processing script. Run the case using the Allrun script. [contactAngleCavity.tar] What is the current *bug* behaviour? The contact angle at the specified locations is not equal to the set contact angle. What is the expected *correct* behavior? The contact angle should be 70°, as specified in `0.orig/contactAngle`. Relevant logs and/or images The Figure below shows the calculated angle between the first two and the last two mesh points. Note that the angle does not change with a longer run time. The tutorial runs only for 0.2s. (contactAngle.png) Environment information OpenFOAM version : v2212 Operating system : ubuntu docker container on windows Possible fixes The calculations regarding the contact angle can be found in the file Code:
$FOAM_SRC/dynamicFaMesh/interfaceTrackingFvMesh/freeSurfacePointDisplacement.C Code:
patchMirrorPoints[patchI] = peCentres + ((I - 2*N*N)&delta); |
|
March 21, 2023, 06:23 |
|
#2 |
New Member
Iason Tsiapkinis
Join Date: Jul 2019
Posts: 8
Rep Power: 7 |
Any ideas or tips?
|
|
Tags |
contactangle, dynamicfvmesh, dynamicmesh, multiphase, tutorial |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Virtualization] OpenFOAM oriented tutorial on using VMware Player - support thread | wyldckat | OpenFOAM Installation | 2 | July 11, 2012 16:01 |
How to install CGNS under windows xp? | lzgwhy | Main CFD Forum | 1 | January 11, 2011 18:44 |
ParaView for OF-1.6-ext | Chrisi1984 | OpenFOAM Installation | 0 | December 31, 2010 06:42 |
CGNS lib and Fortran compiler | manaliac | Main CFD Forum | 2 | November 29, 2010 06:25 |
FoamToEnsight/EnsightParts/VTK not working with motorbike tutorial | madad2005 | OpenFOAM Bugs | 5 | October 18, 2009 04:48 |