EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

ECCOMAS NEWSLETTER

DECEMBER 2012

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

CONTENTS FORWARD...... 3 MANOLIS PAPADRAKAKIS

ECCOMAS ΜEDALS AND AWARDS 2012 ...... 4 MANOLIS PAPADRAKAKIS

ECCOMAS 2012 CONGRESS REPORT……...... 7 JOSEF EBERHARDSTEINER

INTERNATIONAL RESEARCH TRAINING GROUPS, A FRUITFUL MODEL OF COLLABORATION FOR PHD STUDENTS …...... 9 PETER WRIGGERS, OLIVIER ALLIX

COMPUTATIONAL MODELING OF DAMAGE IN CONCRETE DAMS DUE TO ALKALI-AGGREGATE REACTION...... 13 CLAUDIA COMI, UMBERTO PEREGO, R. FEDELE, R. PIGNATELLI

HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES FOR COMPRESSIBLE FLOWS ON UNSTRUCTURED GRIDS...... 18 REMI ABGRALL, DANTE DE SANTIS, MARIO RICCHIUTO

ECCOMAS THEMATIC AND REGIONAL CONFERENCES 2013...... 22

ND 2 ECCOMAS YOUNG INVESTIGATORS CONFERENCE 2013 - PHD OLYMPIAD 2013…………………….23

MEMBER ASSOCIATIONS FORUM…………………………………………………………………………………………….……...... 24

ECCOMAS CONGRESS 2016…………………………………………………………………………….………………...... 25 TH TH 5 EUROPEAN ON COMPUTATIONAL MECHANIC - 6 EUROPEAN ON ...... 26

ALL RIGHTS RESERVED BY AUTHORS. ECCOMAS NEWSLETTER — DECEMBER 2012

EDITORS TECHNICAL EDITOR

MANOLIS PAPADRAKAKIS EKKEHARD RAMM NATIONAL TECHNICAL UNIVERSITY OF ATHENS UNIVERSITÄT STUTTGART PANAGIOTA KOUNIAKI PRESIDENT OF ECCOMAS VICE PRESIDENT OF ECCOMAS NATIONAL TECHNICAL UNIVERSITY OF ATHENS

PEKKA NEITTAANMÄKI PEDRO DÍEZ UNIVERSITY OF JYVÄSKYLÄ UNIVERSITAT POLITÉCNICA DE CATALUNYA VICE PRESIDENT OF ECCOMAS SECRETARY OF ECCOMAS

2

FORWARD

2012 was undoubtedly an extremely The Newsletter also contains a schemes can be successfully applied successful year for ECCOMAS, which presentation by PETER WRIGGERS to perform the simulation of was marked by the VI European and OLIVIER ALLIX on a fruitful complex flows, like three Congress on Computational collaboration model of PhD students dimensional turbulent flows. Methods in Applied Sciences and under the initiative of the Despite the availability of powerful (ECCOMAS 2012), International Research Training computers, which has made hosted on the premises of the Groups (IRTG) for providing possible to perform simulations University of Vienna, , opportunities for joint doctoral with an extremely large numbers of September 10-14, 2012, and the degrees of freedom, the problem of th training programs between German celebration of the 20 Anniversary Universities and Universities from achieving accurate simulations of of the creation of ECCOMAS. The foreign countries. The research and complex flows has been addressed Congress was a great success thanks study programs are jointly in the last years with a new strategy to the thoughtful and tireless efforts developed and supervised by the which is reviewed in this article. of the organizers. At his stage I partner universities, while the Instead of using finer grids, higher would like to take the opportunity doctoral students in the program order methods are implemented to express my deepest appreciation complete a six-month research stay (pages18-21). to the Chairmen of the ECCOMAS at the respective partner 2012 Congress, HERBERT MANG, institutions (pages 9-12). A list of approved ECCOMAS JOSEF EBERHARDSTEINER, HELMUT Thematic Conferences and Special Interest Conferences for 2013 is OEHM and RANZ LAUDIA OMI and MBERTO B F C C U provided on page 22. Announce- RAMMERSTORFER for their PEREGO, in collaboration with R. ments for the forthcoming Young impeccable and professional FEDELE and R. PIGNATTELLI, report Investigators Conference (YIC 2013) organization. Last but not least, I on the computational modeling of and the ECCOMAS Olympiad 2013 would like to thank the Secretary damage in concrete dams due to can be found on page 23. General, MARTINA POELL, the heart alkali-aggregate reaction which Announcements of the joint and soul of the Congress, who plays a major role for the safety of a ECCOMAS Conferences: The 5th contributed immensely to its large number of dams built all over European Conference on Computa- success. The Opening Ceremony the world. In the ICOLD tional (ECCM V) and the and the Plenary lectures took place (International Commission of Large 6th European Conference on in the Golden Hall of the famous Dams) more than one hundred large Computational Fluid Dynamics Vienna’s Musikverein where the concrete dams have been already (ECFD VI) to be held in Barcelona, ECCOMAS Medals and Awards for identified as being seriously affected Spain in July 2013 can be found on 2012 were also conferred. by the structural consequences of page 26, while the announcement concrete chemical expansion, due to of the next ECCOMAS Congress in This electronic edition of the the Alkali-Aggregate Reduction 2016, which was decided to take ECCOMAS Newsletter contains the (AAR). The authors discuss reliable place in Crete, Greece in June 2016, presentation of the winners of the mathematical models for the can be found on page 25. ECCOMAS Medals and Awards for assessments of the dam structural 2012 as well as the winners of the safety and the design of a possible two best PhD Theses for 2011 in the rehabilitation actions (pages 13-17). field of Computational Methods in

Applied Sciences and Engineering REMI ABGRALL, DANTE DE SANTIS

(pages 4-6), while a Report on the and MARIO RUCCHINTO present ECCOMAS 2012 Congress is recent achievements in high order MANOLIS PAPADRAKAKIS provided by the co-Chairman of the residual distribution schemes for PRESIDENT OF ECCOMAS Congress, JOSEF EBERHARDSTEINER compressible flows on unstructured [email protected] (pages 7-8). grids, showing that such kind of

3

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

ECCOMAS MEDALS AND AWARDS FOR 2012

