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

Time step help in Explicit method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 13, 2013, 13:06
Default Time step help in Explicit method
  #1
Member
 
sud
Join Date: May 2013
Posts: 65
Rep Power: 4
archeoptyrx is on a distinguished road
Hello everyone ,

I have wrote a code for solving compressible navier stokes code over a flat plate (same as given in the Anderson CFD Book) . when i tried simulating it with free stream conditions at higher altitude for example

Plate length of 1e-5 encounters slip regime above 10km altitude , which means that slip boundary conditions should be incorporated to get proper results , and the code blows up for free stream conditions above 20km (because transition regime comes into play )

Am i correct guys ? Navier stokes equation works only for continuum regime right ?

So i am using MACCORMACK Method with the time step given in book . Its because of the time step , which brings up complex numbers and its causes the program to stop .

What should i do now ?

Any suggestion will be of great help

Thanks in advance
archeoptyrx is offline   Reply With Quote

Old   August 13, 2013, 13:28
Default
  #2
Senior Member
 
flotus1's Avatar
 
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,104
Rep Power: 19
flotus1 will become famous soon enoughflotus1 will become famous soon enough
While the NS equations are only valid in the continuum regime, a NS solver should still be able to produce a (wrong) solution in the rarefied regime.
If your NS code diverges for higher Knudsen numbers, there must be something wrong with the code.
flotus1 is offline   Reply With Quote

Old   August 13, 2013, 13:31
Default
  #3
Member
 
sud
Join Date: May 2013
Posts: 65
Rep Power: 4
archeoptyrx is on a distinguished road
Alex,

I checked the code and i find that i implemented the time step (stability criteria) wrong . This produces the complex number first which makes the solver to blow up for higher Knudsen numbers
archeoptyrx is offline   Reply With Quote

Old   August 13, 2013, 13:35
Default
  #4
Member
 
sud
Join Date: May 2013
Posts: 65
Rep Power: 4
archeoptyrx is on a distinguished road
Take a look at my time step , I implemented it wrong .

Code:
function delt=Tstep(u,v,T,rho,mu)
% K         -Fudge factor 0.5<K<0.8
% delx      - xinterval
% dely      - y interval
% u         - horizontal velocity at (i,j)
% v         - Vertical velocity at (i,j)
% a         - Speed of sound at (i,j)
% mu        - Dynamic viscosity at (i,j)
% Pr        - Prandtl Number
% rho       - Density at (i,j)
% gamma     - iscentropic constant
global K Pr gamma imax jmax delx dely;
% global u v T rho mu;

a=Sound(T);
vda=zeros(imax,jmax);
deltcfl=zeros(imax,jmax);
for i=1:imax
    for j=1:jmax
%                mu(i,j)=Dynvis(T(i,j),2);
        vda(i,j)=((4*mu(i,j)*gamma*mu(i,j))/(3*rho(i,j)*Pr));
    end
end
vdas=max(max(vda(:)));
for i=1:imax
    for j=1:jmax
        A=abs(u(i,j))/delx;
        B=abs(v(i,j))/dely;
        C=a(i,j)*sqrt((1/(delx^2))+(1/(dely^2)));
        D=2*vdas*((1/(delx^2))+(1/(dely^2)));
        deltcfl(i,j)=(A+B+C+D)^(-1);
    end
end

delt=real(min(min(K*deltcfl(:))));
end
archeoptyrx is offline   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
Unstabil Simulation with chtMultiRegionFoam mbay101 OpenFOAM Running, Solving & CFD 13 December 28, 2013 14:12
mixerVesselAMI2D's mass is not balancing sharonyue OpenFOAM Running, Solving & CFD 6 June 10, 2013 09:34
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34
AMI speed performance danny123 OpenFOAM 19 October 24, 2012 07:44
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM 6 April 12, 2011 11:24


All times are GMT -4. The time now is 19:13.