Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F84205230
plot_P_vs_I.py
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Sat, Sep 21, 08:48
Size
1015 B
Mime Type
text/x-python
Expires
Mon, Sep 23, 08:48 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
20957610
Attached To
R4670 PySONIC (old)
plot_P_vs_I.py
View Options
#!/usr/bin/env python
# -*- coding: utf-8 -*-
# @Author: Theo Lemaire
# @Date: 2017-07-17 11:47:50
# @Email: theo.lemaire@epfl.ch
# @Last Modified by: Theo Lemaire
# @Last Modified time: 2018-08-21 16:07:32
''' plot profile of acoustic Intensity (in W/cm^2) vs Pressure (in kPa) '''
import
numpy
as
np
import
matplotlib.pyplot
as
plt
from
PySONIC.utils
import
Pressure2Intensity
rho
=
1075
# kg/m3
c
=
1515
# m/s
P
=
np
.
logspace
(
np
.
log10
(
1e1
),
np
.
log10
(
1e7
),
num
=
500
)
# Pa
Int
=
Pressure2Intensity
(
P
,
rho
,
c
)
# W/m2
fig
,
ax
=
plt
.
subplots
()
ax
.
set_xlabel
(
'$Pressure\ (kPa)$'
)
ax
.
set_ylabel
(
'$I_{SPPA}\ (W/cm^2)$'
)
ax
.
set_xscale
(
'log'
)
ax
.
set_yscale
(
'log'
)
ax
.
plot
(
P
*
1e-3
,
Int
*
1e-4
)
Psnaps
=
np
.
logspace
(
1
,
7
,
7
)
# Pa
for
Psnap
in
Psnaps
:
Isnap
=
Pressure2Intensity
(
Psnap
,
rho
,
c
)
# W/m2
ax
.
plot
(
np
.
array
([
Psnap
,
Psnap
])
*
1e-3
,
np
.
array
([
0.0
,
Isnap
])
*
1e-4
,
'--'
,
color
=
'black'
)
ax
.
plot
(
np
.
array
([
0
,
Psnap
])
*
1e-3
,
np
.
array
([
Isnap
,
Isnap
])
*
1e-4
,
'--'
,
color
=
'black'
)
plt
.
show
()
Event Timeline
Log In to Comment