$(3$ ONLY FOR CONTEXT SOLVE $4$ PLS$)$ Assume H$2$ behaves as an ideal gas with constant specific heat. A $3$ m$3$

rigid tank holds H$2$ at

$250$ kPa and $550$ K$.$ The gas is cooled down to $350$ K$.$ Without using me$40$matlab,

determine the final pressure $($kPa$)$ and the heat transfer of the system $($kJ$).$

$>>[$Pf$\_$hand,Q$\_$hand$]=$ hw$3$partc$3$

$$ Pf$\_$hand $[$kPa$]($double$)$ final pressure

$$ Q$\_$hand $[$kJ$]($double$)$ heat transfer during the process

$----------------------$

$4)$ Repeat Part C$\mathrm{.}3$ but this time with me$40$matlab. Observe the difference between this

approach and the previous manual calculation.

$>>[$Pf$\_$lib,Q$\_$lib$]=$ hw$3$partc$4$

$$ Pf$\_$lib $[$kPa$]($double$)$ final pressure

$$ Q$\_$lib $[$kJ$]($double$)$ heat transfer during the process

$----------------------$

NON EDITABLE CODE PART:

$\%$ Part C

$\%$ NOTE: finish the local functions for each subproblem in the FUNCTION$($s$)$ SECTION

