|
[Sponsors] |
Boundary conditions for vortex particle passing through a square domain |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 18, 2017, 09:47 |
Boundary conditions for vortex particle passing through a square domain
|
#1 |
New Member
Join Date: Aug 2016
Posts: 19
Rep Power: 10 |
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 |
Display Modes | |
|
|
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 |