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

Solve the Navier stokes equations

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

Like Tree2Likes
  • 1 Post By FMDenaro
  • 1 Post By mprinkey

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 8, 2016, 05:51
Question Solve the Navier stokes equations
  #1
New Member
 
Rime
Join Date: Jun 2015
Posts: 13
Rep Power: 10
Rime is on a distinguished road
Hi!
How to solve unsteady N-S equations without using projection method (solve the coupled equation)?

Thank you
Rime is offline   Reply With Quote

Old   March 8, 2016, 07:47
Default
  #2
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
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 carefully--either staggered mesh or Rhie-Chow. 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 easy--so 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 single-phase incompressible flow, this is generally NOT recommended especially for transient flows.

You can also try to use Newton-Krylov methods with some underlying preconditioner (say an NS projection method!) to converge to the full-nonlinear 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 fractional-step non-iterative methods. Note that all SIMPLE/PISO-type 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 iteration-by-iteration or timestep-by-timestep is a central portion of almost every primative-variable incompressible CFD solver that I am aware of.
mprinkey is offline   Reply With Quote

Old   March 8, 2016, 09:35
Default
  #3
Senior Member
 
Jonas T. Holdeman, Jr.
Join Date: Mar 2009
Location: Knoxville, Tennessee
Posts: 128
Rep Power: 18
Jonas Holdeman is on a distinguished road
If you want to try something really different that doesn't involve projection-based mixed methods, look at several FEM codes on the cfd wiki using a "pressureless" method.

cfd-online.com/Wiki/Source_code_archive_-_educational
Jonas Holdeman is offline   Reply With Quote

Old   March 8, 2016, 09:59
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
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
Rime likes this.
FMDenaro is offline   Reply With Quote

Old   March 11, 2016, 05:29
Default
  #5
New Member
 
Rime
Join Date: Jun 2015
Posts: 13
Rep Power: 10
Rime is on a distinguished road
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)?
Rime is offline   Reply With Quote

Old   March 11, 2016, 07:05
Default
  #6
Senior Member
 
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25
mprinkey will become famous soon enough
The problem is generally called "pressure-velocity coupling" or "even-odd decoupling." This is a problem where the pressure and velocity fields can decouple on collated grids and exhibit cell-to-cell oscillations. On 2D Cartesian grids, this results in a checkerboard-type oscillations. It occurs because the central differencing approximation of the pressure gradient is approximated as (p[i+1] - p[i-1])/ (2*dx). The p[i] value does not enter the computation and so the odd-values 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 Rhie-Chow 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.
Rime likes this.
mprinkey 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
Write a solver to solve the Fithigh-Nagumo 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 07: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 10:31
Computational complexity of Navier Stokes equations Marco Ellero Main CFD Forum 5 May 5, 1999 22:07


All times are GMT -4. The time now is 18:39.