The ECCOMAS Awards for 2012 The ECCOMAS winner of the best Award in the field of Computaional were conferred during the opening PhD Thesis for 2011 in the field of . The awards are given ceremony of the Sixth ECCOMAS Computational Solids and every two years at the ECCOMAS Congress on Computational Structures was LENA WIECHERT- Congress or at the ECCOMAS Methods in Applied Sciences and YOSHIHARA for the thesis entitled: Computational Solid and Structural Engineering, in the Golden Hall of “Computational Modeling of Multi- Mechanics and Fluid Dynamics the famous Vienna’s Musikverein Conferences, respectively and are th Field and Multi-Scale Phenomena in on September 10 , 2012. The three Respiratory Mechanics” under the accompanied with a prize of 2.000 Medals to honor eminent supervision of Wolfgang Wall, Euros. researchers: the Ritz-Galerkin conducted at the Technical Medal and the Euler and Prandtl University of Munich, . The ECCOMAS Zienkiewicz Award Medals and the two Awards: the for young scientists in Computation- Lions and Zinkiewicz Awards for The ECCOMAS winner of the best al Engineering Sciences is awarded young investigators. Moreover, the Ph.D. Thesis for 2011 in the field of ex aequo to ELÍAS CUETO and two best PhD Theses for 2011, in Computational Fluid Dynamics was ALESSANDRO REALI. the field of Computational Methods MAX DUARTE for the thesis in Applied Sciences and Elías Cueto obtained his BSc and Engineering, were also announced entitled: ”Αdaptive numerical methods for the simulation of time- MSc in Mechanical Engineering at the opening ceremony. from the University of Oviedo, space multi-scale reacting fronts» under the supervision of Marc Spain and his PhD from the For the two Ph.D. Awards, every University of Zaragoza, on the Member Association proposes one Massot, Frédérique Laurent- Nègre and Stéphane Descombes, natural element method. He is now candidate and the ECCOMAS a Professor of Mechanical - Awards Committee selects two conducted at the Ecole Centrale de Paris, France. ing Applied Mechanics and awardees, one in the field of Bioengineering at the University of “Computational Solids and Both PhD Awards were conferred at Zaragoza and holds also a visiting Structures” and the other in field of professorship at the Ecole Centrale “Computational Fluid Dynamics”. the First Young Investigators Conference, in Aveiro, Portugal de Nantes in France. He is an active The presentation of the awards researcher in a number of topics in takes place annually at an earlier this year. computational mechanics including appropriate thematic event and are meshless methods with emphasis accompanied with a prize of 2.000 The two ECCOMAS Awards for young scientists are the Zienkiewicz on the simulation of forming Euros. processes, model reduction Award in the field of Computational Engineering Sciences and the Lions methods by means of Proper Generalized Decompositions and its

application to real-time simulation in biomedical applications.

Alessandro Reali received his MSc on locking-free mixed-enhanced

finite elements and his PhD on constitutive modeling of shape memory alloys, both from the Lena Wiechert-Yoshihara Max Duarte ECCOMAS winner of the best PhD Thesis for ECCOMAS winner of the best PhD Thesis for University of Pavia, Italy. He is now 2011 in the field of Computational Solids 2011 in the field of Computational Fluid an Associate Professor in the and Structures Dynamics

4

Department of Civil Engineering and and transitional molecular/ methods for dynamic systems, but Architecture, at the University of continuum flows. soon he was devoted to the solution Pavia. His main research topic is of computational fluid dynamics isogeometric analysis, a field he The Euler Medal for outstanding and problems. Since then, his research started to explore in collaboration sustained contribution to the area of activities have spread over a range of with T.J.R. Hughes and his group. He computational solids and structural multidisciplinary fields which he has has made substantial contributions mechanics and coupled problems in contributed relevant theories and to IGA in many fields among which, engineering and the Prandtl Medal methods of scientific and industrial structural vibrations, wave for outstanding and sustained relevance. He was one of the first propagation, structural problems, contribution to the area of who applied Meshless Methods to turbulence in fluids, collocation computational fluid dynamics, are fluid dynamics and in particular the methods. He is Recipient of the ERC conferred every two years at the finite point method in convective “Ideas” Starting Independent ECCOMAS main Conferences. transport and fluid flow. He was the Researcher Grant. main pioneer of the development of Ludwig Prandtl was a pioneer Particle Methods using a Finite The ECCOMAS Lions Award for German scientist who developed a Element mesh. These methods are young scientists in Computational rigorous mathematical analyses revolutionary for the computational Mathematics is awarded to HARALD which he implemented for solutions of VAN BRUMMELEN. establishing the basis of the applied problems with free-surfaces, science of aeronautical engineering. breaking waves, multi-fluids and Harald van Brummelen obtained his In the 1920s he developed the fragmentation problems. In addition MSc in mechanical engineering from mathematical basis for the to his scientific activity, he has Twente University of Technology and fundamental principles of subsonic developed an intensive activity in the his PhD in from aerodynamics, including transonic transfer of the outcome of his the University of Amsterdam and he velocities. His studies identified the research to the industrial sector in is currently a Professor at Eindhoven boundary layer, thin-airfoils, and Europe and in Latin America. University of Technology where he lifting-line theories. He is often holds the chair on Multiscale referred to as the father of modern Leonhard Euler was an outstanding Engineering Fluid Dynamics at the aerodynamics and founder of Swiss and physicist. Department of Mechanical modern Fluid Mechanics. The He made important discoveries in Engineering. He is the Director of the Prandtl number was named after fields as diverse as infinitesimal Eindhoven Multiscale Institute and him, while the crater Prandl on the calculus and graph theory. He is the vice-chairman of the ECCOMAS far side of the Moon is named in his father of much of the modern Committee on Computational and honor. Theodor von Karman was one mathematical terminology and Applied Mathematics. His research of his notable students. The winner notation and is considered to be the focuses on the development and of Ludwig Prandtl Medal is SERGIO preeminent mathematician of the analysis of advanced numerical IDELSOHN. 18th century, and one of the techniques for multiscale flow greatest of all time. He is also one of problems, with emphasis on free- Sergio Idelsohn started his research the most prolific boundary and interface problems, endeavors on the finite element ever. fluid-structure interaction problems analysis of shells and reduced order A statement attributed to Pierre- Simon Laplace, quote "Read Euler, read Euler, he is the master of us all", expresses Euler's influence on mathematics. In addition to mathematics he is also renowned for his work in mechanics and particularly in stability theory. The winner of Leonhard Euler Award is HERBERT MANG. Elías Cueto and Alessandro Reali Harald van Brummelen ECCOMAS ex aequo winners of the ECCOMAS winner of the Lions Award for Zienkiewicz Award for young scientists in young scientists in Computational Computational Engineering Sciences Mathematics

5

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

Sergio Idelsohn Herbert Mang Erwin Stein Winner of the ECCOMAS Prandtl Medal Winner of the ECCOMAS Euler Medal Winner of the ECCOMAS Ritz-Galerkin Medal Herbert Mang has devoted a lot of his disciplines. The medal is conferred encountered in engineering practice. scientific potential in computational every four years at the ECCOMAS His fields of activity are extensive and stability analysis throughout his Congress. different. His scientific achievements professional career. The awardee is a are reflected in numerous scholar of highest distinction and Walter Ritz, who was born a few publications on physical and mathe- personality and definitely is one of hundred kilometers from Vienna in matical modelling, mathematical and the most prominent scientists in our southern , began his , numerical community. His research covers a studies at ETH Zurich in engineering and software engineering broad area in structural and solid where he soon found that he could in continuum and structural mechanics including modeling, not live with the approximations and mechanics. His work has been computational aspects and compromises involved with honoured by a number of applications in engineering practice. engineering of that period I would organizations and universities both His trademark is his ability to adjust say, and switched to more national and international. In addition and apply advanced computational mathematically exact studies in to his research achievements he had a methods, deeply rooted in a scientific physics, with Albert Einstein being remarkable contribution in inspiring basis, to real life applications. He is one of his classmates. Boris Galerkin, his students to pursue research in the initiator of a ground-breaking new who was born in , part of the Computational Mechanics and research direction in computational at that time, studied becoming a mentor for a large structural stability aimed at improving at St. Petersburg Technological number of doctoral students, many of the postbuckling behavior of thin Institute mathematics. But as he them became university professors. elastic shells by small changes of the needed funds to survive, he then The leitmotivs of Leibniz “theoria cum original design. His activities are studied engineering and worked for praxis” and “commune bonum” is documented in more than 20 books many years as a structural engineer guiding his journey through his entire and over 440 papers in scientific before entering the academic life. Ritz life, while the motto used for the journals and books. He is a member and Galerkin common characteristic promotion of the Leibniz exhibition: of several National Academies of was that they both combined quote “Gottfried Wilhelm Leibniz – Sciences or Engineering and is mathematical theory and engineering way ahead of his time – Philosopher, decorated with many awards and practice to develop independently, Mathematician, Physicist and recognitions for his work. among other, the Ritz-Galerkin Engineer” could also be applied to variational approximation methods of him. The Ritz-Galerkin Medal is the highest boundary and initial value problems. Award given by ECCOMAS. It honours The winner of the Ritz-Galerkin Award On behalf of ECCOMAS my warmest individuals who have made is ERWIN STEIN. congratulations to all recipients of the outstanding and sustained 2012 ECCOMAS Awards and Medals. contributions in the field of Erwin stein studied engineering and computational methods in applied mathematics and has a long career sciences and engineering. The medal both as a researcher and as a carries the images of Ritz and Galerkin structural engineer. Being a strong MANOLIS PAPADRAKAKIS in recognition of the synergy between advocate of the interaction of RESIDENT OF mathematics, numerical analysis, P ECCOMAS Mathematics and Mechanics he mechanics of solids and structures, implemented this school of thought [email protected] fluid dynamics, and other engineering for solving real-world problems 6

