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 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.

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.

It is of great significance to study the influence of subgrade filling on permafrost temperature field in permafrost area for the smooth construction and safe operation of railway.

The paper builds up the model for the hydrothermal coupling calculation of permafrost using finite element software COMSOL to study how permafrost temperature field changes in the short term after subgrade filling, on which basis it proposes the method of calculation for the concave distortion of freezing front in the subgrade-covered area.

The results show that the freezing front below the subgrade center sinks due to the thermal effect of subgrade filling, which will trigger hydrothermal erosion in case of sufficient moisture inflows, leading to the thawing settlement or the cracking of the subgrade, etc. The heat output of soil will be hindered the most in case of July filling, in which case the sinking and the distortion of the freezing front is found to be the most severe, which the recovery of the permafrost temperature field, the slowest, constituting the most unfavorable working condition. The concave distortion of the freezing front in the subgrade area increases with the increase in temperature difference between the filler and ground surface, the subgrade height, the subgrade width and the volumetric thermal capacity of filler, while decreases with the increase of the thermal conductivity of filler. Therefore, the filler chose for engineering project shall be of small volumetric thermal capacity, low initial temperature and high thermal conductivity whenever possible.

The concave distortion of the freezing front under different working conditions at different times after filling can be calculated using the method proposed.

This paper aims to analyze the pressure distribution of rectangular aerostatic thrust bearing with a single air supply inlet using the complex potential theory and conformal mapping.

The Möbius transform is used to map the interior of a rectangle onto the interior of a unit circle, from which the pressure distribution and load carrying capacity are obtained. The calculation results are verified by finite difference method.

The constructed Möbius formula is very effective for the performance characteristics researches for the rectangular thrust bearing with a single air supply inlet. In addition, it is also noted that to obtain the optimized load carrying capacity, the square thrust bearing can be adopted.

The Möbius transform is found suitable to describe the pressure distribution of the rectangular thrust bearing with a single air supply inlet.

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 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.

Incineration has become increasingly important in many large cities around the world because of fast urbanization and population growth. The benefits of energy production and large reduction in the waste volume to landfills also contribute to its growing adaptation for solid waste management for these cities. At the same time, the environmental impact of the pollutant gases emitted from the incineration process is a common concern for various stakeholders which must be properly addressed. To minimize the pollutant gas emission levels, as well as maximize the energy efficiency, it is critically important to optimize the combustion performance of an incinerator freeboard which would require the development of reliable approaches based on computational fluid dynamics (CFD) modeling. A critical task in the CFD modeling of an incinerator furnace requires the specification of waste characteristics along the moving grate as boundary conditions, which is not available in standard CFD packages at present. This study aims to address this gap by developing a numerical incinerator waste bed model.

A one-dimensional Lagrangian model for the incineration waste bed has been developed, which can be coupled to the furnace CFD model. The changes in bed mass due to drying, pyrolysis, devolatilization and char oxidation are all included in the model. The mass and concentration of gases produced in these processes through reactions are also predicted. The one-dimensional unsteady energy equations of solid and gas phases, which account for the furnace radiation, conduction, convection and heat of reactions, are solved by the control volume method.

The Lagrangian model is validated by comparing its prediction with the experimental data in the literature. The predicted waste bed height reduction, temperature profile and gas concentration are in reasonable agreement with the observations.

The simplicity and efficiency of the model makes it ideally suitable to be used for coupling with the computational furnace model to be developed in future (so as to optimize incinerator designs).

One of the major challenges in the design of thermal equipment is to minimize the entropy production and enhance the thermal dissipation rate for improving energy efficiency of the devices. In several industrial applications, the structure of thermal device is cylindrical shape. In this regard, this paper aims to explore the impact of isothermal cylindrical solid block on nanofluid (Ag – H₂O) convective flow and entropy generation in a cylindrical annular chamber subjected to different thermal conditions. Furthermore, the present study also addresses the structural impact of cylindrical solid block placed at the center of annular domain.

The alternating direction implicit and successive over relaxation techniques are used in the current investigation to solve the coupled partial differential equations. Furthermore, estimation of average Nusselt number and total entropy generation involves integration and is achieved by Simpson and Trapezoidal’s rules, respectively. Mesh independence checks have been carried out to ensure the accuracy of numerical results.

Computations have been performed to analyze the simultaneous multiple influences, such as different thermal conditions, size and aspect ratio of the hot obstacle, Rayleigh number and nanoparticle shape on buoyancy-driven nanoliquid movement, heat dissipation, irreversibility distribution, cup-mixing temperature and performance evaluation criteria in an annular chamber. The computational results reveal that the nanoparticle shape and obstacle size produce conducive situation for increasing system’s thermal efficiency. Furthermore, utilization of nonspherical shaped nanoparticles enhances the heat transfer rate with minimum entropy generation in the enclosure. Also, greater performance evaluation criteria has been noticed for larger obstacle for both uniform and nonuniform heating.

The current numerical investigation can be extended to further explore the thermal performance with different positions of solid obstacle, inclination angles, by applying Lorentz force, internal heat generation and so on numerically or experimentally.

A pioneering numerical investigation on the structural influence of hot solid block on the convective nanofluid flow, energy transport and entropy production in an annular space has been analyzed. The results in the present study are novel, related to various modern industrial applications. These results could be used as a firsthand information for the design engineers to obtain highly efficient thermal systems.