MATLABAssignment
MATLABAssignment
MATLABAssignment
Page 2 of 5
Dr
A.
Cookson
Department
of
Mechanical
Engineering
2018/2019
University
of
Bath
Your
finite
element
mesh
and
material
parameters
must
represent
the
three-‐
layered
structure
of
skin,
as
follows:
Parameter
values
that
you
should
use
in
each
layer
are
provided
in
Table
1
–
with
the
exception
that
in
this
question
we
assume
zero
blood
flow
i.e.
G=0
everywhere.
The
layers
of
the
tissue
are
defined
at
the
following
x
coordinates:
E=0.00166667,
D=0.005,
B=0.01.
a) Run
your
code
for
the
following
initial
conditions
and
Dirichlet
boundary
conditions,
for
a
maximum
of
50
seconds:
T (x, 0) = 310.15K, T (x = B, t) = 310.15K, T (x = 0, t) = 393.15K
<latexit sha1_base64="S+gsjDMsRLnFMWkh0JHLvCu3rDQ=">AAACMXicdZDLSgMxFIYz9VbrbdSlm2ARKpSSsYi6KJS6Kbip0IvQDiWTpm1o5mKSEcvQV3Ljm4ibLhRx60uYmXahrR4I/PzfOZyc3wk4kwqhqZFaWV1b30hvZra2d3b3zP2DpvRDQWiD+NwXdw6WlDOPNhRTnN4FgmLX4bTljK5j3nqgQjLfq6txQG0XDzzWZwQrbXXNaj33mEensASLFipY5/AmDzv3Ie5BDUqVPFT/MZSfoatigrpmFhVQUnBZWHORBfOqdc2XTs8noUs9RTiWsm2hQNkRFooRTieZTihpgMkID2hbSw+7VNpRcvEEnminB/u+0M9TMHF/TkTYlXLsOrrTxWooF1ls/sXaoepf2hHzglBRj8wW9UMOlQ/j+GCPCUoUH2uBiWD6r5AMscBE6ZAzOgRr8eRl0TwrxHneomy5Mo8jDY7AMcgBC1yAMqiCGmgAAp7AK3gD78azMTU+jM9Za8qYzxyCX2V8fQNUE6CX</latexit>
sha1_base64="/3jbs8nrmRVJwXnYeTlBGQWw8oc=">AAACMXicdZBLS0JBFMfPtZfZ61bLNoMSGIncm0S1EMQ2QhsDX6Aic8dRB+c+mpkbifhd+gRt+ibRxkURbfsSjY9FaR0Y+PP/ncOZ83cCzqSyrLERWVldW9+Ibsa2tnd298z9g4r0Q0FomfjcFzUHS8qZR8uKKU5rgaDYdTitOv3rCa/eUyGZ75XUIKBNF3c91mEEK221zEIp+ZCyTlAWZWwrbZ+jmxRq3IW4jTTI5lNI/ces1AxdZaaoZSastDUttCzsuUjk4o3TRwAotsyXRtsnoUs9RTiWsm5bgWoOsVCMcDqKNUJJA0z6uEvrWnrYpbI5nF48QsfaaaOOL/TzFJq6PyeG2JVy4Dq608WqJxfZxPyL1UPVuWwOmReEinpktqgTcqR8NIkPtZmgRPGBFpgIpv+KSA8LTJQOOaZDsBdPXhaVs/Qkz1udRh5mFYUjiEMSbLiAHBSgCGUg8ASv8AbvxrMxNj6Mz1lrxJjPHMKvMr6+AWaSoiA=</latexit>
sha1_base64="deDdzdjGHZafshfKucR35742ZyY=">AAACMXicdZDLSgMxFIYz9VbrbdSlm9AiVBzKjEXURaHUTcFNhd6gLSWTZtrQzMUkI5bSt/A53Pgm4qYLRdz6EmamXWirBwI//3cOJ+e3A0aFNM2pllhZXVvfSG6mtrZ3dvf0/YO68EOOSQ37zOdNGwnCqEdqkkpGmgEnyLUZadjD64g37gkX1PeqchSQjov6HnUoRlJZXb1czT4Y5gkswLxl5qxzeGPA9l2IelCBQsmA8j9mGjN0lY9RV8+YOTMuuCysucgU0+3Tx2lxVOnqL+2ej0OXeBIzJETLMgPZGSMuKWZkkmqHggQID1GftJT0kEtEZxxfPIHHyulBx+fqeRLG7s+JMXKFGLm26nSRHIhFFpl/sVYoncvOmHpBKImHZ4uckEHpwyg+2KOcYMlGSiDMqforxAPEEZYq5JQKwVo8eVnUz3JRnrcqjRKYVRIcgTTIAgtcgCIogwqoAQyewCt4A+/aszbVPrTPWWtCm88cgl+lfX0Db/yjpg==</latexit>
Create
a
plot
showing
the
spatial
temperature
distribution
in
the
tissue
for
various
points
in
time
between
0
and
50
seconds.
Explain
the
shape
of
the
curves
as
they
change
through
time
based
on
your
knowledge
of
the
physics,
initial
conditions,
and
boundary
conditions.
b)
Once
you
have
solved
the
temperature
distribution
in
space
and
time,
use
this
information
to
determine
the
level
of
tissue
damage
that
will
result.
To
do
this,
evaluate
the
following
integral
numerically
at
the
mesh
node
for
x=E:
t Z ✓ ◆
98 12017
= 2 ⇥ 10 exp dt
tburn (T 273.15)
<latexit sha1_base64="42sVyH86fqDkPq4E9nzrj4BTk4A=">AAACRHicbZDNbhMxFIU95a+EvwBLNhYRUrpoZKdCKYtKFSxgWaSmrRSnkce5k1i1PSP7DmpkzcN10wdgxxOwYQFCbBFOmgW0XNnS0Tn3ytdfXhkdkLEv2cat23fu3tu833rw8NHjJ+2nz45CWXsFQ1Wa0p/kMoDRDoao0cBJ5UHa3MBxfvZumR9/Ah906Q5xUcHYypnThVYSkzVpj8R7aa2ke1Roh5OYTl571zSnERvapwK1hUA5O41vdhsq4LwSBgrsbovCSxV5n/FBE7uH2/3BTo+/3mqE17M5btEpTtod1mOrojcFX4sOWdfBpP1ZTEtVW3CojAxhxFmF4yg9amWgaYk6QCXVmZzBKEkn02rjuILQ0FfJmdKi9Ok6pCv374kobQgLm6dOK3EermdL83/ZqMZidxy1q2oEp64eKmpDsaRLonSqPSg0iySk8jrtStVcJjiYuLcSBH79yzfFUb/HWY9/ZJ39t2scm+QFeUm6hJMB2ScfyAEZEkUuyFfynfzILrNv2c/s11XrRraeeU7+qez3H0PHrxg=</latexit>
sha1_base64="2KX0IRJtk+S+0uqZGcHmcGCNlxs=">AAACRHicbZDNThsxFIXvUNrS9IfQLruxQJXCgmgchAKLSogu2iWVCCDFIfI4nsTC9ozsO1Ujax6ODQ/Ark/QTRetKraoTsKCn17Z0tE598rXX1Zq5TFNfyRLT5afPnu+8qLx8tXrN6vNtbfHvqickD1R6MKdZtxLrazsoUItT0snucm0PMnOP83yk2/SeVXYI5yWcmD42KpcCY7RGjb77DM3hpOPhCmLwxBPVjlb12cBa9IhDJWRntD0LOzt1oTJ7yXTMsfWFssdF4F2UtqtQ+toq9PdbtOdzZo5NZ7gJhnhsLmRttN5kceC3oqN/YO8IgBwOGxesVEhKiMtCs2979O0xEHgDpXQsm6wysuSi3M+lv0oLY+rDcIcQk0+RGdE8sLFa5HM3bsTgRvvpyaLnYbjxD/MZub/sn6F+e4gKFtWKK1YPJRXmmBBZkTJSDkpUE+j4MKpuCsREx7hYOTeiBDowy8/FsedNk3b9OuMBixqBd7DOrSAQhf24QscQg8EXMBP+A1/ksvkV/I3uV60LiW3M+/gXiU3/wBG5LCV</latexit>
sha1_base64="tzOz+MPQxdZIpKPAb0paT/Q1CBQ=">AAACRHicbZDPbhMxEMa9pdAS/jSFIxysVkjpodE6VZVyQKrgQI9FatpKcRp5ndnEqu1d2bOIyNpX4J248AAckLhy4dJDK8QV4SQ9QMvIlj5934w8/mWlVh7T9FuydGf57r2V1fuNBw8fPV5rrj859kXlJPRkoQt3mgkPWlnooUINp6UDYTINJ9n5m1l+8h6cV4U9wmkJAyPGVuVKCozWsNnnb4Uxgr6iXFkchniyytm6PgtY0w7lqAx4ytKz8HKvphw+lFxDjq1tnjshA+ukrFuH1tF2p7vTZrtbNXdqPMEtOsJhczNtp/OitwW7Fpv7r/Pq+fePXw+HzS98VMjKgEWphfd9lpY4CMKhkhrqBq88lEKeizH0o7QirjYIcwg1fRGdEc0LF69FOnf/ngjCeD81Wew0Aif+ZjYz/5f1K8z3BkHZskKwcvFQXmmKBZ0RpSPlQKKeRiGkU3FXKiciwsHIvREhsJtfvi2OO22Wttm7GQ2yqFXyjGyQFmGkS/bJATkkPSLJJ/KDXJKr5HNykfxMfi1al5Lrmafkn0p+/wFLQbLX</latexit>
Integrate
between
the
time
points,
tburn,
at
which
the
temperature
T
becomes
greater
than
317.15
Kelvin,
and
time,
t
(which
will
be
50
seconds
for
all
these
questions).
Note
that
the
Matlab
function,
trapz(x),
uses
the
trapezium
rule
to
integrate
the
vector
x,
assuming
an
interval
of
1.
Therefore
multiply
the
output
by
the
size
of
your
timestep,
Δt,
to
obtain
the
final
value
of
this
integral.
Page 3 of 5
Dr
A.
Cookson
Department
of
Mechanical
Engineering
2018/2019
University
of
Bath
If
Γ >
1
at
x=E,
this
produces
a
second-‐degree
burn.
State
the
value
of
Γ
you
have
calculated
at
this
node
(don’t
be
alarmed
if
it
is
very
large
or
small),
and
hence
state
whether
you
expect
a
second-‐degree
burn
to
occur.
2. Now
use
your
code
to
determine
the
minimum
temperature
reduction
(to
the
nearest
0.5
degree
C)
that
must
be
achieved
by
the
protective
clothing
at
the
boundary
x=0,
in
order
that
second-‐degree
burns,
as
measured
at
x=E,
are
not
caused
by
t=50s.
State
the
final
Dirichlet
boundary
condition
at
x=0
that
achieves
this,
and
explain
how
you
used
your
code,
and
any
other
calculations,
to
estimate
this
value.
Include
any
data
or
figures
that
you
feel
are
relevant
in
order
to
demonstrate
your
approach
to
estimating
this.
3. Re-‐run
your
results
for
question
2,
but
now
including
the
blood
flow
related
reaction
&
source
terms
–
use
the
values
in
Table
1
for
each
region
of
tissue.
How
much
does
this
change
the
temperature
distribution
in
space
and
time?
Do
you
consider
this
to
be
an
important
effect
to
consider
for
future
modelling
of
this
problem?
4. Other
than
changing
the
geometry
to
a
3D
representation
of
the
body,
comment
briefly
on
the
accuracy
of
the
modelling
assumptions
and
suggest
how
these
could
be
modified
in
order
to
provide
more
realistic
results?
Parameter
Epidermis
Dermis
Sub-‐cutaneous
k
25
40
20
G
0
0.0375
0.0375
ρ
1200
1200
1200
c
3300
3300
3300
ρb
-‐
1060
1060
cb
-‐
3770
3770
Tb
-‐
310.15
310.15
Table
1:
Parameter
values
for
questions
2,
3
and
4.
Note
that
some
are
realistic
and
others
are
not,
having
been
chosen
to
allow
you
to
solve
your
model
in
a
reasonable
time
frame.
Page 4 of 5
Dr
A.
Cookson
Department
of
Mechanical
Engineering
2018/2019
University
of
Bath
SUBMISSION
GUIDELINES
Structure
your
report
as
a
set
of
answers
to
these
questions
–
there
is
no
requirement
to
write
this
in
a
lab
report
format.
However
your
report
must
be
self-‐contained
and
therefore
must
not
assume
that
the
reader
knows
the
content
in
this
document.
• You
must
include
all
your
Matlab
source
code
as
text
in
the
Appendices.
Do
not
paste
your
code
into
the
document
as
an
OLE
or
as
an
image.
• Do
not
upload
zipped
or
compressed
folders
of
these
source
files.
• Your
code
should
use
meaningful
variable
names
and
include
comments,
in
line
with
good
practice.
• Word
limit
of
2000
words
(not
including
source
code)
Submit
your
work
using
the
online
submission
function
on
the
unit’s
Moodle
page.
Deadline:
4pm
on
Monday,
3rd
December
2018.
Marks
will
be
awarded
based
on
the
following
criteria:
• Correctness
of
numerical
results
• Demonstration
of
knowledge
of
material
presented
in
the
course
including
more
advanced
features
of
the
finite
element
method
• Clear
presentation
of
results
• Correctness
&
readability
of
Matlab
code
• Level
of
insight
&
clarity
of
explanation
of
the
meaning
of
the
results
Page 5 of 5