/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1606+                                  |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

// NACA 4314 - http://www.ppart.de/aerodynamics/profiles/NACA4.html

Lu -0.50;         //Upstream length from tip of leading edge
Ld  2.50;
Hy  0.50;         //Height of computational domain perpendicular to flow
Y1 #calc "$Hy*0.20";
Y2 #calc "$Hy*0.75";
Y3 #calc "$Hy*0.05";
Zw 0.025;
Y4 #calc "-1.00*$Y1";
Y5 #calc "-1.00*$Y2";
Y6 #calc "-1.00*$Y3";
Hd #calc "-1.00*$Hy";

X1 -0.20;
X2 0.50;
X3 1.25;

nx1  20;
nx2  40;
nx3  40;

ny1  50;
ny2  20;
nz   1;
nd   100;
NLE  10;
NTE  3;
n1   10;

vertices
(
  // Upper surface
  (0.000023   0.011008    0.000)        //  0
  (0.025207   0.042111    0.000)        //  1
  (0.423108   0.105449    0.000)        //  2
  (0.999494   0.001650    0.000)        //  3
  // Lower Surface
  (0.005456   -0.009554    0.000)       //  4
  (0.041212   -0.025379    0.000)       //  5
  (0.420457   -0.027870    0.000)       //  6 
  (0.999833   -0.001460    0.000)       //  7

  ($Lu   $Y1   0.000)                   //  8
  ($X1   $Y2   0.000)                   //  9
  ($X2   $Y2   0.000)                   //  10
  ($X3   $Y2   0.000)                   //  11
  ($Ld   $Y2   0.000)                   //  12
  ($Ld   $Y3   0.000)                   //  13
  ($Lu   $Hy   0.000)                   //  14
  ($X1   $Hy   0.000)                   //  15
  ($X2   $Hy   0.000)                   //  16
  ($X3   $Hy   0.000)                   //  17
  ($Ld   $Hy   0.000)                   //  18
  
  ($Lu   $Y4   0.000)                   //  19
  ($X1   $Y5   0.000)                   //  20
  ($X2   $Y5   0.000)                   //  21
  ($X3   $Y5   0.000)                   //  22
  ($Ld   $Y5   0.000)                   //  23
  ($Ld   $Y6   0.000)                   //  24
  ($Lu   $Hd   0.000)                   //  25
  ($X1   $Hd   0.000)                   //  26
  ($X2   $Hd   0.000)                   //  27
  ($X3   $Hd   0.000)                   //  28
  ($Ld   $Hd   0.000)                   //  29
  
  (0.000023   0.011008     $Zw )        // 30
  (0.025207   0.042111     $Zw )        // 31
  (0.423108   0.105449     $Zw )        // 32
  (0.999494   0.001650     $Zw )        // 33

  (0.005456   -0.009554    $Zw )        // 34
  (0.041212   -0.025379    $Zw )        // 35
  (0.420457   -0.027870    $Zw)         // 36
  (0.999833   -0.001460    $Zw )        // 37
  
  ($Lu   $Y1   $Zw)                     // 38
  ($X1   $Y2   $Zw)                     // 39
  ($X2   $Y2   $Zw)                     // 40
  ($X3   $Y2   $Zw)                     // 41
  ($Ld   $Y2   $Zw)                     // 42
  ($Ld   $Y3   $Zw)                     // 43
  ($Lu   $Hy   $Zw)                     // 44
  ($X1   $Hy   $Zw)                     // 45
  ($X2   $Hy   $Zw)                     // 46
  ($X3   $Hy   $Zw)                     // 47
  ($Ld   $Hy   $Zw)                     // 48
  
  ($Lu   $Y4   $Zw)                     // 49
  ($X1   $Y5   $Zw)                     // 50
  ($X2   $Y5   $Zw)                     // 51
  ($X3   $Y5   $Zw)                     // 52
  ($Ld   $Y5   $Zw)                     // 53
  ($Ld   $Y6   $Zw)                     // 54
  ($Lu   $Hd   $Zw)                     // 55
  ($X1   $Hd   $Zw)                     // 56
  ($X2   $Hd   $Zw)                     // 57
  ($X3   $Hd   $Zw)                     // 58
  ($Ld   $Hd   $Zw)                     // 59
);

