In this paper the one-dimensional two-phase Stefan problem is studied analytically leading to a system of non-linear Volterra-integral-equations describing the heat distribution in each phase. For this the unified transform method has been employed which provides a method via a global relation, by which these problems can be solved using integral representations. To do this, the underlying partial differential equation is rewritten into a certain divergence form, which enables to treat the boundary values as part of the integrals. Classical analytical methods fail in the case of the Stefan problem due to the moving interface. From the resulting non-linear integro-differential equations the one for the position of the phase change can be solved in a first step. This is done numerically using a fix-point iteration and spline interpolation. Once obtained, the temperature distribution in both phases is generated from their integral representation.
In direct numerical simulations (DNS) of turbulent Couette flow, the observation has been made that the long streamwise rolls increase in length with the Reynolds number (Lee & Moser, J. Fluid Mech., vol. 842, 2018, pp. 128–145). To understand this, we employ both linear stability theory and its extension to resolvent analysis. For this, we emphasise the high Reynolds number ( $Re \rightarrow \infty$ ) and small streamwise wavenumbers ( $\alpha \rightarrow 0$ ) limit, imposing the distinguished limit $Re_{\alpha }=Re \, \alpha = O(1)$ . We find that in case of linear stability theory, $Re_{\alpha }$ acts as a global invariant in the resulting eigenvalue problem, while in case of resolvent analysis, $Re_{\alpha }$ acts as a local invariant in the behaviour of the energy of the system characterised through the first singular value $\sigma _1$ of the resolvent operator within the investigated asymptotic limit. In order to obtain constant streamwise structures for increasing Reynolds numbers, the respective streamwise wavenumber has to decrease, which verifies the observations from DNS studies of an increasing length of the streamwise structures with the Reynolds number. In linear stability theory, a parameter reduction is achieved for the above asymptotic limit, resulting in the modified Orr–Sommerfeld and Squire equations being dependent only on $Re_{\alpha }$ . The behaviour of both the coherent structures obtained from linear stability theory and resolvent analysis are compared with each other and show similar behaviours over $Re_{\alpha }$ .
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.