Documentation du code de simulation numérique SUNFLUIDH

Cette page est en lecture seule. Vous pouvez afficher le texte source, mais ne pourrez pas le modifier. Contactez votre administrateur si vous pensez qu'il s'agit d'une erreur.
===== Examples of wall boundary conditions  =====

The user finds here some examples illustrating the construction of the wall boundary conditions :
   * Heated enclosed cavity : Temperature imposed on sides walls & adiabatic conditions on upper and lower walls
   * Two bars in an enclosed cavity : One bar with imposed temperature, the second one with imposed heat flux 
   * An enclosed cavity with a thermally conductive bar
   * A heated bar in a lid-driven cavity

<WRAP info>
Only 2D geometrical configurations are considered. So, data associated to the 3rd direction K could be removed. They are still kept in these examples even if they are useless.
 
</WRAP>

Some abbreviations considered here :
   * **__WBCT__** (Wall Boundary Condition Type) refers any namelist describing the type of wall boundary conditions applied to different physical quantities (temperature,velocity, species) :
     * The enthalpy or temperature : [[sunfluidh:heat_wall_boundary_condition_setup_namelist|Namelist "Heat_Wall_Boundary_Condition_Setup"]]
     * The velocity components : [[sunfluidh:velocity_wall_boundary_condition_setup_namelist|Namelist "Velocity_Wall_Boundary_Condition_Setup"]]
     * The species mass fractions : [[sunfluidh:species_wall_boundary_condition_setup_namelist|Namelist "Species_Wall_Boundary_Condition_Setup"]]

   * **__IBD__** (Immersed Body Descriptor) defines the geometry and the location of the immersed bodies

-----
==== Heated enclosed cavity ====
-----
<WRAP info>
   * The domain is a square cavity (2D geometry).
   * Bottom and top walls are insulated. The heat flux is null.
   * The left wall is set to the upper constant and uniform temperature : T= 305K.\\ The right wall is set to the lower constant and uniform "temperature : T= 295K.
   * The velocity boundary conditions are no-slip and impermeability conditions.
   * There is no immersed body. 
   * Walls of the cavity are located at the ends of the domain. This is the basic domain configuration initialized by the code if any other boundary condition is not placed at the ends (like inlet, outlet, symmetrical or periodical boundary conditions)
</WRAP>
{{ :sunfluidh:example_heat_driven_cavity.jpg?direct&600 |}}