edges
(
    // Upper body (upper circular arcs).
    arc 0   4  (0 0 0)
    arc 3   7  (1 0 0)
    arc 30  34 (0 0 $Zw)
    arc 33  37 (1 0 $Zw)
	
    // Upper body longitudinal splines.
    spline 0 1
    (
     (0.002175   0.016867    0.000)       
     (0.005767   0.022940    0.000)       
     (0.010804   0.029196    0.000)       
     (0.017286   0.035601    0.000)       
    )
    spline 30 31
    (
     (0.002175   0.016867     0.025 )       
     (0.005767   0.022940     0.025 )      
     (0.010804   0.029196     0.025 )       
     (0.017286   0.035601     0.025 )      
    )
	
    spline 1 2
    (
     (0.034559   0.048674    0.000)      
     (0.045326   0.055235    0.000)       
     (0.057488   0.061732    0.000)      
     (0.071018   0.068099    0.000)      
     (0.085883   0.074269    0.000)       
     (0.102044   0.080171    0.000)  
     (0.119458   0.085732    0.000) 
     (0.138071   0.090882    0.000) 
     (0.157827   0.095551    0.000)   
     (0.178662   0.099672    0.000)   
     (0.200507   0.103182    0.000)   
     (0.223287   0.106020    0.000)    
     (0.246922   0.108135    0.000)  
     (0.271330   0.109479    0.000)  
     (0.296422   0.110012    0.000) 
     (0.321054   0.109872    0.000) 
     (0.346008   0.109333    0.000) 
     (0.371384   0.108409    0.000) 
     (0.397108   0.107110    0.000)  
    )
	
    spline 31 32
    (
     (0.034559   0.048674     $Zw ) 
     (0.045326   0.055235     $Zw ) 
     (0.057488   0.061732     $Zw ) 
     (0.071018   0.068099     $Zw )  
     (0.085883   0.074269     $Zw )  
     (0.102044   0.080171     $Zw ) 
     (0.119458   0.085732     $Zw ) 
     (0.138071   0.090882     $Zw ) 
     (0.157827   0.095551     $Zw ) 
     (0.178662   0.099672     $Zw ) 
     (0.200507   0.103182     $Zw ) 
     (0.223287   0.106020     $Zw ) 
     (0.246922   0.108135     $Zw )  
     (0.271330   0.109479     $Zw ) 
     (0.296422   0.110012     $Zw ) 
     (0.321054   0.109872     $Zw ) 
     (0.346008   0.109333     $Zw )  
     (0.371384   0.108409     $Zw ) 
     (0.397108   0.107110     $Zw )  
    )
	
    spline 2 3
    (
     (0.449309   0.103443    0.000) 
     (0.475636   0.101108    0.000) 
     (0.502016   0.098465    0.000) 
     (0.528373   0.095535    0.000) 
     (0.554634   0.092341    0.000)  
     (0.580726   0.088904    0.000) 
     (0.606576   0.085252    0.000)  
     (0.632113   0.081408    0.000)
     (0.657267   0.077398    0.000) 
     (0.681969   0.073249    0.000) 
     (0.706152   0.068989    0.000) 
     (0.729751   0.064643    0.000) 
     (0.752702   0.060241    0.000) 
     (0.774943   0.055808    0.000) 
     (0.796416   0.051375    0.000)  
     (0.817063   0.046969    0.000) 
     (0.836830   0.042617    0.000) 
     (0.855665   0.038349    0.000) 
     (0.873519   0.034193    0.000) 
     (0.890346   0.030175    0.000) 
     (0.906103   0.026324    0.000)  
     (0.920749   0.022666    0.000)  
     (0.934246   0.019227    0.000)  
     (0.946562   0.016031    0.000)   
     (0.957664   0.013102    0.000)  
     (0.967526   0.010462    0.000)   
     (0.976122   0.008131    0.000)  
     (0.983432   0.006127    0.000)  
     (0.989437   0.004465    0.000) 
     (0.994123   0.003159    0.000) 
     (0.997478   0.002218    0.000)  
    )
	
    spline 32  33
    (
     (0.449309   0.103443     $Zw ) 
     (0.475636   0.101108     $Zw )  
     (0.502016   0.098465     $Zw )  
     (0.528373   0.095535     $Zw )   
     (0.554634   0.092341     $Zw )    
     (0.580726   0.088904     $Zw ) 
     (0.606576   0.085252     $Zw ) 
     (0.632113   0.081408     $Zw ) 
     (0.657267   0.077398     $Zw ) 
     (0.681969   0.073249     $Zw ) 
     (0.706152   0.068989     $Zw ) 
     (0.729751   0.064643     $Zw )  
     (0.752702   0.060241     $Zw )  
     (0.774943   0.055808     $Zw )  
     (0.796416   0.051375     $Zw )
     (0.817063   0.046969     $Zw )
     (0.836830   0.042617     $Zw ) 
     (0.855665   0.038349     $Zw )
     (0.873519   0.034193     $Zw ) 
     (0.890346   0.030175     $Zw )
     (0.906103   0.026324     $Zw )
     (0.920749   0.022666     $Zw ) 
     (0.934246   0.019227     $Zw ) 
     (0.946562   0.016031     $Zw ) 
     (0.957664   0.013102     $Zw ) 
     (0.967526   0.010462     $Zw )  
     (0.976122   0.008131     $Zw )  
     (0.983432   0.006127     $Zw )  
     (0.989437   0.004465     $Zw ) 
     (0.994123   0.003159     $Zw ) 
     (0.997478   0.002218     $Zw ) 
    )
	
    spline 4 5
    (
     (0.010137   -0.013618    0.000) 
     (0.016086   -0.017218    0.000)  
     (0.023270   -0.020368    0.000) 
     (0.031658   -0.023082    0.000) 
    )

	spline 34 35
    (
     (0.010137   -0.013618     $Zw )  
     (0.016086   -0.017218     $Zw ) 
     (0.023270   -0.020368     $Zw )  
     (0.031658   -0.023082     $Zw )  
    )
	
	spline 5 6
    (
     (0.051895   -0.027280    0.000) 
     (0.063667   -0.028810    0.000)
     (0.076486   -0.029994    0.000)  
     (0.090312   -0.030862    0.000) 
     (0.105100   -0.031446    0.000) 
     (0.120810   -0.031779    0.000)  
     (0.137397   -0.031898    0.000)
     (0.154822   -0.031841    0.000)  
     (0.173042   -0.031647    0.000)  
     (0.192017   -0.031359    0.000)  
     (0.211708   -0.031018    0.000)   
     (0.232074   -0.030669    0.000)    
     (0.253078   -0.030357    0.000)     
     (0.274680   -0.030127    0.000)  
     (0.296841   -0.030022    0.000)   
     (0.320578   -0.029942    0.000)  
     (0.344975   -0.029671    0.000)   
     (0.369797   -0.029222    0.000)  
     (0.394980   -0.028616    0.000)  
    )

	spline 35 36
    (
     (0.051895   -0.027280     $Zw )  
     (0.063667   -0.028810     $Zw )  
     (0.076486   -0.029994     $Zw )   
     (0.090312   -0.030862     $Zw )  
     (0.105100   -0.031446     $Zw )  
     (0.120810   -0.031779     $Zw )  
     (0.137397   -0.031898     $Zw ) 
     (0.154822   -0.031841     $Zw )  
     (0.173042   -0.031647     $Zw )  
     (0.192017   -0.031359     $Zw )  
     (0.211708   -0.031018     $Zw )  
     (0.232074   -0.030669     $Zw )  
     (0.253078   -0.030357     $Zw )  
     (0.274680   -0.030127     $Zw )  
     (0.296841   -0.030022     $Zw ) 
     (0.320578   -0.029942     $Zw ) 
     (0.344975   -0.029671     $Zw )   
     (0.369797   -0.029222     $Zw )   
     (0.394980   -0.028616     $Zw )   
    )
	
	spline 6 7
    (
     (0.446163   -0.027006    0.000)  
     (0.472028   -0.026042    0.000)  
     (0.497984   -0.024996    0.000)   
     (0.523963   -0.023887    0.000)   
     (0.549895   -0.022730    0.000)  
     (0.575709   -0.021542    0.000)  
     (0.601336   -0.020336    0.000)   
     (0.626706   -0.019124    0.000)    
     (0.651750   -0.017917    0.000)    
     (0.676399   -0.016724    0.000)   
     (0.700585   -0.015553    0.000)  
     (0.724240   -0.014411    0.000) 
     (0.747298   -0.013302    0.000) 
     (0.769696   -0.012231    0.000) 
     (0.791370   -0.011200    0.000)  
     (0.812258   -0.010214    0.000) 
     (0.851442   -0.008377    0.000)   
     (0.869626   -0.007531    0.000)    
     (0.886800   -0.006733    0.000)  
     (0.902914   -0.005986    0.000) 
     (0.917922   -0.005290    0.000)  
     (0.931779   -0.004648    0.000) 
     (0.944445   -0.004059    0.000)  
     (0.955881   -0.003527    0.000)  
     (0.966055   -0.003052    0.000)  
     (0.974935   -0.002636    0.000) 
     (0.982494   -0.002280    0.000)  
     (0.988711   -0.001987    0.000) 
     (0.993566   -0.001758    0.000)  
     (0.997044   -0.001593    0.000)   
     (0.999135   -0.001494    0.000)    
    )

	spline 36 37
    (
     (0.446163   -0.027006     $Zw )  
     (0.472028   -0.026042     $Zw ) 
     (0.497984   -0.024996     $Zw )   
     (0.523963   -0.023887     $Zw )   
     (0.549895   -0.022730     $Zw )  
     (0.575709   -0.021542     $Zw )   
     (0.601336   -0.020336     $Zw )   
     (0.626706   -0.019124     $Zw )   
     (0.651750   -0.017917     $Zw )  
     (0.700585   -0.015553     $Zw )  
     (0.724240   -0.014411     $Zw )   
     (0.747298   -0.013302     $Zw )  
     (0.769696   -0.012231     $Zw )  
     (0.791370   -0.011200     $Zw )  
     (0.812258   -0.010214     $Zw ) 
     (0.832301   -0.009272     $Zw )  
     (0.851442   -0.008377     $Zw ) 
     (0.869626   -0.007531     $Zw )  
     (0.886800   -0.006733     $Zw )  
     (0.902914   -0.005986     $Zw ) 
     (0.917922   -0.005290     $Zw )  
     (0.931779   -0.004648     $Zw ) 
     (0.944445   -0.004059     $Zw )  
     (0.955881   -0.003527     $Zw )  
     (0.966055   -0.003052     $Zw )  
     (0.974935   -0.002636     $Zw )  
     (0.982494   -0.002280     $Zw )   
     (0.988711   -0.001987     $Zw )  
     (0.993566   -0.001758     $Zw )  
     (0.997044   -0.001593     $Zw )   
     (0.999135   -0.001494     $Zw ) 
    )
);

