CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Circular initial condition on a cartesian mesh

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 28, 2008, 08:35
Default Circular initial condition on a cartesian mesh
  #1
jinwon park
Guest
 
Posts: n/a
I am using the smoothed heaviside function to impose the circular I.C. on a cartesian mesh. There are a little stair-case effects in simulations. Could anyone advice me to prevent unphysical stair-case effects(zig-zag) in the case where a circular I.C. is embedded on a cartesian mesh?
  Reply With Quote

Old   March 28, 2008, 11:54
Default Re: Circular initial condition on a cartesian mesh
  #2
Praveen. C
Guest
 
Posts: n/a
If you are using a finite volume method, you should set the IC using cell averages of the initial condition. THat is the correct thing to do and will give you a better representation of the IC. FOr other schemes I dont know.
  Reply With Quote

Old   March 28, 2008, 12:19
Default Re: Circular initial condition on a cartesian mesh
  #3
jinwon park
Guest
 
Posts: n/a
Could you tell me how to use cell averages of the initial condition? I am confusion of that term. Thanks in advance!
  Reply With Quote

Old   March 28, 2008, 12:30
Default Re: Circular initial condition on a cartesian mesh
  #4
Praveen. C
Guest
 
Posts: n/a
For any conserved quantity u(x,y,z) compute cell average for cell c as

U_c = Integral_over_c u(x,y,z) dx dy dz / volume(c)

and use this as initial condition in cell c. If initial condition is a simple function, then try to compute the integral analytically. Otherwise use some quadrature rule.
  Reply With Quote

Old   March 29, 2008, 09:38
Default Re: Circular initial condition on a cartesian mesh
  #5
jinwon park
Guest
 
Posts: n/a
As you mentioned as "If initial condition is a simple function, then try to compute the integral analytically", the initial condition in my cases is simple. I do not understand how to compute the integral analytically. Could you teach me how to do? Thanks for your help.
  Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SLTS+rhoPisoFoam: what is rDeltaT??? nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM 6 April 12, 2011 11:24
Low Mach number Compressible jet flow using LES ankgupta8um OpenFOAM Running, Solving & CFD 7 January 15, 2011 14:38
MRFSimpleFOAM goes divergenced! renyun0511 OpenFOAM Running, Solving & CFD 0 November 19, 2009 03:11


All times are GMT -4. The time now is 00:44.