Re: [Nektar-users] Flow around a sphere
Hi Silvia, Your setup seems to run fine for me on master. Could you please let me know the exact command you are running and the seg fault it produces and I will try to figure out what is going on? The composite ID in the .xml corresponds to the physical surface/volume tag in the .msh file, which in turn is filled from the .geo file if a numeric tag is specified. In your case as you have provided a char-expression instead, the .msh file is assigned a unique numeric tag automatically in the order of the physical surface/volume definitions in the .geo file. Here that is from 1-5 and the corresponding part in the .msh file is: $PhysicalNames 5 2 1 "inlet" 2 2 "outlet" 2 3 "walls" 2 4 "sphere" 3 5 "domain" $EndPhysicalNames This then corresponds to the C ID="1" to C ID="5" you see in the .xml file. You can assign a specific numeric tag in the .geo file by adding it after the char-expression e.g. Physical Surface("inlet", 101) = {2}; would correspond to 101 in the .msh file and then also in the .xml file as C ID="101". Hopefully that helps explain it! Thanks, Ed From: <nektar-users-bounces@imperial.ac.uk> on behalf of "Silvia Ceccacci (HDR)" <silvia.ceccacci@hdr.mq.edu.au> Date: Friday, 10 April 2020 at 03:09 To: nektar-users <nektar-users@imperial.ac.uk> Subject: [Nektar-users] Flow around a sphere This email from silvia.ceccacci@hdr.mq.edu.au originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list<https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fspam.ic.ac.uk%2FSpamConsole%2FSenders.aspx&data=02%7C01%7Cel326%40exeter.ac.uk%7C9859c2f302544ccd927508d7dcf43169%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C637220813712022384&sdata=zGn4Jw6E2G0mhsf0Xd%2BUqTArs12FH3kokVt0URTMwTo%3D&reserved=0> to disable email stamping for this address. Dear all, I hope this email finds you all well. I am trying to simulate flow around a sphere. I get an error of segmentation fault. I am generating a mesh from gmsh, converting it into Nektar format, then using the IncNavierStokesSolver. I am wondering whether it may be a problem with the association of composite IDs in boundary conditions. However I do not seem to find the problem. I attach the files I am using. If anyone could have a look and help me understand what I am doing wrong I would be very grateful. Anyways, I also have another question that may be related with the issue above. What is the way to know which composite ID corresponds to the related physical surface of your geometry? In other words, is there a way to associate a composite ID to a specific surface in order to assign boundary conditions? For example, in my case, in file s_session.xml, I say C[1] (ID="0") is the inlet surface because in file s.xml it is the first composite, which, I am assuming, corresponds to the first physical surface I have defined in file s.geo (Physical Surface("inlet")={2};). Is this correct? I hope my questions make sense. Thank you in advance for your help, and Happy Easter! Kind regards, Silvia
Hi Ed, I have sorted out the issue on the simulation. There was a mistake in the definition of the parameter NumSteps, as I had not defined the FinalTime. Thank you for the explanation on the composite IDs and boundary conditions. I was very helpful. Kind regards, Silvia ________________________________ From: Laughton, Edward <el326@exeter.ac.uk> Sent: Friday, 10 April 2020 9:33 PM To: Silvia Ceccacci (HDR) <silvia.ceccacci@hdr.mq.edu.au>; nektar-users <nektar-users@imperial.ac.uk> Subject: Re: [Nektar-users] Flow around a sphere Hi Silvia, Your setup seems to run fine for me on master. Could you please let me know the exact command you are running and the seg fault it produces and I will try to figure out what is going on? The composite ID in the .xml corresponds to the physical surface/volume tag in the .msh file, which in turn is filled from the .geo file if a numeric tag is specified. In your case as you have provided a char-expression instead, the .msh file is assigned a unique numeric tag automatically in the order of the physical surface/volume definitions in the .geo file. Here that is from 1-5 and the corresponding part in the .msh file is: $PhysicalNames 5 2 1 "inlet" 2 2 "outlet" 2 3 "walls" 2 4 "sphere" 3 5 "domain" $EndPhysicalNames This then corresponds to the C ID="1" to C ID="5" you see in the .xml file. You can assign a specific numeric tag in the .geo file by adding it after the char-expression e.g. Physical Surface("inlet", 101) = {2}; would correspond to 101 in the .msh file and then also in the .xml file as C ID="101". Hopefully that helps explain it! Thanks, Ed From: <nektar-users-bounces@imperial.ac.uk> on behalf of "Silvia Ceccacci (HDR)" <silvia.ceccacci@hdr.mq.edu.au> Date: Friday, 10 April 2020 at 03:09 To: nektar-users <nektar-users@imperial.ac.uk> Subject: [Nektar-users] Flow around a sphere This email from silvia.ceccacci@hdr.mq.edu.au originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list<https://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fspam.ic.ac.uk%2FSpamConsole%2FSenders.aspx&data=02%7C01%7Cel326%40exeter.ac.uk%7C9859c2f302544ccd927508d7dcf43169%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C637220813712022384&sdata=zGn4Jw6E2G0mhsf0Xd%2BUqTArs12FH3kokVt0URTMwTo%3D&reserved=0> to disable email stamping for this address. Dear all, I hope this email finds you all well. I am trying to simulate flow around a sphere. I get an error of segmentation fault. I am generating a mesh from gmsh, converting it into Nektar format, then using the IncNavierStokesSolver. I am wondering whether it may be a problem with the association of composite IDs in boundary conditions. However I do not seem to find the problem. I attach the files I am using. If anyone could have a look and help me understand what I am doing wrong I would be very grateful. Anyways, I also have another question that may be related with the issue above. What is the way to know which composite ID corresponds to the related physical surface of your geometry? In other words, is there a way to associate a composite ID to a specific surface in order to assign boundary conditions? For example, in my case, in file s_session.xml, I say C[1] (ID="0") is the inlet surface because in file s.xml it is the first composite, which, I am assuming, corresponds to the first physical surface I have defined in file s.geo (Physical Surface("inlet")={2};). Is this correct? I hope my questions make sense. Thank you in advance for your help, and Happy Easter! Kind regards, Silvia
participants (2)
-
Laughton, Edward
-
Silvia Ceccacci (HDR)