[bgulzar22]

Hello everyone, I need guidance on applying a temperature distribution on a wall. Could someone please provide the procedure or point me to relevant documentation?

[TKatt]

You want to use the python wrapper


Please take a look at TestCases/py_wrapper/custom_heat_flux/run_ad.py and probably better TestCases/py_wrapper/custom_load_fea/run.py and take it from there. The code is from the custom_load_fea case and does eseentially what you want - which is to apply a spatially varying temperature BC.

  # Apply a load based on the coordinates.
  if nVertex > 0:
    MarkerCoords = SU2Driver.MarkerCoordinates(MarkerID)
    L = 0.5
    dx = L / 16 # known from mesh settings in this case.
    for iVertex in range(nVertex):
      x = MarkerCoords(iVertex, 0)
      nodalForce = (2 * x / L) * dx
      # Half load due to half dx on first and last node.
      if abs(x) < 1e-6 or abs(x - L) < 1e-6:
        nodalForce = nodalForce / 2
      SU2Driver.SetMarkerCustomFEALoad(MarkerID, iVertex, (0, nodalForce))

[bgulzar22]

Originally Posted by TKatt

You want to use the python wrapper


Please take a look at TestCases/py_wrapper/custom_heat_flux/run_ad.py and probably better TestCases/py_wrapper/custom_load_fea/run.py and take it from there. The code is from the custom_load_fea case and does eseentially what you want - which is to apply a spatially varying temperature BC.

  # Apply a load based on the coordinates.
  if nVertex > 0:
    MarkerCoords = SU2Driver.MarkerCoordinates(MarkerID)
    L = 0.5
    dx = L / 16 # known from mesh settings in this case.
    for iVertex in range(nVertex):
      x = MarkerCoords(iVertex, 0)
      nodalForce = (2 * x / L) * dx
      # Half load due to half dx on first and last node.
      if abs(x) < 1e-6 or abs(x - L) < 1e-6:
        nodalForce = nodalForce / 2
      SU2Driver.SetMarkerCustomFEALoad(MarkerID, iVertex, (0, nodalForce))

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
Could you please look into it and please lend some help. I will be grateful.