[alvora]

Hi,

I think the fluid density (rho) is not at all required to solve the current continuity equation (charge conservation equation).

Let's say,
rhoE - Space charge density [C/m³ or A-s/m³],
mu - electrical mobility [A-s²/kg]
ePot - electrical potential [kg-m²/(A-s³)]
De - Diffusion coefficient [m²/s]
U - fluid velocity [m/s]

then
Current continuity equation:

fvScalarMatrix rhoEEqn
(
fvm::ddt(rhoE) //Time derivative
+ fvm::div(rhoFlux, rhoE) //Convective term - rhoFlux is surfaceScalarField
- fvm::laplacian(D, rhoE) //Diffusive Term, -/+ sign for +/- space charge density
);
rhoEEqn.solve();

where,
surfaceScalarField rhoFlux = - mu*mesh.magSf()*fvc::snGrad(ePot)

If the space charge density is also dragged by the fluid velocity, then
surfaceScalarField rhoFlux = - mu*mesh.magSf()*fvc::snGrad(ePot) + fvc::linearInterpolate(U)&mesh.Sf()


Isn't it this way? correct me if I am wrong or missing something.


Kind regards
Alpesh

[damon707]

Dear all,

I would like to do a simulation of boat propeller with AMI solid body rotating motion and I am interested in the far-wake characteristics of the propeller. Im guessing that buoyancy forces will be important for the far-wake of a propeller thus I need to use a solver which takes buoyancy into account. Since there will be very small stratification of the density in such a problem, will I be able to use buoyantboussinesqpimplefoam for this simulation? Also are there any tutorials or some papers available in this kind of propeller simulations? Are there any known problems of the buoyantboussinesqpimplefoam solver?

Kind regards,
George