![]() |
Sapphire++
v1.3.0-121-g2682132 |
Inflow boundary condition. More...
#include <config.h>
Public Member Functions | |
| BoundaryValueFunction (const PhysicalParameters &physical_parameters, const unsigned int system_size) | |
| Constructor. | |
| void | bc_vector_value_list (const std::vector< dealii::Point< dim > > &points, const unsigned int boundary_id, std::vector< dealii::Vector< double > > &bc_values) const |
| Values of the distribution function \( f \) at the boundary specified by boundary_id. | |
Inflow boundary condition.
| dim | Dimension of the reduced phase space \( (\mathbf{x},p) \) |
|
inline |
Constructor.
| physical_parameters | User defined runtime parameters |
| system_size | Number of expansion coefficients, normally \( (l_{\mathrm{max}} + 1)^2 \) |
|
inline |
Values of the distribution function \( f \) at the boundary specified by boundary_id.
Return a vector containing the values of the expansion coefficients \( f_{i(l,m,s)}(t, \mathbf{x}, p) \) representing the distribution function \( f \) at points on the boundary with boundary_id.
The zeroth component of the vector corresponds to \( f_{000} \), the first component to \( f_{110} \) etc. The mapping between the system index \( i \) and the spherical harmonic indices \( (l,m,s) \) is given by lms_indices:
The points in reduced phase space are handed over using the vector points. Each of its elements has dim components. The first dim_cs components correspond to the spatial coordinates \( \mathbf{x} \) and if the momentum term is activated the last component to the momentum coordinate \( \ln p \):
where points[0] picks out the first point in points.
The boundary_id is used to specify the values of \( f \) at a specific boundary of the computational domain.The relation between the possible boundaries and the boundary ids depends on the dimension dim. For a 1D spatial problem with momentum the dependence, the boundary ids are,
| Boundary | ID |
|---|---|
| lower x | 0 |
| upper x | 1 |
| lower p | 2 |
| upper p | 3 |
Only if in the parameter file the corresponding boundary is set to inflow, the function bc_vector_value will be called.
| points | Points in reduced phase space on the boundary |
| boundary_id | ID of the boundary of the computational domain |
| bc_values | Return vector \( f_{i(l,m,s)}(t, \mathbf{x}, p) \) |
[Boundary value]
[Boundary value]