New answers tagged

3

A search for the specific coefficients listed led me to the method ROS3PRL from J. Sieber, Konvergenzanalyse und Numerische Tests für die Prothero–Robinson–Gleichung (Master thesis), TU Darmstadt, 2014. I can't seem to find this thesis online, but the method is mentioned in the following which may be of interest. Rang, Joachim. "Improved traditional ...


4

Crank-Nicolson is a very good classical approach for parabolic PDE like the heat transfer PDE to which it was originally applied. It is relatively easy to understand and implement so it is often presented in basic courses on numerical methods for PDE. pdepe is also very well-suited to this class of PDE (the second "p" in pdepe stands for parabolic). It has ...


2

Since this is a question with pretty subjective answers, I'll add a couple to Prof. Bangerth's very good list. the theorem of adjoint/dual operators and spaces is pretty crucial to Computational Science. We know that dual-consistent discretizations of the PDEs can obtain superconvergence which is a nice property. But I think the more commonly used outcomes ...


1

So I ran with my adjoint solver, but without a singularity on a sequence of refined meshes, and didn't see that behavior, as expected. My guess is as follows: If you think of the adjoint vector as a green's function, relating perturbations in the residual operator to a delta in the objective function (in this case weighted projections of pressure), I ...


10

You'll get everyone to give different answers to this question, and maybe that's alright. Here are some of my favorite ones: Taylor's theorem that a function (of sufficient smoothness) equals its Taylor expansion plus a remainder term. One can consider the Bramble-Hilbert lemma as a variation of Taylor's theorem, but it has different applications and is ...


26

There are a few ways to conserve energy during ODE integration. Method 1: Symplectic Integration The cheapest way that is to use a symplectic integrator. A symplectic integrator solves the ODE on a symplectic manifold if it comes from one, and so if the system comes from a Hamlitonian system, then it will solve on some perturbed Hamiltonian trajectory. ...


Top 50 recent answers are included