Skip to content

Commit

Permalink
Deploying to gh-pages from @ ebddc1f 🚀
Browse files Browse the repository at this point in the history
  • Loading branch information
CyrilleBonamy committed Apr 16, 2024
1 parent c0d9423 commit 81a4559
Show file tree
Hide file tree
Showing 5 changed files with 9 additions and 7 deletions.
Binary file added bed_profiles_kOmega.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file added maximum_depth_kOmega.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
5 changes: 3 additions & 2 deletions tutorials_LES.html
Original file line number Diff line number Diff line change
Expand Up @@ -284,7 +284,7 @@ <h3>Contents</h3>
maxAlphaCo 1;

favreAveraging false;
</pre></section><section id="lauch_Chan"><h3><a href="#lauch_Chan">Computation launching</a></h3><p>As for the previous case, you can launch the computation by executing <em>Allrun</em> for turbulence initialization. The simulations is run in parallel on 16 cores. Once the end time has been modified and favreAveraging keyword set to true, you can launch the computation of time average variables by executing <em>AllrunAverage</em>.</p></section><section id="post_Chan"><h3><a href="#post_Chan">Post-processing using python</a></h3><p>The post-processing python scripts writedata_tuto3DChannel560.py and plot_tuto3DChannel560.py are located in the folder tutorials/Py. To run this script, the latest output (16s) should be reconstructed.</p><p>The script writedata_tuto3DChannel560.py reads time averaged OpenFoam data, perform a spatial averaging operation and stores the 1D vertical profiles in a netCDF file located in the postProcessing directory of the case.</p><p>The script plot_tuto3DChannel560.py reads the netCDF file and plots vertical profile of concentration, velocities and Reynolds stresses</p><figure class="m-figure"><img src="tuto3DChannel560.png" alt="Image" style="width: 700px;" /><figcaption>Figure 2: Concentration, velocity and Reynolds stress profiles for the 3DChannel560 tutorial.</figcaption></figure></section></section><section id="OscillSheetFlow_sec"><h2><a href="#OscillSheetFlow_sec">3DOscillSheetFlow: Oscillatory sheet flow</a></h2><p>In this tutorial, the oscillatory sheet flow experimental configuration with fine sand from O&#x27;Donoghue and Wright (2004) is reproduced using Large-Eddy Simulation (LES).</p><section id="preprocOdo"><h3><a href="#preprocOdo">Pre-processing</a></h3><p>This tutorial is distributed with SedFoam under the folder <em>sedFoamDirectory/tutorials/3DOscillSheetFlow</em>.</p></section><section id="meshgen_Odo"><h3><a href="#meshgen_Odo">Mesh and boundary conditions</a></h3><p>The numerical domain is a periodic box presented in figure 3 with <svg class="m-math" style="width: 6.744em; height: 1.494em; vertical-align: -0.380em;" viewBox="0 -10.711001 64.74148 14.346348">
</pre></section><section id="lauch_Chan"><h3><a href="#lauch_Chan">Computation launching</a></h3><p>As for the previous case, you can launch the computation by executing <em>Allrun</em> for turbulence initialization. The simulations is run in parallel on 16 cores. Once the end time has been modified and favreAveraging keyword set to true, you can launch the computation of time average variables by executing <em>AllrunAverage</em>.</p></section><section id="post_Chan"><h3><a href="#post_Chan">Post-processing using python</a></h3><p>The post-processing python scripts <em>writedata_tuto3DChannel560.py</em> and <em>plot_tuto3DChannel560.py</em> are located in the folder <em>tutorials/Py</em>. To run this script, the latest output (16s) should be reconstructed.</p><p>The script <em>writedata_tuto3DChannel560.py</em> reads time averaged OpenFoam data, perform a spatial averaging operation and stores the 1D vertical profiles in a netCDF file located in the postProcessing directory of the case.</p><p>The script <em>plot_tuto3DChannel560.py</em> reads the netCDF file and plots vertical profile of concentration, velocities and Reynolds stresses</p><figure class="m-figure"><img src="tuto3DChannel560.png" alt="Image" style="width: 700px;" /><figcaption>Figure 2: Concentration, velocity and Reynolds stress profiles for the 3DChannel560 tutorial.</figcaption></figure></section></section><section id="OscillSheetFlow_sec"><h2><a href="#OscillSheetFlow_sec">3DOscillSheetFlow: Oscillatory sheet flow</a></h2><p>In this tutorial, the oscillatory sheet flow experimental configuration with fine sand from O&#x27;Donoghue and Wright (2004) is reproduced using Large-Eddy Simulation (LES).</p><section id="preprocOdo"><h3><a href="#preprocOdo">Pre-processing</a></h3><p>This tutorial is distributed with SedFoam under the folder <em>sedFoamDirectory/tutorials/3DOscillSheetFlow</em>.</p></section><section id="meshgen_Odo"><h3><a href="#meshgen_Odo">Mesh and boundary conditions</a></h3><p>The numerical domain is a periodic box presented in figure 3 with <svg class="m-math" style="width: 6.744em; height: 1.494em; vertical-align: -0.380em;" viewBox="0 -10.711001 64.74148 14.346348">
<title>
$\delta=\sqrt{2\nu^f/\omega}$
</title>
Expand Down Expand Up @@ -415,7 +415,8 @@ <h3>Contents</h3>
maxAlphaCo 0.3;

