|
[Sponsors] |
Writing a Streamfunction-vorticity simulation in MatLab |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 11, 2017, 09:57 |
Writing a Streamfunction-vorticity simulation in MatLab
|
#1 |
New Member
Ed
Join Date: Nov 2017
Posts: 2
Rep Power: 0 |
I've been tasked with modelling this 2D incompressible viscous flow using these boundary conditions and for these values of U₁ and U₂ on a rectangular square grid with 193 by 129 grid points. I have been asked to model the steady-state flow at Re=200.
I’m not sure how I should go about this. Does anyone recommend a comprehensive guide or perhaps a similar example I’d be able to read through and reverse engineer? |
|
November 11, 2017, 11:38 |
|
#2 |
Member
AGN
Join Date: Dec 2011
Posts: 70
Rep Power: 14 |
U could refer this link
https://curiosityfluids.com/2016/03/...n-cavity-flow/ https://www.iist.ac.in/sites/default.../psi-omega.pdf https://in.mathworks.com/matlabcentr...en-cavity-flow and change that Boundary conditions as per ur requirement. |
|
November 11, 2017, 12:47 |
Vorticity equation
|
#3 |
Senior Member
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 11 |
There are two great books for this problem:
1. Computational Techniques for Fluid Dynamics 2: Specific Flow Problems, Fletcher 2. Computational Fluid Mechanics and Heat Transfer, Anderson With that said most literature references will use SOR as the iterative solver for to find the streamfunction. I would avoid using such a method as it takes much longer for a convergence solution. BiCGStab or GMRES are more appropriate solvers . |
|
November 11, 2017, 13:11 |
|
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
just a comment, out of the topic of this post...SOR can be slower in terms of number of iterations but if you consider the real computational time there is almost nothing fast as a SOR (that can be also easily parallelized) programmed for a standard second order discrete Laplace equation. You need to work on 3D problems with some hundreds of millions of nodes to see GMRES/BiCGStab being more convenient.
|
|
November 11, 2017, 16:49 |
BiCGStab
|
#5 |
Senior Member
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 11 |
Dr. Denaro,
I was speaking in terms of number of iterations for convergence, otherwise you are correct. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] refineWallLayer Error | Yuby | OpenFOAM Meshing & Mesh Conversion | 2 | November 11, 2021 12:04 |
writing a udf to obtain velocity and vorticity at each point | reza.erfani | Fluent UDF and Scheme Programming | 1 | June 23, 2017 18:44 |
Writing patches during dynamic mesh motion simulation | eelcovv | OpenFOAM Post-Processing | 1 | January 16, 2016 16:21 |
writing out pressure values from OF Simulation | gillers18 | OpenFOAM Post-Processing | 1 | January 15, 2016 20:36 |
restarting paused transient simulation using reactingFoam | JMDag2004 | OpenFOAM Running, Solving & CFD | 1 | August 10, 2015 11:15 |