We need just to set one **WPCP**.


   &Heat_Wall_Boundary_Condition_Setup
      Wall_BC_DataSetName     ="Set1",
      West_Heat_BC_Option    = 0 ,  East_Heat_BC_Option    = 0 ,  Back_Heat_BC_Option    = 1  ,  Front_Heat_BC_Option    = 1 ,
      South_Heat_BC_Option    = 0 , North_Heat_BC_Option    = 0,
      West_Heat_Function_Type= 0 ,  East_Heat_Function_Type= 0 ,  Back_Heat_Function_Type= 0  ,  Front_Heat_Function_Type= 0 ,   
      South_Heat_Function_Type= 0 , North_Heat_Function_Type= 0,
      West_Wall_BC_Value= 305.   ,  East_Wall_BC_Value= 295.  ,  Back_Wall_BC_Value= 0.0   ,  Front_Wall_BC_Value= 0.0   ,   
      South_Wall_BC_Value= 0.0   , North_Wall_BC_Value= 0.0  /
      
      
      &Velocity_Wall_Boundary_Condition_Setup
      Wall_BC_DataSetName     ="Set1",
      West_Velocity_I_BC_Option= 0 ,  East_Velocity_I_BC_Option= 0 ,  Back_Velocity_I_BC_Option= 0 , Front_Velocity_I_BC_Option= 0, 
      South_Velocity_I_BC_Option= 0 , North_Velocity_I_BC_Option= 0 ,
      West_Velocity_J_BC_Option= 0 ,  East_Velocity_J_BC_Option= 0 ,  Back_Velocity_J_BC_Option= 0 , Front_Velocity_J_BC_Option= 0, 
      South_Velocity_J_BC_Option= 0 , North_Velocity_J_BC_Option= 0 ,
      West_Velocity_K_BC_Option= 0 ,  East_Velocity_K_BC_Option= 0 ,  Back_Velocity_K_BC_Option= 0 , Front_Velocity_K_BC_Option= 0, 
      South_Velocity_K_BC_Option= 0 , North_Velocity_K_BC_Option= 0 ,
      West_Velocity_I_Function_Type= 0,  East_Velocity_I_Function_Type= 0, Back_Velocity_I_Function_Type= 0,  
      Front_Velocity_I_Function_Type= 0, South_Velocity_I_Function_Type= 0, North_Velocity_I_Function_Type= 0 ,
      West_Velocity_J_Function_Type= 0,  East_Velocity_J_Function_Type= 0, Back_Velocity_J_Function_Type= 0, 
      Front_Velocity_J_Function_Type= 0, South_Velocity_J_Function_Type= 0, North_Velocity_J_Function_Type= 0 ,
      West_Velocity_K_Function_Type= 0,  East_Velocity_K_Function_Type= 0, Back_Velocity_K_Function_Type= 0, 
      Front_Velocity_K_Function_Type= 0, South_Velocity_K_Function_Type= 0, North_Velocity_K_Function_Type= 0 ,
      West_Wall_Velocity_I= 0.0   ,   East_Wall_Velocity_I= 0.0   ,   Back_Wall_Velocity_I= 0.0    , Front_Wall_Velocity_I= 0.0   ,
      South_Wall_Velocity_I= 0.0    , North_Wall_Velocity_I= 0.0 ,
      West_Wall_Velocity_J= 0.0   ,   East_Wall_Velocity_J= 0.0   ,   Back_Wall_Velocity_J= 0.0    , Front_Wall_Velocity_J= 0.0   ,
      South_Wall_Velocity_J= 0.0    , North_Wall_Velocity_J= 0.0 ,
      West_Wall_Velocity_K= 0.0   ,   East_Wall_Velocity_K= 0.0   ,   Back_Wall_Velocity_K= 0.0    , Front_Wall_Velocity_K= 0.0  ,
      South_Wall_Velocity_K= 0.0    , North_Wall_Velocity_K= 0.0 /

      
<note important>
   * There is just one WBCT. This name is "Set1". It is automatically linked with the walls belonging to the domain ends.
   * The  no-slip and impermeable boundary conditions are actually used by default in the code. For this reason, this kind of boundary conditions could be not explicitly set.
</note>
==== Two bars in an enclosed cavity ====
<WRAP info>
   * The domain is a square cavity with a length of $L= 1.0$ (2D geometry).
   * Walls of the cavity are set to the constant and uniform temperature : $T= 300 K$.
   * Walls of the first bar are set to the constant and uniform heat flux : $Q= 5W/m^2$. Keep in mind that :
     * heating the fluid from walls is defined with a negative heat flux
     * cooling the fluid from walls is defined with a positive heat flux
   * Walls of the second bar are set to the same boundary conditions as the cavity : $T=300 K$.
   * The velocity boundary conditions are no-slip and impermeability conditions for all the walls (cavity and bars). As they correspond to the default boundary conditions, they are not explicitly set here.
   * There are two immersed bodies : the couple of bars. Each bar is a rectangle with a size of $0.1$x$0.2$. The center of the first bar is $0.3$x$0.3$, the second one is $0.7$x$0.7$
</WRAP>
{{ :sunfluidh:example_enclosed_cavity_with_2_bars.jpg?600direct |}}
We need two **WBCT** for setting the two different sets of wall boundary conditions, the first one for the imposed temperature, the second one for the imposed heat flux. We also need to set two ** IBD** in order to define the positions of the couple of bars.

