OpenMP + MPI parallel implementation of a numerical method for solving a kinetic equation

A two-level OpenMP + MPI parallel implementation is used to numerically solve a model kinetic equation for problems with complex three-dimensional geometry. The scalability and robustness of the method are demonstrated by computing the classical gas flow through a circular pipe of finite length and the flow past a reentry vehicle model. It is shown that the two-level model significantly speeds up the computations and improves the scalability of the method.     

Performance analysis of a hybrid MPI/OpenMP application on multi-core clusters

The mixing of shared memory and message passing programming models within a single application has often been suggested as a method for improving scientific application performance on clusters of shared memory or multi-core systems. DL_POLY, a large scale molecular dynamics application programmed using message passing programming has been modified to add a layer of shared memory threading and the performance analysed on two multi-core clusters. At lower processor numbers, the extra overheads from shared memory threading in the hybrid code outweigh performance benefits gained over the pure MPI code. On larger core counts the hybrid model performs better than pure MPI, with reduced communication time decreasing the overall runtime.     

Flow and heat transfer of a nanofluid over a nonlinearly stretching sheet: A numerical study

Steady, laminar boundary fluid flow which results from the non-linear stretching of a flat surface in a nanofluid has been investigated numerically. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis. The resulting non-linear governing equations with associated boundary conditions are solved using variational finite element method (FEM) with a local non-similar transformation. The influence of Brownian motion number (Nb), thermophoresis number (Nt), stretching parameter (n) and Lewis number (Le) on the temperature and nanoparticle concentration profiles are shown graphically. The impact of physical parameters on rate of heat transfer (−θ′(0)) and mass transfer (−?′(0)) is shown in tabulated form. Some of results have also been compared with explicit finite difference method (FDM). Excellent validation of the present numerical results has been achieved with the earlier nonlinearly stretching sheet problem of Cortell [16] for local Nusselt number without taking the effect of Brownian motion and thermophoresis.     

An efficient hybrid numerical method based on an additive scheme for solving coupled systems of singularly perturbed linear parabolic problems

We construct an efficient hybrid numerical method for solving coupled systems of singularly perturbed linear parabolic problems of reaction-diffusion type. The discretization of the coupled system is based on the use of an additive or splitting scheme on a uniform mesh in time and a hybrid scheme on a layer-adapted mesh in space. It is proven that the developed numerical method is uniformly convergent of first order in time and third order in space. The purpose of the additive scheme is to decouple the components of the vector approximate solution at each time step and thus make the computation more efficient. The numerical results confirm the theoretical convergence result and illustrate the efficiency of the proposed strategy.     

A numerical solution for the heat transfer from a wedge with non-isothermal surfaces

Summary In this paper, we present results obtained using a numerical method for calculating the development of the thermal boundary layer on a non-porous wedge following a step change in surface temperature. The method is shown to be very accurate in comparison with previous analytical solutions.

---

Résumé Cet article présente les résultats obtenus à l'aide d'une méthode numérique pour calculer le développement de la couche limite thermique sur des dièdres non poreux avec une température pariétale en échelon. La méthode se révèle très précise par comparaison avec les solutions analytiques jusqu'alors employées.

     

A fast running numerical model based on the implementation of volume forces for prediction of pressure drop in a fin tube heat exchanger

Numerical based design of geometrical structures is common when studying systems involving heat exchangers, a central component in several fields, such as industrial, vehicle and household systems. The geometrical structure of heat exchangers is generally comprised by closely placed fins and tube bundles. The creation of a mesh grid for a geometrically compact heat exchanger will result in a dense structure, which is not feasible for personal computer usage. Hence, volume forces were created based on Direct Numerical Simulations (DNS) on a Flow Representative Volume (FRV) of a tube fin heat exchanger in an internal duct system of a heat pump tumble dryer. A relation of the volume averaged velocity and the volume averaged force was established in two different FRV models with a finite element simulation in COMSOL. This relation was subsequently used to create flow resistance coefficients based on volume averaged expressions of fluid velocity and volume forces. These flow resistance coefficients were implemented in two respective porous models, which represent the entire heat exchanger except the interior arrangements of fins and tube bundles. Hence, the computation time was reduced thanks to the absence of a dense mesh grid. Experimental results of the entire heat exchanger showed good agreement with the second porous model in terms of pressure drop and volume flow rate.     

A new relaxation framework for quadratic assignment problems based on matrix splitting

Quadratic assignment problems (QAPs) are known to be among the hardest discrete optimization problems. Recent study shows that even obtaining a strong lower bound for QAPs is a computational challenge. In this paper, we first discuss how to construct new simple convex relaxations of QAPs based on various matrix splitting schemes. Then we introduce the so-called symmetric mappings that can be used to derive strong cuts for the proposed relaxation model. We show that the bounds based on the new models are comparable to some strong bounds in the literature. Promising experimental results based on the new relaxations are reported.     

Boundary control of heat transfer in a three-dimensional material: A hyperbolic model

We consider a boundary value problem modeling heat propagation in a homogeneous body in the framework of a hyperbolic heat conduction law. We construct a class of controls by temperature modes on the boundary that provide the desired body temperature at a given time.     

A discrete eigenfunctions method for computing mixed hyperbolic problems based on an implicit difference scheme

In this paper numerical solutions of mixed hyperbolic problems are computed using a discrete eigenfunctions method combined with an implicit difference scheme. This new numerical technique preserves the qualitative properties of the analytic solution due to the Sturm-Liouville structure of the underlying discrete linear boundary-value problem and has computational stability advantages vs other methods. Illustrative examples are included.     

Numerical study of Sivashinsky equation using a splitting scheme based on Crank‐Nicolson method

