Simulation of the heat accumulator operation of the internal combustion engine preheating system

We considered the heat accumulator with a phase-transfer heat-accumulating material, which serves for pre-heating of the car’s internal combustion engine. Simulation of the heat accumulator operation allows to build calculated graphs of temperature change of the heat-accumulating material in time, and afterwards to determine the charging time of the heat accumulator depending on its design features, thus, by modelling the most optimal design solution. We performed numerical computations of the system engine – circulating fluid – heat storage material – environment in two stages. In the first stage, we calculated the parameters of thermal resistance in the engine system and pipe manifold for different coolant temperatures according to the method of finite volume in the CFD system. In the second stage the problem was solved numerically by the method of equivalent thermal circuit. We carried out phase transition simulation using the Stefan condition, based on the thermal balance for the phase separation surface. We constructed numerical algorithmic models for calculations of temperature change of heat-accumulating material in time. Such calculations allowed determining the optimal number of U-shaped tubes based on which we proposed the heat accumulator design. We manufactured the heat accumulator, tested, and proved its efficiency and positive effect on the engine warm-up time and the passenger compartment.

A free boundary problem with a Stefan condition for a ratio-dependent predator-prey model

<abstract><p>In this paper we study a ratio-dependent predator-prey model with a free boundary caused by predator-prey interaction over a one dimensional habitat. We study the long time behaviors of the two species and prove a spreading-vanishing dichotomy; namely, as $ t $ goes to infinity, both prey and predator successfully spread to the whole space and survive in the new environment, or they spread within a bounded area and eventually die out. The criteria governing spreading and vanishing are obtained. Finally, when spreading occurs we provide some estimates to the asymptotic spreading speed of the moving boundary $ h(t) $.</p></abstract>

Revisiting the Fisher–Kolmogorov–Petrovsky–Piskunov equation to interpret the spreading–extinction dichotomy

The Fisher–Kolmogorov–Petrovsky–Piskunov model, also known as the Fisher–KPP model, supports travelling wave solutions that are successfully used to model numerous invasive phenomena with applications in biology, ecology and combustion theory. However, there are certain phenomena that the Fisher–KPP model cannot replicate, such as the extinction of invasive populations. The Fisher–Stefan model is an adaptation of the Fisher–KPP model to include a moving boundary whose evolution is governed by a Stefan condition. The Fisher–Stefan model also supports travelling wave solutions; however, a key additional feature of the Fisher–Stefan model is that it is able to simulate population extinction, giving rise to a spreading–extinction dichotomy . In this work, we revisit travelling wave solutions of the Fisher–KPP model and show that these results provide new insight into travelling wave solutions of the Fisher–Stefan model and the spreading–extinction dichotomy. Using a combination of phase plane analysis, perturbation analysis and linearization, we establish a concrete relationship between travelling wave solutions of the Fisher–Stefan model and often-neglected travelling wave solutions of the Fisher–KPP model. Furthermore, we give closed-form approximate expressions for the shape of the travelling wave solutions of the Fisher–Stefan model in the limit of slow travelling wave speeds, c ≪1.

Revisiting the Fisher-KPP equation to interpret the spreading-extinction dichotomy

The Fisher-KPP model supports travelling wave solutions that are successfully used to model numerous invasive phenomena with applications in biology, ecology, and combustion theory. However, there are certain phenomena that the Fisher-KPP model cannot replicate, such as the extinction of invasive populations. The Fisher-Stefan model is an adaptation of the Fisher-KPP model to include a moving boundary whose evolution is governed by a Stefan condition. The Fisher-Stefan model also supports travelling wave solutions; however, a key additional feature of the Fisher-Stefan model is that it is able to simulate population extinction, giving rise to aspreading-extinction dichotomy.In this work, we revisit travelling wave solutions of the Fisher-KPP model and show that these results provide new insight into travelling wave solutions of the Fisher-Stefan model and the spreading-extinction dichotomy. Using a combination of phase plane analysis, perturbation analysis and linearisation, we establish a concrete relationship between travelling wave solutions of the Fisher-Stefan model and often-neglected travelling wave solutions of the Fisher-KPP model. Furthermore, we give closed-form approximate expressions for the shape of the travelling wave solutions of the Fisher-Stefan model in the limit of slow travelling wave speeds,c≪ 1.

An integral relationship for a fractional one-phase Stefan problem

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.

THE QUENCHING OF SOLUTIONS OF A REACTION–DIFFUSION EQUATION WITH FREE BOUNDARIES

This paper concerns the quenching phenomena of a reaction–diffusion equation $u_{t}=u_{xx}+1/(1-u)$ in a one dimensional varying domain $[g(t),h(t)]$, where $g(t)$ and $h(t)$ are two free boundaries evolving by a Stefan condition. We prove that all solutions will quench regardless of the choice of initial data, and we also show that the quenching set is a compact subset of the initial occupying domain and that the two free boundaries remain bounded.

A three-phase free boundary problem with melting ice and dissolving gas

We develop a mathematical model for a three-phase free boundary problem in one dimension that involves interactions between gas, water and ice. The dynamics are driven by melting of the ice layer, while the pressurized gas also dissolves within the meltwater. The model incorporates the Stefan condition at the water–ice interface along with Henry's law for dissolution of gas at the gas–water interface. We employ a quasi-steady approximation for the phase temperatures and then derive a series solution for the interface positions. A non-standard feature of the model is an integral free boundary condition that arises from mass conservation owing to changes in gas density at the gas–water interface, which makes the problem non-self-adjoint. We derive a two-scale asymptotic series solution for the dissolved gas concentration, which because of the non-self-adjointness gives rise to a Fourier series expansion in eigenfunctions that do not satisfy the usual orthogonality conditions. Numerical simulations of the original governing equations are used to validate series approximations.

Calculation of Phase-Change Boundary Position in Continuous Casting

The problem of determination of the phase-change boundary position at the mathematical modeling of continuous ingot temperature field is considered. The description of the heat transfer process takes into account the dependence of the thermal physical characteristics on the temperature, so that the mathematical model is based on the nonlinear partial differential equations. The boundary position between liquid and solid phase is given by the temperatures equality condition and the Stefan condition for the two-dimensional case. The new method of calculation of the phase-change boundary position is proposed. This method based on the finite-differences with using explicit schemes and on the iteration method of solving of non-linear system equations. The proposed method of calculation is many times faster than the real time. So that it amenable to be used for model predictive control of continuous semifinished product solidification.

Ice ripple formation at large Reynolds numbers

A free-surface-induced morphological instability is studied in the laminar regime at large Reynolds numbers ($\mathit{Re}= 1\text{{\ndash}} 1{0}^{3} $) and on sub-horizontal walls ($\vartheta \lt 3{0}^{\ensuremath{\circ} } $). We analytically and numerically develop the stability analysis of an inclined melting–freezing interface bounding a free-surface laminar flow. The complete solution of both the linearized flow field and the heat conservation equations allows the exact derivation of the upper and lower temperature gradients at the interface, as required by the Stefan condition, from which the dispersion relationship is obtained. The eigenstructure is obtained and discussed. Free-surface dynamics appears to be crucial for the triggering of upstream propagating ice ripples, which grow at the liquid–solid interface. The kinematic and the dynamic conditions play a key role in controlling the formation of the free-surface fluctuations; these latter induce a streamline distortion with an increment of the wall-normal velocities and a destabilizing phase shift in the net heat transfer to the interface. Three-dimensional effects appear to be crucial at high Reynolds numbers. The role of inertia forces, vorticity, and thermal boundary conditions are also discussed.