Thank you Denis. I tried that. The function log(c/(1-c)) looks like:
<https://lh3.googleusercontent.com/-98tPVjPUANE/V9KVjOfyXdI/AAAAAAAAAAM/OMFk_X10nP82jlxIfN1cEjtDx0Fvdf9lQCLcB/s1600/log.png> and it is difficult to linearise in the whole range. And the assembled matrices are solved out of deal.II so I definitely need to a weak form. Regards Dragan On Friday, September 9, 2016 at 6:19:15 AM UTC+1, Denis Davydov wrote: Hi Dragan, > > Did you do a literature review to see what researchers are doing in your > field? > Most likely you need to linearize it w.r.t. to "c" so that you will be > solving for \Delta c around some points c_0 employing something like > Newton-Raphson method. > > Regards, > Denis. > > On Thursday, September 8, 2016 at 10:48:43 PM UTC+2, dragn...@gmail.com > wrote: >> >> Hello guys, >> >> I would need a small help with deriving the weak form. The equation >> without the unnecessary terms is the following: >> log(c/(1-c)) - k*Laplacian(c) = 0 >> I am not sure how to treat the term log(c/(1-c)). >> Thank you. >> >> Regards, >> Dragan >> > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
