This package contains a block that computes the wave strengths.
Release Notes:
| Name | Description |
|---|---|
PDE.FiniteVolume.Fluxes.Roe.WaveStrength.a
Takes the averages rho, a and differences deltarho, deltav and deltap as input and computes the three wave strengths
Release Notes:
| Type | Name | Default | Description |
|---|---|---|---|
| Integer | n | worldModel1.n |
| Type | Name | Description |
|---|---|---|
| input RealInput | u[worldModel1.n + 1] | |
| input RealInput | u1[worldModel1.n + 1] | |
| input RealInput | u2[worldModel1.n + 1] | |
| input RealInput | u3[worldModel1.n + 1] | |
| input RealInput | u4[worldModel1.n + 1] | |
| output RealOutput | y[worldModel1.n + 1] | |
| output RealOutput | y1[worldModel1.n + 1] | |
| output RealOutput | y2[worldModel1.n + 1] |
block a extends Icons.BlockIcon; outer PDE.World.worldModel worldModel1; parameter Integer n = worldModel1.n; equation for j in 1:n+1 loop //Inserting 0.0001 to avoid division by zero y[j] = (1/(2*((u1[j]+0.0001)^2)))*(u4[j] - u[j]*u1[j]*u3[j]); y1[j] = (1/(((u1[j]+0.0001)^2)))*((u1[j]^2)*u2[j] - u4[j]); y2[j] = (1/(2*((u1[j]+0.0001)^2)))*(u4[j] + u[j]*u1[j]*u3[j]); end for; public Modelica.Blocks.Interfaces.RealInput u[worldModel1.n + 1]; Modelica.Blocks.Interfaces.RealInput u1[worldModel1.n + 1]; Modelica.Blocks.Interfaces.RealInput u2[worldModel1.n + 1]; Modelica.Blocks.Interfaces.RealInput u3[worldModel1.n + 1]; Modelica.Blocks.Interfaces.RealInput u4[worldModel1.n + 1]; Modelica.Blocks.Interfaces.RealOutput y[worldModel1.n + 1]; Modelica.Blocks.Interfaces.RealOutput y1[worldModel1.n + 1]; Modelica.Blocks.Interfaces.RealOutput y2[worldModel1.n + 1]; end a;