luchete80 / EPSol

0 stars 0 forks source link

Bad F derivative in 1 Radial Flow element example #85

Closed luchete80 closed 3 years ago

luchete80 commented 3 years ago

At FIRST GAUSS POINT F derivative is null. But for this example F is varying in both directions.

rg,sg -0.5773502691896258 -0.5773502691896258

UV [[ 0.04980974] [ 0.00435778] [ 0.09961947] [ 0.00871557] [ 0.09961947] [-0.00871557] [ 0.04980974] [-0.00435778]]

v [[ 0.08909343] [-0.00450025]]

dVxy [[-0.05 ] [ 0. ] [ 0.00625492] [ 0.07383136]]

for l in range(4): for m in range(16):
for n in range(2): temp[l,n]=temp[l,n]+BsigF[l,m,n]*UF[m] print("temp[l,n]",temp) F derivative temp[l,n] [[0. 0.] [0. 0.] [0. 0.] [0. 0.]]

luchete80 commented 3 years ago

at first gauss point: v [[ 0.08909343] [-0.00450025]] F [[0.85207259] [0. ] [0. ] [1.84529946]]