This package contains deltaQ block that computes the jumps at the cell boundaries.
Release Notes:
| Name | Description |
|---|---|
PDE.FiniteVolume.FluxLimiter.DeltaQ.deltaQ
Takes the average matrix Q and computes the jumps
giving matrix deltaQ as output. The entry delta_q[i, j] of the matrix contains the jump of the i-th unknown in the system at the interface j.
Release Notes:
| Type | Name | Default | Description |
|---|---|---|---|
| Integer | n | worldModel1.n | |
| Integer | m | worldModel1.m | |
| Integer | gcl | worldModel1.gcl | |
| Integer | gcr | worldModel1.gcr | |
| Integer | p | 1 |
| Type | Name | Description |
|---|---|---|
| input RealInput | u[worldModel1.m, worldModel1.n + worldModel1.gcl + worldModel1.gcr] | |
| output RealOutput | y[worldModel1.m, worldModel1.n + worldModel1.gcl + worldModel1.gcr - 1] |
block deltaQ
extends Icons.BlockIcon;
outer PDE.World.worldModel worldModel1;
parameter Integer n = worldModel1.n;
parameter Integer m = worldModel1.m;
parameter Integer gcl = worldModel1.gcl;
parameter Integer gcr = worldModel1.gcr;
parameter Integer p = 1;
equation
for i in 1:m loop
for j in 1:n+gcl+gcr-1 loop
y[i, j] = u[i, j+1] - u[i, j];
end for;
end for;
public
Modelica.Blocks.Interfaces.RealInput u[worldModel1.m,worldModel1.n +
worldModel1.gcl + worldModel1.gcr];
Modelica.Blocks.Interfaces.RealOutput y[worldModel1.m,worldModel1.n
+ worldModel1.gcl + worldModel1.gcr - 1];
end deltaQ;