Outils pour utilisateurs

Outils du site


datasetmeca:flow2dpastheadtedcylinder1

Différences

Ci-dessous, les différences entre deux révisions de la page.

Lien vers cette vue comparative

Prochaine révision
Révision précédente
datasetmeca:flow2dpastheadtedcylinder1 [2021/12/15 11:51] – créée yanndatasetmeca:flow2dpastheadtedcylinder1 [2023/05/26 18:54] (Version actuelle) yann
Ligne 1: Ligne 1:
 ===== 2D flow around a heated square-cylinder - $Re_{h}= 50$, $Ra= 5\; 10^6$ ===== ===== 2D flow around a heated square-cylinder - $Re_{h}= 50$, $Ra= 5\; 10^6$ =====
 +
 <note> <note>
-Numerical study based on the publication of Pham et al., Journal of Marine Science and Technology, Vol. 18, No. 3, pp. 361-368 (2010)\\ +A heated square-cylinder is placed in channel flow
-The motion of the cylinder is vertical and forced with sinusoïdal function\\+
 </note> </note>
  
-**Author** :  Y. Fraigneau  CNRS-LIMSI (yann.fraigneau@limsi.fr) \\ +**Author** :  Y. Fraigneau, LISN, UMR CNRS 9015 (yann.fraigneau@lisn.fr) \\ 
-**Date**   :  July 2016\\ +**Date**   :  June 2019\\ 
-**Simulation type** : DNS ([[https://sunfluidh.limsi.fr|Sunfluidh code]]) \\ +**Simulation type** : DNS (code Sunfluidh -- [[https://sunfluidh.lisn.upsaclay.fr|Sunfluidh code]]) \\ 
-**Location** : /DATABASE_2DFLOW_AROUND_OSCILLATING_CYLINDER_DNS \\ +**Location** : /DATABASE_2DFLOW_AROUND_HEATED_SQUARECYLINDER_DNS \\ 
-**Status** : Free access \\ +**Status** : Free access for laboratory members, on request for external members \\ 
-**Data size** : ~ Gb \\+**Data size** : ~ Gb \\
  
-{{ :datasetmeca:snaphot_oscyl_pham.png?900 |}}+{{ :datasetmeca:snapshot_heatcyl_snap.png?900 |}} 
  
-[[https://sunfluidh.limsi.fr/sunfluidh:sunfluidh_gallery| A video is available here ]]+[[https://sunfluidh.lisn.upsaclay.fr/doku.php?id=sunfluidh:sunfluidh_gallery#d_flow_past_a_heated_square-cylinder| A video is available here ]]
  
 [[:start|Come back to first page]] [[:start|Come back to first page]]
Ligne 23: Ligne 23:
 === 2D sketch === === 2D sketch ===
  
-{{ :datasetmeca:sketch2d_config_oscyl.jpg?800 |}} +{{ :datasetmeca:sketch_config_heatcyl.jpg?600 |}}
-\\ +
- +
  
 === Referential : cartesian geometry === === Referential : cartesian geometry ===
Ligne 32: Ligne 29:
       * x(i) : downstream direction       * x(i) : downstream direction
       * y(j) : normal direction       * y(j) : normal direction
-    - origin : +    - origin : lower left corner of the computational domain 
-      * $x_0= 0: centre of the cylinder +      * $x_0= 10h$  
-      * $y_0= 0$ : mid-height+      * $y_0= -6h$ 
  
 ===  Reference scales === ===  Reference scales ===
  
    * Density           : mass density of the fluid ($\rho_0$)    * Density           : mass density of the fluid ($\rho_0$)
-   Length            : cylinder diameter ($D$)+   Lengths           : cylinder size ($h$) and $H_u= 10h$ the distance between the domain's top and the square-cylinder's top.
    * Velocity          : velocity at inlet ($U_0$)    * Velocity          : velocity at inlet ($U_0$)
    * Dynamic viscosity : dynamic viscosity of the fluid ($\mu_0$)    * Dynamic viscosity : dynamic viscosity of the fluid ($\mu_0$)
-   * Reynolds number : $Re_H= \frac{\rho_0.U_0.D}{\mu_0}= 185+   * thermal diffusivity :  ($\kappa_0$) 
-   Strouhal number : $St_0=\frac{f}{f_0}$ (oscillating-cylinder frequency $f$ over vortex-shedding frequency $f_0$ when the cylinder is at rest) +   * Reynolds number   : $Re_h= \frac{\rho_0.U_0.h}{\mu_0}= 50
-   * motion magnitude : $A \times D(where A is the motion factor)  +   Rayleigh number   : $Ra_{Hu} \frac{\rho_0.\beta.g\Delta T.H_u^3}{\mu_0 \kappa_0}= 5\;10^6$
 === Non-dimensionalised data === === Non-dimensionalised data ===
  
    * velocity   : $U^*=\frac{U}{U_0}$     * velocity   : $U^*=\frac{U}{U_0}$ 
    * density    : $\rho^*= \rho/\rho_0= 1$    * density    : $\rho^*= \rho/\rho_0= 1$
-   * coordinates : $x*=\frac{x}{D}$, $y^*=\frac{y}{D}$+   * coordinates : $x*=\frac{x}{h}$, $y^*=\frac{y}{h}$ 
 +   * temperature : $T^*=\frac{T-T_c}{T_h-T_c}$ ($T_h$ and $T_c$ are respectively the temperature imposed at the cylinder's walls and the temperature at the inlet and channel's walls)
  
 === Computational domain === === Computational domain ===
  
    - **Domain scope**    - **Domain scope**
-           - computational domain size +           - computational domain size (channel flow) 
-                * Downstream direction(x) : $L_x^*= 90.0$ (upward part $L_u= 40$, downward part $L_d=50$) +                * Downstream direction(x) : $L_x^*= 50.0$ (upward part $L_u= -10$, downward part $L_d=40$) 
-                * Normal direction    (y) : $L_y^*= 80.0$ +                * Normal direction    (y) : $H*= 20.0$ 
-          - Oscillating cylinder (vertical oscillating motion around the mid-height of the domain) +          - heated square cylinder  
-                * centre position at rest : $x_c= 0.0$ and $y_c= 0.0$ +                * size: $h^*= 1$ 
-                * diameter: $D^*= 1$ +                * wall temperature T_h^*= 1 
-                * oscillation frequency : $St= 0.95$ or $St= 1.2$ +
-                * The cylinder is modeled with a pseudo-penalisation method (Pasquetti et al., Applied Numerical Mathematics, 2008)+
   - **Boundary conditions**   - **Boundary conditions**
       * Inlet           : imposed pressure uniform velocity ($U_0=1$)       * Inlet           : imposed pressure uniform velocity ($U_0=1$)
       * Outlet          : Orlansky type       * Outlet          : Orlansky type
-      * Bottom and top conditions symmetry+      * Bottom and top conditions walls at (T_c^*=0)
   - **Spatial resolution**   - **Spatial resolution**
-       * mesh size : $480 \times 480 $  (576.000 cells)+       * mesh size : $512 \times 256 $  (131.072 cells)
        * About cell-size           * About cell-size   
-          * $\Delta x*$ : from $0.0125$ to $0.78785$$ (downstream direction) +          * $\Delta x*$ : from $0.03125$ to $0.723$ (downstream direction) 
-          * $\Delta y*$ : from  $0.0125$ to $0.65924$ (normal direction)+          * $\Delta y*$ : from  $0.03125$ to $0.187$ (normal direction)
  
 [[:start|Come back to first page]] [[:start|Come back to first page]]
Ligne 79: Ligne 75:
  
    * **Time series **    * **Time series **
-       * from probes 
           * Physical quantities : velocity components along x and y directions (u,v) and pressure (p)           * Physical quantities : velocity components along x and y directions (u,v) and pressure (p)
           * 1 probe           * 1 probe
-          * Time step= $\,10^{-2}$ time unit +          * Time step= $\approx 6 \,10^{-2}$ time unit (irregular time step) 
-          * Time range : 400 to 600 time units +          * Time range : 250 time units 
-          * Location : Xi= 3., Xj=0.0  +          * Location : Xi= 20., Xj=2.0  
-          * File name (per physical quantity): x_ins_00000.d   with x= u,v,p +          * File name (per physical quantity): x_ins_00000.d   with x= u,v,t 
-      * about cylinder motion +
-          * center position as a function of time (ASCII file ibm_position01.dat) +
-          * cylinder velocity as a function of time (ASCII file ibm_velocity01.dat) +
-          * global force components applied to the cylinder (ASCII file  ibm_force01.dat) +
-          * force contributions exerted on the moving cylinder (ASCII file  ibm_force_contribution01.dat, __the four last columns can be ignored__)+
  
    * **3D snapshots**    * **3D snapshots**
          * Instantaneous fields : velocity components in x, y and z directions (U,V), the pressure (P) and the phase function related to the body motion (TRACE)          * Instantaneous fields : velocity components in x, y and z directions (U,V), the pressure (P) and the phase function related to the body motion (TRACE)
-         * Recording rate : 0.4 time unit +         * Recording rate : time unit 
-         * Time range from from 200.0 to 600.0 time units+         * Time range from from 20.0 to 250.0 time units
          * File name : res_xxxxx_yyyyyyy.d           * File name : res_xxxxx_yyyyyyy.d 
-               * MPI subdomain ID: 0 +               * MPI subdomain ID: from to 3 (case with MPI domain decomposition) 
-               * Time ID : from to 502+               * Time ID : from to 250
  
    
Ligne 106: Ligne 97:
 ---- ----
  
-**Data size** : ~ Gb \\ +**Data size** : ~ Gb \\ 
-**Main directory** : /vol/DATABASE_MECA/DATABASE_2DFLOW_AROUND_OSCILLATING_CYLINDER_DNS \\ +**Main directory** : /vol/DATABASE_MECA/DATABASE_2DFLOW_HEATED_SQUARRECYLINDER_DNS \\
-For more details about files, see the [[https://sunfluidh.limsi.fr| wiki doc of Sunfluidh]]+
  
-Intermediate directories (four cases of study)  +Directories & files
-    * $St= 0.95$ and $A= 0.2$ - directory : CAS_F0.95_A0.2 +
-    * $St= 0.95$ and $A= 0.5$ - directory : CAS_F0.95_A0.5 +
-    * $St= 1.20$ and $A= 0.2$ - directory : CAS_F1.2_A0.2 +
-    * $St= 1.20$ and $A= 0.5$ - directory : CAS_F1.2_A0.5 +
- +
-Endpoint directories & files+
  
                                                                    
     /DATASETUP       : ASCII files     /DATASETUP       : ASCII files
-                       input data file for sunfluidh : input3d.dat+                       input data file for sunfluidh : input3d.dat (and input3d.dat_mpi4 for the case with MPI domain decomposition)
          
     /GRID            : ASCII files     /GRID            : ASCII files
                        input data file         : data_meshgen.d                        input data file         : data_meshgen.d
-                       grid files for sunfluidh: maillx.d, mailly.d, maillz.d+                       grid files for sunfluidh: maillx.d, mailly.d, maillz.d (sequential case) 
 +                                                 mesh.tar (archive containing grid files for a domain decomposition with 4 MPI processes)
            
     /SNAPSHOTS       : snapshots binary files res_xxxxx_yyyyyyy.d     /SNAPSHOTS       : snapshots binary files res_xxxxx_yyyyyyy.d
 +    /SNAPSHOTS_MPI4  : snapshots resulting from a domain decomposition with 4 MPI processes
          
     /TIMESERIES      : ASCII files     /TIMESERIES      : ASCII files
-                       u_ins_00000.d , v_ins_00000.d and p_ins_00000.d (timeseries of the velocity components and pressure from probe) +                       u_ins_00000.d , v_ins_00000.d and t_ins_00000.d (timeseries of the velocity components and temperature from probe) 
-                       ibm_position01.dat, ibm_velocity01.dat (about cylinder motion) +
-                       ibm_force01.dat, ibm_force_contribution01.dat (about forces exerted on the cylinder) +
-                       +
                                    
 [[:start|Come back to first page]] [[:start|Come back to first page]]
datasetmeca/flow2dpastheadtedcylinder1.1639565465.txt.gz · Dernière modification : 2021/12/15 11:51 de yann

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki