{"id":716,"date":"2016-05-16T15:45:33","date_gmt":"2016-05-16T20:45:33","guid":{"rendered":"http:\/\/rpchurchill.com\/?p=716"},"modified":"2019-01-24T10:52:32","modified_gmt":"2019-01-24T15:52:32","slug":"a-specific-example-of-equations-to-be-solved-by-matrix","status":"publish","type":"post","link":"https:\/\/rpchurchill.com\/wordpress\/posts\/2016\/05\/16\/a-specific-example-of-equations-to-be-solved-by-matrix\/","title":{"rendered":"A Specific Example of Equations To Be Solved By Matrix"},"content":{"rendered":"<p>Dividing a cross-section of a steel billet into nodes, we can start by looking at the energy balance for each node:<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.rpchurchill.com\/images\/articles\/ThermalMatrixEqn1.png\" width=\"500px\" \/><\/p>\n<p>All this equation is saying is that the change in the quantity of energy in a given node per unit time is equal to the sum of the energy flows into and out of the node.  The quantity to the left is the accumulation terms while the quantities on the right side of the equation are terms for radiative heat transfer from the environment to the surface of the steel, conductive heat transfer between nodes within the steel, convective heat transfer between the external atmosphere and the steel, and inductive heat transfer into the steel.<\/p>\n<p>Let\u2019s start by breaking the mass in the accumulator term down as follows:<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.rpchurchill.com\/images\/articles\/ThermalMatrixEqn2.png\" width=\"500px\" \/><\/p>\n<p>We continue by eliminating the convection and induction terms because they are comparatively small or nonexistent, and then by dividing all the remaining terms by the unit area:<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.rpchurchill.com\/images\/articles\/ThermalMatrixEqn3.png\" width=\"500px\" \/><\/p>\n<p>Next we can expand terms:<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.rpchurchill.com\/images\/articles\/ThermalMatrixEqn4.png\" width=\"500px\" \/><\/p>\n<p>Steal a first order dependent variable from the fourth order term:<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.rpchurchill.com\/images\/articles\/ThermalMatrixEqn5.png\" width=\"500px\" \/><\/p>\n<p>Substituting back into the previous equation yields:<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.rpchurchill.com\/images\/articles\/ThermalMatrixEqn6.png\" width=\"500px\" \/><\/p>\n<p>Gathering terms we get:<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/www.rpchurchill.com\/images\/articles\/ThermalMatrixEqn7.png\" width=\"500px\" \/><\/p>\n<p>The coefficients for the T<sub>n<\/sub> terms and the T<sub>c<\/sub> terms on the left side of the equation fill the initialized coefficient matrix while the terms on the right side of the equation fill the initialized source matrix.  There is one T term for each node.  The radiation transfer terms drop out for any node that isn\u2019t on the surface of the piece of steel being heated.  The conduction terms drop out between any two nodes that aren\u2019t internally connected.<\/p>\n<p>This equation gets written for each node individually and then the system of equations gets solved.  All of the T<sub>n<\/sub> terms for the new time step end up in the source matrix at the end.  Note that the temperatures are all represented in absolute degrees.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dividing a cross-section of a steel billet into nodes, we can start by looking at the energy balance for each node: All this equation is saying is that the change in the quantity of energy in a given node per &hellip; <a href=\"https:\/\/rpchurchill.com\/wordpress\/posts\/2016\/05\/16\/a-specific-example-of-equations-to-be-solved-by-matrix\/\">Continue reading <span class=\"meta-nav\">&rarr;<\/span><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[14],"tags":[98,223,77],"_links":{"self":[{"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/posts\/716"}],"collection":[{"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/comments?post=716"}],"version-history":[{"count":3,"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/posts\/716\/revisions"}],"predecessor-version":[{"id":1506,"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/posts\/716\/revisions\/1506"}],"wp:attachment":[{"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/media?parent=716"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/categories?post=716"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/rpchurchill.com\/wordpress\/wp-json\/wp\/v2\/tags?post=716"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}