Menu

Discussion Forum

New Discussion
Filter

Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Diffusion in two layers!

Posted Feb 17, 2015, 7:43 a.m. EST Version 4.4 4 Replies

Mikael Hedenqvist
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Hi, I am pretty new on Multiphysics. I want to make a solid sphere with two layers. I want to calculate the mass uptake in the inner layer (from origo to the border between the layers). As a first approximation there will be two diffusion coefficients but the same solubility in both. The boundary will be a continous flux (in the next case two different solubilities will be used). How do I set up the system. As I understand there are various ways of doing this, two components, two physics and two domains, please advise how to make sucg systems.
4 Replies Last Post Feb 24, 2015, 4:08 a.m. EST
Lasse Murtomäki
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 10 years ago Feb 17, 2015, 9:12 a.m. EST
Hi

1. Make an axisymmetric 2D-model.
2. Draw a semicircle with a desired layer thickness.
3. Mesh across the boundary very densely (attached).
4. No boundary condition needed between the domains.
5. Initial concentration in the inner sphere = 0.
6. Compute with logarithmic time step, because initially there is an enormous concentration gradient across the boundary.
7. Surface integration in the inner sphere, activate "Compute volume integral".
8. Check that the total amount = volume integral over entire geometry = constant = 4/3·pi.

Thus, only one physics (Transport of Dilutes Species), one component. I make always dimensionless models, if possible, like this:

R = r/r0 (r0 = outer sphere radius)
D = 1, because
T = Dt/r0² (dimensionless time)
C = c/c_init (c_init = initial conc. in the outer sphere)

At T = 0.1 80% of the material is in the inner sphere.

br
Lasse
Hi 1. Make an axisymmetric 2D-model. 2. Draw a semicircle with a desired layer thickness. 3. Mesh across the boundary very densely (attached). 4. No boundary condition needed between the domains. 5. Initial concentration in the inner sphere = 0. 6. Compute with logarithmic time step, because initially there is an enormous concentration gradient across the boundary. 7. Surface integration in the inner sphere, activate "Compute volume integral". 8. Check that the total amount = volume integral over entire geometry = constant = 4/3·pi. Thus, only one physics (Transport of Dilutes Species), one component. I make always dimensionless models, if possible, like this: R = r/r0 (r0 = outer sphere radius) D = 1, because T = Dt/r0² (dimensionless time) C = c/c_init (c_init = initial conc. in the outer sphere) At T = 0.1 80% of the material is in the inner sphere. br Lasse

Attachments:

Lasse Murtomäki
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 10 years ago Feb 18, 2015, 7:04 a.m. EST
Hi

And the model file.

Lasse
Hi And the model file. Lasse
Lasse Murtomäki
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 10 years ago Feb 24, 2015, 1:46 a.m. EST
I have only version 5.0.

br

Lasse
I have only version 5.0. br Lasse
Mikael Hedenqvist
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 10 years ago Feb 24, 2015, 4:08 a.m. EST
OK, thanks anyway,

Regards

Mikael
OK, thanks anyway, Regards Mikael

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.

Suggested Content