CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

obtaining significantly different Flow profiles after only changing Mesh size

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 30, 2020, 13:06
Default obtaining significantly different Flow profiles after only changing Mesh size
  #1
New Member
 
berat
Join Date: Jun 2020
Posts: 6
Rep Power: 6
beratcagan is on a distinguished road
Hello,
so i am new to OpenFOAM and simulation in general and i am currently trying to simulate a two dimentional pipe flow. I managed to obtain a satisfiying result/flow profile in my simulation.
Now I am looking to making the simualtion bigger by simply changing the value form the "converttoMeter" command in the BlockMesh file. When i do so i am getting a completly unrealistic flow profile, even though every other part of the Simualtion is exactly the same. I allready tried increasing/decreacing the velocity at the input patch, this didnt change anything though.

BlockMesh of the working Simulation:

convertToMeters 0.000001;

vertices
(
(0 0 0)
(35000 0 0)
(35000 3500 0)
(0 3500 0)
(0 0 65)
(35000 0 65)
(35000 3500 65)
(0 3500 65)
);


blockMesh of the faulty simulation

convertToMeters 1;

vertices
(
(0 0 0)
(35000 0 0)
(35000 3500 0)
(0 3500 0)
(0 0 65)
(35000 0 65)
(35000 3500 65)
(0 3500 65)
);

as i already mentioned the rest of the simulation is exactly the same, despite me tryed to change up the velocity which didnt change anything though.

now i am curious what i am doing wrong. i would appreate the sightest hints which would help me solve my problem.
Attached Images
File Type: jpg correct.jpg (66.0 KB, 7 views)
File Type: jpg faulty.jpg (61.8 KB, 9 views)
beratcagan is offline   Reply With Quote

Old   June 30, 2020, 16:07
Default
  #2
Senior Member
 
Troy Snyder
Join Date: Jul 2009
Location: Akron, OH
Posts: 220
Rep Power: 19
tas38 is on a distinguished road
In your faulty simulation, the length scales are increased dramatically. Correspondingly, your Reynolds number has increased significantly given the same U and fluid properties rho and mu.

Is the solution converging for the high-Re case?

Is the flow modeled as laminar or turbulent?
tas38 is offline   Reply With Quote

Old   July 1, 2020, 12:28
Default
  #3
New Member
 
berat
Join Date: Jun 2020
Posts: 6
Rep Power: 6
beratcagan is on a distinguished road
thank you, I never concidered the reynolds number. But i calculated the Reynolds number by hand and get 2.17. so i should still have laminar flow.

Quote:
Is the solution converging for the high-Re case?
to be hornest i didnt know what convergence in fluid dynamics are up till now, after looking it up. i can tell you that i didn't implement anything to help with convergence. so my simulation probably doesn't converge for high-Re cases.
Do i have any other options of assuring convergence besides implementing relaxationFactors?

Quote:
Is the flow modeled as laminar or turbulent?
i am using the piso solver and i have the simulationType set to laminar in the turbulenceProperties file.
beratcagan is offline   Reply With Quote

Old   July 1, 2020, 13:14
Default
  #4
Senior Member
 
Troy Snyder
Join Date: Jul 2009
Location: Akron, OH
Posts: 220
Rep Power: 19
tas38 is on a distinguished road
Using incompressible piso in openfoam, the only fluid property to specify is kinematic viscosity. What value did you prescibe for this and the inlet velocity for the faulty mesh case? I
tas38 is offline   Reply With Quote

Old   July 1, 2020, 13:22
Default
  #5
New Member
 
berat
Join Date: Jun 2020
Posts: 6
Rep Power: 6
beratcagan is on a distinguished road
Code:
transportModel  Newtonian;

nu              nu [ 0 2 -1 0 0 0 0 ] 1;

CrossPowerLawCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1;
    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1;
    m               m [ 0 0 1 0 0 0 0 ] 1;
    n               n [ 0 0 0 0 0 0 0 ] 1;
}

BirdCarreauCoeffs
{
    nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1;
    nuInf           nuInf [ 0 2 -1 0 0 0 0 ] 1;
    k               k [ 0 0 1 0 0 0 0 ] 0;
    n               n [ 0 0 0 0 0 0 0 ] 1;
}
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
	value		uniform (0.62e-3 0 0);

    }
    outlet
    {
        type            zeroGradient;
    }
    walls
    {
        type		noSlip;
    }
}
but i should mention here that i am am also using this line of code here. so the units are in micro. e.g μg, μm etc.

Code:
DimensionedConstants
{
    unitSet             micro;

microCoeffs
{
universal
{
c c [ 0 1 -1 0 0 0 0 ] 2.99792e+8; // speed of light in vacuum (µm/µs)
G G [ -1 3 -2 0 0 0 0 ] 6.67429e-20; // gravitational constant (µm^3/(pgµs^2))
h h [ 1 2 -1 0 0 0 0 ] 6.62607e-13; // Planck's constant (pgµm^2/µs)
}
electromagnetic
{
e e [ 0 0 1 0 0 1 0 ] 1.602176634e-4; // elementary charge (pC)
}
atomic
{
me me [ 1 0 0 0 0 0 0 ] 9.10938e-16; // electron mass (pg)
mp mp [ 1 0 0 0 0 0 0 ] 1.67262e-12; // proton mass (pg)
}
physicoChemical
{
mu mu [ 1 0 0 0 0 0 0 ] 1.66054e-12; // atomic mass unit (pg)
k k [ 1 2 -2 -1 0 0 0 ] 1.38065e-8; // Boltzman constant (pgµm^2/µs^2K)
}
standard
{
//- Standard pressure [pg/µmµs^2 ]
Pstd Pstd [ 1 -1 -2 0 0 0 0 ] 1e2; // 1 pg/µmµs^2
//- Standard temperature [degK]
Tstd Tstd [ 0 0 0 1 0 0 0 ] 298.15; // should be same as in SI unit system
}
}

}
beratcagan is offline   Reply With Quote

Old   July 1, 2020, 13:28
Default
  #6
New Member
 
berat
Join Date: Jun 2020
Posts: 6
Rep Power: 6
beratcagan is on a distinguished road
kinematic viscosity is 1
velocity is (0.62e-3 0 0)

i should also metion, that i have been using the "unitSet micro" command which is supposed to change the units into micro. so its not kg,m etc anymore, but rather μg, μm etc.
beratcagan is offline   Reply With Quote

Reply

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
foam-extend-4.1 release hjasak OpenFOAM Announcements from Other Sources 19 July 16, 2021 06:02
[snappyHexMesh] Layers:problem with curvature giulio.topazio OpenFOAM Meshing & Mesh Conversion 10 August 22, 2012 10:03
OF 1.6 | Ubuntu 9.10 (64bit) | GLIBCXX_3.4.11 not found piprus OpenFOAM Installation 22 February 25, 2010 14:43
mass flow rate and size of mesh zandi CFX 9 February 21, 2010 05:31
Changing mesh size for independency test az_f ANSYS Meshing & Geometry 1 August 27, 2009 14:55


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