Solving system of equations: long evaluation time with no results

Please help in solving the following system of equations. Working on non-Newtonian and Newtonian immiscible fluids (fluid mechanics)

immnnewt 2.nb (77.1 KB)

It is a tough system. I would do away with InverseFunction this way, but it is stille hard:

oldSystem = eq1 && eq2 && eq3 && eq4;
invFuncInstances = Union[Cases[oldSystem,
     InverseFunction[_][_], All, Heads -> True]] /.
   InverseFunction[f_][x_] :> {InverseFunction[f][x] ->
      If[x === B1, g1, g2],
     f[If[x === B1, g1, g2]] == x};
{substs, newEqs} = Transpose[invFuncInstances];
newSystem = (oldSystem /. substs) && And @@ newEqs // FullSimplify
sol0 = Solve[newSystem[[{2, 4}]], {A, B}][[1]]
newSystem0 = FullSimplify[Rest@newSystem /. sol0]
FindRoot[newSystem0, {{g1, -3}, {A1, -1}, {B1, I}},
 WorkingPrecision -> 50]

Thank you for your efforts, but I couldn’t solve my set of equations using the set of commands that you sent.
I have attached a set of differential equations with boundary conditions. Is it possible to solve analytically? using Mathematica. Please suggest me the method.

Immiscible fluidnonnewtonian.nb (2.72 KB)

It is very difficult to interpret the solution formula, because there are so many ambiguities: inverse functions and noninteger powers, which have more branches in the complex plane.

Thank you for your response. Ys I tried this but it did not work.

Ok. Thank you