Search results

This paper aims to use the Lattice Boltzmann method (LBM) to numerically simulate the natural convection heat transfer of Cu-water nanofluid in an L-shaped enclosure with curved boundaries.

LBM on three different models of curved L-shape cavity using staircase approach is applied to perform a comparative investigation for the effects of curved boundary on fluid flow and heat transfer. The staircase approximation is a straightforward and efficient approach to simulating curved boundaries in LBM.

The effect of curved boundary on natural convection in different parameter ranges of Rayleigh number and nanoparticle volume fraction is investigated. The curved L-shape results are also compared to the rectangular L-shape results that were also achieved in this study. The curved boundary LBM simulation is also validated with existing studies, which shows great accuracy in this study. The results show that the top curved boundary in curved L-shape models causes a notable increase in the Nusselt number values.

Based on existing literature, there is a lack of comparative studies which would specifically examine the effects of curved boundaries on natural convection in closed cavities. Particularly, the application of curved boundaries to an L-shape cavity has not been examined. In this study, curved boundaries are applied to the sharp corners of the bending section in the L-shape and the results of the curved L-shape models are compared to the simple rectangular L-shape model. Hence, a comparative evaluation is performed for the effect of curved boundaries on fluid flow in the L-shape enclosure.

An extension of the Schwarz‐Christoffel transformation is described to formally map polygons which contain curved boundaries. The curved boundaries are divided into small ‘curved elements’ and each element is approximated by a second degree polynomial (higher degree polynomials can also be used). The iterative algorithm of evaluating the unknown constants of the basic S‐C transformation described in a companion paper is applied to the extended S‐C transformation to compute its unknown constants, including the coefficients of the polynomials. Excellent results are achieved as far as accuracy and convergence are concerned. Examples including a practical application, are provided. The mapping of curved polygons is important because they provide a better model of a physical device.

Present study aims to extend the lattice Boltzmann method (LBM) to simulate radiation in geometries with curved boundaries, as the first step to simulate radiation in complex porous media. In recent years, researchers have increasingly explored the use of porous media to improve the heat transfer processes. The lattice Boltzmann method (LBM) is one of the most effective techniques for simulating heat transfer in such media. However, the application of the LBM to study radiation in complex geometries that contain curved boundaries, as found in many porous media, has been limited.

The numerical evaluation of the effect of the radiation-conduction parameter and extinction coefficient on temperature and incident radiation distributions demonstrates that the proposed LBM algorithm provides highly accurate results across all cases, compared to those found in the literature or those obtained using the finite volume method (FVM) with the discrete ordinates method (DOM) for radiative information.

For the case with a conduction-radiation parameter equal to 0.01, the maximum relative error is 1.9% in predicting temperature along vertical central line. The accuracy improves with an increase in the conduction-radiation parameter. Furthermore, the comparison between computational performances of two approaches reveals that the LBM-LBM approach performs significantly faster than the FVM-DOM solver.

The difficulty of radiative modeling in combined problems involving irregular boundaries has led to alternative approaches that generally increase the computational expense to obtain necessary radiative details. To address the limitations of existing methods, this study presents a new approach involving a coupled lattice Boltzmann and first-order blocked-off technique to efficiently model conductive-radiative heat transfer in complex geometries with participating media. This algorithm has been developed using the parallel lattice Boltzmann solver.

A scheme is presented for the automatic generation of triangular meshes over general analytical curved surfaces with explicit control on the discretization error over the interior of the domain. The element size on the curved surface is estimated based on the given allowable discretization error and the surface curvature. The mesh generation starts with the subdivision of the curved boundary of the domain into straight line segments in compliance with the given discretization error constraint. Elements are generated directly over the surface in such a way that the discretization error is kept within the required geometrical tolerance. Those elements violating the given constraint are noted and post‐processed by either element subdivision or node repositioning. Various schemes as to how improvement should be made are proposed, which reduce the discretization error by different strategies depending on the position of the point where the maximum discretization error occurs.

This paper aims to develop a new isogeometric boundary element formulation based on non-uniform rational basis splines (NURBS) curves for solving Reissner’s shear-deformable plates.

The generalized displacements and tractions along the problem boundary are approximated as NURBS curves having the same rational B-spline basis functions used to describe the geometrical boundary of the problem. The source points positions are determined over the problem boundary by the well-known Greville abscissae definition. The singular integrals are accurately evaluated using the singularity subtraction technique.

Numerical examples are solved to demonstrate the validity and the accuracy of the developed formulation.

This formulation is considered to preserve the exact geometry of the problem and to reduce or cancel mesh generation time by using NURBS curves employed in computer aided designs as a tool for isogeometric analysis. The present formulation extends such curves to be implemented as a stress analysis tool.

To study and to analyze a second order finite‐element boundary‐flux approximation using isoparametric numerical integration.

The numerical finite‐element integration is the main method used in this research. Since a domain with curved boundary is considered we apply an isoparametric approach. The lumped flux formulation is another method of approach in this paper.

This research study presents a careful analysis of the combined effect of the numerical integration and isoparametric FEM on the boundary‐flux error. Some L₂‐norm estimates are proved for the approximate solutions of the problem under consideration.