blocks
(
 hex (0	 1	9	8	30	31	39	38) ($nx1 $ny1 $nz) simpleGrading (1 1 1)
 hex (1	 2	10	9	31	32	40	39) ($nx2 $ny1 $nz) simpleGrading (1 1 1) 
 hex (2	 3	11	10	32	33	41	40) ($nx3 $ny1 $nz) simpleGrading (1 1 1)
 hex (8	 9	15	14	38	39	45	44) ($nx1 $ny2 $nz) simpleGrading (1 1 1)
 hex (9	 10	16	15	39	40	46	45) ($nx2 $ny2 $nz) simpleGrading (1 1 1)
 hex (10	11	17	16	40	41	47	46) ($nx3 $ny2  $nz) simpleGrading (1 1 1)
 
 hex (4  0	8   19	34	30	38	49) ($NLE $ny1  $nz) simpleGrading (1 1 1)
 
 hex (3	7	24	13	33	37	54	43) ($NTE $nd  $nz) simpleGrading (1 1 1)
 
 hex (3	 13	12	11	33	43	42	41) ($nd $ny1 $nz) simpleGrading (1 1 1)
 hex (11 12	18	17	41	42	48	47) ($nd $ny2  $nz) simpleGrading (1 1 1)
 hex (7	 22	23	24	37	52	53	54) ($ny1 $nd  $nz) simpleGrading (1 1 1)
 hex (22  28  29	23	52	58	59	53) ($ny2 $nd  $nz) simpleGrading (1 1 1)
 
 hex (5	4	19	20	35	34	49	50) ($nx1 $ny1  $nz) simpleGrading (1 1 1)
 hex (6	5	20	21	36	35	50	51) ($nx2 $ny1  $nz) simpleGrading (1 1 1)
 hex (7	6	21	22	37	36	51	52) ($nx3 $ny1  $nz) simpleGrading (1 1 1)
 hex (20	19	25	26	50	49	55	56) ($nx1 $ny2  $nz) simpleGrading (1 1 1)
 hex (21	20	26	27	51	50	56	57) ($nx2 $ny2  $nz) simpleGrading (1 1 1)
 hex (22	21	27	28	52	51	57	58) ($nx3 $ny2  $nz) simpleGrading (1 1 1)
);

boundary (
 Inlet  {
  type patch;
   faces (
    (25	19	55	49) (19	 8	49	38) (8	14	38	44) 
   );
 }
 
 Outlet {
  type patch;
   faces (
    (29	23	59	53) (23	24	53	54) (24	13	54	43) (13	12	43	42) (12	18	42	48) 
   );
 }
 
 Aerofoil {
  type wall;
   faces (
    (0	1	30	31) (1	2	31	32) (2	3	32	33) (3	7	33	37) 
	(7	6	37	36) (6	5	36	35) (5	4	35	34) (4	0	34	30)
   );
 }
 
 Upper {
  type symmetry;
   faces (
    (14	15	44	45) (15	16	45	46) (16	17	46	47) (17	18	47	48) 
   );
 }
 
 Lower {
  type symmetry;
   faces (
    (25	26	55	56) (26	27	56	57) (27	28	57	58) (28	29	58	59)
   );
 }
);
 
mergePatchPairs ();


// ************************************************************************* //