ADRIEN SCHEUER
MIT Mathematics Postdoctoral Associate
Computational mechanics
Multiscale modelling
Model order reduction

Massachusetts Institute of Technology
MIT Mathematics
Earth Resources Laboratory  Prof. Laurent DEMANET
 ascheuer [at] mit.edu
8 Research

Postdoctoral Associate in Mathematics2019  Present
Massachusetts Institute of Technology, USA
BAEF & Gustave Boël  Sofina Fellow 2019
Prof. Laurent Demanet
Development of a hybrid framework combining datadriven and modelbased methods for the study of corrosionerosion phenomena

FNRS Research Fellow2015  19
Université catholique de Louvain, Belgium  Ecole Centrale de Nantes, France
Prof. Roland Keunings (UCL) and Prof. Francisco Chinesta (ECN / ENSAM ParisTech)
PhD thesis: Multiscale modelling of fibre suspensions: Particle inertia, confined flows and datadriven approach
PhD thesis in cotutelle between UCL and ECN1112/2018: Research stay with Prof. A. Ammar at ENSAM Angers, France
0205/2019: Research stay at ESI Group Tokyo, Japan

Visiting FellowSummer 2014 / Summer 2015 / Summer 2016
Massachusetts Institute of Technology, USA
Prof. Laurent Demanet and Dr. Leonardo ZepedaNúñez
Implementation of a high performance solver for the Helmholtz equation with applications to seismic imaging and wave propagation.
Education

Ph.D. in Engineering and Technology2015  18
Université catholique de Louvain, Belgium  Ecole Centrale de Nantes, France

Master of Engineering (M.Eng.) in Applied Mathematics2013  15
Université catholique de Louvain, Belgium
Summa cum laude
Specialization in "Modelling and simulation of physical systems" and "Numerical algorithms" with a focus on mathematical modelling of complex fluids and flows.
Master thesis: "Multiscale modelling of dilute suspensions: from confined fibres to rigid clusters composed of rods"
Supervisors: Prof. Roland Keunings and Prof. Francisco Chinesta 
Bachelor of Applied Science  Engineering (B.Eng.)2010  13
Université catholique de Louvain, Belgium
Summa cum laude
Major: Applied Mathematics, minor: Computer Science
Programming Skills

Matlab

C

Python

Java
Languages

French

English

