From Plastic Pipes and Bottles to Bioengineering Applications: Fluid-Structure Interaction Procedures for Flexible Systems

Transcription

1 University College Dublin UNIVERZITET U ZENICI Mašinski fakultet From Plastic Pipes and Bottles to Bioengineering Applications: Fluid-Structure Interaction Procedures for Flexible Systems Aleksandar Karac Zagreb, June 2007

2 FV method Coupling schemes in OpenFOAM Two-system approach: Fluid Solid Analysis method Strongly coupled interacting Solids relatively flexible: plastic pipelines, containers, tanks/reservoir sloshing, airbags, cardiovascular flows One-system approach: Fluid + Solid Analysis method Fully coupled coupled Solid-fluid phase transformations: material processing/forming (moulding, casting, extrusion, welding, ) 2/29

3 FV method C. Greenshields H. Jasak H. Weller Two-system approach: START OF TIME STEP Fluid Solid Analysis method 1. Solve fluid (or 3.) 2. Pass infromation from the fluid to solid (or 4.) 3. Solve structure (or 1.) 4. Pass infromation from the solid to fluid (or 2.) END OF TIME STEP (implicit scheme only) 3/29

4 for ( FV method ) { runtime++, runtimes++;!runtime.end() &&!runtimes.end(); runtime++, runtimes++ # include "movesolidmesh.h" # include "movefluidmesh.h" # include "updateinletdata.h" do { # include "fluid.h" # include "fluidtosolid.h" # include "structure.h" # include "solidtofluid.h" # include "updateproperties.h" Fluid Solid Analysis method } while (initialresidualsolid > convergencetolerancesolid && ++icorrs < ncorrs); # include "calculatefields.h" } } return(0); 4/29

5 FV method C. Greenshields H. Jasak H. Weller Two-system approach: START OF TIME STEP Fluid Solid Analysis method 1. Solve fluid (or 3.) 2. Pass infromation from the fluid to solid (or 4.) 3. Solve structure (or 1.) 4. Pass infromation from the solid to fluid (or 2.) - incompressible flow - compressible flow - multiphase flow - turbulance flows - fluid models - mesh motion possibility -... END OF TIME STEP (implicit scheme only) 5/29

6 FV method C. Greenshields H. Jasak H. Weller Two-system approach: START OF TIME STEP Fluid Solid Analysis method 1. Solve fluid (or 3.) 2. Pass infromation from the fluid to solid (or 4.) - pass tractions (pressure + shear stresses) - one-to-one exchange - exchange via interpolation 3. Solve structure (or 1.) 4. Pass infromation from the solid to fluid (or 2.) END OF TIME STEP (implicit scheme only) 6/29

7 FV method C. Greenshields H. Jasak H. Weller Two-system approach: START OF TIME STEP Fluid Solid Analysis method 1. Solve fluid (or 3.) 2. Pass infromation from the fluid to solid (or 4.) 3. Solve structure (or 1.) 4. Pass infromation from the solid to fluid (or 2.) - small-strains (total Lagrangean, incremental Lagrangean) - large-strains (updated Lagrangean; mesh motion) - material models (elastic, non-linear elastic, elasto-plastic,...) - cohesive boundaries to account for fracture END OF TIME STEP (implicit scheme only) 7/29

8 FV method C. Greenshields H. Jasak H. Weller Two-system approach: START OF TIME STEP Fluid Solid Analysis method 1. Solve fluid (or 3.) 2. Pass infromation from the fluid to solid (or 4.) 3. Solve structure (or 1.) - pass displacement (increments), velocity field - one-to-one exchange - exchange via interpolation 4. Pass infromation from the solid to fluid (or 2.) END OF TIME STEP (implicit scheme only) 8/29

9 FV method C. Greenshields H. Jasak H. Weller Two-system approach: START OF TIME STEP Fluid Solid Analysis method 1. Solve fluid (or 3.) 2. Pass infromation from the fluid to solid (or 4.) 3. Solve structure (or 1.) STABILITY ISSUES 4. Pass infromation from the solid to fluid (or 2.) END OF TIME STEP (implicit scheme only) 9/29

10 FV method C. Greenshields H. Weller One-system approach: Fluid + Solid Analysis method Fluid: t 2 ρvdv+ ρvv nds= 2ηε ηtr( ε) I pi ds 3 n V S S Solid: 2 ρ dv ρ ds 2Nε Ntr( ε) p ds ds t v + vv n = I I 3 n + Σ n V S S S t t 1 2 N= µ t, p= trσ = Kstr ε, 2 t ( ) Σ= µε µ r ε dt 3 I /29

11 Falling bottles and containers Problem description H p =p 0 v v =0? v 1 a) b) c) 11/29