The authors offer a general study within the framework of the boundary‐flux approximation theory, which completes the results of published works in this scientific field of research.

A useful application is to employ appropriate quadrature formulae without violating the precision of the boundary‐flux FEM. The lumped mass approximation is also an important practical approach to the problem in question.

The paper presents an entire investigation in FE boundary‐flux approximation theory, in particular, elements of arbitrary degree and domains with curved boundaries. The work is addressed to the possible related fields of interest of postgraduate students and specialists in fluid mechanics and numerical analysis.

The purpose of this paper is to extend the penalty concept to treat partial slip, free surface, contact and related boundary conditions in viscous flow simulation.

The penalty partial‐slip formulation is analysed and related to the classical Navier slip condition. The same penalty scheme also allows partial penetration through a boundary, hence the implementation of porous wall boundaries. The finite element method is used for investigating and interpreting penalty approaches to boundary conditions.

The generalised penalty approach is verified by means of a novel variant of the circular‐Couette flow problem, having partial slip on one of the cylindrical boundaries, for which an analytic solution is derived. Further verificationis provided by consideration of viscous flow over a sphere with partial slip on the surface, and comparison of numerical and classical solutions. Numerical studies illustrate the versatility of the approach.

The penalty approach is applied to some different boundaries: partial slip and partial penetration with no/full slip/penetration as limiting cases; free surface; space‐ and time‐varying boundary conditions which allow progressive contact over time. Application is made to curved and inclined boundaries. Sensitivity of flow to penalty parameters is an avenue for continued research, as is application of the penalty approach for non‐Newtonian flows.

This is the first work to show the relation between penalty formulation of boundary conditions and physical boundary conditions. It provides a method that overcomes past difficulties in implementing partial slip on boundaries of general shape, and which handles progressive contact. It also provides useful benchmark problems for future studies.

This work is concerned with the development and the numerical investigation of a hybridizable discontinuous Galerkin (HDG) method for the simulation of two‐dimensional time‐harmonic electromagnetic wave propagation problems.

The proposed HDG method for the discretization of the two‐dimensional transverse magnetic Maxwell equations relies on an arbitrary high order nodal interpolation of the electromagnetic field components and is formulated on triangular meshes. In the HDG method, an additional hybrid variable is introduced on the faces of the elements, with which the element‐wise (local) solutions can be defined. A so‐called conservativity condition is imposed on the numerical flux, which can be defined in terms of the hybrid variable, at the interface between neighbouring elements. The linear system of equations for the unknowns associated with the hybrid variable is solved here using a multifrontal sparse LU method. The formulation is given, and the relationship between the considered HDG method and a standard upwind flux‐based DG method is also examined.

The approximate solutions for both electric and magnetic fields converge with the optimal order of p+1 in L₂ norm, when the interpolation order on every element and every interface is p and the sought solution is sufficiently regular. The presented numerical results show the effectiveness of the proposed HDG method, especially when compared with a classical upwind flux‐based DG method.

The work described here is a demonstration of the viability of a HDG formulation for solving the time‐harmonic Maxwell equations through a detailed numerical assessment of accuracy properties and computational performances.

In the newly released ASTM standard specification for additive manufacturing file (AMF) format – version 1.1 – Hermite curve-based interpolation is used to refine input triangles to generate denser mesh with smoother geometry. This paper aims to study the problems of constructing smooth geometry based on Hermite interpolation on curves and proposes a solution to overcome these problems.

A formulation using triangular Bézier patch is proposed to generate smooth geometry from input polygonal models. Different configurations on the boundary curves in the formulation are analyzed to further enrich this formulation.

The study shows that the formulation given in the AMF format (version 1.1) can lead to the problems of inconsistent normals and undefined end-tangents.

The scheme has requirements on the input normals of a model, only C0 interpolation can be generated on those cases with less-proper input.

To overcome the problems of smooth geometry generation in the AMF format, the authors propose an enriched scheme for computing smooth geometry by using triangular Bézier patch. For the configurations with less-proper input, the authors adopt the Boolean sum and the Nielson’s point-opposite edge interpolation for triangular Coons patch to generate the smooth geometry as a C0 interpolant.

The structural design problem can be viewed as an iterative design loop with each iteration involving two stages for topology and shape designs with genetic algorithm (GA) as the optimization tool for both.

The topology optimization problem, which is ill posed, is regularized using a constraint on perimeter and solved using GA. The problem is formulated as one of compliance minimization subject to volume constraint for the single loading case. A dual formulation of this has been used for the multiple loading cases resulting in as many behavioral constraints as there are loading cases. The tentative topology given by the topology optimization module is taken and the domain boundary is approximated using straight lines, B‐splines and cubic spline curves and design variables are selected among the boundary defining points. Optimum boundary shape of the problem has been obtained using GA in two different ways: without stress constraints; and with stress constraints.

The proposed two stage strategy has been tested on benchmark structural optimization problems and its performance is found to be extremely good.

The strategy appears to be eminently suitable for implementation in a general purpose FE software as an add‐on module for structural design optimization.

It has been observed that the integrated topology and shape design method is robust and easy to implement in comparison with other techniques. The computing time requirements for the GA does not appear daunting in the present scenario of high performance parallel computing and improved GA techniques.