Dear Nicolò,
The data has to be a vector field on the domain. The brick is
equivalent to the term (which can be added equivalently with
'add_linear_term')
'(V.Normal).Test_u'
For 'V' the data and 'u' the vector variable. So the dimension is 3xN
with N the number of degrees of freedom of a scalar field. If you use a
vector field, this should be just N, a priori.
Best regards,
Yves
Le 22/10/2019 à 11:48, Nicolò Cambiaso a écrit :
Sorry, I forgot to define you the region, that is a side of the cube,
gf_mesh_set( Mesh_M , 'region', RIGHT_BOUND, fb1);
fb1 = gf_mesh_get( Mesh_M , 'outer faces with direction', [ 1 0 0],
0.01);
Thank you very much,
Have a nice day
Nicolò Cambiaso
Il giorno mar 22 ott 2019 alle ore 11:43 Nicolò Cambiaso
<[email protected] <mailto:[email protected]>> ha scritto:
Dear GetFem users,
I'm trying to apply a normal boundary condition, using the command
gf_model_set(model M, 'add normal source term brick', mesh_im
mim,string varname, string dataname, int region),
using matlab interface.
and I previously defined the components as follows:
mim = gf_mesh_im( mf, IT);
IT = gf_integ('IM_TETRAHEDRON(1)')
mf = gf_mesh_fem(Mesh_M, 3);
Mesh_M is a cube, [0,1]x [0,1]x [0,1]
string varname is 'ut', and it was defined as: gf_model_set(md,
'add fem variable', 'ut', mf); (md is the model).
but the problem is the "string dataname". The dataname was
initialized by
gf_model_set(md, 'add initialized data', 'dataname', Data);
and previously Data was defined as a 3xN matrix, where 3 is the
dimension of the vectors (3D). Hoping right now that everything is
correct, what is N ??
Because it's not the number of degrees of freedom of mf, as i
thought at the beginning.
Thank you very much,
Have a nice day
Nicolò Cambiaso
--
Yves Renard ([email protected]) tel : (33) 04.72.43.87.08
INSA-Lyon
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard
---------
