
[Sponsors] 
March 8, 2016, 05:51 
Solve the Navier stokes equations

#1 
New Member
Rime
Join Date: Jun 2015
Posts: 28
Rep Power: 4 
Hi!
How to solve unsteady NS equations without using projection method (solve the coupled equation)? Thank you 

March 8, 2016, 07:47 

#2 
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 322
Rep Power: 16 
You can use fully coupled methods that linearize the mass and momentum components together and put those entries into a single sparse matrix and solve that repeatedly. Construction of particularly the mass conservation and pressure terms needs to be handled carefullyeither staggered mesh or RhieChow. That linearized system may need to be iterated a few times to converge. This leads to a very large sparse badly conditioned linear system and it generally takes a lot more work to solve than the individual fields in segregated methods. This is because the pressure equation for mass conservation is VERY hard to solve while the momentum equations are fairly easyso the full system ends up stagnating because of the twitchy pressure behavior. It can be very useful for some problems with complicated physics, but for singlephase incompressible flow, this is generally NOT recommended especially for transient flows.
You can also try to use NewtonKrylov methods with some underlying preconditioner (say an NS projection method!) to converge to the fullnonlinear system solution at each timestep. This may be recommended for some difficult problems or for very large timesteps, but again, probably not for transient incompressible NS. Projection methods are popular because they are fast and accurate enough, especially if the timestep must be small for accuracy reasons (as in LES/DNS).. And by "projection methods" here, I mean general fractionalstep noniterative methods. Note that all SIMPLE/PISOtype schemes for incompressible flow involve a projection step of sorts that solves for a pressure(like) variable and uses it to make the mass flux (aka velocity) field conservative. The projection process of rendering the mass flux field conservative iterationbyiteration or timestepbytimestep is a central portion of almost every primativevariable incompressible CFD solver that I am aware of. 

March 8, 2016, 09:35 

#3 
Senior Member
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 109
Rep Power: 11 
If you want to try something really different that doesn't involve projectionbased mixed methods, look at several FEM codes on the cfd wiki using a "pressureless" method.
cfdonline.com/Wiki/Source_code_archive__educational 

March 8, 2016, 09:59 

#4 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 3,250
Rep Power: 37 
If you want to have some idea of the general coupled problem that can be then decomposed in fractional steps, I suggest
https://www.researchgate.net/profile...3fd6000000.pdf http://faculty.nps.edu/fxgirald/Home...t_JCP_2002.pdf 

March 11, 2016, 05:29 

#5 
New Member
Rime
Join Date: Jun 2015
Posts: 28
Rep Power: 4 
thank you for the interesting answers and documents.
I have a question, what's the problem that i can have if I use for spatial discretization "collocation in the same node" (velocity and pressure)? 

March 11, 2016, 07:05 

#6 
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 322
Rep Power: 16 
The problem is generally called "pressurevelocity coupling" or "evenodd decoupling." This is a problem where the pressure and velocity fields can decouple on collated grids and exhibit celltocell oscillations. On 2D Cartesian grids, this results in a checkerboardtype oscillations. It occurs because the central differencing approximation of the pressure gradient is approximated as (p[i+1]  p[i1])/ (2*dx). The p[i] value does not enter the computation and so the oddvalues of i and the even values i can (and do) drift apart, often catastrophically.
The solution to this is either staggered mesh, whereby the pressure terms across a momentum cell can be computed as (p[east]  p[west])/dx where east and west of the pressures on the faces of the momentum cell. Or, by using RhieChow to interpolate the discrete momentum relation (not including the pressure term) from the two neighboring cells to their shared face and evaluating the face velocity/massflux using that momentum interpolation and the pressure gradient term as (p[1]  p[0])/dn where 0 and 1 are the cells on the positive and negative side of the face and dn is the distance between the centroids of those two cells. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Write a solver to solve the FithighNagumo equations  Alex1418  OpenFOAM Programming & Development  0  December 14, 2014 18:24 
MATLAB Finite Volume Unstructured Triangular Navier stokes  Mh.R  Main CFD Forum  0  October 18, 2011 06:06 
Navier Stokes equations in rotation frame..?  vinayender  Main CFD Forum  2  December 1, 2009 01:12 
Navier Stokes Equations  J.J.  Main CFD Forum  2  June 29, 2000 09:31 
Computational complexity of Navier Stokes equations  Marco Ellero  Main CFD Forum  5  May 5, 1999 21:07 