The mathematical modeling of a planar solid‐liquid interface in the solidification of a dilute binary alloy is formulating by one of nonintegrable, nonlinear evolution equation known as Sivashinsky equation. In the first part of this paper, the mathematical modeling of Sivashinsky equation is briefly discussed. Since, the exact solutions of this equation is yet unknown, obtaining its numerical solution plays an important role to simulate its behavior. Therefore, in the second part, a second‐order splitting finite difference scheme, based on Crank‐Nicolson method, is investigated to approximate the solution of the Sivashinsky equation with homogeneous boundary conditions. We prove the solvability of the present scheme and establish the error estimate of the numerical scheme.     

A study of splitting scheme for hyperbolic conservation laws with source terms

This study deals with the convergence of a numerical scheme for conservation laws including source terms. A splitting method for source term integration is presented. More precisely, the convergence of the numerical solution towards the entropy solution is proved in the scalar case. Because of the effect of source term, the constructed scheme is total variation bounded. Numerical experiments for one-dimensional shallow water equation are presented to demonstrate the performance of the scheme.     

A numerical scheme based on a solution of nonlinear advection–diffusion equations

A mathematical formulation of the two-dimensional Cole–Hopf transformation is investigated in detail. By making use of the Cole–Hopf transformation, a nonlinear two-dimensional unsteady advection–diffusion equation is transformed into a linear equation, and the transformed equation is solved by the spectral method previously proposed by one of the authors. Thus a solution to initial value problems of nonlinear two-dimensional unsteady advection–diffusion equations is derived. On the base of the solution, a numerical scheme explicit with respect to time is presented for nonlinear advection–diffusion equations. Numerical experiments show that the present scheme possesses the total variation diminishing properties and gives solutions with good quality.     

Optimal boundary control of heat transfer in a three-dimensional material: A hyperbolic model

We consider a boundary value problem that describes heat propagation in a homogeneous isotropic body in the framework of the hyperbolic heat conduction model. We construct a class of boundary data (controls) depending on a functional parameter and providing a given distribution of the body temperature at a given time. By using the Lagrange method, from the constructed class, we choose a subclass of controls minimizing a given loss function.     

A numerical scheme for a class of sweeping processes

The aim of this paper is to study a whole class of first order differential inclusions, which fit into the framework of perturbed sweeping process by uniformly prox-regular sets. After obtaining well-posedness results, we propose a numerical scheme based on a prediction-correction algorithm and we prove its convergence. Finally we apply these results to a problem coming from the modelling of crowd motion.     

A numerical investigation of blow-up in reaction-diffusion problems with traveling heat sources

This paper studies the numerical solution of a reaction-diffusion differential equation with traveling heat sources. According to the fact that the locations of heat sources are known, we add auxiliary mesh points exactly at heat sources and present a novel moving mesh algorithm for solving the problem. Several examples are provided to demonstrate the efficiency of the new moving mesh method, especially in the case of two or three traveling heat sources. Moreover, numerical results illustrate that the speed of the movement of the heat source is critical for blow-up when there is only one traveling heat source. For the case of two traveling heat sources, blow-up depends not only on the speed but also on the distance between the two traveling heat sources.     

Complexity analysis and numerical implementation of a short-step primal-dual algorithm for linear complementarity problems

In this paper we deal with the study of the polynomial complexity and numerical implementation for a short-step primal-dual interior point algorithm for monotone linear complementarity problems LCP. The analysis is based on a new class of search directions used by the author for convex quadratic programming (CQP) [M. Achache, A new primal-dual path-following method for convex quadratic programming, Computational and Applied Mathematics 25 (1) (2006) 97-110]. Here, we show that this algorithm enjoys the best theoretical polynomial complexity namely , iteration bound. For its numerical performances some strategies are used. Finally, we have tested this algorithm on some monotone linear complementarity problems.     

A numerical scheme based on mean value solutions for the helmholtz equation on triangular grids

A numerical treatment for the Dirichlet boundary value problem on regular triangular grids for homogeneous Helmholtz equations is presented, which also applies to the convection-diffusion problems. The main characteristic of the method is that an accuracy estimate is provided in analytical form with a better evaluation than that obtained with the usual finite difference method. Besides, this classical method can be seen as a truncated series approximation to the proposed method. The method is developed from the analytical solutions for the Dirichlet problem on a ball together with an error evaluation of an integral on the corresponding circle, yielding accuracy. Some numerical examples are discussed and the results are compared with other methods, with a consistent advantage to the solution obtained here.

     

Analytic solutions to heat transfer problems on a basis of determination of a front of heat disturbance

With the use of additional boundary conditions in integral method of heat balance, we obtain analytic solution to nonstationary problem of heat conductivity for infinite plate. Relying on determination of a front of heat disturbance, we perform a division of heat conductivity process into two stages in time. The first stage comes to the end after the front of disturbance arrives the center of the plate. At the second stage the heat exchange occurs at the whole thickness of the plate, and we introduce an additional sought-for function which characterizes the temperature change in its center. Practically the assigned exactness of solutions at both stages is provided by introduction on boundaries of a domain and on the front of heat perturbation the additional boundary conditions. Their fulfillment is equivalent to the sought-for solution in differential equation therein. We show that with the increasing of number of approximations the accuracy of fulfillment of the equation increases. Note that the usage of an integral of heat balance allows the application of the given method for solving differential equations that do not admit a separation of variables (nonlinear, with variable physical properties etc.).     

Finite volume ADI scheme for hybrid dimension heat conduction problems set in a cross-shaped domain

Lithuanian Mathematical Journal - In this paper, we construct an alternating direction implicit (ADI) type finite volume numerical scheme to solve a nonclassical nonstationary heat conduction...