Пример. Тепловой поток на квадратной плите

Функции > Решение и оптимизация > Решатели дифференциальных уравнений > Пример. Тепловой поток на квадратной плите - II

Пример. Тепловой поток на квадратной плите - II

Найдите решение для стационарного распределения температуры квадратной плиты с помощью решателя дифференциальных уравнений в частных производных relax.

Решение уравнения Пуассона

Решите уравнение теплопроводности, где значения функции источника известны и граничные условия являются ненулевыми.

Функция relax базируется на совершенно другом методе решения и, следовательно, требует другой набор аргументов.

1. Определите пять квадратных матриц a, b, c, d и e, которые будут содержать коэффициенты для аппроксимации лапласиана:

Эти массивы могут иметь любой указанный вами размер. Чем больше размер массивов, тем мельче будет сетка в решении.

2. Определите размер квадратной плиты:

3. Определите коэффициенты:

4. Определите мощность и позицию постоянного источника.

5. Определите квадратную матрицу f, чей размер равен размеру сетки, которая будет содержать известные граничные значения функции F(x,y) и значения начального приближения для неизвестных внутренних значений.

◦ Граничное условие вдоль верхнего края:

◦ Граничное условие вдоль нижнего края:

◦ Граничное условие боковых краев:

6. Определите переменную r, представляющую Jacobi spectral radius, как вещественное число со значением от 0 до 1.

Этот параметр управляет сходимостью алгоритма. Если появляется сообщение об ошибке "слишком много итераций", попытайтесь уменьшить r.

7. Вызовите функцию relax:

<region id="ID0ERU15" top="1516.8000000000002" left="38.400000000000006" height="25.6" width="204.10333333333335" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math resultRef="30">
    <define xmlns="http://schemas.mathsoft.com/math50">
      <id xml:space="preserve" labels="VARIABLE">F</id>
      <apply>
        <id xml:space="preserve" labels="FUNCTION" label-is-contextual="true">relax</id>
        <sequence>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">a</id>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">b</id>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">c</id>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">d</id>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">e</id>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">S</id>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">f</id>
          <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">r</id>
        </sequence>
      </apply>
    </define>
  </math>
</region>

8. Создайте 3D-график, чтобы показать распределение тепла по квадратной плите.

Щелкните для копирования этого выражения

<region id="ID0EBW15" top="1747.2000000000003" left="38.400000000000006" height="301.40000715255735" width="297.6000071525574" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white">
    <plot3D>
      <traces>
        <trace3D resultRef="32" num-of-points="41" SurfaceFill="#00FFFFFF">
          <traceStyle color="#FFED1D2F">surface</traceStyle>
        </trace3D>
      </traces>
      <graph-size width="240.00000715255737" height="240.00000715255737" />
      <axes>
        <xAxis rank="1" start="0" end="32">
          <axisLine />
          <axisGrid>
            <gridFrequency>5</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <numberFormat>
            <general precision="3" />
          </numberFormat>
          <xyDomain auto-scale="true" allow-negative-values="true" scale-type="linear">
            <startValue />
            <endValue />
          </xyDomain>
        </xAxis>
        <yAxis rank="1" start="0" end="32">
          <axisLine />
          <axisGrid>
            <gridFrequency>5</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <numberFormat>
            <general precision="3" />
          </numberFormat>
          <xyDomain auto-scale="true" allow-negative-values="true" scale-type="linear">
            <startValue />
            <endValue />
          </xyDomain>
        </yAxis>
        <zAxis rank="1" start="-3" end="4.5">
          <axisLine />
          <axisGrid>
            <gridFrequency>6</gridFrequency>
            <gridLabels display="true" />
            <gridLines />
            <tickMarks display="true" />
          </axisGrid>
          <axisLabel />
          <numberFormat>
            <general precision="3" />
          </numberFormat>
          <zDomain auto-scale="true" allow-negative-values="true" scale-type="linear">
            <startValue />
            <endValue />
          </zDomain>
        </zAxis>
      </axes>
      <plotEquations>
        <plotEquation>
          <math>
            <apply xmlns="http://schemas.mathsoft.com/math50">
              <neg />
              <id xml:space="preserve">F</id>
            </apply>
          </math>
          <math>
            <placeholder xmlns="http://schemas.mathsoft.com/math50" />
          </math>
        </plotEquation>
      </plotEquations>
      <CameraTransformMatrix Columns="4" Rows="4">
        <Cells>
          <Cell>-0.052623873051966662</Cell>
          <Cell>0.043294250885890617</Cell>
          <Cell>0.99767546618388891</Cell>
          <Cell>0</Cell>
          <Cell>0.095704097834471175</Cell>
          <Cell>0.99467978821191028</Cell>
          <Cell>-0.038116198399097093</Cell>
          <Cell>0</Cell>
          <Cell>-0.99401783366431629</Cell>
          <Cell>0.093475808436936883</Cell>
          <Cell>-0.056487340124674094</Cell>
          <Cell>0</Cell>
          <Cell>0</Cell>
          <Cell>0</Cell>
          <Cell>0</Cell>
          <Cell>1</Cell>
        </Cells>
      </CameraTransformMatrix>
      <CameraType>Orthographic</CameraType>
    </plot3D>
  </plot>
</region>

9. Создайте контурный график, чтобы показать линии постоянной температуры.

Было ли это полезно?