Search results

We consider, as discretization in space, the nonconforming mesh developed in SUSHI (Scheme Using Stabilization and Hybrid Interfaces) developed in Eymard et al. (2010) for a semi-linear heat equation. The time discretization is performed using a uniform mesh. We are concerned with a nonlinear scheme that has been studied in Bradji (2016) in the context of the general framework GDM (Gradient Discretization Method) (Droniou et al., 2018) which includes SUSHI. We provide sufficient conditions on the size of the spatial mesh and the time step which allow to prove a W1,∞(L2)-error estimate. This error estimate can be viewed as an improvement for the W1,2(L2)-error estimate proved in Bradji (2016). The W1,∞(L2)-error estimate we want to prove in this note was stated without proof in Bradji (2016, Remark 7.2, Page 1302). Its proof is based on a comparison with an appropriately chosen auxiliary finite volume scheme along with the derivation of some new estimates on its solution.

The main purpose of this paper is to introduce the gradient discretisation method (GDM) to a system of reaction diffusion equations subject to non-homogeneous Dirichlet boundary conditions. Then, the authors show that the GDM provides a comprehensive convergence analysis of several numerical methods for the considered model. The convergence is established without non-physical regularity assumptions on the solutions.

In this paper, the authors use the GDM to discretise a system of reaction diffusion equations with non-homogeneous Dirichlet boundary conditions.

The authors provide a generic convergence analysis of a system of reaction diffusion equations. The authors introduce a specific example of numerical scheme that fits in the gradient discretisation method. The authors conduct a numerical test to measure the efficiency of the proposed method.

This work provides a unified convergence analysis of several numerical methods for a system of reaction diffusion equations. The generic convergence is proved under the classical assumptions on the solutions.

This study aims to use new formula derived based on the shifted Jacobi functions have been defined and some theorems of the left- and right-sided fractional derivative for them have been presented.

In this article, the authors apply the method of lines (MOL) together with the pseudospectral method for solving space-time partial differential equations with space left- and right-sided fractional derivative (SFPDEs). Then, using the collocation nodes to reduce the SFPDEs to the system of ordinary differential equations, which can be solved by the ode45 MATLAB toolbox.

Applying the MOL method together with the pseudospectral discretization method converts the space-dependent on fractional partial differential equations to the system of ordinary differential equations.

This paper contributes to gain choosing the shifted Jacobi functions basis with special parameters a, b and give the authors this opportunity to obtain the left- and right-sided fractional differentiation matrices for this basis exactly. The results of the examples are presented in this article. The authors found that the method is efficient and provides accurate results, and the authors found significant implications for success in the science, technology, engineering and mathematics domain.

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.

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.

Confronting the unveiled sophisticated structural and physical characteristics of permanent magnets, notably the samarium–cobalt (Sm-Co) alloy, This work aims to introduce a simulation scheme that can link physics-based micromagnetics on the nanostructures and magnetostatic homogenization on the mesoscale polycrystalline structures.

The simulation scheme is arranged in a multiscale fashion. The magnetization behaviors on the nanostructures examined with various orientations are surrogated as the micromagnetic-informed hysterons. The hysteresis behavior of the mesoscale polycrystalline structures with micromagnetic-informed hysterons is then evaluated by computational magnetostatic homogenization.

The micromagnetic-informed hysterons can emulate the magnetization reversal of the parameterized Sm-Co nanostructures as the local hysteresis behavior on the mesostructures. The simulation results of the mesoscale polycrystal demonstrate that the demagnetization process starts from the grain with the largest orientation angle (a) and then propagates to the surrounding grains.

The presented scheme depicts the demand for integrating data-driven methods, as the parameters of the surrogate hysteron intrinsically depend on the nanostructure and its orientation. Further hysteron parameters that help the surrogate hysteron emulate the micromagnetic-simulated magnetization reversal should be examined.

This work provides a novel multiscale scheme for simulating the polycrystalline permanent magnets’ hysteresis while recapitulating the nanoscale mechanisms, such as the nucleation of domains, and domain wall migration and pinning. This scheme can be further extended to simulate the part-level hysteresis considering the mesoscale features.

The purpose of this study is to analyze the nonhomogeneous model on the mixed convection of Al2O3–Fe3O4 Bingham plastic hybrid nanofluid in a ventilated enclosure subject to an externally imposed uniform magnetic field. Entropy generation and the pressure drop are determined to analyze the performance of the heat transfer. The significance of Joule heating arising due to the applied magnetic field on the heat transfer of the yield stress fluid is described.

The ventilation in the enclosure of heated walls is created by an opening on one vertical wall through which cold fluid is injected and another opening on the opposite vertical wall through which fluid can flow out.

This study finds that the inclusion of Fe₃O₄ nanoparticles with the Al₂O₃-viscoplastic nanofluid augments the heat transfer. This rate of enhancement in heat transfer is higher than the rate by which the entropy generation is increased as well as the enhancement in the pressure drop. The yield stress has an adverse effect on the heat transfer; however, it favors thermal mixing. The magnetic field, which is acting opposite to the direction of the inlet jet, manifests heat transfer of the viscoplastic hybrid nanofluid. The horizontal jet of cold fluid produces the optimal heat transfer.

The objective of this study is to analyze the impact of the inclined cold jet of viscoplastic electrically conducting hybrid nanofluid on heat transfer from the enclosure in the presence of a uniform magnetic field. The combined effect of hybrid nanoparticles and a magnetic field to enhance heat transfer of a viscoplastic fluid in a ventilated enclosure has not been addressed before.

The paper targets on providing new experimental data for validation of the well-established mathematical models within the framework of the lattice Boltzmann method (LBM), which are applied to problems of casting processes in complex mould cavities.

An experimental campaign aiming at the free-surface flow within a system of narrow channels is designed and executed under well-controlled laboratory conditions. An in-house lattice Boltzmann solver is implemented. Its algorithm is described in detail and its performance is tested thoroughly using both the newly recorded experimental data and well-known analytical benchmark tests.

The benchmark tests prove the ability of the implemented algorithm to provide a reliable solution when the surface tension effects become dominant. The convergence of the implemented method is assessed. The two new experimentally studied problems are resolved well by simulations using a coarse computational grid.

A detailed set of original experimental data for validation of computational schemes for simulations of free-surface gravity-driven flow within a system of narrow channels is presented.

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.