ECCOMAS 2012 CONGRESS REPORT

7

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

JOSEF EBERHARDSTEINER CO-CHAIRMEN OF ECCOMAS CONGRESS

[email protected]

8

International Research Training Groups a FruitfulM odel of Collaboration for PHD Students

The funding of International Research Foundation) while the French collaboration in engineering sciences. Training Groups (IRTG) is an initiative partners used their own funds to launched by DFG in 1997 in order to support the scientific exchange. The Several procedures were installed to provide opportunities for joint training group is a collaboration of ensure that the goals of the IRTG doctoral training programs between researchers representing faculties of 1627 are in the core of its activities. A German universities and universities civil and mechanical engineering as strong cooperation for years with the from foreign countries. The research well as mathematics at LUH and ENS LMT of ENS Cachan, established by and study programs are jointly Cachan. Prof. P. Wriggers, was instrumental developed and supervised by the for the establishment of the IRTG’s partner universities. Doctoral The training group started October functional structure and its working students in the program complete a 2010 with 12 German PhD students. program for. The training group is six-month research stay at the Later 6 PhD students from ENS managed by P. Wriggers and O. Allix, respective partner institution. The Cachan joined the research training from the Institute of Continuum key factors of the evaluation of the group that is now fully operational Mechanics at Leibniz Universität IRTG program by an international having in total 20 principal Hannover and the LMT at ENS committee are: investigators from both universities. Cachan, respectively. The PhD  Excellent, coherent and multidis- The seminars organized every year students are supervised using a ciplinary research program are open to other PhD student of the double coaching concept. Within this  Specific qualification program two universities which enables the approach a German and a French PI  Innovative supervision concept student to have a large view on are jointly responsible for assisting  Promotion of early independence research activities on both sides. the students and for evaluating and international mobility Even though international success and progress. As part of the cooperation and multidisciplinary educational program each of the In 2010 a multidisciplinary group of work is pursued in many aspects and German PhD students e.g. has to scientists from the École Normale at different levels in the scientific spend a minimum of 6 months of his/ Supérieure (ENS) Cachan, France, and community, the group is one of the her training period at ENS Cachan to the Leibniz Universität Hannover few research training groups that work on their project. (LUH), Germany, applied successfully establishes a true European for a grant of an international research training group (IRTG 1627) with the purpose to jointly educate PhD students. The research topics and general scientific vision of the training group is “Virtual Materials and Structures and their Validation” which has wide spread applications in many areas of engineering. The German partners received a grant from the DFG (German Science

Official starting of the IRTG Hanover September 2010

9

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

Typical working atmosphere of "our meetings” Participants of the student's workshop in Mehringen

At the heart of the training group are Dmytro Rodman from the Institute of quenching conditions on pitting the PHD students who are supposed Material Science in Mechanical resistance. Staying at the LMT new to actively take part in the Engineering provided the following research approaches and modelling organization of the IRTG. Hence they statement: methods were learned in a friendly determine the best period for their “In November 2012 I completed my working atmosphere. The everyday stay abroad, they organize seminars research stay of six month at the LMT contact with colleagues from the and workshops for which they have Cachan. During this time he worked LMT helped to improve the French the opportunity to invite senior with Prof. Sylvie Pommier on the skills and to gain insight in the way of researchers to give a course on a modelling of the fatigue strength life in France.” given subject of interest. We asked behaviour of induction hardened two of the PhD students from the gearwheels by employing water-air André Hürkamp from the Institute of Leibniz Universität Hannover, who spray cooling. The purpose of the Mechanic and Computational already stayed at the LMT, to give a cooperation was the development of Mechanic in Civil Engineering short report of their involvement in a numerical model regarding the provided insights with respect to his the research training group and superposition of contact stresses due scientific interests: about their experiences during their to loading and residual stresses due “As a member of the IRTG 1627, I had stay at the LMT in Cachan. to heat treatment in order to the opportunity to participate in the determine the influence of different research exchange program between the Leibniz Universität Hannover and the LMT-Cachan near Paris. I lived and worked in Cachan from May 7th to November 30th. During that period, I collaborated with Christophe Cluzel of the Department of Structures and Systems. Together we worked on the modeling of damage behavior of short fiber reinforced concrete. The interdisciplinary structure of the LMT and the pleasant working atmosphere made it very easy to get in contact with researchers from different fields. This provided not only an insight into topics out of the own field of research but was also very fruitful for myself. Altogether, It has been a Photo of the participants to the Porquerolles workshop

10

Vyacheslav Boyarkin (National Metallurgical Academy of Ukraine (NMetAU), Dnepropetrovsk) and Prof. Dr. Anton Krysko (Saratov State Technical University, Russia) held courses or presented their highly appreciated research results. In the meantime, the close contacts that were established between all scientists lead to several exchange visits abroad by PhD students. With this the PhD students were able to broaden their international Seminar of IRTG 1627 in Mehringen experience on top of the “basic” stay at ENS Cachan. valuable experience to stay for A series of special seminars, almost 7 month in this research organized by the IRTG 1627, assisted With respect to 2013, it is environment and also in the to bulid up special knowledge. Well- anticipated that the IRTG 1627 is metropolis region of Paris.” known researchers like Prof. Joze well prepared to meet its goals, and Korelc, (University of Ljubljana, that the first cohort of PhD students Different activities of the group in Slovenia), Prof. Laura de Lorenzis will successfully finish their the past were organized such that (University of Salento, Lecce), Prof. dissertation. the group could grow together and Marco Paggi (University of Torino), an effective collaboration between Dr. Guillaume Pacquaut (Université Of course such cooperation is very all members could be established. de Lyon, INSA-Lyon), Prof. Daya demanding from the PhD point of Several meetings and workshops Reddy (South African Research Chair view but it is our feeling that they took place, e.g. on the island of in Computational Mechanics, benefit a lot of this early Porquerolles (Hyères/France), and at University of Cape Town), Dr.-Ing. international experience. In the past the island of Wangerooge (Germany), where some PhD tried to symbolize both the pleasure and the difficult equilibrium that characterize an IRTG.

At these workshops the PhD students presented the current status of their work. Thus it was possible to find for each PhD student the best possible co-advisor from the partner university. Furthermore topics of interest such as the importance of the validation of computational models, general approaches to multiscale analyses, material models for composites, modeling and simulation of coupled problems and optimization of materials as well as modeling and simulation of uncertainties have been discussed. These are of great interest for the group of PhD students who work on the topics summarized in Box 1.

11

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

Majeed Bishara, M.Sc. Microbuckling and Fiber Kinking in Composites Using Multilevel Simulation

Ing. Matteo Capaldo PGD Approximation for Virtual Structural Testing

Ing. Amaury Courard Virtual Chart Concept in Computational Structural Mechanics

Nicolas Courrier Coupling partially converged data and a multiparametric strategy for the optimisation of assemblies

Dipl.-Ing. Daniel Gottschalk Multiscale modelling and characterization of heterogeneous media exhibiting size effects

Dipl.-Math. Barbara Development of a friction approach for the FE method of sheet metal forming based on Homann multi-scale modelling

Dipl.-Ing. Malte Hothan Development of non-convex shaped particles for DEM – A Multiscale Approach to Granu- lar Materials

André Hürkamp, M.Sc. Multi-Scale modelling of damage in laminated composites

Dipl.-Math. A. Issaoui Stochastic BEM and Multiscale-Modeling

Andreas Kempe, M.Sc. Virtual Development of Nanocomposites using the Molecular-Dynamic Dipl.-Ing. Jocelyn Lopez Towards robust and accurate simulation of composites shock-absorbers

Daniel Nolte, M.Sc. Three-dimensional XFEM analysis for cracks in composite materials

Ing. Enrico Obert Damage modeling of woven composite materials on the micro scale and the meso scale

Dipl.-Wirtsch.-Ing. Dmytro Application of spray cooling for surface hardening Rodman Guillaume Rosenzveig Modelisation and simulation in structures computation for tolerancing and assistance in assembling of large and flexible aeronautical parts

Dipl.-Ing.(FH) Ole Stegen Multi-scale constitutive modeling of carbon black reinforced rubber in a finite strain framework for FE analysis

Dipl.-Ing.(FH) Martin Enhancement of magnesium structures to support cardiovascular tissue scaffolds in the Weidling, M.Sc. high-pressure system

Michael Wenzel, M.Sc. Non-Newtonian Fluids in Fractured Media

Box1: Current PhD students and topics of IRTG 1627

such experience was often acquired later in a career For more details and further information see: and initiative like this and others like the newly http://www.irtg1627.uni-hannover.de launched ECCOMAS Young Investigators Conferences, will have a beneficial impact on international cooperation very soon. Additionally we would like to PETER WRIGGERS stress that this experience is positive not only for the PhD students but also for the involved senior OLIVIER ALLIX researchers, as was expressed by all the members of the IRTG during our annual workshop. [email protected]

12

COMPUTATIONAL MODELING OF DAMAGE IN CONCRETE DAMS DUE TO ALKALI-AGGREGATE REACTION

INTRODUCTION The concrete volume overall diagnosis and control and for expansion, may lead to macroscopic practical methods for determination The ICOLD (International Commission structural deformation preventing of model parameters. on Large Dams) Register of Dams the correct service operation. currently indexes 37641 large dams Moreover, selfstresses induced by (i.e. more than 15 meters high) differential swelling in the dam body EVIDENCES OF DAM through-out the world, out of which easily exceed the concrete strength, about nine thousands are concrete causing the propagation of BEHAVIOR dams. Among these, the majority of macroscopic cracks, usually visible Several case histories have been dams located in the more advanced from inspection galleries. In some extensively documented in the litera- countries have been built decades cases, selfstresses arising from ture. Typical reported macroscopic ago and more than one hundred of structurally constrained expansion evidences are: increase in crest them have been already identified as may cause indirect macroscopic elevation (from 5 μm/m/year to 150 being seriously affected by the cracks, typically in correspondence of μm/m/year); usually not uniform structural consequences of concrete dam geometry discontinuities. expansion in the vertical direction, chemical expansion, due in most this being higher in the dam upper cases to the so called Alkali- Even though, in view of their size and part; in arch dams, long term Aggregate Reaction (AAR). economical impact, dams are upstream drift, produced by volume probably the concrete structures increase (in gravity dams this is At the micro-scale, AAR is a chemical which have attracted more attention usually less visible, since it can be reaction which occurs in concrete be- from the scientific community, the produced only by non-uniform tween the highly alkaline cement problem of AAR is of relevance also through-the-thickness expansion or paste and non-crystalline silicon di- for other types of structures, such by a particular dam geometry); visible oxide (silica), which is found in many e.g. concrete bridges. macroscopic cracks in inspection gal- common aggregates and/or dissolved leries or at geometry discontinuities. within aggregates. The term AAR Expensive rehabilitation interven- Figure 1a shows the Pian Telessio includes three types of reaction: the tions on affected dams, aimed at arch dam, built in 1955 in Italy, where most common alkali-silica reaction, recreating suitable service conditions, an expansion rate of 10 μm/m/year the alkali-silicate reaction and the have been undertaken in many cases. was estimated at the dam top, less frequent alkali-carbonate These interventions are usually producing a considerable upstream reaction. The main product of the finalized to extend the dam service drift, clearly visible in Figure 1b, reaction between alkali and silica is a life well beyond the initial design du- reaching 55 mm at the center of the gel which, in the presence of water, ration of few decades. The correct crest in 2005, when rehabilitation expands initially filling up the pre- assessment of the dam structural works began. Several horizontal existing concrete pores. Further safety and the design of a possible cracks were also visible in the upper expansion may induce significant rehabilitation action require there- inspection gallery *1+. stresses within the concrete matrix, fore the possibility to predict the dam causing local micro-cracking with long term behavior under the effects consequent degradation of elastic of AAR, which in turn demands for properties and material strength. reliable mathematical models for

13

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

Fig. a) Rehabilitation works at Pian Telessio dam (after *2+); b) radial crest displacement (negative values denote upstream direction) (after *1+).

