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

results of my MAC solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 5, 2000, 09:42
Default results of my MAC solver
  #1
Maciej Matyka
Guest
 
Posts: n/a
hello everybody!

at http://panoramix.ift.uni.wroc.pl/~maq/

you have a couple of gif images generated by my MAC solver. as you see - its not a good result. i dont know what happens, whad im doing wrong. size of the mesh contains 32x16 cells. fluid have 5*5 particles in water basin and 10*10 particles in drop. viscosity 0.01, g=-0.98, v0=-1.0.

what im doing wrong?

im using The MAC Method described by Harlow and Welch 1963 (Los Alamos paper)

/// maq

ps. maybe put here more pictures?

ps2. help!!!
  Reply With Quote

Old   December 5, 2000, 10:33
Default Re: results of my MAC solver
  #2
Maciej Matyka
Guest
 
Posts: n/a
sorry. changes:

now files on panoramix (see prev. mail) are generated for:

width = 400 height = 400

nx = 40 (40 cells in x dir.) ny = 40

gy = -0.1 gx = 0

vis = 0.01

and still it is nothing good...

buuuuuuuuuu...

/// maq
  Reply With Quote

Old   December 29, 2000, 15:36
Default Re: results of my MAC solver
  #3
Krash
Guest
 
Posts: n/a
Hey Matyka, I am also working on solving a problem via MAC method. The problem I have experienced are stability related: I had to decrease time step, pick a fluid with a high viscosity,and let it run for many iterations. then, it finally converges...kind of. The convergence is "iffy" at best. I found this little tidbit on the Internet after struggling for weeks:

MAC: Marker-And-Cell

To treat incompressible, free surface flows, the MAC method was developed by Harlow and Welch [Harlow et al, 1965] as a variation of the PIC method but treating applications that extended beyond those addressed by the vorticity*stream-function method. The MAC method was the first successful technique for incompressible flows. Particles were used as markers to locate the material in the mesh and, consequently, to define the location of the free-surface. The MAC method had the advantage of a more compact finite difference stencil and tight coupling between the pressure and velocity fields. To treat the fluid incompressibility, a solution to the Poisson equation for the pressure was used. This was in contrast to later methods that solved the coupled velocity-pressure equations, as discussed by Viecelli [Viecelli, 1969]. Although the solution of Poisson's equation was numerically simple, the specification of the velocity boundary conditions were not straightforward. There was some controversy at the time about the relative stability of the MAC method, and this was resolved in the now-classic paper by Hirt [Hirt, 1968Hirt, C.W., J. Comp. Phys., Vol. 2, p. 339, 1968.], in which he showed that the MAC method is unstable with centered momentum advection unless the viscosity is sufficiently large. This work was the precursor of the modern truncation error subtraction analysis. This controversy illustrated the T-3 approach: the development was always on the physics, with limited application of mathematical analysis of, e.g., convergence and stability properties. The MAC method is still in use and has profited from the added efficiency of modern conjugate gradient schemes for solving the Poisson equation.

Krash
  Reply With Quote

Old   January 5, 2001, 15:24
Default Re: results of my MAC solver
  #4
Arturo
Guest
 
Posts: n/a
Hi, I'm a newbie and I would like to ask you about the existing codes for the MAC solver.

My advisor thinks that this method can be useful for my project about porous media.

Thanks

Arturo
  Reply With Quote

Old   January 5, 2001, 17:38
Default Re: results of my MAC solver
  #5
Krash
Guest
 
Posts: n/a
Hey Arturo, I recently coded (via Matlab) a program to solve the flow between two solid flat plates. The technique of "artificial compressibility" and the MAC method were used. I don't know if it will help with your project, but I can definitely send it to you.

The big points using this method as I see them are the damn staggered grid and deciding whether to use Poissons equation for pressure or another explicit method. Also, make sure all of the dang values are getinng calculated at the correct times.

Let me know if I can help.

Krash
  Reply With Quote

Reply

Thread Tools
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
OpenCL linear solver for OpenFoam 1.7 (alpha) will come out very soon qinmaple OpenFOAM Announcements from Other Sources 4 August 10, 2012 11:00
chtMultiregionFoam Solver msarkar OpenFOAM 1 February 24, 2010 02:56
Getting too many iterations by velocity solving (aborting). Changing U - Solver? suitup OpenFOAM Running, Solving & CFD 0 January 20, 2010 07:45
switching from coupled solver to segregated Oz FLUENT 2 November 8, 2006 16:02
MAC Solver Source Code Page Updated Maciej Matyka Main CFD Forum 0 July 6, 2001 06:41


All times are GMT -4. The time now is 11:31.