{ Th2LS_1.PDE  - program for FlexPDE }

{русский язык не поддерживается для заголовков - названия на английском}

  title "Heat flow in 2D region" {2D задача теплопроводности}

  select
     errlim = 1e-3	{точность вычислений}

  Variables
     Temp 	     {неизвестная ф-я - температура}

  definitions
     hl=1.4               {высота области}
     wl=0.8               {ширина области}
     wwl=0.2            {ширина стенок}
     hal=0.7               {дополнительный размер}
     K = 46.7             {коэффициент теплопроводности}
     Source = 0         {интенсивность источников тепла}
     hf=60	    {коэффициент теплообмена film coefficient}
     Text = 0              {температура окужающей среды - bulk temperature of the adjacent medium}
     Tinp1 = 90          {подаваемая температура на нижнюю левую грань}
     Tinp2 = 180          {подаваемая температура на нижнюю правую грань}
     flux = -K*grad(Temp) {вектор потока тепла}

    pel=(hal/2)/(wl-wwl)  {параметр эллиптичности}

  equations
     div(flux) = Source    {стационарное уравнение теплопроводности}
		   { (div(flux) = 0, в отсутствии источников тепла Source = 0) }

  boundaries               {расчетная область и граничные условия}
     Region 1              {регион 1 - будет одна расчетная область}
	start (0,0)
	value(Temp)=Tinp1  {главное граничное условие: температура на на нижнюю левую грань}
	line to (wwl,0)
	natural(Temp)=0       {естественное граничное условие: тепловой поток равен нулю - теплоизолированная граница}
	line to (wwl,hl-hal-wwl) {далее повторяется предыдущее условие, если не написано новое}
	line to (wl-wwl,0)
	value(Temp)=Tinp2  {главное граничное условие: температура на нижнюю правую грань}
	line to (wl,0)
	natural(Temp)=0 {естественное граничное условие: тепловой поток равен нулю - теплоизолированная граница}
	line to (wwl,hl-hal)
	arc (center=wwl,hl-hal/2) to (wl,hl-hal/2)   {Дуга эллипса: лучше строить не более 90 градусов}
	arc (center=wwl,hl-hal/2) to (wwl,hl)
	natural(Temp)=-hf*(Text-Temp)   {естественное граничное условие: конвективный теплообмен на верхней грани}
	line to (0,hl)
	natural(Temp)=0 {естественное граничное условие: тепловой поток равен нулю - теплоизолированная граница}
	line to close {возврат в стартовую точку для замыкания границы}

	natural(Temp)=0
	start (wwl,hl-hal/2-(wl-2*wwl)*pel)     {новая стартовая точка для построения новой замкнутой границы - моделирование отверстия}
	arc (center=wwl,hl-hal/2) to (wl-wwl,hl-hal/2)
	arc (center=wwl,hl-hal/2) to (wwl,hl-hal/2+(wl-2*wwl)*pel)
	line to close

  monitors
     contour(Temp) {показ распределения температуры в процессе решения}

  plots            {графический вывод результатов}
     grid(x,y) as "Finite element mesh" {русский язык не поддерживается для заголовков - названия на английском}
     contour(Temp) painted as "Temperature" {Картина распределения температуры  - градиентая заливка с опцией painted}
     contour(magnitude(flux))  painted  as "Heat flow vector magnitude" {Картина распределения модуля вектора потока тепла}

      vector(flux) as "Heat Flow vector" {Картина распределения вектора потока тепла}



  end