Analytical-Numerical Treatment of the One-Phase Stefan Problem with Constant Applied Heat Flux

Author(s):  
Otto G. Ruehr
Filomat ◽  
2018 ◽  
Vol 32 (3) ◽  
pp. 985-990
Author(s):  
Merey Sarsengeldin ◽  
Stanislav Kharin ◽  
Samat Kassabek ◽  
Zamanbek Mukambetkazin

Exact solution of inverse one phase Stefan problem is represented in the form of linear combination of integral error functions. Heat flux function is reconstructed and coefficients of solution function are found exactly. Test problem was considered for engineering purposes and it was shown that by collocation method the error for three points does not exceed 0:01%. Error estimate was calculated by maximum principle.


2006 ◽  
Vol 5 (1) ◽  
pp. 03 ◽  
Author(s):  
A. Boucíguez ◽  
R. Lozano ◽  
M. A. Lara

The one phase Stefan problem in a semi - infinite slab with heat flux boundary condition  proportional  to  t½   and  with  constant  temperature  boundary condition are presented here. In these two cases the exact solution exists, the relation  between  the  two  boundary  conditions  is  presented  here,  and  the equivalence between the two problems is demostrated.


2020 ◽  
Vol 20 (2) ◽  
pp. 437-458 ◽  
Author(s):  
Félix del Teso ◽  
Jørgen Endal ◽  
Juan Luis Vázquez

AbstractThe classical Stefan problem is one of the most studied free boundary problems of evolution type. Recently, there has been interest in treating the corresponding free boundary problem with nonlocal diffusion. We start the paper by reviewing the main properties of the classical problem that are of interest to us. Then we introduce the fractional Stefan problem and develop the basic theory. After that we center our attention on selfsimilar solutions, their properties and consequences. We first discuss the results of the one-phase fractional Stefan problem, which have recently been studied by the authors. Finally, we address the theory of the two-phase fractional Stefan problem, which contains the main original contributions of this paper. Rigorous numerical studies support our results and claims.


1971 ◽  
Vol 26 (1) ◽  
pp. 10-17 ◽  
Author(s):  
A. R. Allnatt

AbstractA kinetic equation is derived for the singlet distribution function for a heavy impurity in a lattice of lighter atoms in a temperature gradient. In the one dimensional case the equation can be solved to find formal expressions for the jump probability and hence the heat of transport, q*. for a single vacancy jump of the impurity, q* is the sum of the enthalpy of activation, a term involving only averaging in an equilibrium ensemble, and two non-equilibrium terms in­volving time correlation functions. The most important non-equilibrium term concerns the cor­relation between the force on the impurity and a microscopic heat flux. A plausible extension to three dimensions is suggested and the relation to earlier isothermal and non-isothermal theories is indicated


2021 ◽  
Author(s):  
Xueyun Wang ◽  
Xueqiao Xu ◽  
Philip B Snyder ◽  
Zeyu Li

Abstract The BOUT++ six-field turbulence code is used to simulate the ITER 11.5MA hybrid scenario and a brief comparison is made among ITER baseline, hybrid and steady-state operation (SSO) scenarios. Peeling-ballooning instabilities with different toroidal mode numbers dominate in different scenarios and consequently yield different types of ELMs. The energy loss fractions (ΔWped/Wped) caused by unmitigated ELMs in the baseline and hybrid scenarios are large (~2%) while the one in the SSO scenario is dramatically smaller (~1%), which are consistent with the features of type-I ELMs and grassy ELMs respectively. The intra ELM divertor heat flux width in the three scenarios given by the simulations is larger than the estimations for inter ELM phase based on Goldston’s heuristic drift model. The toroidal gap edge melting limit of tungsten monoblocks of divertor targets imposes constraints on ELM energy loss, giving that the ELM energy loss fraction should be smaller than 0.4%, 1.0%, and 1.2% for ITER baseline, hybrid and SSO scenarios, correspondingly. The simulation shows that only the SSO scenario with grassy ELMs may satisfy the constraint.


2018 ◽  
Vol 21 (4) ◽  
pp. 901-918 ◽  
Author(s):  
Sabrina Roscani ◽  
Domingo Tarzia

Abstract A one-dimensional fractional one-phase Stefan problem with a temperature boundary condition at the fixed face is considered by using the Riemann–Liouville derivative. This formulation is more convenient than the one given in Roscani and Santillan (Fract. Calc. Appl. Anal., 16, No 4 (2013), 802–815) and Tarzia and Ceretani (Fract. Calc. Appl. Anal., 20, No 2 (2017), 399–421), because it allows us to work with Green’s identities (which does not apply when Caputo derivatives are considered). As a main result, an integral relationship between the temperature and the free boundary is obtained which is equivalent to the fractional Stefan condition. Moreover, an exact solution of similarity type expressed in terms of Wright functions is also given.


2009 ◽  
Vol 20 (2) ◽  
pp. 187-214 ◽  
Author(s):  
WAN CHEN ◽  
MICHAEL J. WARD

The dynamics and oscillatory instabilities of multi-spike solutions to the one-dimensional Gray-Scott reaction–diffusion system on a finite domain are studied in a particular parameter regime. In this parameter regime, a formal singular perturbation method is used to derive a novel ODE–PDE Stefan problem, which determines the dynamics of a collection of spikes for a multi-spike pattern. This Stefan problem has moving Dirac source terms concentrated at the spike locations. For a certain subrange of the parameters, this Stefan problem is quasi-steady and an explicit set of differential-algebraic equations characterizing the spike dynamics is derived. By analysing a nonlocal eigenvalue problem, it is found that this multi-spike quasi-equilibrium solution can undergo a Hopf bifurcation leading to oscillations in the spike amplitudes on an O(1) time scale. In another subrange of the parameters, the spike motion is not quasi-steady and the full Stefan problem is solved numerically by using an appropriate discretization of the Dirac source terms. These numerical computations, together with a linearization of the Stefan problem, show that the spike layers can undergo a drift instability arising from a Hopf bifurcation. This instability leads to a time-dependent oscillatory behaviour in the spike locations.


Sign in / Sign up

Export Citation Format

Share Document