My Math Forum  

Go Back   My Math Forum > College Math Forum > Applied Math

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}$

$\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)?
Benit13 is offline  
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 series-based 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_1-x_0) + \frac{1}{2} f''(x_0)(x_1-x_0)^2 + O((x_1-x_0)^3)$
$\displaystyle f(x_2) = f(x_0) + f'(x_0)(x_2-x_0) + \frac{1}{2} f''(x_0)(x_2-x_0)^2 + O((x_2-x_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_2-x_1} + O((x_2-x_1)^3), ~if~ x_2-x_0=x_0-x_1$ (central difference), or $\displaystyle +O((x_2-x_1)^2)$ otherwise.
Definite integrals work similarly, though I'd have to sit down with one for a few minutes to figure it out (Runge-Kutta methods are developed this way, though).

Chebyshev approximations, if I'm not mistaken, approximate the function not as a superposition of Nth-order 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.
DarnItJimImAnEngineer is offline  
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 user-input 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(2j-1)}{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.
Attached Images
File Type: jpg table.jpg (15.1 KB, 5 views)

Last edited by Benit13; September 6th, 2019 at 02:02 AM.
Benit13 is offline  

  My Math Forum > College Math Forum > Applied Math

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 aaron-math Calculus 7 September 22nd, 2011 11:25 AM
Chebyshev's theorem finitehelp Algebra 6 July 31st, 2009 02:50 PM
Understanding Integral Calculus W-Unit Calculus 8 May 7th, 2009 10:18 PM
Chebyshev's theorem sangfroid Advanced Statistics 0 January 17th, 2008 11:38 AM

Copyright © 2019 My Math Forum. All rights reserved.