Add temperature field variable to vertical component of force vector NS
Asked by
Joakim Bø
Hi!
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.
Thanks!
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
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.