
Applied Math Applied Math Forum 
 LinkBack  Thread Tools  Display Modes 
September 3rd, 2019, 04:54 AM  #1 
Senior Member Joined: Apr 2014 From: Glasgow Posts: 2,161 Thanks: 734 Math Focus: Physics, mathematical modelling, numerical and computational solutions  Understanding a specific Chebyshev integral
I have a situation where I'm trying to understand the physics of a cooling tower. Part of the solution for the model I'm working with involves an integral: $\displaystyle C = \int_{T_{low}}^{T_{high}} \frac{1}{\left(h'  h_a\right)} dT$ where $\displaystyle h'$ is a slowly varying, monotonically increasing function with temperature, T, and $\displaystyle h_a$ is linearly increasing with a gradient $\displaystyle m$. Both functions for $\displaystyle h'$ and $\displaystyle h_a$ can be evaluated for any known T and $\displaystyle h' > h_a$ for all $\displaystyle T_{low} \le T \le T_{high}$. Now... according to a particular reference (I don't have a copy, but it's known as British Standard 4485... catchy!) the solution to this integral involves performing a Chebyshev approximation and then performing a rectangular integration on the pieces: $\displaystyle C = \frac{\Delta T}{4}\sum_{i=1}^{4} \frac{1}{\left(h'(T_i)  h_a(T_i)\right)}$ $\displaystyle \Delta T = T_{high}  T_{low}$ where $\displaystyle T_1 = T_{high}  0.1026728 \Delta T$ $\displaystyle T_2 = T_{high}  0.4062038 \Delta T$ $\displaystyle T_3 = T_{low} + 0.4062038 \Delta T$ $\displaystyle T_4 = T_{low} + 0.1026728 \Delta T$ I don't understand this approximation at all. Could someone please explain this to me (or point out some links or references that explain this outcome)? 
September 5th, 2019, 10:41 AM  #2 
Senior Member Joined: Jun 2019 From: USA Posts: 213 Thanks: 90 
Since no one's answered, let me say that I can't tell you specifically where those numbers came from, but I can maybe give you a general idea. If these were Taylor seriesbased integrals, I could probably derive it for you on the spot. With those, you are using Taylor series expansion to approximate derivatives (or integrals) of the function using only a finite footprint (such as the four points you have here). For example: $\displaystyle f(x_1) = f(x_0) + f'(x_0)(x_1x_0) + \frac{1}{2} f''(x_0)(x_1x_0)^2 + O((x_1x_0)^3)$ $\displaystyle f(x_2) = f(x_0) + f'(x_0)(x_2x_0) + \frac{1}{2} f''(x_0)(x_2x_0)^2 + O((x_2x_0)^3)$ Then by choosing a certain number of points and/or a certain desired accuracy (leading order of the dropped terms), you can combine these expansions to approximate any of the derivatives $\displaystyle f^{(n)}(x_0)$ as a weighted average of neighbouring points $\displaystyle c_1 f(x_1) + c_2 f(x_2) + ...$, e.g., $\displaystyle f'(x_0) = \frac{f(x_2)f(x_1)}{x_2x_1} + O((x_2x_1)^3), ~if~ x_2x_0=x_0x_1$ (central difference), or $\displaystyle +O((x_2x_1)^2)$ otherwise. Definite integrals work similarly, though I'd have to sit down with one for a few minutes to figure it out (RungeKutta methods are developed this way, though). Chebyshev approximations, if I'm not mistaken, approximate the function not as a superposition of Nthorder polynomials, but as polynomials in sines and cosines. You'd have to go googling for derivations to find the exact source of those coefficients, though. Out of curiosity, are you sure British Standard 4485 is still the industrial standard? (Not that it would surprise me terribly.) Assuming your Rayleigh/Prandtl numbers aren't changing wildly, it's not expensive to calculate the convection coefficients for dozens or even hundreds of temperatures, and get a much better approximation of the integral than just relying on four points. 
September 6th, 2019, 01:56 AM  #3 
Senior Member Joined: Apr 2014 From: Glasgow Posts: 2,161 Thanks: 734 Math Focus: Physics, mathematical modelling, numerical and computational solutions 
Thanks for your reply. I've been researching Chebyshev approximation techniques and integrals now for several days and I wrote a Python script that manages to perform function approximations for any userinput function and solve the integral between two arbitrary limits. It's been a good exercise. Numerical recipes sort of came to the rescue in this regard. I'm not certain whether BS4485 is *the* current standard used for cooling tower performance calculations, but it seems vestiges of it are still sitting in old codes for modelling cooling towers. It's adoption is probably due to the simplicity of the method for hand calculations in the 80s. The closest I've found to finding out the origin of those particular numbers are these two characteristics: 1. The roots of Chebyshev polynomials of the first kind are: $\displaystyle x_j = \cos\left(\frac{\pi(2j1)}{2N}\right)$ where N is the order of the the polynomial and j = 1, ..., N. These are the points which should be used to evaluate integrals and consequently those numbers should probably be associated with cosines of some angle. However, it seems no combination of j or N can give the values specified in the table. 2. I found a table in a book on thermal heat transfer which provides weights and ordinate directions for a discretized ordinate method for radiative heat transfer (see attached file). The methods described in the book talk about approximation methods like Chebyshev approximation, but applied to radiative transfer problems. Those S4 numbers are the same ones! Unfortunately, I can't find any references that explain where those numbers come from and how they are calculated; the original reference in the table is hidden behind a paywall. Last edited by Benit13; September 6th, 2019 at 02:02 AM. 

Tags 
chebyshev, integral, specific, understanding 
Thread Tools  
Display Modes  

Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Chebyshev's inequality  kameleon_99  Advanced Statistics  3  January 24th, 2014 05:50 AM 
Chebyshev polynomials  aaronmath  Calculus  7  September 22nd, 2011 11:25 AM 
Chebyshev's theorem  finitehelp  Algebra  6  July 31st, 2009 02:50 PM 
Understanding Integral Calculus  WUnit  Calculus  8  May 7th, 2009 10:18 PM 
Chebyshev's theorem  sangfroid  Advanced Statistics  0  January 17th, 2008 11:38 AM 