Kwestia Stefana

W matematyce i jej zastosowaniach, zwłaszcza do przejść fazowych w materii, problem Stefana jest szczególnym rodzajem problemu wartości brzegowych dla układu równań różniczkowych cząstkowych (PDE), w którym granica między fazami może przesuwać się w czasie. Klasyczny problem Stefana ma na celu opisanie ewolucji granicy między dwiema fazami materiału przechodzącego przemianę fazową , na przykład topnienie ciała stałego, takiego jak lód w wodę . Osiąga się to poprzez rozwiązanie równań ciepła w obu obszarach, z zastrzeżeniem danych warunków brzegowych i początkowych. Na styku faz (w klasycznym problemie) temperatura jest ustawiana na temperaturę przemiany fazowej. Aby zamknąć system matematyczny, wymagane jest dalsze równanie, warunek Stefana . Jest to bilans energetyczny, który określa położenie poruszającego się interfejsu. Zauważ, że ta ewoluująca granica jest nieznaną (hiper-)powierzchnią ; stąd problemy Stefana są przykładami problemów z brzegami swobodnymi .

Analogiczne problemy pojawiają się na przykład w badaniu przepływu mediów porowatych, matematycznych finansów i wzrostu kryształów z roztworów monomerów.

Uwaga historyczna

Problem został nazwany na cześć Josefa Stefana (Jožefa Stefana), słoweńskiego fizyka , który około 1890 roku przedstawił ogólną klasę takich problemów w serii czterech artykułów dotyczących zamarzania gruntu i tworzenia się lodu morskiego . Jednak około 60 lat wcześniej, w 1831 roku, podobny problem, dotyczący formowania się skorupy ziemskiej, badali Lamé i Clapeyron . Problem Stefana dopuszcza rozwiązanie podobieństwa , często nazywane Neumannem rozwiązanie, które rzekomo zostało przedstawione w serii wykładów na początku lat sześćdziesiątych XIX wieku.

Obszerny opis historii problemów Stefana można znaleźć u Rubinsteina.

Przesłanki opisu matematycznego

Z matematycznego punktu widzenia fazy są jedynie regionami, w których rozwiązania bazowego PDE są ciągłe i różniczkowalne aż do rzędu PDE. W problemach fizycznych takie rozwiązania reprezentują właściwości ośrodka dla każdej fazy. Ruchome granice (lub interfejsy ) to nieskończenie cienkie powierzchnie , które oddzielają sąsiednie fazy; w związku z tym rozwiązania podstawowego środowiska PDE i jego pochodnych mogą wykazywać nieciągłości między interfejsami.

Bazowe PDE nie są ważne na interfejsach zmiany fazy; dlatego do uzyskania domknięcia potrzebny jest dodatkowy warunek — warunek Stefana . Warunek Stefana wyraża lokalną prędkość poruszającej się granicy jako funkcję wielkości ocenianych po obu stronach granicy faz i jest zwykle wyprowadzany z ograniczenia fizycznego. Na przykład w problemach wymiany ciepła ze zmianą fazy zasada zachowania energii nakazuje, aby nieciągłość strumienia ciepła na granicy była uwzględniana przez szybkość ciepła utajonego (które jest proporcjonalne do lokalnej prędkości interfejsu).

Regularność równania została zbadana głównie przez Luisa Caffarelliego i udoskonalona przez prace Alessio Figalli , Xaviera Ros-Otona i Joaquima Serry

Sformułowanie matematyczne

Jednowymiarowy jednofazowy problem Stefana

Jednofazowy problem Stefana opiera się na założeniu, że można pominąć jedną z faz materiału. Zwykle osiąga się to zakładając, że faza ma temperaturę zmiany fazy, a zatem wszelkie odchylenia od niej prowadzą do zmiany fazy. Jest to matematycznie wygodne przybliżenie, które upraszcza analizę, jednocześnie demonstrując podstawowe idee stojące za procesem. Kolejnym standardowym uproszczeniem jest praca w bezwymiarowym , tak że temperatura interfejsie może być ustawiona na zero, a wartości pola dalekiego na lub .

