CG  Version 25
commonMacros.h
Go to the documentation of this file.
1 c **********************************************************************
2 c This file contains some commonly used macros.
3 c **********************************************************************
4 
5 
6 c Define macros for the derivatives based on the dimension, order of accuracy and grid-type
7 #beginMacro defineDerivativeMacros(DIM,ORDER,GRIDTYPE)
8 
9 #defineMacro U(cc) u(i1,i2,i3,cc)
10 #defineMacro UU(cc) uu(i1,i2,i3,cc)
11 
12 #If #DIM == "2"
13  #If #ORDER == "2"
14  #If #GRIDTYPE == "rectangular"
15  #defineMacro UX(cc) ux22r(i1,i2,i3,cc)
16  #defineMacro UY(cc) uy22r(i1,i2,i3,cc)
17  #defineMacro UXX(cc) uxx22r(i1,i2,i3,cc)
18  #defineMacro UXY(cc) uxy22r(i1,i2,i3,cc)
19  #defineMacro UYY(cc) uyy22r(i1,i2,i3,cc)
20  #defineMacro ULAP(cc) ulaplacian22r(i1,i2,i3,cc)
21  #Elif #GRIDTYPE == "curvilinear"
22  #defineMacro UX(cc) ux22(i1,i2,i3,cc)
23  #defineMacro UY(cc) uy22(i1,i2,i3,cc)
24  #defineMacro UXX(cc) uxx22(i1,i2,i3,cc)
25  #defineMacro UXY(cc) uxy22(i1,i2,i3,cc)
26  #defineMacro UYY(cc) uyy22(i1,i2,i3,cc)
27  #defineMacro ULAP(cc) ulaplacian22(i1,i2,i3,cc)
28 
29  #defineMacro RXX() rxx22(i1,i2,i3)
30  #defineMacro RXY() rxy22(i1,i2,i3)
31  #defineMacro RYY() ryy22(i1,i2,i3)
32  #defineMacro SXX() sxx22(i1,i2,i3)
33  #defineMacro SXY() sxy22(i1,i2,i3)
34  #defineMacro SYY() syy22(i1,i2,i3)
35  #Else
36  stop 888
37  #End
38  #Else
39  #If #GRIDTYPE == "rectangular"
40  #defineMacro UX(cc) ux42r(i1,i2,i3,cc)
41  #defineMacro UY(cc) uy42r(i1,i2,i3,cc)
42  #defineMacro UXX(cc) uxx42r(i1,i2,i3,cc)
43  #defineMacro UXY(cc) uxy42r(i1,i2,i3,cc)
44  #defineMacro UYY(cc) uyy42r(i1,i2,i3,cc)
45  #defineMacro ULAP(cc) ulaplacian42r(i1,i2,i3,cc)
46  #Elif #GRIDTYPE == "curvilinear"
47  #defineMacro UX(cc) ux42(i1,i2,i3,cc)
48  #defineMacro UY(cc) uy42(i1,i2,i3,cc)
49  #defineMacro UXX(cc) uxx42(i1,i2,i3,cc)
50  #defineMacro UXY(cc) uxy42(i1,i2,i3,cc)
51  #defineMacro UYY(cc) uyy42(i1,i2,i3,cc)
52  #defineMacro ULAP(cc) ulaplacian42(i1,i2,i3,cc)
53 
54  #defineMacro RXX() rxx42(i1,i2,i3)
55  #defineMacro RXY() rxy42(i1,i2,i3)
56  #defineMacro RYY() ryy42(i1,i2,i3)
57  #defineMacro SXX() sxx42(i1,i2,i3)
58  #defineMacro SXY() sxy42(i1,i2,i3)
59  #defineMacro SYY() syy42(i1,i2,i3)
60  #Else
61  stop 888
62  #End
63  #End
64 #Else
65  #If #ORDER == "2"
66  #If #GRIDTYPE == "rectangular"
67  #defineMacro UX(cc) ux23r(i1,i2,i3,cc)
68  #defineMacro UY(cc) uy23r(i1,i2,i3,cc)
69  #defineMacro UZ(cc) uz23r(i1,i2,i3,cc)
70  #defineMacro UXX(cc) uxx23r(i1,i2,i3,cc)
71  #defineMacro UXY(cc) uxy23r(i1,i2,i3,cc)
72  #defineMacro UXZ(cc) uxz23r(i1,i2,i3,cc)
73  #defineMacro UYY(cc) uyy23r(i1,i2,i3,cc)
74  #defineMacro UYZ(cc) uyz23r(i1,i2,i3,cc)
75  #defineMacro UZZ(cc) uzz23r(i1,i2,i3,cc)
76  #defineMacro ULAP(cc) ulaplacian23r(i1,i2,i3,cc)
77  #Elif #GRIDTYPE == "curvilinear"
78  #defineMacro UX(cc) ux23(i1,i2,i3,cc)
79  #defineMacro UY(cc) uy23(i1,i2,i3,cc)
80  #defineMacro UZ(cc) uz23(i1,i2,i3,cc)
81  #defineMacro UXX(cc) uxx23(i1,i2,i3,cc)
82  #defineMacro UXY(cc) uxy23(i1,i2,i3,cc)
83  #defineMacro UXZ(cc) uxz23(i1,i2,i3,cc)
84  #defineMacro UYY(cc) uyy23(i1,i2,i3,cc)
85  #defineMacro UYZ(cc) uyz23(i1,i2,i3,cc)
86  #defineMacro UZZ(cc) uzz23(i1,i2,i3,cc)
87  #defineMacro ULAP(cc) ulaplacian23(i1,i2,i3,cc)
88 
89  #defineMacro RXX() rxx23(i1,i2,i3)
90  #defineMacro RXY() rxy23(i1,i2,i3)
91  #defineMacro RXZ() rxz23(i1,i2,i3)
92  #defineMacro RYY() ryy23(i1,i2,i3)
93  #defineMacro RYZ() ryz23(i1,i2,i3)
94  #defineMacro RZZ() rzz23(i1,i2,i3)
95 
96  #defineMacro SXX() sxx23(i1,i2,i3)
97  #defineMacro SXY() sxy23(i1,i2,i3)
98  #defineMacro SXZ() sxz23(i1,i2,i3)
99  #defineMacro SYY() syy23(i1,i2,i3)
100  #defineMacro SYZ() syz23(i1,i2,i3)
101  #defineMacro SZZ() szz23(i1,i2,i3)
102 
103  #defineMacro TXX() txx23(i1,i2,i3)
104  #defineMacro TXY() txy23(i1,i2,i3)
105  #defineMacro TXZ() txz23(i1,i2,i3)
106  #defineMacro TYY() tyy23(i1,i2,i3)
107  #defineMacro TYZ() tyz23(i1,i2,i3)
108  #defineMacro TZZ() tzz23(i1,i2,i3)
109 
110  #Else
111  stop 888
112  #End
113 
114  #Else
115 
116  #If #GRIDTYPE == "rectangular"
117  #defineMacro UX(cc) ux43r(i1,i2,i3,cc)
118  #defineMacro UY(cc) uy43r(i1,i2,i3,cc)
119  #defineMacro UZ(cc) uz43r(i1,i2,i3,cc)
120  #defineMacro UXX(cc) uxx43r(i1,i2,i3,cc)
121  #defineMacro UXY(cc) uxy43r(i1,i2,i3,cc)
122  #defineMacro UXZ(cc) uxz43r(i1,i2,i3,cc)
123  #defineMacro UYY(cc) uyy43r(i1,i2,i3,cc)
124  #defineMacro UYZ(cc) uyz43r(i1,i2,i3,cc)
125  #defineMacro UZZ(cc) uzz43r(i1,i2,i3,cc)
126  #defineMacro ULAP(cc) ulaplacian43r(i1,i2,i3,cc)
127  #Elif #GRIDTYPE == "curvilinear"
128  #defineMacro UX(cc) ux43(i1,i2,i3,cc)
129  #defineMacro UY(cc) uy43(i1,i2,i3,cc)
130  #defineMacro UZ(cc) uz43(i1,i2,i3,cc)
131  #defineMacro UXX(cc) uxx43(i1,i2,i3,cc)
132  #defineMacro UXY(cc) uxy43(i1,i2,i3,cc)
133  #defineMacro UXZ(cc) uxz43(i1,i2,i3,cc)
134  #defineMacro UYY(cc) uyy43(i1,i2,i3,cc)
135  #defineMacro UYZ(cc) uyz43(i1,i2,i3,cc)
136  #defineMacro UZZ(cc) uzz43(i1,i2,i3,cc)
137  #defineMacro ULAP(cc) ulaplacian43(i1,i2,i3,cc)
138 
139  #defineMacro RXX() rxx43(i1,i2,i3)
140  #defineMacro RXY() rxy43(i1,i2,i3)
141  #defineMacro RXZ() rxz43(i1,i2,i3)
142  #defineMacro RYY() ryy43(i1,i2,i3)
143  #defineMacro RYZ() ryz43(i1,i2,i3)
144  #defineMacro RZZ() rzz43(i1,i2,i3)
145 
146  #defineMacro SXX() sxx43(i1,i2,i3)
147  #defineMacro SXY() sxy43(i1,i2,i3)
148  #defineMacro SXZ() sxz43(i1,i2,i3)
149  #defineMacro SYY() syy43(i1,i2,i3)
150  #defineMacro SYZ() syz43(i1,i2,i3)
151  #defineMacro SZZ() szz43(i1,i2,i3)
152 
153  #defineMacro TXX() txx43(i1,i2,i3)
154  #defineMacro TXY() txy43(i1,i2,i3)
155  #defineMacro TXZ() txz43(i1,i2,i3)
156  #defineMacro TYY() tyy43(i1,i2,i3)
157  #defineMacro TYZ() tyz43(i1,i2,i3)
158  #defineMacro TZZ() tzz43(i1,i2,i3)
159  #Else
160  stop 888
161  #End
162  #End
163 #End
164 #endMacro
165 
166 
167 c================================================================
168 c Input:
169 c ADTYPE: AD2, AD4, AD24
170 c TURBULENCE_MODEL: INS, SPAL
171 c Output:
172 c artificialDissipation(cc) : inline macro for u,v,w
173 c artificialDissipationTM(cc) : inline macro for "n" or k or epsilon
174 c================================================================
175 #beginMacro defineArtificialDissipationMacro(ADTYPE,DIM,TURBULENCE_MODEL)
176 
177  ! By default there is no AD:
178  #defineMacro artificialDissipation(cc)
179  #defineMacro artificialDissipationTM(cc)
180 
181  ! 2nd-order artificial dissipation
182  #If #ADTYPE == "AD2"
183  #If #DIM == "2"
184  #defineMacro artificialDissipation(cc) +adCoeff2*delta22(cc)
185  #Else
186  #defineMacro artificialDissipation(cc) +adCoeff2*delta23(cc)
187  #End
188 
189  #End
190 
191  ! 4th-order artficial dissipation **todo** implicit-line, self-adjoint versions
192  #If #ADTYPE == "AD4" || #ADTYPE == "AD24"
193  #If #DIM == "2"
194  #defineMacro artificialDissipation(cc) +adCoeff4*delta42(cc)
195  #Else
196  #defineMacro artificialDissipation(cc) +adCoeff4*delta43(cc)
197  #End
198  #End
199 
200  ! Both 2nd and 4th order dissipation
201  #If #ADTYPE == "AD24"
202  #If #DIM == "2"
203  #defineMacro artificialDissipation(cc) +adCoeff2*delta22(cc) + adCoeff4*delta42(cc)
204  #Else
205  #defineMacro artificialDissipation(cc) +adCoeff2*delta23(cc) + adCoeff4*delta43(cc)
206  #End
207  #End
208 
209  #If #TURBULENCE_MODEL == "SPAL"
210  ! Define the ad macro for the turbulent eddy viscosity equation
211  ! Just base the coefficient on the derivatives of that component.
212  #If #ADTYPE == "AD2"
213  #If #DIM == "2"
214  #defineMacro artificialDissipationTM(cc) +( ad21n+cd22n*( abs(UX(cc))+abs(UY(cc)) ) )*delta22(cc)
215  #Else
216  #defineMacro artificialDissipationTM(cc) +( ad21n+cd22n*( abs(UX(cc))+abs(UY(cc))+abs(UZ(cc)) ) )*delta23(cc)
217  #End
218  #Elif #ADTYPE == "AD4"
219  #If #DIM == "2"
220  #defineMacro artificialDissipationTM(cc) +( ad41n+cd42n*( abs(UX(cc))+abs(UY(cc)) ) )*delta42(cc)
221  #Else
222  #defineMacro artificialDissipationTM(cc) +( ad41n+cd42n*( abs(UX(nc))+abs(UY(nc))+abs(UZ(nc)) ) )*delta43(cc)
223  #End
224  #Elif #ADTYPE == "AD24"
225  #If #DIM == "2"
226  #defineMacro artificialDissipationTM(cc) +( ad21n+cd22n*( abs(UX(cc))+abs(UY(cc)) ) )*delta22(cc)\
227  +( ad41n+cd42n*( abs(UX(cc))+abs(UY(cc)) ) )*delta42(cc)
228  #Else
229  #defineMacro artificialDissipationTM(cc) +( ad21n+cd22n*( abs(UX(cc))+abs(UY(cc))+abs(UZ(cc)) ) )*delta23(cc)\
230  +( ad41n+cd42n*( abs(UX(nc))+abs(UY(nc))+abs(UZ(nc)) ) )*delta43(cc)
231  #End
232  #Else
233  stop 4444
234  #End
235 
236  #End
237 #endMacro
238 
239 
240 c================================================================
241 c Input:
242 c ADTYPE: AD2, AD4, AD24
243 c TURBULENCE_MODEL: INS, SPAL
244 c Output:
245 c artificialDissipation(cc) : inline macro for u,v,w
246 c artificialDissipationTM(cc) : inline macro for "n" or k or epsilon
247 c================================================================
248 #beginMacro getArtificialDissipationCoeff(ADTYPE,DIM,TURBULENCE_MODEL)
249 
250  #If #ADTYPE == "AD2" || #ADTYPE == "AD24"
251  #If #DIM == "2"
252  adCoeff2 = ad21+cd22*( abs(UX(uc))+abs(UY(uc))+abs(UX(vc))+abs(UY(vc)) )
253  #Else
254  adCoeff2 = ad21+cd22*( abs(UX(uc))+abs(UY(uc))+abs(UZ(uc))+\
255  abs(UX(vc))+abs(UY(vc))+abs(UZ(vc))+\
256  abs(UX(wc))+abs(UY(wc))+abs(UZ(wc)))
257  #End
258 
259  #End
260 
261  ! Artficial Dissipation **todo** implicit-line, self-adjoint versions
262  #If #ADTYPE == "AD4" || #ADTYPE == "AD24"
263  #If #DIM == "2"
264  adCoeff4 = ad41+cd42*( abs(UX(uc))+abs(UY(uc))+abs(UX(vc))+abs(UY(vc)) )
265  #Else
266  adCoeff4 = ad41+cd42*( abs(UX(uc))+abs(UY(uc))+abs(UZ(uc))+\
267  abs(UX(vc))+abs(UY(vc))+abs(UZ(vc))+\
268  abs(UX(wc))+abs(UY(wc))+abs(UZ(wc)))
269  #End
270  #End
271 
272 #endMacro
273 
274 
Generated on Fri Jan 4 2013 10:19:17 for CG by   doxygen 1.8.3