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.

This study aims to develop an experimentally validated three-dimensional numerical model for predicting different flow patterns produced with a gas dynamic virtual nozzle (GDVN).

The physical model is posed in the mixture formulation and copes with the unsteady, incompressible, isothermal, Newtonian, low turbulent two-phase flow. The computational fluid dynamics numerical solution is based on the half-space finite volume discretisation. The geo-reconstruct volume-of-fluid scheme tracks the interphase boundary between the gas and the liquid. To ensure numerical stability in the transition regime and adequately account for turbulent behaviour, the k-ω shear stress transport turbulence model is used. The model is validated by comparison with the experimental measurements on a vertical, downward-positioned GDVN configuration. Three different combinations of air and water volumetric flow rates have been solved numerically in the range of Reynolds numbers for airflow 1,009–2,596 and water 61–133, respectively, at Weber numbers 1.2–6.2.

The half-space symmetry allows the numerical reconstruction of the dripping, jetting and indication of the whipping mode. The kinetic energy transfer from the gas to the liquid is analysed, and locations with locally increased gas kinetic energy are observed. The calculated jet shapes reasonably well match the experimentally obtained high-speed camera videos.

The model is used for the virtual studies of new GDVN nozzle designs and optimisation of their operation.

To the best of the authors’ knowledge, the developed model numerically reconstructs all three GDVN flow regimes for the first time.

The purpose of this study is to develop stable, convergent and accurate numerical method for solving singularly perturbed differential equations having both small and large delay.

This study introduces a fitted nonpolynomial spline method for singularly perturbed differential equations having both small and large delay. The numerical scheme is developed on uniform mesh using fitted operator in the given differential equation.

The stability of the developed numerical method is established and its uniform convergence is proved. To validate the applicability of the method, one model problem is considered for numerical experimentation for different values of the perturbation parameter and mesh points.

In this paper, the authors consider a new governing problem having both small delay on convection term and large delay. As far as the researchers' knowledge is considered numerical solution of singularly perturbed boundary value problem containing both small delay and large delay is first being considered.

This paper aims to study qualitative properties and approximate solutions of a thermostat dynamics system with three-point boundary value conditions involving a nonsingular kernel operator which is called Atangana-Baleanu-Caputo (ABC) derivative for the first time. The results of the existence and uniqueness of the solution for such a system are investigated with minimum hypotheses by employing Banach and Schauder's fixed point theorems. Furthermore, Ulam-Hyers (UH) stability, Ulam-Hyers-Rassias UHR stability and their generalizations are discussed by using some topics concerning the nonlinear functional analysis. An efficiency of Adomian decomposition method (ADM) is established in order to estimate approximate solutions of our problem and convergence theorem is proved. Finally, four examples are exhibited to illustrate the validity of the theoretical and numerical results.

This paper considered theoretical and numerical methodologies.

This paper contains the following findings: (1) Thermostat fractional dynamics system is studied under ABC operator. (2) Qualitative properties such as existence, uniqueness and Ulam–Hyers–Rassias stability are established by fixed point theorems and nonlinear analysis topics. (3) Approximate solution of the problem is investigated by Adomain decomposition method. (4) Convergence analysis of ADM is proved. (5) Examples are provided to illustrate theoretical and numerical results. (6) Numerical results are compared with exact solution in tables and figures.

The novelty and contributions of this paper is to use a nonsingular kernel operator for the first time in order to study the qualitative properties and approximate solution of a thermostat dynamics system.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

This study aims to simulate the dendritic growth in Stokes flow by iteratively coupling a domain and boundary type meshless method.

A preconditioned phase-field model for dendritic solidification of a pure supercooled melt is solved by the strong-form space-time adaptive approach based on dynamic quadtree domain decomposition. The domain-type space discretisation relies on monomial augmented polyharmonic splines interpolation. The forward Euler scheme is used for time evolution. The boundary-type meshless method solves the Stokes flow around the dendrite based on the collocation of the moving and fixed flow boundaries with the regularised Stokes flow fundamental solution. Both approaches are iteratively coupled at the moving solid–liquid interface. The solution procedure ensures computationally efficient and accurate calculations. The novel approach is numerically implemented for a 2D case.

The solution procedure reflects the advantages of both meshless methods. Domain one is not sensitive to the dendrite orientation and boundary one reduces the dimensionality of the flow field solution. The procedure results agree well with the reference results obtained by the classical numerical methods. Directions for selecting the appropriate free parameters which yield the highest accuracy and computational efficiency are presented.

A combination of boundary- and domain-type meshless methods is used to simulate dendritic solidification with the influence of fluid flow efficiently.

