# Error near boundary layer as 1st grid decrease

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

 LinkBack Thread Tools Display Modes
 October 3, 2007, 00:32 Error near boundary layer as 1st grid decrease #1 CF Guest   Posts: n/a Hello, I'm simulating a flow past cylinder using an O grid, in FVM. I normally intialize solution everywhere to u=1. and do the simulation. As time step increases, the velocity near the boundary will decrease since u=v=0 on its surface. It works ok when 1st grid size is at 5.e-3. When I refine my grid further to 1.e-3, 2 strange things happen to the grids near the boundary. The velocity decreases as it gets nearer the boundary but for the grid cells just next to the boundary, velocity suddenly gets larger and points in the opposite direction. Hence the resulting drag is negative. Interestingly the solution did not explode and other parts of the flow looks fine (separation etc). At 5.e-3, there's no such problem. I also tried to initialize the solution such that near the boundary, the velocity is tangential to it and it gets smaller as it get nearer but the problem is still there. Moreover, I've used the code on a C grid simulating flow past an airfoil and there's no such problem even as grid refines. Hope someone can suggest some ways to solve this problem. Thank you.

 October 3, 2007, 07:04 Re: Error near boundary layer as 1st grid decrease #2 Jake Guest   Posts: n/a What precision are you using? Although 1e-3 doesn't seem too small, it might be worthwhile to try double precision to see the effect.

 October 3, 2007, 09:37 Re: Error near boundary layer as 1st grid decrease #3 CF Guest   Posts: n/a I'm using double precision already. Btw, I wonder if it is because for initial velocity field is not solenoid. I'm not very sure how to do it. Can someone enlighten me? Thanks.

 October 3, 2007, 10:21 Re: Error near boundary layer as 1st grid decrease #4 Jake Guest   Posts: n/a Why not just use a uniform flow field as the initial condition? This is divergence free.

 Thread Tools 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 On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Hoshang Garda FLUENT 1 November 27, 2013 11:24 sina Main CFD Forum 0 February 3, 2011 10:26 ivan_cozza OpenFOAM Native Meshers: snappyHexMesh and Others 0 October 6, 2010 13:47 scottneh STAR-CCM+ 3 September 30, 2010 14:21 Mark CFX 6 November 15, 2004 16:55

All times are GMT -4. The time now is 16:57.

 Contact Us - CFD Online - Privacy Statement - Top