I tried to do this, but it was still wrong. I draw the interface of the
fluid-solid coupling without imposing any physical labels. It still doesn't
pass the Gridin
On Friday, June 14, 2024 at 10:28:41 AM UTC+8 hkch...@gmail.com wrote:
> You should not add Physical Surface("inner wall surface", 776) of
> interface in gmsh,
>
> This interface can be determined by the step-46, this is the drawback of
> dealii as I know.
>
>
>
> 王艺康 <yikangw...@gmail.com> 于2024年6月13日周四 14:27写道：
>
>> Hello，
>> Sure. I just use GridIN to read the grid.
>>
>> using namespace dealii;
>> void run()
>> {
>>
>> Triangulation<3> triangulation;
>> GridIn<3> grid_in;
>> grid_in.attach_triangulation(triangulation);
>>
>> {
>> std::string filename = "heart_have_thickness(1).msh";
>> std::ifstream file(filename);
>> Assert(file, ExcFileNotOpen(filename.c_str()));
>> grid_in.read_msh(file);
>> }
>>
>> }
>>
>> int main()
>> {
>>
>> run();
>>
>> }
>>
>> On Wednesday, June 12, 2024 at 10:46:04 PM UTC+8 blais...@gmail.com
>> wrote:
>>
>>> Can you provide the GMSH code and a small executable that tries to
>>> achieve what you are doing? I can take a look and figure out what is the
>>> issue :)
>>>
>>>
>>> On Wednesday, June 12, 2024 at 3:08:19 a.m. UTC-4 yikangw...@gmail.com
>>> wrote:
>>>
>>>> Hollo,
>>>> I tried your suggestion but dealii still reports an error. I don't
>>>> think it has anything to do with the grid type, it also errors when I only
>>>> use Gridin.
>>>>
>>>> Wang
>>>>
>>>> On Friday, June 7, 2024 at 11:32:16 PM UTC+8 blais...@gmail.com wrote:
>>>>
>>>>> Loading a tet mesh works if you set-up your Mapping/FE accordingly,
>>>>> this is a feature we use extensively.
>>>>>
>>>>>
>>>>> On Wednesday, June 5, 2024 at 1:28:20 p.m. UTC+1 bruno.t...@gmail.com
>>>>> wrote:
>>>>>
>>>>>> Wang,
>>>>>>
>>>>>> Are you loading a hex mesh or a tet mesh? I don't know if loading a
>>>>>> tet mesh works.
>>>>>>
>>>>>> Best,
>>>>>>
>>>>>> Bruno
>>>>>>
>>>>>> Le mar. 4 juin 2024 à 21:12, 王艺康 <yikangw...@gmail.com> a écrit :
>>>>>>
>>>>>>> Bruno,
>>>>>>> In fact that's what I did. The error is the same for regenerated
>>>>>>> meshes entered into dealii. I redrew a simple geometry again without
>>>>>>> defining any physical information still the same error. I don't
>>>>>>> understand
>>>>>>> why this is wrong, if I generate only the fluid-solid coupling
>>>>>>> interface to
>>>>>>> the external contour (solid domain) mesh there is no problem, but the
>>>>>>> same
>>>>>>> error occurs when I generate only the internal (fluid domain).
>>>>>>>
>>>>>>> Wang
>>>>>>>
>>>>>>>
>>>>>>> [image: 捕获1.PNG][image: 捕获2.PNG]
>>>>>>> [image: 捕获3.PNG]
>>>>>>> [image: 捕获4.PNG]
>>>>>>>
>>>>>>>
>>>>>>> On Tuesday, June 4, 2024 at 10:28:21 PM UTC+8 bruno.t...@gmail.com
>>>>>>> wrote:
>>>>>>>
>>>>>>>> Wang,
>>>>>>>>
>>>>>>>> Do you have a geo file that describes the geometry? Can you remove
>>>>>>>> the "Physical Surface" lines in the geo file and regenerate the mesh.
>>>>>>>>
>>>>>>>> Best,
>>>>>>>>
>>>>>>>> Bruno
>>>>>>>>
>>>>>>>> Le mar. 4 juin 2024 à 10:13, 王艺康 <yikangw...@gmail.com> a écrit :
>>>>>>>>
>>>>>>>>> Hello,
>>>>>>>>> Thank you very much for your reply. That's what I thought at
>>>>>>>>> first, but I was wrong, when I removed the corresponding physical
>>>>>>>>> group it
>>>>>>>>> was still the same error. I think gmsh automatically sets the
>>>>>>>>> coupling
>>>>>>>>> interface to solid when drawing the model body mesh, but I don't know
>>>>>>>>> how
>>>>>>>>> to change it. I didn't learn gmsh deep enough.
>>>>>>>>>
>>>>>>>>> Wang
>>>>>>>>>
>>>>>>>>> On Tuesday, June 4, 2024 at 8:28:07 PM UTC+8 bruno.t...@gmail.com
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>> Hello,
>>>>>>>>>>
>>>>>>>>>> Deal.II associates the physical group of a face with a boundary
>>>>>>>>>> id. You need to remove the physical group that are associated with
>>>>>>>>>> interior
>>>>>>>>>> faces.
>>>>>>>>>>
>>>>>>>>>> Best,
>>>>>>>>>>
>>>>>>>>>> Bruno
>>>>>>>>>>
>>>>>>>>>> On Tuesday, June 4, 2024 at 3:28:05 AM UTC-4 yikangw...@gmail.com
>>>>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>>> I drew the model in gmsh but importing to dealii is giving me an
>>>>>>>>>>> error.This step error occurs when importing a model.
>>>>>>>>>>>
>>>>>>>>>>> The violated condition was:
>>>>>>>>>>> boundary_id != numbers::internal_face_boundary_id
>>>>>>>>>>> Additional information:
>>>>>>>>>>> The input arguments for creating a triangulation specified a
>>>>>>>>>>> boundary
>>>>>>>>>>> id for an internal face. This is not allowed.
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> Here only a part of the mesh is shown (the outer contour also
>>>>>>>>>>> generates the mesh),the error should occur at the interface of the
>>>>>>>>>>> two
>>>>>>>>>>> bodies, but I don't know how to fix it. This model will be used for
>>>>>>>>>>> fluid-solid coupling, where the interior is the fluid domain and
>>>>>>>>>>> the
>>>>>>>>>>> contours are the solid domain.
>>>>>>>>>>> [image: 捕获.PNG]
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> --
>>>>>>>>> The deal.II project is located at http://www.dealii.org/
>>>>>>>>> For mailing list/forum options, see
>>>>>>>>> https://groups.google.com/d/forum/dealii?hl=en
>>>>>>>>> ---
>>>>>>>>> You received this message because you are subscribed to a topic in
>>>>>>>>> the Google Groups "deal.II User Group" group.
>>>>>>>>> To unsubscribe from this topic, visit
>>>>>>>>> https://groups.google.com/d/topic/dealii/WD4qeXgUf0E/unsubscribe.
>>>>>>>>> To unsubscribe from this group and all its topics, send an email
>>>>>>>>> to dealii+un...@googlegroups.com.
>>>>>>>>> To view this discussion on the web visit
>>>>>>>>> https://groups.google.com/d/msgid/dealii/6658e3cb-5f21-46a4-ab4f-08bbed3fa975n%40googlegroups.com
>>>>>>>>>
>>>>>>>>> <https://groups.google.com/d/msgid/dealii/6658e3cb-5f21-46a4-ab4f-08bbed3fa975n%40googlegroups.com?utm_medium=email&utm_source=footer>
>>>>>>>>> .
>>>>>>>>>
>>>>>>>> --
>>>>>>> The deal.II project is located at http://www.dealii.org/
>>>>>>> For mailing list/forum options, see
>>>>>>> https://groups.google.com/d/forum/dealii?hl=en
>>>>>>> ---
>>>>>>> You received this message because you are subscribed to a topic in
>>>>>>> the Google Groups "deal.II User Group" group.
>>>>>>> To unsubscribe from this topic, visit
>>>>>>> https://groups.google.com/d/topic/dealii/WD4qeXgUf0E/unsubscribe.
>>>>>>> To unsubscribe from this group and all its topics, send an email to
>>>>>>> dealii+un...@googlegroups.com.
>>>>>>>
>>>>>> To view this discussion on the web visit
>>>>>>> https://groups.google.com/d/msgid/dealii/12ce6fac-41c7-430e-935b-6a2326083850n%40googlegroups.com
>>>>>>>
>>>>>>> <https://groups.google.com/d/msgid/dealii/12ce6fac-41c7-430e-935b-6a2326083850n%40googlegroups.com?utm_medium=email&utm_source=footer>
>>>>>>> .
>>>>>>>
>>>>>> --
>> The deal.II project is located at http://www.dealii.org/
>> For mailing list/forum options, see
>> https://groups.google.com/d/forum/dealii?hl=en
>> ---
>>
> You received this message because you are subscribed to the Google Groups
>> "deal.II User Group" group.
>> To unsubscribe from this group and stop receiving emails from it, send an
>> email to dealii+un...@googlegroups.com.
>> To view this discussion on the web visit
>> https://groups.google.com/d/msgid/dealii/16856316-a225-4d92-9049-e48ffc6c81b4n%40googlegroups.com
>>
>> <https://groups.google.com/d/msgid/dealii/16856316-a225-4d92-9049-e48ffc6c81b4n%40googlegroups.com?utm_medium=email&utm_source=footer>
>> .
>>
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/12abdfe1-f16e-4c50-a9ee-28b5465ffaa8n%40googlegroups.com.
