background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

1

METODY MATEMATYCZNE I 

STATYSTYCZNE 

W INŻYNIERII CHEMICZNEJ

Wykład 3

Elementy analizy pól skalarnych, 

wektorowych i tensorowych

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

2

Analiza pól  – Uwagi wstępne

Bardzo ważnym pojęciem w fizyce a także inżynierii chemicznej jest 

pojęcie pola.

Polem nazywamy pewną funkcję wielu zmiennych, w której 

argumentami są położenie i czas. 

Za pomocą pól opisuje się różne procesy zachodzące w przestrzeni 

i czasie.

pole

(po o enie, czas)

f

ł ż

W zależności od rodzaju wielkości jaką opisuje dane pole rozróżniamy

pola skalarne, wektorowe i tensorowe.

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

3

Analiza pól  – Uwagi wstępne

Stosunkowo prostą zmienną jest czas, który często identyfikowany 

jest ze zbiorem nieujemnych liczb rzeczywistych R

+

. Najczęściej czas 

jest oznaczany literą t.

Zatem tєR+tzn. 0≤t<∞.

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

4

Analiza pól  – Uwagi wstępne

Znacznie bardziej złożone jest zagadnienie opisu położenia. 
Przyjmuje się, że procesy zachodzą w przestrzeni trójwymiarowej 
a zatem do opisu położenia potrzebne są 3 składowe. W zależności 
od geometrii opisywanego zjawiska stosowane mogą być różne 
układy współrzędnych przestrzennych. 

Najczęściej stosowane są trzy rodzaje układów współrzędnych 
przestrzennych:

-

kartezjański układ prostokątny,

-

układ cylindryczny,

-

układ sferyczny.

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

5

Analiza pól  – Układy 

współrzędnych przestrzennych

1. Układ kartezjański.

Ustalenie położenia polega na wyborze w przestrzeni trzech 

wzajemnie prostopadłych  i przecinających się w jednym 

punkcie prostych określanych tradycyjnie jako osie 

współrzędnych x,y,z. 

Położenie „u” danego punktu w przestrzeni określa trójka 

liczb x,y,z będących rzutami punktu u na odpowiednie osie. 

Możemy to zapisać: u=[u

x

,u

y

,u

z

]

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

6

Analiza pól  – Układy 

współrzędnych przestrzennych

x

y

z

u

u

xy

u

x

u

y

u

z

O

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

7

Analiza pól  – Układy 

współrzędnych przestrzennych

1. Układ kartezjański.

Współrzędne układu kartezjańskiego są dowolnymi liczbami 

rzeczywistymi x,y,z:

-∞<x<∞

-∞<y<∞

-∞<z<∞

x

y

z

u

u

xy

u

x

u

y

u

z

O

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

8

Analiza pól  – Układy 

współrzędnych przestrzennych

2. Układ cylindryczny.

Podstawową rolę w układzie cylindrycznym odgrywa płaszczyzna z wyróżnioną

półprostą Ox oraz prostopadła do tej płaszczyzny i przechodząca przez punkt O

oś z. Położenie punktu u określamy za pomocą trójki liczb rzeczywistych 

u=[r

c

,φ,u

z

]=[r

c

,φ,z] gdzie:

x

z

u

u’

u

z

O

φ

r

c

r

c

– odległość punktu u’ będącego

rzutem punktu u na płaszczyznę

od punktu O, czyli długość 

odcinka u’O. r

c

≥0

φ – kąt między półprostą Ox 

a odcinkiem Ou’. 0≤φ<2π

u

z

=z – współrzędna rzutu punktu

u na oś z. -∞<z<∞

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

9

Analiza pól  – Układy 

współrzędnych przestrzennych

3. Układ sferyczny.

Podstawową rolę w układzie sferycznym odgrywa płaszczyzna z wyróżnioną

półprostą Ox oraz prostopadła do tej płaszczyzny i przechodząca przez punkt O

oś z. Położenie punktu u określamy za pomocą trójki liczb rzeczywistych 

u=[r

s

,φ,θ] gdzie:

x

z

u

u’

O

φ

r

s

r

s

– odległość punktu u będącego     

od punktu O, czyli długość 

odcinka uO. r

s

≥0

φ – kąt między półprostą Ox 

