|
[Sponsors] |
|
December 5, 2000, 09:42 |
results of my MAC solver
|
#1 |
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!!! |
|
December 5, 2000, 10:33 |
Re: results of my MAC solver
|
#2 |
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 |
|
December 29, 2000, 15:36 |
Re: results of my MAC solver
|
#3 |
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 |
|
January 5, 2001, 15:24 |
Re: results of my MAC solver
|
#4 |
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 |
|
January 5, 2001, 17:38 |
Re: results of my MAC solver
|
#5 |
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 |
|
|
|
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 |