|
[Sponsors] |
June 1, 2011, 09:16 |
flux source at particular location(boundary)
|
#1 |
New Member
Churamani Prasad
Join Date: May 2011
Posts: 16
Rep Power: 14 |
Hi.
I need to find the derivative of a variale (reactant.molconc) at particular location(x=0.002) which i m able to do it and then specify two sources, NEGATIVE SOURCE at left of the location and POSITIVE SOURCE at the right of the location, where source is some simple function of the gradient of the reactant.molconc . The geometry is simple reactangle with a boundary at x=0.002 I thougt of a way but don't know(90% sure that it won't work) if it will work the part of the code i written is : DO ILOC = 1, NLOC IF (X(ILOC).GE.0.002 .AND. X(ILOC).LE.0.002 .AND. & Y(ILOC).GE.0.00 .AND. Y(ILOC).LE.0.001) THEN SOURCE(ILOC-1) = -(1e-9)*GRAD_PHI(1,ILOC) SOURCE(ILOC+1) = (1e-9)*GRAD_PHI(1,ILOC) END IF END DO but i m not sure if ILOC-1 will be just left of the location and ILOC+1 at the right of the boundary. The problem is how the vertices are assigned to NLOC.. Plz help on this...Suggest some other and better ways to do the same. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] swak4foam building problem | GGerber | OpenFOAM Community Contributions | 54 | April 24, 2015 16:02 |
pisoFoam compiling error with OF 1.7.1 on MAC OSX | Greg Givogue | OpenFOAM Programming & Development | 3 | March 4, 2011 17:18 |
OpenFOAM on MinGW crosscompiler hosted on Linux | allenzhao | OpenFOAM Installation | 127 | January 30, 2009 19:08 |
DxFoam reader update | hjasak | OpenFOAM Post-Processing | 69 | April 24, 2008 01:24 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 17:51 |