a odcinkiem Ou’. 0≤φ<2π

θ – kąt między dodatnim kierunkiem 

osi z a odcinkiem Ou.  0≤θ≤π

θ

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

10

Analiza pól  – Przeliczanie 

współrzędnych

cos

sin

c

c

x

r

y

r

z

z

Współrzędne w różnych układach można wzajemnie przeliczać. Szczególnie ważne

są wzory przeliczeniowe między układem kartezjańskim a cylindrycznym i sferycznym.

Wzory te można otrzymać na podstawie elementarnych relacji geometryczno –

trygonometrycznych.

x

z

u

u’

u

z

O

φ

r

c

y

2

2

0

0

arctan

0

0

0

0

2

0

0

c

r

x

y

y

x

z

z

gdzie

dla x

i y

dla x

dla x

i y

Układ cylindryczny –> układ kartezjański:

.

Układ kartezjański –> układ cylindryczny:

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

11

Analiza pól  – Przeliczanie 

współrzędnych

cos sin

sin sin

cos

s

s

s

x

r

y

r

z

r

Układ sferyczny –> układ kartezjański:

x

z

u

u’

u

z

O

φ

r

s

y

2

2

2

0

2

2

0

0

0

arctan

arctan

0

0

0

0

2

0

0

0

0

0

s

r

x

y

z

y

x

x

y

z

gdzie

dla x

i y

dla x

dla x

i y

dla z

dla z

θ

r

s

.

.

Układ kartezjański –> układ sferyczny:

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

12

Analiza pól  – Pole skalarne

R

s

t

u

f

s

)

,

(

Polem skalarnym nazywamy funkcję rzeczywistą położenia i czasu.

Pole skalarne jest więc funkcją 4 zmiennych: 3 przestrzennych i czasu. Zmienne 

przestrzenne zależą od stosowanego układu współrzędnych:

( , , , )

( , , , )

( , , , )

K

C

c

S

s

s

f

x y z t

s

f

r

z t

s

f r

t

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

13

Pole skalarne – powierzchnia 

ekwiskalarna

R

s

t

u

f

s

)

,

(

Powierzchnią ekwiskalarną nazywamy zbiór punktów przestrzennych

dla których wartości funkcji s w określonym czasie są stałe.

Przykład

Niech s będzie polem określonym za pomocą wzoru:

( )

( , )

s

u

t

s

f u t

const

A

2

2

2

s

x

y

z

Powierzchniami ekwiskalarnymi dla tego pola są powłoki kuliste (sfery) dla których 

odległości od początku układu są stałe. 

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

14

Analiza pól  – Pole wektorowe

3

)

,