__The first **WBCT**__

    &Heat_Wall_Boundary_Condition_Setup
     Wall_BC_DataSetName     ="Set1",
      West_Heat_BC_Option    = 0 ,  East_Heat_BC_Option    = 0 ,  Back_Heat_BC_Option    = 0  ,  Front_Heat_BC_Option    = 0 , South_Heat_BC_Option    = 0 , North_Heat_BC_Option    = 0,
      West_Heat_Function_Type= 0 ,  East_Heat_Function_Type= 0 ,  Back_Heat_Function_Type= 0  ,  Front_Heat_Function_Type= 0 , South_Heat_Function_Type= 0 , North_Heat_Function_Type= 0,
      West_Wall_BC_Value= 300.   ,  East_Wall_BC_Value= 300.  ,  Back_Wall_BC_Value= 300.0   ,  Front_Wall_BC_Value= 300.0   , South_Wall_BC_Value= 0.0   , North_Wall_BC_Value= 0.0  /

__The second **WBCT**__

    &Heat_Wall_Boundary_Condition_Setup
     Wall_BC_DataSetName     ="Set2",
      West_Heat_BC_Option    = 1 ,  East_Heat_BC_Option    = 1 ,  Back_Heat_BC_Option    = 1  ,  Front_Heat_BC_Option    = 1 , South_Heat_BC_Option    = 0 , North_Heat_BC_Option    = 0,
      West_Heat_Function_Type= 0 ,  East_Heat_Function_Type= 0 ,  Back_Heat_Function_Type= 0  ,  Front_Heat_Function_Type= 0 , South_Heat_Function_Type= 0 , North_Heat_Function_Type= 0,
      West_Wall_BC_Value= -5.0   ,  East_Wall_BC_Value= -5.0  ,  Back_Wall_BC_Value= -5.0  ,  Front_Wall_BC_Value= -5.0  , South_Wall_BC_Value= 0.0   , North_Wall_BC_Value= 0.0  /
      
__The first **IBD** (the first bar)__

   &Polyhedral_Immersed_Bodies
       Xi_1= 0.25 , Xj_1= 0.20 , Xk_1= 0.0 ,
       Xi_2= 0.35 , Xj_2= 0.20 , Xk_2= 0.0 ,
       Xi_3= 0.35 , Xj_3= 0.40 , Xk_3= 0.0 ,
       Xi_4= 0.25 , Xj_4= 0.40 , Xk_4= 0.0 ,
       Wall_BC_DataSetName     ="Set2"/
       
__The second **IBD** (the second bar)__

   &Polyhedral_Immersed_Bodies
       Xi_1= 0.65 , Xj_1= 0.60 , Xk_1= 0.0 ,
       Xi_2= 0.75 , Xj_2= 0.60 , Xk_2= 0.0 ,
       Xi_3= 0.75 , Xj_3= 0.80 , Xk_3= 0.0 ,
       Xi_4= 0.65 , Xj_4= 0.80 , Xk_4= 0.0 ,
       Wall_BC_DataSetName     ="Set1"/
       
<note important>
   * The first immersed body is linked to the second **WBCT**(imposed heat flux) and the second body is linked to the first **WBCT** (imposed temperature) by means of the variable "Wall_BC_DataSetName".
   * The **WBCT** named "Set1" is also automatically linked to the walls of the domain end's (cavity walls).
   * The rank order of the boundary conditions is unimportant.
   * **About imposed heat flux**, keep in mind that :
     * heating the fluid from walls is defined with a negative heat flux
     * cooling the fluid from walls is defined with a positive heat flux
</note>

==== Enclosed cavity with a thermally conductive bar ====
<WRAP info>
   * The domain is a square cavity with a length of $L= 1.0$ (2D geometry).
   * Bottom and top walls of the cavity are insulated. The heat flux is null.
   * The left and right walls of the cavity are set to 305K and 295K, respectively.
   * The bar is thermally conductive. Its size is $0.2$x$0.2$ and it is centered on the cavity.
   * The velocity boundary conditions are no-slip and impermeable conditions for all walls.
</WRAP>
We need to set two **WBCT** to define the two different sets of wall boundary conditions, the first one for the imposed temperature and heat flux at the walls of the cavity, the second one to define the physical properties of the bar. We also need to set one ** IBD** in order to define the position of the bar.
{{ :sunfluidh:example_cavity_with_bar.jpg?direct&600 |}}

