Using the $H^{\infty }$-functional calculus for quaternionic operators, we show how to generate the fractional powers of some densely defined differential quaternionic operators of order $m\geq 1$, acting on the right linear quaternionic Hilbert space $L^{2}(\Omega,\mathbb {C}\otimes \mathbb {H})$. The operators that we consider are of the type

\begin{align*} T=i^{m-1}\left(a_1(x) e_1\partial_{x_1}^{m}+a_2(x) e_2\partial_{x_2}^{m}+a_3(x) e_3\partial_{x_3}^{m}\right), \quad x=(x_1,\, x_2,\, x_3)\in \overline{\Omega}, \end{align*}

$\overline {\Omega }$

$\Omega$

$C^{1}$

$\Omega$

$\mathbb {R}^{3}$

$(R)$

$e_1,\, e_2,\, e_3\in \mathbb {H}$

$a_1,\,a_2,\, a_3: \overline {\Omega } \subset \mathbb {R}^{3}\to \mathbb {R}$

$T$

$T$

$T$

$P_{\alpha }(T)$

$\alpha \in (0,1)$

$T$

$T_l:=a_l\partial _{x_l}^{m}$

Steady-state diffusion in long axisymmetric structures is considered. The goal is to assess one-dimensional approximations by comparing them with axisymmetric eigenfunction expansions. Two problems are considered in detail: a finite tube with one end that is partly absorbing and partly reflecting; and two finite coaxial tubes with different cross-sectional radii joined together abruptly. Both problems may be modelled using effective boundary conditions, containing a parameter known as the trapping rate. We show that trapping rates depend on the lengths of the finite tubes (and that they decay slowly as these lengths increase) and we show how trapping rates are related to blockage coefficients, which are well known in the context of potential flow along tubes of infinite length.

The fine description of complex fluids can be carried out by describing the evolution of each individual constituent (e.g. each particle, each macromolecule, etc.). This procedure, despite its conceptual simplicity, involves many numerical issues, the most challenging one being that related to the computing time required to update the system configuration by describing all the interactions between the different individuals. Coarse grained approaches allow alleviating the just referred issue: the system is described by a distribution function providing the fraction of entities that at certain time and position have a particular conformation. Thus, mesoscale models involve many different coordinates, standard space and time, and different conformational coordinates whose number and nature depend on the particular system considered. Balance equation describing the evolution of such distribution function consists of an advection-diffusion partial differential equation defined in a high dimensional space. Standard mesh-based discretization techniques fail at solving high-dimensional models because of the curse of dimensionality. Recently the authors proposed an alternative route based on the use of separated representations. However, until now these approaches were unable to address the case of advection dominated models due to stabilization issues. In this paper this issue is revisited and efficient procedures for stabilizing the advection operators involved in the Boltzmann and Fokker-Planck equation within the PGD framework are proposed.