|
|
- Homer Booth
- 5 years ago
- Views:
Transcription
1
2
3
4
5
6
7 ME201/MTH21/ME400/CHE400 ASSIGNMENT #7 PROBLEM 1 Newton's Law of Cooling in Transient Conduction 1. Introduction This notebook uses Mathematica to solve problem 1 of Assignment #7, in which we analyze transient heat conduction in a slab of width L. The boundary condition on the right face of the slab is a specified constant heat flux F 0, and the boundary condition on the left face of the slab is a convective flux, described by Newton's law of cooling, with an ambient temperature T A. The initial temperature is constant and is equal to the ambient temperature T A. The relevant parameters are the thermal conductivity k, the thermal diffusivity D f, and the heat transfer coefficient h, all assumed constant. The specific numerical expressions for all functions and parameters are given in section 5. The mathematical formulation of the problem is given below. with k T T t = D f 2 T x H0, tl = h@t H0, tl - T AD, k T, 0 < x < L, t > 0, 2 x x HL, tl = F 0, and T Hx, 0L = T A. (1) 2. Steady State Solution As in other problems of this sort, we must first find the steady-state solution T s (x), and then reformulate the problem in terms of the transient solution T r. The steady-state solution is obtained by setting the time derivative to zero in the original equation for T. The solution of the resulting equation is a linear function of x, and we require that particular linear function which satisfies the boundary conditions at x = 0 and L. The result is easily shown to be In[1]:= Ts@x_D := TA + F0 h + F0 k x We check this. Because it is linear, it clearly satisfies the equation. Now the boundary conditions: In[2]:= Out[2]= k D@Ts@xD, xd ê. x Ø L F0 In[3]:= Hk D@Ts@xD, xd - h HTs@xD - TALL ê. x Ø 0 Out[3]= 0 It checks. 3. Formulation of Problem for Transient Now that we have the steady-state solution, we decompose the full solution into a steady-state part T s (x) and a transient part T r (x,t):
8 2 sol7.nb T Hx, tl = T s HxL + T r Hx, tl. (2) By substituting this decomposition into the original equation for T, we find the following problem for the transient: T r t k T r x H0, tl - h T r H0, tl = 0, This is the problem that we solve by separation of variables. 2 T r = D f, 0 < x < L, t > 0, 2 x T r x HL, tl = 0, and T r Hx, 0L = T 0 HxL - T s HxL. (3) 4. Separation of Variables in Transient Problem Following the standard approach in separation of variables, we seek functions which satisfy the equation for T r and the homogeneous boundary conditions at x = 0, L. We assume a form F(x)G(t) for such solutions. The two separated equations then have the form F ² HxL + lf HxL = 0, with kf H0L - hf H0L = 0, F HLL = 0 and and G HtL + ldg HtL = 0. (4) With the Rayleigh quotient, one may show that all of the eigenvalues are positive. Then the general solution of the equation for F(x) is a linear combination of sin( l x) and cos( l x). Applying the two homogeneous boundary conditions to the general solution, we get a solution which may be put into the form F(x) = cos( l H1 -x)). Simplifying the resulting eigenvalue equation, we get tan HzL = B i êz, where z = L l, and where B i = hlêk. (5) We define the relevant quantities for Mathematica, denoting the nth value of z by z[n], and the nth eigenvalue by lam[n]. In[4]:= lam@n_d := Hz ê LL^2 ê. z -> z@nd where the values z[1], z[2],.. z[n],.. are the roots of eig1 = eig2, where In[5]:= eig1 = Tan@zD; In[]:= eig2 = Bi ê z; with In[7]:= Bi := Hh * LL ê k The generic eigenfunction is In[]:= geneig@x_d := Cos@z H1 - x ê LLD and the nth eigenfunction is given by substituting for z the nth eigenvalue z(n): In[9]:= Feig@x_, n_d := geneig@xd ê. z -> z@nd The generic normalization integral is In[]:= Out[]= nor = Integrate@Hgeneig@xDL ^ 2, x, 0, L<D L Hz + Cos@zD Sin@zDL 2 z
9 sol7.nb 3 and the nth normalization integral is In[11]:= norm@n_d := nor ê. z -> z@nd 5. Parameter Values In this section, we define the parameter values and the initial function. This is the only place in the notebook where these quantities are given specifically. This makes it possible to change a value for the entire calculation by just changing it here. The material properties are appropriate for granite. The initial temperature is taken here to be a constant. We also calculate the value of the Biot number Bi. In[]:= In[13]:= In[]:= In[15]:= In[1]:= In[17]:= h = 22.4; H** Wëm 2 ÿk **L L = 0.75; H** m **L k = 1.2; H** WêmÿK **L Df = 4.1 * - ;H** m 2 ës **L TA =.0; H** C **L T0@x_D =.0; H** C **L In[1]:= F0 = 170.2; H** Wëm 2 **L In[19]:= Bi Out[19]= Determination of Eigenvalues We now determine numerically the eigenvalues. We begin by plotting the expressions eig1 and eig2. At each crossing of these, there is an eigenvalue. In[20]:= SetOption@Plot, ImageSize Ø 250D; In[21]:= Plot@eig1, eig2<, z, 0, <, PlotRange -> -, <, PlotStyle -> RGBColor@1, 0, 0D, RGBColor@0, 0, 1D<, AxesLabel -> "z", "eig1 HredL and eig2 HblueL"<, Ticks Ø 0, Pi ê 2, Pi, 3 Pi ê 2, 2 Pi, 5 Pi ê 2, 3 Pi<<D eig1 HredL and eig2 HblueL 5 Out[21]= p 2 p 3 p 2 2 p 5 p 2 3 p z -5 -
10 4 sol7.nb We see that the first root is between 0 and p/2, and each subsequent root is in an interval of length p. It is also clear that for large n, the nth root is very close to (n-1)p. To find accurate numerical values of the z's, we use the Mathematica command FindRoot. FindRoot requires an initial guess. The nth value of z, as we can see from the graph, is between (n - 1) and (n -0.5 )p. We take as an initial guess (n - 0.5)p The command is then of the form FindRoot[eig1 == eig2,{z, (n - 0.5)p -0.1}]. We now use this command in a Do loop to calculate and display in a table the first 40 roots. This will suffice for any calculation requiring up to 40 terms in the eigenfunction series. We also display the eigenvalue l, and the approximate value z = (n-1)p. In[22]:= In[23]:= In[24]:= zroot@n_d := FindRoot@eig1 == eig2, z, Hn - 0.5L * Pi - 0.1<D zasymp@n_d := N@Hn - 1L * PiD Do@z@nD = z ê. zroot@nd, n, 1, 40<D;
11 sol7.nb 5 In[25]:= TableForm@Table@n, PaddedForm@z@nD,, 4<D, PaddedForm@zasymp@nD,, 4<D, PaddedForm@lam@nD,, 4<D<, n, 1, 40<D, TableHeadings -> None, "n", " z@nd", " zasymp@nd", " lam@nd"<<d Out[25]//TableForm= n z@nd zasymp@nd lam@nd We see that the asymptotic formula for z[n] has an error of less than 0.5% for any n equal to or greater than. It is possible to develop even more accurate, but still simple, asymptotic formulas. 7. Representation of the Initial Condition We are finally ready to begin solving the boundary value for the transient T r (x,t). The form of the solution for the transient is
12 sol7.nb T r Hx, tl = CHnL ExpI-lHnL D f tm cos@hzhnl H1 - xêlld, n=1 () where the coefficients C(n) are determined by the initial conditions satisfied by T r. The formula for the nth coefficient is In[2]:= coeff@n_d := Integrate@Feig@x, nd * HT0@xD - Ts@xDL, x, 0, L<D ê norm@nd We now use a Do loop to construct the first 40 coefficients, and then print them out in a table. The nth numerical value is assigned to c[n]. In[27]:= Do@c@nD = Re@coeff@nDD, n, 1, 40<D; In[2]:= TableForm@ Table@n, PaddedForm@c@nD,, 4<D, PaddedForm@n + 20, D, PaddedForm@c@n + 20D,, 4<D<, n, 1, 20<D, TableHeadings -> None, "n", " c@nd", " n", " c@nd"<<d Out[2]//TableForm= n c@nd n c@nd As a check on all that we have done, we see how well our series represents the initial condition, by plotting both the initial condition (in blue) and its representation by the first 40 terms of our series (in red). In[29]:= exactinit@x_d := T0@xD - Ts@xD In[30]:= seriesinit@x_d := Sum@c@nD * Feig@x, nd, n, 1, 40<D
13 sol7.nb 7 In[31]:= seriesgraph = Plot@exactinit@xD, seriesinit@xd<, x, 0, L<, PlotRange -> -30, 0<, PlotStyle -> RGBColor@0, 0, 1D, Thickness@0.004D<, RGBColor@1, 0, 0D, Thickness@0.004D<<, AxesLabel -> "x", "Exact HblueL, Series HredL"<D Exact HblueL, Series HredL x - Out[31]= We see that the graphical agreement is excellent. We also check a few values. In[32]:= Module@x, n<, x := 0.05 * n; TableForm@ Table@x, PaddedForm@exactinit@xD,, 4<D, PaddedForm@seriesinit@xD,, 4<D<, n, 0, 15<D, TableHeadings -> None, "", " Exact IC", " Series IC"<<DD Out[32]//TableForm= Exact IC Series IC Again excellent agreement, and we conclude that we are using enough terms (40) in the series.. Solution of the Initial-Boundary Value Problem Now that we have the eigenvalues and series coefficients, we can construct the series solution of the problem. We define Tran[x,t,n] as the nth partial sum of the series solution, and we define term[x,t,k] as the kth term in the series. As a special case of importance, we define firstterm[x,t] to be the first term in the series. In[33]:= In[34]:= In[35]:= term@x_, t_, m_d := c@md * Feig@x, md * Exp@-lam@mD * Df * td Tran@x_, t_, n_d := Sum@term@x, t, rd, r, 1, n<d firstterm@x_, t_d := term@x, t, 1D
14 sol7.nb Let's look at the first term: In[3]:= td Out[3]= µ- t Cos@0.717 H xld We may calculate the e-folding time for this mode -- call it t1 -- as the reciprocal of the coefficient of t in the exponent: In[37]:= t1 = 1 ê Hlam@1D * DfL Out[37]= This time, which is in seconds, is a little over 50 hours. We compare this with the diffusion time L 2 ëip 2 D f Mwe derived earlier for zero boundary conditions on both sides of the slab. In hours, the time is In[3]:= L 2 ë I300 Ip 2 DfMM Out[3]= This is very much less than the 50 hours we calculated for the first mode in our problem. There are two reasons for this. The first reason is that our slab is insulated on one end, and this will reduce the cooling rate. One can show that the diffusion time for a slab of width L which has zero temperature on one face and which is insulated on the other face is (4 L 2 MëIp 2 D f ). For our present parameters, the value of this in hours is In[39]:= 4 L 2 ë I300 Ip 2 DfMM Out[39]= This is somewhat closer to our calculated value of 50 hours, but still considerably smaller. There is an another effect in the present problem, and that is an additional thermal resistance at the boundary (characterized by the reciprocal of the heat transfer coefficient h). This extra resistance increases the decay time to the 50 hours that we calculated. If you have studied heat transfer, you will know that because the Biot number is around 1, the thermal resistance at the boundary is comparable with that internal to the slab. Let's look at the decay times of the second and third modes, converting them to hours as we calculate them In[40]:= t2 = 1 ê Hlam@2D * Df * 300L Out[40]= In[41]:= t3 = 1 ê Hlam@3D * Df * 300L Out[41]= Knowledge of these decay times will help us choose time values for plotting the solution. To make this a little easier, we construct a table of the first 20 decay times. In[42]:= t@n_d := 1 ê Hlam@nD * DfL
15 sol7.nb 9 In[43]:= TableForm@Table@n, PaddedForm@t@nD,, 2<D, PaddedForm@t@nD ê 300,, 5<D<, n, 1, 20<D, TableHeadings -> None, "Mode", "Decay Time HsL", "Decay Time HhrL"<<D Out[43]//TableForm= Mode Decay Time HsL Decay Time HhrL Graphs of the Solution We define here graphs of the solution for T versus x, for various values of time. We test the code with a few graphs, and then we construct a sequence to be shown in a Manipulate panel. The function which produces the graph at time t is named snapshot[t]. We use all the 40 terms we have calculated in the series to compute the solution, although this is somewhat extravagent in that fewer terms would be sufficient for the longer times. In[44]:= In[45]:= Trans@x_, t_d := Tran@x, t, 40D Temp@x_, t_d := Trans@x, td + Ts@xD Now we are ready to define snapshot[t], which gives us a plot of temperature versus x at time t. From the ambient temperature and the initial distribution, we can get a range of variation of temperature. We extend this range slightly on either side, and call the extended range plrange: In[4]:= plrange = 0, 30<; Now we define snapshot[t], giving a special definition for snapshot[0] -- namely the initial distribution. For convenience, the function snapshot takes an argument in hours, but converts it internally to seconds, to be consistent with the units of the problem. In[47]:= snapshot@t_d := Module@time<, time = 300 * t; Plot@Temp@x, timed, x, 0, L<, PlotStyle Ø Thickness@0.004D, PlotRange -> plrange, AxesLabel -> "x", "Temperature"<, PlotLabel Ø Row@"t = ", PaddedForm@t, 4, 2<D, " hr"<ddd In[4]:= snapshot@0d := Plot@T0@xD, x, 0, L<, PlotStyle Ø Thickness@0.004D, PlotRange -> plrange, AxesLabel -> "x", "Temperature"<, PlotLabel Ø Row@"t = ", PaddedForm@0, 4, 2<D, " hr"<dd We try it for the initial time and for 1,, 20, 30, 40 and 50 hr.
16 sol7.nb In[49]:= Temperature 30 t = 0.00 hr Out[49]= 15 5 x In[50]:= snapshot@1d Temperature 30 t = 1.00 hr Out[50]= 15 5 x In[51]:= snapshot@d Temperature 30 t =.00 hr Out[51]= 15 5 x
17 sol7.nb 11 In[52]:= Temperature 30 t = hr Out[52]= 15 5 x In[53]:= snapshot@30d Temperature 30 t = hr Out[53]= 15 5 x In[54]:= snapshot@40d Temperature 30 t = hr Out[54]= 15 5 x
18 sol7.nb In[55]:= Temperature 30 t = hr Out[55]= 15 5 x As an aid to visualizing the transient process, we produce a sequence of 201 graphs that can be animated to show the dynamics of the heat loss process. The graphs are at 0.5 hr intervals, running from 0 hr to 0 hr. The graphs are output into a Manipulate panel for ease of viewing. In[5]:= DynamicModule@mangraph, i, s<, Do@mangraph@iD = snapshot@0.5 * id, i, 0, 200<D; Manipulate@mangraph@sD, s, 0, 200, 1<DD s Temperature 30 t = 0.00 hr 25 Out[5]= x Have we reached steady-state at 0 hours? Not quite. Here are the endpoint values at 0 hours and in steady-state. In[57]:= Temp@0, 300 * 0D Out[57]= In[5]:= Temp@L, 300 * 0D Out[5]=
19 sol7.nb 13 In[59]:= Out[59]= In[0]:= Out[0]= Let's look at 300 hours. In[1]:= Temperature 30 t = hr Out[1]= 15 5 x In[2]:= Temp@0, 300 * 300D Out[2]= In[3]:= Temp@L, 300 * 300D Out[3]= Now we are off only in the fourth digit.. The One Term Approximation to the Transient The first term in the transient series has a decay time of over 50 hours, whereas the next term has a decay time of only 3.2 hours. Thus the first term will be dominant for times longer than several hours. We now make a graphical comparison of the series solution (in red) and the approximate solution retaining only one term of the transient series (in blue). The basic graphing command is compgraph[t], which takes a time argument in hours and plots both the series solution and the one-term approximate solution. In[4]:= Tapprox@x_, t_d := Ts@xD + firstterm@x, td In[5]:= compgraph@t_d := Module@time<, time = 300 * t; Plot@Temp@x, timed, Tapprox@x, timed<, x, 0, L<, PlotRange Ø 5, 15<, PlotStyle Ø Red, Blue<, FrameLabel -> "", ""<, Frame Ø True, PlotLabel Ø Row@"Exact T HBlueL, Approx T HRedL; ", "t = ", PaddedForm@t, 4, 2<D, " hr"<ddd We try this for 1 hour.
20 sol7.nb In[]:= Exact T HBlueL, Approx T HRedL; t = 1.00 hr Out[]= We look at graphs every half-hour up to 5 hours. In[7]:= Do@Print@compgraph@0.5 * idd, i, 0, <D; Exact T HBlueL, Approx T HRedL; t = 0.00 hr Exact T HBlueL, Approx T HRedL; t = 0.50 hr
21 sol7.nb 15 Exact T HBlueL, Approx T HRedL; t = 1.00 hr Exact T HBlueL, Approx T HRedL; t = 1.50 hr Exact T HBlueL, Approx T HRedL; t = 2.00 hr
22 1 sol7.nb Exact T HBlueL, Approx T HRedL; t = 2.50 hr Exact T HBlueL, Approx T HRedL; t = 3.00 hr Exact T HBlueL, Approx T HRedL; t = 3.50 hr
23 sol7.nb 17 Exact T HBlueL, Approx T HRedL; t = 4.00 hr Exact T HBlueL, Approx T HRedL; t = 4.50 hr Exact T HBlueL, Approx T HRedL; t = 5.00 hr The approximation is good by 5 hours. We also look at hours.
24 1 sol7.nb In[]:= Exact T HBlueL, Approx T HRedL; t =.00 hr Out[]= Now the agreement is excellent.
25 sol7.nb 19 ME201/MTH21/ME400/CHE400 ASSIGNMENT #7 PROBLEM 2 In[9]:= Clear@L, DfD; We use Mathematica to evaluate the integral occurring in the solution for the coefficient C2[t]. In[70]:= Out[70]= C2@t_D = g0 -b2 t IntegrateA b2 x Cos@w xd, x, 0, t<e -t b2 g0 I-b2 + t b2 Hb2 Cos@t wd + w Sin@t wdlm b2 2 + w 2 We now use Mathematica to check our solution, which we must first define. In[71]:= T@x_, t_d = C2@tD SinB 3 2 p x L F; In[72]:= g@x_, t_d = g0 Cos@w td SinB We now substitute this into the equation. 3 2 p x L F; In[73]:= Out[73]= 0 I 3 2 pm2 check = SimplifyBD@T@x, td, td - Df D@T@x, td, x, 2<D - g@x, td ê. b2 Ø Df F The solution satisfies the equation. Now we check the boundary and initial conditions. L 2 In[74]:= T@0, td Out[74]= 0 In[75]:= D@T@x, td, xd ê. x Ø L Out[75]= 0 In[7]:= T@x, 0D Out[7]= 0 Everything checks.
ME201/MTH281/ME400/CHE400 Newton's Law of Cooling
ME1/MTH281/ME0/CHE0 Newton's Law of Cooling 1. Introduction This notebook uses Mathematica to solve the problem presented in class, in section 5.1 of the notes. The problem is one of transient heat conduction
More informationME201/MTH281ME400/CHE400 Convergence of Bessel Expansions
ME201/MTH281ME400/CHE400 Convergence of Bessel Expansions 1. Introduction This notebook provides a general framework for the construction and visualization of partial sums of Fourier-Bessel series. The
More informationME201/MTH281/ME400/CHE400 Zeros of Bessel Function Derivatives
ME201/MTH281/ME400/CHE400 Zeros of Bessel Function Derivatives In[42]:= We write code here to find the n th zero of the derivative of the Bessel function J m. Here n is a positive integer, and m is any
More informationConcepts of Approximate Error: Approximate Error, Absolute Approximate Error, Relative Approximate Error, and Absolute Relative Approximate Error
Concepts of Approximate Error: Approximate Error, Absolute Approximate Error, Relative Approximate Error, and Absolute Relative Approximate Error Sean Rodby, Luke Snyder, Autar Kaw University of South
More informationPhysics 229& 100 Homework #1 Name: Nasser Abbasi
Physics 229& 100 Homework #1 Name: Nasser Abbasi (converted to 6.0) 1. Practice in transcribing expressions into Mathematica syntax a) Find the determinant of this matrix: In[94]:= a = {{1, 2, 3, 4}, {1,
More informationDifferentiation of Continuous Functions Central Difference Approximation of the First Derivative
Differentiation of Continuous Functions Central Difference Approximation of the First Derivative Ana Catalina Torres, Autar Kaw University of South Florida United States of America kaw@eng.usf.edu Introduction
More informationIn[0]:= Dateist@D Out[0]= 82009, 0, 9, 2, 35, 3.457034< ü Bo 0 < < -> 2 Consider the case when we have a particle in the ground state of PIB of length. In[3]:= Y@_, _D := Definitions: 2 SinB p F In[4]:=
More informationThe wave equation. The basic equation and it's fundamental solutions on bounded domains. u t t = c 2 u x x. An amalgam of sections 1.5 and 2.2.
The wave equation An amalgam of sections.5 and.. The basic equation and it's fundamental solutions on bounded domains Suppose that uhx, tl describes the displacement from equlibrium of a vibrating string
More informationME 406 Bifurcations VII Subcritical Hopf Bifurcation
ME 406 Bifurcations VII Subcritical Hopf Bifurcation sysid Mathematica 4.1.2, DynPac 10.66, 3ê5ê2002 intreset; plotreset; 1. Introduction In this notebook, the seventh in a series of notebooks on bifurcations,
More informationMathematica examples relevant to Legendre functions
Mathematica eamples relevant to Legendre functions Legendre Polynomials are built in Here is Legendre s equation, and Mathematica recognizes as being solved by Legendre polynomials (LegendreP) and the
More informationSample Mathematica code for the paper "Modeling a Diving Board"
Sample Mathematica code for the paper "Modeling a Diving Board" Michael A. Karls and Brenda M. Skoczelas, Department of Mathematical Sciences, Ball State University, Muncie, IN 47306 Beam Data Original
More informationME 201/MTH 281/ME400/CHE400 Bessel Functions
ME 01/MTH 81/ME00/CHE00 Bessel Functions 1. Introduction This notebook has two purposes: to give a brief introduction to Bessel functions, and to illustrate how Mathematica can be used in working with
More informationStudents should read Sections of Rogawski's Calculus [1] for a detailed discussion of the material presented in this section.
Chapter 3 Differentiation ü 3.1 The Derivative Students should read Sections 3.1-3.5 of Rogawski's Calculus [1] for a detailed discussion of the material presented in this section. ü 3.1.1 Slope of Tangent
More informationLecture 22 Appendix B (More on Legendre Polynomials and Associated Legendre functions)
Lecture Appendix B (More on Legendre Polynomials and Associated Legendre functions) Ex.7: Let's use Mathematica to solve the equation and then verify orthogonality. First In[]:= Out[]= DSolve y'' x n^y
More informationHandbook of Computational Analytical Heat Conduction
Filippo de Monte, James V. Beck, et al. January 7, 2013 1 Contents 1. Problem description 2 2. Dimensional governing equations. 2 3. Dimensionless variables. 3 4. Dimensionless governing equations.. 3
More informationA Review of Trigonometry
A Review of Trigonometr 0 Kenneth Levasseur Mathematical Sciences UMass Lowell Kenneth_Levasseur@uml.edu Introduction Trigonometr is often introduced as a sstem of ratios of sides of right triangles. Although
More informationME 406 Using Eigenvector Methods with Mathematica to Solve Linear Autonomous Systems of First Order Differential Equations
ME 406 Using Eigenvector Methods with Mathematica to Solve Linear Autonomous Systems of First Order Differential Equations 1. Introduction In this notebook, we use the methods of linear algebra -- specifically
More informationClass 4: More Pendulum results
Class 4: More Pendulum results The pendulum is a mechanical problem with a long and interesting history, from Galileo s first ansatz that the period was independent of the amplitude based on watching priests
More informationMSE 310/ECE 340. Instructor: Bill Knowlton. Initialize data Clear somedata, someplot. Place data in list somedata
MSE 3/ECE 340 Instructor: Bill Knowlton Examples of Plotting Functions ü Example 1: Using the command ListPlot[ ] to plot data in a scatter plot In[1]:= Initialize data Clear somedata someplot Place data
More informationENSC 388. Assignment #8
ENSC 388 Assignment #8 Assignment date: Wednesday Nov. 11, 2009 Due date: Wednesday Nov. 18, 2009 Problem 1 A 3-mm-thick panel of aluminum alloy (k = 177 W/m K, c = 875 J/kg K, and ρ = 2770 kg/m³) is finished
More informationCalculus of Vector-Valued Functions
Chapter 3 Calculus of Vector-Valued Functions Useful Tip: If you are reading the electronic version of this publication formatted as a Mathematica Notebook, then it is possible to view 3-D plots generated
More informationAPPM 1360 Final Exam Fall 2013
APPM 6 Final Exam Fall On the front of our bluebook, please write: a grading ke, our name, student ID, and section and instructor. This exam is worth 5 points and has 9 questions. Show all work! Answers
More informationProject 2 by Kayli-Anna Barriteau
Project 2 by Kayli-Anna Barriteau Decay of quinine in a malaria patient Data for the decay of a mg dose of quinine in a patient suffering from malaria is given in Mathematica format below. The first number
More informationCalculus I Homework: The Derivatives of Polynomials and Exponential Functions Page 1
Calculus I Homework: The Derivatives of Polynomials and Exponential Functions Page 1 Questions Example Differentiate the function y = ae v + b v + c v 2. Example Differentiate the function y = A + B x
More informationDecember 18, 2017 Section W Nandyalam
Mathematica Gift #11 Problem Complete the following 15 exercises in a Mathematica notebook. You should write the problem in a text cell and show your work below the question. Questions Question 1 Compute
More informationPartial Differential Equations Summary
Partial Differential Equations Summary 1. The heat equation Many physical processes are governed by partial differential equations. temperature of a rod. In this chapter, we will examine exactly that.
More informationSolution to Problem Set 5 Physics 480 / Fall 1999 Prof. Klaus Schulten / Prepared by Pinaki Sengupta & Ioan Kosztin
Solution to Problem Set 5 Physics 480 / Fall 1999 Prof. Klaus Schulten / Prepared by Pinaki Sengupta & Ioan Kosztin Problem 1: Energies of the Bound States of the Morse Potential in the Semiclassical Approximation
More informationA small note on the statistical method of moments for fitting a probability model to data
A small note on the statistical method of moments for fitting a probability model to data by Nasser Abbasi, Nov 16, 2007 Mathematics 502 probability and statistics, CSUF, Fall 2007 f x x Μ 2 2 Σ 2 2 Π
More informationPhysicist's Introduction to Mathematica
Physicist's Introduction to Mathematica Physics 200 Physics 50 Lab Revised for V6.0 Fall Semester 2000 Spring Semester 2006 Summer 2007 Nicholas Wheeler John Boccio John Boccio REED COLLEGE Swarthmore
More informationExam Review 2 nd Semester 6-1 Operations on Functions
NAME DATE PERIOD Exam Review 2 nd Semester 6-1 Operations on Functions Find (f + g)(x), (f g)(x), (f g)(x), and (x) for each f(x) and g(x). 1. f(x) = 8x 3; g(x) = 4x + 5 2. f(x) = + x 6; g(x) = x 2 If
More informationMathematica for Calculus II (Version 9.0)
Mathematica for Calculus II (Version 9.0) C. G. Melles Mathematics Department United States Naval Academy December 31, 013 Contents 1. Introduction. Volumes of revolution 3. Solving systems of equations
More informationIntroduction II. accompany our textbook",john H. Mathews, and Russell W.
Introduction II "Essential Mathematica for Students of Science",James J. Kelly, 2006 http://www.physics.umd.edu/courses/courseware "Mathematica 4.1 Notebooks -Complimentary software to accompany our textbook",john
More informationCHAPTER 4. Introduction to the. Heat Conduction Model
A SERIES OF CLASS NOTES FOR 005-006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 4 A COLLECTION OF HANDOUTS ON PARTIAL DIFFERENTIAL EQUATIONS
More informationUsing Mathematica to study series (part II)
Using Mathematica to study series (part II) Truncation errors Example of exponential; first consider x>0 Remainder if use only first 3 terms in power series for exponential: PlotExpx1xx ^,x, 0,, PlotStyleThick,
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationProblem Set 3 Physics 480 / Fall 1999 Professor Klaus Schulten
Problem Set 3 Physics 480 / Fall 1999 Professor Klaus Schulten Problem 1: Wave Packet in 1-Dimensional Box Using Mathematica evaluate the time-dependence of a Gaussian wave packet moving in a 1-dimensional
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationLIMITS AND DERIVATIVES
2 LIMITS AND DERIVATIVES LIMITS AND DERIVATIVES 2.2 The Limit of a Function In this section, we will learn: About limits in general and about numerical and graphical methods for computing them. THE LIMIT
More information3.2. The Ramsey model - linearization, backward integration, relaxation (complete)
3.. The Ramsey model - linearization, backward integration, relaxation (complete) (Quantitative Dynamic Macroeconomics, Lecture Notes, Thomas Steger, University of Leipzig, winter term /3) Ramsey Model:
More informationa) An even function is symmetric with respect to the y-axis. An odd function is symmetric with respect to the origin.
Course Review Course Review Question 1 Page 479 a) An even function is symmetric with respect to the y-axis. An odd function is symmetric with respect to the origin. b) Substitute x for x in f(x). If f(
More informationSection 1.4 Tangents and Velocity
Math 132 Tangents and Velocity Section 1.4 Section 1.4 Tangents and Velocity Tangent Lines A tangent line to a curve is a line that just touches the curve. In terms of a circle, the definition is very
More informationNon-Linear Problems. Almost Linear Systems Consider the system of equations: (1) Y ' = AY + ghyl
Non-Linear Problems Almost Linear Systems Consider the system of equations: (1) Y ' = AY + ghyl (here Y is a matrix). Suppose that the origin is an isolated equilibrium point (that is, there is some circle
More informationThe particular integral :
Second order linear equation with constant coefficients The particular integral : d f df Lf = a + a + af = h( x) Solutions with combinations of driving functions d f df Lf = a + a + af = h( x) + h( x)
More informationExACT Exact Analytical Conduction Toolbox
ExACT Exact Analytical Conduction Toolbox X23B60T0 case Kevin D. Cole and Lukas G. Brettmann April 28, 2014 "Copyright (c) by Kevin D. Cole and Lukas Brettmann" 1 Description Slab body initially at ambient
More informationMECH 375, Heat Transfer Handout #5: Unsteady Conduction
1 MECH 375, Heat Transfer Handout #5: Unsteady Conduction Amir Maleki, Fall 2018 2 T H I S PA P E R P R O P O S E D A C A N C E R T R E AT M E N T T H AT U S E S N A N O PA R T I - C L E S W I T H T U
More informationQuantum Mechanics using Matrix Methods
Quantum Mechanics using Matrix Methods Introduction and the simple harmonic oscillator In this notebook we study some problems in quantum mechanics using matrix methods. We know that we can solve quantum
More informationThe Ramsey model - linearization, backward integration, relaxation (complete)
The Ramsey model - linearization, backward integration, relaxation (complete) (Quantitative Dynamic Macroeconomics, Lecture Notes, Thomas Steger, University of Leipzig) Ramsey Model: linearized version
More informationTRANSIENT HEAT CONDUCTION
TRANSIENT HEAT CONDUCTION Many heat conduction problems encountered in engineering applications involve time as in independent variable. This is transient or Unsteady State Heat Conduction. The goal of
More informationMA S3-MID-Answers Page 1 of 7. Solve the following Wave Equation. Solution by Fourier Series. Let, where is a constant
MA203-3S3-MID-Answers-2050829-Page of 7 Solve the following Wave Equation., 0, 0 0, 0,, 0,, 0 sin,, 0 Solution by Fourier Series Let, So, and, We have,, or or Let 0(see next page for other cases) Then
More informationShowing the shooting method in detail for n = 0
11.2.7 Finding Eigenvalues by the Shooting Method The shooting method is a means to compute eigenvalues of a boundary problem defined by a onedimensional differential equation and boundary conditions fixed
More informationA Poor Example of a Project - But Gives Some Direction
A Poor Example of a Project - But Gives Some Direction Bill Knowlton Materials Science & Engineering Electrical & Computer Engineering Boise State Unversity The idea behind this example is two-fold: 1.
More information2 Discrete Dynamical Systems (DDS)
2 Discrete Dynamical Systems (DDS) 2.1 Basics A Discrete Dynamical System (DDS) models the change (or dynamics) of single or multiple populations or quantities in which the change occurs deterministically
More informationTOPIC 3. Taylor polynomials. Mathematica code. Here is some basic mathematica code for plotting functions.
TOPIC 3 Taylor polynomials Main ideas. Linear approximating functions: Review Approximating polynomials Key formulas: P n (x) =a 0 + a (x x )+ + a n (x x ) n P n (x + x) =a 0 + a ( x)+ + a n ( x) n where
More informationPhysics with Matlab and Mathematica Exercise #1 28 Aug 2012
Physics with Matlab and Mathematica Exercise #1 28 Aug 2012 You can work this exercise in either matlab or mathematica. Your choice. A simple harmonic oscillator is constructed from a mass m and a spring
More informationOrdinary Differential Equations
Instructor: José Antonio Agapito Ruiz UCSC Summer Session II, 00 Math Ordinar Differential Equations Final Solution This report does not contain full answers to all questions of the final. Instead, ou
More informationINTRO TO LIMITS & CALCULUS MR. VELAZQUEZ AP CALCULUS
INTRO TO LIMITS & CALCULUS MR. VELAZQUEZ AP CALCULUS WHAT IS CALCULUS? Simply put, Calculus is the mathematics of change. Since all things change often and in many ways, we can expect to understand a wide
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationProfessor Faith Morrison
CM3215 Fundamentals of Chemical Engineering Laboratory Professor Faith Morrison Department of Chemical Engineering Michigan Technological University 1 Objective: Measure the heat-transfer coefficient to
More informationChapter 11 Parametric Equations, Polar Curves, and Conic Sections
Chapter 11 Parametric Equations, Polar Curves, and Conic Sections ü 11.1 Parametric Equations Students should read Sections 11.1-11. of Rogawski's Calculus [1] for a detailed discussion of the material
More informationSolution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)
Math 380 Practice Final Solutions This is longer than the actual exam, which will be 8 to 0 questions (some might be multiple choice). You are allowed up to two sheets of notes (both sides) and a calculator,
More informationMathematica Project, Math21b Spring 2008
Mathematica Project, Math21b Spring 2008 Oliver Knill, Harvard University, Spring 2008 Support Welcome to the Mathematica computer project! In case problems with this assignment, please email Oliver at
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationDepending on your Mathematica installation, you might see the solution curve corresponding to q
Problem 5 (a) Clear@pendD; pend@t_, x0_, y0_d := x@td, y@td
More informationNumerical Heat and Mass Transfer
Master degree in Mechanical Engineering Numerical Heat and Mass Transfer 02-Transient Conduction Fausto Arpino f.arpino@unicas.it Outline Introduction Conduction ü Heat conduction equation ü Boundary conditions
More information3.3 Unsteady State Heat Conduction
3.3 Unsteady State Heat Conduction For many applications, it is necessary to consider the variation of temperature with time. In this case, the energy equation for classical heat conduction, eq. (3.8),
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More informationME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included
ME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included sysid Mathematica 6.0.3, DynPac 11.01, 1ê13ê9 1. Introduction Description of the Model In this notebook, we include births and deaths in the
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechanics for Scientists and Engineers David Miller Particles in potential wells The finite potential well Insert video here (split screen) Finite potential well Lesson 7 Particles in potential
More informationFactors of Polynomials Factoring For Experts
Factors of Polynomials SUGGESTED LEARNING STRATEGIES: Shared Reading, Activating Prior Knowledge, Discussion Group, Note-taking When you factor a polynomial, you rewrite the original polynomial as a product
More informationRelativistic Motion of a Charged Particle and Pauli Algebra Jan Vrbik
The Mathematica Journal Relativistic Motion of a Charged Particle and Pauli Algebra Jan Vrbik We introduce some key formulas of special relativity and apply them to the motion of a spinless, charged point
More informationStochastic Integrals and Their Expectations
The Mathematica Journal Stochastic Integrals and Their Expectations Wilfrid S. Kendall This article explains how the Itovsn3 package can be extended to add various properties and rules for ItoIntegral,
More informationHW #6. 1. Inflaton. (a) slow-roll regime. HW6.nb 1
HW6.nb HW #6. Inflaton (a) slow-roll regime In the slow-roll regime, we neglect the kinetic energy as well as f ÿÿ term in the equation of motion. Then H = ÅÅÅ 8 p 3 G N ÅÅÅ m f, 3 H f ÿ + m f = 0. We
More informationLogarithmic Functions
Name Student ID Number Group Name Group Members Logarithmic Functions 1. Solve the equations below. xx = xx = 5. Were you able solve both equations above? If so, was one of the equations easier to solve
More informationIntroduction to Techniques for Counting
Introduction to Techniques for Counting A generating function is a device somewhat similar to a bag. Instead of carrying many little objects detachedly, which could be embarrassing, we put them all in
More informationSO x is a cubed root of t
7.6nth Roots 1) What do we know about x because of the following equation x 3 = t? All in one.docx SO x is a cubed root of t 2) Definition of nth root: 3) Study example 1 4) Now try the following problem
More informationGlobal Warming? Effects on Glaciers
This project has been funded with support from the European Commission. This publication reflects the views only of the author, and the Commission cannot be held responsible for any use which may be made
More informationMidterm Review. Name: Class: Date: ID: A. Short Answer. 1. For each graph, write the equation of a radical function of the form y = a b(x h) + k.
Name: Class: Date: ID: A Midterm Review Short Answer 1. For each graph, write the equation of a radical function of the form y = a b(x h) + k. a) b) c) 2. Determine the domain and range of each function.
More informationMathematica and the Glass Tempering Process
Mathematica and the Glass Tempering Process Jennifer Voitle Introduction Mathematica is being used in industry to solve difficult and highly complex problems, formerly insoluble or requiring the power
More informationPhysics 212E Spring 2004 Classical and Modern Physics. Computer Exercise #2
Physics 212E Spring 2004 Classical and Modern Physics Chowdary Computer Exercise #2 Launch Mathematica by clicking on the Start menu (lower left hand corner of the screen); from there go up to Science
More informationhp calculators HP 33S Using the formula solver part 2 Overview of the formula solver Practice Example: A formula with several variables
Overview of the formula solver Practice Example: A formula with several variables Practice Example: A direct solution Practice Example: A formula with a pole Practice Example: Where two functions intersect
More informationHomework 7 Solutions
Homework 7 Solutions # (Section.4: The following functions are defined on an interval of length. Sketch the even and odd etensions of each function over the interval [, ]. (a f( =, f ( Even etension of
More informationMATH 1231 MATHEMATICS 1B Calculus Section 2: - ODEs.
MATH 1231 MATHEMATICS 1B 2007. For use in Dr Chris Tisdell s lectures: Tues 11 + Thur 10 in KBT Calculus Section 2: - ODEs. 1. Motivation 2. What you should already know 3. Types and orders of ODEs 4.
More informationQuantum Wells The Shooting Method (application to the simple harmonic oscillator and other problems of a particle in a potential)
quantumwellsho.nb Quantum Wells The Shooting Method (application to the simple harmonic oscillator and other problems of a particle in a potential) Introduction In the previous handout we found the eigenvalues
More information() Chapter 8 November 9, / 1
Example 1: An easy area problem Find the area of the region in the xy-plane bounded above by the graph of f(x) = 2, below by the x-axis, on the left by the line x = 1 and on the right by the line x = 5.
More informationDifferentiation 1. The METRIC Project, Imperial College. Imperial College of Science Technology and Medicine, 1996.
Differentiation 1 The METRIC Project, Imperial College. Imperial College of Science Technology and Medicine, 1996. 1 Launch Mathematica. Type
More informationPDE and Boundary-Value Problems Winter Term 2014/2015
PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 6 Saarland University 17. November 2014 c Daria Apushkinskaya (UdS) PDE and BVP lecture 6 17. November 2014 1 / 40 Purpose of Lesson To show
More informationf[x_, μ_] := 4. μ... Nest[f[#,...] &,...,...] Plot[{x, f[...]}, {x, 0, 1}, AspectRatio Equal]
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 In this exercise, we use renormalization-group and scaling methods to study the onset of chaos. There are several routes by which a dynamical system can start exhibiting chaotic
More informationfirst define colors for future graphics red RGBColor 1, 0, 0 ; green RGBColor 0, 1, 0 ; blue RGBColor 0, 0, 1 ; black RGBColor 0, 0, 0 ;
In[1]:= note: to simplifynotations all periodic functions of x are expected to have a period of 2 Pi, as Ch. 11.1 by KR. If a function has a period "2L" as in his later sections, one can always switch
More informationx 4 D: (4, ); g( f (x)) = 1
Honors Math 4 Describing Functions One Giant Review Name Answer Key 1. Let f (x) = x, g(x) = 6x 3, h(x) = x 3 a. f (g(h(x))) = 2x 3 b. h( f (g(x))) = 1 3 6x 3 c. f ( f ( f (x))) = x 1 8 2. Let f (x) =
More information2 = = 0 Thus, the number which is largest in magnitude is equal to the number which is smallest in magnitude.
Limits at Infinity Two additional topics of interest with its are its as x ± and its where f(x) ±. Before we can properly discuss the notion of infinite its, we will need to begin with a discussion on
More information10.7 Trigonometric Equations and Inequalities
0.7 Trigonometric Equations and Inequalities 79 0.7 Trigonometric Equations and Inequalities In Sections 0., 0. and most recently 0., we solved some basic equations involving the trigonometric functions.
More informationThe connected locus for complex cubic iteration
The connected locus for complex cubic iteration A preprint version of a Mathematical graphics column from Mathematica in Education and Research. Mark McClure Department of Mathematics University of North
More informationMATH 114 Calculus Notes on Chapter 2 (Limits) (pages 60-? in Stewart)
Still under construction. MATH 114 Calculus Notes on Chapter 2 (Limits) (pages 60-? in Stewart) As seen in A Preview of Calculus, the concept of it underlies the various branches of calculus. Hence we
More informationProblem Set 1 Solution
Problem Set 1 Solution Friday, 9 September 2011 Physics 111 Exercise 1 Watch the 20-minute Wolfram tutorial introducing Mathematica at http://url. wolfram.com/fgkkpc/. Even if you already know how to use
More informationChapter 4: Transient Heat Conduction. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University
Chapter 4: Transient Heat Conduction Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Objectives When you finish studying this chapter, you should be able to: Assess when the spatial
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES
More informationApplications of integrals
ApplicationsofIntegrals.nb Applications of integrals There are many applications of definite integrals and we cannot include all of them in a document intended to be a review. However there are some very
More informationME 406 Assignment 6 Solutions
ME 406 Assignment 6 Solutions sysid Mathematica 6.0.3, DynPac 11.02, 3ê12ê2009 inreset; plotreset; Problem 1 ü Part a We look for equilibria. The first equation requires y = 0, and the second equation
More informationMath 41 Final Exam December 9, 2013
Math 41 Final Exam December 9, 2013 Name: SUID#: Circle your section: Valentin Buciumas Jafar Jafarov Jesse Madnick Alexandra Musat Amy Pang 02 (1:15-2:05pm) 08 (10-10:50am) 03 (11-11:50am) 06 (9-9:50am)
More information