Search results
1 – 10 of 233A. Huerta and F. Casadei
The arbitrary Lagrangian—Eulerian (ALE)formulation, which is already well established in the hydrodynamics andfluid‐structure interaction fields, is extended to materials…
Abstract
The arbitrary Lagrangian—Eulerian (ALE) formulation, which is already well established in the hydrodynamics and fluid‐structure interaction fields, is extended to materials with memory, namely, non‐ linear path‐dependent materials. Previous attempts to treat non‐ linear solid mechanics with the ALE description have, in common, the implicit interpolation technique employed. Obviously, this implies a numerical burden which may be uneconomical and may induce to give up this formulation, particularly in fast‐transient dynamics where explicit algorithms are usually employed. Here, several applications are presented to show that if adequate stress updating techniques are implemented, the ALE formulation could be much more competitive than classical Lagrangian computations when large deformations are present. Moreover, if the ALE technique is interpreted as a simple interpolation enrichment, adequate—in opposition to distorted or locally coarse—meshes are employed. Notice also that impossible computations (or at least very involved numerically) with a Lagrangian code are easily implementable in an ALE analysis. Finally, it is important to observe that the numerical examples shown range from a purely academic test to real engineering simulations. They show the effective applicability of this formulation to non‐linear solid mechanics and, in particular, to impact, coining or forming analysis.
Details
Keywords
Ali Ayyed Abdul-Kadhim, Fue-Sang Lien and Eugene Yee
This study aims to modify the standard probabilistic lattice Boltzmann methodology (LBM) cellular automata (CA) algorithm to enable a more realistic and accurate computation of…
Abstract
Purpose
This study aims to modify the standard probabilistic lattice Boltzmann methodology (LBM) cellular automata (CA) algorithm to enable a more realistic and accurate computation of the ensemble rather than individual particle trajectories that need to be updated from one time step to the next (allowing, as such, a fraction of the collection of particles in any lattice grid cell to be updated in a time step, rather than the entire collection of particles as in the standard LBM-CA algorithm leading to a better representation of the dynamic interaction between the particles and the background flow). Exploitation of the inherent parallelism of the modified LBM-CA algorithm to provide a computationally efficient scheme for computation of particle-laden flows on readily available commodity general-purpose graphics processing units (GPGPUs).
Design/methodology/approach
This paper presents a framework for the implementation of a LBM for the simulation of particle transport and deposition in complex flows on a GPGPU. Towards this objective, the authors have shown how to map the data structure of the LBM with a multiple-relaxation-time (MRT) collision operator and the Smagorinsky subgrid-scale turbulence model (for turbulent fluid flow simulations) coupled with a CA probabilistic method (for particle transport and deposition simulations) to a GPGPU to give a high-performance computing tool for the calculation of particle-laden flows.
Findings
A fluid-particle simulation using our LBM-MRT-CA algorithm run on a single GPGPU was 160 times as computationally efficient as the same algorithm run on a single CPU.
Research limitations/implications
The method is limited by the available computational resources (e.g. GPU memory size).
Originality/value
A new 3D LBM-MRT-CA model was developed to simulate the particle transport and deposition in complex laminar and turbulent flows with different hydrodynamic characteristics (e.g. vortex shedding, impingement, free shear layer, turbulent boundary layer). The solid particle information is encapsulated locally at the lattice grid nodes, allowing for straightforward mapping of the datastructure onto a GPGPU enabling a massive parallel execution of the LBM-MRT-CA algorithm. The new particle transport algorithm was based on the local (bulk) particle density and velocity and provides more realistic results for the particle transport and deposition than the standard LBM-CA algorithm.
Details
Keywords
Seth Dillard, James Buchholz, Sarah Vigmostad, Hyunggun Kim and H.S. Udaykumar
The performance of three frequently used level set-based segmentation methods is examined for the purpose of defining features and boundary conditions for image-based Eulerian…
Abstract
Purpose
The performance of three frequently used level set-based segmentation methods is examined for the purpose of defining features and boundary conditions for image-based Eulerian fluid and solid mechanics models. The focus of the evaluation is to identify an approach that produces the best geometric representation from a computational fluid/solid modeling point of view. In particular, extraction of geometries from a wide variety of imaging modalities and noise intensities, to supply to an immersed boundary approach, is targeted.
Design/methodology/approach
Two- and three-dimensional images, acquired from optical, X-ray CT, and ultrasound imaging modalities, are segmented with active contours, k-means, and adaptive clustering methods. Segmentation contours are converted to level sets and smoothed as necessary for use in fluid/solid simulations. Results produced by the three approaches are compared visually and with contrast ratio, signal-to-noise ratio, and contrast-to-noise ratio measures.
Findings
While the active contours method possesses built-in smoothing and regularization and produces continuous contours, the clustering methods (k-means and adaptive clustering) produce discrete (pixelated) contours that require smoothing using speckle-reducing anisotropic diffusion (SRAD). Thus, for images with high contrast and low to moderate noise, active contours are generally preferable. However, adaptive clustering is found to be far superior to the other two methods for images possessing high levels of noise and global intensity variations, due to its more sophisticated use of local pixel/voxel intensity statistics.
Originality/value
It is often difficult to know a priori which segmentation will perform best for a given image type, particularly when geometric modeling is the ultimate goal. This work offers insight to the algorithm selection process, as well as outlining a practical framework for generating useful geometric surfaces in an Eulerian setting.
Details
Keywords
Michel Bellet, Olivier Jaouen and Isabelle Poitrault
The present paper addresses the computer modelling of pipe formation in metal castings.
Abstract
Purpose
The present paper addresses the computer modelling of pipe formation in metal castings.
Design/methodology/approach
As a preliminary, a brief review of the current state‐of‐the‐art in pipe shrinkage computation is presented. Then, in first part, the constitutive equations that have to be considered in thermomechanical computations are presented, followed by the main lines of the mechanical finite element resolution. A detailed presentation of an original arbitrary Lagrangian‐Eulerian (ALE) formulation is given, explaining the connection between the Lagrangian and the quasi Eulerian zones, and the treatment of free surfaces.
Findings
Whereas most existing methods are based on thermal considerations only, it is demonstrated in the current paper that this typical evolution of the free surface, originated by shrinkage at solidification front and compensating feeding liquid flow, can be effectively approached by a thermomechanical finite element analysis.
Research limitations/implications
Future work should deal with the following points: identification of thermo‐physical and rheological data, automatic and adaptive mesh refinement, calculation of the coupled deformation of mold components, development of a two‐phase solid/liquid formulation.
Practical implications
An example of industrial application is given. The proposed method has been implemented in the commercial software THERCAST® dedicated to casting simulation.
Originality/value
The proposed numerical methods provide a comprehensive approach, capable of modelling concurrently all the main phenomena participating in pipe formation.
Details
Keywords
F. H. Hamdan and P. J. Dowling
This paper, which is concerned with fluid‐structure interactionanalysis, is a sequel to our earlier paper which gave an introduction to thenumerical treatment of such systems. The…
Abstract
This paper, which is concerned with fluid‐structure interaction analysis, is a sequel to our earlier paper which gave an introduction to the numerical treatment of such systems. The paper is divided into five main sections. In the first two, a state‐of‐the‐art review on near‐field and far‐field fluid structure interaction is presented. In attempting to highlight where current research should be directed, only the most widely used computer codes are reviewed in the third section. Conclusions are presented in the fourth section.
Details
Keywords
M. Vahdati and M. Imregun
Presents a finite element/volume method for non‐linear aeroelasticity analyses of turbomachinery blades. The method uses an Arbitrary Lagrangian‐Eulerian (ALE) kinematical…
Abstract
Presents a finite element/volume method for non‐linear aeroelasticity analyses of turbomachinery blades. The method uses an Arbitrary Lagrangian‐Eulerian (ALE) kinematical description of the fluid domain, in which the grid points can be displaced independently of the fluid motion. In addition, it employs an iterative implicit formulation similar to that of the Implicit‐continuous Eulerian (ICE) technique, making it applicable to flows at all speeds. A deforming mesh capability that can move the grid to conform continuously to the instantaneous shape of an aeroelastically deforming body without excessive distortion is also included in the algorithm. The unsteady aerodynamic loads are obtained using inviscid Euler equations. The model for the solid is general and can accommodate any spatial or modal representation of the structure. Determines the flutter stability of the system by studying the aeroelastic time response histories which are obtained by integration of the coupled equations of motion for both the fluid and the structure. Develops and demonstrates in 2D the formulation, which includes several corrections for better numerical stability. The cases studied include NACA64A006 and NACA0012 aerofoils and the EPFL Configuration 4 cascade. Finds the results from the numerical indicate good overall agreement with other published work and hence demonstrates the suitability of an ICED‐ALE formulation for turbomachinery applications.
Details
Keywords
This paper gives a review of the finite element techniques (FE)applied in the area of material processing. The latest trends in metalforming, non‐metal forming and powder…
Abstract
This paper gives a review of the finite element techniques (FE) applied in the area of material processing. The latest trends in metal forming, non‐metal forming and powder metallurgy are briefly discussed. The range of applications of finite elements on the subjects is extremely wide and cannot be presented in a single paper; therefore the aim of the paper is to give FE users only an encyclopaedic view of the different possibilities that exist today in the various fields mentioned above. An appendix included at the end of the paper presents a bibliography on finite element applications in material processing for the last five years, and more than 1100 references are listed.
Details
Keywords
S.M. Hosseini, M.T. Manzari and S.K. Hannani
This paper sets out to present a fully explicit smoothed particle hydrodynamics (SPH) method to solve non‐Newtonian fluid flow problems.
Abstract
Purpose
This paper sets out to present a fully explicit smoothed particle hydrodynamics (SPH) method to solve non‐Newtonian fluid flow problems.
Design/methodology/approach
The governing equations are momentum equations along with the continuity equation which are described in a Lagrangian framework. A new treatment similar to that used in Eulerian formulations is applied to viscous terms, which facilitates the implementation of various inelastic non‐Newtonian models. This approach utilizes the exact forms of the shear strain rate tensor and its second principal invariant to calculate the shear stress tensor. Three constitutive laws including power‐law, Bingham‐plastic and Herschel‐Bulkley models are studied in this work. The imposition of the incompressibility is fulfilled using a penalty‐like formulation which creates a trade‐off between the pressure and density variations. Solid walls are simulated by the boundary particles whose positions are fixed but contribute to the field variables in the same way as the fluid particles in flow field.
Findings
The performance of the proposed algorithm is assessed by solving three test cases including a non‐Newtonian dam‐break problem, flow in an annular viscometer using the aforementioned models and a mud fluid flow on a sloping bed under an overlying water. The results obtained by the proposed SPH algorithm are in close agreement with the available experimental and/or numerical data.
Research limitations/implications
In this work, only inelastic non‐Newtonian models are studied. This paper deals with 2D problems, although extension of the proposed scheme to 3D is straightforward.
Practical implications
This study shows that various types of flow problems involving fluid‐solid and fluid‐fluid interfaces can be solved using the proposed SPH method.
Originality/value
Using the proposed numerical treatment of viscous terms, a unified and consistent approach was devised to study various non‐Newtonian flow models.
Details
Keywords
A. Mitter, J.P. Malhotra and H.T. Jadeja
A modelling approach of gas solid flow, considering different physical phenomenon such as fluid turbulence, particle turbulence and interparticle collision effects are presented…
Abstract
A modelling approach of gas solid flow, considering different physical phenomenon such as fluid turbulence, particle turbulence and interparticle collision effects are presented. The approach is based on the two‐fluid model formulation where both phases are treated as continuum. This implies that the gas phase as well as the particle phase are weighted by their separate volumetric fractions. According to the experimental results and numerical simulations, the inter‐particle collision possesses a significant influence of turbulence level on particle transport properties in gas solid turbulent flow even for dispersed phase volume fraction (α<0.01). Comparisons in predictions have been depicted with inclusion of interparticle collision effect in the equation of particle turbulent kinetic energy and with exclusion of this effect. Experimental research has been conducted in a thermal power plant depicting higher erosion resistance of noncircular square sectioned coal pipe bends in comparison with those with circular cross section, the salient features of the experimental work are presented in this paper. Experiments have been conducted to determine, pressure drop in straight and curved portions of conduits conveying air coal mixtures in a thermal power plant. Validation of this experimental data with numerical predictions have been presented.
Details
Keywords
Marianne Francois, Eray Uzgoren, Jelliffe Jackson and Wei Shyy
Multiphase flow computations involve coupled momentum, mass and energy transfer between moving and irregularly shaped boundaries, large property jumps between materials, and…
Abstract
Multiphase flow computations involve coupled momentum, mass and energy transfer between moving and irregularly shaped boundaries, large property jumps between materials, and computational stiffness. In this study, we focus on the immersed boundary technique, which is a combined Eulerian‐Lagrangian method, to investigate the performance improvement using the multigrid technique in the context of the projection method. The main emphasis is on the interplay between the multigrid computation and the effect of the density and viscosity ratios between phases. Two problems, namely, a rising bubble in a liquid medium and impact dynamics between a liquid drop and a solid surface are adopted. As the density ratio increases, the single grid computation becomes substantially more time‐consuming; with the present problems, an increase of factor 10 in density ratio results in approximately a three‐fold increase in CPU time. Overall, the multigrid technique speeds up the computation and furthermore, the impact of the density ratio on the CPU time required is substantially reduced. On the other hand, the impact of the viscosity ratio does not play a major role on the convergence rates.
Details