A parallel and efficient multi-split XFEM for 3-D analysis of heterogeneous materials

in Computer Methods in Applied Mechanics and Engineering (2018)

We propose a parallel and computationally efficient multi-split XFEM approach for 3-D analysis of heterogeneous materials. In this approach, multiple discontinuities (pores and reinforcement particles ... [more ▼]

We propose a parallel and computationally efficient multi-split XFEM approach for 3-D analysis of heterogeneous materials. In this approach, multiple discontinuities (pores and reinforcement particles) may intersect any given element (we call those elements multi-split elements). These discontinuities are modeled by imposing additional degrees of freedom at the nodes. The main advantage of the proposed scheme is that the mesh size remains independent of the relative distance among the heterogeneities/discontinuities. The pores and reinforcement particles are assumed to be spherical. The simulations are performed for uniform and non-uniform heterogeneity distribution. The Young’s modulus of the heterogeneous material is evaluated for different amount of pores and reinforcement particles. To demonstrate the computational efficiency of the multi-split XFEM, elastic damage analysis is performed for the unit cell with 5% pores and 5% reinforcement particles under uniaxial tensile loading. These simulations show that the Young’s modulus decreases linearly with the increase in the volume fraction of the pores and increases linearly with the increase in volume fraction of reinforcement particles. The multi-split XFEM is found to be at least 1.8 times computationally efficient than standard XFEM and at least 6.7 times computationally efficient than FEM. [less ▲]

A Parallel Multiscale DEM-VOF Method For Large-Scale Simulations Of Three-Phase Flows

in Proceedings of ECCM-ECFD 2018 (2018)

A parallel dual-grid multiscale DEM-VOF coupling is here investigated. Dual- grid multiscale couplings have been recently used to address different engineering problems involving the interaction between ... [more ▼]

A parallel dual-grid multiscale DEM-VOF coupling is here investigated. Dual- grid multiscale couplings have been recently used to address different engineering problems involving the interaction between granular phases and complex fluid flows. Nevertheless, previous studies did not focus on the parallel performance of such a coupling and were, therefore, limited to relatively small applications. In this contribution, we propose an insight into the performance of the dual-grid multiscale DEM-VOF method for three- phase flows when operated in parallel. In particular,we focus on a famous benchmark case for three-phase flows and assess the influence of the partitioning algorithm on the scalability of the dual-grid algorithm. [less ▲]

Parameter identification problem in bimaterial human skin and sensitivity analysis : Uncertainties in biomechanics of skin

Scientific Conference (2019, December 09)

The proposed paper concerns the prediction of the numerical response of a biomechanical structure submitted to an unknown external loading state. The methodology is based on homogeneous and then ... [more ▼]

The proposed paper concerns the prediction of the numerical response of a biomechanical structure submitted to an unknown external loading state. The methodology is based on homogeneous and then heterogeneous structures such as healthy or pathological cutaneous tissues that can be mechanically tested in vivo under a patchy knowledge of boundary conditions. Experimental data corresponding to the extension of a piece of skin located between two pads with displacement enslavement, represent input data to the numerical model. Data are reaction force on one pad and displacement field between the two pads and all around. The numerical model consists of a representation of the bi-material domain geometry with neo-hookean behaviors. The boundary conditions and loadings of the experimental extension test are imposed. The materials parameters have been identified by inverse method starting from a constrained cost function minimizing the difference between the calculated displacements field and experimental displacements field obtained by digital image correlation and taking into account the reaction force as a constraint. An analysis of the model sensitivity to material parameters is presented. [less ▲]

A PAS-state constrained optimization formulation for the anticipatory network traffic control problem

Scientific Conference (2016, September)

A perspective on urban economics and geographical agent-based models. (Keynote Speech)

Scientific Conference (2015, September)

In this presentation, we will argue that there are many benefits to further liaise urban economic models and geographically -aware agent-based simulation models for further understanding the structure and ... [more ▼]

In this presentation, we will argue that there are many benefits to further liaise urban economic models and geographically -aware agent-based simulation models for further understanding the structure and dynamics of cities, especially at the scale of city regions but also, eventually, for contributing micro understanding of changes in a system of cities. At first sight it appears that these two fields of urban research are very separated: urban economics relies mostly on mathematical proofs, a parsimonious set of interactions, static equilibria and a rather homogeneous space and set of agents; while geographical agent-based models emphasize the dynamics and complexity of cities, the heterogeneity of agents and space, and rely solely on computational methods. The presentation will showcase a series of theoretical and applied research to stress the idea that the two can fruitfully be combined and that results are then palatable to both urban economists and geographers and can therefore participate to consolidating a wider and formalised urban theory. The geographical turn towards agent-based approaches, compared to more aggregate models and behaviourally-blind geocomputational methods, is actually a chance for geographers to formalise market as well as non-market interactions in a manner that is consistent with urban economic theory and enhance the treatment of space in this theory. Relying on a formal description of agents utility permits to discuss normative instruments where environmental and welfare impacts can be traded-off by policy makers. Furthermore,, land or housing market outcomes of the models and the analytical efforts actually help to reduce parametric space and calibrate or validate models against independent data. [less ▲]

PIV study of the effect of piston position on the in-cylinder swirling flow during the scavenging process in large two-stroke marine diesel engines

in Journal of Marine Science and Technology (2013), 18(1), 133-143

A simplified model of a low speed large two- stroke marine diesel engine cylinder is developed. The effect of piston position on the in-cylinder swirling flow during the scavenging process is studied ... [more ▼]

A simplified model of a low speed large two- stroke marine diesel engine cylinder is developed. The effect of piston position on the in-cylinder swirling flow during the scavenging process is studied using the stereo- scopic particle image velocimetry technique. The mea- surements are conducted at different cross-sectional planes along the cylinder length and at piston positions covering the air intake port by 0, 25, 50 and 75 %. When the intake port is fully open, the tangential velocity profile is similar to a Burgers vortex, whereas the axial velocity has a wake- like profile. Due to internal wall friction, the swirl decays downstream, and the size of the vortex core increases. For increasing port closures, the tangential velocity profile changes from a Burgers vortex to a forced vortex, and the axial velocity changes correspondingly from a wake-like profile to a jet-like profile. For piston position with 75 % intake port closure, the jet-like axial velocity profile at a cross-sectional plane close to the intake port changes back to a wake-like profile at the adjacent downstream cross-sectional plane. This is characteristic of a vortex breakdown. The non-dimensional velocity profiles show no significant variation with the variation in Reynolds number. [less ▲]

POD-based model order reduction for the simulation of strong nonlinear evolutions in structures: Application to damage propagation

in IOP Conference Series: Materials Science and Engineering (2014), 10(1),

In this paper, we develop a bridge between POD-based model order reduction techniques and the classical Newton-Krylov solvers to derive an efficient solution procedure for highly nonlinear problems ... [more ▼]

In this paper, we develop a bridge between POD-based model order reduction techniques and the classical Newton-Krylov solvers to derive an efficient solution procedure for highly nonlinear problems undergoing strong topological changes. [less ▲]

POD-based reduction methods, the Quasicontinuum Method and their Resemblance

Scientific Conference (2016, June 27)

Poroelastic model parameter identification using artificial neural networks: on the effects of heterogeneous porosity and solid matrix Poisson ratio

in Computational Mechanics (2020), 66

Predictive analysis of poroelastic materials typically require expensive and time-consuming multiscale and multiphysics approaches, which demand either several simplifications or costly experimental tests ... [more ▼]

Predictive analysis of poroelastic materials typically require expensive and time-consuming multiscale and multiphysics approaches, which demand either several simplifications or costly experimental tests for model parameter identification. This problem motivates us to develop a more efficient approach to address complex problems with an acceptable computational cost. In particular, we employ artificial neural network (ANN) for reliable and fast computation of poroelastic model parameters. Based on the strong-form governing equations for the poroelastic problem derived from asymptotic homogenisation, the weighted residuals formulation of the cell problem is obtained. Approximate solution of the resulting linear variational boundary value problem is achieved by means of the finite element method. The advantages and downsides of macroscale properties identification via asymptotic homogenisation and the application of ANN to overcome parameter characterisation challenges caused by the costly solution of cell problems are presented. Numerical examples, in this study, include spatially dependent porosity and solid matrix Poisson ratio for a generic model problem, application in tumour modelling, and utilisation in soil mechanics context which demonstrate the feasibility of the presented framework. [less ▲]

A posteriori error estimation for extended finite elements by an extended global recovery

in International Journal for Numerical Methods in Engineering (2008), 76(8), 1123-1138

This contribution presents an extended global derivative recovery for enriched finite element methods (FEMs), such as the extended FEM along with an associated error indicator. Owing to its simplicity ... [more ▼]

This contribution presents an extended global derivative recovery for enriched finite element methods (FEMs), such as the extended FEM along with an associated error indicator. Owing to its simplicity, the proposed scheme is ideally suited to industrial applications. The procedure is based on global minimization of the L2 norm of the difference between the raw strain field (C-1) and the recovered (C0) strain field. The methodology engineered in this paper extends the ideas of Oden and Brauchli (Int. J. Numer. Meth. Engng 1971; 3) and Hinton and Campbell (Int. J. Numer. Meth. Engng 1974; 8) by enriching the approximation used for the construction of the recovered derivatives (strains) with the gradients of the functions employed to enrich the approximation employed for the primal unknown (displacements). We show linear elastic fracture mechanics examples, both in simple two-dimensional settings, and for a three-dimensional structure. Numerically, we show that the effectivity index of the proposed indicator converges to unity upon mesh refinement. Consequently, the approximate error converges to the exact error, indicating that the error indicator is valid. Additionally, the numerical examples suggest a novel adaptive strategy for enriched approximations in which the dimensions of the enrichment zone are first increased, before standard h- and p-adaptivities are applied; we suggest to coin this methodology e-adaptivity. Copyright © 2008 John Wiley & Sons, Ltd. [less ▲]