__The first **WBCT**__

    &Heat_Wall_Boundary_Condition_Setup
      Wall_BC_DataSetName     ="Set1",
      West_Heat_BC_Option    = 0 ,  East_Heat_BC_Option    = 0 ,  Back_Heat_BC_Option    = 1  ,  Front_Heat_BC_Option    = 1 , South_Heat_BC_Option    = 0 , North_Heat_BC_Option    = 0,
      West_Heat_Function_Type= 0 ,  East_Heat_Function_Type= 0 ,  Back_Heat_Function_Type= 0  ,  Front_Heat_Function_Type= 0 , South_Heat_Function_Type= 0 , North_Heat_Function_Type= 0,
      West_Wall_BC_Value= 305.   ,  East_Wall_BC_Value= 295.  ,  Back_Wall_BC_Value= 0.0   ,  Front_Wall_BC_Value= 0.0   , South_Wall_BC_Value= 0.0   , North_Wall_BC_Value= 0.0  /
 
__The second **WBCT**__
     
    &Heat_Wall_Boundary_Condition_Setup
      Wall_BC_DataSetName     ="Set2",
      West_Heat_BC_Option    = 2 ,  East_Heat_BC_Option    = 2 ,  Back_Heat_BC_Option    = 2  ,  Front_Heat_BC_Option    = 2 ,
      South_Heat_BC_Option    = 0 , North_Heat_BC_Option    = 0,
      Material_Thermal_Conductivity= 1.00 , Material_Mass_Heat_Capacity= 3000.00 , Material_Density= 900.00 /
   
      
__The first **IBD** (the bar) __

   &Polyhedral_Immersed_Bodies
       Xi_1= 0.40 , Xj_1= 0.40 , Xk_1= 0.0 ,
       Xi_2= 0.40 , Xj_2= 0.40 , Xk_2= 0.0 ,
       Xi_3= 0.60 , Xj_3= 0.60 , Xk_3= 0.0 ,
       Xi_4= 0.60 , Xj_4= 0.60 , Xk_4= 0.0 ,
       Wall_BC_DataSetName     ="Set2"/


<note important>
   * The **WBCT** named "Set1" is automatically linked to the walls of the domain end's (cavity walls).
   * The bar (immersed body) is linked to the second **WBCT**.
</note>  
     
==== Heated bar in a lid-driven cavity ====

The example is very similar to the previous one but the fluid is driven by the top wall of the cavity.
<WRAP info>
   * The domain is a square cavity with a length of $L= 1.0$ (2D geometry).
   * Bottom and top walls of the cavity are insulated. The heat flux is null.
   * The left and right walls of the cavity are set to 305K and 295K, respectively.
   * The bar is thermally conductive. Its size is $0.2$x$0.2$ and it is centered on the cavity.
   * The velocity boundary conditions are no-slip and impermeability conditions for all walls **except at the top of the cavity** where $V= 1.5m.s^{-1}$. For these reason, the wall boundary conditions for the velocity must be set.
</WRAP>
We need to set two **WBCT** to define the two different sets of wall boundary conditions, the first one for the imposed temperature and heat flux at the walls of the cavity, the second one to define the physical properties of the bar. We also need to set one ** IBD** in order to define the position of the bar.

{{ :sunfluidh:example_liddrivenheatedcavity.jpg?direct&600 |}}

