Comments (4)
Hi
I believe it was drawn from the Simson KTH solver. But it's been a while, would be great if you could check it out😀
Mikael
from shenfun.
The scheme is indeed published (without reference) in the 2007 SIMSON technical report from KTH. The overall scheme is semi-implicit, paring an explicit Low-Storage RK33 (3-stage, 3rd-order) for nonlinear terms with Crank-Nicolson for linear terms, so rather than being an "IRK" or "SIRK" scheme, wherein any implicit part is baked into the RK, it would get labeled (by a method paper for RK schemes) something like RK(3,3)-2N + CN, where the 2N refers to the kind of low-storage scheme it seems to be.
I worked backwards to some early SIMSON method papers but never saw an explicit reference for the exact LSRK Butcher coefficients used, but at this point any modern LSRK reference, such as those pointed to by Nodepy, describe how to construct an RK(3,3)-2N scheme, so the only open question was why those particular coefficients were chosen over other examples.
The original scheme paper for all 2N LSRK schemes seems to be Williamson 1980.
Thanks for the help!
from shenfun.
Thanks a lot for the clarification! I have spent far too little time on the temporal integrators. Any upgrade or additions would be greatly appreciated😀
from shenfun.
I found a lead in Appendix D of Canuto, Quarteroni, Hussaini, and Zang (new editions, Vol. 2 on CFD methods), where they reference a Spalart and Watmuff 1993 JFM publication as providing "another common LSRK+CN" method due to Wray (unpublished), but in fact I finally found the SIMSON RK(3,3)-2N coefficients in a 1991 JCP paper from Spalart, Moser, and Rogers (https://doi.org/10.1016/0021-9991(91)90238-G), in which they say the coefficients are only "almost 3rd order", so the original Williamson 1980 (3,3) coefficients or the (4,3) coefficients from Carpenter and Kennedy (1994) might provide better results just by swapping coefficients (and bumping the rk loops by 1 stage), but obviously the best person to ask about optimal coefficients would be an integrator expert like David Ketcheson. I'm just a turbulence modeling/theory guy 😃.
Anyway, that concludes the hunt for where the coefficients came from.
I am curious if you'd ever tested any of the IRK3/ETD/ETDRK4 integrators from shenfun using a CFL-based adaptive time step, and thus had to call the setup() function for the solver over and over? Would this be crazy inefficient, or would the alltoall cost of a large-scale 3D HIT or channel flow make the solver matrix-assembly steps seem utterly cheap in comparison? If you haven't ever tried then I'll definitely run some tests to see what I could get away with. I'm running decaying HIT cases using a slight modification of your shenfun_NS example, basically the spectralLES module I developed as a PhD student (https://github.com/ctowery/teslapy/tree/master/spectralLES) but upgraded to shenfun, and since the maximum time step increases dramatically over time, there will be a time interval where re-running the setup for a bigger dt will be worth it. (Question is if that interval is closer to every step or 1000 steps.)
from shenfun.
Related Issues (20)
- Dx not working with 1D mixed function space HOT 4
- Issue on Stokes equations with doubly periodic boundary conditions. HOT 2
- Error when evaluating derivatives for nonstandard domain HOT 7
- Installation problem: error in "cimport fastgl_wrap" HOT 2
- Old shenfun script for 2D Poisson equation dosn't work with new version HOT 2
- 1d poisson HOT 1
- Bug in FunctionSpace with bc as dictionary HOT 4
- Solve DGL with non-constant (discrete) variable? HOT 5
- AttributeError: 'tuple' object has no attribute 'copy' HOT 1
- Install problem HOT 1
- problems in solving 2D Allen Cahn equation HOT 2
- Two questions in the demo MixedPoisson.py HOT 7
- How to compute the numerical solution on the general points HOT 4
- 1. Solver for multidimensional problem 2. Product with variable coefficients HOT 4
- Two questions about the Legendre-Galerkin method for the Cahn-Hilliard equation HOT 2
- How to compute the integral about the nonlinear term in the Navier-Stokes equation? HOT 2
- How to assemble the right-hand side vector? HOT 3
- Understanding the chebyshev solver in poisson2ND.py HOT 3
- Time dependent boundary conditions in segments for Rayleigh-Benard 2D HOT 2
- Implementing the one-dimensional Kuramoto-Sivashinsky system HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from shenfun.