In this work, we applied the Markov chain Monte Carlo (MCMC) method for the estimation of parameters appearing in the Pennes’ formulation of the bioheat transfer equation. The inverse problem of parameter estimation was solved with the simulated transient temperature measurements. A one-dimensional (1D) test case was used to explore the capabilities of using the MCMC method in bioheat transfer problems, specifically for the detection of skin tumors by using surface temperature measurements. The analysis of the sensitivity coefficients was performed in order to examine linear dependence and low sensitivity of the model parameters. The solution of the direct problem was verified with a commercial code. The results obtained in this work show the ability of using inverse heat transfer analysis for the detection of skin tumors.
W artykule przedstawiono wyniki obliczeń przeprowadzonych dla mokrej chłodni kominowej o ciąg naturalnym ze zrzutem spalin w chłodni. Obliczenia przeprowadzono z wykorzystaniem kodu numerycznej mechaniki płynów ANSYS Fluent 13.0. Celem obliczeń było sprawdzenie możliwości omywania powłoki chłodni przez wprowadzane spaliny, które mogą powodować korozję powłoki chłodni. Przeprowadzone obliczenia wykazały, że przy pogodzie bezwietrznej i o słabym wietrze omywanie powłoki przez spaliny nie występuje.
Praca dotyczy optymalizacji kształtu ciała sprężystego i przewodzącego ciepło. Ciało to wymienia ciepło z otoczeniem na drodze promieniowania i konwekcji, może być wieło-spójne a jego powierzchnie mogą się wzajemnie opromieniowywać. Zadanie promieniowania uwzględnienia występowanie stref zacienionych. W zadaniu termospręźystości uwzględnia się tytko wpływ poła temperatury na połę naprężeń. Zarówno zadanie przewodzenia ciepła, sprzężone z rozwiązaniem problemu promieniowania, jak i zadanie sprężystości są rozwiązywane metodą elementów brzegowych. Do rozwiązywania zadania optymalizacji kształtu użyto algorytmu ewolucyjnego. Łatwość modyfikacji siatek w MEB pozwała na prostą automatyzację obliczeń ewolucyjnych. Sama procedura optymalizacyjna cechuje się dużą elastycznością i odpornością na obecność łokałnych ekstremów.
EN
Shape optimization of heat conducting, elastic bodies subjected to thermal and standard loads is considered. Interaction of stress and temperature fields is modelled using the formulation of uncoupled steady state thermoelasticity. The presence of heat radiation with mutual irradiation of the boundaries and the presence of shadow zones is taken into account. To find the optimal shape evolutionary algorithm is used. The boundary element method is applied to discretize the thermoelasticity, conduction and radiation problems.
Two applications of the bonduary element technique to radiative analysis are discussed: radiation in participating media and the presence of concave, self irradiating cavities in heat conducting bodies. When compared with classic Hottel's zoning approach, BEM leads to shorter computing times. This is achieved because the most time consuming procedure of Hottel's approach, volume integration, can be avoided. Additionally, instead of multiple integrals only a single integration ought to be performed. The paper presents a general, self adaptive quadrature capable of evaluating integrals over the elements located close to the observation point. The paper shows also an original technique of coupling heat radiation and conduction. Both problems are solved by using BEM. Static condensation of linear unknowns before entering the iterative loop, makes this approach very efficient. Numerical examples are included.
PL
Rozpatrzono dwa zastosowania metody elementów brzegowych (MEB) w zagadnieniach promieniowania ciepła: transmisja promieniowania w ośrodku optycznie czynnym i przewodzenia ciepła w ciałach zawierających promieniujące wnęki. Czas obliczeń przy zastosowaniu MEB jest krótszy niż w klasycznej metodzie strefowej Hottela. W MEB unika się całkowania po obiętości, czyli najbardziej czasochłonnego etapu metody Hottela. Dodatkowo, w miejsce całek wielokrotnych obliczanych w metodzie strefowej, w MEB oblicza się całki jednokrotne, W artykule zaprezentowano ogólną, adaptacyjną klawiaturę do obliczeń całek po elementach leżących blisko punktu obserwacji. Przedstawiono także oryginalną procedurę sprzęgania zadań promieniowania we wnękach i przewodzenia w ściankach tworzących tę wnękę. Oba problemy cząstkowe rozwiązywane są metodą elementów brzegowych. Przed rozpoczęciem procesu iteracyjnego rozwiązywania równań bilansu ciepła, eliminowane są liniowe niewiadome. Wykorzystuje przy tym technikę statystycznej kondensacji. Procedura taka prowadzi do znacznego skrócenia czasu obliczeń. Artykuł zawiera przykłady numeryczne.
JavaScript jest wyłączony w Twojej przeglądarce internetowej. Włącz go, a następnie odśwież stronę, aby móc w pełni z niej korzystać.