How to get partial surface integration ? Log Out | Topics | Search
Moderators | Register | Edit Profile

FlexPDE User's Forum » User Postings » How to get partial surface integration ? « Previous Next »

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

Liem Peng Hong (liemph)
Member
Username: liemph

Post Number: 8
Registered: 09-2004
Posted on Sunday, June 25, 2006 - 09:43 pm:   

Dea Forum,

I'd like to evaluate the leakage at the boundary of a 2D (RZ geometry) diffusion problem. I can get the total leakage by the following command:

report (surf_integral(normal(Fluxd), "Outer") / vol_integral(Production, "Reactor")) as "Lkg"

However, how can I get the leakage breakdown, i.e. only from the top, bottom and lateral surface ? I tried to divide the "Outer" surface into 3 surfaces (top, bottom and lateral surfaces of the cylinder) but failed.

Thanks in advance.

Liem.
Top of pagePrevious messageNext messageBottom of page Link to this message

Robert G. Nelson (rgnelson)
Moderator
Username: rgnelson

Post Number: 651
Registered: 06-2003
Posted on Tuesday, July 11, 2006 - 12:58 pm:   

Overlay named features on the segments of the boundary you want to integrate, and then specify those boundary names in the surf_integral() call.
Top of pagePrevious messageNext messageBottom of page Link to this message

Liem Peng Hong (liemph)
Member
Username: liemph

Post Number: 9
Registered: 09-2004
Posted on Thursday, July 20, 2006 - 10:20 am:   

Dear Robert,
Thank you for your reply, however still I can not understand the meaning of overlay. Here is my domain definition (R-Z geometry, right cylinder).

BOUNDARIES
REGION 1 "Reactor"
start "Outer"
( 0., Height/2) Natural (Flux) = -Flux/(0.7104*Lamda_Tr)
line to ( Rad, Height/2) Natural (Flux) = -Flux/(0.7104*Lamda_Tr)
line to ( Rad, 0.) Natural (Flux) = 0.0
line to ( 0., 0.) Natural (Flux) = 0.0
line to close

If I want to evaluate the surf_integral() for just the segment of (0, Height/2) -> (Rad, Height/2), what should I do ?

Thanks in advance,

Liem
Top of pagePrevious messageNext messageBottom of page Link to this message

Robert G. Nelson (rgnelson)
Moderator
Username: rgnelson

Post Number: 661
Registered: 06-2003
Posted on Thursday, July 20, 2006 - 04:07 pm:   

Add to BOUNDARIES section:
FEATURE START "IntegralPath" (0,Height/2) LINE TO (Rad,Height/2)

Add wherever you want to use it:
... surf_integral(<expression>,"IntegralPath")

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