CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
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

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

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.

Approach

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.

Question

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

Notes

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

Reply

Tags
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 08:30
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 06:15
Low torque values on Screw Turbine Shaun Waters CFX 34 July 23, 2015 09:16
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 08:00
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55


All times are GMT -4. The time now is 01:04.