MSC Nastran; Patran Tutorial - Transient Heat Transfer with Phase Change

in this example we'll be performing a
transient heat transfer analysis with a
phase change full details of this
exercise are on page 450 of the PDF
linked in the video description below
I'll make a new folder on my desktop
called problem 28 I'll start Patrat the
units for this exercise are meters watts
degrees Celsius joules and kilograms
click new
I move down and select that folder I
just made
call it prop 28 and here I have to
change the analysis type to thermal I'll
click OK the first thing I'll do after
that I'll create a cube so here under
geometry go to solids and hit XYZ this
cube is 0.1 in every direction after you
hit after you enter this vector use the
same after you enter this factor use
this origin point and hit apply so
here's your cube if I go to the Home tab
and click shaded here I can review the
cube shaded now let's go ahead and
define our material by going to the
properties tab and clicking isotropic
here I'll call this material matte under
input properties I'll give it a thermal
conductivity of 0.6 a specific heat of
four-to-one 7 and density of 1000 click
OK and hit apply go back to input
properties switch this to phase changes
here under phase change temperature
range type in 2 and my lane heat value
is 3 3 4 0 0 0 click OK hit apply
I'll go ahead and assign this material
to the solid by clicking 3d property
solid call this prop under input
properties also like the materia just
made click OK and notice it's a entered
here click OK for application region
select this solid one added ok and apply
now going on to our boundary conditions
I'll apply an initial temperature by
going to initial conditions temperature
I'll call this 20 20 C actually for my
input data initial temperature will be
20 click OK or your application region
will use the point filter tool here make
sure the cursor is here and click pick
also it selects all the points add it ok
and apply now I'll apply element uniform
convection condition here I'll just type
in convection and this will eventually
be a 3d element so click input data here
for convection coefficient type in 100
for ambient temperature type in negative
20 click OK and for your application
region this will be applied to all the
phases that are free to the bare
essentially so here click basal solid
and select or pickle add it ok and apply
you have the markers here next we'll
mesh this so select the solid mushroom
here you have a few options you can
select Ted wedge or hex
if I select at and I select the solid
this is a what the mesh will look like
if I do a wedge this is what that will
look like but in this case we just want
X with the global edge length of 0.1 or
if I were to make it point zero five you
will see eight elements one two three
four five six seven eight but in this
case we're just doing one element so
give it a global edge link that matches
the width length and I know this cube
now we just have one element and one
thing I neglected to do is for my load
bc's I needed to define a time-dependent
load case I'll do that now so create a
low case call it time underscore LC make
it time-dependent for your input data
just add the convection then the initial
condition we created a moment ago hit
apply here you can minimize this and go
back home and click model tree I'll just
go over the boundary conditions making
sure they're properly defined and show
temperature that is correct
convection that is correct so now I can
move on to my analysis click analyze
entire model your solution types should
be selected to steady state click OK
here we'll go ahead and modify one of
our sub cases will modify the time LC
sub case or load case we just created
under sub case parameters the number of
increments will be
or actually cancel this go back to
solution type select transient analysis
click OK and your sub cases modify the
time LC sub case on your sub case
parameters give it an initial time step
of 5 the number of steps this 1000 click
ok here for output request we want two
things we want the temperatures which
already selected here you'll see a
thermal and we want the rate of exchange
of enthalpies so here at the bottom you
select that in it includes each dot now
you can click OK here at the bottom or
just click apply out here click cancel
you'll notice the we get feedback at the
sub case has been modified for sub case
select make sure you select time L see
if you select default it's a empty sub
case with no boundary conditions in its
respective load case click OK here and
select apply once that's done and put
the results and click apply the results
tab select the graph tool click this
icon select all the results and now what
we want to view here are our our
temperatures so here for the node we
want node 1 and to make this easier on
this we'll go to the Home tab turn on
label control and turn on the node
numbers we want to see it for node 1 and
hit apply
I'll get this graph when I compare it to
what we should get and it's actually
what we do get so this is the
temperature response of the freezing
cube for our next graph we'll view the
rate of enthalpy change so back here in
the results tab
we'll see rates of the enthalpy change
here simply select that and hit apply OC
it modify and you'll see that the ring
or the the range is from thirty to zero
so it goes from thirty to zero so we've
obtained the same answers make sure to
save and one last tip you can hide the
plot by just unchecking here and that
concludes this example