CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[swak4Foam] Set surface normal manually in groovyBC

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By Mitchell Baum

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 25, 2016, 20:22
Default Set surface normal manually in groovyBC
  #1
New Member
 
Mitchell Baum
Join Date: Sep 2015
Location: Australia
Posts: 7
Rep Power: 10
Mitchell Baum is on a distinguished road
Hi all,

I am using groovyBC to define a power law velocity profile at a pipe outlet. According to the attached figure I have a nozzle in a rectangular domain which I have meshed using the snappyHexMesh utility. Due to the limitations of snappyHexMesh I have some surface irregularities at the edge of the orifice. When I implement the line:
Code:
valueExpression "-Un*pow((1-(r/R)),(1/n))*normal()";
I understand that when I define the according velocity normal in the valueExpression I am defining velocity vectors to the surface of the constituent cells of that patch, however due to my surface irregularities at the orifice edge I am getting some vectors pointing away from my desired direction. Is it possible to set the direction using the surface vector say (5.000000e-001 3.041372e-014 8.660254e-001) (which corresponds directly from my patch .stl file) as opposed to using normal()?

Thanks in advance for your help.
Mitchell
Attached Images
File Type: jpg PowerLawNozzleVelocityProfile.JPG (77.9 KB, 18 views)

Last edited by Mitchell Baum; March 1, 2016 at 00:18.
Mitchell Baum is offline   Reply With Quote

Old   March 15, 2016, 00:25
Default
  #2
New Member
 
Mitchell Baum
Join Date: Sep 2015
Location: Australia
Posts: 7
Rep Power: 10
Mitchell Baum is on a distinguished road
While I didn't manage to attain a solution to my original question in this post, I achieved a much more sound solution by adjusting my snappyHexMesh processes. Prior to running the snappyHexMesh command I ran the surfaceFeatureExtract command while I also changed around some settings in snappyHexMeshDict. I have attached my most recent mesh.

Although I have solved my original issue in this case, I feel it would be beneficial if anyone could offer guidance in assigning a predefined vector to a patch.

Regards,
Mitchell
Attached Images
File Type: jpg UpdatedPowerLawVelocityProfile.JPG (89.4 KB, 12 views)
Thamali likes this.
Mitchell Baum is offline   Reply With Quote

Reply

Tags
groovybc, surface vector


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
Normal vector of the free surface kknv001 FLOW-3D 1 February 11, 2021 06:56
Ansys Licence Serve on Ubuntu 16.04 LTS david.pasquale ANSYS 2 January 20, 2017 11:52
Set surface normal component of a vector to zero ChrisA OpenFOAM Running, Solving & CFD 0 July 19, 2013 14:27
How to show the transient case? H.P.LIU Phoenics 7 July 13, 2010 04:31
[GAMBIT] Construction of surface by set of points pankos ANSYS Meshing & Geometry 0 June 19, 2010 08:25


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