Concept explainers
To calculate: The temperature on the heated plate with fixed temperature at the boundaries and round corners at upper right and lower left corners; considering the value of weighing factor for overrelaxation as 1.5 and performing iteration up to
Answer to Problem 7P
Solution:
Heated plate having rounded upper right and lower left corners has the temperature distribution with the desired level of accuracy as,
Explanation of Solution
Given Information:
The fixed temperature at top and right boundary is
For overrelaxation the weighing factor is 1.5.
Limiting condition for iteration is
Formula used:
For a square grid, that is,
With
Overrelaxation is applied using the following expression after each iteration,
In the above expression,
The percent relative error is calculated using the following expression,
Laplacian equation for a node adjacent to an irregular boundary is solved for temperature distribution as,
Calculation:
Consider,
Consider the value of the parameters as,
The temperature at the node
Thus,
Further, the temperature at the node
Thus,
Liebmann’s method is used to solve the temperature on the heated plate iteratively in MATLAB, until the percent relative error falls below 1%.
Use the following MATLAB code to implement the Liebmann’s method and the equation at the nodes
Execute the above code to obtain the temperature distribution and the number of iterations performed, in the output as shown below,
Hence, the temperature on the heated plate with desired accuracy is achieved after 9 iterations as obtained in the above MATLAB output.
Want to see more full solutions like this?
Chapter 29 Solutions
EBK NUMERICAL METHODS FOR ENGINEERS
- 1.6. By manipulating Taylor series, determine the constant C for an error expansion of (1.3) of the form wj−u' (xj) ~ Ch¼u (5) (x;), where u (5) denotes the fifth derivative. Based on this value of C and on the formula for u(5) (x) with u(x) = esin(x), determine the leading term in the expansion for w; - u'(x;) for u(x) = esin(x). (You will have to find maxε[-T,T] |u(5) (x)| numerically.) Modify Program 1 so that it plots the dashed line corresponding to this leading term rather than just N-4. This adjusted dashed line should fit the data almost perfectly. Plot the difference between the two on a log-log scale and verify that it shrinks at the rate O(h6).arrow_forwardDefine sinc(x) = sin(x)/x, except with the singularity removed. Differentiate sinc(x) once and twice.arrow_forward1.4. Run Program 1 to N = 216 instead of 212. What happens to the plot of error vs. N? Why? Use the MATLAB commands tic and toc to generate a plot of approximately how the computation time depends on N. Is the dependence linear, quadratic, or cubic?arrow_forward
- Show that the function f(x) = sin(x)/x has a removable singularity. What are the left and right handed limits?arrow_forward18.9. Let denote the boundary of the rectangle whose vertices are -2-2i, 2-21, 2+i and -2+i in the positive direction. Evaluate each of the following integrals: (a). 之一 dz, (b). dz, (b). COS 2 coz dz, dz (z+1) (d). z 2 +2 dz, (e). (c). (2z+1)zdz, z+ 1 (f). £, · [e² sin = + (2² + 3)²] dz. (2+3)2arrow_forward18.10. Let f be analytic inside and on the unit circle 7. Show that, for 0<|z|< 1, f(E) f(E) 2πif(z) = --- d.arrow_forward
- 18.4. Let f be analytic within and on a positively oriented closed contoury, and the point zo is not on y. Show that L f(z) (-20)2 dz = '(2) dz. 2-20arrow_forward18.9. Let denote the boundary of the rectangle whose vertices are -2-2i, 2-21,2+i and -2+i in the positive direction. Evaluate each of the following integrals: (a). rdz, (b). dz (b). COS 2 coz dz, (z+1) (d). 之一 z 2 +2 dz, (e). dz (c). (2z + 1)2dz, (2z+1) 1 (f). £, · [e² sin = + (2² + 3)²] dz. z (22+3)2arrow_forward18.8. (a). Let be the contour z = e-≤0≤ traversed in the า -dz = 2xi. positive direction. Show that, for any real constant a, Lex dzarrow_forward
- f(z) 18.7. Let f(z) = (e² + e³)/2. Evaluate dz, where y is any simple closed curve enclosing 0.arrow_forward18. If m n compute the gcd (a² + 1, a² + 1) in terms of a. [Hint: Let A„ = a² + 1 and show that A„|(Am - 2) if m > n.]arrow_forwardFor each real-valued nonprincipal character x mod k, let A(n) = x(d) and F(x) = Σ : dn * Prove that F(x) = L(1,x) log x + O(1). narrow_forwardarrow_back_iosSEE MORE QUESTIONSarrow_forward_ios
- Advanced Engineering MathematicsAdvanced MathISBN:9780470458365Author:Erwin KreyszigPublisher:Wiley, John & Sons, IncorporatedNumerical Methods for EngineersAdvanced MathISBN:9780073397924Author:Steven C. Chapra Dr., Raymond P. CanalePublisher:McGraw-Hill EducationIntroductory Mathematics for Engineering Applicat...Advanced MathISBN:9781118141809Author:Nathan KlingbeilPublisher:WILEY
- Mathematics For Machine TechnologyAdvanced MathISBN:9781337798310Author:Peterson, John.Publisher:Cengage Learning,