Hello Apurva, I noticed that you are encountering some challenges when it comes to defining the boundary conditions for your heat transfer problem. While you have successfully defined the 'Convection coefficient' as a function handle, it seems that you are facing difficulty in doing the same for the 'AmbientTemperature'. According to the documentation, the 'AmbientTemperature' is expected to be of type 'double', which means it cannot be defined as a function handle.
To work around this limitation, one approach you can consider is reapplying the boundary conditions individually. By doing so, you can achieve non-constant ambient temperatures on the faces of your 3D geometry. Here's a suggested modification to your approach:
% Define the convection coefficient function handle
h_c = [10, 15, 20, 25]; % Example convection coefficient values
cold_coefficient = @(location, state) getcoldcoefficient(location, state, h_c);
% Define the ambient temperature vector
Tc = [25, 30, 35, 40]; % Example ambient temperature values
% Create the thermal model
model2 = createpde('thermal', 'transient');
% Define your 3D geometry (example: a cube)
L = 1; % Length of the cube
geometry = [3; 4; 0; L; 0; L; 0; L; 0; L; 0; L];
% Create a geometry object and assign it to the model
g = geometryFromEdges(model2, 'Edges', geometry);
model2.Geometry = g;
% Generate a mesh
generateMesh(model2);
% Apply the non-constant boundary conditions
faceIndices = [7, 8, 9, 10];
thermalBC(model2, 'Face', faceIndices, 'ConvectionCoefficient', cold_coefficient);
% Set the ambient temperature for each face individually
for k = 1:numel(faceIndices)
faceID = faceIndices(k);
thermalBC(model2, 'Face', faceID, 'AmbientTemperature', Tc(k));
end
In this modified code, after defining the geometry and generating the mesh, we apply the non-constant convection coefficient boundary condition using the “cold_coefficient” function handle.
Then, we iterate over each face index and apply the individual ambient temperatures using a loop. The “thermalBC” function is called for each face separately, specifying the face ID and the corresponding ambient temperature from the “Tc” vector.
Hope It helps!
