[bgulzar22]

Dear SU2 Community,

I am currently working on implementing a simulation using SU2 and have encountered a challenge that I hope some of you might be able to assist with. Below is a brief overview of my problem:

Problem Overview:
I am working on a simulation that involves setting custom temperature boundary condition on a specific marker using Python scripting. However, I'm facing difficulties in properly setting this boundary condition based on temperature data from an external file. I have written a Python script to read the temperature data, assign it to the appropriate marker, and run the simulation. However, despite my efforts, I am encountering issues with the implementation.

To provide more context and allow for a better understanding of the problem, I have prepared a comprehensive package containing the following:

Configuration file
Mesh file
Temperature data file
Python script for setting the boundary condition and running the simulation
You can access the files package through this link:https://drive.google.com/drive/folde...usp=drive_link

Request for Assistance:
I am reaching out to the SU2 community in the hopes that someone with expertise in Python scripting for SU2 simulations can review my implementation and provide insights into resolving the issues I'm facing. Specifically, I would greatly appreciate any guidance or suggestions on how to properly set the custom temperature boundary conditions using Python scripting.

If anyone is willing to lend their expertise and offer assistance, I would be extremely grateful. Your support would not only help me overcome this obstacle but also contribute to advancing my understanding and proficiency in utilizing SU2 for my simulations.

Thank you in advance for your time and assistance.

Best regards,

[bigfootedrockmidget]

I haven't looked at your setup, but the way I would have done it is take this case as a starting point:
https://github.com/su2code/SU2/blob/...T_FlatPlate.py


Then change the time-dependent temperature BC for a spatially dependent BC using something like this:

marker_coords = driver.MarkerCoordinates(marker_id) 
x = marker_coords(i_vertex, 0) 
y = marker_coords(i_vertex, 1)
temperature = temperature_function(x,y)

And in your case read your file in temperature_function and interpolate the values at (x,y)

[bigfootedrockmidget]

It looks like you are doing something similar so first take take the working testcase from the repository, verify that it works, then modify it so it uses spatially varying temperature. Just T=x+y or something simple like that. Then get temperature from a file. Verify all intermediate steps.

[bgulzar22]

Originally Posted by bigfootedrockmidget

It looks like you are doing something similar so first take take the working testcase from the repository, verify that it works, then modify it so it uses spatially varying temperature. Just T=x+y or something simple like that. Then get temperature from a file. Verify all intermediate steps.

Thank you for your response , I have run the test case for unsteady one, test case works fine but when I try to run my case I face a segmentation error. Could you please look into it and please let me know where does the fault lie , and what mistake I am making. Please I would be grateful.

[bigfootedrockmidget]

Does your code run fine when you just set temperature=300 in your python script?
Does your code run fine when you set temperature=300+x+y ?

[bgulzar22]

Originally Posted by bigfootedrockmidget

Does your code run fine when you just set temperature=300 in your python script?
Does your code run fine when you set temperature=300+x+y ?

No, I get the same segmentation error when I run the python script, but when I run the simulation without the python script , it runs fine with isothermal boundary condition.

[bgulzar22]

Originally Posted by bigfootedrockmidget

Does your code run fine when you just set temperature=300 in your python script?
Does your code run fine when you set temperature=300+x+y ?

Now the simulation works when I use the updated python script https://drive.google.com/file/d/13M6...usp=drive_link
but I am not sure if it is using temperature distribution correctly. As in my simulation I have given temperature at WALL as 421K(ISOTHERMAL BC). Could you suggest me a way, so that I can check whether the updated temperature distribution is applied on the WALL.
Thank YOU in advance.