maxDeltaT 1e-3;
</pre></section><section id="lauch_Odo"><h3><a href="#lauch_Odo">Computation launching</a></h3><p>As for the previous case, you can launch the computation by executing <em>Allrun</em>. The simulations is run in parallel on 16 cores.</p></section><section id="post_Chan"><h3><a href="#post_Chan">Post-processing using python</a></h3><p>The post-processing python scripts writedata_tuto3DOscillSheetFlow.py and plot_tuto3DOscillSheetFlow.py are located in the folder tutorials/Py. To run this script, all the simulation outputs after 20s (and only them) should be reconstructed using the following command:</p><p>reconstructPar -time &#x27;20:&#x27;*</p><p>The script writedata_tuto3DOscillSheetFlow.py reads OpenFoam data, performs a spatial averaging operation, then performs a phase averaging operation taking into account four wave periods and stores the 1D vertical profiles in netCDF files located in the postProcessing directory of the case.</p><p>The script plot_tuto3DOscillSheetFlow.py reads the netCDF files and plots vertical profile of concentration, velocities and Reynolds stresses at different moments of the wave period.</p><figure class="m-figure"><img src="tuto3DOscillSheetFlow.png" alt="Image" style="width: 600px;" /><figcaption>Figure 4: Concentration and velocity profiles for the 3DOscillSheetFlow tutorial.</figcaption></figure></section></section>
</pre></section><section id="lauch_Odo"><h3><a href="#lauch_Odo">Computation launching</a></h3><p>As for the previous case, you can launch the computation by executing <em>Allrun</em>. The simulations is run in parallel on 16 cores.</p></section><section id="post_Chan"><h3><a href="#post_Chan">Post-processing using python</a></h3><p>The post-processing python scripts <em>writedata_tuto3DOscillSheetFlow.py</em> and <em>plot_tuto3DOscillSheetFlow.py</em> are located in the folder tutorials/Py. To run this script, all the simulation outputs after 20s (and only them) should be reconstructed using the following command:</p><pre>reconstructPar -time &#x27;20:&#x27;
</pre><p>The script <em>writedata_tuto3DOscillSheetFlow.py</em> reads OpenFoam data, performs a spatial averaging operation, then performs a phase averaging operation taking into account four wave periods and stores the 1D vertical profiles in netCDF files located in the postProcessing directory of the case.</p><p>The script <em>plot_tuto3DOscillSheetFlow.py</em> reads the netCDF files and plots vertical profile of concentration, velocities and Reynolds stresses at different moments of the wave period.</p><figure class="m-figure"><img src="tuto3DOscillSheetFlow.png" alt="Image" style="width: 600px;" /><figcaption>Figure 4: Concentration and velocity profiles for the 3DOscillSheetFlow tutorial.</figcaption></figure></section></section>
</div>
</div>
</div>
Expand Down
2 changes: 1 addition & 1 deletion tutorials_RAS.html
Original file line number Diff line number Diff line change
Expand Up @@ -1612,7 +1612,7 @@ <h3>Contents</h3>
maxAlphaCo 0.3;

