Numerical Modeling of Wave Propagation in Particulate ...

18
Modeling of Wave Propagation in Particulate Composites Phani Mylavarapu, Srivathsa Boddapati Defence Metallurgical Research Laboratory, Kanchanbagh 1 17-11-2011

Transcript of Numerical Modeling of Wave Propagation in Particulate ...

Page 1: Numerical Modeling of Wave Propagation in Particulate ...

Modeling of Wave Propagation in Particulate

Composites

Phani Mylavarapu, Srivathsa Boddapati

Defence Metallurgical Research Laboratory,

Kanchanbagh

1 17-11-2011

Page 2: Numerical Modeling of Wave Propagation in Particulate ...

Outline • Introduction • Problem Definition • FEM Modeling • Approach • Results and Discussions • Conclusions and Future work

2 17-11-2011

Page 3: Numerical Modeling of Wave Propagation in Particulate ...

Syntactic foams (SF)

Two phase Three phase

microballoons

matrix

void

SEM Image of Syntactic Foams

ηi= ri/ro

t d0

Underwater explorer

Courtesy: Oceanexplorer 3 17-11-2011

Page 4: Numerical Modeling of Wave Propagation in Particulate ...

Ultrasonic Testing

• Attenuation

• Ultrasonic Velocity – Longitudinal Velocity (𝑉𝑙) – Shear Velocity (𝑉𝑠)

E. Segal and J. L. Rose, in: Research Techniques in Non-Destructive Testing, Vol. IV,R. S. Sharpe (Ed.), pp. 275–316. Academic Press, London (1980).

xeA

A 0

ras

2* lVL 2* sVG GL

GLGE

)43(

)(22GL

GL

4 17-11-2011

Page 5: Numerical Modeling of Wave Propagation in Particulate ...

Modulus Predicted from UT Sample Name

Predicted UI Modulus, MPa

SHPB Modulus, MPa

Quasi -Static Modulus,

Mpa

Pure Epoxy 4476±183 4272-5084 2320±40

S2210 3667±89

3310-3903 2156

S2250 2173±41

1701-2244 1512

S3210 4071±45

3048-4099 2310

S3250 2800±88

2146-2808 2052

S3840 3481±25

3244-3664 2351

K4650 3385±109

2901-3446 2473

5 17-11-2011

Page 6: Numerical Modeling of Wave Propagation in Particulate ...

6

Cluster Interaction With Ultrasonic Energy

Page 7: Numerical Modeling of Wave Propagation in Particulate ...

Problems with Ultrasonic Testing

• Difficult to distinguish the reflections originating from epoxy or particle or epoxy/particle interface

• Waves due to multiple reflections

S2230

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0 5000 10000 15000 20000 25000 30000 35000

Time, ns

Am

pli

tud

e, V

olt

s

First back

wall reflection Second back

wall reflection Third back

wall reflection

Front wall reflection

7 17-11-2011

1a

1a

1a

Typically UT performed

& Experimentally difficult

Page 8: Numerical Modeling of Wave Propagation in Particulate ...

Modeling using COMSOL

• BC’s available for source of excitation: – Gaussian point source – Pressure – Plane wave radiation

• BC’s available for reflection/transmission between interfaces – Impedance

17-11-2011 8

Qcp

pqp

p oo

2

21. 01

qpn t

Sound Hard Boundary

t

p

zqpn t

i

t

11.

Page 9: Numerical Modeling of Wave Propagation in Particulate ...

Mesh Details and CFL • 0.285 element size • ƛ/10 • Same element size for all models • CFL Criteria: 5e-9 s

17-11-2011 9

Page 10: Numerical Modeling of Wave Propagation in Particulate ...

Case 1: 1a

10 17-11-2011

Page 11: Numerical Modeling of Wave Propagation in Particulate ...

Case 1 (contd…)

11 17-11-2011

Page 12: Numerical Modeling of Wave Propagation in Particulate ...

Case 1 (contd…)

12 17-11-2011

Page 13: Numerical Modeling of Wave Propagation in Particulate ...

Case 1 (contd…)

13 17-11-2011

Page 14: Numerical Modeling of Wave Propagation in Particulate ...

17-11-2011 14

1a

Page 15: Numerical Modeling of Wave Propagation in Particulate ...

17-11-2011 15

Comparison of pressure in the particle

1a

1

a

Page 16: Numerical Modeling of Wave Propagation in Particulate ...

17-11-2011 16

1a

Page 17: Numerical Modeling of Wave Propagation in Particulate ...

Conclusions and Future Work • Effect of ratio of wavelength to particle size on

reflection studied. • Ratios of less than, greater than and equal to 1. • Multiple reflections from particles

17-11-2011 17

Page 18: Numerical Modeling of Wave Propagation in Particulate ...

Acknowledgements

• Director, DMRL • Dr. Nikhil Gupta, Brooklyn University

17-11-2011 18