PHYSICAL EVIDENCES FROM the amount of available water, the MATHEMATICAL MODELS greater is the final volumetric LABORATORY TESTS expansion and the faster is the AND APPLICATIONS A necessary requisite for gel reaction; Mathematical models for the production and growth is the the AAR expansion may lead to simulation of the structural availability of reactive ingredients. micro-cracking and to a consequences of AAR need to These are the reactive aggregates, deterioration of Young’s modulus consider the interaction between the alkalis pre-sent in the cement and peak stress, even with the concrete matrix and the paste, and the water present in vanishing macroscopic stresses; for expanding gel. This is a concrete pores. When one of the the range of compressive stresses phenomenon evolving in time, ingredients is depleted, the reaction considered in the experiments, the governed by the kinetics of the comes to a rest. Recent studies chemical kinetics is not affected by chemical reaction, which strongly have however shown that in many the stress state; however, directions depends on cur-rent temperature cases alkalis become available with and intensity of local principal and moisture availability. Different time from certain types of reactive stresses strongly influence the models have been presented in the aggregates. In these cases the expansion redistribution within the literature at various levels of reaction can continue also after material. complexity. The simplest approach, depletion of the alkalis contained in consisting in modeling the gel the cement paste.

Several experimental results on AAR are available in the literature *3-5+, providing the basic understanding of the reaction mechanisms required for the formulation of a mathematical model. The main outcomes are: at the macro-scale, the AAR determines an expansion of the concrete mass, growing until a condition of chemical equilibrium is reached; the rate of the chemical reaction significantly depends on environmental conditions, such as temperature and moisture. In particular: the higher is the temperature, the faster is the reaction development; the larger is Figure 2. Schema of multiphase models hierarchy: Vs= skeleton volume, Vv=pore volume, Vg= gel volume, Vw= water volume, Va= pore air volume .

