CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[mesh manipulation] manipulation of an unstructerd mesh

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 17, 2020, 14:19
Question manipulation of an unstructerd mesh
New Member
Join Date: Sep 2019
Posts: 2
Rep Power: 0
HandsomePotato is on a distinguished road
hey foamers,

I'm doing an exhaust gas dispersion simulation with OF4.1. I have recieved an unstructered mesh from an external party. This is a circular domain with a ship in the middle. I want to check the effect of a modified velocity profile on the exhaust gas dispersion. on the domain boundary i am using groovyBC with tabulated values:
/0 -> U
type groovyBC;
fields (z U);
name Uval;
outOfBounds clamp;
fileName "$FOAM_CASE/0/";

valueExpression "(vector(cos(a*pi/180),sin(a*pi/180),0)*Uval(pos().z)-ship)";
gradientExpression "vector(0,0,0)";
fractionExpression "(phi > 0 ) ? 0 : 1";

which is working !!
However the domain has a radius of 2500 meter, when the velocity has reached the ship the profile has already grown towards an atmospheric boundary layer.
To solve this i want to cut this mesh on the windward side such that it only has a 100 meter to reach the ship.
I tried to do this using topoSet and than subset mesh
/system -> topoSetDict:
name a;
type cellSet;
action new;
source boxToCell;
box (3000 3000 500)(-3000 100 0);

name a;
type cellSet;
action invert;

>>>topoSet &>log.topoSet
Time = 0
mesh not changed.
Created cellSet a
Applying source boxToCell
Adding cells with center within boxes 1((3000 3000 500) (-3000 100 0))
cellSet a now size 0
Time = 0
mesh not changed.
Read set cellSet a with size 0
Inverting cellSet
cellSet a now size 37343105

>>>subsetMesh -overwrite a -patch Domain &>log.subsetMesh
Subsetting field nut
Subsetting field k
Subsetting field S
Subsetting field p_rgh
Subsetting field alphat
Subsetting field p
Subsetting field T
Subsetting field omega
Subsetting field U
--> FOAM Warning :
From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(c onst fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
in file groovyBCFvPatchField.C at line 131
No value defined for U on Domain therefore using the internal field next to the patch
Writing subsetted mesh and fields to time 0

I think the FOAM warning comes from the fact that the table is only used when actually simulating, thus it does not have a value when manipulating the mesh. When launching the simulation this warning is found in every log-file until the solver is employed and just runs smoothly.

i tried out topoSet with subsetMesh command on a simple blockmesh which works as intended, but when i try it on this unstructured mesh it does not change the mesh when i review the result in paraview.

Does anybody know what i am doing wrong, or does someone have a different method for manipulation a unstructured mesh? Furthermore another method on how to impose a modified velocity profile in the internalfield of this mesh is more than welcome.

I know openFOAM able to do this!!!
Thanks in advance.
HandsomePotato is offline   Reply With Quote


mesh manipulation, subsetmesh, toposet, unstructured mesh

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
[snappyHexMesh] High quality mesh for wind in complex urban environment ziboaa OpenFOAM Meshing & Mesh Conversion 1 January 12, 2021 15:33
[snappyHexMesh] Creating multiple multiple cell zones with snappyHexMesh - a newbie in deep water! divergence OpenFOAM Meshing & Mesh Conversion 0 January 23, 2019 04:17
[ICEM] 2D hybrid mesh (unstructured mesh highly dependent on structured mesh parameters) shubham jain ANSYS Meshing & Geometry 1 April 10, 2017 05:03
[ICEM] Generating Mesh for STL Car in Windtunnel Simulation tommymoose ANSYS Meshing & Geometry 48 April 15, 2013 04:24
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 11:55

All times are GMT -4. The time now is 19:07.