Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.
Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.
discretize block into smaller blocks
Posted 12.06.2012, 14:25 GMT-4 Geometry Version 4.2a 6 Replies
Please login with a confirmed email address before reporting spam
I have a block domai (500x5000x500), which I want to discretize it into smaller blocks of 25x25x25.
I tried to create 19 layers of 25m width (left,front,top), but after 20 min of running, I gave up. Is there a faster way to do this?
P.S.
The commands I used to do this
model.geom('geom1').feature.create('blk1', 'Block');
model.geom('geom1').feature('blk1').set('pos', {'0' '0' '-250'});
model.geom('geom1').feature('blk1').set('base', 'center');
model.geom('geom1').feature('blk1').set('size', {'500' '500' '500'});
thickness=num2str(25);
for i=1:20-1
model.geom('geom1').feature('blk1').setIndex('layer', thickness, i-1);
end
Please login with a confirmed email address before reporting spam
not sure why but that makes some 20*200*20 = 80000 blocks, that is almost a mesh ;)
I can believe it takes somet time, still it should be possible, but you have to deal with 6 time as many boundaries (in asembly something less in Union mode). You should at least turn off the automatic node update and plot update in the preferences ifyou navigate with such a model. With our largest model and some 1200 COMSOL generated domains we needed 5 minutes to travel up and down in the model tree
Are you sure you cannot use a mapped mesh instead and distrbute your variablesin the BCs and the physics as interpoaltion functions over the spatial variables ?
--
Good luck
Ivar
Please login with a confirmed email address before reporting spam
Thanks again.
The main concerain I am doing this, is to be able to refer to a spefic block easily, like block number 111.
With mapped mesh, I haven't found an easy way to assign attributes to blocks, especially when I am using parametric solver.
Any ideas/suggestions would be nice.
Please login with a confirmed email address before reporting spam
In COMSOL you do not access the mesh, you access the domains and the boundaries, based on their coordinates (x,y,z,t). Most physics and BC entries are expecting fields V(x,y,z,t) and not just constant or Parameter values. All physics apply to the element dx*dy*dz which might map to a mesh element or to several. The entity is not the mesh, that is "just" a discretization of the domain and boundaries for the maths of the numerical processing. And a mesh element is further discretized y higher order functions (instead of having 3-4-6-9 ... node elemets)
This is a fundamental difference between COMSOL and most older FEM programmes.
And a source of confusion and frustration for those of us having dealt too long with older programmes. This is so until you catch the new approach, then you ask yourself how you could have been working in the old way for so long ...
--
Good luck
Ivar
Please login with a confirmed email address before reporting spam
I need your advise again.
In the attached model, if I used a mapped mesh, (mesh1), I fail to find a solution (convergence of CG).
If I use tetrahedra (mesh2), I can find a solution.
Why is that?
Thanks in advance.
Attachments:
Please login with a confirmed email address before reporting spam
first of all onthe principle, you can leave the mapped sweep mesh for the cube, only convert the 6 BOUNDARIES of the cube to tri, and then finish with a Thets
Apart from that I do not see why the mesh should matter, but i might be missing somethingt too ;)
But why do you not trust thet mesh, if it solves ?
It could be that your source term is to abrupt and that you should "smooth" it with some rectangle functions
--
Good luck
Ivar
Please login with a confirmed email address before reporting spam
I guess it is some kind of a bug.
Thank you for your help.