Rocstar  1.0
Rocstar multiphysics simulation application
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RFLO_RoeFluxSecond.F90
Go to the documentation of this file.
1 ! *********************************************************************
2 ! * Rocstar Simulation Suite *
3 ! * Copyright@2015, Illinois Rocstar LLC. All rights reserved. *
4 ! * *
5 ! * Illinois Rocstar LLC *
6 ! * Champaign, IL *
7 ! * www.illinoisrocstar.com *
8 ! * sales@illinoisrocstar.com *
9 ! * *
10 ! * License: See LICENSE file in top level of distribution package or *
11 ! * http://opensource.org/licenses/NCSA *
12 ! *********************************************************************
13 ! *********************************************************************
14 ! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, *
15 ! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES *
16 ! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND *
17 ! * NONINFRINGEMENT. IN NO EVENT SHALL THE CONTRIBUTORS OR *
18 ! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER *
19 ! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, *
20 ! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE *
21 ! * USE OR OTHER DEALINGS WITH THE SOFTWARE. *
22 ! *********************************************************************
23 !******************************************************************************
24 !
25 ! Purpose: compute convective fluxes based on 2nd-order Roe upwind scheme.
26 !
27 ! Description: none.
28 !
29 ! Input: region = data of current region.
30 !
31 ! Output: region%levels%mixt%rhs = convective fluxes added to the residual.
32 !
33 ! Notes: uses MUSCL scheme with kappa=1/3.
34 !
35 !******************************************************************************
36 !
37 ! $Id: RFLO_RoeFluxSecond.F90,v 1.3 2008/12/06 08:44:28 mtcampbe Exp $
38 !
39 ! Copyright: (c) 2003 by the University of Illinois
40 !
41 !******************************************************************************
42 
43 SUBROUTINE rflo_roefluxsecond( region )
44 
45  USE moddatatypes
46  USE moddatastruct, ONLY : t_region
47  USE modinterfaces, ONLY : rflo_getdimensphys, rflo_getcelloffset, &
48  rflo_getnodeoffset, rflo_roefluxpatch
49  USE moderror
50  USE modparameters
51  IMPLICIT NONE
52 
53 #include "Indexing.h"
54 
55 ! ... parameters
56  TYPE(t_region) :: region
57 
58 ! ... loop variables
59  INTEGER :: i, j, k, ipatch
60 
61 ! ... local variables
62  INTEGER :: ilev, ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend
63  INTEGER :: icoff, ijcoff, inoff, ijnoff, ijkc0, ijkc1, ijkc2, ijkc3, ijkn
64  INTEGER :: indcp, indmol, indsvel
65 
66  REAL(RFREAL) :: limfac, limfac3, rvolref, vola, eps2(3), dvar(5), dvarm(5), &
67  dvarp(5), rhl, rhr, qsl, qsr, pav, deltl(5), deltr(5), fc(5),&
68  eps2n, rgas, gam, ggm1, rl, rr, ul, ur, vl, vr, wl, wr, &
69  pl, pr, hl, hr, qsrl, qsrr
70  REAL(RFREAL), POINTER :: cv(:,:), dv(:,:), rhs(:,:), si(:,:), sj(:,:), sk(:,:)
71  REAL(RFREAL), POINTER :: gv(:,:), vol(:), sivel(:), sjvel(:), skvel(:)
72 
73 ! ... functions
74  REAL(RFREAL) :: muscl3, af, bf, eps
75 
76 !******************************************************************************
77 ! limiter function
78 
79  muscl3(af,bf,eps) = (bf*(2._rfreal*af*af+eps)+af*(bf*bf+2._rfreal*eps))/ &
80  (2._rfreal*af*af+2._rfreal*bf*bf-af*bf+ &
81  3._rfreal*eps+1.e-30_rfreal)
82 
83  CALL registerfunction( region%global,'RFLO_RoeFluxSecond',&
84  'RFLO_RoeFluxSecond.F90' )
85 
86 ! get dimensions and pointers -------------------------------------------------
87 
88  ilev = region%currLevel
89 
90  CALL rflo_getdimensphys( region,ilev,ipcbeg,ipcend, &
91  jpcbeg,jpcend,kpcbeg,kpcend )
92  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )
93  CALL rflo_getnodeoffset( region,ilev,inoff,ijnoff )
94 
95  cv => region%levels(ilev)%mixt%cv
96  dv => region%levels(ilev)%mixt%dv
97  gv => region%levels(ilev)%mixt%gv
98  rhs => region%levels(ilev)%mixt%rhs
99  vol => region%levels(ilev)%grid%vol
100  si => region%levels(ilev)%grid%si
101  sj => region%levels(ilev)%grid%sj
102  sk => region%levels(ilev)%grid%sk
103  sivel => region%levels(ilev)%grid%siVel
104  sjvel => region%levels(ilev)%grid%sjVel
105  skvel => region%levels(ilev)%grid%skVel
106  indsvel = region%levels(ilev)%grid%indSvel
107  indcp = region%levels(ilev)%mixt%indCp
108  indmol = region%levels(ilev)%mixt%indMol
109  limfac = region%mixtInput%limFac
110 
111 ! normalise epsilon^2 for all limited variables (rho, u, v, w, p) -------------
112 
113  limfac3 = limfac*limfac*limfac
114  rvolref = 1._rfreal/region%global%limVolRef**1.5_rfreal
115  eps2(1) = limfac3*region%global%limRef(1)*region%global%limRef(1)*rvolref
116  eps2(2) = limfac3*region%global%limRef(2)*region%global%limRef(2)*rvolref
117  eps2(3) = limfac3*region%global%limRef(3)*region%global%limRef(3)*rvolref
118 
119 ! flux in i-direction (except through boundary) -------------------------------
120 
121  DO k=kpcbeg,kpcend
122  DO j=jpcbeg,jpcend
123  DO i=ipcbeg+1,ipcend
124  ijkc0 = indijk(i ,j,k,icoff,ijcoff)
125  ijkc1 = indijk(i-1,j,k,icoff,ijcoff)
126  ijkc2 = indijk(i-2,j,k,icoff,ijcoff)
127  ijkc3 = indijk(i+1,j,k,icoff,ijcoff)
128  ijkn = indijk(i ,j,k,inoff,ijnoff)
129 
130  vola = (0.5_rfreal*(vol(ijkc0)+vol(ijkc1)))**1.5_rfreal
131 
132 ! ----- limited extrapolations
133 
134  dvar(1) = cv(cv_mixt_dens,ijkc0) - cv(cv_mixt_dens,ijkc1)
135  dvar(2) = dv(dv_mixt_uvel,ijkc0) - dv(dv_mixt_uvel,ijkc1)
136  dvar(3) = dv(dv_mixt_vvel,ijkc0) - dv(dv_mixt_vvel,ijkc1)
137  dvar(4) = dv(dv_mixt_wvel,ijkc0) - dv(dv_mixt_wvel,ijkc1)
138  dvar(5) = dv(dv_mixt_pres,ijkc0) - dv(dv_mixt_pres,ijkc1)
139 
140  dvarm(1) = cv(cv_mixt_dens,ijkc1) - cv(cv_mixt_dens,ijkc2)
141  dvarm(2) = dv(dv_mixt_uvel,ijkc1) - dv(dv_mixt_uvel,ijkc2)
142  dvarm(3) = dv(dv_mixt_vvel,ijkc1) - dv(dv_mixt_vvel,ijkc2)
143  dvarm(4) = dv(dv_mixt_wvel,ijkc1) - dv(dv_mixt_wvel,ijkc2)
144  dvarm(5) = dv(dv_mixt_pres,ijkc1) - dv(dv_mixt_pres,ijkc2)
145 
146  dvarp(1) = cv(cv_mixt_dens,ijkc3) - cv(cv_mixt_dens,ijkc0)
147  dvarp(2) = dv(dv_mixt_uvel,ijkc3) - dv(dv_mixt_uvel,ijkc0)
148  dvarp(3) = dv(dv_mixt_vvel,ijkc3) - dv(dv_mixt_vvel,ijkc0)
149  dvarp(4) = dv(dv_mixt_wvel,ijkc3) - dv(dv_mixt_wvel,ijkc0)
150  dvarp(5) = dv(dv_mixt_pres,ijkc3) - dv(dv_mixt_pres,ijkc0)
151 
152  eps2n = eps2(1)*vola
153  deltl(1) = 0.5_rfreal*muscl3(dvar(1) ,dvarm(1),eps2n)
154  deltr(1) = 0.5_rfreal*muscl3(dvarp(1),dvar(1) ,eps2n)
155  eps2n = eps2(2)*vola
156  deltl(2) = 0.5_rfreal*muscl3(dvar(2) ,dvarm(2),eps2n)
157  deltr(2) = 0.5_rfreal*muscl3(dvarp(2),dvar(2) ,eps2n)
158  deltl(3) = 0.5_rfreal*muscl3(dvar(3) ,dvarm(3),eps2n)
159  deltr(3) = 0.5_rfreal*muscl3(dvarp(3),dvar(3) ,eps2n)
160  deltl(4) = 0.5_rfreal*muscl3(dvar(4) ,dvarm(4),eps2n)
161  deltr(4) = 0.5_rfreal*muscl3(dvarp(4),dvar(4) ,eps2n)
162  eps2n = eps2(3)*vola
163  deltl(5) = 0.5_rfreal*muscl3(dvar(5) ,dvarm(5),eps2n)
164  deltr(5) = 0.5_rfreal*muscl3(dvarp(5),dvar(5) ,eps2n)
165 
166 ! ----- left and right states
167 
168  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc1*indmol)
169  gam = gv(gv_mixt_cp,ijkc1*indcp)/(gv(gv_mixt_cp,ijkc1*indcp)-rgas)
170  ggm1 = gam/(gam-1._rfreal)
171  rl = cv(cv_mixt_dens,ijkc1) + deltl(1)
172  ul = dv(dv_mixt_uvel,ijkc1) + deltl(2)
173  vl = dv(dv_mixt_vvel,ijkc1) + deltl(3)
174  wl = dv(dv_mixt_wvel,ijkc1) + deltl(4)
175  pl = dv(dv_mixt_pres,ijkc1) + deltl(5)
176  hl = ggm1*pl/rl + 0.5_rfreal*(ul*ul+vl*vl+wl*wl)
177  qsrl = (ul*si(xcoord,ijkn)+vl*si(ycoord,ijkn)+ &
178  wl*si(zcoord,ijkn)-sivel(ijkn*indsvel))*rl
179 
180  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc0*indmol)
181  gam = gv(gv_mixt_cp,ijkc0*indcp)/(gv(gv_mixt_cp,ijkc0*indcp)-rgas)
182  ggm1 = gam/(gam-1._rfreal)
183  rr = cv(cv_mixt_dens,ijkc0) - deltr(1)
184  ur = dv(dv_mixt_uvel,ijkc0) - deltr(2)
185  vr = dv(dv_mixt_vvel,ijkc0) - deltr(3)
186  wr = dv(dv_mixt_wvel,ijkc0) - deltr(4)
187  pr = dv(dv_mixt_pres,ijkc0) - deltr(5)
188  hr = ggm1*pr/rr + 0.5_rfreal*(ur*ur+vr*vr+wr*wr)
189  qsrr = (ur*si(xcoord,ijkn)+vr*si(ycoord,ijkn)+ &
190  wr*si(zcoord,ijkn)-sivel(ijkn*indsvel))*rr
191 
192 ! ----- fluxes
193 
194  pav = 0.5_rfreal*(pl+pr)
195  fc(1) = 0.5_rfreal*(qsrl +qsrr )
196  fc(2) = 0.5_rfreal*(qsrl*ul+qsrr*ur) + si(xcoord,ijkn)*pav
197  fc(3) = 0.5_rfreal*(qsrl*vl+qsrr*vr) + si(ycoord,ijkn)*pav
198  fc(4) = 0.5_rfreal*(qsrl*wl+qsrr*wr) + si(zcoord,ijkn)*pav
199  fc(5) = 0.5_rfreal*(qsrl*hl+qsrr*hr) + sivel(ijkn*indsvel)*pav
200 
201  rhs(cv_mixt_dens,ijkc0) = rhs(cv_mixt_dens,ijkc0) + fc(1)
202  rhs(cv_mixt_xmom,ijkc0) = rhs(cv_mixt_xmom,ijkc0) + fc(2)
203  rhs(cv_mixt_ymom,ijkc0) = rhs(cv_mixt_ymom,ijkc0) + fc(3)
204  rhs(cv_mixt_zmom,ijkc0) = rhs(cv_mixt_zmom,ijkc0) + fc(4)
205  rhs(cv_mixt_ener,ijkc0) = rhs(cv_mixt_ener,ijkc0) + fc(5)
206 
207  rhs(cv_mixt_dens,ijkc1) = rhs(cv_mixt_dens,ijkc1) - fc(1)
208  rhs(cv_mixt_xmom,ijkc1) = rhs(cv_mixt_xmom,ijkc1) - fc(2)
209  rhs(cv_mixt_ymom,ijkc1) = rhs(cv_mixt_ymom,ijkc1) - fc(3)
210  rhs(cv_mixt_zmom,ijkc1) = rhs(cv_mixt_zmom,ijkc1) - fc(4)
211  rhs(cv_mixt_ener,ijkc1) = rhs(cv_mixt_ener,ijkc1) - fc(5)
212  ENDDO ! i
213  ENDDO ! j
214  ENDDO ! k
215 
216 ! flux in j-direction (except through boundary) -------------------------------
217 
218  DO k=kpcbeg,kpcend
219  DO j=jpcbeg+1,jpcend
220  DO i=ipcbeg,ipcend
221  ijkc0 = indijk(i,j ,k,icoff,ijcoff)
222  ijkc1 = indijk(i,j-1,k,icoff,ijcoff)
223  ijkc2 = indijk(i,j-2,k,icoff,ijcoff)
224  ijkc3 = indijk(i,j+1,k,icoff,ijcoff)
225  ijkn = indijk(i,j ,k,inoff,ijnoff)
226 
227  vola = (0.5_rfreal*(vol(ijkc0)+vol(ijkc1)))**1.5_rfreal
228 
229 ! ----- limited extrapolations
230 
231  dvar(1) = cv(cv_mixt_dens,ijkc0) - cv(cv_mixt_dens,ijkc1)
232  dvar(2) = dv(dv_mixt_uvel,ijkc0) - dv(dv_mixt_uvel,ijkc1)
233  dvar(3) = dv(dv_mixt_vvel,ijkc0) - dv(dv_mixt_vvel,ijkc1)
234  dvar(4) = dv(dv_mixt_wvel,ijkc0) - dv(dv_mixt_wvel,ijkc1)
235  dvar(5) = dv(dv_mixt_pres,ijkc0) - dv(dv_mixt_pres,ijkc1)
236 
237  dvarm(1) = cv(cv_mixt_dens,ijkc1) - cv(cv_mixt_dens,ijkc2)
238  dvarm(2) = dv(dv_mixt_uvel,ijkc1) - dv(dv_mixt_uvel,ijkc2)
239  dvarm(3) = dv(dv_mixt_vvel,ijkc1) - dv(dv_mixt_vvel,ijkc2)
240  dvarm(4) = dv(dv_mixt_wvel,ijkc1) - dv(dv_mixt_wvel,ijkc2)
241  dvarm(5) = dv(dv_mixt_pres,ijkc1) - dv(dv_mixt_pres,ijkc2)
242 
243  dvarp(1) = cv(cv_mixt_dens,ijkc3) - cv(cv_mixt_dens,ijkc0)
244  dvarp(2) = dv(dv_mixt_uvel,ijkc3) - dv(dv_mixt_uvel,ijkc0)
245  dvarp(3) = dv(dv_mixt_vvel,ijkc3) - dv(dv_mixt_vvel,ijkc0)
246  dvarp(4) = dv(dv_mixt_wvel,ijkc3) - dv(dv_mixt_wvel,ijkc0)
247  dvarp(5) = dv(dv_mixt_pres,ijkc3) - dv(dv_mixt_pres,ijkc0)
248 
249  eps2n = eps2(1)*vola
250  deltl(1) = 0.5_rfreal*muscl3(dvar(1) ,dvarm(1),eps2n)
251  deltr(1) = 0.5_rfreal*muscl3(dvarp(1),dvar(1) ,eps2n)
252  eps2n = eps2(2)*vola
253  deltl(2) = 0.5_rfreal*muscl3(dvar(2) ,dvarm(2),eps2n)
254  deltr(2) = 0.5_rfreal*muscl3(dvarp(2),dvar(2) ,eps2n)
255  deltl(3) = 0.5_rfreal*muscl3(dvar(3) ,dvarm(3),eps2n)
256  deltr(3) = 0.5_rfreal*muscl3(dvarp(3),dvar(3) ,eps2n)
257  deltl(4) = 0.5_rfreal*muscl3(dvar(4) ,dvarm(4),eps2n)
258  deltr(4) = 0.5_rfreal*muscl3(dvarp(4),dvar(4) ,eps2n)
259  eps2n = eps2(3)*vola
260  deltl(5) = 0.5_rfreal*muscl3(dvar(5) ,dvarm(5),eps2n)
261  deltr(5) = 0.5_rfreal*muscl3(dvarp(5),dvar(5) ,eps2n)
262 
263 ! ----- left and right states
264 
265  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc1*indmol)
266  gam = gv(gv_mixt_cp,ijkc1*indcp)/(gv(gv_mixt_cp,ijkc1*indcp)-rgas)
267  ggm1 = gam/(gam-1._rfreal)
268  rl = cv(cv_mixt_dens,ijkc1) + deltl(1)
269  ul = dv(dv_mixt_uvel,ijkc1) + deltl(2)
270  vl = dv(dv_mixt_vvel,ijkc1) + deltl(3)
271  wl = dv(dv_mixt_wvel,ijkc1) + deltl(4)
272  pl = dv(dv_mixt_pres,ijkc1) + deltl(5)
273  hl = ggm1*pl/rl + 0.5_rfreal*(ul*ul+vl*vl+wl*wl)
274  qsrl = (ul*sj(xcoord,ijkn)+vl*sj(ycoord,ijkn)+ &
275  wl*sj(zcoord,ijkn)-sjvel(ijkn*indsvel))*rl
276 
277  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc0*indmol)
278  gam = gv(gv_mixt_cp,ijkc0*indcp)/(gv(gv_mixt_cp,ijkc0*indcp)-rgas)
279  ggm1 = gam/(gam-1._rfreal)
280  rr = cv(cv_mixt_dens,ijkc0) - deltr(1)
281  ur = dv(dv_mixt_uvel,ijkc0) - deltr(2)
282  vr = dv(dv_mixt_vvel,ijkc0) - deltr(3)
283  wr = dv(dv_mixt_wvel,ijkc0) - deltr(4)
284  pr = dv(dv_mixt_pres,ijkc0) - deltr(5)
285  hr = ggm1*pr/rr + 0.5_rfreal*(ur*ur+vr*vr+wr*wr)
286  qsrr = (ur*sj(xcoord,ijkn)+vr*sj(ycoord,ijkn)+ &
287  wr*sj(zcoord,ijkn)-sjvel(ijkn*indsvel))*rr
288 
289 ! ----- fluxes
290 
291  pav = 0.5_rfreal*(pl+pr)
292  fc(1) = 0.5_rfreal*(qsrl +qsrr )
293  fc(2) = 0.5_rfreal*(qsrl*ul+qsrr*ur) + sj(xcoord,ijkn)*pav
294  fc(3) = 0.5_rfreal*(qsrl*vl+qsrr*vr) + sj(ycoord,ijkn)*pav
295  fc(4) = 0.5_rfreal*(qsrl*wl+qsrr*wr) + sj(zcoord,ijkn)*pav
296  fc(5) = 0.5_rfreal*(qsrl*hl+qsrr*hr) + sjvel(ijkn*indsvel)*pav
297 
298  rhs(cv_mixt_dens,ijkc0) = rhs(cv_mixt_dens,ijkc0) + fc(1)
299  rhs(cv_mixt_xmom,ijkc0) = rhs(cv_mixt_xmom,ijkc0) + fc(2)
300  rhs(cv_mixt_ymom,ijkc0) = rhs(cv_mixt_ymom,ijkc0) + fc(3)
301  rhs(cv_mixt_zmom,ijkc0) = rhs(cv_mixt_zmom,ijkc0) + fc(4)
302  rhs(cv_mixt_ener,ijkc0) = rhs(cv_mixt_ener,ijkc0) + fc(5)
303 
304  rhs(cv_mixt_dens,ijkc1) = rhs(cv_mixt_dens,ijkc1) - fc(1)
305  rhs(cv_mixt_xmom,ijkc1) = rhs(cv_mixt_xmom,ijkc1) - fc(2)
306  rhs(cv_mixt_ymom,ijkc1) = rhs(cv_mixt_ymom,ijkc1) - fc(3)
307  rhs(cv_mixt_zmom,ijkc1) = rhs(cv_mixt_zmom,ijkc1) - fc(4)
308  rhs(cv_mixt_ener,ijkc1) = rhs(cv_mixt_ener,ijkc1) - fc(5)
309  ENDDO ! i
310  ENDDO ! j
311  ENDDO ! k
312 
313 ! flux in k-direction (except through boundary) -------------------------------
314 
315  DO k=kpcbeg+1,kpcend
316  DO j=jpcbeg,jpcend
317  DO i=ipcbeg,ipcend
318  ijkc0 = indijk(i,j,k ,icoff,ijcoff)
319  ijkc1 = indijk(i,j,k-1,icoff,ijcoff)
320  ijkc2 = indijk(i,j,k-2,icoff,ijcoff)
321  ijkc3 = indijk(i,j,k+1,icoff,ijcoff)
322  ijkn = indijk(i,j,k ,inoff,ijnoff)
323 
324  vola = (0.5_rfreal*(vol(ijkc0)+vol(ijkc1)))**1.5_rfreal
325 
326 ! ----- limited extrapolations
327 
328  dvar(1) = cv(cv_mixt_dens,ijkc0) - cv(cv_mixt_dens,ijkc1)
329  dvar(2) = dv(dv_mixt_uvel,ijkc0) - dv(dv_mixt_uvel,ijkc1)
330  dvar(3) = dv(dv_mixt_vvel,ijkc0) - dv(dv_mixt_vvel,ijkc1)
331  dvar(4) = dv(dv_mixt_wvel,ijkc0) - dv(dv_mixt_wvel,ijkc1)
332  dvar(5) = dv(dv_mixt_pres,ijkc0) - dv(dv_mixt_pres,ijkc1)
333 
334  dvarm(1) = cv(cv_mixt_dens,ijkc1) - cv(cv_mixt_dens,ijkc2)
335  dvarm(2) = dv(dv_mixt_uvel,ijkc1) - dv(dv_mixt_uvel,ijkc2)
336  dvarm(3) = dv(dv_mixt_vvel,ijkc1) - dv(dv_mixt_vvel,ijkc2)
337  dvarm(4) = dv(dv_mixt_wvel,ijkc1) - dv(dv_mixt_wvel,ijkc2)
338  dvarm(5) = dv(dv_mixt_pres,ijkc1) - dv(dv_mixt_pres,ijkc2)
339 
340  dvarp(1) = cv(cv_mixt_dens,ijkc3) - cv(cv_mixt_dens,ijkc0)
341  dvarp(2) = dv(dv_mixt_uvel,ijkc3) - dv(dv_mixt_uvel,ijkc0)
342  dvarp(3) = dv(dv_mixt_vvel,ijkc3) - dv(dv_mixt_vvel,ijkc0)
343  dvarp(4) = dv(dv_mixt_wvel,ijkc3) - dv(dv_mixt_wvel,ijkc0)
344  dvarp(5) = dv(dv_mixt_pres,ijkc3) - dv(dv_mixt_pres,ijkc0)
345 
346  eps2n = eps2(1)*vola
347  deltl(1) = 0.5_rfreal*muscl3(dvar(1) ,dvarm(1),eps2n)
348  deltr(1) = 0.5_rfreal*muscl3(dvarp(1),dvar(1) ,eps2n)
349  eps2n = eps2(2)*vola
350  deltl(2) = 0.5_rfreal*muscl3(dvar(2) ,dvarm(2),eps2n)
351  deltr(2) = 0.5_rfreal*muscl3(dvarp(2),dvar(2) ,eps2n)
352  deltl(3) = 0.5_rfreal*muscl3(dvar(3) ,dvarm(3),eps2n)
353  deltr(3) = 0.5_rfreal*muscl3(dvarp(3),dvar(3) ,eps2n)
354  deltl(4) = 0.5_rfreal*muscl3(dvar(4) ,dvarm(4),eps2n)
355  deltr(4) = 0.5_rfreal*muscl3(dvarp(4),dvar(4) ,eps2n)
356  eps2n = eps2(3)*vola
357  deltl(5) = 0.5_rfreal*muscl3(dvar(5) ,dvarm(5),eps2n)
358  deltr(5) = 0.5_rfreal*muscl3(dvarp(5),dvar(5) ,eps2n)
359 
360 ! ----- left and right states
361 
362  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc1*indmol)
363  gam = gv(gv_mixt_cp,ijkc1*indcp)/(gv(gv_mixt_cp,ijkc1*indcp)-rgas)
364  ggm1 = gam/(gam-1._rfreal)
365  rl = cv(cv_mixt_dens,ijkc1) + deltl(1)
366  ul = dv(dv_mixt_uvel,ijkc1) + deltl(2)
367  vl = dv(dv_mixt_vvel,ijkc1) + deltl(3)
368  wl = dv(dv_mixt_wvel,ijkc1) + deltl(4)
369  pl = dv(dv_mixt_pres,ijkc1) + deltl(5)
370  hl = ggm1*pl/rl + 0.5_rfreal*(ul*ul+vl*vl+wl*wl)
371  qsrl = (ul*sk(xcoord,ijkn)+vl*sk(ycoord,ijkn)+ &
372  wl*sk(zcoord,ijkn)-skvel(ijkn*indsvel))*rl
373 
374  rgas = 8314.3_rfreal/gv(gv_mixt_mol,ijkc0*indmol)
375  gam = gv(gv_mixt_cp,ijkc0*indcp)/(gv(gv_mixt_cp,ijkc0*indcp)-rgas)
376  ggm1 = gam/(gam-1._rfreal)
377  rr = cv(cv_mixt_dens,ijkc0) - deltr(1)
378  ur = dv(dv_mixt_uvel,ijkc0) - deltr(2)
379  vr = dv(dv_mixt_vvel,ijkc0) - deltr(3)
380  wr = dv(dv_mixt_wvel,ijkc0) - deltr(4)
381  pr = dv(dv_mixt_pres,ijkc0) - deltr(5)
382  hr = ggm1*pr/rr + 0.5_rfreal*(ur*ur+vr*vr+wr*wr)
383  qsrr = (ur*sk(xcoord,ijkn)+vr*sk(ycoord,ijkn)+ &
384  wr*sk(zcoord,ijkn)-skvel(ijkn*indsvel))*rr
385 
386 ! ----- fluxes
387 
388  pav = 0.5_rfreal*(pl+pr)
389  fc(1) = 0.5_rfreal*(qsrl +qsrr )
390  fc(2) = 0.5_rfreal*(qsrl*ul+qsrr*ur) + sk(xcoord,ijkn)*pav
391  fc(3) = 0.5_rfreal*(qsrl*vl+qsrr*vr) + sk(ycoord,ijkn)*pav
392  fc(4) = 0.5_rfreal*(qsrl*wl+qsrr*wr) + sk(zcoord,ijkn)*pav
393  fc(5) = 0.5_rfreal*(qsrl*hl+qsrr*hr) + skvel(ijkn*indsvel)*pav
394 
395  rhs(cv_mixt_dens,ijkc0) = rhs(cv_mixt_dens,ijkc0) + fc(1)
396  rhs(cv_mixt_xmom,ijkc0) = rhs(cv_mixt_xmom,ijkc0) + fc(2)
397  rhs(cv_mixt_ymom,ijkc0) = rhs(cv_mixt_ymom,ijkc0) + fc(3)
398  rhs(cv_mixt_zmom,ijkc0) = rhs(cv_mixt_zmom,ijkc0) + fc(4)
399  rhs(cv_mixt_ener,ijkc0) = rhs(cv_mixt_ener,ijkc0) + fc(5)
400 
401  rhs(cv_mixt_dens,ijkc1) = rhs(cv_mixt_dens,ijkc1) - fc(1)
402  rhs(cv_mixt_xmom,ijkc1) = rhs(cv_mixt_xmom,ijkc1) - fc(2)
403  rhs(cv_mixt_ymom,ijkc1) = rhs(cv_mixt_ymom,ijkc1) - fc(3)
404  rhs(cv_mixt_zmom,ijkc1) = rhs(cv_mixt_zmom,ijkc1) - fc(4)
405  rhs(cv_mixt_ener,ijkc1) = rhs(cv_mixt_ener,ijkc1) - fc(5)
406  ENDDO ! i
407  ENDDO ! j
408  ENDDO ! k
409 
410 ! fluxes through boundaries ---------------------------------------------------
411 
412  DO ipatch=1,region%nPatches
413  CALL rflo_roefluxpatch( region,region%levels(ilev)%patches(ipatch) )
414  ENDDO
415 
416 ! finalize --------------------------------------------------------------------
417 
418  CALL deregisterfunction( region%global )
419 
420 END SUBROUTINE rflo_roefluxsecond
421 
422 !******************************************************************************
423 !
424 ! RCS Revision history:
425 !
426 ! $Log: RFLO_RoeFluxSecond.F90,v $
427 ! Revision 1.3 2008/12/06 08:44:28 mtcampbe
428 ! Updated license.
429 !
430 ! Revision 1.2 2008/11/19 22:17:38 mtcampbe
431 ! Added Illinois Open Source License/Copyright
432 !
433 ! Revision 1.1 2004/11/29 20:51:40 wasistho
434 ! lower to upper case
435 !
436 ! Revision 1.6 2003/11/20 16:40:40 mdbrandy
437 ! Backing out RocfluidMP changes from 11-17-03
438 !
439 ! Revision 1.2 2003/10/01 23:52:10 jblazek
440 ! Corrected bug in moving noslip wall BC and grid speeds.
441 !
442 ! Revision 1.1 2003/05/29 17:28:43 jblazek
443 ! Implemented Roe scheme.
444 !
445 !******************************************************************************
446 
447 
448 
449 
450 
451 
452 
k
j indices k indices k
Definition: Indexing.h:6
rhs
NT rhs
Definition: rational_rotation.h:86
kpcbeg
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE kpcbeg
Definition: RFLO_GetDimensPhys.h:27
moderror::registerfunction
subroutine registerfunction(global, funName, fileName)
Definition: ModError.F90:449
jpcbeg
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE jpcbeg
Definition: RFLO_GetDimensPhys.h:27
ipcend
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE ipcend
Definition: RFLO_GetDimensPhys.h:27
rflo_getnodeoffset
subroutine rflo_getnodeoffset(region, iLev, iNodeOffset, ijNodeOffset)
Definition: RFLO_GetDimensions.F90:267
ipcbeg
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE ipcbeg
Definition: RFLO_GetDimensPhys.h:27
i
blockLoc i
Definition: read.cpp:79
rflo_getcelloffset
subroutine rflo_getcelloffset(region, iLev, iCellOffset, ijCellOffset)
Definition: RFLO_GetDimensions.F90:235
j
j indices j
Definition: Indexing.h:6
jpcend
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE jpcend
Definition: RFLO_GetDimensPhys.h:27
rflo_roefluxpatch
subroutine rflo_roefluxpatch(region, patch)
Definition: RFLO_RoeFluxPatch.F90:45
moderror::deregisterfunction
subroutine deregisterfunction(global)
Definition: ModError.F90:469
rflo_roefluxsecond
subroutine rflo_roefluxsecond(region)
Definition: RFLO_RoeFluxSecond.F90:43
rflo_getdimensphys
subroutine rflo_getdimensphys(region, iLev, ipcbeg, ipcend, jpcbeg, jpcend, kpcbeg, kpcend)
Definition: RFLO_GetDimensions.F90:109