Steady State Heat Transfer: Input File Settings
Morris Jowas
Hello all,
I am a new deal.ii user pursuing a PhD in thermal management of printed circuit boards .I am trying to analyse a steady state heat transfer problem on a 3d rectangular domain. The domain has a convective heat flux boundary condition on the top and bottom surfaces. There is also a discrete heat source on the top surface.
I am inputting my geometry using as a vtk file using the read_vtk() function. My questions:
- Is there a way to include the boundary conditions in the input file?
- Is it possible to include the material properties of the cells in the input file. Each cell can potentially be anisotropic and have different thermal conductivity in the x, y and z directions.
From reading the documentation it seems that the vtk input file can include MaterialID and BoundaryID data, but it would be useful for me to also include the accompanying properties of the identified materials and boundaries.
I appreciate your assistance in advance, especially an example implementation of such a file.
Kind regards,
Morris Jowas
Wolfgang Bangerth
Morris,
> I am a new deal.ii user pursuing a PhD in thermal management of printed
> circuit boards .I am trying to analyse a steady state heat transfer
> problem on a 3d rectangular domain. The domain has a convective heat
> flux boundary condition on the top and bottom surfaces. There is also a
> discrete heat source on the top surface.
>
>
> - Is there a way to include the boundary conditions in the input file?
>
> - Is it possible to include the material properties of the cells in the
> input file. Each cell can potentially be anisotropic and have different
> thermal conductivity in the x, y and z directions.
>
> From reading the documentation it seems that the vtk input file can
> include MaterialID and BoundaryID data, but it would be useful for me to
> also include the accompanying properties of the identified materials and
> boundaries.
to include this sort of information into the mesh file. But, in
practice, that's not how people do things -- instead, one just puts
material and boundary ids into the mesh file, and then make the
association with material properties and what a particular boundary id
means in the solver.
In other words, in your assembly loop, you'd just query the material_id
of a cell, and based on that decide what the material properties of that
cell should be.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/