Work Packages Control of parameter dependent problems

Greedy algorithm for Parametric Vlasov-Fokker-Planck System

Authors: - 10 January 2018

Download Code
1. Numerical experiments

Consider the one dimensional linear Vlasov-Fokker-Planck (VPFP) as following.

In the numerical experiments, we compare the results from the greedy algorithm with the results from the polynomial interpolation of the quadrature points. we set $\epsilon = 1, T = 0.5$ for all examples. More specifically, we set

is the approximation solution obtained from either method, $f(T, x,v,z_i)$ is the exact solution. $z_i$ is uniformly chosen from $[a,b]$. ${\textbf{Dimension}}$ refers to the number of snap shots.

Observe Figure 1-6, we found that the Greedy algorithm converges to the solution faster than the other method. Especially, when there is discontinuity, the polynomial interpolation doesn’t converge as dimension increases, however, the Greedy algorithm still converges and almost keeps the same convergent rate.

Algorithm 1.1. The algorithm for the polynomial interpolation:

For Dimension $= n$, ${\alpha_i}_{i=1}^n$ is the n roots of $n$-th order Legendre polynomial, and let

Then by Lagrange interpolation,

1.1 Linear VFP with parametric initial data

In the first experiment, we consider the case where the parameter $z\in[1,5]$ only involved in the initial data,

Figure 1:

Figure 2:

Algorithm 1.2. The algorithm for PDE with only parameter in initial data: Offline:

Step 1: Set $N = {z_i: z_i = i\delta_z, \delta_z = \frac{b-a}{N_z}, i = 0,1, \dots, N_z}$, we set $N_z = 100$ in this case.

where $ \vert f \vert^2 = \int f^2 dxdv$.

Step 2: After obtaining $z_i^\ast, f_i = f_0(x,v,z_i^\ast)$, $i =1, \dots, n-1$,

where $P^{n-1} f_0$ is the projection onto the subspace spanned by $f_1, \dots, f_{n-1}$

Step 3: Stops when

Online: After we obtained the snap shots ${z_i^\ast}_{i=1}^n$,

Step 1: Orthogonalize ${f_i}_{i=1}^n$, and calculate the corresponding

Step 2: for any given $z$, one can approximate $f(T, x,v,z)$ by

figure1Figure 1: Continuous initial data

figure2Figure 2: Discontinuous initial data

1.2 Linear parametric VFP with deterministic initial data

In the second experiment, we consider the case where the parameter $z\in[1,5]$ only involved in the PDE,

Figure 3:

Figure 4:

Algorithm 1.3. The algorithm for PDE with only parameters in PDE: Offline: Set $N = {z_i: z_i = i\delta_z, \delta_z = \frac{b-a}{N_z}, i = 0,1, \dots, N_z}$, we set $N_z = 100$ in this case.

Step 1: choose

Step 2: After obtaining $z_i^\ast, f_i = f(T,x,v,z_i^\ast)$, $i =1, \dots, n-1$,

Step 3: Stops when the above maximum residual smaller than $\epsilon$.

Online: After we obtained the snap shots ${z_i}_{i=1}^n$, and the corresponding $f_i = f(T, x,v,z_i)$,

figure3Figure 3: Continuous parameter in PDE

figure4Figure 4: Discontinuous parameter in PDE

1.3 Linear parametric VFP with parametric initial data

In the third experiment, we consider the case where the parameter $z\in[1,5]$ involved in both the initial data and the PDE

Figure 5:

Figure 6:

Algorithm 1.4. The algorithm for PDE with parameters in PDE and initial data:

Offline: Set $N = {z_i: z_i = i\delta_z, \delta_z = \frac{b-a}{N_z}, i = 0,1, \dots, N_z}$, we set $N_z = 10000$ in this case.

Step 1: Using Algorithm 1.4 to obtain 50 snap shots out of $N$ to form a new set $N^\ast$. Step 2: Using Algorithm 1.3 to get snap shot ${ z_i^\ast }_{i=1}^n$

Online: Same as Algorithm 1.3.

Remark 1.5. In the experiments, the exact solution $f(T,x,v,z)$ refers to the numerical solution to VFP, that is, $f(T,x,v,z)$ is a $N_x\times N_v$ dimensional vector obtained by the following scheme

where $P(f)$ is the discretization of $\delta_v (M\delta_v (\frac{f}{M}))$,

Remark 1.6. For the matlab coding, the algorithm for the examples is in the “graph.m” file of the folder “eq1”, “eq2”, “eq3” respectively. You can change the coefficients in “sig_1”, “sig_2”, “sig_3”, and the electric field in “fcn_phi_x”, and initial data in “fcn_phi_x_0”, “fcn_rho_0”.

P.S. Remember to change the path of other functions at the beginning of “graph.m”.

figure5Figure 5: Continuous parameter in PDE and initial data

figure6Figure 6: Discontinuous parameter in PDE and initial data