Flow Square
Flow Square
1. Non-reactive flows
u i 0,
t x i
(1)
u j uiu j
t
xi
ij
g j ,
x j xi
(2)
u j
t
ij
u i u j
x i
x i
g j ,
(4)
(5)
p *
u j
.
t
x j
x i2
u
u j
u k
ij i
x k
x j x i
2 p
2
3
ij
u j
x j
d
t
1
. ,
t
(6)
u j u j
p *
t.
x j
(7)
Yi
u j Yi
t
x j
x j
Di Yi
x j
i ,
(8)
T
T T
u jT
,
t
x j
x j c p x j c p
(9)
F O P,
cYF / O 1
cT
(10)
YF / O
Y
, cYP P , and
YF / O, u
YP , b
T Tu ,
Tb Tu
(11)
c
c
u j c
D
c , (12)
t
x j
x j
x j
c k exp
n
T (1 c) c ,
(13)
sY F
u j sY F
t
x j
x j
O
s
P
1 s
T
Q
D sY F
x j
s F ,
YO
u j YO
t
x j
x j
D YO
x j
O .
(20)
(21)
(14)
F ,O
u j F ,O
x j
x j
D F ,O
x j
(22)
F ,O sY F YO ,
F , P YF
O , P
(15)
YP
,
1 s
(16)
Y
YO P ,
1 s
F ,T YF
c pT
F ,O YO
sc p T
(23)
(17)
OX
,
FU OX
(18)
(19)
As shown in Fig. 2b, the flame location in nonpremixed systems is often assumed as the
location where st . Here, st is the
4 / 10
st
OX
.
FU OX
(24)
u j
D
.
t
x j
x j
x j
(25)
p RT .
For st
YF 0 ,
YO YO0 1
st
YP (1 s)YF0 ,
T (Tb TO0 )
TO0 .
st
For st
(26)
YF YF0
st
,
1 st
n1 n 1 (1 ) n ,
YO 0 ,
YP
(1 s)YF0
(27)
(28)
1 st
Tb TF0
T 0 T
F st b ,
st 1
st 1
5 / 10
u
v
2
u
u p
uv
.
q , F
,G 2
v
uv
v p
E
( E p)u
( E p )v
(33)
3. Subsonic/supersonic flows
In high speed flows, the viscous term in the
momentum transport equation (Eq. 2) may be
negligible:
p
.
u j
u i u j
t
xi
x j
(29)
They are solved (also with Eq 31) for fullycompressible subsonic/ supersonic (inviscid)
flows.
( E p)u i 0 ,
t x i
4. Numerical methods
In Flowsquare ver 4, there are a variety of
choices of numerical methods. Each method
has pros and cons, and they need to be
considered for the choice. In this section,
these methods are summarised.
4.1. Differentiations
(30)
(31)
q F G
0,
t x y
x ( f j 1 f j 1 ) /( 2xi ) ,
i j
(34)
x
i
( f 1 f 2 ) / x i ,
(35)
x
i
( f N f N 1 ) / x i .
(36)
x ( f j 2 8 f j 1 8 f j 1 fi 2 ) /(12xi )
i j
(36)
where
f 25 f1 48 f 2 36 f3 16 f 4 3 f5
x
12xi
i 1
(37)
6 / 10
(38)
x
12xi
i 2
f
3 f 10 f N 1 18 f N 2 6 f N 3 f N 4
N
x
12xi
i N 1
(39)
x
i N
25 f N 48 f N 1 36 f N 2 16 f N 3 3 f N 4
12xi
(40)
f n t
f n
,
t
(44)
f n f n f n ,
f n
(45)
f
15 5
t f n n ,
16 9
t
(46)
f n f n f n ,
f n
(47)
f
8
153
t
f n n ,
15 128
t
(48)
f n1 f n f n .
(49)
f n
t .
t
(41)
2 f ,
i 1, j 2 i , j i 1, j
x 2
i , j 1 2 i , j i , j 1
f n t
,
t 2
f n 1 f n
f n
t .
t
(50)
y 2
(51)
f i, j
(42)
(43)
7 / 10
in, j 1 (1 ) in, j
in1, j in11,j in, j 1 in, j 11
f ij
2
2
x
y
x 2 y 2
,
2(x 2 y 2 )
(53)
The
relaxation
parameter fil
ranges
recommended.
4.5. Boundary conditions
Various boundary conditions can be used in
Flowsquare. Boundary conditions are roughly
classified into two types: Dirichlet boundary
condition and Neumann boundary condition.
Dirichlet BC is also called fixed BC, because it
specifies the values that a solution takes on
the boundary of the computational domain. In
contrast, Neumann BC specifies values that
the derivative of a solution is to take on the
boundary of the domain. The value is
determined to make the derivative zero
( f 1 f 2 and f N f N 1 ). In Flowsquare, these
8 / 10
x ( f 2 f N ) /( 2xi ) ,
i 1
(54)
x
i
(55)
( f 1 f N 1 ) /(2x i ) .
i 1
u p u mag sin 2
n wave
Nx 1
j 1
cos 2
n wave
Nx 1
i 1
u mag cos 2
n wave
Nx 1
(56)
j 1
sin 2
n wave
Nx 1
i 1
v p u mag cos 2
n wave
Nx 1
j 1
sin 2
n wave
Nx 1
i 1
u mag sin 2
n wave
Nx 1
j 1
cos 2
n wave
Nx 1
(57)
i 1 n wave
u mag sin
N x 1 2 n 1
j 1 n wave
cos
M
N x 1 2 n 1
,(58)
up
i 1 n wave
n 0
u mag cos
N x 1 2 n 1
j 1 n wave
sin
N x 1 2 n 1
9 / 10
u mag
M
vp
n 0
u mag
i 1 n wave
cos
N x 1 2 n 1
j 1 n wave
sin
N x 1 2 n 1
,(59)
i 1 n wave
sin
N x 1 2 n 1
j 1 n wave
cos
N x 1 2 n 1
u rnd
M
up
n 0
u rnd
i 1 n wave
sin
N x 1 2 n 1
j 1 n wave
cos
N x 1 2 n 1
,(60)
i 1 n wave
cos
N x 1 2 n 1
j 1 n wave
sin
N x 1 2 n 1
i 1 n wave
v rnd cos
N x 1 2 n 1
j 1 n wave
sin
M
N x 1 2 n 1
,(61)
up
i 1 n wave
n 0
v rnd sin
N x 1 2 n 1
j 1 n wave
cos
N x 1 2 n 1
and
M log 2 nwave .
The
10 / 10