|
[Sponsors] |
August 21, 2000, 11:25 |
Mooving Grids
|
#1 |
Guest
Posts: n/a
|
Hi, I implemented mooving grids in my flow solver and chose as test case a circular cylinder with a constant inflowing velocity field. The grid is circular and symmetric. The wall does not moove. I supposed that the result is the same if I let the grid rotate or not, in particular I expected a symmetric flow field. However, if the grid mooves, the flow field is not symmetric any more. The asymmetry is the stronger the larger is the grid speed. I did a benchmark against TascFlow and found a similar result. These calculation were steady state. If I perform transient calculations, this effect occures if the time step is rather large, e.g. if I resolve one grid revolution with less than 10 time steps and the grid velocity is of the same order as the inflow velocity, irrespective if I use the space conservation law or calculate the grid velocities by pure kinematics. What might be the reason, and what other test cases can be recommended? In the following some implementation details: - Absolute, cartesian velocity components as dependent variables - Co-located cell-centred variable arrangement - The only difference to fix grids is the treatment of the convective fluxes, namely a "relative mass flux" for the convective transport is introduced by the grid movement - The Conti- or pressure correction equation is kept unchanged,except that the divergence is evaluated by the "relative mass flux"; should have no influence for pure rotation of the circular grid. Thank you very much for any hint!
With best regards from Munich Romuald Skoda |
|
August 21, 2000, 11:47 |
Re: Mooving Grids
|
#2 |
Guest
Posts: n/a
|
(1). If you rotate the grid at a constant speed, then it becomes very difficult to compute the flow using the fixed frame of reference. (2). If you sit on the rotating frome and write the governing equations on the rotating frame of reference, then the grid will be stationary relative to it. (3). The problem is: rotation is not translation, constant rotation will produce acceleration, while constant translation will not produce acceleration. As a result, the governing equations will be different in the rotating frame of reference. (4). You also need to take care of the relative motion of the boundary conditions, if you are using the rotating frame of reference. (5). If on the other hand, you are using fixed reference frame, then the moving grid is basically a transient flow problem, there will be contributions from the motion of the grid. So, you should be very careful in deriving the governing equations.
|
|
August 22, 2000, 00:36 |
Re: Mooving Grids
|
#3 |
Guest
Posts: n/a
|
is your flow laminar or turbelent
|
|
August 22, 2000, 00:55 |
Re: Moving Grids
|
#4 |
Guest
Posts: n/a
|
The test case is laminar (Re=1). In general, I am going to use wall functions later. In response to Dr. Chien: what in particular is the problem about staying in the absolute frame and let the grid move? I chose the absolute frame because I am going to couple the moving domain with a fixed one. By the way: Sorry for the ugly spelling error in "moving"
|
|
August 22, 2000, 06:18 |
Re: Moving Grids
|
#5 |
Guest
Posts: n/a
|
The best way to test your implmentation of moving grid is to have no physical flow - in your case use zero inflow velocity. Even better is to use inviscid so you only have convection term to worry about since the grid velocity only affects the convection.
Now use the same case, the mesh movement itself should not create any flow at all (in reality due to numerical precision, you may have velocity in the order of 10E-5 or so). |
|
August 22, 2000, 07:44 |
Re: Moving Grids
|
#6 |
Guest
Posts: n/a
|
I would use the following simple test case: Set up a regular structured grid, let the outer boundaries fixed and then just contract the inner mesh. E.g. x and y components do decrease and then increase again, such that one gets a volume change. (Peric test case) If the fluid is at rest, nothing should happen. If something happens there is something wrong with the implementation of the GCL, I would guess.
Formulating both, the volumes of the cells and grid velocities on geometric quantities of course does give wrong results due to violation of the GCL. (I assume that this is almost not worth mentioning anymore) I would first look for the implementation of the GCL, if that is ok. Frank |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] OpenFOAM - structured or unstructured Grids ? | thomasduerr | OpenFOAM Meshing & Mesh Conversion | 22 | November 24, 2018 23:18 |
A question about dynamic overset grids | JXIA | Main CFD Forum | 4 | January 10, 2008 03:50 |
Reloading grids in Pre | Pascale Fonteijn | CFX | 1 | July 7, 2003 13:33 |
UDF Face Fluxes & Partitioned Grids | Greg Perkins | FLUENT | 0 | April 6, 2001 01:50 |
Can I define grids motion as a time functions | Darcy | CFX | 4 | March 29, 2001 21:14 |