What I am trying to solve right now is that I have a magnetic particle (nanoparticle to be exact) floating in liquid that is viscoelastic and apply sinusoidal magnetic field. The applied magnetic field will create torque that aligns the magnetic moment of the particle to the field. Also, the particle is large enough to be thermally blocked meaning that when the particle is rotated, it will physically rotate the particle not just the internal magnetic moment.
So, when the particle is rotated, it will drag along the fluid at the surface creating shear stress.
As of now I am interested in maximum force that the particle drags the fluid so I will assume that the internal magnetic moment and applied magnetic field are perpendicular to each other.
My goal is calculate shear stress so I start with calculating force. Suppose that the applied torque is in +y-axis so we have only forces in x- and z-axis direction. $\vec{F} = F_x\hat{a_x}+F_z\hat{a_z}$ and $\vec{r} = r_x\hat{a_x} + r_y\hat{a_y} + r_z\hat{a_z}$. And then I apply $\vec{\tau} = \vec{r} \times \vec{F}$ and get $\vec{\tau} = -(r_xF_z - r_zF_x)\hat{a_y}$.
Decomposing components, and calculate only on y-axis we get
$$ \tau_y = -|\vec{r}|\cos{\theta}\cos{\phi}\cdot |\vec{F}|\cos{\theta} + |\vec{r}|\sin{\theta}\cdot |\vec{F}|\sin{\theta} $$
Integrating from $\theta = -\frac{\pi}{2}$ to $\theta = \frac{\pi}{2}$ and from $\phi = 0$ to $\phi = 2\pi$,
$$\tau_y = |\vec{r}||\vec{F}|\int_{\theta = -\frac{\pi}{2}}^{\theta = \frac{\pi}{2}}\int_{\phi = 0}^{\phi = 2\pi} -\cos^2{\theta}\cos{\phi}\cdot + \sin{\theta}\cdot \sin{\theta} d\phi d\theta $$
I got this
$$\tau_y = |\vec{r}||\vec{F}|\pi^2$$
Now comes where I get stuck when I want to apply $\sigma = G\gamma$ linear relationship (in reality, this is not always the case but I want to test it first with linear assumption) where $\sigma$ is shear stress, $\gamma$ is shear strain, and $G$ is shear modulus.
From my understanding, shear stress is $\frac{F}{A}$, but the force is not constant so I have to calculate force acted on infinitsimal surface area $dA = r^2\sin{\theta}d\phi d\theta$. However, I feel something is wrong here and cannot proceed forward. I found the result $ d\sigma = \frac{F}{dA}$ very strange and wrong but I have no idea what is the better way.
I also welcome other solution or equations with more complicated or rigorous approach if there is any.