Inverse problems in electromagnetism, namely, the recovery of sources (currents or charges) or system data from measured effects, are usually ill-posed or, in the numerical formulation, ill-conditioned and require suitable regularization to provide meaningful results. To test new regularization methods, there is the need of benchmark problems, which numerical properties and solutions should be well known. Hence, this study aims to define a benchmark problem, suitable to test new regularization approaches and solves with different methods.

To assess reliability and performance of different solving strategies for inverse source problems, a benchmark problem of current synthesis is defined and solved by means of several regularization methods in a comparative way; subsequently, an approach in terms of an artificial neural network (ANN) is considered as a viable alternative to classical regularization schemes. The solution of the underlying forward problem is based on a finite element analysis.

The paper provides a very detailed analysis of the proposed inverse problem in terms of numerical properties of the lead field matrix. The solutions found by different regularization approaches and an ANN method are provided, showing the performance of the applied methods and the numerical issues of the benchmark problem.

The value of the paper is to provide the numerical characteristics and issues of the proposed benchmark problem in a comprehensive way, by means of a wide variety of regularization methods and an ANN approach.

The purpose of this paper is to examine a solution strategy for coupled nonlinear magnetic-thermal problems and apply it to the heating process of a thin moving steel sheet. Performing efficient numerical simulations of induction heating processes becomes ever more important because of faster production development cycles, where the quasi steady-state solution of the problem plays a pivotal role.

To avoid time-consuming transient simulations, the eddy current problem is transformed into frequency domain and a harmonic balancing scheme is used to take into account the nonlinear BH-curve. The thermal problem is solved in steady-state domain, which is carried out by including a convective term to model the stationary heat transport due to the sheet velocity.

The presented solution strategy is compared to a classical nonlinear transient reference solution of the eddy current problem and shows good convergence, even for a small number of considered harmonics.

Numerical simulations of induction heating processes are necessary to fully understand certain phenomena, e.g. local overheating of areas in thin structures. With the presented approach it is possible to perform large 3D simulations without excessive computational resources by exploiting certain properties of the multiharmonic solution of the eddy current problem. Together with the use of nonconforming interfaces, the overall computational complexity of the problem can be decreased significantly.

Joule heating effect is a pervasive phenomenon in electro-osmotic flow because of the applied electric field and fluid electrical resistivity across the microchannels. Its effect in electro-osmotic flow field is an important mechanism to control the flow inside the microchannels and it includes numerous applications.

This research article details the numerical investigation on alterations in the profile of stream wise velocity of simple Couette-electroosmotic flow and pressure driven electro-osmotic Couette flow by the dynamic viscosity variations happened due to the Joule heating effect throughout the dielectric fluid usually observed in various microfluidic devices.

The advantages of the Joule heating effect are not only to control the velocity in microchannels but also to act as an active method to enhance the mixing efficiency. The results of numerical investigations reveal that the thermal field due to Joule heating effect causes considerable variation of dynamic viscosity across the microchannel to initiate a shear flow when EDL (Electrical Double Layer) thickness is increased and is being varied across the channel.

This research work suggest how joule heating can be used as en effective mechanism for flow control in microfluidic devices.

This study aims to analyse the non-linear losses of a porous media (stack) composed by parallel plates and inserted in a resonator tube in oscillatory flows by proposing numerical correlations between pressure gradient and velocity.

The numerical correlations origin from computational fluid dynamics simulations, conducted at the microscopic scale, in which three fluid channels representing the porous media are taken into account. More specifically, for a specific frequency and stack porosity, the oscillating pressure input is varied, and the velocity and the pressure-drop are post-processed in the frequency domain (Fast Fourier Transform analysis).

It emerges that the viscous component of pressure drop follows a quadratic trend with respect to velocity inside the stack, while the inertial component is linear also at high-velocity regimes. Furthermore, the non-linear coefficient b of the correlation ax + bx² (related to the Forchheimer coefficient) is discovered to be dependent on frequency. The largest value of the b is found at low frequencies as the fluid particle displacement is comparable to the stack length. Furthermore, the lower the porosity the higher the Forchheimer term because the velocity gradients at the stack geometrical discontinuities are more pronounced.

The main novelty of this work is that, for the first time, non-linear losses of a parallel plate stack are investigated from a macroscopic point of view and summarised into a non-linear correlation, similar to the steady-state and well-known Darcy–Forchheimer law. The main difference is that it considers the frequency dependence of both Darcy and Forchheimer terms. The results can be used to enhance the analysis and design of thermoacoustic devices, which use the kind of stacks studied in the present work.