Hi R,

I tried out modifying the thermo electric example, but I was getting

convdiff: left over: ['verbose', '__builtins__', 'n_step', '__file__', '__doc__', '__name__', 't1', 'sys', 'data_dir', 't0', '__package__', '_filename', 'main', 'os', 'cwd']
convdiff: reading mesh (/usr/local/lib/python2.7/dist-packages/sfepy/meshes/2d/rectangle_fine_quad.mesh)...
convdiff: ...done in 0.02 s
convdiff: creating regions...
convdiff:     Wall
convdiff:     Top
convdiff:     Surface
convdiff:     Driven
convdiff:     Entry
convdiff:     Omega
convdiff: ...done in 0.03 s
convdiff: using solvers:
                ts: ts
               nls: newton
                ls: ls
convdiff: equation "balance":
convdiff: + dw_div_grad.5.Omega(m.viscosity, v, u)
       + dw_convect.5.Omega(v, u)
       - dw_stokes.5.Omega(v, p) = 0
convdiff: equation "incompressibility":
convdiff: dw_stokes.5.Omega(u, q) = 0
convdiff: setting up dof connectivities...
convdiff: ...done in 0.00 s
convdiff: updating variables...
convdiff: ...done
convdiff: matrix shape: (44949, 44949)
convdiff: assembling matrix graph...
convdiff: ...done in 0.15 s
convdiff: matrix structural nonzeros: 1747752 (8.65e-04% fill)
convdiff: updating materials...
convdiff:     m
convdiff: ...done in 0.01 s
convdiff: nls: iter: 0, residual: 2.002082e-02 (rel: 1.000000e+00)
convect_build_vtg(): ERR_Switch
Traceback (most recent call last):
File "examples/phase_change/convective_diffusive/convective_diffusive.py", line 139, in <module>
    main()
File "examples/phase_change/convective_diffusive/convective_diffusive.py", line 125, in main
    state_el = problem.solve()
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/problemDef.py", line 933, in solve
    vec = solvers.nls(vec0)
File "/usr/local/lib/python2.7/dist-packages/sfepy/solvers/nls.py", line 345, in __call__
    mtx_a = fun_grad(vec_x)
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/evaluate.py", line 66, in eval_tangent_matrix
    mtx = pb.equations.eval_tangent_matrices(vec, mtx)
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 640, in eval_tangent_matrices
    self.evaluate(mode='weak', dw_mode='matrix', asm_obj=tangent_matrix)
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 526, in evaluate
    asm_obj=asm_obj)
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/equations.py", line 766, in evaluate
    ret_status=True)
File "/usr/local/lib/python2.7/dist-packages/sfepy/terms/terms.py", line 1474, in evaluate
    vals.append(self.sign * val)
RuntimeError: ccore error (see above)

Then I modified it to solve without any timestep and still it is throwing errors:

convdiff: left over: ['os', '__builtins__', '__doc__', '__name__', '__package__', 'verbose', 'sys', 'data_dir', 'cwd', '_filename', 'main', '__file__']
convdiff: reading mesh (/usr/local/lib/python2.7/dist-packages/sfepy/meshes/2d/rectangle_fine_quad.mesh)...
convdiff: ...done in 0.02 s
convdiff: creating regions...
convdiff:     Wall
convdiff:     Top
convdiff:     Surface
convdiff:     Driven
convdiff:     Entry
convdiff:     Omega
convdiff: ...done in 0.03 s
convdiff: using solvers:
                ts: no ts
               nls: newton
                ls: ls
convdiff: equation "balance":
convdiff: + dw_div_grad.5.Omega(m.viscosity, v, u)
       + dw_convect.5.Omega(v, u)
       - dw_stokes.5.Omega(v, p) = 0
convdiff: equation "incompressibility":
convdiff: dw_stokes.5.Omega(u, q) = 0
convdiff: setting up dof connectivities...
convdiff: ...done in 0.00 s
convdiff: updating materials...
convdiff:     m
convdiff: ...done in 0.01 s
Traceback (most recent call last):
File "examples/phase_change/convective_diffusive/convective_diffusive.py", line 126, in <module>
    main()
File "examples/phase_change/convective_diffusive/convective_diffusive.py", line 112, in main
    flow = problem.solve()
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/problemDef.py", line 929, in solve
    state0.apply_ebc(force_values=force_values)
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/state.py", line 107, in apply_ebc
    self.variables.apply_ebc(self.vec, force_values=force_values)
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/variables.py", line 308, in apply_ebc
    var.apply_ebc(vec, self.di.indx[var.name].start, force_values)
File "/usr/local/lib/python2.7/dist-packages/sfepy/fem/variables.py", line 1780, in apply_ebc
    eq_map = self.eq_map
AttributeError: 'FieldVariable' object has no attribute 'eq_map'

I am attaching the code.
I have spent lot of hours trying to find my error, but to no avail.
Kindly help.

Regards.

On Wednesday, 3 July 2013 17:25:42 UTC+5:30, Ankit Mahato wrote:

On Wednesday, 3 July 2013 17:19:57 UTC+5:30, Robert Cimrman wrote:
On 07/03/2013 01:38 PM, Ankit Mahato wrote:
>
>
> On Wednesday, 3 July 2013 13:14:37 UTC+5:30, Robert Cimrman wrote:
>>
>> On 07/03/2013 01:51 AM, Ankit Mahato wrote:
>>> R,
>>>
>>> Attached is the Temperature distribution along the length with Peclet
>>> number graph which you wanted to have a look at.
>>> I used ParaView to generate it.
>>> It is validated by the Fig 6.2 of Computer Simulation of Flow and Heat
>>> Transfer, P S Ghoshdastidar, Tata McGraw-Hill.
>>
>> Good! Btw. are there some analytical relations that can be verified, for
>> example on some special simple geometries/boundary conditions? It would be
>> interesting to see the comparison of analytical/numerical curves in a
>> single
>> figure (I have no access to the book).
>>
>
> Yes this one can be verified with the analytical solution.
> I was about to ask you that the analytical vs simulation curve should be
> plotted for how many test cases?

As you wish (if it is ok :)) BTW. http://terri.toybox.ca/python-soc/ still does
not show your blog posts - is the subscription correct?

Yes R I saw it yesterday and I thought maybe it is not updated.
The subscription link points to http://ankitmahato.blogspot.in/search/label/Python
which is correct.
I think I should write to Terri

r.

>>
>> Similar figures could be also done in sfepy directly, check [1] - it is
>> mostly
>> useful when one tweaks and runs a simulation repeatedly, as the setup
>> needs
>> some effort...
>>
>
> Okie :)
>
>
>>
>> r.
>>
>> [1] http://sfepy.org/doc-devel/primer.html#probing
>>
>