__The first **WBCT**__

    &Heat_Wall_Boundary_Condition_Setup
      West_Heat_BC_Option    = 0 ,  East_Heat_BC_Option    = 0 ,  Back_Heat_BC_Option    = 1  ,  Front_Heat_BC_Option    = 1 , South_Heat_BC_Option    = 0 , North_Heat_BC_Option    = 0,
      West_Heat_Function_Type= 0 ,  East_Heat_Function_Type= 0 ,  Back_Heat_Function_Type= 0  ,  Front_Heat_Function_Type= 0 , South_Heat_Function_Type= 0 , North_Heat_Function_Type= 0,
      West_Wall_BC_Value= 305.   ,  East_Wall_BC_Value= 295.  ,  Back_Wall_BC_Value= 0.0   ,  Front_Wall_BC_Value= 0.0   , South_Wall_BC_Value= 0.0   , North_Wall_BC_Value= 0.0  /
      
     &Velocity_Wall_Boundary_Condition_Setup
      West_Velocity_I_BC_Option= 0 ,  East_Velocity_I_BC_Option= 0 ,  Back_Velocity_I_BC_Option= 0 , Front_Velocity_I_BC_Option= 0, 
      South_Velocity_I_BC_Option= 0 , North_Velocity_I_BC_Option= 0 ,
      West_Velocity_J_BC_Option= 0 ,  East_Velocity_J_BC_Option= 0 ,  Back_Velocity_J_BC_Option= 0 , Front_Velocity_J_BC_Option= 0, 
      South_Velocity_J_BC_Option= 0 , North_Velocity_J_BC_Option= 0 ,
      West_Velocity_K_BC_Option= 0 ,  East_Velocity_K_BC_Option= 0 ,  Back_Velocity_K_BC_Option= 0 , Front_Velocity_K_BC_Option= 0, 
      South_Velocity_K_BC_Option= 0 , North_Velocity_K_BC_Option= 0 ,
      West_Velocity_I_Function_Type= 0,  East_Velocity_I_Function_Type= 0, Back_Velocity_I_Function_Type= 0,  
      Front_Velocity_I_Function_Type= 0, South_Velocity_I_Function_Type= 0, North_Velocity_I_Function_Type= 0 ,
      West_Velocity_J_Function_Type= 0,  East_Velocity_J_Function_Type= 0, Back_Velocity_J_Function_Type= 0, 
      Front_Velocity_J_Function_Type= 0, South_Velocity_J_Function_Type= 0, North_Velocity_J_Function_Type= 0 ,
      West_Velocity_K_Function_Type= 0,  East_Velocity_K_Function_Type= 0, Back_Velocity_K_Function_Type= 0, 
      Front_Velocity_K_Function_Type= 0, South_Velocity_K_Function_Type= 0, North_Velocity_K_Function_Type= 0 ,
      West_Wall_Velocity_I= 0.0   ,   East_Wall_Velocity_I= 0.0   ,   Back_Wall_Velocity_I= 0.0    , Front_Wall_Velocity_I= 1.5   ,
      South_Wall_Velocity_I= 0.0    , North_Wall_Velocity_I= 0.0 ,
      West_Wall_Velocity_J= 0.0   ,   East_Wall_Velocity_J= 0.0   ,   Back_Wall_Velocity_J= 0.0    , Front_Wall_Velocity_J= 0.0   ,
      South_Wall_Velocity_J= 0.0    , North_Wall_Velocity_J= 0.0 ,
      West_Wall_Velocity_K= 0.0   ,   East_Wall_Velocity_K= 0.0   ,   Back_Wall_Velocity_K= 0.0    , Front_Wall_Velocity_K= 0.0  ,
      South_Wall_Velocity_K= 0.0    , North_Wall_Velocity_K= 0.0 /
 
