CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Boundary conditions for vortex particle passing through a square domain

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   August 18, 2017, 08:47
Default Boundary conditions for vortex particle passing through a square domain
New Member
Join Date: Aug 2016
Posts: 19
Rep Power: 9
mtgoncalves is on a distinguished road

I am trying to simulate a single vortex particle which convects along with the freestream velocity through a free-surface square domain. Initially the particle is outside of the domain, so I simulate the particle approaching, entering, travelling through, exiting and distancing from the domain.


For this case the analytical solution for the velocity and pressure are known in closed-form. So I can compute the velocity and pressure at any given point in space and instance in time. I would like to use the analytical solution to impose boundary conditions.

Results and problem

Previously, I had good results for a case where the particle passes the domain from the outside (so it never enters or exits the domain). In this case I used the analytical solution to impose a non-uniform Dirichlet b.c. for the velocity on all four sides. This worked well, but for the case where the particle passes through the domain I get wiggles along the side at which the particle exits the domain.


What type of boundary condition(s) could I use for this case?


Analytical solution:

The velocity equals the freestream velocity plus the induction by the particle. The pressure is calculated via Bernoulli's principle.

u = uInf + Gamma/(2*pi*r^2)*(1 - exp(-0.5*(r/e)^2))*ry
v = vInf - Gamma/(2*pi*r^2)*(1 - exp(-0.5*(r/e)^2))*rx
p = pInf + 0.5*[(uInf*uInf + vInf*vInf) - (u*u + v*v)]

uInf = freestream velocity in x direction
vInf = freestream velocity in y direction
pInf = freestream pressure
Gamma = particle strength
e = particle core radius
r = Euclidean distance between particle and point of evaluation
rx = x component of r
ry = y component of r

ps: The analytical solution is NaN at r = 0. I have ensured that r never becomes 0 by positioning the particle such that it never coincides with a cell center or face center.
mtgoncalves is offline   Reply With Quote


boundary conditions, non-uniform bc, potential flow, unsteady b.c.

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 07:30
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 05:15
Low torque values on Screw Turbine Shaun Waters CFX 34 July 23, 2015 08:16
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 15:55

All times are GMT -4. The time now is 11:49.