I’m working with a channel geometry that has periodic variations. My goal is to solve the Navier–Stokes equations to compute the velocity field, and then use this field to advect a passive scalar.
Currently, I compute the flow over the entire domain and then use it to advect the scalar field, but this approach is inefficient. Since the flow field is periodic, but the passive scalar is not, a more efficient strategy would be to solve for the flow in a single periodic cell and then tile this solution across the full domain. This way, the velocity field can be extended to the whole domain efficiently and used for scalar advection.
I’m wondering: how can this kind of periodic flow solution and domain extension be implemented in FreeFEM? Thanks for your time and suggestions in advance.
The attached file is my current implementation.
full curved channelS.edp (3.9 KB)