Hi Wei,

thanks for getting back to me. 

I found the error in the code of PulseWaveSystem.cpp and in the analytical solution of the inverse Jacobian in BetaPressureArea.cpp.  

I also found that the solver is hardcoded for junctions of 3 vessels or less and I am confident that the example Merging.xml is incorrect. 

I added trifurcations and confluences and currently I am trying to add junctions with multiple inlet/outlets. 

BW,

Kai 
From: nektar-users-bounces@imperial.ac.uk <nektar-users-bounces@imperial.ac.uk> on behalf of Wei He <weiyh.he@gmail.com>
Sent: Wednesday, July 3, 2024 15:33
To: 🌸那慢慢🌸 <2644878715@qq.com>
Cc: nektar-users <nektar-users@imperial.ac.uk>
Subject: Re: [Nektar-users] Adaptive SFD
 
You don't often get email from weiyh.he@gmail.com. Learn why this is important
I can see the symmetry bcs for u is defined as D u=1, maybe you can set it as N u=0.

Regards,
Wei

On Tue, 2 Jul 2024 at 12:31, 🌸那慢慢🌸 <2644878715@qq.com> wrote:

This email from 2644878715@qq.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 to disable email stamping for this address.

 

Dear nektar++ users and developers,

I am using the Adaptive SFD method to calculate the flow field information when a square column experiences vortex shedding. But I have calculated many times according to the prompts in the user guide document, and the following two results are displayed:

Vectors are lineary independent

NaN found during time integration

How should I set up an XML file to ensure correct calculation? I will attach my calculation file.

I run the command ：IncNavierStokesSolver Cylinder-4.xml.gz Cylinder-4.xml

Hongtao Cheng

_______________________________________________
Nektar-users mailing list
Nektar-users@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/nektar-users