
[Sponsors] 
July 7, 2010, 00:17 
Oscillating Cylinder in Parallel Flow: Analytical Solution

#1 
New Member
Join Date: Jul 2010
Posts: 8
Rep Power: 9 
I need a paper (lecture or class notes) showing a basic analytical solution to parallel flow past an oscillating cylinder. The oscillation can be transverse or parallel to the flow. The flow is laminar, incompressible, inviscid, and unsteady. So basically, I want to solve Laplace's equation for a moving cylinder. I believe the unsteady nature of the flow appears in the momentum equation, which for this problem reduces to an unsteady Bernoulli's equation.
Ostensibly, the potential flow solution around a stationary cylinder is quite simple and easily accessible as contained in any introductory fluids text. But I've had no luck whatsoever finding the analytical solution for the moving cylinder. I have found a bunch of papers with a numerical solution though. Feel free to send those along as well. Related question just to make sure: For the problem stated (inviscid, incompressible), I won't be solving for vortex shedding? Need to include the viscous term along with the convection term to do so? Thanks for the help. [I need it for background of my research.] 

July 7, 2010, 08:30 

#2 
Senior Member
Join Date: Jul 2009
Posts: 247
Rep Power: 12 
According to theory, there is no vortex shedding if the flow is inviscid. The flow all separates at the rear stagnation point. Conversely, there is no general theoretical solution for viscous flow over a cylinder where separation and shedding occur.


July 7, 2010, 09:33 

#3 
Super Moderator

If you assume potential flow, at any instant, you have Laplace equation for the potential with some Neumann conditions on the cylinder and Dirichlet conditions in the farfield. The potential adjusts instantaneously to the boundary motion. To find the solution, you can transform to a frame in which the cylinder is at rest at that instant, use the standard solution for steady cylinder, and then transform back to your original frame. Once you have the potential, you get the pressure from the Bernoulli equation.
If you get the solution, please write about it here. 

July 7, 2010, 17:05 

#4  
New Member
Join Date: Jul 2010
Posts: 8
Rep Power: 9 
Quote:
dphi/dt + p/pho + 1/2V^2 + gz = constant So you can't just solve for the potential at each time step because of the term dphi/dt. The change in the potential from the previous time step matters for the solution. I've looked more and I still can't find an analytical solution for this problem. Anyone else have luck? Also, how exactly would I approach this numerically? If I'm using a Galerkin/FEM approach, should my test functions be a seperation of variables type function, i.e. X(x)Y(y)T(t)? The continuity equation reduces to Laplace's equation which is only a function of the 2D spatial variables. But then the unsteady Bernoulli's equation (reduction of the momentum equation) is a function of both space and time. I've only previously solved steady heat conduction and stress strain problems that included just one equation to satisfy. A little unsure about taking it to 2 equations and temporal discretization in an FEM context (I've done finite difference but that was rather trivial). 

July 8, 2010, 00:32 

#5 
Super Moderator


July 8, 2010, 12:14 

#6  
New Member
Join Date: Jul 2010
Posts: 8
Rep Power: 9 
Quote:
This unsteady "flow" is governed by the derivative with respect to time. By doing it your way (solving the Laplace equation for the instantaneous steady state according to the current BC's), you can't know the dphi/dt term. You've ignored everything in the past and you can't do that when a time derivative is present. The two equations (Laplace and unsteady Bernoulli) must be solved simultaneously. 

July 8, 2010, 14:50 

#7 
Senior Member
Join Date: Jul 2009
Posts: 247
Rep Power: 12 
Praveen is correct  the Laplace equation can be solved for the velocity potential, and the only time dependence enters in through the boundary conditions. Bernoulli's equation is then used to compute the pressure once you have the potential. The first thing I would try is to simply separate the variables  assume phi(x,t) = f(x)g(t). Then clearly Lap(f) = 0, and if the boundary conditions allow separation, you should be able to solve this equation for f, with g(t) being determined by the timedependent part of the BC. Then use Bernoulli's equation to find the pressure field p(x,t) for any value of time.
Note that this is not a formal solution procedure, but simply a suggestion of one possible approach. Also note that this does not ignore the fact that phi depends on time, but makes use of the fact that the velocity potential can only have a certain form of dependence on time based on the fact that it satisfies Lap(phi) = 0 (no time derivative in this equation). Last edited by agd; July 8, 2010 at 15:05. 

March 23, 2012, 06:49 
unsteady bernoulli equation/ natural convection open channel

#8 
New Member
Join Date: Dec 2010
Posts: 3
Rep Power: 8 
Hi every body,
I'm studying natural convection in vertical open channel (only two vertical walls without extended volumes). Unsteady NavierStokes is used to describe the flow inside the channel. For pressure boundary conditions, I'm using steady Bernoulli equation to describe the drop of pressure between freestream and the channel inlet. When using steady bernoulli I got every time a steady state motion which is not very realistic. Now, I'm trying to use unsteady instead of steady bernoulli equation. There is actually an additional term in Unsteady Bernoulli equation which is \int_{1}^{2} \frac{\partial V}{\partial t}ds. As I don't have any information about the variation of velocity along a streamline between freestream (point 1) and the channel inlet (point 2), I would like to know if there is possibility to model this additional term. Thank you 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Incorrect Drag and Drag Coefficient for flow over a cylinder  ozzythewise  Main CFD Forum  8  June 13, 2012 06:24 
bubbly flow and parallel processing  mvee  FLUENT  0  September 12, 2007 05:08 
IcoFoam parallel woes  msrinath80  OpenFOAM Running, Solving & CFD  9  July 22, 2007 02:58 
Analytical Solution for a Free Surface Flow  Ferreira, VG  Main CFD Forum  0  February 25, 2007 05:56 
Turbulent steady flow around a circular cylinder  Mirek Kabacinski  FLUENT  0  July 23, 2003 18:40 