
[Sponsors] 
November 27, 2007, 08:58 
Non dimensionalisation Problems  Please Advise

#1 
Guest
Posts: n/a

I am trying to solve a problem involving Marangoni convection in which an air bubble is surrounded by a liquid oil layer. The governing equations are NavierStokes and The energy Eqn. I am trying to benchmark against data in the literature but my data points are all out by a constant factor, so I am assuming I have made some error in the nondimensionalisation of the problem, or rather, in the implementation of the nondimensional problem when using a commercial CFD code.
The original eqns were nondimensionalised using the same reference scales as the author to which I am comparing results, and I get the same nondimensional form of the eqns. Original NS: (ρ_&alpha*(∂u/∂t + (u.V)u) = V[pI + (μ_&alpha*(Vu + (Vu)^t)]  (ρ_op)*(β_&alpha*g*(TT_op) Non_Dimensional NS: (ρ_α/ρ_oil)*(∂u'/∂t' + (u'.V')u') = V'[p'I +(μ_α/μ_oil)* (V'u' + (V'u')^T)]  (ρ_α/ρ_oil)*(β_α/β_oil)*(Ra/Pr)*T'z Original Energy: (ρ_&alpha*(Cp_&alpha*(∂T/∂t + u.VT) = (k_&alpha*(V^2)T Non_Dimensional Energy: (ρ_α/ρ_oil)*(Cp_α/Cp_oil)*Pr*(∂T'/∂t' + u'.V'T') = (k_α/k_oil)*(V'^2)T' The important boundary condition (BC) is at the liquidvapor interface, where the shear stress is related to surface tension effects. Dimensional BC: μ(∂u_tang/∂n) = (∂σ/∂T)*(∂T/∂s) NonDimensional BC: (∂u_tang'/∂n') = (Ma/Pr)*(∂T'/∂s') Apologies for the appearance of the eqns. The prim symbol denotes dimensionless quantities. ρ,μ,Cp,k,β refer to density, dynamic viscosity, specific heat, thermal conductivity and volume expansion coefficient respectively. u is the velocity, and V is the differential operator. g is gravitational acceleration, T is temperature, and Ra, Pr, Ma are the Rayleigh, Prandtl and Marangoni numbers. ∂σ/∂T is the tempeture derivative of surface tension. The subscript alpha refers to the reference fluid, so if the reference fluid is chosen to be the oil (as properties of liquid generally much greater in magnitude to properties of vapor/air), should all those ratios in the nondimensional form should go to unity? I have tried to construct this problem using both FLUENT and COMSOL codes, but with no success so far. In both cases, I have nondimensionlised the original geometry (0.5mm radius bubble, 10mm liquid layer depth) using the reference length scale, the radius. Hence I have a dimensionless geometry with bubble radius of 1 and liquid layer depth 20. I must set the following parameters: density, dynamic viscosity, specific heat, thermal conductivity, volume expansion coefficient, gravitational acceleration, and temperature derivative of surface tension. My problem is what values do I give these properties??? I am also confused by the appearance of the Pr term in my nondimensional energy eqn. Any help will be greatly appreciated. Thanks you. 

November 30, 2007, 16:30 
Re: Non dimensionalisation Problems  Please Advis

#2 
Guest
Posts: n/a

First, be sure whethere the commercial software solver is based on dimension or nondimensional....i hope Fluent is a dimensional software..it cannot take values nondimensonally....just check....i am also not sure ...


November 30, 2007, 19:00 
Re: Non dimensionalisation Problems  Please Advis

#3 
Guest
Posts: n/a

I think FLUENT is dimensional, but I thought it shouldn't matter if the ratios of the fluid properties and the geometry are such that the relevant dimensionless properties like Re, Pr, Ra, Ma are the same for both the dimensional and dimensionless cases?


February 28, 2011, 13:23 

#4 
New Member
lyes
Join Date: Feb 2011
Posts: 24
Rep Power: 8 
fluent trat only dimensional equations but we can use dimensionless equations and change something so that Fluent will understand them like dimensional equation when u write the dimensional and non diemnsonal u see that if we set ro=1 and viscosity + 1/Re they apear the same and like that everything is ok but results of Fluent are dmensionless


February 28, 2011, 13:29 

#5 
New Member
lyes
Join Date: Feb 2011
Posts: 24
Rep Power: 8 
hello i'm working on simulate a turbulent flow in a rectangularr channel were heat blocks are amounted i want have a Gr of 2.5 E +7 but i dont know how to set the thermal expansion coefficient and delta T in Fluent ?


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
[ICEM] Problems with coedge curves and surfaces  tommymoose  ANSYS Meshing & Geometry  0  August 5, 2011 16:02 
Needed Benchmark Problems for FSI  Mechstud  Main CFD Forum  4  July 26, 2011 12:13 
Nondimensionalisation Problems  Please Advise  Bren  Main CFD Forum  19  August 21, 2007 06:19 
Some problems with Star CD  Micha  Siemens  0  August 6, 2003 13:55 
Airpak: Problems with heat flow contour conditions  Ivan de Miguel Urain  FLUENT  1  November 16, 2001 10:36 