fprintf$(\text{'}$

Part C

$\text{'})$;

fprintf$(\text{'}$

$1)$

$\text{'})$;

$[$Ti$,$m$\_$liq$2,$Tf$,$Wsys$]=$ hw$3$partc$1$;

fprintf$(\text{'}$Initial temperature Ti $=\%$g K

$\text{'},$Ti$)$;

fprintf$(\text{'}$Final temperature Tf $=\%$g K

$\text{'},$Tf$)$;

fprintf$(\text{'}$Mass of liquid of state $2$ m$\_$liq$2=\%$g kg

$\text{'},$m$\_$liq$2)$;

fprintf$(\text{'}$Total work Wsys $=\%$g kJ

$\text{'},$Wsys$)$;

fprintf$(\text{'}$

$2)$

$\text{'})$;

units $=$ hw$3$partc$2$;

fprintf$(\text{'}$# of required units $=\%$g

$\text{'},$units$)$;

fprintf$(\text{'}$

$3)$

$\text{'})$;

$[$Pf$\_$hand,Q$\_$hand$]=$ hw$3$partc$3$;

fprintf$(\text{'}$Final pressure Pf$\_$hand $=\%$g kPa

$\text{'},$Pf$\_$hand$)$;

fprintf$(\text{'}$Transferred heat Q$\_$hand $=\%$g kJ

$\text{'},$Q$\_$hand$)$;

fprintf$(\text{'}$

$4)$

$\text{'})$;

$[$Pf$\_$lib,Q$\_$lib$]=$ hw$3$partc$4$;

fprintf$(\text{'}$Final pressure Pf$\_$lib $=\%$g kPa

$\text{'},$Pf$\_$lib$)$;

fprintf$(\text{'}$Transferred heat Q$\_$lib $=\%$g kJ

$\text{'},$Q$\_$lib$)$;

$-----------------------------$

EDITABLE CODE PART:$($ ADD ONTO IT PLS THANKS$)$

$\%4)$

function $[$Pf$\_$lib,Q$\_$lib$]=$ hw$3$partc$4$

Pf$\_$lib$=0$; Q$\_$lib$=0$;

$\%$ Start writing your code under here

$----------------------$

attached is the me$40$library to be consulted for the problem

ME$40$ MATLAB SOFTWARE

The software package contains classes state and process and function thplot. The

following is the description of their use. You can also get the built$-$in information by typing help

state, help process, help process$/$get$,$ and help thplot in the MATLAB

command line.

Function thplot

Function thplot displays saturated lines in a graphical form. The following options are

available:

thplot default to 'water'

thplot$($substance$)$ substance: 'water',$\text{'}$r$134$a$\text{'}$

The properties displayed on the axes can be chosen from the corresponding selection boxes, and

the units can be chosen by clicking on the units with the right$-$hand button of the mouse.

Class state

A call

st $=$ state$($substance$,\{$property$1,$value$1,$units$1$;property$2,$value$2,$units$2\})$

creates an object of class state. Each such object contains the following fields:

st$.$p pressure $($in units of kPa$)$

st$.$t temperature $($in units of K$)$

st$.$v specific volume $($in units of m$3/$kg$)$

st$.$u specific energy $($in units of kJ$/$kg$)$

st$.$h specific enthalpy $($in units of kJ$/$kg$)$

st$.$s specific entropy $($in units of kJ$/$kg K$)$

st$.$x quality $($fraction$)$

st$.$phys 'liquid','sat liquid','twophase','sat vapor','vapor','ideal gas'

st$.$substance 'water', 'air', 'nitrogen', $...$

st$.$molW molecular weight $($kg$/$kmol$)$

st$.$R gas constant

Examples are:

s$1=$ state$(\text{'}$water$\text{'},\{\text{'}$p$\text{'},1\mathrm{.}23,$'MPa';$\text{'}$x$\text{'},0\mathrm{.}55,\text{'}\text{'}\})$

s$2=$ state$(\text{'}$water$\text{'},\{\text{'}$v$\text{'},0\mathrm{.}0011,\text{'}\text{'}$ ;'sat',$\text{'}\text{'},$'liquid'$\})$

s$3=$ state$(\text{'}$air$\text{'},\{\text{'}$p$\text{'},1\mathrm{.}23,$'bar';$\text{'}$h$\text{'},2345,\text{'}\text{'}\})$

The empty units field, $\text{'}\text{'},$ implies base units $($in addition to help state, typing

$>>$Units.base

displays the base units$).$

Object of class state can be a single state or a list of states, i$.$e$.,$

st $=[$ st$1$ st$2]$

or$,$ equivalently,

UCB, Spring $2024$

M$.$ Frenklach Page $2$ of $3$

ME$40,$ Thermodynamics Last updated $1/9/24$

st$(1)=$ st$1$;

st$(2)=$ st$2$;

creates an object with two states, st$1$ and st$2,$ each defined through a call to state.

To get the property values of state st$,$ use the dot syntax, for example,

st$.$p or st$.(\text{'}$p$\text{'}),$

which returns the value of p in base units $($i$.$e$.,$ kPa$).$

Once you created state object st$,$ you can display its states on a thermodynamic diagram via

plot$($st$)$ or st$.$plot

This will display the same diagram as calling thplot, but with states marked with red $$dots$.$

Clicking with right$-$hand button of the mouse on a state dot displays the state number; selecting

this box displays all the property values of the state. Selecting one of them draws a constant$-$

property line; selecting it again removes this line. Clicking with the right$-$hand button on a line

connecting states, displays the constant property of the corresponding process if exists.

Class process

Having a set of states, st $=[$st$1$ st$2...$ stN$],$ whose length$($st$)>1,$ a call

pr $=$ process$($st$)$

creates an object of class process, i$.$e$.,$ a set of connected states,

st$(1)->$ st$(2)->...->$ st$($end$).$

When the initial state is the same as the final state, i$.$e$.,$ st $=[$st$1$ st$2...$ st$1],$ the call to

process makes a cycle,

st$(1)->$ st$(2)->...->$ st$(1)$

A call with multiple state sets, e$.$g$.,$

pr $=$ process$([$st$1$ st$2$ st$3$ st$4$ st$5$ st$6],[$st$4$ st$2])$

creates a split$-$flux process,

st$(1)->$ st$(2)->...->$ st$(6)$ with st$4$ also connected to st$2.$

You can get process properties by the following calls to get$($pr$,..)$ :