COMSOL Conference 2008

Similar documents
In this section is given an overview of the common elasticity models.

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

Airflow and Contaminant Simulation with CONTAM

Introduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012

Effects of Polymer Concentration and Molecular Weight on the Dynamics of Visco-Elasto- Capillary Breakup

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

NUMERICAL SIMULATION OF FLOW AROUND A NACA 0012 AIRFOIL IN TRANSIENT PITCHING MOTION USING IMMERSED BOUNDARY METHOD WITH VIRTUAL PHYSICAL MODEL

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Turbulent Flow. Turbulent Flow

A Frequency-Domain Approach for Transient Dynamic Analysis using Scaled Boundary Finite Element Method (I): Approach and Validation

Study Guide For Exam Two

Lecture Note 3. Eshelby s Inclusion II

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

COMPUTATIONAL FLUID DYNAMICS ME G515

One Dimensional Axial Deformations

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Introduction to Turbulence Modeling

Tensor Smooth Length for SPH Modelling of High Speed Impact

Simulation of Turbulent Flow Using FEM

Numerical Simulation of Lid-Driven Cavity Flow Using the Lattice Boltzmann Method

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Chapter 02: Numerical methods for microfluidics. Xiangyu Hu Technical University of Munich

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Physics 181. Particle Systems

Macroscopic Momentum Balances

Survey of applications of discrete vortex method in civil engineering

A Cartesian-grid integrated-rbf method for viscoelastic flows

NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD

CONTROLLED FLOW SIMULATION USING SPH METHOD

MEEM 3700 Mechanical Vibrations

Second Order Analysis

11. Dynamics in Rotating Frames of Reference

Optimal Control of Temperature in Fluid Flow

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Flow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved

REAL TIME AIRFLOW SIMULATION IN BUILDINGS

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

Chapter 3. r r. Position, Velocity, and Acceleration Revisited

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

NUMERICAL SIMULATION OF FLOW OVER STEPPED SPILLWAYS

Numerical Heat and Mass Transfer

Indeterminate pin-jointed frames (trusses)

Computational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH

THE NEAR-WALL INFLUENCE ON THE FLOW AROUND A SINGLE SQUARE CYLINDER.

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P.

Please initial the statement below to show that you have read it

Simulation of Heat Transfer during Artificial Ground Freezing Combined with Groundwater Flow

CFD simulation with multiphase flows in porous media and open mineral storage pile

Lecture # 02: Pressure measurements and Measurement Uncertainties

A STRONG COUPLING SCHEME FOR FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMICALLY MOVING BOUNDARIES IN VISCOUS INCOMPRESSIBLE FLOWS

Solution of the Navier-Stokes Equations

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Simulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM

A new integrated-rbf-based domain-embedding scheme for solving fluid-flow problems

Introduction to Compact Dynamical Modeling. I.2.a Conservations and Constitute Laws. Luca Daniel Massachusetts Institute of Technology

NUMERICAL APPLICATIONS AND VERIFICATION OF AN INTEGRATED FLOW-STRESS MODEL IN PROCESSING OF THERMOSET COMPOSITES

Modeling of Dynamic Systems

1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #3: Hydraulic Head and Fluid Potential. p o. p o + p

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Structural Dynamics and Earthquake Engineering

LAB 4: Modulus of elasticity

Finite Element Modelling of truss/cable structures

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

GENERAL THEORY FOR THE MODELLING OF LAMINATED MULTI-LAYER PLATES AND ITS APPLICATION FOR THE ANALYSIS OF SANDWICH PANELS WITH PLY DROPS

Lecture 8 Modal Analysis

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

The Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially

Higher Order Wall Boundary Conditions for Incompressible Flow Simulations

Turbulent Flow in Curved Square Duct: Prediction of Fluid flow and Heat transfer Characteristics

16 Circular Footing on a Semi-infinite Elastic Medium

The Finite Element Method

One-sided finite-difference approximations suitable for use with Richardson extrapolation

DEFORMABLE GROUND CONTACT MODELING: AN OPTIMIZATION APPROACH. Jeff Reinbolt. EML 5595 Mechanics of the Human Locomotor System Final Project Report

TURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March

Introduction to Computational Fluid Dynamics

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Axial Turbine Analysis

Energy configuration optimization of submerged propeller in oxidation ditch based on CFD

THREE-DIMENSION DYNAMIC SOIL-STRUCTURE INTERACTION ANALYSIS USING THE SUBSTRUCTURE METHOD IN THE TIME DOMAIN

Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Mechanics Physics 151

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

Thermodynamics General

HYBRID LBM-FVM AND LBM-MCM METHODS FOR FLUID FLOW AND HEAT TRANSFER SIMULATION

PHYS 705: Classical Mechanics. Newtonian Mechanics

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Physics 207 Lecture 6

Numerical Nonlinear Analysis with the Boundary Element Method

The Dissolution and Transport of Radionuclides from Used Nuclear Fuel in an Underground Repository

Fictitious Boundary and Moving Mesh Methods for the Numerical Simulation of Particulate Flow

Equation of State Modeling of Phase Equilibrium in the Low-Density Polyethylene Process

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems:

Transcription:

