Скачиваний:
34
Добавлен:
15.06.2014
Размер:
483.33 Кб
Скачать

Пример точного решения ДУ (продолжение3)

u

 

1;

u

 

x 1 0;

 

g ( x ) 1;

f ( x ) 1 .

 

0

0; q

0

1;

1

1; q

1

0

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем сразу

c1

g ( 0 )

 

 

1;

p ( x )

 

f ( x ) d x x

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p (1)

 

1

 

 

 

1

• Из

q

1

(

c1

)

1

( c2

 

 

 

 

 

 

 

 

 

 

 

 

 

g (1)

 

 

 

 

 

 

 

g (1)

 

 

0

 

 

 

 

 

 

 

 

 

 

1

p ( x ) d x c1 g ( x )

0

1 d x ) 1 ; g ( x )

Получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 (

1

1

1

) 1( c2

1

x

d x 1

1

 

1

d x ) 0; c2 1 .5

 

 

 

 

 

 

 

1

 

1

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

Окончательно

 

x

p ( x )

 

 

 

 

 

x

1

 

x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u ( x ) c2

 

 

 

d x c1

 

d x 1 .5

 

x

 

 

 

g ( x )

 

g ( x )

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

05.01.2011

11

Постановка задач для ДУ в частных производных (ДУЧП)

Большинство окружающих нас физических полей описывается ДУЧП второго порядка, или системами таких ДУЧП.

Если процессы не слишком интенсивны, то можно ограничиться линейными дифференциальными уравнениями второго порядка

aij

2 u

bi

u

cu f

x1 x ,

x 2 y ,

x3 z ,

x 4 t .

xi x j

 

ij

i

xi

 

 

 

 

Свойства решений существенно зависят от коэффициентов, стоящих при старших производных. С помощью соответствующей замены независимых переменных уравнение может быть приведено к одному из трех типов

Параболические

Гиперболические

Эллиптические

05.01.2011

12

Параболические

• Канонический вид

 

u

 

2

u

 

2

u

 

2

u

f

 

 

 

x 2

y 2

z 2

 

t

 

 

 

Уравнение

 

 

 

 

 

 

 

 

 

теплопроводности

 

u

 

 

 

u

 

 

a

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

t

 

x

x

 

y

u

g y

 

 

u

f

 

g

 

 

 

z

z

 

 

распространение возмущения

 

 

 

u

 

u

 

2 u

 

t=0

 

 

 

t

x 2

 

 

 

 

t=t1>0

 

 

 

 

 

 

 

 

t=t2>t1

 

 

 

0

 

 

 

 

0

z

 

 

 

 

05.01.2011

13

 

 

 

 

 

Гиперболические

 

 

 

 

 

Канонический вид

 

2 u

 

 

2 u

 

2 u

 

2 u

f

 

 

 

 

 

 

 

 

 

 

 

t 2

x 2

y 2

z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Волновое

 

2 u

 

 

 

 

u

 

 

 

 

 

u

 

 

u

 

 

уравнение

 

a

2

 

g

 

 

y

g

 

g

 

f

Д’Аламбера

t

 

 

x

 

x

 

 

 

 

y

 

z

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Распространение одномерного возмущения

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

t=0

 

 

 

 

 

2 u

 

2 u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

2

z

2

 

 

 

 

 

t=t2

 

t=t1

 

 

 

 

t=t1

t=t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

Точное решение

 

u

 

f ( t

 

 

x )

 

 

 

 

 

 

05.01.2011

 

 

 

 

 

 

 

 

 

 

 

 

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эллиптические (стационарные процессы)

• Канонический вид

2 u

 

2 u

 

2 u

f

 

 

 

 

 

x 2

y 2

z 2

 

 

 

 

Уравнение

Пуассона

(Лапласса)

 

 

u

 

 

u

 

 

 

u

f

 

g

 

 

 

g

 

 

 

g

 

 

 

 

 

 

 

 

x

x

 

y

y

 

 

z

z

 

При постановке задач для ДУ в частных производных обычно

требуется найти распределение , удовлетворяющее ДУ в некоторой области пространства с границей Г.

y

Г

x

05.01.2011

15

Граничные условия

Общее решение ДУЧП содержит произвольные дифференцируемые функции, например

Решением

 

2

u

 

2

u

является

u f (t z )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

2

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Начальные условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

t 0 U

0

( x , y , z )

 

u

 

 

U 0 ( x , y , z )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

t 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Граничные условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

первого рода

 

 

 

второго рода

 

 

 

третьего рода

 

 

Дирихле

 

 

 

 

 

Неймана

 

 

 

 

Ньютона

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

Г ( Г )

 

 

 

 

 

u

 

 

( Г )

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

( Г )

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

 

 

 

 

 

 

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

05.01.2011

16

Пример решения

 

 

 

u

 

2 u

1

 

 

 

 

 

u

 

 

 

 

 

 

 

 

t

z 2

 

 

 

 

 

 

 

 

 

 

 

1;

u

 

x 1 0;

u

 

t 0 0

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

05.01.2011

17

Как получают дифференциальные уравнения

ДУ являются следствием фундаментальных законов природы и применения интегральных теорем

Например, Французский математик и физик Жан Батист Жозеф Фурье в 1822 г. установил закон теплопроводности (закон Фурье):

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

 

 

 

 

 

 

 

 

q T n

T

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

T гр а д

dS

В т

 

2

 

В т

 

n

м

 

 

 

 

 

 

 

 

 

 

м г р а д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

05.01.2011

18

Пример получения стационарного двумерного уравнения теплопроводности

используя закон Фурье

Выберем элементарный

 

y/2

объем V= x y

 

 

 

в центре которого находится

 

 

 

точка x,y .

- x/2

 

x/2

 

 

 

 

 

 

 

 

- y/2

Допустим, что в таком кубике в

 

 

единицу времени выделяется

Q V

qV ( x , y ) x y

количество теплоты, равное

05.01.2011

19

Составим уравнение баланса для элементарного объема

 

 

 

 

T

qV

x y

 

 

 

x

 

 

 

 

 

 

 

 

 

x x

2

T

x

x x

2

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

y

y

y

 

 

 

 

2

 

 

T

 

 

 

 

 

 

y

y

 

 

y

 

 

 

 

 

 

 

 

2

x

Разделим на x y:

 

 

 

 

T

 

 

 

T

 

 

 

 

 

 

T

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

y

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x

x

 

 

x x

x

 

 

 

y

y

 

 

 

y

qV

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Перейдем к пределу 0

T

xx

 

 

 

 

 

T

qV ( x , y )

 

 

 

 

 

 

 

 

y

 

 

 

y

 

 

 

 

05.01.2011

20