(

R

w

t

u

f

w

Polem wektorowym nazywamy funkcję wektorową położenia i czasu.

Ponieważ wektor w przestrzeni R

3

ma 3 składowe, więc funkcja f musi w sposób

niezależny określać te składowe. Rodzaj tych składowych zależy od rodzaju

stosowanego układu współrzędnych:

[

,

,

]

[

,

,

]

[

,

,

]

c

s

x

y

z

r

z

r

w

w w w

w

w w w

w

w w w


background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

15

Analiza pól  – Pole wektorowe

Pole wektorowe jest zatem równoważne trzem funkcjom rzeczywistym 

czterech zmiennych. 

W zależności od rodzaju układu funkcje te mają różne postacie:

( , )

[

( , , , ),

( , , , ),

( , , , )]

( , )

[

( , , , ),

( , , , ),

( , , , )]

( , )

[

( , , , ),

( , , , ),

( , , , )]

K

x

y

z

C

r

c

c

z

c

S

r

s

s

s

w

f

u t

w x y z t w x y z t w z y z t

w

f u t

w r

z t w r

z t w r

z t

w

f u t

w r

t w r

t w r

t



background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

16

Analiza pól  – Pole tensorowe

( , )

T

f u t



Polem tensorowym nazywamy funkcję tensorową położenia i czasu.

W przypadku tensorów drugiego rzędu określonych w przestrzeni trójwymiarowej

tensor ma 9 składowych, więc funkcja f musi w sposób niezależny określać te 

składowe. Każda ze składowych jest funkcją 4 zmiennych.

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

17

Analiza pól  – Pole tensorowe

Pole tensorowe jest zatem równoważne 9 - ciu funkcjom rzeczywistym 

czterech zmiennych.

W układzie kartezjańskim popularny jest macierzowy zapis tensorów:

( , )

( , )

( , )

( , )

( , )

( , )

( , )

( , )

( , )

xx

xy

xz

yx

yy

yz

zx

zy

zz

f

u t

f

u t

f

u t

T

f

u t

f

u t

f

u t

f

u t

f

u t

f

u t



background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

18

Analiza pól  – Różniczkowanie

Pola jako funkcje można różniczkować. Ponieważ jednak są to funkcje 

dosyć złożone więc również różniczkowanie nie jest proste. 

Ogólnie różniczkowanie jest operatorem.

Istnieją różne operatory różniczkowania pól noszące różne nazwy. 

Najważniejsze z nich to:

-gradient pola skalarnego

dywergencja pola wektorowego

rotacja pola wektorowego 

laplasjan (operator Laplace’a – nie mylić z transformatą Laplace’a) 

pola skalarnego

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

19

Analiza pól  – Różniczkowanie

Formalne definicje poszczególnych operatorów są niezależne od rodzaju 

układu współrzędnych przestrzennych, jednakże nie nadają się one do 

obliczeń praktycznych (podobnie jak formalna definicja zwykłej pochodnej 

czy całki). 

W praktyce stosuje się wzory, których postać zależy od rodzaju układu 

współrzędnych. Z reguły najprostsze są wzory dla układu kartezjańskiego 

i te wzory przedstawię poniżej.

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

20

Analiza pól  – Gradient pola 

skalarnego

Operację gradientu wykonuje się na polu skalarnym. 

Wynikiem jest pole wektorowe. Dla układu kartezjańskiego mamy:

z

f

y

f

x

f

s

grad

t

z

y

x

f

s

,

,

)

(

)

,

,

,

(

Odpowiednie wzory określające gradient w innych układach 

współrzędnych można znaleźć w podręcznikach lub poradnikach 

matematycznych. 

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

21

Właściwości gradientu pola 

skalarnego

Można wykazać że: 

1. Operator gradientu jest liniowy tzn. obowiązuje wzór

gdzie – α i β dowolne liczby rzeczywiste, s

1

i s

2

dowolne pola skalarne

2. Wektor gradientu wskazuje kierunek, w którym wartość pola 

rośnie najszybciej.

3. Wektor gradientu jest prostopadły do powierzchni ekwiskalarnej. 

1

2

1

2

(

)

( )

( )

grad

s

s

grad s

grad s

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

22

Analiza pól  – Dywergencja pola 

wektorowego

Operację dywergencji wykonuje się na polu wektorowym. 

Wynikiem jest pole skalarne. Inne określenie tego operatora

to rozbieżność.

[

( , , , ),

( , , , ),

( , , , )]

( )

x

y

z

y

x

z

w

w x y z t w x y z t w x y z t

f

f

f

div w

x

y

z

Odpowiednie wzory określające dywergencję w innych układach 

współrzędnych można znaleźć w podręcznikach lub poradnikach 

matematycznych. 

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

23

Niektóre właściwości operatora 

dywergencji pola wektorowego

Można wykazać że:

1

2

1

2

(

)

(

)

(

)

div

w

w

div w

div w

1. Dywergencja jest operatorem liniowym tzn.

2. Zachodzi wzór:

(

)

( )

( )

div s w

s div w

grad s

w

gdzie α i β – dowolne liczby rzeczywiste, 

w

1

i w

2

dowolne pola wektorowe

gdzie s – dowolne pole skalarne, 

w - dowolne pole wektorowe

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

24

Analiza pól  – Rotacja pola 

wektorowego

Operację rotacji wykonuje się na polu wektorowym. Wynikiem jest pole 

wektorowe. Polskie określenie tego operatora to wirowość

[

( , , , ),

( , , , ),

( , , , )]

( )

[ , , ]

x

y

z

x

y

z

y

z

x

x

z

y

y

x

z

w

w x y z t w x y z t w x y z t

rot w

r r r

gdzie

w

w

r

y

z

w

w

r

z

x

w

w

r

x

y

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

25

Analiza pól  – Laplasjan pola 

skalarnego

( , , , )

( )

[

( )]

s

f x y z t

s

div grad s

Laplasjan jest operatorem złożonym składającym się z operatorów gradientu 

i diwergencji.

Laplasjanem działamy na pole skalarne w wyniku otrzymując inne pole skalarne:

Dla układu kartezjańskiego operator Laplace’a jest sumą drugich pochodnych 

cząstkowych pola skalarnego względem współrzędnych przestrzennych:

2

2

2

2

2

2

)

