
[Sponsors] 
April 15, 2008, 21:20 
Hi Patrick,
Yes, it is a bi

#21 
New Member
Richard Jones
Join Date: Mar 2009
Location: Adelaide, South Australia, Australia
Posts: 22
Rep Power: 10 
Hi Patrick,
Yes, it is a bit of an overkill, but I will be modifying the trailing edge to have wiggles in it's profile (both chord and span wise)  I'm new to panel methods, but according to Martin Hepperle they aren't appropriate for very wavy surfaces, so for now I'm soldiering on with fvm. I'll give the zero velocity patch a try and see how I go. Thanks Pat Richard 

April 16, 2008, 23:55 
I modified the two TE cells an

#22 
New Member
Richard Jones
Join Date: Mar 2009
Location: Adelaide, South Australia, Australia
Posts: 22
Rep Power: 10 
I modified the two TE cells and it just makes a zero velocity point at the tip of the TE around which is a colorful wagonwheel of velocity!
I had more of a look at panel methods, and I think I was wrong (or my interpretation of Martin Hepperle's words was wrong). I'll try it and if I run into snags with my wavy profile I will increase the number of panels and go from there. If anyone else has suggestions for making the FVM potential solution on an airfoil work, let me know, I'm willing to experiment. Cheers, Richard. 

April 17, 2008, 03:20 
You could try to split the mes

#23 
Super Moderator
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 18 
You could try to split the mesh along faces which more or less follow the expected streamlines (using createBaffles). You'll have to program a bit to construct a faceSet with all those faces in it.
Or just allow some more physics and run simpleFoam ... 

April 17, 2008, 04:23 
oops,
forgot that the governi

#24 
Member
Patrick Bourdin
Join Date: Mar 2009
Posts: 40
Rep Power: 10 
oops,
forgot that the governing equation in potentialFoam is the pressure equation. So the zerovelocity boundary condition at the TE, should be indirectly enforced through a pressure boundary condition. Try to enforce at the TE panels (on top of the zerovelocity bc) a pressure equal to the total pressure in the freestream (roughly 0.5*U_inlet^2 + p_outlet  rho is missing cause p is normalized by rho in incompressible solvers). Otherwise, as Mattijs said, try simpleFoam. Use a zero viscosity, and slip bc at the airfoil surface to simulate an inviscid flow, along with a bounded convection scheme (the amount of artificial viscosity produced by these schemes should enable the Kutta condition 'automagically' like in a real viscous flow, and provide a stabilizing effect on the computation  hopefully) As regards the panel method, any higher order one (quadratic doublet or linear vortex strengths) should do the trick (with an appropriate panel distribution around the bends, and provided the bends on your profile are more of the rounded types than of the sawtoothed type). At some point, you should however give in to viscous flow simulations, because the waviness of your airfoil surface may have a significant effect on the boundary layer behavior. 

April 17, 2008, 05:08 
Thanks guys,
Yes they'll de

#25 
New Member
Richard Jones
Join Date: Mar 2009
Location: Adelaide, South Australia, Australia
Posts: 22
Rep Power: 10 
Thanks guys,
Yes they'll definitely be the rounded type bends  and yes I'm doing viscous solutions too (actually I'm comparing them). What I'm aiming for is cheap (potential + integral BL) and expensive (rans) methods of finding the boundary layer properties at the TE  and hopefully have both in OpenFOAM. I just tried the pressure condition and I get rather high velocities (1290m/s!) around the TE. I'll try simpleFoam to see if it works, but I would assume the run time would be in the same order of magnitude as the ones with viscosity, which makes the whole exercise useless for me. Perhaps OpenFOAM would be good at inverting a matrix from a panel method formulation.. 

October 17, 2008, 03:27 
Hi,
I am new to the movin

#26 
Member
John Wang
Join Date: Mar 2009
Location: Singapore
Posts: 72
Rep Power: 10 
Hi,
I am new to the moving mesh option in OpenFoam, and I'm wondering how to setup the boundary condition for the "tail" motion? I have managed to use splitmesh to generate the tail as internal surface, but I am lost in setting up the boundary condition so that the root of the tail could not move. Can someone help me on this topic? Thanks 

February 23, 2009, 05:56 
Hey..Guys.
I am a new user of

#27 
Member
Join Date: Mar 2009
Posts: 39
Rep Power: 10 
Hey..Guys.
I am a new user of OF. I want to simulate subsonic flow over NACA0012.I am using O topology.I want to learn, to make mesh such that according to angle of attack, boundary should be change.I have not much idea about it. Pls. guide me. Give me some reference. 

April 23, 2009, 04:05 

#28  
Member
John Wang
Join Date: Mar 2009
Location: Singapore
Posts: 72
Rep Power: 10 
Quote:


July 21, 2009, 02:38 
Solver fails

#29 
Member
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 10 
Hi mesh movers!
First of all thanks for providing the tools to the responsible humans! The tetFem solver (cellDecomp.) is quite slow but that doesn't matter. My problem now is the following: The solver works fine and solves with 300  700 iterations; tolerance = 1e9. Mesh stays fine. UNTIL the INITIAL residual falls to ~ 1e8, only 1020 iterations are done and the mesh gets destroyed beautifully... 1. Can this be caused by skew faces (I've got 50 out of 500.000; Max skewness 77) ??? 2. Or due to a too high timestep ??? 3. Is there a possibility to define a minimum number of iteration like for other solvers ??? Sorry: How to do it? BTW: I'm using 1.4.1dev. Thanks a lot! Marco 

November 18, 2013, 08:06 
Moving mesh with boundary layer

#30 
New Member
Join Date: Oct 2013
Posts: 2
Rep Power: 0 
Hi guys, i am simulating FSI on a zerothickness membrane. But the mesh motion solver will destroy the boudary layer, please see my pictures. What i have figured out is that, the problem happens due to the laplacian equation is solved on the cell displacement instead of point displacement, and when doing volumnpoint interpolation it gives inverted cells on the boundary layer.
From google i read the vortexbased mesh motion is only in the OFdev, but in standard OF there is only fv based mesh motion. Therefore, i can not use Mr. Jasak's motion solver, although it looks able to solve my problem. Any ideas? before inverted: beforeInverting.png after inverted: afterInverting.png 

April 14, 2015, 23:10 

#31 
Member
Hua
Join Date: May 2012
Posts: 30
Rep Power: 7 
Dear Foamer,
Now,I want to do FSI with openFOAM to solve fluid and other third FEM open source codes to solve solid.But I am not sure how to invoke FEM solver in openFOAM? Since it seems that you have coupled a FEM code (Feap) and OpenFOAM to solve FSI. Any suggestions will be appreciated. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Moving Mesh Problem  Rashad  FLUENT  0  August 28, 2006 04:31 
moving mesh problem  eos  Siemens  4  January 27, 2006 21:42 
moving mesh problem  walid  Siemens  6  April 1, 2005 17:31 
Moving mesh problem  Samir  Siemens  0  November 10, 2004 14:35 
Moving mesh problem  Jan Majer  Siemens  2  June 10, 2003 04:19 