TELEMAC-MASCARET Forum: Initialization of k-epsilon model (1/1)

Welcome, Guest

TOPIC: Initialization of k-epsilon model

Initialization of k-epsilon model 12 years 5 months ago #4689

qilong
OFFLINE
Expert Boarder
Posts: 340
Thank you received: 33

Hello,

I'm looking at the subroutines of k-epsilon model. In the subroutine KEPINI.f, I found a code block which was commented.It's the initialization of k and epsilon according to a theory by Vincent Boyer modified by Martin Ferrand.

! A THEORY BY VINCENT BOYER MODIFIED BY MARTIN FERRAND
!
! DO IPOIN2 = 1,NPOIN2
! DO IPLAN = 1,NPLAN
!
! ARBITRARY COMPUTATION OF K EXPRESSED AS A PERCENTAGE OF SPEED
!
! AK(IPOIN2,IPLAN) = 1.D-3*U(IPOIN2,IPLAN)**2
! AK(IPOIN2,IPLAN) = MAX(AK(IPOIN2,IPLAN),KMIN)
!
! COMPUTATION OF EPSILON
!
! EP INITIALISED ACCORDING TO UETOIL**3/KAPPA/Y
! WHERE UETOIL IS CALCULATED FROM THE VALUE OF K AT THE WALL
!
! IF(IPLAN.EQ.1) THEN
! DIST = (Z(IPOIN2,2)-ZF(IPOIN2))/FICTIFEPS
! ELSE
! DIST = Z(IPOIN2,IPLAN)-ZF(IPOIN2)
! ENDIF
! EP(IPOIN2,IPLAN)=CMU**0.75*SQRT(AK(IPOIN2,1)**3)/KARMAN/DIST
! EP(IPOIN2,IPLAN)=MAX(EP(IPOIN2,IPLAN),EMIN)
! ENDDO
! ENDDO

Why not using this theory to initialize the k and epsilon? Thanks!

Best regards,
Qilong

The administrator has disabled public write access.

Re: Initialization of k-epsilon model 12 years 5 months ago #4693

jmhervouet	Hello, You can use this theory but we cannot put it as a standard one. As a matter of fact it gave too high values in the open sea in one of our applications. Other subject: I am trying to get references for RNG and YAP, but am still waiting for them. With best regards, Jean-Michel Hervouet
	The administrator has disabled public write access.

Re: Initialization of k-epsilon model 12 years 5 months ago #4695

qilong OFFLINE Expert Boarder Posts: 340 Thank you received: 33	Thank you very much for your kind help! Best regards, Qilong
	The administrator has disabled public write access.