(

z

f

y

f

x

f

s

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

26

Analiza pól  – Operator Hamiltona 

z

y

x

,

,

W układzie kartezjańskim jest to symboliczny wektor, którego składowymi

są operatory różniczkowania względem zmiennych przestrzennych:

W celu uporządkowania i łatwego zapisu powyższych pojęć czasami

jest stosowany tzw. operator Hamiltona oznaczony symbolem „nabla”.

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

27

Analiza pól  – Operator Hamiltona 

Za pomocą operatora „nabla” można otrzymać wszystkie do tej pory 

zdefiniowane pojęcia:

2

( )

( )

,

( )

( )

(

)

( )

grad s

s

s

div w

w

w

s

div grad s

s

s

rot w

w

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

28

Analiza pól  – Pola potencjalne

Operacją odwrotną do różniczkowania jest całkowanie. W przypadku pola

wektorowego odpowiednikiem całki nieoznaczonej (czyli tzw. funkcji

pierwotnej) może być pojęcie potencjału. Pojęcie to można definiować

tylko dla tzw. pól potencjalnych.

Dane pole wektorowe w nazywamy polem potencjalnym, jeżeli istnieje

pole skalarne U, którego gradientem jest dane pole w. Czyli

( )

,

,

U

U

U

w

grad U

x

y

z

Potencjał danego pola (podobnie jak całka nieoznaczona) jest określony

z dokładnością do pewnej stałej addytywnej.

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

29

Analiza pól  – Pola potencjalne

Warunkiem koniecznym i wystarczającym aby dane pole wektorowe w

było polem potencjalnym jest jego bezwirowość tzn. że jego rotacja

musi być równa 0:

( )

[

( )]

0

rot w

rot grad U

Warunek ten w układzie kartezjańskim można zapisać za pomocą 3 równań:

y

y

x

x

z

z

w

w

w

w

w

w

y

z

z

x

x

y

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

30

Analiza pól  – Całkowanie

W analizie pól rozpatruje się również operacje analogiczne do całki

oznaczonej. Najważniejsze są tzw. całki objętościowe całki

powierzchniowe

Dla pola skalarnego s w zamkniętym obszarze przestrzennym Ω

definiuje się całkę objętościową jako granicę:

1

lim

( ) ( )

( )

0

n

i

i

i

j

i

n

i

i

i

i

sdV

s u V

u

V

W dowolnym układzie przestrzennym dla odpowiednio zdefiniowanego

zbioru Ω można całkę objętościową zapisać za pomocą całki potrójnej

W odpowiednich granicach. 

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

31

Analiza pól  – Całkowanie

Dla pola wektorowego w i ograniczonej powierzchni zorientowanej A
definiuje się całkę powierzchniową drugiego rodzaju jako granicę:

1

lim

( ( )

) ( )

( )

0

n

i

i

i

j

i

n

i

i

i

i

wd

w u

n A a

a

a

a

u

a

A

A

A

A

W dowolnym układzie przestrzennym dla odpowiednio zdefiniowanej

powierzchni A można całkę powierzchniową zapisać za pomocą całki 
podwójnej w odpowiednich granicach. 

background image

© Prof. Antoni Kozioł, Wydział Chemiczny Politechniki Wrocławskiej

32

Analiza pól  – Twierdzenie 

Gaussa

W szczególnym przypadku, gdy mamy dane pole wektorowe w 

zdefiniowane w zamkniętym obszarze Ω, którego brzegiem jest 

powierzchnia zorientowana na zewnątrz A całkę powierzchniową tego
pola po tej powierzchni można wyrazić za pomocą całki objętościowej.

Umożliwia to tzw. twierdzenie Gaussa:

( )

wd

div w dV

A

A