German
9 Research publications
Google Scholar record
Below is a list of publications sorted out by topic:
Feel free to contact me for more information concerning my past and current research!
Material on this web site is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. Copyright holders claiming that the material available below is not in accordance with copyright terms and constraints are invited to contact the author by email and ask him to remove the links to specific manuscripts.
PhD thesis
A. Scheuer
PDF
Multiscale modelling of fibre suspensions: particle inertia, confined flows and datadriven approach
Université catholique de Louvain / Ecole Centrale de Nantes, October 2018
ABSTRACT Suspensions of fibres and nonspherical particles are encountered in many fields ranging from engineering to biology, e.g. papermaking, composite manufacturing, pharmaceutical applications, red blood cells, foodprocessing and cosmetics industries, etc. Predicting the evolution of the orientation state of the particles is crucial to estimate the rheology of the suspension, that is its flow behaviour, as well as the final properties of the material. Jeffery's theory, describing the kinematics of a single particle immersed in an homogeneous flow of Newtonian fluid, lays the foundation for almost every models used today. Coarser representations, built upon this theory, have been introduced later to describe statistically the orientation state of the particles, either using a probability density function, or even moments of this function (AdvaniTucker orientation tensors). The assumptions underlying Jeffery's model are however quite restrictive to predict reliably what happens in fibre suspensions flows encountered in industrial processes. In this thesis, we first revisit this model, studying the impact of particle inertia and of confinement (wall effects) on the particle kinematics. In each case, we propose a multiscale approach, but given the challenges to upscale the microscopic description to the macroscopic scale, we then came up with an innovative approach based on datadriven simulations to circumvent upscaling issues and inaccuracies introduced by macroscopic closure approximations. Finally, we developed efficient numerical methods to simulate fluid flows in thin geometries, considering, within the Proper Generalized Decomposition (PGD) framework, an inplane/outofplane separated representation of the solutions of the incompressible NavierStokes equations.
Multiscale modelling of fibre suspensions
Orientation of a fibre in a shear flow  confined (red) vs. unconfined (blue)
M. Perez, A. Scheuer, E. AbissetChavanne, F. Chinesta, R. Keunings
PDF
A multiscale description of orientation in simple shear flows of confined rod suspensions
Journal of NonNewtonian Fluid Mechanics, 233, 6174 (2016)
ABSTRACT The multiscale description of dilute or semidilute suspensions involving rods has been successfully accomplished and applied in many scenarios of industrial interest. Many processes involve, however, the flow of rod suspensions in very narrow gaps whose thickness is much smaller than the rod length. In these conditions, the evolution of rod orientation is expected to be affected by confinement effects. In the present work, we propose a multiscale description of rod orientation in confined conditions and simple shear flows.
A. Scheuer, E. AbissetChavanne, F. Chinesta, R. Keunings
PDF
Secondgradient modelling of orientation development and rheology of dilute confined suspensions
Journal of NonNewtonian Fluid Mechanics, 237, 5464 (2016)
ABSTRACT We address the extension of Jeffery’s model, governing the orientation of rods im mersed in a Newtonian fluid, to confined regimes occurring when the thickness of the flow domain is narrower than the rod length. The main modelling ingredients concern: (i) the consideration of the rod interactions with one or both gap walls and their effects on the rod orientation kinematics; and (ii) the consideration of non uniform strain rates at the scale of the rod, requiring higherorder descriptions. Such scenarios are very close to those encountered in real composites forming processes and have never been appropriately addressed from a microstructural point of view. We also show that confinement conditions affect the rheology of the suspension.
R. Mezher, M. Perez, A. Scheuer, E. AbissetChavanne, F. Chinesta, R. Keunings
PDF
Analysis of the Folgar & Tucker model for concentrated fibre suspensions in unconfined and confined shear flows via direct numerical simulation
Composites Part A, 91, 388397 (2016)
ABSTRACT The classical Jeffery model allows for the prediction of the flowinduced orientation in dilute fibre suspensions. In most industrial applications, however, fibre suspensions are concentrated and fibrefibre interactions cannot be ignored any longer. These interactions have been traditionally modelled at the mesoscopic and macroscopic scales by introducing a phenomenological diffusion term inducing a randomizing effect. In the socalled Folgar & Tucker (F&T) model, widely used in applications, the diffusion coefficient is assumed to scale linearly with the flow intensity, the latter being described by the second invariant of the rate of strain tensor. Modifications and alternatives to the F&T model have been proposed in view of the difficulty for the F&T model to explain an apparent orientation delay observed experimentally in injectionmoulded parts. The noticed deviations were attributed to the intense fibrefibre interactions, thus pointing to the limitations of a phenomenological diffusion term for describing them. In the present work, we revisit the F&T model and compare its predictions with those obtained by simplified stateoftheart direct numerical simulation (DNS) in unconfined and confined simple shear flows for a range of shear rates and concentrations, the latter ensuring intense fibrefibre interactions. In unconfined flows, we find that the F&T model agrees quantitatively with the DNS results once an adequate closure relation is considered for approximating the fourthorder orientation tensor involved in the F&T model. Thus, the results seem to confirm, at least in simple shear flows, the F&T assumption for the form of the isotropic rotary diffusion function scaling linearly with the magnitude of the scalar rate of deformation. Also, a linear scaling of the diffusivity with the fibre concentration is observed. This conclusion remains unexpectedly valid under moderatelyconfined flow conditions as soon as an advanced fitted closure, like the IBOF, is considered within the F&T model. Other simpler closures (e.g. quadratic or hybrid), however, definitively fail to address confinement issues as also reported in our former work for the dilute regime. Obviously the validity of these conclusions depends on the validity of the considered stateoftheart DNS, issue that remains at present an open question.
R. Ibañez, A. Scheuer, E. Lopez, E. AbissetChavanne, F. Chinesta, R. Keunings
PDF
From elastic homogenization to upscaling of nonNewtonian fluid flows in porous media
International Journal of Material Forming, 11, 607617 (2018)
ABSTRACT Upscaling behaviors of heterogeneous microstructures to define macroscopic effective media is of major interest in many areas of computational mechanics, in particular those related to materials and processes engineering. In this paper, we explore the possibility of defining a macroscopic behavior manifold from microscopic calculations, and then use it directly for efficiently performing manifoldbased simulations at the macroscopic scale. We consider in this work upscaling of nonNewtonian flows in porous media, and more particularly the ones involving shortfibre suspensions.
A. Scheuer, E. AbissetChavanne, F. Chinesta, R. Keunings
PDF
Microscopic modelling of orientation kinematics of nonspherical particles suspended in confined flows using unilateral mechanics
Comptes Rendus Mécanique, 346, 4856 (2018)
ABSTRACT The properties of reinforced polymers strongly depend on the microstructural state, that is, the orientation state of the fibres suspended in the polymeric matrix, induced by the forming process. Understanding flowinduced anisotropy is thus a key element to optimize both materials and process. Despite the important progresses accomplished in the modelling and simulation of suspensions, few works addressed the fact that usual processing flows evolve in confined configurations, where particles characteristic lengths may be greater than the thickness of the narrow gaps in which the flow takes place. In those circumstances, orientation kinematics models proposed for unconfined flows must be extended to the confined case. In this short communication, we propose an alternative modelling framework based on the use of unilateral mechanics, consequently exhibiting a clear analogy with plasticity and contact mechanics. This framework allows us to revisit the motion of confined particles in Newtonian and nonNewtonian matrices. We also prove that the confined kinematics provided by this model are identical to those derived from microstructural approaches (Perez et al. (2016)).
A. Scheuer, A. Ammar, E. AbissetChavanne, E. Cueto, F. Chinesta, R. Keunings, S.G. Advani
PDF
Datadriven upscaling of orientation kinematics in suspensions of rigid fibres
Computer Modeling in Engineering & Sciences, 117, 367386 (2018)
ABSTRACT Describing the orientation state of the particles is often critical in fibre suspension applications. Macroscopic descriptors, the socalled secondorder orientation tensor (or moment) leading the way, are often preferred due to their low computational cost. Closure problems however arise when evolution equations for the moments are derived from the orientation distribution functions and the impact of the chosen closure is often unpredictable. In this work, our aim is to provide macroscopic simulations of orientation that are cheap, accurate and closurefree. To this end, we propose an innovative databased approach to the upscaling of orientation kinematics in the context of fibre suspensions. Since the physics at the microscopic scale can be modelled reasonably enough, the idea is to conduct accurate offline direct numerical simulations at that scale and to extract the corresponding macroscopic descriptors in order to build a database of scenarios. During the online stage, the macroscopic descriptors can then be updated quickly by combining adequately the items from the database instead of relying on an imprecise macroscopic model. This methodology is presented in the wellknown case of dilute fibre suspensions (where it can be compared against closurebased macroscopic models) and in the case of suspensions of confined or electricallycharged fibres, for which stateoftheart closures proved to be inadequate or simply do not exist.
M. Perez, A. Scheuer, E. AbissetChavanne, A. Ammar, F. Chinesta, R. Keunings
PDF
On the multiscale description of microstructured fluids composed of aggregating rods
Continuum Mechanics and Thermodynamics, 31, 955967 (2019)
ABSTRACT When addressing the flow of concentrated suspensions composed of rods, dense clusters are observed. Thus, the adequate modelling and simulation of such a flow requires addressing the kinematics of these dense clusters and their impact on the flow in which they are immersed. In a former work, we addressed a first modelling framework of these clusters, assumed so dense that they were considered rigid and their kinematics (flowinduced rotation) were totally defined by a symmetric tensor $\mathbf c$ with unit trace representing the cluster conformation. Then, the rigid nature of the clusters was relaxed, assuming them deformable, and a model giving the evolution of both the cluster shape and its microstructural orientation descriptor (the socalled shape and orientation tensors) was proposed. This paper compares the predictions coming from those models with finerscale discrete simulations inspired from molecular dynamics modelling.
A. Scheuer, G. Gregoire, E. AbissetChavanne, F. Chinesta, R. Keunings
PDF
Modelling the effect of particle inertia on the orientation kinematics of fibres and spheroids immersed in a simple shear flow
Computers and Mathematics with Applications, In press (2019)
ABSTRACT Simulations of flows containing nonspherical particles (fibres or ellipsoids) rely on the knowl edge of the equation governing the particle motion in the flow. Most models used nowadays are based on the pioneering work of Jeffery (1922), who obtained an equation for the motion of an ellipsoidal particle immersed in a Newtonian fluid, despite the fact that this model relies on strong assumptions: negligible inertia, unconfined flow, dilute regime, flow unperturbed by the presence of the suspended particle, etc. In this work, we propose a dumbbellbased model aimed to describe the motion of an inertial fibre or ellipsoid suspended in a Newtonian fluid. We then use this model to study the orientation kinematics of such particle in a linear shear flow and compare it to the inertialess case. In the case of fibres, we observe the appearance of periodic orbits (whereas inertialess fibres just align in the flow field). For spheroids, our model predicts an orbit drift towards the flowgradient plane, either gradually (slight inertia) or by first rotating around a moving oblique axis (heavy particles). MultiParticle Collision Dynamics (MPCD) simulations were carried out to assess the model predictions in the case of inertial fibres and revealed similar behaviours.
3D Helmholtz equation
Solution of the Helmholtz equation for the SEAM model
L. Zepeda Núñez, A. Scheuer, R. J. Hewett, L. Demanet
The method of polarized traces for the 3D Helmholtz equation
Geophysics, 84, 313333 (2019)
ABSTRACT We present a fast solver for the 3D Helmholtz equation, in heterogeneous, constant density, acoustic media, in the highfrequency regime. The solver is based on the method of polarized traces, a layered domaindecomposition method, where the subdomains are connected via transmission conditions prescribed by the discrete Green's representation formula and artificial reflections are avoided by enforcing nonreflecting boundary conditions between layers. The method of polarized traces allows us to consider only unknowns at the layer interfaces, reducing the overall cost and memory footprint of the solver. We demonstrate that polarizing the wavefields in this manner yields an efficient preconditioner for the reduced system, whose rate of convergence is \textit{independent} of the problem frequency. The resulting preconditioned system is solved iteratively using GMRES, where we never assemble the reduced system or preconditioner, rather we implement them via solving the Helmholtz equation locally within the subdomains. The method is parallelized using MPI and coupled with a distributed linear algebra library and pipelining to obtain an empirical online runtime $\cO(\max(1,R/L) N \log N)$ where $N = n^3$ is the total number of degrees of freedom, $L$ is the number of subdomains, and $R$ is the number of righthand sides. This scaling is favorable for regimes in which the number of sources (distinct righthand sides) is large, for example enabling largescale implementations of \LZ{frequencydomain} full waveform inversion (FWI).
Unclassified
E. Lopez, A. Scheuer, E. AbissetChavanne, F. Chinesta
PDF
On the effect of phase transition on the manifold dimensionality: Application to the Ising model
Mathematics and Mechanics of Complex Systems, 6, 251265 (2018)
ABSTRACT Fields can be represented in a discrete manner from their values at some locations, the nodes when considering finite elements descriptions. Thus, each discrete scalar solution can be considered as a point in $\mathbb R^N$ ($N$ being the number of nodes used for approximating the scalar field). Most manifold learning techniques (linear and nonlinear) are based on the fact that those solutions define a slow manifold of dimension $n \ll N$ embedded in the space $\mathbb R^N$. This paper explores such a behaviour in systems exhibiting phase transitions in order to analyse the evolution of the local dimensionality $n$ when the system moves from one side of the critical behaviour to the other. For that purpose we consider the Ising model.
R. Ibañez, A. Scheuer, E. AbissetChavanne, F. Chinesta, A. Huerta, R. Keunings
PDF
A simple microstructural viscoelastic model for flowing foams
International Journal of Material Forming, 12, 295306 (2019)
ABSTRACT The numerical modelling of forming processes involving the flow of foams requires taking into account the different problem scales. Thus, in industrial applications a macroscopic approach is suitable, whereas the macroscopic flow parameters depend on the cellular structure: cell size, shape, orientation, etc. Moreover, the shape and orientation of the cells are induced by the flow. A fully microscopic description remains useful to understand the foam behaviour and the topological changes induced by the cell elongation or distortion, however, from an industrial point of view, microscopic simulations remain challenging to address practical applications involving flows in complex 3D geometries. In this paper, we propose a viscoelastic flow model where the foam microstructure is represented from suitable microstructure descriptors whose evolution is governed by the macroscopic flow kinematics.
s Music
I enjoy playing the cello and am a former member of the LouvainlaNeuve Student Symphony Orchestra (OSEL).
` Contact info
 Massachusetts Institute of Technology
MIT Mathematics / Earth Resources Laboratory
Green Building  Office 54323
77 Massachusetts Avenue
02139 Cambridge, MA
United States  Email: ascheuer [at] mit.edu
 http://math.mit.edu/~ascheuer/