Anemometer

Description

An anemometer, a flow sensing device, consists of a heater and temperature sensors before and after the heater, placed either directly in the flow or in its vicinity. They are located on a membrane to minimize heat dissipation through the structure. Without any flow, the heat dissipates symmetrically into the fluid. This symmetry is disturbed if a flow is applied to the fluid, which leads to a convection on the temperature field and therefore to a difference between the temperature sensors (see Fig.1 below) from which the fluid velocity can be determined.

The physical model can be expressed by the convection-diffusion partial differential equation[4]:

$\rho c \frac{\partial T}{\partial t} = \nabla \cdot (\kappa \nabla T ) - \rho c v \nabla T + \dot q,$

where $\rho$ denotes the mass density, $c$ is the specific heat, $\kappa$ is the thermal conductivity, $v$ is the fluid velocity, $T$ is the temperature and $\dot q$ the heat flow into the system caused by the heater.

The solid model has been generated and meshed in ANSYS. Triangular PLANE55 elements have been used for meshing and discretizing by the finite element method. The order of the system is $n = 29008$.

Example with 1 parameter:

The $n$ dimensional ODE system has the following transfer function $G(p) = c((sE - A_{v0}- p(A_{v1} - A_{v0}))^{-1}b)$

with the fluid velocity $p(=v) \in [0, 1]$ as single parameter. Here $E$ is the heat capacitance matrix, and Failed to parse (unknown function "\jmat"): \jmat M_{d} incorporates the heat dissipation and thus contains the convection and conduction terms. $b$ is the load vector which is derived from separating the spatial and temporal variables in $\dot{q}$ and the FEM discretization w.r.t. the spatial variables. $A$ are the stiffness matrices with $\_v0$ for pure diffusion, $\_v1$ diffusion and convection. Thus, for obtaining pure convection you have to compute $A_{v1} - A_{v0}$.