@@ -40,10 +40,6 @@ def main(inflow: 'inflow velocity' = 10,
40
40
tθ = lambda f : theta * f + (1 - theta ) * t0 (f )
41
41
# 1st order FD
42
42
δt = lambda f : (f - t0 (f )) / function .Argument ('dt' , ())
43
- # extrapolation for pressure
44
- tp = lambda f : (1.5 * f - 0.5 * t0 (f ))
45
- # 2nd order FD
46
- tt = lambda f : tp (δt (f ))
47
43
48
44
# Nutils namespace
49
45
ns = function .Namespace ()
@@ -53,7 +49,7 @@ def main(inflow: 'inflow velocity' = 10,
53
49
ns .x0 = geom # reference geometry
54
50
ns .dbasis = domain .basis ('std' , degree = 1 ).vector (2 )
55
51
ns .d_i = 'dbasis_ni ?meshdofs_n'
56
- ns .umesh = tt (ns .d ) # mesh velocity
52
+ ns .umesh = δt (ns .d ) # mesh velocity
57
53
ns .x_i = 'x0_i + d_i' # moving geometry
58
54
ns .ubasis , ns .pbasis = function .chain ([domain .basis ('std' , degree = 2 ).vector (2 ), domain .basis ('std' , degree = 1 ), ])
59
55
ns .F_i = 'ubasis_ni ?F_n' # stress field
@@ -118,7 +114,7 @@ def main(inflow: 'inflow velocity' = 10,
118
114
timestep = 0
119
115
arguments = dict (lhs = numpy .zeros (len (ns .ubasis )), F = numpy .zeros (len (ns .ubasis )), meshdofs = numpy .zeros (len (ns .dbasis )), dt = timestepsize )
120
116
121
- nhist = 4 # history length required by the formulation
117
+ nhist = 2 # history length required by the formulation
122
118
arguments .update ((k + '0' * i , v ) for k , v in tuple (arguments .items ()) for i in range (nhist ))
123
119
124
120
while interface .is_coupling_ongoing ():
0 commit comments