 So the next session will be taken by Mr. Krishnakant who is a research scholar from IIT Hyderabad and presently he is working with us in the open form team of 4C IIT Bombay. The session will be on Laplacian form plus conduction, ICO form plus scalar transport form. So Krishnakant over to you. So okay, I welcome you all in this day 3 session 5. This session will be on Laplacian form and it is used for conduction problems. And what Laplacian form actually is like it solves a simple one equation which is shown here. So a temperature equation with just the Laplacian term and a short term. So it can be used for a conduction problem. It is both steady, it can solve both steady and transient problems. It implies the non-orthogonal corrections but turbulence models are not incorporated in this. Okay, so the main junk of code is shown in here like this in the colored box. So this is the main part of the code. And so you might be like from our previous sessions you are very well familiar with the discretization and all. But I am still showing you for the sake of completeness. And then so the case we will be studying is using this Laplacian form we will be studied in a cavity problem where you can say it is a material slab metal slab. And the domain and the wall and all those will be same as the grid driven cavity problem as you are very familiar with this grid driven cavity case. And we will be specifying the time gradient and Laplacian schemes as shown here. Like time will be this far as we are going earlier backwards scheme and grid T will be Gaussian here Laplacian will be the corrected schemes. The material diffusion property will be 2 10 to the power minus 3 meters per second. And the boundary condition will also be able to take very simple boundary conditions. Will be grid state condition. But later on I will be showing you how you can add a simple source term inside this slab in the half of the domain slab. So without further ado, let's start the process. Okay. So I am using WSL. And so this is my first one to spin. So let's start with I need to source my open pool. You guys may not need to because you have already in your bash RC file. So can I get a yes or something like all are you are with me till here and started with the terminal. It's a hands on session. So you can type in the chat box why yes or no you can type. So okay. So this in this case we will be showing you how to build a case file from scratch. So let's start. So but we will not write or obviously all the files of ourselves. So just let us copy a simple file from the home tutorials for the laplacian phone. So I think you all are familiar with this process like you can use the tab of button so that you can get like the whole part of the code. So I am copying this file from home tutorials. A basic laplacian form. And here it is the flange. And we are copying it at the run folder. Okay. And let's say we rename it to a cavity. It's fine. Are you all done till here. So let us simultaneously see something like so that you can also check it. You can check these things whether it is working or not in this. Let's go to explorer.exe so that you know where what is going on. So let's say I am here. Side by side I will be showing you that in my open form this is my user directory. This is my run directory. So I have copied this file. Okay. So now let's go inside. Yeah. Hello. Sir. No such files. No such file directory. Maybe you have some problem with the like the upper command cp command with the location. You can, I mean are you all getting this problem? Most of you are getting this problem? Yes sir. I am also getting the same problem that file is not showing. Okay. So you are in the home directory right? You all are in the home directory right? If you do this ls you are in the home directory where you have you might be having this open form folder here. Yes sir. Okay. So why don't you one of you just share the screen so that you are facing the problem I will see. Yes sir. Okay. Someone share. So yeah it worked right? It worked. Very, I mean you don't see it but I mean we will let it go there where you can see the file. It worked. I don't see that there any error. So okay it's showing. So stop sharing. I will share my screen and continue. So let's go where you can see that you have copied just now. So we can go to the form run directory where you have copied your file. So in the run directory you can see your integrity that you have just now copied. We are copying it in the home run directory. Are you all able to see? Can I get like some of you are able to see some of you are in here. Okay. So let's go. Hello sir. Yes. Sir actually we don't have it in file. You don't you might have not copied. I mean the copy thing you need to do the upper upper file the second line in the screen. And I think Abhushan has shared the command so you can copy it from the text box. Okay. So let's go into the heated cavity and let me make sure that you have these files. Zero all clean constant flag answers system. So you might be seeing this dot answers file first time. So it's a I mean if you want to run this kind of. So you can create convert the answers files to open form mesh using just answers to form command. And we won't be dealing I mean if time permits we will deal with it later. So but in here for our file for our case we will be we need a block mesh file for from the lead driven cavity. So let's copy this block mesh file from your from the from the lead driven cavity. So for that let's go to the your deep driven cavity problems. The forms for tutorials and the lead driven cavity problem is in the compressible cavity cavity. And in their system directly and inside that block mesh. Since I am already inside this heated cavity. So I will copy it in the system. And we are I'm not renaming it so it will copy it with the name of blogmas itself. And how to state. Callity poverty blogmas didn't know such file about it. The icofoam is missing here. In front of cavity. So. Are you guys able to follow up this command line. I copied the block much dick from the lead driven cavity case. Done till here. Yes, no. So okay let's then copy. So we need to edit this new block for the new block mesh in it all the pressure all the boundary conditions. Mainly the bonding conditions are for the previous case for the flange case which is that is what we have copied in the tutorial. So we need to change the bonding conditions accordingly for our our case. So for that we will go to. I did the file in the. Zero folder. So here you can see the boundary. Are you all able to see it. So the boundary field are for the previous case of the flange case. So but we are doing the lead driven cavity problem. So our boundary will be for the lead driven cavity case which is moving on. And she was. And then by default let us say I keep on the file faces. Empty. So we will give the fixed. Condition that is fixed by permission and the value let's say keep it for the internal field as. 300. For the moving wall that is the top wall let's say. High temperature 350 for all the sidewalls let's say. 300 again. Okay. So I said this. Please can you show the file again. Define. So we are changing the face name like the boundary face in the block magic we defined the boundary patch name. So in the boundary condition then it will ask for the body condition for that patch. So that is what we are defining here. Since in the lead driven cavity you are familiar that there are. There are three kind of faces that is moving walls, face walls and front and back. So in front and back I mean one way to define I mean. Other than all the moving wall and fix walls whatever with the other patch fields patch walls for all those it will treat as as empty. Just a generic way I have defined it otherwise you can in general people front and back also you can give the name. So it will be good to go from here further. People followed. Yes. Yes. Form tutorial incompressible from there it is showing missing destination file. Mission destination file. Use a library open. Can I share my screen. Okay. Mission destination file. Follow the whole I mean I cannot see that you have type the whole command. You are just in the CPE phone tutorials in principle the whole command you have to type isn't it. And there is no destination. So that's why this is giving that. Okay. Okay. So let's go ahead we have modified our temperature field according to our case. So we are we are no need to change the schemes and solvers so that's why we are good to go. So can you again just show the temperature. Okay thank you. Okay. So we are good to go. We can do the block mesh and we can do on the Laplacian. So one of you if you can tell me like what is happening in here. Why temperature is solved so many times. You can type in the chat box. And why it is taking 00 iterations. And what is it necessary or I can undo it. So let's say I want to I mean it's unnecessary why I need to solve it every time right. So let's say I will change. So what I need to do. Not bad. It's down in the solution if I go. So in here you see that they have turned on the non orthogonal correction thing. So make it zero. Then. I mean it will it will not do the non orthogonal correction it will just solve one time. So let's say to say some computation time I'm making it zero. Okay. So locate this thing is done. So let's see our result. I open. So I think you are by now you are familiar how to open your. Result home so. People go ahead and take your result home file. So at the end it should look like something like this. Okay. Although some of you might have got it right. Yes. No. Got the result. So okay. So the next step we will do what we will add the source term some source term in here. Okay. So for that what I will do I will copy one more thing which is my every option scheme. So for that maybe copy it from home. Tutorials. Compressible. Home. And then temperature. So. The wine setting. Constant. And that system. Okay. So I'm copying one every option five. This new file. I'm copying in my system directory. Any issue they are copying this new file. Everyone copied. What did you modify in FV solution. FV solution I modified the non-orthogonal correction to zero. But I have not run the case again. So it doesn't matter. And for this this is a simple case it won't change the solution also. As you can see it is taking zero iteration also. Since it was orthogonal mesh, block mesh. It doesn't matter. Non-orthogonal correction is not needed. People copied the FV options file. So let us see what this new file we have copied. So it's actually set up the source term in your solutions. So this is one way to set your source terms. So what is the case that we have caused the file that we have copied is for different case. So for our case I'm going to modify it. Okay. So I'm going to modify it also in a very generic manner. So this is very generic way to explain you what is inside this. So for you let's replace this file with this new lines. I think you will get the commands and all the text files. I mean someone will forward you this command line. You can copy paste from the chat box. This new line of code you need to add. Replace it in the FV option files. So is the code forwarded? Line's forwarded. I forwarded it. I think forwarded now. So okay then save this file. So what this line that we are coding we are saying is we are adding this scalar coded source term because it's temperature is scalar term. We are adding a scalar source term. And the field that we are modifying it is t name you can give anything. Let's say I'm giving xx. And so in the all the source term you need to add it is a volume which is a finite volume code. So we need volume also. And the source term we are defining the scalar field t source equation will repeat this equation for this field. It will modify it will take the source it will take the pointer of this equation source terms and that will be saved in the t source. And later on I'm adding some terms in this pointer source pointer for the temperature equation. So where I'm adding so in this you can see this is the cell loop that is present here. So for all the cells so here you can see the selection mode is all so for all the cells and in that cells where this this is the for xx square plus y square minus r square if less than 0 less than 1 then I'm adding the source term. Okay. So this is x square. So x is here msc component msc for all say cell level I this will go look over all the cells and for the I itself it has the coordinate vector xyz. So the component 0 is the xx component component 1 is the y component 2d case it is so we are adding x square plus y square minus r square less than 1 for all the cells satisfying that condition we are adding the source term. Okay. So let us say I save this. So we will rerun the case so what we will do we will clean the tutorial or cleaning it we have the command phone clean tutorials. So we need to do the block mesh again. People are following. Can you show that clean tutorial command once. Thank you. So the application phone if I run now. So it is reading that the option file. And now you can see there is no non orthogonal directions just out one time with one iteration and let's see the new results. Sir. Yes. After the command phone clean tutorials. We have to do block mesh. Yes. Sir but after doing block mesh it is showing the form existing. Exiting. Yes. You are in the heated cavity file folder only. And have you previously done the block mesh thing previously without this idea of something. At that time you have run the case. Yes, previously I run. Okay. Share yours. Do this phone clean tutorials once again for me. Do the block mesh. It's been right. I think you previously missed that. So okay. We are good to check it now. Our new solution. Anyone till here followed up. I mean run the cases all the cases. See the new result. So you can see the difference right here because of this. And we have run till time three so and 50 is the source time added. So around 150 temperature profile is developing. So all of you. Done. We still have any doubts regarding the solver or a scheme and anything. Any question. Of the solver. Sir. Yeah. Okay. So. Okay. Let's do this one more case of this. We can see the plan and this right. So we I will show you how you can do that. But we have already modified the temperature file and all. So we will what we will do come out of this. I will copy that case once again. Okay. Cp phone tutorials. Basic. Flash. Flash. Flash. Flash. Okay. Copy this file in here. The destination location in this. So. Yeah. Okay. So let's go to the plan. Okay. So let's go to the plan. Okay. And you will. We have again the same case files. We don't need to take any important condition and schemes because it's already set. Case file so. We will do that. This is for two poem and flange. We are we can scale our mesh also. So I'm scaling the. Ansys file to let's say. Money minus three. So it run properly. Applation form to you need to capital. To be the small. Capital. So let's say touch. For the new file I will load now. So this is how the new. Just answers file you that you have run. So this is one tutorial for how to you can convert the mesh of answers. No one having any doubt this. So what do you mean by scale. Scaling like the while making the this geometry may not one may not be taking care of like the dimension maybe meter or anything like something like that. So if if suppose while making the this geometry I have kept let's say this distance as one meter. So I have if I do the scaling then it will become one mm. Like if I do the scaling with one e-minus three. Actually you have the. All the yeah all the geometry will be multiplied with one e-minus three. You you you have made you have made this geometry by by different different dimension in answers and you. No actually it's come with the tutorial case. I have not. I mean in the open-phone tutorial case it comes with that. Okay. If you want to see the dimensions and all you can check here is what I see. This will tell you the dimensions all those dimensions are previously if you don't do the dimensions and you consider it might be 25 or something like that. One e-minus three. Okay. One thing I missed actually one of you can redo it and just tell me what happens. I haven't changed the transfer properties. So keep it to e-minus three. Because in the later case in the next session I will be I mean we can analyze our results like in the test in the previous case. You do without the source term with the same diffusivity to e-minus three to minus three. You can see how the results vary with the source term and then later when we will add the velocity field properly with iphone how our results will modify because of the velocity field. So it will be a good comparison. But so we miss it. So you can do it later on for your own sake. I don't think anyone any doubt. I think if no doubt then after answers to form you can take your results. If no doubt then maybe I can start the another session so that there we will have the spare time. It's a good idea. You can start the next one.