14

expansion as a temperature-like Since water is consumed in the and 1, expressing the degree of deformation, has proven to be chemical reaction, the loss of water reaction advancement as a inadequate for prediction purposes. mass should enter in the gel mass percentage of its asymptotic value. More realistic models should balance. This effect is however The evolution of the reaction extent consider the coexistence of different negligible, the degree of local water is expressed through empirical phases in a material representative saturation and the temperature level relations in terms of the previously volume element (RVE). A hierarchy being the key factors influencing the determined evolutions of tempera- of different approaches is reaction kinetics. A simple three- ture and moisture. schematically summarized in Figure phase model can there-fore be 2. The three different phases, i.e. formulated as in Figure 3, where the The two-phase model then accounts concrete skeleton, expanding gel and effect of local water pressure on the for the coupled response of the water, are assumed to be RVE static equilibrium is explicitly homogenized solid skeleton and gel homogenized within the RVE, in the accounted for. Even though this to the assigned histories of spirit of Biot’s theory of multiphase effect is negligible in quantitative temperature, of moisture and of models (see e.g. *6+). In the most terms when compared to the gel external loads. An isotropic damage general scheme, the three phases pressure, it is a necessary ingredient model has been used in *7+ for the are independently modeled, and for the simulation of concrete shrink- solid skeleton. Damage evolves as a balance equations are solved for age, where water pressure attains con-sequence of both the each of them. Different levels of negative values. microscopic gel-skeleton interaction coupling can however be considered. and the macroscopic selfstresses due The effect of concrete shrinkage is to inhomogeneous and/or restrained however negligible for the prediction concrete swelling. Anisotropic of the long term dam behavior, so damage models (see e.g. *8+) have that two phase models seem to rep- shown to better describe the local resent a viable alternative. In this skeleton-gel inter-action, at the cost case, temperature and moisture however of signifi-cantly higher distributions in space and time, complexity, so that iso-tropic models under assigned boundary conditions, seem at the moment to be more are preliminarily and independently suitable for practical applications. computed solving their diffusion equations. The gel content is then Figure 6 shows the evolution of taken to be proportional to the reaction extent and concrete Figure 3. Three phase simplified model. reaction extent, a phenomenological damage in a model of Beauharnois scalar variable varying between 0 gravity dam (Canada), built in 1932. It can be noted that damage starts developing near the dam outer skin, which is the only part affected by moisture con-tent and in the dam upper part, which is exposed to higher tempera-tures, in accordance with the usually observed real dam behavior. A macroscopic crack, due to inhomogeneous swelling, starts developing after about six years, in accordance with reported data.

Consideration of the solutions of heat and moisture transport equations in a half space under uniform and constant boundary conditions, allows to estimate the characteristic temperature and moisture lengths, i.e. the lengths where combined effects of diffusion Figure 4. Patterns of the reaction extent and damage after (a) 3, (b) 6 and (c) 60 years *7+. (of heat or moisture) and ASR occur.

15

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

Figure 5. Simultaneous influence of temperature and degree of saturation on the temperature(a) and moisture (b) characteristic lengths *7+.

Figure 5 shows the evolution of magnitude of the moisture the inspection gallery. This temperature (a) and moisture (b) characteristic length. additional crack has not been characteristic lengths as functions of documented. However, it is worth temperature and water saturation. It An application of the two phase noting that this crack, embedded appears that the temperature model, with the assumption of inside the dam body can be hardly characteristic length is two/three uniform water saturation, is shown detected and that important orders of magnitude larger than the in Figure 6. Reference is made to a retrofitting provisions, not moisture one. In concrete dams, the well-known case history in dam accounted for by the model, have structural dimensions are engineering, concerning Fontana been implemented very soon after comparable with the heat gravity dam (North Carolina, United discovering the first severe crack, characteristic length but are much States). Fontana Dam (maximum altering the pre-existing stress state. bigger than the moisture one. height, length, and thickness at the Therefore, while an important basis equal to 146, 720 and 114 m, CONCLUSIONS temperature gradient can be respectively) was completed in 1946. observed, the moisture diffusive Only 3 years later, in 1949, a pattern The evolution of the AAR in concrete process affects the external of cracking was first observed, depends on chemo-physical boundary of the structure only for a together with an upstream processes at the micro-scale. Though depth of about 40 cm and the movement of the structure. In 1972, difficult to model and to quantify, degree of saturation is constant a large longitudinal crack near the the important negative consequenc- elsewhere, equal to the imposed left abutment, in both the upstream es on concrete structures, initial value. This remark, together and downstream walls of the particularly concrete dams, that the with the outcome of Figure 4, foundation drainage gallery inside reaction may produce, justify the showing that the most relevant dam- the dam was found *10+. Figure 7 development of simplified models, age phenomena take place in the shows the obtained results in terms which could be readily employed in dam core, far from the outer skin of reaction extent and tensile engineering applications. where complete water saturation damage pattern at different time of can be conveniently assumed, the analysis. At the third year a The proposed multi-phase chemo- suggest to simplify the problem highly localized damage zone starts elastic-damage models have been formulation assuming uniform water forming at the high-left corner of the implemented as user subroutines in saturation throughout the dam body drain-age gallery and it quickly the commercial finite-element code (two-phase model *9+), thus reducing develops reaching the downstream Abaqus. Good agreement of the the problem complexity and the face and prosecuting from the other predicted damage distributions with need for parameter identification. It comer. These fully damaged bands the observed damage patterns has must be noted, however, that this well correspond to the real cracks been obtained in the applications, assumption may be inaccurate for detected by in situ inspections, see confirming the practical validity of other problems, such as, e.g. AAR in Figure 6. At a later stage, the finite element computational bridges, where the geometrical adopted model predicts localization modeling for this type of problems. dimensions are of the same order of of damage also few meters above 16

Figure 6. a) Fontana Dam: crack trajectory inside dam body (after *10+). b) Finite element mesh *9+.

