2D temperature distribution Log Out | Topics | Search
Moderators | Register | Edit Profile

FlexPDE User's Forum » User Postings » 2D temperature distribution « Previous Next »

Author Message
Top of pagePrevious messageNext messageBottom of page Link to this message

abdulwhab (tolana)
New member
Username: tolana

Post Number: 1
Registered: 03-2009
Posted on Sunday, March 22, 2009 - 09:41 pm:   

hope that someone will help me with this, here is the problem statment.

Use FlexPDE to solve for the temperature distribution in a wall of thickness b
extending from 0 to 10b in the y direction, and from - to + in the direction
perpendicular to the x and y directions. The surfaces at x= +/- b/2 and the top surface at
y=10b are held at temperature To, whereas the bottom surface of the wall at y=0 is
maintained at temperature T1.

and here what I did.
application/octet-stream
My_1.pde (1.0 k)
Top of pagePrevious messageNext messageBottom of page Link to this message

Robert G. Nelson (rgnelson)
Moderator
Username: rgnelson

Post Number: 1228
Registered: 06-2003
Posted on Monday, March 23, 2009 - 07:02 pm:   

You haven't put any boundary conditions on the problem. You need to use the VALUE(temp)=T0 and VALUE(temp)=T1 commands to place the boundary conditions.

However, adjacent boundaries with different VALUE conditions are a physical impossibility. It implies an infinite gradient, therefore an infinite flux at the corner point. But go ahead and do it. Just realize that the corner is unphysical.
Top of pagePrevious messageNext messageBottom of page Link to this message

abdulwhab (tolana)
New member
Username: tolana

Post Number: 2
Registered: 03-2009
Posted on Monday, March 23, 2009 - 11:25 pm:   

Thaks for the response, it was very helpfull. I have attached the file to make sure that I got it write.
application/octet-stream
My_1.pde (1.0 k)
Top of pagePrevious messageNext messageBottom of page Link to this message

abdulwhab (tolana)
Junior Member
Username: tolana

Post Number: 3
Registered: 03-2009
Posted on Monday, March 23, 2009 - 11:28 pm:   

I also have assumed my own geometry and solved for it, I am just wondering if this goemetery applicable or not.
application/octet-stream
2.pde (1.0 k)
Top of pagePrevious messageNext messageBottom of page Link to this message

Marek Nelson (mgnelson)
Moderator
Username: mgnelson

Post Number: 116
Registered: 07-2007
Posted on Tuesday, March 24, 2009 - 09:14 pm:   

Your script "my_1" has the boundary conditions that you wanted.

The script "2" has a problem with the region definition. You have made a "bow-tie" not a rectangle. The FlexPDE gridder can handle this (you get two triangles with their tips touching), but now you have two boundaries crossing, one with a VALUE(temp)=70 the other is -15. Check your region definition. Your last point is probably in error.
Top of pagePrevious messageNext messageBottom of page Link to this message

abdulwhab (tolana)
Member
Username: tolana

Post Number: 4
Registered: 03-2009
Posted on Wednesday, March 25, 2009 - 02:52 am:   

Thank you very much for being generous in answering my question.

One last question, if the hight of the wall is infinity, then would FlexPDE handle the problem and if not what are my alternatives.

Note: I tride to change the boundry condition for y to be a huge number, but it did not work it out.
Top of pagePrevious messageNext messageBottom of page Link to this message

Robert G. Nelson (rgnelson)
Moderator
Username: rgnelson

Post Number: 1232
Registered: 06-2003
Posted on Thursday, March 26, 2009 - 06:35 pm:   

1.
Computers cannot compute with infinities, and it would be hard to draw the picture on your screen.

2.
Use a height that is great enough that a variation in the height makes negligible effect on the solution at the bottom.

3.
Use a less huge number.

Add Your Message Here
Post:
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Password:
Options: Enable HTML code in message
Automatically activate URLs in message
Action:

Topics | Last Day | Last Week | Tree View | Search | Help/Instructions | Program Credits Administration