__The second **WBCT**__
     
    &Heat_Wall_Boundary_Condition_Setup
    West_Heat_BC_Option    = 2 ,  East_Heat_BC_Option    = 2 ,  Back_Heat_BC_Option    = 2  ,  Front_Heat_BC_Option    = 2 , 
    South_Heat_BC_Option    = 0 , North_Heat_BC_Option    = 0,
      Material_Thermal_Conductivity= 1.00 , Material_Mass_Heat_Capacity= 3000.00 , Material_Density= 900.00 /
      
      &Velocity_Wall_Boundary_Condition_Setup
      West_Velocity_I_BC_Option= 0 ,  East_Velocity_I_BC_Option= 0 ,  Back_Velocity_I_BC_Option= 0 , Front_Velocity_I_BC_Option= 0, 
      South_Velocity_I_BC_Option= 0 , North_Velocity_I_BC_Option= 0 ,
      West_Velocity_J_BC_Option= 0 ,  East_Velocity_J_BC_Option= 0 ,  Back_Velocity_J_BC_Option= 0 , Front_Velocity_J_BC_Option= 0, 
      South_Velocity_J_BC_Option= 0 , North_Velocity_J_BC_Option= 0 ,
      West_Velocity_K_BC_Option= 0 ,  East_Velocity_K_BC_Option= 0 ,  Back_Velocity_K_BC_Option= 0 , Front_Velocity_K_BC_Option= 0, 
      South_Velocity_K_BC_Option= 0 , North_Velocity_K_BC_Option= 0 ,
      West_Velocity_I_Function_Type= 0,  East_Velocity_I_Function_Type= 0, Back_Velocity_I_Function_Type= 0,  
      Front_Velocity_I_Function_Type= 0, South_Velocity_I_Function_Type= 0, North_Velocity_I_Function_Type= 0 ,
      West_Velocity_J_Function_Type= 0,  East_Velocity_J_Function_Type= 0, Back_Velocity_J_Function_Type= 0, 
      Front_Velocity_J_Function_Type= 0, South_Velocity_J_Function_Type= 0, North_Velocity_J_Function_Type= 0 ,
      West_Velocity_K_Function_Type= 0,  East_Velocity_K_Function_Type= 0, Back_Velocity_K_Function_Type= 0, 
      Front_Velocity_K_Function_Type= 0, South_Velocity_K_Function_Type= 0, North_Velocity_K_Function_Type= 0 ,
      West_Wall_Velocity_I= 0.0   ,   East_Wall_Velocity_I= 0.0   ,   Back_Wall_Velocity_I= 0.0    , Front_Wall_Velocity_I= 0.0   ,
      South_Wall_Velocity_I= 0.0    , North_Wall_Velocity_I= 0.0 ,
      West_Wall_Velocity_J= 0.0   ,   East_Wall_Velocity_J= 0.0   ,   Back_Wall_Velocity_J= 0.0    , Front_Wall_Velocity_J= 0.0   ,
      South_Wall_Velocity_J= 0.0    , North_Wall_Velocity_J= 0.0 ,
      West_Wall_Velocity_K= 0.0   ,   East_Wall_Velocity_K= 0.0   ,   Back_Wall_Velocity_K= 0.0    , Front_Wall_Velocity_K= 0.0  ,
      South_Wall_Velocity_K= 0.0    , North_Wall_Velocity_K= 0.0 ,
     End_of_Data_Block= .true. /
      
__The first **IBD** (the bar) __

   &Polyhedral_Immersed_Bodies
       Xi_1= 0.40 , Xj_1= 0.40 , Xk_1= 0.0 ,
       Xi_2= 0.40 , Xj_2= 0.40 , Xk_2= 0.0 ,
       Xi_3= 0.60 , Xj_3= 0.60 , Xk_3= 0.0 ,
       Xi_4= 0.60 , Xj_4= 0.60 , Xk_4= 0.0 ,
       Wall_BC_DataSetName     ="Set2"/

<note important>
   * The first **WBCT** is automatically linked to the walls of the domain end's (cavity walls). The wall temperatures are imposed as in the previous case. Moreover the I-velocity component is set to 1.5 at the top wall (Front_Wall_Velocity_I= 1.5).
   * The bar (immersed body) is linked to the second **WBCT**. The usual no-slip and impermeability boundary conditions are explicitly set here because the default values are no longer available due to the previous **WBCT**.
   * When several **WBCT** are present and at least one of them contains a specific set of velocity boundary conditions, which differ from the usual ones, the namelist [[sunfluidh:velocity_wall_boundary_condition_setup_namelist| "Velocity_Wall_Boundary_Condition_Setup"]] must then be explicitly set for each **WBCT**, even for the usual velocity boundary conditions.
</note>  

[[sunfluidh:heat_wall_boundary_condition_setup_namelist|Click here to come back to the "Heat_Wall_Boundary_Condition_Setup"]] \\
[[sunfluidh:velocity_wall_boundary_condition_setup_namelist|Click here to come back to the "Velocity_Wall_Boundary_Condition_Setup"]] \\
[[sunfluidh:species_wall_boundary_condition_setup_namelist|Click here to come back to the "Species_Wall_Boundary_Condition_Setup"]] \\
[[sunfluidh:polyhedral_immersed_bodies_setup_namelist|Click here to come back to the "Polyhedral_Immersed_Bodies"]] \\
[[sunfluidh:cylindrical_immersed_bodies_setup_namelist|Click here to come back to the "Cylindrical_Immersed_Bodies"]] \\