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

Writing a Streamfunction-vorticity simulation in MatLab

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 11, 2017, 09:57
Default Writing a Streamfunction-vorticity simulation in MatLab
  #1
New Member
 
Ed
Join Date: Nov 2017
Posts: 2
Rep Power: 0
Woody27327 is on a distinguished road
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?
Woody27327 is offline   Reply With Quote

Old   November 11, 2017, 11:38
Default
  #2
Member
 
AGN
Join Date: Dec 2011
Posts: 70
Rep Power: 14
arungovindneelan is on a distinguished road
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.
arungovindneelan is offline   Reply With Quote

Old   November 11, 2017, 12:47
Default Vorticity equation
  #3
Senior Member
 
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 11
selig5576 is on a distinguished road
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 .
selig5576 is offline   Reply With Quote

Old   November 11, 2017, 13:11
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by selig5576 View Post
I would avoid using such a method as it takes much longer for a convergence solution. BiCGStab or GMRES are more appropriate solvers .
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.
selig5576 likes this.
FMDenaro is offline   Reply With Quote

Old   November 11, 2017, 16:49
Default BiCGStab
  #5
Senior Member
 
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 11
selig5576 is on a distinguished road
Dr. Denaro,

I was speaking in terms of number of iterations for convergence, otherwise you are correct.
selig5576 is offline   Reply With Quote

Reply


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
[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


All times are GMT -4. The time now is 14:15.