
[Sponsors] 
April 23, 2018, 08:40 

#21 
Member
Peng Liang
Join Date: Mar 2014
Posts: 59
Rep Power: 12 
Quote:
Bests, Peng 

April 23, 2018, 08:49 

#22  
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18 
Quote:
is the normal vector to the gas/liquid interface. In your problem, I will say it depends. If you have an immiscible interface between your two fluids. Yes, you can go for a Volume of Fluid solver as a starting point. 

July 2, 2018, 14:29 

#23 
New Member
Rajesh Singh
Join Date: Jun 2010
Posts: 15
Rep Power: 16 
Dear FOAMer,
I have imported StarCCM+ mesh to openfoam 5.0 and quality of mesh was also good. The CCM mesh is trimmed mesh with polyhedral cells. The flow simulation for interFoam with species transport equation (Haroun formulation)diverges after some iterations. Below is report of the checkMesh. Help for resolving this issue would be appreciated. Thanks in advance ******************************* Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 2709670 faces: 7273254 internal faces: 6550582 cells: 2217107 faces per cell: 6.23508 boundary patches: 10 point zones: 0 face zones: 1 cell zones: 1 Overall number of cells of each type: hexahedra: 1558763 prisms: 17795 wedges: 1051 pyramids: 73 tet wedges: 51 tetrahedra: 184093 polyhedra: 455281 Breakdown of polyhedra by number of faces: faces number of cells 4 4 5 43 6 76272 7 247451 8 4264 9 70800 10 3318 11 356 12 16481 13 700 14 204 15 35328 16 60 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology floweriod1 110302 113851 ok (nonclosed singly connected) flow:sheet1 229147 166554 ok (nonclosed singly connected) flow:sheet2 228457 166118 ok (nonclosed singly connected) flowutg 9532 9665 ok (nonclosed singly connected) flow:wallt 9576 9998 ok (nonclosed singly connected) floweriod2 110292 113851 ok (nonclosed singly connected) flow:inl 5071 5594 ok (nonclosed singly connected) flow:wallb 10551 11107 ok (nonclosed singly connected) flowutl 5792 6040 ok (nonclosed singly connected) flow:ing 3952 4069 ok (nonclosed singly connected) Checking geometry... Overall domain bounding box (0.0133625 0.013 0.0851) (0.0133628 0.013 0.007) Mesh (nonempty, nonwedge) directions (1 1 1) Mesh (nonempty) directions (1 1 1) Boundary openness (4.18881e14 2.88815e15 1.82866e17) OK. Max cell openness = 2.78469e16 OK. Max aspect ratio = 9.20971 OK. Minimum face area = 1.10391e11. Maximum face area = 6.4e07. Face area magnitudes OK. Min volume = 5.01959e16. Max volume = 5.12e10. Total volume = 3.4427e05. Cell volumes OK. Mesh nonorthogonality Max: 57.8894 average: 11.7876 Nonorthogonality check OK. Face pyramids OK. Max skewness = 3.2775 OK. Coupled point location match (average 0) OK. Mesh OK. End 

June 7, 2022, 07:55 

#24  
Member
sadra mahmoudi
Join Date: Feb 2021
Location: Austria
Posts: 39
Rep Power: 5 
Quote:
Hello Cyp, I am trying to simulate a single bubble movement in a solution of water and sugar with interIsoFoam solver, OF2112. I modified the solver and coupled the density, surface tension and viscosity of solution to the concentration of sugar which is different in various parts of the domain. In order to solve the distribution of sugar (a passive scalar) in the geometry, I added a new equation to the solver as below: fvScalarMatrix CEqn ( fvm::ddt(C) + fvm::div(phi, C)  fvm::laplacian(dc,C) == fvOptions(C) ); CEqn.relax(); fvOptions.constrain(CEqn); CEqn.solve(); fvOptions.correct(C); Now, the problem is, sugar concentration penetrates inside the bubble which is not correct. I would like to know how can I prevent sugar entering the bubble? I would be more than happy if you share your opinion with me. 

June 7, 2022, 07:57 
a question

#25  
Member
sadra mahmoudi
Join Date: Feb 2021
Location: Austria
Posts: 39
Rep Power: 5 
Quote:
Hello Astrodan, I am trying to simulate a single bubble movement in a solution of water and sugar with interIsoFoam solver, OF2112. I modified the solver and coupled the density, surface tension and viscosity of solution to the concentration of sugar which is different in various parts of the domain. In order to solve the distribution of sugar (a passive scalar) in the geometry, I added a new equation to the solver as below: fvScalarMatrix CEqn ( fvm::ddt(C) + fvm::div(phi, C)  fvm::laplacian(dc,C) == fvOptions(C) ); CEqn.relax(); fvOptions.constrain(CEqn); CEqn.solve(); fvOptions.correct(C); Now, the problem is, sugar concentration penetrates inside the bubble which is not correct. I would like to know how can I prevent sugar entering the bubble? I would be more than happy if you share your opinion with me. 

Tags 
drop, interface position, interfoam, openfoam 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Wind turbine simulation  Saturn  CFX  58  July 3, 2020 01:13 
An error has occurred in cfx5solve:  volo87  CFX  5  June 14, 2013 17:44 
RPM in Wind Turbine  Pankaj  CFX  9  November 23, 2009 04:05 
Convective Heat Transfer  Heat Exchanger  Mark  CFX  6  November 15, 2004 15:55 
Replace periodic by inletoutlet pair  lego  CFX  3  November 5, 2002 20:09 