Given a positive integer $n$, a finite field $\mathbb{F}_q$ of $q$ elements ($q$ odd), and a non-degenerate symmetric bilinear form $B$ on $\mathbb{F}_q^n$, we determine the largest possible cardinality of pairwise $B$-orthogonal subsets $\mathcal{E} \subseteq \mathbb{F}_q^n$, that is, for any two vectors $\mathbf{x}, \mathbf{y} \in \mathcal{E}$, one has $B (\mathbf{x}, \mathbf{y}) = 0$.