Mathematics Asked by user541686 on December 25, 2020

I’m a bit confused how the Backward-Euler implicit integration method works for multiple variables (i.e., systems of differential equations).

For single variables, we solve the implicit equation

$$x(t + Delta t) = f(Delta t, x(t), dot{x}(t + Delta t))$$

But when I have (say) two variables, I’m not sure which variables’ future values I should use.

Do I solve

$$x_1(t + Delta t) = f_1(Delta t, x_1(t), dot{x}_1(t + Delta t), x_2(t), dot{x}_2(t + Delta t))$$

$$x_2(t + Delta t) = f_2(Delta t, x_2(t), dot{x}_2(t + Delta t), x_1(t), dot{x}_1(t + Delta t))$$

or do I solve

$$x_1(t + Delta t) = f_1(Delta t, x_1(t), dot{x}_1(t + Delta t), x_2(t), dot{x}_2(t))$$

$$x_2(t + Delta t) = f_2(Delta t, x_2(t), dot{x}_2(t + Delta t), x_1(t), dot{x}_1(t))$$

and why?

You are not even close. Please try to follow standard notation.

The differential equation is $dot x(t)=f(t,x(t))$. The implicit Euler method has the step $$ x(t+Δt)simeq x(t)+Δt·f(t+Δt, x(t+Δt)) $$ or in a more numerically oriented manner $$ x_{j+1}=x_j+Δt·f(t_{j+1}, x_{j+1}),qquad t_j=t_0+j·Δt, $$ or in correspondence with the Butcher tableau $$ k=f(t_{j+1},x_j+Δt·k)\ x_{j+1}=x_j+Δt·k $$ This formula is valid for scalar ODE as well as for vector ODE. If $L$ is the Lipschitz constant of $f$ in the $x$ direction and $|Δt|·L<1$ then the fixed point iteration of the first equation can be used as is. If this condition is not satisfied, for instance in the interesting case that the ODE is stiff, then one has to use the (simplified) Newton method. That is, with $J=partial_xf(t+Δt,x_j)$ you iterate $$ k_+=k-(I-Δt·J)^{-1}(k-f(t+Δt,x_j+Δt·k)) $$ If you stop this iteration after one step, you get the most simple Rosenbrock method.

Correct answer by Lutz Lehmann on December 25, 2020

2 Asked on February 22, 2021

algebraic geometry algebraic number theory commutative algebra integral extensions modules

2 Asked on February 22, 2021 by wyatt

1 Asked on February 22, 2021

2 Asked on February 22, 2021 by user330477

1 Asked on February 22, 2021

5 Asked on February 21, 2021 by arshiya

1 Asked on February 21, 2021 by michael-maier

1 Asked on February 21, 2021 by user160

1 Asked on February 21, 2021 by mashed-potato

0 Asked on February 21, 2021 by melisa-karimi

2 Asked on February 21, 2021 by annie-marie-heart

group theory lie algebras lie groups representation theory semi riemannian geometry

0 Asked on February 21, 2021 by tuong-nguyen-minh

1 Asked on February 20, 2021 by stanislas-castellana

1 Asked on February 20, 2021 by canine360

0 Asked on February 20, 2021 by sam-blitz

3 Asked on February 20, 2021

0 Asked on February 20, 2021 by red_trumpet

1 Asked on February 20, 2021 by rubenscube

1 Asked on February 19, 2021

Get help from others!

Recent Questions

- How Do I Get The Ifruit App Off Of Gta 5 / Grand Theft Auto 5
- Iv’e designed a space elevator using a series of lasers. do you know anybody i could submit the designs too that could manufacture the concept and put it to use
- Need help finding a book. Female OP protagonist, magic
- Why is the WWF pending games (“Your turn”) area replaced w/ a column of “Bonus & Reward”gift boxes?
- Does Google Analytics track 404 page responses as valid page views?

Recent Answers

- Peter Machado on Why fry rice before boiling?
- Lex on Does Google Analytics track 404 page responses as valid page views?
- Joshua Engel on Why fry rice before boiling?
- haakon.io on Why fry rice before boiling?
- Jon Church on Why fry rice before boiling?

© 2023 AnswerBun.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP