# -- Poiseuille flow in a tube, 3D cylindrical solution, BC driven.
##############################################################################
# 4 element axisymmetric 3D flow in a tube, parabolic input profile.

<FIELDS>
	u v w p
</FIELDS>

<USER>
	u = 1-y*y
	v = 0
	w = 0
	p = 0
</USER>

<TOKENS>
	CYLINDRICAL = 1
	N_Z         = 8
	N_TIME      = 2
	N_P      = 9
	N_STEP      = 20
	D_T         = 0.01
	KINVIS      = 0.5
	TOL_REL     = 1e-15
</TOKENS>

<GROUPS NUMBER=4>
	1	v	value
	2	w	wall
	3	o	outflow
	4	a	axis
</GROUPS>

<BCS NUMBER=4>
	1	v	4
			<D>	u = 1.0-y*y	</D>
			<D>	v = 0		</D>
			<D>     w = 0		</D>
			<H>	p = 0		</H>
	2	w	4
			<D>	u = 0		</D>
			<D>	v = 0		</D>
			<D>	w = 0		</D>
			<H>	p = 0		</H>
	3	o	4
			<N>	u = 0		</N>
			<N>	v = 0		</N>
			<N>	w = 0		</N>
			<D>	p = 0		</D>
	4	a	4
			<A>	u = 0		</A>
			<A>	v = 0		</A>
			<A>	w = 0		</A>
			<A>	p = 0		</A>
</BCS>

<NODES NUMBER=9>
	1	0.0	0.0	0.0
	2	2.0	0.0	0.0
	3	4.0	0.0	0.0
	4	0.0	0.5	0.0
	5	2.0	0.5	0.0
	6	4.0	0.5	0.0
	7	0.0	1.0	0.0
	8	2.0	1.0	0.0
	9	4.0	1.0	0.0
</NODES>

<ELEMENTS NUMBER=4>
	1	<Q>	1 2 5 4		</Q>
	2	<Q>	2 3 6 5		</Q>
	3	<Q>	4 5 8 7		</Q>
	4	<Q>	5 6 9 8		</Q>
</ELEMENTS>

<SURFACES NUMBER=8>
	1	1	1	<B>	a	</B>
	2	2	1	<B>	a	</B>
	3	2	2	<B>	o	</B>
	4	4	2	<B>	o	</B>
	5	4	3	<B>	w	</B>
	6	3	3	<B>	w	</B>
	7	3	4	<B>	v	</B>
	8	1	4	<B>	v	</B>
</SURFACES>
