Search results

This paper focuses on the unconditionally optimal error estimates of a linearized second-order scheme for a nonlocal nonlinear parabolic problem. The first step of the scheme is based on Crank–Nicholson method while the second step is the second-order BDF method.

A rigorous error analysis is done, and optimal L² error estimates are derived using the error splitting technique. Some numerical simulations are presented to confirm the study’s theoretical analysis.

Optimal L² error estimates and energy norm.

The goal of this research article is to present and establish the unconditionally optimal error estimates of a linearized second-order BDF finite element scheme for the reaction-diffusion problem. An optimal error estimate for the proposed methods is derived by using the temporal-spatial error splitting techniques, which split the error between the exact solution and the numerical solution into two parts, that is, the temporal error and the spatial error. Since the spatial error is not dependent on the time step, the boundedness of the numerical solution in L∞-norm follows an inverse inequality immediately without any restriction on the grid mesh.

The purpose of this paper is to provide an insight and to solve numerically the identification of timewise terms and free boundaries coefficient appearing in the heat equation from over-determination conditions.

The formulated coefficient identification problem is inverse and ill-posed, and therefore, to obtain a stable solution, a nonlinear Tikhonov regularization least-squares approach is used. For the numerical discretization, the finite difference method combined with a regularized nonlinear minimization is performed using the MATLAB subroutine lsqnonlin.

The numerical results presented for two examples show the efficiency of the computational method and the accuracy and stability of the numerical solution even in the presence of noise in the input data.

The mathematical formulation is restricted to identify coefficients in unknown components dependent on time, and this may be considered as a research limitation. However, there is no research implication to overcome this, as the known input data is also limited to single temperature in heat equation with Stefan conditions, and the first- and second-order heat moments measurements at a particular time location.

As noisy data are inverted, the study models real situations in which practical measurements are inherently contaminated with noise.

The identification of the timewise terms and free boundaries will be of great interest in the heat transfer community and related fluid flow applications.

The current investigation advances previous studies, which assumed that the coefficient multiplying the lower order temperature term depends on time. The knowledge of this physical property coefficient is very important in heat transfer and fluid flow. The originality lies in the insight gained by performing for the numerical simulations of inversion to find the timewise terms and free boundaries coefficient dependent on time in the heat equation from noisy measurements.

The purpose of the study is to solve numerically the inverse problem of determining the time-dependent convection coefficient and the free boundary, along with the temperature in the two-dimensional convection-diffusion equation with initial and boundary conditions supplemented by non-local integral observations. From the literature, there is already known that this inverse problem has a unique solution. However, the problem is still ill-posed by being unstable to noise in the input data.

For the numerical discretization, this paper applies the alternating direction explicit finite-difference method along with the Tikhonov regularization to find a stable and accurate numerical solution. The resulting nonlinear minimization problem is solved computationally using the MATLAB routine lsqnonlin. Both exact and numerically simulated noisy input data are inverted.

The numerical results demonstrate that accurate and stable solutions are obtained.

The inverse problem presented in this paper was already showed to be locally uniquely solvable, but no numerical solution has been realized so far; hence, the main originality of this work is to attempt this task.

This paper aims to review the existing bond-based peridynamic (PD) and state-based PD heat conduction models, and further propose a refined bond-based PD thermal conduction model by using the PD differential operator.

The general refined bond-based PD is established by replacing the local spatial derivatives in the classical heat conduction equations with their corresponding nonlocal integral forms obtained by the PD differential operator. This modeling approach is representative of the state-based PD models, whereas the resulting governing equations appear as the bond-based PD models.

The refined model can be reduced to the existing bond-based PD heat conduction models by specifying particular influence functions. Also, the refined model does not require any calibration procedure unlike the bond-based PD. A systematic explicit dynamic solver is introduced to validate 1 D, 2 D and 3 D heat conduction in domains with and without a crack subjected to a combination of Dirichlet, Neumann and convection boundary conditions. All of the PD predictions are in excellent agreement with the classical solutions and demonstrate the nonlocal feature and advantage of PD in dealing with heat conduction in discontinuous domains.

The existing PD heat conduction models are reviewed. A refined bond-based PD thermal conduction model by using PD differential operator is proposed and 3 D thermal conduction in intact or cracked structures is simulated.

The purpose of this paper is to demonstrate how anomalous diffusion behaviors can be manifest in physically realizable phase change systems.

In the presence of heterogeneity the exponent in the diffusion time scale can become anomalous, exhibiting values that differ from the expected value of 1/2. Here the author investigates, through directed numerical simulation, the two-dimensional melting of a phase change material (PCM) contained in a pattern of cavities separated by a non-PCM matrix. Under normal circumstances we would expect that the progress of melting F(t) would exhibit the normal diffusion time exponent, i.e., F∼t1/2. The author’s intention is to investigate what features of the PCM cavity pattern might induce anomalous phase change, where the progress of melting has a time exponent different from n=1/2.

When the PCM cavity pattern has an internal length scale, i.e., when there is a sub-domain pattern which, when reproduced, gives us the full domain pattern, the direct simulation recovers the normal ∼t1/2 phase change behavior. When, however, there is no internal length scale, e.g., the pattern is a truncated fractal, an anomalous super diffusive behavior results with melting going as t n; n > 1/2. By studying a range of related fractal patterns, the author is able to relate the observed sub-diffusive exponent to the cavity pattern’s fractal dimension. The author also shows, how the observed behavior can be modeled with a non-local fractional diffusion treatment and how sub-diffusion phase change behavior (F∼t n; n < 1/2) results when the phase change nature of the materials in the cavity and matrix are inverted.

Although the results clearly demonstrate under what circumstances anomalous phase change behavior can be practically produced, the question of an exact theoretical relationship between the cavity pattern geometry and the observed anomalous time exponent is not known.

The clear role of the influence of heterogeneity on heat flow behavior is illustrated. Suggesting that modeling heat and fluid flow in heterogeneous systems requires careful consideration.

The novel direct simulation of melting in a two-dimensional PCM cavity pattern provides a clear illustration of anomalous behavior in a classic heat and fluid flow system and by extension provides motivation to continue the numerical investigation of anomalous and non-local behaviors and fractional calculus tools.

In this work, a modified thermoelastic model of heat conduction, including higher order of time derivative, is constructed by extending the Roychoudhuri model (TPL) (Choudhuri, 2007). In this new model, Fourier’s law of heat conduction is replaced by using Taylor series expansions, including three different phase lags for the heat flux, the thermal displacement and the temperature gradient. The generalized thermoelasticity models of Lord–Shulman (Lord and Shulman, 1967), Green and Naghdi (1991), dual-phase lag (Tzou, 1996) and three-phase lag (TPL) (Choudhuri, 2007) are obtained as special cases. The paper aims to discuss these issues.

The aim of this work is to establish a new generalized mathematical model of thermoelasticity that includes TPL in the vector of heat flux, and in the thermal displacement and temperature gradients extending TPL model (Li et al., 2019e). In this model, Fourier law of heat conduction is replaced by using Taylor series expansions to a modification of the Fourier law with introducing three different phase lags for the heat flux vector, the temperature gradient, and the thermal displacement gradient and keeping terms up with suitable higher orders.

The established high-order three-phase-lag heat conduction model reduces to the previous models of thermoelasticity as special cases.

In this paper, a TPL thermoelastic model is developed by extending the Roychoudhuri (Sherief and Raslan, 2017) model (TPL) considering the Taylor series approximation of the equation of heat conduction. This model is an alternative construction to the TPL model. The new model includes high order of TPL in the vector of heat flux, and in the thermal displacement and temperature gradients.