JuliaIBPM / ViscousFlow.jl

A framework for simulating viscous incompressible flows about arbitrary body shapes.
https://JuliaIBPM.github.io/ViscousFlow.jl/latest
Other
100 stars 28 forks source link

Explicit RHS #34

Open mkalia94 opened 4 years ago

mkalia94 commented 4 years ago

Hi,

Firstly, thanks for a very efficient high-level setup for Navier Stokes simulations!

I'm looking at the fluid flow past a cylinder example, and need an explicit right hand side.

Am I correct in assuming that for vorticity w for a particular time t, the corresponding right hand side is given by:

rhs = 1/Re*sys.L*w .+ ConstrainedSystems.r₁(w,t,sys) .+ plan_constraints(w,t)[1](f)

where f is f = VectorData(collect(body)), body is the cylinder.

Note that in my case I do not apply the PointForce as done in the examples.

Thanks!