Add temperature field variable to vertical component of force vector NS

Asked by Joakim Bø


I'm currently constructing a Chorin/IPCS Navier Stokes solver using Fenics. The temperature/energy equation is also incorporated.

So, my question is, I'm going to use the Boussinesq bouyancy approximation for incompressible fluids, but how do I apply this to my force vector?

The temperature field T1 in the last time step is known, and I want to add it to my force vector so it becomes f = 0i + T0[i,j]j in the spatial point i, j for the next time step. Any tips on the easiest way of doing this?

This is just for testing purposes, the full model is going to be implemented at a later stage.


Question information

English Edit question
FEniCS Project Edit question
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Johannes Ring (johannr) said :

FEniCS no longer uses Launchpad for Questions & Answers. Please consult the documentation on the FEniCS web page for where and how to (re)post your question:

Can you help with this problem?

Provide an answer of your own, or ask Joakim Bø for more information if necessary.

To post a message you must log in.