This package contains the DiffusionFlux block for the computation of fluxes at the cell boundaries.
Release Notes:
| Name | Description |
|---|---|
PDE.FiniteVolume.Fluxes.DiffusionFlux.DiffusionFlux
Implements the diffusion flux
By using this flux, the average update rule becomes:
Release Notes:
| Type | Name | Default | Description |
|---|---|---|---|
| Integer | n | worldModel1.n | |
| Integer | gcl | worldModel1.gcl | |
| Integer | gcr | worldModel1.gcr | |
| Real | beta | 1.0 | |
| Real | deltax | 1/n |
| Type | Name | Description |
|---|---|---|
| input RealInput | u[worldModel1.n + worldModel1.gcl + worldModel1.gcr] | |
| output RealOutput | y[worldModel1.n + 1] |
block DiffusionFlux
extends Icons.BlockIcon;
outer PDE.World.worldModel worldModel1;
parameter Integer n = worldModel1.n;
parameter Integer gcl = worldModel1.gcl;
parameter Integer gcr = worldModel1.gcr;
parameter Real beta = 1.0;
parameter Real deltax = 1/n;
equation
// for i in 1:gcl loop
// y[i] = u[i];
// end for;
for i in gcl:gcl+n loop
y[i-gcl+1] = (-beta/(deltax))*(u[i+1] - u[i]);
end for;
// for i in gcl+n+1:gcl+n+gcr loop
// y[i] = u[i];
// end for;
public
Modelica.Blocks.Interfaces.RealInput u[worldModel1.n + worldModel1.
gcl + worldModel1.gcr];
Modelica.Blocks.Interfaces.RealOutput y[worldModel1.n + 1];
end DiffusionFlux;