==================================================== [Example 1] Tracer transport in a straight channel ==================================================== Flow calculation by Nays2DH ================================= Select a solver ------------------- In the [iRIC start page] , select [Create New Project], and when the [Select Solver] screen appears, choose [Nays2DH iRIC 4.x 1.0 64bit] and click [OK] button. .. figure:: images/01/Select_Nays2dh.png :align: center :width: 600pt : Select Solver A windows with "Untitled - iRIC 4.x.x.xxxx [Nays2DH]" appears as :numref:`01_mudai`. .. _01_mudai: .. figure:: images/01/mudai.png :align: center :width: 100% : Untitled .. _01_lavel_koshi: Grid Generation ------------------- From the main menu of the screen, :numref:`01_mudai`, choose [Grid]->[Select Algorithm to Create Grid] as :numref:`Select_Alg`. .. _Select_Alg: .. figure:: images/01/Select_Alg.png :align: center :width: 100% : Select Algorithm to Create Grid In the [Select Grid Creating Algorithm] window, select [Simple Straight and Meandering Channel Creator] and click [OK] (:numref:`01_kanni`). .. _01_kanni: .. figure:: images/01/kanni.png :align: center :width: 600pt : Select Grid Creating Algorithm In the window of :numref:`01_koushi_1` , click "Channel Shape" and set [Select Channel Shape of the Main Part] as [straight channel], and other values as shown in :numref:`01_koushi_1`, then click [Create Grid]. .. _01_koushi_1: .. figure:: images/01/koushi_1.png :align: center :width: 600pt :Setting Channel Shape When the confirmation window appears as :numref:`01_koushi_3`, click [Yes] to generate the grid, then the computational grid is generated as :numref:`01_koushi_4` . .. _01_koushi_3: .. figure:: images/01/koushi_3.png :align: center :width: 400pt :Confirmation of mapping .. _01_koushi_4: .. figure:: images/01/koushi_4.png :align: center :width: 100% :Grid Generation Compete Setting of calculation conditions for flow by Nays2DH ------------------------------------------------------- The next step is to set the calculation conditions. From the menu bar, select [Calculation Conditions]->[Settings], then the [Calculation condition setting window] as :numref:`01_joken_1` appears. .. _01_joken_1: .. figure:: images/01/joken_1.png :align: center :width: 600pt :Calculation Condition Window As :numref:`01_joken_2`, in the [Group] of the [Boundary Condition], click [Edit] at the [Time series of discharge at upstream and water level at downstream]. Then the [Time series of discharge at upstream and water level at downstream] appears as :numref:`01_joken_3` . .. _01_joken_2: .. figure:: images/01/joken_2.png :align: center :width: 600pt : Boundary Condition .. _01_joken_3: .. figure:: images/01/joken_3.png :align: center :width: 600pt : Time series of discharge at upstream settings In :numref:`01_joken_3`, input [Time] and [Discharge] values, and click [OK] when you finish, and close this window. .. _01_joken_4: .. figure:: images/01/joken_4.png :align: center :width: 600pt :Time parameters Select [Time] and set parameters as :numref:`01_joken_4` and click [OK]. .. _res_Nays2DH: Flow calculation run by Nays2DH ---------------------------------- From the main menu, when you select [Simulation]->[Run], you will get the message like :numref:`01_warning` . Then, select [OK] and save the project with an appropriate name. At this time, do not save the project as an ipro file, but save it as a project. .. _01_warning: .. figure:: images/01/warning.png :align: center :width: 400pt :warning | A window as :numref:`01_jikko` is shown during the computation, and :numref:`01_keisan` appears when the computation is finished. | Then press [OK], and the computation is completed. .. _01_jikko: .. figure:: images/01/jikko.png :align: center :width: 100% :Window when the solver is running .. _01_keisan: .. figure:: images/01/keisan.png :align: center :width: 250pt :Computation completed .. note:: Whenever you finished the computation, select [File]->[Save] from the menu bar to save the results as :numref:`01_hozon` . This result is important for later analysis by GELATO. .. _01_hozon: .. figure:: images/01/hozon.png :align: center :width: 100% :Saving computational results Visualization of the calculated results ---------------------------------------------- After the calculation, select [Calculation Result] -> [Open New 2D Post-processing Window] to open the visualization window. .. _01_kekka_0: .. figure:: images/01/kekka_0.png :align: center :width: 100% : 2D Post-processing Window Velocity Vectors ^^^^^^^^^^^^^^^^^^^^ In the [Object Browser], put check marks in the boxes by [Arrow] and [Velocity], click Focus on [Arrow] and click the right mouse button [Properties]. Vector setting" window as :numref:`01_kekka_2` appears. Set the values in the red line and click [OK]. :numref:`01_kekka_6` is the depth-averaged velocity vector. Here, the velocity distribution is uniform under the constant flow condition. .. _01_kekka_2: .. figure:: images/01/kekka_2.png :align: center :width: 600pt : Vector Settings .. _01_kekka_6: .. figure:: images/01/kekka_6.png :align: center :width: 100% : Depth averaged velocity vectors Display Particle Movement ^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Uncheck "Vectors" in the Object Browser, and put check marks in "Particles" and "Velocity" ( :numref:`01_kekka_9` ) .. _01_kekka_9: .. figure:: images/01/kekka_9.png :align: center :width: 100% : Particles(1) Right click [Particle] and select [Properties] as :numref:`01_kekka_10` . .. _01_kekka_10: .. figure:: images/01/kekka_10.png :align: center :width: 100% : Particles(2) Set parameters for particle injection as shown in red box in :numref:`01_kekka_11` . .. _01_kekka_11: .. figure:: images/01/kekka_11.png :align: center :width: 250pt : Set particle parameters As shown in :numref:`01_kekka_12` , set time bar back to zero, and select [Animation]->[Start/Stop Animation] rom the main menu bar. Then the particle animation starts. .. _01_kekka_12: .. figure:: images/01/kekka_12.png :align: center :width: 100% : Start Particle Animation .. _01_kekka_13: .. figure:: images/01/nays2d_particle.gif :align: center :width: 100% : Particle animation by Nays2DH As can be seen in :numref:`01_kekka_13`, since the sub-grid scale turbulence is not included in the output velocity from the solver. It only shows very simple steady and uniform movement. Tracer Tracking by GELATO =========================== Starting GELATO ---------------- From the iRIC startup screen, select [New Project], and in the solver selection screen appears. Select "GELATO" and click "OK" ( :numref:`01_GELATO_kido` ). .. _01_GELATO_kido: .. figure:: images/01/GELATO_kido.png :align: center :width: 600pt : Selecting GELATO and Starting A window with [Untitled -iRIC 4.x.xxxx] [GELATO] appears, and the GELATO session is started. (:numref:`01_GELATO_openning` ) .. _01_GELATO_openning: .. figure:: images/01/GELATO_openning.png :align: center :width: 100% : Opening GELATO At this stage, the [Grid] in the [Object Browser] shows [No data] as shown in :numref:`01_GELATO_openning` , we will first import the grid data created in :ref:`01_lavel_koshi` session. .. _01_GELATO_import: .. figure:: images/01/GELATO_import.png :align: center :width: 100% : Grid data import Right click [Grid(No Data)] and select [Import] as (:numref:`01_GELATO_import` ). .. _01_GELATO_koshi_1: .. figure:: images/01/GELATO_koshi_1.png :align: center :width: 600pt : Select CGNS file contains grid data As shown in :numref:`01_GELATO_koshi_1`, select [Case1.cgn] which contains the grid data used in the previous section of [Computational Results of Nays2DH], and click [Open]. .. _01_GELATO_wng: .. figure:: images/01/GELATO_wng.png :align: center :width: 400pt : Warning Message A warning message is coming out as :numref:`01_GELATO_wng` , Just click [Yes] without worry, and the grid import is completed as :numref:`01_GELATO_grid` . .. _01_GELATO_grid: .. figure:: images/01/GELATO_grid.png :align: center :width: 100% : Grid import completed Single Tracer Tracking(Without Turbulent Diffusivity) -------------------------------------------------------- Condition Settings ^^^^^^^^^^^^^^^^^^^^^ Choose [Calculation Condition]->[Setting] as :numref:`01_joken_0` .. _01_joken_0: .. figure:: images/01/joken_0.png :align: center :width: 100% : Calculation Condition Settings(0) Set parameters as follows. [Flow information file name] is Locat of the CGNS file to read the calculation result of the flow field. Here, the CGNS file produced by the Nays2DH computation.( :ref:`res_Nays2DH` ). .. _01_GELATO_joken_1: .. figure:: images/01/GELATO_joken_1.png :align: center :width: 600pt : Basic Settings .. _01_GELATO_joken_2: .. figure:: images/01/GELATO_joken_2.png :align: center :width: 600pt : Primary Tracers Supplying Condition .. _01_GELATO_joken_3: .. figure:: images/01/GELATO_joken_3.png :align: center :width: 600pt : Secondary Tracers Supplying Condition .. _01_GELATO_joken_4: .. figure:: images/01/GELATO_joken_4.png :align: center :width: 600pt : Time Settings for Normal Tracers .. _01_GELATO_joken_5: .. figure:: images/01/GELATO_joken_5.png :align: center :width: 600pt : Diffusion Condition Launch GELATO ^^^^^^^^^^^^^^^ From the main menu bar, select [Simulation]->[Run], then you are asked as :numref:`01_GELATO_jikko_0`. When you click [OK] and save project, the computation starts as :numref:`01_GELATO_jikko_1`. .. _01_GELATO_jikko_0: .. figure:: images/01/warning.png :align: center :width: 400pt : warning .. _01_GELATO_jikko_1: .. figure:: images/01/GELATO_jikko_1.png :align: center :width: 100% : Launch GELATO When the computation finishes, :numref:`01_GELATO_jikko_2` appears, and click [OK] for confirmation. .. _01_GELATO_jikko_2: .. figure:: images/01/GELATO_jikko_2.png :align: center :width: 250pt : Computation finished Visualization of Computational Results ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ From the main menu, select [Calculation Result]->[Open ne 2D Post-processing Window] as :numref:`01_GELATO_kekka_0`, then [2D Post Processing Window] will appear. .. _01_GELATO_kekka_0: .. figure:: images/01/GELATO_kekka_0.png :align: center :width: 100% : Open 2D Post Processing Window | Right-click [Primary Nomal Tracers] and [Secondary Nomal Tracers] in the [Object Browser] and click [propertie]. | Then [Particles Scalar Setting] that appears,and you can set the primary and secondary have different colors by setting like :numref:`01_GELATO_kekka_2`. .. _01_GELATO_kekka_2: .. figure:: images/01/GELATO_kekka_2.png :align: center :width: 100% : Setting particles colors From the main menu, select [Animation]->[Start/Stop] as :numref:`01_GELATO_kekka_1`, animation starts ( :numref:`01_GELATO_00` ). .. _01_GELATO_kekka_1: .. figure:: images/01/GELATO_kekka_1.png :align: center :width: 100% : Visualization of computational results .. _01_GELATO_00: .. figure:: images/01/GELATO_00.gif :align: center :width: 70% : Tracer movement(No diffusivity) It is obviously very simple because it doesn't including any turbulent effect (:numref:`01_GELATO_00`). Single Tracer Tracking(With Turbulent Diffusivity) ----------------------------------------------------- Setting Computational Condition ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Change the calculation conditions to take into account for the effect of turbulent diffusion. From the main menu, select [Calculation Conditions] → [Setting], and show the :numref:`01_GELATO_joken_6`. check the box of [Diffusion Condition]->[Diffusivity Correction] , set the parameter [A Value] to [1], and then click "OK". .. _01_GELATO_joken_6: .. figure:: images/01/GELATO_joken_6.png :align: center :width: 600pt : Calculation Condition (Diffusion Condition) Launch GELATO and the Results Visualization ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Computation can be conducted through the same procedure as previous example, the animation becomes as :numref:`01_GELATO_01`. .. _01_GELATO_01: .. figure:: images/01/GELATO_01.gif :align: center :width: 70% : Tracer Movement(With Turbulent Diffusivity A=1) When the value of A is set as [10], the results become as :numref:`01_GELATO_10`, the effect of the turbulent becomes stronger. .. _01_GELATO_10: .. figure:: images/01/GELATO_10.gif :align: center :width: 70% : Tracer Movement(With Turbulent Diffusivity A=10)