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.
total weight calculation for a domain with variable density
Posted Dec 28, 2015, 7:36 a.m. EST Materials, Parameters, Variables, & Functions, Studies & Solvers 2 Replies
Please login with a confirmed email address before reporting spam
I’m trying to model a polymer degradation process with Comsol in which I introduced degradation to my model by defining variable concentration of one of its ionomer blocks.
Then I defined its Density as: “(MW1[g/mol]*cIon1[mol/m^3])+const”.
And for calculating polymer weight loss during degradation process: I defined integration over Density.
cIon1=-R1 (which is degradation reaction)
I made sure that degradation reaction is positive by (R1=max(esp^2,r1))
However as I integrate Density on my domain it decreases and then increases after certain time!!!! Please help!!
Then I defined its Density as: “(MW1[g/mol]*cIon1[mol/m^3])+const”.
And for calculating polymer weight loss during degradation process: I defined integration over Density.
cIon1=-R1 (which is degradation reaction)
I made sure that degradation reaction is positive by (R1=max(esp^2,r1))
However as I integrate Density on my domain it decreases and then increases after certain time!!!! Please help!!
2 Replies Last Post Jan 3, 2016, 2:59 a.m. EST
Please login with a confirmed email address before reporting spam
Posted:
9 years ago
Dear Aida,
You probably mean “clon1t=-R1”. I suggest plotting clon1 to check that it remains positive, and also using a tighter convergence tolerance in your solver.
Nagi Elabbasi
Veryst Engineering
You probably mean “clon1t=-R1”. I suggest plotting clon1 to check that it remains positive, and also using a tighter convergence tolerance in your solver.
Nagi Elabbasi
Veryst Engineering
Please login with a confirmed email address before reporting spam
Posted:
9 years ago
Dear Aida,
You probably mean “clon1t=-R1”. I suggest plotting clon1 to check that it remains positive, and also using a tighter convergence tolerance in your solver.
Nagi Elabbasi
Veryst Engineering
Dear Nagi,
Thank you for your reply.I checked cIon1 and as you said it becomes negative so I modified my model:
CIon1=max(eps^2,cIon1)
W=CIon1*MW1+const
cIon1t=-R1
R1=max(eps^2,r1)
However the problem is still there, W increases after a certain time. I probe max value for CIon1 in the domain and the confusing fact is the it increases more than the initial value.
Please help.
Aida