Example 03: Dropping object inside a tank
To calculate the gravity driven flow when a density difference occurred in vertical direction such as falling a heavy liquid from the top.
Creation of the calculation grid and setting the initial conditions
As explained in the other examples and the introduction, create the grid using, [Grid], [Select Algorithm to Create Grid] and then select [Grid Generator for Nays3DV]. Then the grid creation window will appear.
In grid creation window, give channel shape parameters as shown in Figure 48.
Then we can give Chanel bed condition. As here we use the default condition flat(no bar), no modifications are needed.
If new grids are added or width is varied it is possible to set them. As in this example no grids added and no width variations, no modifications are needed in them.
Initial water surface profile tab is used to give downstream depth, water surface slope and initial water surface purtavation. It can be seen as shown in Figure 49. After giving all the parameters click on [Create Grid].
Then the grid is created and a confirmation message box will appear asking to map the geographic data as shown in Figure 50. and click on [Yes].
This will map the geographic data to the grid and the mapped grid can be seen as shown in Figure 51.
Save the project with [File], [Save as .ipro] or [Save as Project] as in Figure 52.
Setting the calculation conditions and simulation
Set the calculation conditions with, [Calculation Condition], [Setting].
Calculation condition window will open.
Set computational parameters as shown in Figure 53.
Then give hydraulic boundary conditions. Since the boundary conditions are closed boundaries , boundary condition window is inactive as shown in Figure 54.
Then give initial and Boundary concentrations as shown in Figure 55.
Here background concentration is the concentration inside tank and perturbed concentration is concentration of the dropping item. To drop and go down the tank, the purturbed concentration should be higher than the background concentration.
The i, j, k start and end locations of the purturbed concentration should be given as shown in above figure.
Then the time and iteration parameters are give as shown in Figure 56.
Then give the physical parameters as given in Figure 57.
After setting the calculation conditions, save and close the calculation condition window.
Save the project again. Now start the simulation by [Simulation], [Run].
The end of simulation can be seen as shown in Figure 58.
Visualization of results
After calculation solver stopped, go to [Calculation Results], [Open new 3D Post-Processing Window].
Here it is easy to visualize as isosurfaces. Go to [Object Browser] - [iRIC Zone] - [Isosurfaces] and right click on [Isosurfaces]. Isosurface setting window will appear as shown in Figure 59.
Iso value need to be set and the min and max value can be used as a guide to set Iso value.
Colour also need to be set to be seen nicely. Default colour is white and if the colour setting is forgetted, isosurface can’t be seen.
The movement of concentration can be seen as shwon in Figure 60.
If need to visualize the concentration as countours, goto [Object Browser] - [iRICZone] - [Contours]. Right click on contours and [Add]. Contour Setting window will appear as shown in Figure 61.
Contours of concentration at the start and end can be seen as shown in Figure 62.