March 2021 - Nektar-users - imperial.ac.uk

Rotational arguments for periodic faces
by Vidhi Zala 28 Mar '21

28 Mar '21

******************* This email 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. ******************* Hi everyone, I am trying to implement a 3D advection example on a cube with composite elements (hex+tets+pyramids+quads+tri). I am using periodic boundary conditions. I am getting the following error: "Unable to determine face orientation". I referred to the relevant section of user guide (please see screenshot). I am unable to understand how should I change the rotation argument such that I can simulate a periodic flow. It would really help if I get some more details about this. Many thanks, V ,-~~-.___._. / | ' \ | |--------. Vidhi Zala ( ) 0 | | | School of Computing \_/-, ,----' | | | Scientific Computing and Imaging Institute ==== !_!-v---v--. University of Utah / \-'~; .--------. / __/~| ._-""|| | vidhi.zala(a)utah.edu =( _____|_|____||________| http://sci.utah.edu/people/vidhi.html I support Ruff Haven Crises Shelter for Homeless Pets <https://www.ruffhaven.org/donate>. Consider a donation today :)

2 3

Outflow Boundary Condition for ShallowWaterSolver
by Boyuan Yu 27 Mar '21

27 Mar '21

******************* This email 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. ******************* Hello all, I am using ShallowWaterSolver and my application requires an outflow boundary condition. However, there seems no appropriate outflow boundary condition for ShallowWaterSolver (the "HOutflow" and "PressureOutflow" boundary conditions seem to only work for IncNavierStokesSolver). So how could the outflow boundary condition be implemented for ShallowWaterSolver? The "outflow" boundary condition here only requires that the field variables advected out of the domain without reflection using some kind of extrapolation. Thank you. Boyuan Yu

2 1

Mesh growth rate control for NekMesh
by Shujie Jiang 27 Mar '21

27 Mar '21

******************* This email 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 Nektar users, I am using NekMesh to generate mesh for my simulation. I found some reflection on the edge of the coarse mesh and the refine mesh (like the picture below). [cid:3879835e-1dea-4ad8-952e-7a020d496fdd][cid:f2252f49-e539-4496-aafd-4b0319898a27] Can I control the growth rate of mesh to give a smooth transition from a refined mesh to a coarse mesh? Thank you very much! Kind regards, Shujie University of southampton

2 1

A Question Concerning how to Include Source Terms in ShallowWaterSolver
by Boyuan Yu 27 Mar '21

27 Mar '21

******************* This email 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. ******************* Hello all, I would like to express my appreciation for the efforts made by the developer team and community to create such a great open-sourced software. Currently, I am working with shallow water equations and my application requires the inclusion of the source terms of: (1) bed-friction term; (2) channel slope term. However, these particular source terms are not included in ShallowWaterSolver of Nektar++ yet and it seems to be not possible to use any forcing terms in ShallowWaterSolver. As far as I am concerned, the method to use these source terms is to modify the solver ShallowWaterSolver source code and recompile all the library. But this method is quite time-consuming. So I am writing to ask is there any simpler way to implement this functionality. Thank you in advance for your answer. Boyuan Yu

2 1

Re: [Nektar-users] Residuals Stability Simulation
by Gao, Ankang 15 Mar '21

15 Mar '21

Hi Matt, If you are using the ModifiedArnoldi driver, the Residual after each iteration is shown in the *.evl file, for example, -- Iteration = 10, H(k+1, k) = 1.2317e-01 Magnitude Angle Growth Frequency Residual EV: 0 1.0000e+00 0.0000e+00 1.0701e-08 0.0000e+00 2.8913e-07 EV: 1 7.1131e-01 0.0000e+00 -3.4064e-01 0.0000e+00 6.8432e-05 Cheers Ankang 发件人:"Matt Duran" <matt.duran60(a)gmail.com> 发送时间:2021-03-15 08:33:45 (星期一) 收件人: Nektar-users(a)imperial.ac.uk 抄送: 主题: [Nektar-users] Residuals Stability Simulation This email from matt.duran60(a)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. Hello, I was running a stability simulation using the IncNavierStokesSolver and was wondering if there was some way to check or display the residuals after each iteration. Thank you for your time, Matt Duran Mechanical and Aerospace Engineering University of Central Florida

1 0

white noise
by Alessandro Alberto de Lima 15 Mar '21

15 Mar '21

******************* This email 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. ******************* Hi Guys. Can you help me setting up a white noise during a few steps on m simulation? On the section 3.5.5 of the user guide I found this: <FORCE TYPE="Noise"> <WHITENOISE> [VALUE] <WHITENOISE/>  <UPDATEFREQ> [VALUE] <UPDATEFREQ/> <NSTEPS> [VALUE] <NSTEPS/> </FORCE> What's the meaning of those white noises? I need to specify one for each direction on three dimension problems? I need to assign that noise to a specific boundary region (a wall, for example)? Thank you. Alessandro

2 1

Fieldconvert question
by Vidhi Zala 15 Mar '21

15 Mar '21

******************* This email 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. ******************* Hi all, I have limited experience with FieldConvert utility and need help with the module interppointdatatofld. Here is the way I am using it ./FieldConvert-g -m interppointdatatofld:frompts=tmp.pts compositemesh2d.xml file1.dat Is it possible to have the output as vtk/vtu file instead of .dat file? Most other modules in FieldConvert allow output as both .dat and .vtu. Any help will be appreciated! Thanks, V ,-~~-.___._. / | ' \ | |--------. Vidhi Zala ( ) 0 | | | School of Computing \_/-, ,----' | | | Scientific Computing and Imaging Institute ==== !_!-v---v--. University of Utah / \-'~; .--------. / __/~| ._-""|| | vidhi.zala(a)utah.edu =( _____|_|____||________| http://sci.utah.edu/people/vidhi.html I support Ruff Haven Crises Shelter for Homeless Pets <https://www.ruffhaven.org/donate>. Consider a donation today :)

2 1

Residuals Stability Simulation
by Matt Duran 15 Mar '21

15 Mar '21

******************* This email 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. ******************* Hello, I was running a stability simulation using the IncNavierStokesSolver and was wondering if there was some way to check or display the residuals after each iteration. Thank you for your time, Matt Duran Mechanical and Aerospace Engineering University of Central Florida

2 1

Roots of Jocobi or Legendre polynomial?
by Dongmyung Suh 07 Mar '21

07 Mar '21

******************* This email 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 All, I am a user of Nektar ++, and learning its foundation with Sherwin's book, Spectral hp element method for CFD [2ed]. I have a small question, while I study the book. When the book explains the Gauss-Lobatto-Legendre formula for integration (pg.61) and differentiation (pg 68), xi are the roots of Jocobi polynomial. However, in other references ( https://en.wikipedia.org/wiki/Gaussian_quadrature) they are the roots of Legendre polynomial. Which one is right? or both are correct? Please let me be clear about that. Thank you. Dongmyung Suh

2 2

setting-up bc's for Compressible Navier-Stokes solver
by wasim sarwar 07 Mar '21

07 Mar '21

******************* This email 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 All, Good Morning, and wishing all of you very good day! I am trying to solve the problem of convergent-divergent nozzle with pressure at inlet and outlet. Does anyone have the idea how to setup the pressure inlet and pressure outlet boundary conditions in compressible Navier-Stokes solver. Also, give an idea of how-to setup the symmetry bc about x-axis. The input file is also attached herein for reference. Mant thanks in advance Wasim Sarwar UPC, Barcelona

2 1