12 Falling bottles and containers Eksperimetnal set-up Upper cap Water-filled bottle Strain gauges Tie bars PT holder Rubber O-ring Pressure transducers PT3 PT1 PT2 SG1 SG2 Jubilee clip Lower cap Rubber cushion SG1 SG2 PT3 12/29

13 Falling bottles and containers Numerical set-up Simple geometry Flat base Free surface (p=const=1bar) Solid Fluid 50 cells Impact end ( V=0, U=0) Curved base 20 cells 3 cells 13/29

14 Falling bottles and containers OpenFOAM simulation simple geometry 14/29

15 Falling bottles and containers OpenFOAM simulation base shape effect 15/29

16 Falling bottles and containers Some quantitative results Pressure, bar experiment FV, two-system FV, one-system Time, ms 16/29

17 Falling bottles and containers Some quantitative results 6 10 Pressure, bar Pressure, bar experiment FV, two-system FV, one-system Time, ms experiment FV, two-system FV, one-system Time, ms 17/29

18 Falling bottles and containers Simulations with fracture numerical s Cohesive zone model Symmetry plane Free surface (fixed pressure and zero velocity gradient) Interface solid Symmetry plane Outside traction free Interface fluid Crack interface solid Base direction mixed Crack interface fluid 18/29

19 Falling bottles and containers Simulations with fracture OpenFOAM simulation 19/29

20 Fracture of PE pipes EPRSC C. Greenshields H. Jasak P. S. Leevers K. C. Pandya V. Tropsa G. Venizelos J. G. Williams Problem description - Small scale steady-state (S4) Video: A. Karac, A. Paizis 20/29

21 Fracture of PE pipes H. Jasak V. Tropsa OpenFOAM simulation 21/29

22 Fracture of PE pipes H. Jasak V. Tropsa OpenFOAM simulation 22/29

23 Bioengineering blood flow related problems Atherosclerosis Abdominal Aortic Aneurysm Fluid-Structure Interaction in Bioengineering by Valentine Kanyanta 23/29

24 Bioengineering - blast trauma to lungs MoD O. Alakija Transport Sport & horse Racing Military/police Experiments on live animals 24/29

25 Bioengineering - blast trauma to lungs MoD O. Alakija Combined numerical/experimental study σ Stress Strain ε Alveoli (70-350µm) σ Bulk Alveolar wall (5µm thick) Damage δ Scale bridging 25/29

26 Bioengineering - blast trauma to lungs MoD O. Alakija v alveoli wall (tissue) air space symmetry on both sides and fixed at the bottom (a) (b) (a) Micrograph of a cross-section through an alveoli duct and sac system (b) Two dimensional mesh of the alveoli stack Wave propagation v=20 m/s 26/29

27 ? Bioengineering - blast trauma to lungs H. Parsa Construction of the surrogate lung model Polymeric foam filled with microcapsules Optical microscope view of a microcapsule 27/29

28 Bioengineering - blast trauma to lungs H. Parsa CT Scan of the surrogate lung * Before St. Vincent s University Hospital CT scanner After * Walsh, C., Final year project, UCD, /29

29 Bioengineering other ongoing projects A. Safari detachment of bacterial biofilms fracture properties of adhesives fracture properties of diamonds OF SCIENCE, TECHNOLOGY for cutting AND MEDICINE tools simulation of 1 st metatarsalphalangeal joint simulation of hip joint 29/29

30 THANK YOU ALEKSANDAR KARAČ Mašinski fakultet Fakultetska ZENICA Bosnia and Herzegovina TEL: +387 (0) ext. 140 FAX: +387 (0) s: akarac@mf.unze.ba aleksandar.karac@ucd.ie A. KARAC, From Plastic Pipes and Bottles to Bioengineering Applications. 3 rd Progress