CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Solve Coupled Equations

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 12, 2008, 01:08
Default Hi I would like to solve c
  #1
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 206
Rep Power: 9
nandiganavishal is on a distinguished road
Hi

I would like to solve coupled equations (using relaxation or Newtons method, but first would like to see for relaxation method). Can you help me which tutorial or solver should i follow in order to solve this problem.

Further I have solved a poisson equation adding source terms that are functions of x and y successfully. I had incorporated Dirichlet B.C's for a square 1*1 block (the b.c's are also functions of x and y). May I know how do we incorporate neumann b.cs for the same problem (they are also functions of x and y).

Say the four sides are inlet,outlet, upperwall and lowerwall.

My B.C's are

inlet (AT x =0) - Dirichlet B.C's = (y*(1-y)) (something of this
type).
Similarly (at y=0) (lowerwall) = -x^3 + sin (x)..

Neumann B.C's at x =1 (outlet) = -3-2*sin(x)(some arbitary function)
Neumann B.C's at y =1 (upperwall) = 4+sin(y)

For the inlet case I have my B.C statements like these... Which works fine..
label inletPatchID = mesh.boundaryMesh().findPatchID("inlet");
fvPatchScalarField& inletT = T.boundaryField()[inletPatchID];
const vectorField& faceCentresInlet =
mesh.Cf().boundaryField()[inletPatchID];
forAll(inletT, pointI)
{
inletT[pointI] = faceCentresInlet[pointI].y()*(1.0 -
faceCentresInlet[pointI].y());
}

Can anyone let me know how do we incorporate the neumann b.cs for this problem..

Kindly help me out..

Regards

Vishal
nandiganavishal is offline   Reply With Quote

Old   November 12, 2008, 23:47
Default Hi, I need to solve a coupl
  #2
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 206
Rep Power: 9
nandiganavishal is on a distinguished road
Hi,

I need to solve a coupled equation which is of the form
laplacian(phi) = K * (C1-C2) ( K is some constant)
ddt(C1) = laplacian (phi,C1)+ extra terms;
ddt(C2) = laplacian(phi,C2)+ extra terms;

Could anyone let me know how exactly should i got go about solving this equation. Basically i want to solve a PNP equation, and I don't want it to be solved using relaxation method, I want to solve it using a single matrix (using some newton method or any other iteration techniques).

Can anyone throw some light with regard to this, as to what example i have to follow and how I should go about doing this.

Kindly help.

Thanks

Vishal
nandiganavishal is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Solve coupled equationss using direct solvers nandiganavishal OpenFOAM Running, Solving & CFD 0 November 23, 2008 18:21
How to solve coupled transport equations nandiganavishal OpenFOAM Running, Solving & CFD 4 November 22, 2008 17:43
Coupled equations solution litonx OpenFOAM Running, Solving & CFD 7 February 14, 2007 19:31
How to solver coupled transport equations chris1980 OpenFOAM Running, Solving & CFD 7 August 25, 2006 08:18


All times are GMT -4. The time now is 07:29.