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

The problem of bed load transport rate based on FE3.2-famesh

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2018, 00:48
Smile The problem of bed load transport rate based on FE3.2-famesh
  #1
New Member
 
LXJ
Join Date: Apr 2015
Posts: 12
Rep Power: 5
liuxiaojian2015 is on a distinguished road
Hi all,
when i calculated the bed load transport rate (qbi=q0*τ/τi - C*|q0|dη/dxi) component in the i direction by using OpenFOAM famesh, i found the result of dη/dxi in two sides in y direction (two sides patch: symmetryPlane) is wrong by simulating the scour of pileline under current and using five layer grids. The result show the value is zero in internal face ceters in y direction, however it is very big in two sides. Please help me if someone know how to slove it. Thank you.

The code:
//########################### bed load transport rate ###########################//
hA=mesh.boundaryField()[patchID].faceCentres.z();
edgeScalarField absWLEdge = fac::interpolate(hA);
areaVectorField gha = fac::edgeIntegrate(absWLEdge.mesh().Le() * absWLEdge);
// or //areaVectorField gha = fac::grad(hA);

areaVectorField qbA = q0A*shearN - 1.5 * mag(q0A) * gha;

Best,

LXJ
Attached Images
File Type: jpg grad(hA).y().jpg (77.1 KB, 14 views)
File Type: jpg pipeline.jpg (85.9 KB, 15 views)
File Type: jpg result.jpg (114.9 KB, 15 views)
File Type: jpg ???????.jpg (41.8 KB, 12 views)

Last edited by liuxiaojian2015; March 13, 2018 at 03:04.
liuxiaojian2015 is offline   Reply With Quote

Old   March 13, 2018, 06:01
Default
  #2
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 179
Rep Power: 8
Santiago is on a distinguished road
Quote:
Originally Posted by liuxiaojian2015 View Post
Hi all,
when i calculated the bed load transport rate (qbi=q0*τ/τi - C*|q0|dη/dxi) component in the i direction by using OpenFOAM famesh, i found the result of dη/dxi in two sides in y direction (two sides patch: symmetryPlane) is wrong by simulating the scour of pileline under current and using five layer grids. The result show the value is zero in internal face ceters in y direction, however it is very big in two sides. Please help me if someone know how to slove it. Thank you.

The code:
//########################### bed load transport rate ###########################//
hA=mesh.boundaryField()[patchID].faceCentres.z();
edgeScalarField absWLEdge = fac::interpolate(hA);
areaVectorField gha = fac::edgeIntegrate(absWLEdge.mesh().Le() * absWLEdge);
// or //areaVectorField gha = fac::grad(hA);

areaVectorField qbA = q0A*shearN - 1.5 * mag(q0A) * gha;

Best,

LXJ
Which are your bc's for velocity in the spanwise direction?
Santiago is offline   Reply With Quote

Old   March 13, 2018, 06:08
Default
  #3
New Member
 
LXJ
Join Date: Apr 2015
Posts: 12
Rep Power: 5
liuxiaojian2015 is on a distinguished road
Quote:
Originally Posted by Santiago View Post
Which are your bc's for velocity in the spanwise direction?
Hi Santiago,

Yes, you are right. The spanwise direction (i.e., x direction) is set as the velocity direction. Two sides are perpendicular to y direction, and parallel to x direction.

Best,

Liu
liuxiaojian2015 is offline   Reply With Quote

Old   March 13, 2018, 06:10
Default
  #4
Senior Member
 
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 179
Rep Power: 8
Santiago is on a distinguished road
Quote:
Originally Posted by liuxiaojian2015 View Post
Hi Santiago,

Yes, you are right. The spanwise direction (i.e., x direction) is set as the velocity direction. Two sides are perpendicular to y direction, and parallel to x direction.

Best,

Liu
Ok,

Two questions:

1. Which Boundary conditions you set for velocity in the domain?
2. Do you start your simulation from a solution obtained by not moving the grid/potential?
Santiago is offline   Reply With Quote

Old   March 13, 2018, 06:31
Smile
  #5
New Member
 
LXJ
Join Date: Apr 2015
Posts: 12
Rep Power: 5
liuxiaojian2015 is on a distinguished road
Quote:
Originally Posted by Santiago View Post
Ok,

Two questions:

1. Which Boundary conditions you set for velocity in the domain?
2. Do you start your simulation from a solution obtained by not moving the grid/potential?
Hi, Santiago

Thank you for you reply. First, in order to explain this answer, i show a picture for the detailed model setup (plese ignore my clumsy drawing, ), the Boundary conditions for the velocity is set in the left side of model. Second, i attempted the simulation by not moving the mesh motion, the result of grad(hA).y() is acquired in such condition and it is showed again in the attachment.

Best,

Liu
Attached Images
File Type: jpg SETUP.jpg (19.2 KB, 13 views)
File Type: jpg grad(hA).y().jpg (77.1 KB, 14 views)
liuxiaojian2015 is offline   Reply With Quote

Old   March 13, 2018, 06:41
Smile
  #6
New Member
 
LXJ
Join Date: Apr 2015
Posts: 12
Rep Power: 5
liuxiaojian2015 is on a distinguished road
Quote:
Originally Posted by Santiago View Post
Ok,

Two questions:

1. Which Boundary conditions you set for velocity in the domain?
2. Do you start your simulation from a solution obtained by not moving the grid/potential?
I also attempted the simulation by not considering the current velocity, i.e., i set the velocity is about zero. This value will not induce the sediment transport. Thus, no bad results were showed.
liuxiaojian2015 is offline   Reply With Quote

Old   March 13, 2018, 07:05
Smile
  #7
New Member
 
LXJ
Join Date: Apr 2015
Posts: 12
Rep Power: 5
liuxiaojian2015 is on a distinguished road
Quote:
Originally Posted by Santiago View Post
Ok,

Two questions:

1. Which Boundary conditions you set for velocity in the domain?
2. Do you start your simulation from a solution obtained by not moving the grid/potential?
I am so sorry to trouble you again. I give a error introdution several minutes ago. When i use the smaller current velocity, the bed surface gradient is also wrong, please see the attachment. The error also exist in the patch junction. Is this indicate the fac::grad() is wrong in patch junction?

I need emphasize the gha.x() and gha.z() keep same variation for the five layer mesh. Just gha.y() show the error. Thank you.
Attached Images
File Type: jpg grad(hA).y()----small current velocity.jpg (57.3 KB, 9 views)

Last edited by liuxiaojian2015; March 13, 2018 at 09:26.
liuxiaojian2015 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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Multiphase liquid-solid spouted bed: pressure problem ghost82 FLUENT 7 November 10, 2013 13:12
Species transport model problem argeus FLUENT 7 March 12, 2013 04:50
Specifying volume/Mass fraction in Species Transport for problem without inlet tanjinjack FLUENT 1 October 9, 2011 05:42
Fixed Bed Gasifier, Species Transport Problem therandomestname FLUENT 11 May 8, 2011 10:20
moving load problem aditi chauhan ANSYS 0 October 9, 2010 05:47


All times are GMT -4. The time now is 16:05.