Adding a cross product source term
I want to solve for 2-Fluid (ion-electron) steady state flow. Now it seems that between mhdFoam and simpleFoam this should be possible. In essence I need to add a v x B term to the simpleFoam velocity solver. My guess is this should look something like:
Code:
{ |
Followup
Well I got a version to compile anyway, now I just have to figure out how to set my initial condition for my solver (gradient fields, etc). Here's the code:
Code:
volVectorField UxF; |
Hello lazersos,
as far as i get it you add a source term called UxF and then you have to integrate it over the control volume. For that you call "+fvm::Sp(q,UxF)" to integrate UxF over each cell volume, right? But what is q here? thanks for a reply, rupert |
This is what I got to compile:
Code:
UxF = U ^ F; |
All times are GMT -4. The time now is 21:23. |