Hi Everyone, I just wanted to say thank you to those of you who noticed the sign errors on Slides 6 and 7. The coefficients A21 and A23 should be - k L / Delta x. In general for the Laplacian operator, the coefficient of the diagonal (A22) should equal -1 times the sum of the other coefficients, and so the coefficients A21 and A23 should have a - sign! In Equation 1 and 2, Qsource should be negative, as the sign convention here is positive heat fluxes denote heat leaving the cell (so Qsource should be -Qsource as heat is entering the cell). For those of you who support the channel on Patreon, the PDF slides have been corrected there. For everyone else, I have pinned this post so that you are aware of the typos if you go through the derivation in your own work. Thanks again for watching the lectures and supporting the channel :) Aidan
I hope youtube provides a way to re-upload with corrections in slides to the video in a much quicker way. It helps everyone. This is because not everyone sees both videos and comments. @Aidan: It may help everyone.
Yes, me too! I haven't found a better way of doing this yet. It would be really good if there was a way of having a 'living lecture' that we can update with corrections and improvements. The best I can think of so far is to pin a comment with corrections / include the corrections in the description
@@fluidmechanics101 Hello there, I think that the explanation stated above could be rather misleading for someone. There are strict rules that should be followed during the calculation procedure. First, we need to specify a positive direction for the heat flux (let´s assume a positive direction for heat flowing INTO the element). Therefore, any heat flowing out of the element will have a minus sign. Then, we need to stick to the second rule which is the fact that heat flows from a higher temperature to a lower one. So, according to our specified diagram, the temperature differences should always be written as follows: (temperature at the tail of the arrow - temperature at the head of the arrow). The trick is that we can specify the directions of the arrows in the diagram arbitrarily as long as we stick to these rules. Hope it helps! Thanks for the great work!
It’s just great to see how you manage to unfold the less intuitive concepts into a super clear explanation based on understanding rather than learning definitions by heart. Keep it up man!
Your mastery over the topic is so great that everything challenging about CFD becomes easy to grasp thanks to your beautifully chosen examples during the lecture.
Good grief! What a perfectly wonderful explanatory video. After four years of engineering I didn't get it half as well as I got it from your presentation. Thank you!
Hello from Russia. it's really cool that you record such informative and explanatory videos, thank you very much, I really want you to release your videos more often, as they help to better understand the work of cfd code. I'm really looking forward to the second part about the residuals. thank you!
Very good talk, and well explained. I like the physical intuition, it is exactly the way I interpret residuals as well. For a typical incompressible solver, the residual normally represents how 'divergence free' your flow is. For our code, we actually use both L_2 and L_inf and have relations between them to assess convergence. I am sure you will touch on this in the next parts but felt I would just add the below: Another added approach I have seen to residual normals is volume weighting. Here, every cell's residual is multiplied by the volume of the cell, and then the domain volume is divided through instead of N. This allows for small stiff cells to not hold up a simulation. This can be useful on unstructured meshes. Yet another thing to do is to divide your residual by b, this normalises the residual in a way that makes the problem less dependent on the dimensionality of the simulation. In this case, the interpretation is more along the lines of how many orders do you want to drop the residual 'relative to the starting point'. This can provide more consistent convergence tolerances being applied.
This is great! Finally residuals are being addressed..went through some cfd manuals but there was never really a clear explanation about what or how..thanks for the awesome video. Can't wait to see the next one
Q: In practice, what do you think is the acceptable range of representative residual values to say that you have good results in your CFD? Some papers say around 10^-4 and below for all transport parameters. Can't wait for the next parts of your lecture.
Yes, 10^-4 seems to be pretty good in most CFD solvers (just from personal experience and the experience of others). But as always ... please check other quantities to check you have convergence (forces, moments, point monitors etc). Residuals are not enough by themselves!
Your lucid way of explanation is just awesome. When I have question like how can I get just a value for the residual that may represent residuals of each cells.? And there you explain. I do not even recognise that the video has ended. Thank you so much for your hardwork.
Eagerly waiting and excited to hear your talks they are great and help me out very much in my masters studies it an request can you please make some video on Aeroacoustics cfd simulation basics
Wallahi this guy explains things really well! He is a master! Please consider explaining the difference in using fvSolution in openfoam, which is faster GMAC or PCG solvers, smoothers and those complicated topics 😭
Thank you so much for the explanation. I'm used to follow what my lecturer does since he doesn't explain anything abt residual concept. Now I really understand it after watching this vid. Looking forward for Part II. Keep it up man!
Just a few days back was researching the residual outputs from solvers in OF! Thank you. Can you please also consider making videos on linear solvers and compressible solves.
Thanks for this amazing video. Residuals topic can be very confusing at times and it’s difficult to interpret the residual line plots. Eagerly looking forward for the subsequent parts. Also it would be helpful if you could address oscillating residual plots which I had once encountered in my CFD simulation.
Thanks Man for the informative work, this kind of work over youtube may not be benifitial to you compared to the youtubers who provide silly content but, the gold stays gold and whatever else is whatever else. keep the good work up, we are excited for more !
Thank you for these lectures once again, you are such a good pedagog on this subject! There is a slight error on the slide at 4:15 where you go through the numerical calculation. The T2 term should be negative overall, following all the algebra.
It’s just disappointing, that there are not all parts of your already out! Awesome how easy cou can deliver that complicated stuff in an easy manner. Thanks a lot!
Very good video! I am looking forward to watching the rest. Hoping that you will address oscillating residuals for steady-state computations and if it is the end of the world if it happens. I have also read that bad residuals do not always mean a bad solution. It is often better to track variables of interest, such as drag, flow rate, heat transfer rate, etc.
Really usedful video! Thank you. In case of SIMPLE algorithm (say in OpenFOAM), what does the pressure residual physically mean? Is it the imbalance of the pressure Poisson equation? In which case, would it represent an imbalance in force?
This is actually a fairly tricky question. In classical implementations of the SIMPLE algorithm (that solve a pressure correction equation rather than an equation for the pressure), the residual directly represents the mass error in a given cell (i.e the sum of the mass flow rates across the faces of the cell equals the residual which should tend to zero). However in OpenFOAM, an equation is solved for pressure rather than the pressure correction. Hence, it is slightly different but I still tend to think of it as the mass error in a given cell, just because it is easier for my understanding (Notice that it is a scalar equation, not a vector equation, so the residual cannot represent a vector quantity like force. Mass however is a scalar)
As usual, great content! Thank you a lot, I’m a CFD user but thanks to you I’ve been building up some of the concepts that always slipped away from my understanding.
Thank you sou much for explaining these concepts in a simple and understandable way. Also can you please tell ne how you make your slides look so good? How do you make thise images for the cells and plots?
It's was a very useful class to understand the residual physics. You explain clearly all terms in the equations. But I have a question, what happen with the scaled residuals that are used in Fluent? does it have the same meaning?