Całka powierzchniowa jest jednym z tych pojęć matematycznych, które z jednej strony wydają się abstrakcyjne, z drugiej zaś opisują bardzo konkretne zjawiska fizyczne i inżynierskie. Pozwala mierzyć nie tylko pole zakrzywionej powierzchni, ale także rozkład wielkości takich jak temperatura, gęstość masy czy strumień pola wektorowego na skomplikowanych obiektach. Jej pełne zrozumienie wymaga spojrzenia zarówno od strony formalnej analizy matematycznej, jak i przez pryzmat zastosowań w naukach przyrodniczych oraz technice.
Intuicyjna idea całki powierzchniowej
Wyobraźmy sobie cienką, zakrzywioną membranę, na przykład fragment czaszy satelitarnej, na której rozłożona jest masa o zmiennej gęstości. Chcąc obliczyć całkowitą masę, nie wystarczy znać zwykłe pole rzutowane na płaszczyznę – trzeba uwzględnić rzeczywisty kształt powierzchni oraz to, jak gęstość zmienia się w każdym punkcie. To właśnie zadanie dla całki powierzchniowej skalarnej. Z kolei, gdy analizujemy przepływ powietrza lub pola magnetycznego przez powierzchnię, interesuje nas ile wektorowego „strumienia” przechodzi przez daną powierzchnię: tu pojawia się całka powierzchniowa wektorowa, zwana strumieniową.
W klasycznej edukacji matematycznej zaczyna się od całek oznaczonych funkcji jednej zmiennej, następnie przechodzi do całek podwójnych i potrójnych. Całka powierzchniowa jest naturalnym kolejnym krokiem: zamiast całkować po przedziale czy obszarze płaskim, całkujemy po zakrzywionej powierzchni zanurzonej w przestrzeni trójwymiarowej. Jej zrozumienie wymaga połączenia kilku wątków: geometrii analitycznej, rachunku różniczkowego oraz algebry liniowej.
Formalna definicja i klasyfikacja całek powierzchniowych
Formalne zdefiniowanie całki powierzchniowej opiera się na pojęciu parametryzacji powierzchni. Rozpatrzmy gładką powierzchnię S w przestrzeni trójwymiarowej R³, którą można opisać za pomocą odwzorowania
r(u, v) = (x(u, v), y(u, v), z(u, v))
gdzie (u, v) należy do pewnego obszaru D w płaszczyźnie. Parę (u, v) traktujemy jako współrzędne wewnętrzne, a funkcja r opisuje sposób „rozciągnięcia” płaskiego obszaru D na zakrzywioną powierzchnię S. Dzięki temu przenosimy całkowanie z trudnej, trójwymiarowej powierzchni na znacznie prostszy obszar w R².
Wyróżnia się dwa główne typy całek powierzchniowych:
- całka powierzchniowa funkcji skalarnej po powierzchni, oznaczana najczęściej ∫∫S f dS,
- całka powierzchniowa pola wektorowego, tzw. strumieniowa, oznaczana ∫∫S F · n dS.
W pierwszym przypadku f jest funkcją przypisującą każdemu punktowi na powierzchni wartość rzeczywistą (np. temperaturę). W drugim F jest polem wektorowym (np. prędkością przepływu płynu), a n jest wektorem normalnym do powierzchni, wyznaczającym orientację. Iloczyn skalarny F · n opisuje część pola skierowaną „przez” powierzchnię.
Całka powierzchniowa funkcji skalarnej
Dla funkcji skalarnej f zdefiniowanej na powierzchni S, całka powierzchniowa ma postać
∫∫S f(p) dS,
gdzie p należy do S. Definicja ścisła wynika z granicy sum Riemanna. Dzielimy powierzchnię na bardzo małe fragmenty Sk, wybieramy w każdym punkty pk i sumujemy wartości funkcji pomnożone przez pola tych fragmentów:
∑ f(pk) ΔSk.
Całka to granica tej sumy, gdy maksymalna średnica fragmentów dąży do zera. W praktycznych obliczeniach nie korzysta się bezpośrednio z tej definicji granicznej, lecz przechodzi się do współrzędnych (u, v) poprzez parametryzację. Wówczas element powierzchni dS zastępowany jest wyrażeniem opartym na pochodnych cząstkowych ru = ∂r/∂u i rv = ∂r/∂v.
Kluczowym obiektem jest wektor
ru × rv,
który jest prostopadły do powierzchni oraz którego długość równa jest polu równoległoboku rozpiętego na wektorach stycznych w danym punkcie. Wynika stąd znany wzór:
∫∫S f dS = ∫∫D f(r(u, v)) |ru × rv| dudv.
Funkcja podcałkowa jest więc przekształcona składnikiem geometrycznym, który uwzględnia zakrzywienie i „rozciągnięcie” powierzchni względem płaskiej dziedziny D. Ten czynnik ma charakter czysto geometryczny i jest fundamentem w rozumieniu, że całkowanie po powierzchni to w istocie całkowanie po odpowiednio zważonej dziedzinie płaskiej.
Całka powierzchniowa pola wektorowego (strumieniowa)
Całka powierzchniowa pola wektorowego F po powierzchni S ma postać
∫∫S F · n dS,
gdzie n jest jednostkowym wektorem normalnym do S. Ta całka ma bardzo istotne znaczenie fizyczne: opisuje łączny strumień pola przez badaną powierzchnię. Interpretując F jako prędkość przepływu płynu, a n jako kierunek „przebicia” przez powierzchnię, otrzymujemy ilość płynu przepływającą w jednostce czasu przez S.
Podobnie jak w przypadku funkcji skalarnej, stosujemy parametryzację r(u, v). Wektor normalny niejednostkowy możemy zapisać jako
N = ru × rv,
a wtedy element wektorowy powierzchni dSwektor = N dudv. Całka strumieniowa przyjmuje wówczas formę:
∫∫S F · n dS = ∫∫D F(r(u, v)) · (ru × rv) dudv.
W tym zapisie nie ma już osobno czynnika długości i wektora jednostkowego – ich rola zawarta jest w jednym wektorze ru × rv. W zależności od wyboru orientacji powierzchni, wektor ten może być skierowany „na zewnątrz” lub „do wewnątrz” danego obiektu. Zmiana orientacji zmienia znak całki.
W fizyce i technice strumieniowe całki powierzchniowe są wszechobecne. Wystarczy wymienić podstawowe przykłady: strumień pola elektrycznego przez powierzchnię zamkniętą w prawie Gaussa, ilość energii promieniowania przechodząca przez powierzchnię w optyce, czy ilość ciepła transportowanego przez konwekcję w inżynierii cieplnej.
Rola metryk i uogólnienie na powierzchnie zakrzywione
Dotychczas rozpatrywaliśmy powierzchnie zanurzone w klasycznej przestrzeni euklidesowej R³. Jednak pojęcie całki powierzchniowej ma głębszy sens w szerszym kontekście geometrii różniczkowej i struktur metrycznych. W ogólnym przypadku rozważamy powierzchnię jako dwuwymiarową rozmaitość zadaną parametrami (u, v) wyposażoną w metrykę g, która definiuje sposób mierzenia długości i pól na tej powierzchni.
W bazie {∂/∂u, ∂/∂v} tensor metryczny ma współczynniki gij, i, j = 1, 2, a element pola dS wyraża się przez pierwiastek z wyznacznika metryki:
dS = √det(g) dudv.
Dla powierzchni zanurzonej w R³ metryka jest indukowana z metryki euklidesowej przez parametryzację r i można wykazać, że
det(g) = |ru × rv|².
Dzięki temu ogólne ujęcie metryczne zgodne jest z klasycznym. Jednak takie spojrzenie otwiera drogę do analiz na powierzchniach, które nie muszą być „wygodnie” zanurzone w R³, jak choćby sfera w ogólnej teorii względności, gdzie istotne jest zakrzywienie czasoprzestrzeni i lokalna struktura geometryczna. Całki powierzchniowe w tym kontekście służą do definiowania wielkości globalnych, takich jak całkowita krzywizna czy wielkości związane z polem grawitacyjnym.
Przykłady obliczeniowe: sfera i powierzchnia walca
Rozważmy obliczenie pola całkowitego sfery o promieniu R za pomocą całki powierzchniowej. Naturalną parametryzacją jest zastosowanie współrzędnych sferycznych: niech
r(φ, θ) = (R sinφ cosθ, R sinφ sinθ, R cosφ),
gdzie φ należy do [0, π], a θ do [0, 2π]. Wówczas
rφ = (R cosφ cosθ, R cosφ sinθ, −R sinφ),
rθ = (−R sinφ sinθ, R sinφ cosθ, 0).
Obliczamy iloczyn wektorowy:
rφ × rθ = (R² sin²φ cosθ, R² sin²φ sinθ, R² sinφ cosφ).
Jego długość wynosi
|rφ × rθ| = R² sinφ.
Chcąc policzyć pole powierzchni sfery, całkujemy funkcję f ≡ 1:
∫∫S 1 dS = ∫0π ∫02π R² sinφ dθ dφ = 4πR².
Otrzymaliśmy dobrze znaną formułę na pole powierzchni sfery. Jednak wynik ten zyskał interpretację jako szczególny przypadek całki powierzchniowej.
Dla walca o promieniu R i wysokości H możemy wybrać parametryzację bocznej powierzchni w postaci
r(θ, z) = (R cosθ, R sinθ, z),
gdzie θ należy do [0, 2π], a z do [0, H]. Otrzymujemy
rθ = (−R sinθ, R cosθ, 0),
rz = (0, 0, 1).
Iloczyn wektorowy ma postać
rθ × rz = (R cosθ, R sinθ, 0),
a jego długość to R. Całkując funkcję jednostkową:
∫∫S 1 dS = ∫0H ∫02π R dθ dz = 2πRH,
co odpowiada znanemu wzorowi na pole powierzchni bocznej walca. Nawet tak elementarne wzory geometryczne możemy więc interpretować jako rezultat zastosowania rachunku całkowego na powierzchniach.
Zastosowania w fizyce klasycznej
Pojęcie całki powierzchniowej odgrywa centralną rolę w fizyce klasycznej, szczególnie w elektromagnetyzmie, mechanice płynów i teorii potencjału. Jednym z najważniejszych zastosowań jest prawo Gaussa, które łączy strumień pola elektrycznego przez dowolną powierzchnię zamkniętą z ładunkiem elektrycznym wewnątrz. W zapisie matematycznym prawo to przyjmuje formę
∫∫S E · n dS = Qwew/ε₀,
gdzie E jest polem elektrycznym, Qwew łącznym ładunkiem wewnątrz S, a ε₀ przenikalnością elektryczną próżni. To równanie ma głębokie konsekwencje symetrii i umożliwia szybkie wyprowadzenie rozkładów pola dla prostych konfiguracji ładunków, wykorzystując odpowiednio dobrane powierzchnie całkowania.
W mechanice płynów całki powierzchniowe wyrażają ilość materii lub pędu przepływającą przez powierzchnię. Dla pola prędkości v i gęstości ρ, całkowity przepływ masy przez S w jednostce czasu opisuje się jako
∫∫S ρ v · n dS.
Całka ta pozwala formułować równania bilansu masy i pędu w postaci globalnej, a następnie, korzystając z twierdzeń przejścia do postaci lokalnej, prowadzi do równania ciągłości i równań Naviera–Stokesa. Bez aparatu całek powierzchniowych sformułowanie nowoczesnej mechaniki płynów byłoby praktycznie niemożliwe.
Twierdzenie Gaussa–Ostrogradskiego i związek z całkami objętościowymi
Jednym z najbardziej fundamentalnych wyników teorii całek powierzchniowych jest twierdzenie Gaussa–Ostrogradskiego (nazywane też twierdzeniem dywergencji). Łączy ono całkę strumieniową po powierzchni zamykającej pewien obszar z całką objętościową dywergencji pola wektorowego w tym obszarze. W najczęściej spotykanej formie twierdzenie brzmi
∫∫∂V F · n dS = ∭V div F dV,
gdzie V jest obszarem trójwymiarowym, a ∂V jego brzegiem, czyli powierzchnią zamkniętą. Intuicyjnie, lewa strona opisuje całkowity strumień pola „wypływający” z V, podczas gdy prawa strona sumuje lokalne źródła i ubytki w objętości. To globalno-lokalne przejście odgrywa fundamentalną rolę w całej fizyce polowej, od hydrodynamiki po elektromagnetyzm.
Twierdzenie Gaussa–Ostrogradskiego można interpretować jako uogólnienie jednowymiarowego twierdzenia Newtona–Leibniza. Zamiast różnicy wartości funkcji na końcach przedziału, otrzymujemy sumę efektów wewnętrznych (dywergencji) w całym obszarze, równą strumieniowi przez jego brzeg. Uogólnienia tego wynikają również w formalizmie form różniczkowych, gdzie całka po brzegu jest równa całce z pochodnej zewnętrznej po wnętrzu.
Całki powierzchniowe w elektrodynamice i optyce
W elektromagnetyzmie całki powierzchniowe występują niemal na każdym kroku. Prawo Gaussa w elektrodynamice i jego magnetostatyczny odpowiednik (opisujący brak pojedynczych ładunków magnetycznych) to najprostsze przykłady. Równania Maxwella łączą pola elektryczne i magnetyczne, a w ich globalnej formie każdemu równaniu różniczkowemu odpowiada postać całkowa, w której całki powierzchniowe i krzywoliniowe są podstawowymi składnikami.
W optyce falowej pojęcie strumienia energii promieniowania opisuje wektor Poyntinga. Całkując jego składową normalną do powierzchni, otrzymujemy łączną moc przenoszoną przez falę elektromagnetyczną przez daną powierzchnię. Z kolei w optyce geometrycznej rozkład natężenia światła na detektorach może być przybliżany przez całkowanie po ich powierzchni odpowiednio zdefiniowanego pola energii. W projektowaniu układów optycznych czy antenowych takie obliczenia są codziennością.
Co więcej, zastosowanie całek powierzchniowych w fotonice i optyce zintegrowanej obejmuje analizę propagacji fal w strukturach falowodowych, gdzie ważne jest dokładne określenie pola elektromagnetycznego na granicach różnych ośrodków. Tutaj całki powierzchniowe służą nie tylko do obliczania wielkości energetycznych, ale także do formułowania warunków brzegowych i sprzężeń między modami falowymi.
Całki powierzchniowe w mechanice continuum i inżynierii
W mechanice ośrodków ciągłych całki powierzchniowe pojawiają się np. w równaniach równowagi sił działających na elementy materiału. Napięcie mechaniczne w ciele sprężystym opisywane jest przez tensor naprężeń, a siła działająca na mały fragment powierzchni jest wynikiem działania tego tensora na wektor normalny. Całkowita siła na wybranej powierzchni S dana jest wzorem
∫∫S T n dS,
gdzie T jest tensorem naprężeń. W konstrukcjach inżynierskich, od mostów po kadłuby samolotów, analiza rozkładu naprężeń i sił jest kluczowa dla zapewnienia bezpieczeństwa i optymalizacji masy. Te zagadnienia formułuje się w języku równań różniczkowych cząstkowych z naturalnymi warunkami brzegowymi zapisanymi właśnie jako całki po powierzchniach.
W inżynierii cieplnej i procesowej całki powierzchniowe opisują także transport ciepła czy masy przez granice urządzeń, wymienników ciepła, przewodów. Wymiana ciepła przez konwekcję na powierzchni ściany może być modelowana przez całkę z lokalnej gęstości strumienia ciepła. Z kolei w aerodynamice siły nośne i oporu na skrzydłach samolotu są obliczane, w bardziej zaawansowanych modelach numerycznych, poprzez całkowanie rozkładów ciśnień oraz naprężeń lepkościowych na powierzchni skrzydła.
Obliczenia numeryczne i metody dyskretyzacji powierzchni
W realnych zastosowaniach rzadko mamy do czynienia z idealnie prostymi powierzchniami o znanej, eleganckiej parametryzacji. Geometrie struktur inżynierskich, obiektów biologicznych czy krajobrazów geofizycznych są złożone i nieregularne. Dlatego w praktyce dominują numeryczne metody obliczania całek powierzchniowych. Najczęściej stosuje się dyskretyzację powierzchni na elementy skończone, najprościej trójkąty lub czworokąty, na których lokalnie aproksymuje się zarówno powierzchnię, jak i funkcję podcałkową.
W metodzie elementów skończonych (FEM) czy objętości skończonych (FVM) integralne formy równań bilansowych są przekształcane właśnie w sumy po dyskretnych elementach powierzchniowych i objętościowych. W każdej z komórek siatki numerycznej całka powierzchniowa jest przybliżana przez odpowiednie formuły kwadratur, dobierane z uwzględnieniem wymaganej dokładności oraz gładkości funkcji. Przy skomplikowanych geometriach wykorzystuje się siatki nieliniowe, adaptacyjne oraz specjalne techniki rekonstrukcji powierzchni.
Zaawansowane oprogramowanie inżynierskie (np. do symulacji CFD czy MES) zawiera wyspecjalizowane algorytmy tworzenia siatek trójwymiarowych i integralnych formuł dla powierzchni złożonych. Rozwój tych narzędzi silnie wpływa na możliwości analizy i projektowania, pozwalając na modelowanie zjawisk, których opis analityczny jest praktycznie nieosiągalny. W efekcie, całka powierzchniowa, choć definiowana teoretycznie za pomocą granic sum Riemanna, w codziennej praktyce naukowej i przemysłowej istnieje przede wszystkim jako obiekt numeryczny.
Całki powierzchniowe w geofizyce, astronomii i naukach o Ziemi
W naukach o Ziemi i planetach całki powierzchniowe są narzędziem do opisu globalnych wielkości fizycznych. Przykładowo, globalną średnią temperaturę na powierzchni planety czy globalną sumę opadów można wyrazić jako całkę powierzchniową odpowiednich funkcji skalarnych, znormalizowaną przez całkowite pole powierzchni. W meteorologii całki po powierzchni kuli ziemskiej pozwalają definiować wskaźniki klimatyczne, obciążenia radiacyjne czy bilanse energetyczne atmosfery.
W geofizyce grawitacyjnej całki po powierzchni sfery (lub sferoidy obrotowej) służą do obliczania potencjału grawitacyjnego i pola grawitacyjnego generowanego przez rozkład mas wewnątrz planety. W astrofizyce, całki powierzchniowe odgrywają ważną rolę w opisie jasności i rozkładu promieniowania gwiazd, gdzie rozkład temperatury i własności emisyjnych na powierzchni wpływa na obserwowane widmo i całkowitą jasność.
Dodatkowo, w globalnej sejsmologii rozkłady naprężeń i deformacji w skorupie ziemskiej są analizowane również poprzez wielkości zdefiniowane jako całki po powierzchniach granicznych warstw geologicznych. Opis procesów takich jak subdukcja płyt tektonicznych, generowanie trzęsień ziemi czy przepływ ciepła przez powierzchnię Ziemi wymaga bogatego formalizmu całkowego.
Związek z formami różniczkowymi i uogólnienie pojęcia
W nowoczesnym ujęciu matematycznym całka powierzchniowa jest szczególnym przypadkiem całkowania form różniczkowych. W tym formalizmie pola skalarne i wektorowe opisuje się jako 0-, 1-, 2- lub 3-formy definiowane na rozmaitościach różniczkowych. Całkowanie formy k-tego rzędu odbywa się po k-wymiarowych podrozmaitościach. Dla powierzchni dwuwymiarowych naturalny jest więc opis w języku 2-form.
Przykładowo, strumień pola wektorowego F można interpretować jako całkę z odpowiedniej 2-formy ω, zdefiniowanej tak, aby jej działanie na parę wektorów stycznych odpowiadało iloczynowi skalarnemu F z wektorem normalnym. Wówczas całka strumieniowa po S zapisana jest jako
∫S ω.
Taki język nie tylko porządkuje wiele wzorów i upraszcza dowody twierdzeń typu Stokesa czy Gaussa, ale otwiera też drogę do naturalnych uogólnień w geometrii różniczkowej, teorii wiązek czy topologii algebraicznej. Z tej perspektywy całki powierzchniowe nie są odosobnionym konstruktem, lecz częścią spójnego systemu opisującego całkowanie w różnych wymiarach.
Znaczenie dydaktyczne i epistemologiczne
Rozumienie całek powierzchniowych jest ważnym etapem rozwoju myślenia matematycznego studentów nauk ścisłych i inżynierskich. Wymaga ono harmonijnego połączenia intuicji geometrycznej z precyzyjnym rachunkiem analitycznym. Uczy, że pojęcie całki nie jest związane jedynie z „obszarem pod wykresem funkcji”, lecz z ogólniejszym procesem sumowania nieskończenie wielu drobnych wkładów rozłożonych na obiektach o różnej wymiarowości.
Z dydaktycznego punktu widzenia niezwykle ważne jest wyćwiczenie pracy z parametryzacjami, transformacjami współrzędnych, a także z interpretacją fizyczną. Dzięki temu całka powierzchniowa przestaje być zbiorem formalnych wzorów, a staje się narzędziem do modelowania rzeczywistych zjawisk w naukach przyrodniczych i technice. Poznanie jej związku z całkami krzywoliniowymi, objętościowymi i twierdzeniami typu Stokesa pozwala też dostrzec głęboką jedność wielu pozornie odrębnych zagadnień rachunku wektorowego.
Na poziomie epistemologicznym pojęcia takie jak całka powierzchniowa pokazują, w jaki sposób matematyka rozwija uogólnione struktury, zdolne opisywać zarówno idealne konstrukcje geometryczne, jak i niezwykle złożone systemy fizyczne. To uogólnienie – od sumy po prostych przedziałach do całkowania po zakrzywionych powierzchniach – ilustruje, jak abstrakcja zwiększa moc wyjaśniającą teorii naukowych.
FAQ – najczęstsze pytania o całki powierzchniowe
Czym różni się całka powierzchniowa od podwójnej całki po obszarze w płaszczyźnie?
Podwójna całka po obszarze w płaszczyźnie dotyczy funkcji określonej na płaskim regionie R², a element pola ma postać dA = dxdy lub po transformacji |J| dudv. Całka powierzchniowa dotyczy obiektu dwuwymiarowego zanurzonego w trójwymiarowej przestrzeni, często zakrzywionego. Element powierzchni dS zależy od parametryzacji r(u, v) i ma postać |ru × rv| dudv. Dodatkowo, w wersji wektorowej uwzględnia się orientację przez wektor normalny, co pozwala mierzyć strumień pola przez powierzchnię.
Jak intuicyjnie rozumieć całkę powierzchniową pola wektorowego?
Całkę strumieniową można rozumieć jako całkowitą ilość „przepływu” pola przez powierzchnię. Wyobraź sobie pole prędkości płynu: w każdym punkcie masz wektor mówiący, jak i dokąd płyn się porusza. Wektor normalny do powierzchni określa kierunek „przelotu” przez nią. Iloczyn skalarny pola z wektorem normalnym w pojedynczym punkcie daje lokalną gęstość strumienia. Całkując ten iloczyn po całej powierzchni, sumujemy wkład z każdego małego fragmentu i otrzymujemy łączną ilość płynu (lub innej wielkości) przechodzącą w jednostce czasu.
Dlaczego parametryzacja jest konieczna przy obliczaniu całek powierzchniowych?
Parametryzacja umożliwia „spłaszczenie” zakrzywionej powierzchni do prostszego obszaru w płaszczyźnie, na którym łatwiej wykonywać rachunki. Bezpośrednie operowanie na współrzędnych trójwymiarowych jest niewygodne, ponieważ trudno jest wtedy jednoznacznie określić element pola powierzchni. Parametryzacja r(u, v) wprowadza lokalne współrzędne, w których małe fragmenty powierzchni są dobrze przybliżone przez równoległoboki rozpięte na wektorach stycznych ru i rv. Dzięki temu możemy zastosować standardowy rachunek całkowy w R² i automatycznie uwzględnić krzywiznę przez czynnik |ru × rv|.
Czy całka powierzchniowa zawsze zależy od orientacji powierzchni?
W przypadku całki powierzchniowej funkcji skalarnej ∫∫S f dS wynik nie zależy od orientacji, ponieważ integrujemy po prostu wartości funkcji pomnożone przez dodatnią miarę pola. Natomiast przy całce strumieniowej ∫∫S F · n dS wybór orientacji (czyli kierunku wektora normalnego) jest kluczowy. Zmiana orientacji na przeciwną zmienia znak całki, co odpowiada fizycznej różnicy między strumieniem „na zewnątrz” a „do wewnątrz” objętości. Dlatego w zastosowaniach fizycznych orientację trzeba zawsze określać explicite.
Jakie są praktyczne sposoby liczenia całek powierzchniowych w inżynierii?
W zastosowaniach inżynierskich całki powierzchniowe zwykle oblicza się numerycznie. Powierzchnię dzieli się na elementy (np. trójkąty), dla których znana jest lokalna parametryzacja lub przybliżenie geometryczne. Następnie stosuje się formuły kwadratur, wybierając kilka punktów w każdym elemencie i odpowiednie wagi, by aproksymować wartość całki. W metodzie elementów skończonych i objętości skończonych te procedury są zautomatyzowane w oprogramowaniu, które generuje siatkę, wyznacza wektory normalne, uwzględnia zakrzywienie i oblicza całki w ramach globalnego rozwiązania równań fizycznych.
