2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


В этом разделе нельзя создавать новые темы.

Если Вы хотите задать новый вопрос, то не дописывайте его в существующую тему, а создайте новую в корневом разделе "Помогите решить/разобраться (М)".

Если Вы зададите новый вопрос в существующей теме, то в случае нарушения оформления или других правил форума Ваше сообщение и все ответы на него могут быть удалены без предупреждения.

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

Обязательно просмотрите тему Правила данного раздела, иначе Ваша тема может быть удалена или перемещена в Карантин, а Вы так и не узнаете, почему.



Начать новую тему Ответить на тему
 
 Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение12.02.2013, 20:13 


05/02/13
4
Дано двумерное уравнение теплопроводности:
$a \frac{\partial u}{\partial t}$ = $\frac{\partial^2 u}{\partial x^2}$ + $\frac{\partial^2 u}{\partial y^2}$, $a = \frac {\rho C_v} \lambda $
где $a = a_1$ или $a_2$, в зависимости от области (рисунок).
В точках разрыва задается условие идеального контакта по температуре:
$T_l$=$T_r$;
$\lambda_l$ $\frac {\partial u} {\partial x}$=$\lambda_r$ $\frac {\partial u} {\partial x}$.
Нигде не могу найти как апроксимировать идеальный контакт по температуре. Подскажите формулы или в какой литературе их можно поискать?

Для одномерного случая я видел такую формулу:
$u_i = \frac {\lambda_l} {\lambda_l + \lambda_r} u_{i-1} + \frac {\lambda_r} {\lambda_l + \lambda_r} u_{i+1} $
$u_i$, соответственно, точка разрыва.
Но там для решения использовалась явная схема. Если для решения использовать неявную разностную схему (например, схему переменных направлений), то как апроксимировать условие идеального контакта не понятно.
Изображение

 Профиль  
                  
 
 Re: Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение13.02.2013, 11:34 
Заслуженный участник
Аватара пользователя


01/08/06
2401
Уфа
Чтобы с такими вещами не мучаться, можно идти к посконным корням уравнения теплопроводности, а именно — к интегральным уравнениям баланса, от которых получаются как сам дифур, так и условия неразрывности теплового потока (т.е. их не надо дополнительно рассматривать). К слову, если у нас какая-то граница теплоизолирована (т.е. заданы краевые условия II рода с нулём в правой части), то и эти условия получаются "автоматически".

Уравнения баланса аппроксимируются не сказать, чтобы сильно сложнее, чем дифференциальные, хотя повозиться, конечно, придётся. Зато мы получаем (дополнительно) более ясное представление о моделируемом процессе, а также полезное свойство консервативности схемы.
Сетку при этом нужно проектировать так, чтобы внутри ячеек коэффициенты не терпели разрывов, только на границах. Уравнения для узла получаются как сумма интегральных соотношений по примыкающим к узлу ячейкам (эти интегральные соотношения в конечном итоге выражаются линейно через значения в узлах сетки). Интегрирование по времени ведётся от текущего узла к будущему для явной схемы и к прошлому для неявной. Как-то так.

 Профиль  
                  
 
 Re: Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение13.02.2013, 16:58 


05/02/13
4
Я "прямого" способа численного решения задачи нет?

 Профиль  
                  
 
 Re: Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение13.02.2013, 19:28 
Заслуженный участник


11/05/08
31509
Вот Вы жалуетесь, что уравнение многомерно. Между тем разрыв в нём -- абсолютно одномерен. И, следовательно, все рецепты на этот счёт (которые Вы, как сами говорили, знаете) -- действуют.

А по существу -- да. Если коэффициент теплопроводности разрывен, то надо переходить к вариационной постановке задачи. И уже из неё вытягивать разностные уравнения. Они всегда будут так или иначе работать, но эффективно -- лишь в том случае, если узлы придутся на линию разрыва.

 Профиль  
                  
 
 Re: Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение13.02.2013, 19:49 


05/02/13
4
Ок, разрыв одномерен. Но станадртным способом решение двух- и трех-мерного уравнения теплопроводности являются схемы расщепления, которые являются неявными.
Для явной схемы, я понимаю, можно заранее найти $u_{i-1}$ и $u_{i+1}$ и по ним вычислить $u_i$ в точке разрыва. В случае метода переменных напралений (для двумерного случая), значения температуры находятся из системы линейных уранений. Корректно ли будет после решения данной системы, перестичать значения $u_i$ в точках азрыва по формулам, которые я приводил в первом посте?
Может быть вы знаете какую-нибудь литературу, где освещены подобные задачи?
По поводу перехода к вариационной постановке задачи, вопрос тот же - где можно почитать подробнее?

 Профиль  
                  
 
 Re: Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение14.02.2013, 10:00 
Заслуженный участник
Аватара пользователя


