Let $\Omega$ and open subset of $\mathbb{R}^n$. Let us consider the following operator: $$ \Delta_A (u)\, \, \colon= \text{div}(A \nabla u ), \qquad u \in C^{\infty}(\Omega) $$ where $A(x)$ is a matrix with smooth coefficients and uniformly elliptic on $\mathbb{R}^n$, i.e. $A$ is symmetric and there are positive constants $\lambda$ and $\Lambda$ such that for all $x \in \mathbb{R}^n$ and for all vectors $v, w \in \mathbb{R}^n$: $$ \langle v, A(x) w\rangle \le \Lambda |v||w| $$ $$ \langle v, A(x) v\rangle \ge \lambda |v|^2 $$

Consider the distance function from the origin $\rho(x) := dist(x, 0)$.

Is it possible to have an estimate of the kind:

$$ \Delta_A(\rho) \le \frac{K_1}{\rho} + K_2 \qquad \text{on } \Omega = \mathbb{R}^n \setminus\{0\} $$ where $K_1, K_2$ are constants possibly depending on $n$ and on $\lambda$ and $\Lambda$?