Hi All

I am a stranger to mathematica. I would like to calculate the average intensity in a UV reactor. (The reactor is a cylinder with a UV lamp placed parallel to the height of the cylinder.)

cylinder dimensions(say)base=40 mm, height=600mm;

length of UV lamp=450mm).

The Intensity at a point due to the UV lamp is given by:I= P/(2*pi*r*L).

where r=distance of the pt. from the centre of the lamp;

pi=3.14;

L= lamp length

P=power of the lamp.

If i want to calculate the Ave. Intensity in the reactor, how to program this in Mathematica.

Ave. Intensity I_ave=Summation(I) over all the points of the volume.

Please give me a hint how to write a program for this in Mathematica.

Here i am giving my procedure, which doesn't give desired result

x = 0; y = 0;(x,y) is the coordinates of the lamp centre.

P = 10.3; Power of the lamp .

l = 0.45; Length of the lamp.

r0 = 0.025; Radius of the lamp.

r[xk_, yk_] := Sqrt[(x - xk)^2 + (y - yk)^2]; distances from the centre of the lamp to Kth point in the domain.

Intensity[xk_, yk_] := (P/(2*\[Pi]*(Sqrt[(x - xk)^2 + (y - yk)^2])*l));

Integrate[Intensity[xk, yk], {xk, -0.1, 0.1}, {yk, 0, 0.2}]

Actually i need to calculate the ave. intensity in a volume V of a Reactor. For this i have to find the Intensity at a point P(Xk,Yk), which is at a dist. of 'r' from the centre of the lamp ie (0,0).

Please anyone tell me how to do this and how to plot the Intensity field.

Thanks in Advance

Link to the forum page for this post:

http://www.yqcomputer.com/ :Forum_ViewTopic&pid=12752#p12752

Posted through http://www.yqcomputer.com/ [[postId=12752]]

Hi,

you want to calculate the mean value of a function: fun in a cylindrical

region. You can do this by calculating the integral of the function over

the region and divide by the volume. In your case this is best done ijn

cylindrical coordinates. Here, the volume element is:

r dphi dr dz

where r is the distance from the axis, z the heigth and phi the angle.

As the function fun only depends on r an z the angular integral can be

done at once, giving 2 Pi. This leaves:

2 Pi NIntegrate[r fun[r,z],{r,rmin,rmax},{z,0,zmax}] /

((rmax^2-rmin^2)Pi zmaz)

Daniel

1. Beginner: Phong Shading - what is my color from intensity

2. Beginner: Shading models- what is my color from intensity

3. advice to a beginner's beginner?

4. a beginner, beginners question

5. No Class at ALL!!! beginner/beginner question

6. Help for a Beginner Among Beginners

8. Beginner of Beginner getting an error

9. Tutorial for beginner/ Tutorial voor beginner

10. beginner - calculate totals per day

11. Beginner: Calculating with fields

12. beginners problem with decimal point display

13. beginners question: awarding points to a certain amount bought

14. beginner problem - point me in right direction?

15. Beginner question: Precision of floating point arithmetic...

2 post • Page:**1** of **1**