01/08/06
2401
Уфа
Alex2040 в сообщении #683014 писал(а):
Для одномерного случая я видел такую формулу:
$u_i = \frac {\lambda_l} {\lambda_l + \lambda_r} u_{i-1} + \frac {\lambda_r} {\lambda_l + \lambda_r} u_{i+1} $
$u_i$, соответственно, точка разрыва.
Но там для решения использовалась явная схема. Если для решения использовать неявную разностную схему (например, схему переменных направлений), то как апроксимировать условие идеального контакта не понятно.
Похоже не на непонятки с разрывными коэффициентами, а на непонимание, что такое вообще неявные схемы и с чем их едят. Или непонимание метода переменных направлений. Вот Вы выписали же уравнение, связывающее $u_{i-1}$, $u_i$ и $u_{i+1}$. Ну да, в явной схеме оно нам даёт возможность вычислить $u_i$ по простой формуле. А в неявной даёт нужное уравнение в систему. Которое только формулами для коэффициентов отличается от остальных уравнений. Главное — что оно есть, выписано. В чём проблема? Вы не понимаете, что с этой системой дальше делать? Или у Вас трудности с применением хорошо работающего в одномерном случае метода прогонки к многомерному случаю, т.е. с пониманием метода переменных направлений? Или Вы саму схему хотите заранее факторизовать по направлениям введением промежуточного слоя? Последний вопрос от предпоследнего отличается, это разные методы. Но в любом случае на этом этапе вопрос с разрывными коэффициентами уже решён.

 Профиль  
                  
 
 Re: Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение21.02.2013, 14:08 


05/02/13
4
Хорошо, положим $u_k = u_{k-1}d_1 + u_{k+1}d_2$. Тогда матрица системы преобразуется следущим образом:

$\begin{matrix}
\begin{matrix}
 &  &  &  & k &  &  &  & 
\end{matrix}
\\ 
\begin{pmatrix}
\cdots & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots &  0\\ 
\vdots & 0 & a_{k-1} & b_{k-1} & c_{k-1} & 0 & 0 & 0 & \vdots  \\ 
\vdots & 0 & 0 & a_k & b_k & c_k & 0 & 0 & \vdots \\ 
\vdots &0 & 0 & 0 & a_{k+1} & b_{k+1} & c_{k+1} & 0 & \vdots\\ 
0  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots \\
\end{pmatrix}
\end{matrix}
\sim$

$\begin{matrix}
\begin{matrix}
 &  &  &  & k &  &  &  & 
\end{matrix}
\\ 
\begin{pmatrix}
\cdots & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots &  0\\ 
\vdots & 0 & a_{k-1} & b_{k-1}+c_{k-1}d_1 & 0 & c_{k+1}d_2 & 0 & 0 & \vdots  \\ 
\vdots & 0 & 0 & a_k+b_kd_1 & 0 & c_k+b_kd_2 & 0 & 0 & \vdots \\ 
\vdots &0 & 0 & a_{k+1}d_1 & 0 & b_{k+1} + a_{k+1}d_2 & c_{k+1} & 0 & \vdots\\ 
0  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots  & \cdots \\
\end{pmatrix}
\end{matrix}
$

Тождественными преобразованиями преобразуем матрицу к трехдиагональному виду и так же решаем методом прогонки. Я правильно понимаю?

 Профиль  
                  
 
 Re: Многомерное ур-е теплопроводности с разравными коэфф.
Сообщение21.02.2013, 14:51 
Заслуженный участник
Аватара пользователя


23/08/07
4629
Нов-ск
Alex2040 в сообщении #686635 писал(а):
Я правильно понимаю?

В каждой внутренней точке (включая точки разрыва коэффициента) имеем уравнение вида
$A_j u_{j-1}+B_j u_{j}+C_j u_{j+1}=D_j.$
Плюс два уравнения, аппроксимирующие граничные условия справа и слева.

Решайте. Хоть прогонкой, хоть подгонкой. Где тут правильно или неправильно понимать?

-- Чт фев 21, 2013 16:29:01 --

$$\frac{\partial u}{\partial t}=\frac{\partial}{\partial x}\lambda \frac{\partial u}{\partial x}+\frac{\partial}{\partial y}\lambda \frac{\partial u}{\partial y}$$
$$\frac{u^{n+1}_{jk}-u^{n}_{jk}}{\tau}=\frac{1}{h_x}\left(\lambda_{j+1/2}\frac{u^{n+1}_{j+1}-u^{n+1}_{j}}{h_x}-
 \lambda_{j-1/2}\frac{u^{n+1}_{j}-u^{n+1}_{j-1}}{h_x} \right)$$
$$+\frac{1}{h_y} \left( \lambda_{k+1/2}\frac{u^{n+1}_{k+1}-u^{n+1}_{k}}{h_y}-
 \lambda_{k-1/2}\frac{u^{n+1}_{k}-u^{n+1}_{k-1}}{h_y}\right) $$

На линии разрыва $\lambda_{j}=(\lambda_{R}+\lambda_{L})/2$, $\lambda_{j\pm1/2}$ - понятно что

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 8 ] 

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
cron
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group