BIBLIOGRAPHY modeling of concrete affected by alkali–silica reaction under variable *1+ F., Amberg, R., Bremen, N., Brizzo, temperature and humidity Rehabilitation of the Pian Telessio conditions. International Journal of dam (It) affected by AAR Reaction. Solids and Structures, 49: 3367-3380, 23rd Congress on large dams, 2012. Brasilia; Q. 90, ICOLD, 2009. *8+ C., Comi, U., Perego, Anisotropic *2+ F., Amberg, Performance of dams damage model for concrete affected affected by expanding concrete. In: by alkali–silica reaction. International Schleiss, A., Boes, R. M. (Eds.), Dams Journal of Damage Mechanics 20: and reservoirs under changing 598–617, 2011. challenges. ICOLD Meeting, CRC *9+ C., Comi, R., Fedele, U., Perego, A Press/Balkema, 115-122, 2011. chemo–thermo-damage model for *3+ C., Larive, Apports combinés de the analysis of concrete dams l’expérimentation et de la affected by alkali–silica reaction. modélisation la comprehénsion de Mechanics of Materials 41: 210–230, l’alcaliréaction et de ses effets 2009. mécaniques. Ph.D. Thesis, LCPC, *10+ A.R., Ingraffea, Case studies of 1998. simulation of fracture in concrete *4+ S., Multon, Evaluation expéri- dams. Engineering Fracture mentale et théorique des effects Mechanics, 35: 553–564, 1990. mechaniques de l’alcaliréaction sur des structures modéles. Ph.D. Thesis, LCPC, 2003. CLAUDIA COMI *5+ S., Multon, F., Toutlemonde, UMBERTO PEREGO Effect of applied stresses on alkali- silica reaction-induced expansions. DEPARTMENT OF STRUCTURAL Cement and Concrete Research, 6: ENGINEERING 912–920, 2006. POLITECNICO DI MILANO *6+ F.J., Ulm, O., Coussy, L., Kefei, C., ITALY Larive, Thermo-chemo-mechanics of asr expansion in concrete structures. IN COLLABORATION WITH Figure 7. Fontana Dam: reaction extent (first ASCE Journal of Engineering column) and tensile damage (second R. FEDELE Mechanics 126: 233–242, 2000. column) after 3, 7 and 25 years of AAR *7+ C., Comi, B., Kirchmayr, R., R. PIGNATELLI activity *9+. Pignatelli, Two-phase damage [email protected]

17

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

HIGH ORDER RESIDUAL DISTRIBUTION SCHEMES FOR COMPRESSIBLE FLOWS ON UNSTRUCTURED GRIDS

INTRODUCTION grids, higher-order methods (order Because of these difficulties, higher- of accuracy more than two) are used. order schemes are still not popular In the last decades Computational outside academic world, but in the Fluid Dynamics (CFD) has seen an The advantage of high-order recent years their capabilities has impressive development, and is methods is well known. Assuming been pushed further and further. currently used to perform very that the solution is approximated complex simulations in industrial with a polynomial of order k, than an Between the higher-order schemes, applications. For example, in the effective numerical method will probably, one of the most used aerodynamic field, CFD tools are provide a numerical solution with an method is the Discontinuous routinely used in the design and the error e which decrease with a Galerkin (DG) method *1+, which is optimization of aircraft as characteristic mesh size parameter h based on a weak formulation of the complement to the wind tunnel according to e ∞ hk+1, this means governing equations with a discon- experiments. that, since the error decreases with a tinuous polynomial approximation of power k + 1, for a given number of the solution on each element. While Most of the numerical schemes used DOF higher-order methods computationally compact and in industrial applications are based guarantee a smaller level of error probably more flexible, DG schemes on at most second-order accurate respect to low-order scheme. On the suffer from the serious drawback of a Finite Volume (FV) methods other hand, increasing the order of very fast growth of the number of developed in the last twenty years. discretization produces an increase DOF with the cell polynomial degree Despite the great capabilities shown, also of the cost of the numerical and also the limiting strategy is not these schemes still suffer from an scheme, however since the cost completely clear in the contest of the excess of dissipation which make the grows with a fixed factor respect to DG schemes. An interesting simulation of many flow problems the order it appears that the alternative to the DG schemes is very difficult. For example, computational effort required by represented by the Residual aerodynamic flows at high Reynolds higher-order schemes to get a Distribution (RD) method *3+. In RD number with standard FV codes solution with a fixed level of error is schemes the formulation remains would require a such great number smaller respect to that of low-order local, as in DG, but the number of of degrees of freedom (DOF) to make schemes. DOF growths less quickly. The price the simulations of complex aircraft to pay is to impose a continuous configurations very expensive and Despite their advantages, high-order approximation of the solution, alt- time consuming, in particular for methods are not mature yet and, hough a discontinuous formulation is unsteady applications. generally, not as robust as low-order possible in the RD framework *4, 5+. ones. Convergence to the steady Although the progress made in the state maybe more problematic, The history of RD schemes goes back construction of computers has made especially for stiff problems like to the pioneering work of Ni, an possible to perform simulations with turbulent flows. Furthermore, for engineer at Bombardier in the an extremely large number of transonic or supersonic flows with 80’ *10+. He has proposed a degrees of freedom (DOF), in the last discontinuities, higher-order reformulation of the Lax-Wendrof years the problem to perform methods lack the robustness of schemes that can be coded in a way accurate simulations of complex monotone shock capturing very similar to what is done in Finite flows has been tackled down with a techniques developed in the past Elements: a loop over elements, and new strategy: instead of using finer decades for low-order schemes. a series of local operations within

18

each element; all this in RD METHODS The choice of the distribution combination with multigrids. P.L. process is a crucial aspect in the RD Roe’s contributions have been The RD method refers to a class of methods. Some researchers instrumental. In particular his numerical methods, in which the developed upwind schemes which celebrated Roe solver *11+ is indeed numerical solution is used to are particularly suitable for designed as a Fluctuation Splitting construct a cell residual which is advective flows *8+. However such scheme, another name for RD then distributed to all DOF of each schemes have some drawbacks: schemes. He has also, in collabora- element. they are generally formulated for tion with D. Sidilkover, designed and simplectic elements (triangles or analyzed the N scheme, a genuinely The solution on each element is tetrahedral) and their extension to multidimensional extension of the approximated by functions of order high-order formulation is Murman-Cole-Roe scalar scheme, as k, which are assumed continuous on complicated *9+. We focus instead well as many other contributions in each element and on the face of on the construction of central the field. The work has been further two adjacent elements. Standard schemes which can be formulated continued by the VKI group, under Lagrange polynomial are used to on every kind of element and can be H. Deconinck’s leadership, in parallel construct such a continuous easily extended to any order. with and in collaboration with the U. approximation (but also other Another important aspect in the of Michigan at Ann Arbor group approach can be used, e.g., NURBS distribution process is the possibility leaded by P.L Roe. This work has *12+). In the case of steady state to construct a non-linear scheme been given in many VKI lecture simulations, the approximated through a so-called limiting series in the 90’s and early 2000’s. solution is used to define an integral procedure which allows one to The idea was to be able to combine quantity on each element, called construct a high order scheme with ideas coming from finite element total residual, which is given by the monotone shock capturing *2+. methods (the structure of the integral of the governing equations scheme), up-winding and evaluated on the element. In order In order to accelerate the monotonicity preservation as for to obtain an equation for each DOF, convergence to the steady state, finite volume so-called high order the total residual is first locally implicit schemes must be used. We scheme, via non linear limiting. distributed to each DOF of the constructed an effective implicit Other contributors have been H. element and the contributions scheme based on Newton-Krylov Nishikawa (now at NIA), D. Caraeni coming from different elements are method with a preconditioned (then in Lundt University, now in a collected for each DOF. This defines matrix-free technique *7, 6+. private company in the US), a group a set of non-linear equations which in Lisbon, a group in Bari (Italy) is solved via iterative methods. The RESULTS under the supervision of P. de Palma change of the values of the solution and Pr. Napolitano, and the during the iterative process is driven We present in the follows some Bordeaux group. The present paper by the non-zero total residuals on results obtained with our RD is an emanation of this latter group. the elements, at the steady state method in order to show the The RD schemes have very natural the total residual on each element effectiveness of the numerical connections with the streamline vanishes and the steady state scheme. The first test considered is diffusion method, as well as the solution is obtained. laminar flow over a NACA-0012 Residual Based Compact scheme by airfoil at zero angle of incidence, the Lerat and collaborators at ENSAM, Paris, France.

