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

Projection method and mesh free

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 13, 2015, 18:34
Default Projection method and mesh free
  #1
New Member
 
Rime
Join Date: Jun 2015
Posts: 13
Rep Power: 10
Rime is on a distinguished road
I have a problem in my code (diverge)
to solve navier stokes equation, i use mesh free method (uniform) and projection method; in the boundary condition I have Neumann boundary for pressure and for intermediate velocity u^*=2 and u^* =0,I use U^*.\hat{n} = U^{n+1}.\hat{n} but I don't update the velocity in boundary ie i give all time u=2 and 0 is correct???
I don't know what's the problem,
help please (sorry for my english)
Rime is offline   Reply With Quote

Old   June 13, 2015, 19:01
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,756
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rime View Post
I have a problem in my code (diverge)
to solve navier stokes equation, i use mesh free method (uniform) and projection method; in the boundary condition I have Neumann boundary for pressure and for intermediate velocity u^*=2 and u^* =0,I use U^*.\hat{n} = U^{n+1}.\hat{n} but I don't update the velocity in boundary ie i give all time u=2 and 0 is correct???
I don't know what's the problem,
help please (sorry for my english)

Are you using an implicit method for diffusion in the momentum?

if the physical u at the boundaries is steady you can let the same value but u* = u_n+1 is first order accurate in time (read the paper of Kim and Moin on JCP).
However, that should not cause the divergence of the solution therefore check in the code for some bugs.
FMDenaro is offline   Reply With Quote

Old   June 13, 2015, 19:09
Default
  #3
New Member
 
Rime
Join Date: Jun 2015
Posts: 13
Rep Power: 10
Rime is on a distinguished road
I use Crank-Nicolson scheme for convective and diffusive terms
Rime is offline   Reply With Quote

Old   June 14, 2015, 05:15
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,756
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rime View Post
I use Crank-Nicolson scheme for convective and diffusive terms

how do you treat the non-linear term in the matrix?
FMDenaro is offline   Reply With Quote

Old   June 14, 2015, 06:33
Default
  #5
New Member
 
Rime
Join Date: Jun 2015
Posts: 13
Rep Power: 10
Rime is on a distinguished road
thank you for your reply.
I supposed that the term u_n+1 is transported by u_n
Rime is offline   Reply With Quote

Old   June 14, 2015, 06:39
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,756
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Rime View Post
thank you for your reply.
I supposed that the term u_n+1 is transported by u_n

so you are linearizing at first order...but the method should be second order accurate in time..

I suggest:

1) Use CN only for the diffusion and the second order Adams-Bashforth for the convective term
2) Check if you satisfy the continuity constraint in each cell at machine accuracy. You can see that also after 1 time step
FMDenaro is offline   Reply With Quote

Old   June 14, 2015, 06:54
Default
  #7
New Member
 
Rime
Join Date: Jun 2015
Posts: 13
Rep Power: 10
Rime is on a distinguished road
thank you I'll try that.
can you give me more information on suggestion 2)
Rime is offline   Reply With Quote

Old   June 15, 2015, 15:43
Default
  #8
New Member
 
Rime
Join Date: Jun 2015
Posts: 13
Rep Power: 10
Rime is on a distinguished road
Hi Sir
I have a variable density, to solve the pressure equation I use radial basis function,( its the mesh method that I use)
What is your advice?
Thank you
Rime 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
[ICEM] Problem with using the Delaunay volume mesh method manahara ANSYS Meshing & Geometry 12 December 25, 2014 22:54
Issue with arbitrary Lagrangian-Eulerian method and mesh optimization mikolchon Main CFD Forum 0 August 16, 2014 16:50
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 07:42
[snappyHexMesh] Layers:problem with curvature giulio.topazio OpenFOAM Meshing & Mesh Conversion 10 August 22, 2012 10:03
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 19:10


All times are GMT -4. The time now is 06:45.