jednowymiarowy blok lodu początkowo w ; . Najbardziej znana postać problemu Stefana polega na topnieniu przez narzuconą stałą temperaturę na lewej granicy, pozostawiając region zajmowane przez wodę. Stopiona głębokość, oznaczona przez jest nieznaną funkcją czasu. Problem Stefana jest zdefiniowany przez

  • Równanie ciepła:
  • Stała temperatura, powyżej temperatury topnienia, na lewej granicy:
  • Interfejs w temperaturze topnienia jest ustawiony na
  • Warunek Stefana: stosunkiem utajonego do określonego ciepła jawnego (gdzie specyficzna wskazuje, że jest podzielona przez masę). Zauważ, że ta definicja wynika naturalnie z braku wymiarowania i jest używana w wielu tekstach, jednak może być również zdefiniowana jako odwrotność tego .
  • Początkowy rozkład temperatury:
  • Początkowa głębokość stopionej bryły lodu:
Rozwiązanie Neumanna, otrzymane przy użyciu zmiennych samopodobnych, wskazuje, że położenie granicy jest określone przez gdzie spełnia transcendentalne równanie
Temperatura w cieczy jest wtedy dana przez

Aplikacje

Oprócz modelowania topnienia ciał stałych, problem Stefana jest również używany jako model asymptotycznego zachowania się (w czasie) bardziej złożonych problemów. Na przykład Pego używa dopasowanych rozszerzeń asymptotycznych, aby udowodnić, że rozwiązania Cahna-Hilliarda dla problemów separacji faz zachowują się jak rozwiązania nieliniowego problemu Stefana w pośredniej skali czasu. Ponadto rozwiązanie równania Cahna-Hilliarda dla mieszaniny binarnej jest rozsądnie porównywalne z rozwiązaniem problemu Stefana. W tym porównaniu problem Stefana został rozwiązany przy użyciu metody śledzenia przodu, ruchomej siatki z jednorodnymi warunkami brzegowymi Neumanna na zewnętrznej granicy. Problemy Stefana można również zastosować do opisu przemian fazowych innych niż ciało stałe-ciecz lub płyn-ciecz.

Problem Stefana ma również bogatą teorię odwrotną; w takich problemach głębokość spotkania (lub krzywa lub hiperpowierzchnia ) s jest znanym układem odniesienia, a problemem jest znalezienie u lub f .

Zaawansowane formy problemu Stefana

Klasyczny problem Stefana dotyczy materiałów stacjonarnych o stałych właściwościach termofizycznych (zwykle niezależnie od fazy), stałej temperaturze przemiany fazowej oraz, w powyższym przykładzie, chwilowej zmianie temperatury początkowej na wyraźną wartość na granicy. W praktyce właściwości termiczne mogą się różnić, a szczególnie zawsze, gdy zmienia się faza. Skok gęstości przy zmianie fazy indukuje płynny ruch: wypadkowa energia kinetyczna nie figuruje w standardowym bilansie energetycznym. W przypadku chwilowego przełącznika temperatury początkowa prędkość płynu jest nieskończona, co skutkuje początkową nieskończoną energią kinetyczną. W rzeczywistości warstwa cieczy jest często w ruchu, co wymaga adwekcji lub konwekcji w równaniu ciepła . Temperatura topnienia może zmieniać się w zależności od rozmiaru, krzywizny lub prędkości interfejsu. Niemożliwe jest natychmiastowe przełączanie temperatur, a następnie trudne jest utrzymanie dokładnie ustalonej temperatury granicznej. Co więcej, w nanoskali temperatura może nawet nie być zgodna z prawem Fouriera.

Wiele z tych problemów zostało rozwiązanych w ostatnich latach w różnych zastosowaniach fizycznych. W przypadku krzepnięcia przechłodzonych stopów analizę, w której temperatura przemiany fazowej zależy od prędkości międzyfazowej, można znaleźć w Font i in . Modeluje się krzepnięcie w nanoskali ze zmienną temperaturą przemiany fazowej i efektami energii / gęstości. Zbadano krzepnięcie z przepływem w kanale w kontekście lawy i mikrokanalików lub ze swobodną powierzchnią w kontekście zamarzania wody nad warstwą lodu . Ogólny model obejmujący różne właściwości w każdej fazie, zmienną temperaturę przemiany fazowej i równania ciepła oparte na prawie Fouriera lub równaniu Guyera-Krumhansla jest analizowany w.

Zobacz też

Notatki

odniesienia historyczne

Referencje naukowe i ogólne

Linki zewnętrzne