On 07/02/2013 01:58 PM, Ankit Mahato wrote:
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)
Did you rebuild the sources?
Regards,
Ankit
PS: I was not aware that we have to blog weekly. I will blog for week 1 and
2 asap.
It's either weekly or bi-weekly...
r.
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.
>>>