|
[Sponsors] |
![]() |
![]() |
#1 | |
Senior Member
Join Date: Dec 2010
Posts: 135
Rep Power: 15 ![]() |
Hi guys,
i want to create a linear velocity profile at the inlet patch for simulating couette flow in a channel with less computational effort, therefore I tried following code: Quote:
Last edited by eRzBeNgEl; January 9, 2014 at 07:22. |
||
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Join Date: Dec 2010
Posts: 135
Rep Power: 15 ![]() |
any idea? my results are not right, so there must be a mistake
|
|
![]() |
![]() |
![]() |
![]() |
#3 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 ![]() ![]() |
Quote:
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Join Date: Dec 2010
Posts: 135
Rep Power: 15 ![]() |
Hi,
in general I will produce a linear gradient velocity profile at my inlet in y-direction This is my test setup 1. I tried by groovyBC but I failed - 2. I changed the hard coded parabolicVelocity.C File of HJasak and got the profile at figure 1. - it seems to be right BUT if i calculate further I got wiggles at the corner on bottom. How can I avoid them?(see Figure 2 & Figure 3) and in Couette there should be no velocity gradient in y-direction, but there is! and in general, how can I set this profile to whole domain and not just to the inlet? Here is my Member Function of the linear velocity profile: Code:
void linearVelocityFvPatchVectorField::updateCoeffs() { if (updated()) { return; } // Get range and orientation boundBox bb(patch().patch().localPoints(), true); vector ctr = (bb.max() + bb.min()); const vectorField& c = patch().Cf(); // Calculate local 1-D coordinate for the linear profile scalarField coord = ((c - ctr) & y_)/((bb.max() - bb.min()) & y_); //n=flow direction, maxValue=peakVelocity vectorField::operator=(n_*maxValue_*(1.0-sqr(coord))); |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Greetings to all!
@eRzBeNgEl: There are several details here that are either missing or that you are not taking into account:
Bruno
__________________
|
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Velocity inlet BC profile | antonio.c | FLUENT | 1 | November 1, 2017 14:22 |
atmBoundaryLayerInletVelocity - Velocity Profile not continuous through domain | sdfij6354 | OpenFOAM Running, Solving & CFD | 3 | July 26, 2017 16:16 |
extracting outlet velocity profile from one case to another case's inlet | tonggysun | OpenFOAM | 2 | September 13, 2013 04:19 |
Export & Import Velocity profile as Inlet | eRzBeNgEl | STAR-CCM+ | 6 | March 26, 2012 05:16 |
[swak4Foam] Scale discrete inlet velocity profile with groovyBC | cboss | OpenFOAM Community Contributions | 1 | June 20, 2010 13:02 |