Here, we present some of the recent achievements in the RD schemes, showing that such kind of schemes can be successfully applied to perform the simulation of complex flows, like three dimensional turbulent flows.

Figure 1: Mach number contours for the NACA airfoil. Second order (left) and third order (right) scheme on two different grids with approximately the same number of DOF.

19

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

Figure 2: Particular of the stream-trace near the trailing edge. Second order (left) and third order (right) scheme on two different grids with approximately the same number of DOF.

Figure 3: Grid used numerical simulation of the flow around the delta wing

Mach number is take as 0.5 and the the stream-traces near the trailing solution respect to the second order Reynold number considered is 5 edge for the two orders of scheme. 000. One distinctive feature of this approximation are shown. flow is a steady laminar recirculation As last test case we consider a bubble near the trailing edge of the A more complex test case consider turbulent flow over the M6 wing, airfoil. The Mach contours of the here is the laminar flow over a delta the Spalart-Allmaras turbulence ◦ numerical simulation are shown in wing at 12.5 angle of incidence, the model has been used. The angle of ◦ Fig. 1 for the second and third order Mach number is 0.3 and the incidence is 3.06 , the Mach number scheme. The second order Reynolds number is 4 000 based on is 0.8395 and the Reynold number is 6 simulation has been performed on a the cord at the root. The flow is 11.72×10 bases on the mean grid with approximately 30 000, characterized by a main and aerodynamic chord. The pressure nodes while the third order secondary vortex on the wing which coefficient of the third order simulation has been performed on a remain over long distance far simulation is shown in Fig. 5, while coarser grid such that the number of behind the wing. The computational in Fig. 6 is reported the comparison DOF for the two simulation is almost grid is depicted in Fig. 3 while in Fig. between the numerical and the the same. It can be seen that the for 4 are shown the Mach number and experimental results at two a given number of DOF the higher- the vorticity contours at different different stations along the wing. order scheme produce a better stations along and behind the wing. approximation of the solution, this It can be seen how the higher order can be observed also in Fig. 2 where scheme produce a more accurate

Figure 4: Mach number (left) and vorticity magnitude (right) contours at different section along and behind the delta wing. Second and third order solutions are compared.

20

