reaction-diffusion

113 views
Skip to first unread message

Paul Van Liedekerke

unread,
Aug 7, 2013, 12:26:43 PM8/7/13
to dea...@googlegroups.com
Hi

(New here)

I am exploring the possibilities to couple deal to our software which models the growth of biological cells. Therefore we would need to solve a reaction-diffusion equation (of chemicals) with the cells as (point) sources.
The cells are generally closely packed.
 
We would like Deal to solve these equations and return the solution to our software. However, I have no clear view on the capabilities or the efficiency of deal in this regard..; so could anyone give feedback whether this is a valuable solution ?

Many thanks for any input !

Paul

Armin Ghajar Jazi

unread,
Aug 8, 2013, 2:15:05 PM8/8/13
to dea...@googlegroups.com
Hi Paul,

I am not familiar with with particular physics you are talking about, however I have managed to use Dealii in our FluidSolver code. My code deals with partitioned fluid structure interaction. What I am doing essentially is solving the structural part in Dealii (based on step44 tutorial) and the fluid part in our code.

The first difficulty is of course compiling the two together, which is highly dependent on the compiling routine you are using(if it is with cmake , gnu autochain tools etc.). The new version of Dealii uses cmake where as previous versions were autochain based. After that, getting Dealii functionality in your code is just the matter of including appropriate Header files into your code.

The second part is to decide if both programs live on the same mpi-communicator. In my program Dealii lives on another communicator than my code, and I simply send the information via MPI calls back and forth, this of course necessitates a coupling environment which holds both communicators.

I would be glad to help you on this if you need assistance in any part I am familiar with!

cheers,
Armin




Paul

--
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+un...@googlegroups.com.
For more options, visit https://groups.google.com/groups/opt_out.
 
 

Paul Van Liedekerke

unread,
Aug 12, 2013, 4:50:45 PM8/12/13
to dea...@googlegroups.com
Hi Armin

Thanks for your feedback, yet I am not that far  :-) ..I am currently going through the steps of Deal.II. I think a difficulty will be the coupling of discrete concentration sources (cells) with the continuum environment in FEM. If you have any knowledge on this issue, I ll be glad to hear it !

Cheers

Paul



2013/8/8 Armin Ghajar Jazi <arm...@gmail.com>
You received this message because you are subscribed to a topic in the Google Groups "deal.II User Group" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/dealii/QHEXQw_lk_E/unsubscribe.
To unsubscribe from this group and all its topics, send an email to dealii+un...@googlegroups.com.

Wolfgang Bangerth

unread,
Aug 12, 2013, 4:05:54 PM8/12/13
to dea...@googlegroups.com

> Thanks for your feedback, yet I am not that far :-) ..I am currently going
> through the steps of Deal.II. I think a difficulty will be the coupling of
> discrete concentration sources (cells) with the continuum environment in FEM.
> If you have any knowledge on this issue, I ll be glad to hear it !

If these are point sources, then you'll want to know about
VectorTools::create_point_source_vector().

Best
W.


--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@math.tamu.edu
www: http://www.math.tamu.edu/~bangerth/

Reply all
Reply to author
Forward
0 new messages