Dear Professor Kumar, Many thanks for your help. Having obtained Nektar++ 4.5, I can run your example and I get output resembling your graphic. Looking at the xml file, it is a bit of a mystery to me how Nektar++ knows that it has to treat the temperature as a linear advection-diffusion problem (is this a default behaviour?). Anyhow, thanks for the excellent support. Regards, Ed. On Wed, Aug 19, 2020 at 12:39 AM Abhishek Kumar <abhishek.kir@gmail.com> wrote:
Dear Edward,
For Rayleigh-Benard convection, a BodyForce function is applied for the buoyancy term. Presently, it works perfectly with Nektar 4.5, but there is some issue in Nektar 5 for this function.
I have attached a xml file, which solves the Rayleigh-Benard problem for Pr=0.71 and Ra=5000. Please see section 3.2 of Kumar and Potherat JFM 2020 for the details of validation (attached in this email).
Please download this version nektar++ from the following link
https://gitlab.nektar.info/nektar/nektar/-/tree/b9f265928a72241b886a86eaedf4...
This version corresponds to "Commit b9f26592".
Let me know if you have any further queries.
With regards Abhishek
---------------------------------------------------------------------------------------
Abhishek Kumar
Assistant Professor (Research)
Centre for Fluid and Complex Systems
Coventry University, Coventry CV15FB
The United Kingdom
--------------------------------------------------------------------------------------- On Tue, Aug 18, 2020 at 8:12 PM Edward Threlfall <ejthrelfall@gmail.com> wrote:
This email from ejthrelfall@gmail.com originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list <https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address.
Dear users,
I am trying to set up a simulation of Rayleigh-Benard convection in Nektar++.
I have seen a poster by Hossein and Sherwin from last year's workshop which seems to do a similar sort of thing, but, based on my (limited) understanding of the framework, I am not sure how this was done or how to set up something similar myself (e.g. the system is coupled incompressible (?) Navier-Stokes and advection-diffusion; so, was it the case that the two solvers were run with coupling between them?).
I note there is a separate thread running on convection at the moment; I could have emailed the poster's author directly, but perhaps a public discussion of how to set up this type of problem will benefit others too.
Many thanks,
Dr. Ed Threlfall (UKAEA). _______________________________________________ Nektar-users mailing list Nektar-users@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/nektar-users
--