R,
I was getting an error while running the example.
./simple.py examples/navier_stokes/navier_stokes2d.py
sfepy: left over: ['data_dir', 'verbose', '_filename', '__builtins__', '__package__', '__doc__', '__name__', '__file__']
sfepy: reading mesh (/home/ankit/sfepy/meshes/2d/rectangle_fine_quad.mesh)...
sfepy: ...done in 0.02 s
sfepy: creating regions...
sfepy: Right
sfepy: Bottom
sfepy: Top
sfepy: Surface
sfepy: Omega
sfepy: Walls
sfepy: Driven
sfepy: ...done in 0.03 s
sfepy: equation "balance":
sfepy: + dw_div_grad.5.Omega(fluid.viscosity, v, u)
+ dw_convect.5.Omega(v, u)
- dw_stokes.5.Omega(v, p) = 0
sfepy: equation "incompressibility":
sfepy: dw_stokes.5.Omega(u, q) = 0
sfepy: setting up dof connectivities...
sfepy: ...done in 0.00 s
sfepy: using solvers:
ts: no ts
nls: newton
ls: ls
sfepy: updating variables...
sfepy: ...done
sfepy: matrix shape: (44557, 44557)
sfepy: assembling matrix graph...
sfepy: ...done in 0.13 s
sfepy: matrix structural nonzeros: 1728264 (8.71e-04% fill)
sfepy: updating materials...
sfepy: fluid
sfepy: ...done in 0.01 s
sfepy: nls: iter: 0, residual: 2.023986e-02 (rel: 1.000000e+00)
convect_build_vtg(): ERR_Switch
mem_free_mem(): error exit!
mem_free_mem(): error exit!
mem_free_mem(): error exit!
mem_free_mem(): error exit!
mem_free_mem(): error exit!
mem_free_mem(): error exit!
mem_free_mem(): error exit!
mem_free_mem(): error exit!
Traceback (most recent call last):
File "./simple.py", line 146, in <module>
main()
File "./simple.py", line 143, in main
app()
File "/home/ankit/sfepy/sfepy/applications/application.py", line 29, in call_basic
return self.call(**kwargs)
File "/home/ankit/sfepy/sfepy/applications/pde_solver_app.py", line 213, in call
nls_status=nls_status)
File "/home/ankit/sfepy/sfepy/solvers/ts_solvers.py", line 29, in __call__
state = problem.solve(state0=state0, nls_status=nls_status)
File "/home/ankit/sfepy/sfepy/fem/problemDef.py", line 933, in solve
vec = solvers.nls(vec0)
File "/home/ankit/sfepy/sfepy/solvers/nls.py", line 345, in __call__
mtx_a = fun_grad(vec_x)
File "/home/ankit/sfepy/sfepy/fem/evaluate.py", line 66, in eval_tangent_matrix
mtx = pb.equations.eval_tangent_matrices(vec, mtx)
File "/home/ankit/sfepy/sfepy/fem/equations.py", line 640, in eval_tangent_matrices
self.evaluate(mode='weak', dw_mode='matrix', asm_obj=tangent_matrix)
File "/home/ankit/sfepy/sfepy/fem/equations.py", line 526, in evaluate
asm_obj=asm_obj)
File "/home/ankit/sfepy/sfepy/fem/equations.py", line 766, in evaluate
ret_status=True)
File "/home/ankit/sfepy/sfepy/terms/terms.py", line 1473, in evaluate
diff_var, **kwargs)
File "/home/ankit/sfepy/sfepy/terms/terms.py", line 1311, in eval_real
status = self.call_function(out, fargs)
File "/home/ankit/sfepy/sfepy/terms/terms.py", line 1296, in call_function
raise ValueError('term evaluation failed! (%s)' % self.name)
ValueError: term evaluation failed! (dw_convect)
Regards,
Ankit
PS: I was not aware that we have to blog weekly. I will blog for week 1 and 2 asap.
On Tuesday, 2 July 2013 13:58:20 UTC+5:30, Robert Cimrman wrote:
Now it remains to implement a robust flow solver. Even this small example
shows, that the solution is not obtained easily - try decreasing the viscosity,
and/or increase the Dirichlet velocity - the solver would not converge.
Cheers,
r.
PS: As mentioned in Terri Oda's e-mail, you should blog about your work so far
ASAP!
On 07/01/2013 06:12 PM, Ankit Mahato wrote:
> awesome :)
>
> On Monday, 1 July 2013 15:05:16 UTC+5:30, Robert Cimrman wrote:
>>
>> Hi,
>>
>> I have removed the "3d only" restriction from the Navier Stokes and
>> related
>> terms. There is also a new example:
>> examples/navier_stokes/navier_stokes2d.py.
>>
>> r.
>>
>