[5] M.E. Hubbard, A framework for discontinuous fluctuation distribution. International Journal for Numerical Methods in Fluids 2008; 56:1305-1311. [6] R. Abgrall and D. De Santis and M. Ricchiuto, High order residual distribution scheme for RANS equations. Seventh International Conference on Computational Fluid Dynamic, ICCFD7-2802. [7] R. Abgrall and D. De Santis and Figure 5: Pressure coefficient contours aver the M6 wing for the third order turbulent M. Ricchiuto, High order residual simulation. distribution scheme for Navier- Stokes equations. 20th AIAA Computational Fluid Dynamics Conference, AIAA 2011-3231. [8] Struijs R., Deconinck H., Roe P.L., Fluctuation splitting schemes for the 2D Euler equations. Computational Fluid Dynamics. 1991: VKI-LS 1991- 01 [9] N. Villedieu, T. Quintino, M. Ricchiuto, H. Deconinck, Third order residual distribution schemes for the Navier–Stokes equations. Journal of 2011:230 Figure 6: Comparison between the numerical and the experimental values of pressure 4301-4315. coefficient at different station along the M6 wing for the third order turbulent simulation. [10] R.-H. Ni. A multiple grid scheme for solving the Euler equations. AIAA J., 20:1565–1571, 1981. CONCLUSIONS REFERENCES [11] P.L. Roe. Approximate Riemann solvers, parameter vectors, and We presented a high order accurate [1] B. Cockburn, G.E. Karniadakis, difference schemes Journal of and robust Residual Distribution C.W. Shu, Discontinuous Galerkin Computational Physics, 1981:43 357 scheme for the solution of Navier- methods: theory, computation and -372. Stokes equations. The use of a application (Lecture notes in [12] A. Froehly, Couplage d’un continuous approximation of the computational science and schéma aux résidus distribués `a solution guarantees a smaller engineering). (eds), Springer, Berlin, l’analyse isogéométrique : méthode number of DOF respect to 2000. numérique et outils de génération et discontinuous methods like DG. At [2] R. Abgrall, Essentially non- adaptation de maillage Université de the same time the RD method proved oscillatory residual distribution Bordeaux, September 7th, 2012. to be effective in the simulation of schemes for hyperbolic problems. complex test cases, including Journal of Computational Physics turbulent flows. 2006; 214:773:808. REMI ABGRALL [3] R. Abgrall R., P.L. Roe, High-order DANTE DE SANTIS ACKNOWLEDGEMENTS fluctuation schemes on triangular meshes. Journal of Scientific MARIO RICCHIUTO R. Abgrall has been supported by the Computing 2003; 19:3-36. INRIA ET UNIVERSITE DE BORDEAUX European Research Council Advanced [4] R. Abgrall, C.W. Shu, Development FRANCE Grant ADDECCO # 226316. D. de of residual distribution schemes for Santis has been supported by the the discontinuous : [email protected] IDIHOM EC FP7 grant. the scalar case with linear elements. Communications in Computational Physics 2009; 5:376- 390.

21

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

ECCOMAS THEMATIC AND REGIONAL CONFERENCES 2013

8. 4TH INTERNATIONAL CONFERENCE ON 15. XII INTERNATIONAL CONFERENCE ON 1. CONFERENCE ON COMPUTER COMPUTATIONAL METHODS IN COMPUTATIONAL PLASTICITY METHODS IN MATERIALS STRUCTURAL DYNAMICS AND (COMPLAS XII) TECHNOLOGY - KOMPLASTECH 2013 EARTHQUAKE ENGINEERING September 3 - 5 in Barcelona, Spain January 13 - 16 in Zakopane, (COMPDYN 2013) ECCOMAS THEMATIC CONFERENCE ECCOMAS THEMATIC CONFERENCE June 12 - 14 in Kos Island, Greece ECCOMAS THEMATIC CONFERENCE 16. SECOND CONFERENCE ON THE 2. COMPUTATIONAL METHODS IN EXTENDED FINITE ELEMENT METHOD TUNNELLING - EURO:TUN 2013 9. III SOUTH-EAST EUROPEAN - XFEM 2013 April 17 - 19 in Ruhr-University CONFERENCE ON COMPUTATIONAL September 11 - 13 in Lyon, France Bochum, Germany MECHANICS ECCOMAS THEMATIC CONFERENCE ECCOMAS THEMATIC CONFERENCE June 12 - 14 in Kos Island, Greece ECCOMAS SPECIAL INTEREST 17. SIMULATION AND MODELING OF 3. 3RD INTERNATIONAL CONFERENCE CONFERENCE BIOLOGICAL FLOWS - SIMBIO 2013 ON INVERSE PROBLEMS IN September 12 - 14 in Lisbon, Portugal MECHANICS (IPM 2013) 10. V INTERNATIONAL CONFERENCE ON ECCOMAS THEMATIC CONFERENCE April 24 - 27 in Baranow, Poland COUPLED PROBLEMS IN SCIENCE AND ECCOMAS THEMATIC CONFERENCE ENGINEERING (COUPLED 2013). A 18. III INTERNATIONAL CONFERENCE ON CONFERENCE CELEBRATING THE PARTICLE-BASED METHODS 4. V INTERNATIONAL CONFERENCE ON 60TH. BIRTHDAY OF EUGENIO ONATE (PARTICLES 2013) COMPUTATIONAL METHODS IN June 17 - 19 in Ibiza, Spain September 18 - 20 in Stuttgart, MARINE ENGINEERING (MARINE ECCOMAS THEMATIC CONFERENCE Germany 2013) ECCOMAS THEMATIC CONFERENCE May 29 - 31 in Hamburg, Germany 11. CONFERENCE ON SMART ECCOMAS THEMATIC CONFERENCE STRUCTURES AND MATERIALS 19. IV ECCOMAS THEMATIC CONFERENCE (SMART'13) ON MECHANICAL RESPONSE OF 5. VI INTERNATIONAL CONFERENCE ON June 24 - 28 in Torino, Italy COMPOSITES ADAPTIVE MODELING AND ECCOMAS THEMATIC CONFERENCE September 25 - 27 in Azores, Portugal SIMULATION (ADMOS 2013) ECCOMAS THEMATIC CONFERENCE June 3 - 5 in Lisbon, Portugal 12. MULTIBODY DYNAMICS 2013 ECCOMAS THEMATIC CONFERENCE July 1 - 4 in Zagreb, Croatia 20. VI INTERNACIONAL CONFERENCE ON ECCOMAS THEMATIC CONFERENCE TEXTILE COMPOSITES AND 6. COMPUTATIONAL MODELING OF INFLATABLE STRUCTURES FRACTURE AND FAILURE OF 13. III INTERNATIONAL CONFERENCE ON (STRUCTURAL MEMBRANES 2013) MATERIALS AND STRUCTURES (CFRAC COMPUTATIONAL CONTACT October 9 -11 in Munich, Germany -2013) MECHANICS (ICCCM13) ECCOMAS THEMATIC CONFERENCE June 5 - 7 in Prague, Czech Republic July 10 - 12 in Lecce, Italy ECCOMAS THEMATIC CONFERENCE ECCOMAS THEMATIC CONFERENCE 21. VIPIMAGE 2013 - IV ECCOMAS 7. III INTERNATIONAL CONFERENCE ON THEMATIC CONFERENCE ON 14. 20TH INTERNATIONAL CONFERENCE TISSUE ENGINEERING - ICTE2013 COMPUTATIONAL VISION AND ON COMPUTER METHODS IN June 6 - 8 in Leiria, Portugal MEDICAL IMAGE PROCESSING MECHANICS ECCOMAS THEMATIC CONFERENCE October 17 - 19 in Melia Madeira August 27 - 31 in Poznan, Poland Mare, Madeira Island, Portugal ECCOMAS REGIONAL CONFERENCE ECCOMAS THEMATIC CONFERENCE

22

ND 2 ECCOMAS YOUNG INVESTIGATORS CONFERENCE 2013 ECCOMAS PhD OLYMPIAD 2013

The 2nd ECCOMAS Young Investigators Conference (YIC 2013) will be held in September 2-6, 2013 in Bordeaux (France). It is organized by CSMA and SMAI affiliated young French scientists. This new scientific event focuses on applications of mathematical and computational methods and modeling to different areas of simulation in engineering. It mainly aims at bringing together junior researchers and PhD students in the field to encourage discussions and collaborations. It is also a place where senior scientists are invited to share their experience. New creative and innovative formats will be introduced for this conference, such as opening presentations (short courses) by junior scientists for each session, plenary talks with opponents, lightning talk and "software" sessions, and more... About 150 attendees are expected for the conference, and abstract submission is open via the conference website: http://yic2013.sciencesconf.org L. CHAMOIN C. DOBRZYNSKI CHAIRPERSONS OF YIC 2013

[email protected]

TOPICS Among others, those thematics will be particularly high- lighted:  Computational Solids and Structural Mechanics  Computational Fluid Dynamics  Computational Applied Mathematics  Scientific Computing  Industrial Applications and Challenges

IMPORTANT DATES  Abstract submission deadline: March 15, 2013  Notification of abstract acceptance: May 15, 2013  Early registration deadline: May 31, 2013

23

EUROPEAN COMMUNITY ON COMPUTATIONAL METHODS IN APPLIED SCIENCES

MEMBER ASSOCIATIONS FORUM CIRM-FBK, TRENTO-ITALY CALL FOR APPLICATIONS FOR THE 2013 ACTIVITIES: CONFERENCES, POSTDOC POSITION, VISITING POSITIONS, RE- SEARCH IN PAIRS

CIRM (Centro Internazionale per la Ricerca Matematica, see http://cirm.fbk.eu/en/home ) at Trento-Italy, upon pro- posal of its ADVISORY BOARD, is going to continue and expand its activities in the field of Mathematics Research, in collaboration with the Department of Mathematics of the University of Trento, according to the following program.

I) Conferences. This is the traditional activity, which has seen the CIRM organizing more than 250 mathematical Meetings and Conferences in the period 1979-2012. Applications for the year 2013 should be sent before September 15, 2012, by ordinary mail to the address: Fonda- zione Bruno Kessler, Centro Internazionale per la Ricerca Matematica, Via Sommarive n. 14 – Povo, 38123 Trento or via electronic mail to the address [email protected] , following the usual scheme. In particular they must contain: a scientific proposal, with names of possible speakers, a rather detailed financial budget, which CIRM will be asked to finance partially (up to 50%), specification of the other available or planned financial resources.

II) Postdoc Fellowships. The support of younger researchers is one of the vital functions of a research centre; in this direction CIRM opens the application for an annual Postdoc Position, for researches in the field of Mathematics in the year 2012-2013. The fellowship amounts to Euro 23500 gross per year (taxes should be deducted), which will be supplied in 12 monthly instalments, and requires the presence at the CIRM Institute (address Via Sommarive 14, Povo, Trento). The applications must arrive not later than September 15, 2012, by ordinary or electronic mail to the above men- tioned addresses. For further information please visit the website http://cirm.fbk.eu/en/node/134 .

CIRM activities include two more programs, a program of CIRM Visiting Professors and CIRM Visiting Scholars and another one of Research in Pairs, the latter similar to the one of the Mathematisches Forschungsinstitut Oberwol- fach, which expects the presence of two or three partners from Universities located in different towns, who intend to work together at CIRM, on a definite research project for a period of time ranging from one to six weeks. Applications for these two programs must be sent to the address mentioned above by ordinary or electronic mail. They can be submitted at any time and they must contain a specific indication of the proposed dates for the visit; it is suggested to send applications at least three months before the planned stay since decision making might require a long time. Applications must contain a brief research proposal, a C.V. and a list of publications for each applicant. For further information please visit the websites: http://cirm.fbk.eu/en/Visiting+Professors and http://cirm.fbk.eu/ en/node/122 .

The Director of CIRM Prof. Marco Andreatta

24

ECCOMAS CONGRESS 2016

CRETE, GREECE JUNE 2016

Europe Greece

Crete Heraklion

25