CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > SU2

Non-uniform inlet conditions?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 16, 2013, 11:38
Post Non-uniform inlet conditions?
  #1
New Member
 
Magnus W
Join Date: Jul 2013
Location: Munich - Germany
Posts: 5
Rep Power: 3
winter is on a distinguished road
Hi everyone!

I was interested in validating SU2's incompressible solver to the Pitz Daily (backward facing step) experiment. To conduct the simulations it would be beneficial to use a profile for the velocity and the turbulence quantities (k and omega), as the entrance length is quite long.

I haven't, so far, found that SU2 is capable to use profiled inlet conditions. Has anybody found a way to use non-uniform inlet conditions in SU2?
winter is offline   Reply With Quote

Old   October 3, 2013, 14:18
Default
  #2
Super Moderator
 
Francisco Palacios
Join Date: Jan 2013
Location: Stanford, CA
Posts: 301
Rep Power: 5
fpalacios is on a distinguished road
Quote:
Originally Posted by winter View Post
Hi everyone!

I was interested in validating SU2's incompressible solver to the Pitz Daily (backward facing step) experiment. To conduct the simulations it would be beneficial to use a profile for the velocity and the turbulence quantities (k and omega), as the entrance length is quite long.

I haven't, so far, found that SU2 is capable to use profiled inlet conditions. Has anybody found a way to use non-uniform inlet conditions in SU2?
Dear Magnus,
Thanks a lot for your interest in SU2. Most of the work in SU2 has been done for external flows, and so far we don't have a method to introduce an inlet profile.

However, it is not complicate to add this capability, and we are looking forward for contributions. A good starting point would be to modify.

void CEulerSolver::BC_Inlet(CGeometry *geometry, CSolver **solver_container, CNumerics *conv_numerics, CNumerics *visc_numerics, CConfig *config, unsigned short val_marker)

in solver_direct_mean.cpp and instead of using the velocity at the Infinity (GetVelocity_Inf) create a function that reads the velocity profile from a file and interpolate the solution given a (X,YZ) coordinate geometry->node[iPoint]->GetCoord().

As you can see, SU2 use the weak formulation, so a numerical flux will be computed between the value that you provide and the actual value at the node, and then the result will be added to the total residual at the point.

Our recommendation is to fork the SU2 repository in https://github.com/su2code and contact us if you need further support.

Thanks again for your interest,
Francisco
fpalacios is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Channel flow using InterFOAM DanM OpenFOAM Running, Solving & CFD 38 June 19, 2015 09:45
fluentMeshToFoam multidomain mesh conversion problem Attesz OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 12 May 2, 2013 10:52
Initial conditions for uniform flow andreas OpenFOAM 5 November 16, 2012 16:00
Trouble setting groovyBC interFoam Inlet BC mgdenno OpenFOAM 4 July 1, 2012 08:16
chtMultiRegionSimpleFoam javad814 OpenFOAM 1 September 26, 2011 13:30


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