maxDeltaT 15.e-5;
</pre></section><section id="lauch_SP"><h3><a href="#lauch_SP">Computation launching</a></h3><p>As for the previous case, you can launch the computation by executing * Allrun *. The simulations is run in parallel on 16 cores.</p></section><section id="post_SP"><h3><a href="#post_SP">Post-processing using python</a></h3><p>The post-processing python scripts plot_profiles2DPipelineScour.py and plot_depth2DPipelineScour.py are located in the folder tutorials/Py.</p><p>The script _profiles2DPipelineScour.py plots bed profiles at 11s, 18s and 25s compared with experimental data like in the following figure:</p><figure class="m-figure"><img src="bed_profiles_kEpsilon.png" alt="Image" style="width: 600px;" /><figcaption>Figure 8: Sediment bed profiles at 11s (top), 18s (middle) and 25s (bottom) using a \f$ k-\omega \f$ turbulence model.</figcaption></figure><p>The script plot_depth2DPipelineScour.py plots the maximum scour depth compared with experimental data like in the following figure:</p><figure class="m-figure"><img src="maximum_depth_kEpsilon.png" alt="Image" style="width: 600px;" /><figcaption>Figure 9: Maximum scour depth using a k-omega turbulence model.</figcaption></figure></section></section><section id="Scour3DCylinder"><h2><a href="#Scour3DCylinder">3DScourCylinder: Scour around a vertical cylinder</a></h2><p>In this tutorial, the experimental configuration from Roulund et al (2005) is reproduced numerically.</p><section id="preprocS3D"><h3><a href="#preprocS3D">Pre-processing</a></h3><p>This tutorial is distributed with SedFoam under the folder <em>sedFoamDirectory/tutorials/RAS/3DScour</em>.</p><p>A 1D boundary layer simulation needs to be run (see folder 1D and Allrun.pre) for fore information. The result of this simulations is then used to intialize the 3D configuation and to impose the inlet boundary condition.</p></section><section id="meshgen_S3D"><h3><a href="#meshgen_S3D">Mesh generation</a></h3><p>The numerical domain dimensions are presented in figure 10. A cylinder of diameter <svg class="m-math" style="width: 5.497em; height: 0.851em; vertical-align: -0.000em;" viewBox="0 -8.169366 52.774021 8.169366">
</pre></section><section id="lauch_SP"><h3><a href="#lauch_SP">Computation launching</a></h3><p>As for the previous case, you can launch the computation by executing * Allrun *. The simulations is run in parallel on 16 cores.</p></section><section id="post_SP"><h3><a href="#post_SP">Post-processing using python</a></h3><p>The post-processing python scripts plot_profiles2DPipelineScour.py and plot_depth2DPipelineScour.py are located in the folder tutorials/Py.</p><p>The script _profiles2DPipelineScour.py plots bed profiles at 11s, 18s and 25s compared with experimental data like in the following figure:</p><figure class="m-figure"><img src="bed_profiles_kOmega.png" alt="Image" style="width: 600px;" /><figcaption>Figure 8: Sediment bed profiles at 11s (top), 18s (middle) and 25s (bottom) using a \f$ k-\omega \f$ turbulence model.</figcaption></figure><p>The script plot_depth2DPipelineScour.py plots the maximum scour depth compared with experimental data like in the following figure:</p><figure class="m-figure"><img src="maximum_depth_kOmega.png" alt="Image" style="width: 600px;" /><figcaption>Figure 9: Maximum scour depth using a k-omega turbulence model.</figcaption></figure></section></section><section id="Scour3DCylinder"><h2><a href="#Scour3DCylinder">3DScourCylinder: Scour around a vertical cylinder</a></h2><p>In this tutorial, the experimental configuration from Roulund et al (2005) is reproduced numerically.</p><section id="preprocS3D"><h3><a href="#preprocS3D">Pre-processing</a></h3><p>This tutorial is distributed with SedFoam under the folder <em>sedFoamDirectory/tutorials/RAS/3DScour</em>.</p><p>A 1D boundary layer simulation needs to be run (see folder 1D and Allrun.pre) for fore information. The result of this simulations is then used to intialize the 3D configuation and to impose the inlet boundary condition.</p></section><section id="meshgen_S3D"><h3><a href="#meshgen_S3D">Mesh generation</a></h3><p>The numerical domain dimensions are presented in figure 10. A cylinder of diameter <svg class="m-math" style="width: 5.497em; height: 0.851em; vertical-align: -0.000em;" viewBox="0 -8.169366 52.774021 8.169366">
<title>
$ D=10cm $
</title>
Expand Down
Loading

0 comments on commit 81a4559

Please sign in to comment.