Thin membrane using neumann conditions

 Register Blogs Members List Search Today's Posts Mark Forums Read

 October 28, 2013, 14:40 Thin membrane using neumann conditions #1 New Member   Eng. Emmanuel Luján Join Date: Mar 2013 Posts: 5 Rep Power: 6 Sponsored Links Hi!, I have a problem, I need to model a very thin membrane using neumann conditions. I used Salome Meca to create a mesh with a small sphere inside a box. The membrane is a facezone in the surface of the sphere. I created two patches that are together (cell_membrane_m and cell_membrane_s) from this facezone using createBaffles. My problem is that now I have to put a neumann condition like this: flux normal to the surface represented by cell_membrane_m = = flux normal to the surface represented by cell_membrane_s= = km * ( Potential (Vm) in cell_membrane_m - Potential (Vs) in cell_membrane_s ) (last term is the transmembrane potential) Here is what I did: 0/V ... cell_membrane_m { type groovyBC; variables "Vs{cell_membrane_s}=sum(V*mag(Sf()))/sum(mag(Sf()));km=0.01;Vm=sum(V*mag(Sf()))/sum(mag(Sf()));"; gradientExpression "km*(Vm-Vs)"; value uniform 0; fractionExpression "0"; } cell_membrane_s { type groovyBC; variables "Vm{cell_membrane_m}=sum(V*mag(Sf()))/sum(mag(Sf()));km=0.01;Vs=sum(V*mag(Sf()))/sum(mag(Sf()));"; gradientExpression "km*(Vm-Vs)"; value uniform 0; fractionExpression "0"; } As you can see I'm using groovyBC too. I'm using my own solver, based in lapplacian calculations. But I don't think it is the problem because I used it with other examples and worked well. The error I get is: #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/x86_64-linux-gnu/libc.so.6" #3 double Foam::sumProd(Foam::UList const&, Foam::UList const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::PCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::fvMatrix::solveSegregated(Foam::dict ionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #6 Foam::fvMatrix::solve(Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/pnp" #7 Foam::fvMatrix::solve() in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/pnp" #8 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/pnp" #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/pnp" Excepción de coma flotante (`core' generado) Can you help me please? Thank you! Emmanuel PD: This is an image about the boundary condition I want to set: https://www.dropbox.com/s/4ndt1pbm8sa1vtf/membrane.jpg

 Tags groovybc, neuman conditions, thin membrane

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Saturn CFX 45 February 8, 2016 05:42 beeo OpenFOAM Pre-Processing 20 July 17, 2013 08:39 kostas Main CFD Forum 9 August 17, 2012 05:13 s_h Main CFD Forum 0 August 9, 2009 22:15 Hema kothimbare FLUENT 0 August 7, 2009 01:22