CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[mesh manipulation] 2D Axisymmetric Mesh Not Working for Custom STL Multiphase Geometry

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 12, 2024, 11:16
Default 2D Axisymmetric Mesh Not Working for Custom STL Multiphase Geometry
  #1
New Member
 
Muhd Jamil Al Hakim
Join Date: Feb 2015
Posts: 5
Rep Power: 11
Jamill1283 is on a distinguished road
Send a message via Yahoo to Jamill1283 Send a message via Skype™ to Jamill1283
Dear Foamers,

I'm facing an issue with a multiphase flow simulation that involves modeling a liquid ligament in air. The geometry of the ligament, which is custom-designed using CAD software, has been exported to an STL file for this purpose.

My objective is to simulate this specific geometry in a 2D axisymmetric manner. However, I'm encountering a problem when the simulation is set up with a 2D mesh (having a mesh count of 1 in the z-direction); the liquid phase, which has a volume fraction of 1, does not appear correctly in the simulation. Interestingly, when I increase the mesh count in the z-direction to a higher number, the liquid's shape is accurately represented.

My blockmesh is as below;

Code:
scale 0.001;

blockWidth #eval "3";
ligamentSize #eval "0.5";
theta #eval "45.0 / 8";

arcLocZ #eval "($blockWidth) * sin(degToRad($theta))";
arcLocY #eval "($blockWidth) * cos(degToRad($theta))";
nArcLocZ #eval "-1*$arcLocZ";
ligArcLocZ #eval "($ligamentSize) * sin(degToRad($theta))";
ligArcLocY #eval "($ligamentSize) * cos(degToRad($theta))";

smallArcLocZ #eval "($blockWidth) * sin(degToRad($theta/2))";
smallArcLocY #eval "($blockWidth) * cos(degToRad($theta/2))";
smallLigArcLocZ #eval "($ligamentSize) * sin(degToRad($theta/2))";
smallLigArcLocY #eval "($ligamentSize) * cos(degToRad($theta/2))";

vertices
(
    (0  0 0) 
    (10 0 0) 
    (10 $blockWidth 0) 
    (0  $blockWidth 0) 
    (0 $arcLocY $arcLocZ)
    (10 $arcLocY $arcLocZ)
);

blocks
(
    hex (0 3 4 0 1 2 5 1) (100 100 1) simpleGrading (1 1 1) //the issue is here
);

edges
(
);

patches
(
    patch outlet
    (
        (2 5 1 1)
    )    
);

boundary
(
    front
    {
        type wedge;
        faces
        (
            (0 1 5 4)
            (0 1 2 3)
        );
    }
    ligament
    {
        type symmetry;
        faces
        (
            (3 4 0 0)
        );
    }
    external
    {
        type wall;
        faces
        (
            (3 4 5 2)
        );
    }
    axis
    {
        type  empty;
        faces
        (
            (0 1 1 0)
        );
    }
);

mergePatchPairs
(
);
Setfields as below;

Code:
defaultFieldValues
(
    volScalarFieldValue alpha.phase1 0
    volVectorFieldValue U (0 0 0)
);


regions (
   surfaceToCell {
   file           "./constant/fluidStruct/ligament108.stl";
   outsidePoints  ((0.008 0.00025 0));
   includeInside  true;
   includeOutside false;
   includeCut     true;
   nearDistance   -1;
   curvature      -100;
   fieldValues (
   volScalarFieldValue alpha.phase1 1
  );
 }
);
I'm seeking advice or solutions to ensure that the liquid phase is properly visualized while maintaining the simulation in a strictly 2D axisymmetric setup. Any insights or suggestions would be greatly appreciated.

Attached are visuals for 3D, problem for 2D mesh and zipped STL files;
Screenshot from 2024-03-12 14-29-41.png

Screenshot from 2024-03-12 14-20-39.png

ligament.zip

Thank you in advance for your assistance.
Jamil

Last edited by Jamill1283; March 13, 2024 at 08:32.
Jamill1283 is offline   Reply With Quote

Reply

Tags
axisymmetric wedge, mesh, multiphase

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 37 July 20, 2020 06:38
[ANSYS Meshing] Help with Mesh/Geomety of 2-D Axisymmetric Geometry yetlyg ANSYS Meshing & Geometry 1 July 1, 2015 03:01
[ICEM] Working on imported STL mesh file saitej09 ANSYS Meshing & Geometry 0 January 16, 2015 07:11
[ANSYS Meshing] Meshing Strategy for inside geometry powpow ANSYS Meshing & Geometry 6 January 16, 2013 05:32
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 19:10


All times are GMT -4. The time now is 12:02.