CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Strange temperatur flow in buoyantBoussinesqSimpleFoam with viscous heating

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 30, 2012, 09:11
Default Strange temperatur flow in buoyantBoussinesqSimpleFoam with viscous heating
  #1
New Member
 
Join Date: Dec 2011
Posts: 12
Rep Power: 5
SiCaRiUs is on a distinguished road
Hi all,

i need your help, please.

I have added viscous heating to buoyantBoussinesqSimpleFoam. Everythink works fine, but now i have a strange heat flow in my testcase.

I have a modell, similar to the driven cavity case, but with the upper and the lower wall running against each other. Both moving walls have fixed temperatur value of 300K. just to get a steady state solution. The internal field is also at 300K.

here is what i get from my solver, slice through the cavity and a plot over line from top to bottom:




Here is my code for my TEqn.H:

Code:
{
    kappat = turbulence->nut()/Prt;
    kappat.correctBoundaryConditions();

    volScalarField kappaEff("kappaEff", turbulence->nu()/Pr + kappat);

    fvScalarMatrix TEqn
    (
        fvm::div(phi, T)
      - fvm::Sp(fvc::div(phi), T)
      - fvm::laplacian(kappaEff, T)
      - ((U & fvc::grad(p)) / Cp)
      - (((turbulence->nuEff() / Cp) * (fvc::grad(U) + fvc::grad(U)().T())) && fvc::grad(U)) 
    );

    TEqn.relax();
    TEqn.solve();

    rhok = 1.0 - beta*(T - TRef);
}
Why is there a bend in temperatur. I would have expected a parabola between top and bottom.

I hope you can help me. thanks a lot!

Greetings

Sebastian
SiCaRiUs is offline   Reply With Quote

Old   June 30, 2012, 12:15
Default
  #2
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,124
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
you add a "viscous heating" term, so it maybe not unphysical at all, because your mechanical energy turn into heat by that term
as "grad U" is higher near moving wall then the heat is higher there
nimasam 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
flow over a cylinder urgent! kevin FLUENT 8 August 11, 2015 13:00
Simulation of Steady-State Flow in a Constricted Tube with Laminar Viscous Model rezatabe FLUENT 2 September 9, 2011 12:49
Viscous Heating in a single zone??!! pipin FLUENT 0 July 14, 2008 12:08
Viscous heating Pipin FLUENT 3 July 14, 2008 12:05
Strange Solution for a simple pipe flow!! shekharc Main CFD Forum 4 May 9, 2005 09:21


All times are GMT -4. The time now is 20:50.