On the existence of a positive solution to a nonlocal logistic system with nonlinear advection terms

In this paper, we study a nonlocal logistic system with nonlinear advection terms \begin{equation*} \left{ \begin{array}{lcl} -Δu+\vecα(x)\cdot \nabla (|u|^{p-1}u)&=&\left(a-\int_ΩK_1(x,y)f(u,v)dy \right)u+bv\mbox{ in }Ω,\ -Δv+\vecβ(x)\cdot \nabla (|v|^{q-1}v)&=&\left(d-\int_ΩK_2(x,y)g(u,v)dy \right)v+cu\mbox{ in }Ω,\ \qquad \qquad \qquad \qquad u=v&=&0\mbox{ on }\partialΩ, \end{array} \right. \end{equation*} where $\mathrm{\Omega }\subset {\mathbb{R}}^N$, $N\ge 1$, is a bounded domain with a smooth boundary, $\vec{\alpha }(x)=({\alpha }_1(x),\cdots ,{\alpha }_N(x))$ and $\vec{\beta }(x)=({\beta }_1(x),\cdots ,{\beta }_N(x))$ are flows satisfying suitable conditions, $p,q\ge 1$, $a,b,c,d>0$ and $K_1,K_2:\mathrm{\Omega }\times \mathrm{\Omega }\to \mathbb{R}$ are nonnegative functions, with their specific conditions detailed below. The functions $f$ and $g$ satisfy some assumptions which allow us to use bifurcation theory to prove the existence of solution to problem $(P)$. It is important to highlight that the inclusion of the integral nonlocal term on the right-hand side makes the problem more representative of real-world situations.