# Boundary conditions for vortex particle passing through a square domain

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 August 18, 2017, 08:47 Boundary conditions for vortex particle passing through a square domain #1 New Member   Join Date: Aug 2016 Posts: 19 Rep Power: 8 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.

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

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 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 OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post mpeppels CFX 11 August 22, 2019 07:30 MichaelK CFX 12 September 1, 2016 05:15 Shaun Waters CFX 34 July 23, 2015 08:16 sunilpatil CFX 8 April 26, 2013 07:00 Mark CFX 6 November 15, 2004 15:55

All times are GMT -4. The time now is 22:52.

 Contact Us - CFD Online - Privacy Statement - Top