I am confused. I don't understand the difference a ARMA and a GARCH process.. to me there are the same no ?

Here is the (G)ARCH(p, q) process

$$\sigma_t^2 = 
    \underbrace{
    \underbrace{
     \alpha_0 
     + \sum_{i=1}^q \alpha_ir_{t-i}^2}
    _{ARCH}  
     + \sum_{i=1}^p\beta_i\sigma_{t-i}^2}
    _{GARCH}$$

And here is the ARMA($p, q$):

$$ X_t = c + \varepsilon_t +  \sum_{i=1}^p \varphi_i X_{t-i} + \sum_{i=1}^q \theta_i \varepsilon_{t-i}.\,$$

Is the ARMA simply an extension of the GARCH, GARCH being used only for returns and with the assumption $r = \sigma\varepsilon$ where $\varepsilon$ follows a strong white process?

You are conflating the features of a process with its representation. Consider the (return) process $(Y_t)_{t=0}^\infty$. 


An ARMA(p,q) model specifies the conditional mean of the process as


$$
\begin{align}
\mathbb{E}(Y_t \mid \mathcal{I}_t) &= \alpha_0 + \sum_{j=1}^p \alpha_j Y_{t-j}+ \sum_{k=1}^q \beta_k\epsilon_{t-k}\\
\end{align}
$$
Here, $\mathcal{I}_t$ is the information set at time $t$, which is the $\sigma$-algebra generated by the lagged values of the outcome process $(Y_t)$. 


The GARCH(r,s) model specifies the conditional variance of the process
$$
\begin{alignat}{2}
& \mathbb{V}(Y_t \mid \mathcal{I}_t) &{}={}& \mathbb{V}(\epsilon_t \mid \mathcal{I}_t) \\
\equiv \,& \sigma^2_t&{}={}& \delta_0 + \sum_{l=1}^r \delta_j \sigma^2_{t-l} + \sum_{m=1}^s \gamma_k \epsilon^2_{t-m}
\end{alignat}
$$


Note in particular the first equivalence $ \mathbb{V}(Y_t \mid \mathcal{I}_t)= \mathbb{V}(\epsilon_t \mid \mathcal{I}_t)$.

Aside: Based on this representation, you can write
$$
\epsilon_t \equiv \sigma_t Z_t
$$
where $Z_t$ is a strong white noise process, but this follows from the way the process is defined.


The two models (for the conditional mean and the variance) are perfectly compatible with each other, in that the mean of the process can be modeled as ARMA, and the variances as GARCH. This leads to the complete specification of an ARMA(p,q)-GARCH(r,s) model for the  process as in the following representation
$$
\begin{align}
Y_t &= \alpha_0 + \sum_{j=1}^p \alpha_j Y_{t-j} + \sum_{k=1}^q \beta_k\epsilon_{t-k} +\epsilon_t\\
 \mathbb{E}(\epsilon_t\mid \mathcal{I}_t) &=0,\, \forall t \\
\mathbb{V}(\epsilon_t \mid \mathcal{I}_t) &= \delta_0 + \sum_{l=1}^r \delta_l \sigma^2_{t-l} + \sum_{m=1}^s \gamma_m \epsilon^2_{t-m}\, \forall t
\end{align}
$$

