-
Filter by Topic
Most Popular
All Topics
- List all discussions
Facebook
Twitter
Delicious
Digg
Stumbleupon
More services...Flux in Convection and diffusion module : A quick question
|
Thread index | Previous thread | Next thread | Start a new discussion |
RSS feed |
Turn on email notifications |
8 Replies
Last post: April 8, 2010 9:58pm UTC
March 19, 2010 9:27am UTC
Flux in Convection and diffusion module : A quick question
Hi,
Is there such flux Lagrange multiplier in this mode, as there is for drag/lift or reacf(u), or I can use nflux_total accurately?
Many thanks
Reply | Reply with Quote | Send private message | Report Abuse
March 20, 2010 10:01am UTC in response to Danial Tz
Re: Flux in Convection and diffusion module : A quick question
update: the fluxes calculated using ntotalflux are wrong, specially when there is a moving boundary. this is frustrating.
Reply | Reply with Quote | Send private message | Report Abuse
March 21, 2010 10:53am UTC in response to Danial Tz
Re: Flux in Convection and diffusion module : A quick question
Hi
Is it possible to be slightly more explicit about your model ?
such as version, application mode name, minimum model info, BC, ...
It would make life slightly easer for us to understand and perhaps come with a suggestion, or just learn something ourselves.
Have fun with COMSOL
Ivar
Reply | Reply with Quote | Send private message | Report Abuse
March 21, 2010 5:47pm UTC in response to Ivar Kjelberg
Re: Flux in Convection and diffusion module : A quick question
Hi Ivar,
My question is simple and general.
How can I calculate the mass flux going inside a subdomain accurately?
when I activate non-ideal weak mode for "convection and diffusion", say lm3, I can not integrate it, like I do with "drag coefficient" in "cylinder in the flow" case.
I am using comsol 3.5a, and modeling an oscillating cylinder (ALE) in the flow with reaction inside the cylinder, as a benchmark test. It is just the comsol "cylinder in the flow", but oscillating one, and reaction happening inside.
Many thanks in advance
Reply | Reply with Quote | Send private message | Report Abuse
March 22, 2010 8:33am UTC in response to Danial Tz
Re: Flux in Convection and diffusion module : A quick question
Hi
If I record correctly, lagrange multipliers "lm"s are defined on the nodes and should be "summed" and not integrated, have you tried the posprocessing BC "advanced" feature ? I cannot remember now if they are "automatically" summed rather than integrated.
Good luck
Ivar
Reply | Reply with Quote | Send private message | Report Abuse
March 22, 2010 5:27pm UTC in response to Ivar Kjelberg
Re: Flux in Convection and diffusion module : A quick question
Hi
If I record correctly, lagrange multipliers "lm"s are defined on the nodes and should be "summed" and not integrated, have you tried the posprocessing BC "advanced" feature ? I cannot remember now if they are "automatically" summed rather than integrated.
Good luck
Ivar
Dear Ivar,
They are automatically summed. correct. However, Comsol doesn't let me do any summation/integration on the boundaries for lm of convection & diffusion. So, still the question is how to find flux?
Reply | Reply with Quote | Send private message | Report Abuse
March 22, 2010 5:33pm UTC in response to Danial Tz
Re: Flux in Convection and diffusion module : A quick question
Here is an idea that might help under some special circumstances:
You can use the divergence theorem and re-express your boundary flux as a domain integral which you can calculate easily. Catch is, this way you cannot figure out individual boundary contributions; what you get is the total net flux across all boundaries enclosing your domain.
Even in cases when boundary integrals are calculable, domain integral is supposed to be more accurate so this should be the preferred way of calculating total net flux across boundaries.
By the way, I noticed you are talking about flux across moving boundaries. Not sure how COMSOL implementation is but mathematically this requires caution too. For example convective flux of a scalar quantity b is [b*n dot (v-vb)] where n is unit normal and v and vb are velocity vectors of the fluid and moving boundary with respect to a fixed eulerian frame. That is, you should check to see how it is implemented and account for this boundary "sweeping" effect.
Good luck
Ozgur
Reply | Reply with Quote | Send private message | Report Abuse
March 23, 2010 6:42am UTC in response to Ozgur Yildirim
Re: Flux in Convection and diffusion module : A quick question
Hi
Correct to point out the subtilities to consider with moving boundaries, I overlooked that one, some interesting references in this domain are for me:
Principles of computational Fluid Dynamicsby P. Wesseling, Springer 2009 (soft) and
Computational Fluid Mechanics and heat Transfer, by J.C. Tannehill, D.A. Anderson, R.H. Pletcher, Tylor&Fletcher, 1997
Have a nice reading
Ivar
Reply | Reply with Quote | Send private message | Report Abuse
April 8, 2010 9:58pm UTC in response to Ivar Kjelberg
Re: Flux in Convection and diffusion module : A quick question
It seems that nobody, even the support, knows how to calculate the flux using the weak form (lagrange multipliers), nor reacf() form.
Just to make sure one last time, does anybody know how to calculate flux in convection-diffusion using LMs or reacf()?
ps. reacf(c) always gives zero value (weak form is off)
Reply | Reply with Quote | Send private message | Report Abuse
Rules and guidelines