Presented at the COMSOL Conference 008 Boston Sheer, F.J. 1, Ghadal, S.N. 1 Mechancal Engneerng Department, Oho State Unversty, Columbus, OH, sheer.4@osu.edu Bomedcal Engneerng Department, Oho State Unversty, Columbus, OH, ghadal.1@osu.edu COMSOL Conference 008

Background and Motvaton Eustachan Tube (ET) connects the Mddle Ear to the Nasopharyn. ET has three prmary functons: Dranage of ME flud Equlbrate ME pressure Protect ME from Pathogens ET dysfuncton has been drectly lnked to Otts Meda. Otts Meda has an annual healthrelated cost of $4 bllon dollars n the US. [1] Reprnted from Bluestone, Charles D, and Klen, Jerome O., Otts Meda n Infants and Chldren Thrd edton, W.B. Saunders Company, 001

Modelng Goals In the clncs, we can perform a dagnostc procedure called a Forced Response Test (FRT) Patent s hooked up to epermental apparatus shown above and asked to swallow Frst Goal: Create a transent COMSOL model that can replcate FRT results to serve as model valdaton Use valdated model to eplore other parameters nvolved n ET functon

Epermental Data Normal ET Functon: R v (mmh O/cc/mn) 40 30 0 10 0 Resstance P/Q Epermental Rv 1 3 4 5 6 Tme (seconds) Patents wth ET dysfuncton: R v (mmh O/cc/mn) 30 Epermental (Inflamed) 5 0 15 10 5 0.5 3.0 3.5 4.0 4.5 5.0 5.5 Tme (seconds)

Structure of the Eustachan Tube Defne Cross Sectons Loft Cross Sectons Full 3D CAD Model y (Superor) Nasopharyn z (Dstal) (Lateral) TVPM Cartlage Lumen Mddle Ear LVPM Glandular Tssue Hamulus

Fnal Analyzed 3D Geometry Green doman s the Cartlage Red doman s the Fatty Tssue Blue doman s the Lumen Openng

Sold Mechancs Propertes Both Cartlage and Glandular Tssue are represented by Mooney Rvln Hyperelastc materal models. E cart 300kPa [,3] E Gland 50kPa [,3] κ E 3(1 υ ) C 1 E 4(1 +υ) 1 C J X W hyp P where u u Y Z W hyp ( I 3) + κ( ) 01 1

Sold Boundary Condtons Boundary Condtons: Attachment of ET to the cranal base of the skull Attachment of soft tssue to the bony portons of the ET at the promal and dstal ends Appled Loads: Tensor Vel Palatn (TVP) muscle force Levator Vel Palatn (LVP) muscle force Flud forces

Sold Boundary Condtons (cont) All TVPM load vectors pont to the poston of the Pterygod Hamulus All LVPM load vectors are normal to the appled surface

Sold Flud Couplng Arbtrary Lagrange Euleran (ALE) equatons 0 j X Z Y X X j z y where & In the Sold Doman, flud force s appled on the nterface boundares In Flud Doman, veloctes from sold deformaton s appled to nterface boundares T n p n T )) ) ( ( ( u u I + + η σ y T y y n p n T )) ) ( ( ( u u I + + η σ z T z z n p n T )) ) ( ( ( u u I + + η σ t u u t v v t w w where u, v, and w are the sold dsplacements and u, v, and w are the flud veloctes

Flud Doman Due to the complety of the geometry, a coarse, tetrahedral mesh s used wth hgh order elements to accommodate the large dsplacements wthout any element collapse. Incompressble Contnuty and Naver Stokes equatons 0 v ρ j j v p v v t v + + μ ρ ρ w v u v where Flud equatons are solved n the movng reference frame!!

Model Results Sold Dsplacement feld and Flud Velocty vectors 510-4 m YZ ε eff ( ε X + εy + ε Z ) + ( ε XY + ε XZ + ε 3 ) Effectve Stran 0.35 Model Troubles: Statc solutons only Model unstable Fcttously vscous flud doman Sold Flud couplng only 0 m 0

Test Model Anatomcally correct lumen openng Lumen openng surrounded by cylnder of cartlage Cylnder s loaded wth a unform, outward radal force to epand lumen openng.

Test Model Results

Test Model Results

Future Work Incorporate a vscoelastc materal model for the Fatty Tssue E 1 E η Smulate the dranage/flow of hghly vscous mucus through the ET U Through the use of Steered Molecular Dynamcs, create a multscale model whch ncludes molecular adhesons Flud-mucosa layer ME or bronch γ μ F a E F m Soft tssue elements Nasal Cavty or Alveol Resstance to Arflow (mmh O/cc/mn) 100 75 50 5 Irreversble Constrcton!! 0 39.5 40.0 40.5 41.0 41.5 4.0 4.5 Tme (seconds)

Acknowledgments Ths work was supported by NIH 5R01DC00730 and NIH 5P50DC007667. References 1. Bondy, J., et al., Drect ependtures related to Otts meda dagnoses: etrapolatons from a pedatrc medcad cohort. Pedatrcs, 000. 105:p. 7. Hung, C.T., et al., A paradgm for functonal tssue engneerng of artcular cartlage va appled physologc deformatonal loadng. Ann Bomed Eng, 004. 3: p. 35 49. 3. Krouskop, T.A., et al., Elastc modul of breast and prostate tssues under compresson. Ultrason Imagng, 1998. 0: p. 60 37.