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

Simulation of vertical Bridgman,enthalpy-porosity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 24, 2004, 04:20
Default Simulation of vertical Bridgman,enthalpy-porosity
  #1
cloud0574
Guest
 
Posts: n/a
A question about Numerical simulation of vertical Bridgman crystal growth using enthalpy-porosity method I¡¯m writing my code to simulate the CdTe growth process, and yesterday I encountered a hard nut to crack. I employ a traditional method, enthalpy-porosity formulation with fixed-grid, which has been implemented in Fluent. From several prior papers, I knew that the increasing ampoule pulling rates have the impact on the solid/liquid interface shape. But, I have no idea how to incorporate the pulling rate V into the governing equations. The computation is assumed in the steady state where the pulling rate equals to the crystal growth rate. Please give me a hint or some valuable references. Thanks. BTW, my technique comes from Pro. Shyy Wei¡¯s books and papers.
  Reply With Quote

Old   September 24, 2004, 04:39
Default Re: Simulation of vertical Bridgman,enthalpy-poros
  #2
Santhanu
Guest
 
Posts: n/a
My first question is are u modelling a pure substance or a binary substance modelling.

You cannot implement stefan condition along with pulling rate in in Enthalpy porosity technique. Bcos it relaxes some dT over few grid cells and the result is affected in the source term of the Energy equation. In short the interface is not sharp (in numeric sense it does not pass through single cell) to implement Boundary conditions on the interface.

If u are using pure substance then there is a solution. You have to ALE(Arbitrary Eulerian and Lagrangian Techniques) where you can explicitly track the interface using lagrangian technique and solve it in Eulerian fashion.

You may also look up 'Journal of Crystal growth'
  Reply With Quote

Old   September 24, 2004, 08:43
Default Re: Simulation of vertical Bridgman,enthalpy-poros
  #3
cloud0574
Guest
 
Posts: n/a
Thanks for your reply , Santhanu. My concern is an isotherm phase change problem about a pure substance, CdTe.. You said, "You cannot implement stefan condition along with pulling rate in Enthalpy porosity technique." I found some papers using Enthalpy porosity technique to simulate CdTe vertical Bridgman always give a pulling rate. I think it must be an influencing factor. In fact, I want to reveal the pulling rate's impact on the interface shape, however, as you said, the pulling rate generally is implemented by explicitly tracking the interface. I'm afraid I have to change my method.

However, in Shyy Wei¡¯s Computational Fluid Dynamics With Moving Boundaries Taylor & Francis (1996) PP164

¡°Source term in momentum equation is A(u-u(cast)), A=C(1-f^2)/(f^3+q). u(cast) denote the prescribed velocity at which the solid phase is moving, say, for example, in a crystal growth or a casing process.¡±

It appends a u(cast) while the source term in energy equation is common latent heat update. I have ever thought u(cast) is what I want, but the result is disappointing.
  Reply With Quote

Old   September 28, 2004, 07:05
Default Re: Simulation of vertical Bridgman,enthalpy-poros
  #4
Santhanu
Guest
 
Posts: n/a
Hi,

In the momentum equation in the Enthalpy porosity technique uses a darcy's type law to suppress velocities in the solidified part. and a Ucast is supplied there to take care of the pulling velocity(to take care of the relative movement). It can be a way to do it. It is an easier way of treating the interface. But if u are able to get your results closer to experiments then everything is perfect.

I have never done any Vertical bridgeman modelling but as far as i know it is the ampoule which is given a velocity but not the solidified itself as in Czochralski Process
  Reply With Quote

Old   November 5, 2004, 04:41
Default Re: Simulation of vertical Bridgman,enthalpy-poros
  #5
Madhukar M. Rao
Guest
 
Posts: n/a
It is also possible that thermophysical properties such as thermal conductivity of solid vs. the melt, as well as flow conditions and pulling rates affect the interface shape. We did those calculations many years ago, but, I remember thermophysical properties having a significant impact on the interface shape.

Are you trying to compare the computations with measurements?

Best wishes,

Madhukar M. Rao
  Reply With Quote

Old   December 12, 2004, 02:40
Default Re: Simulation of vertical Bridgman,enthalpy-poros
  #6
nouri sabrina
Guest
 
Posts: n/a
Nouri sabrina Attached of research to the unit of search for develop technological silicon in Algerie. Adress: udts 02 Bd frantz fanon, Bp 399 Alger-gare (algerie). Fax: 213/21/43/35/11 Télé: 213/21/27/98/80 research topic:Numerical analysis of the transfer of heat and momentum in a vertical system of growth of Bridgman

I do not have a clear idea in connection with this subject, which relates to: equations which govern the physical phenomenon, the boundary conditions in a vertical growth which occurs in the furnace of Bridgman. I ask you for illustrations is by e-mail or it is possible some articles. I do not have a clear idea on the method of enthalpy

  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
GUI crash and simulation engine still running RPJones FLOW-3D 2 November 9, 2010 08:18
velocity profile export from a simulation onto another sudhirlv STAR-CCM+ 1 September 12, 2010 18:57
boundary conditions for simulation of vertical channel gharek Main CFD Forum 3 April 12, 2010 01:59
strange simulation error Ralf Schmidt FLUENT 2 May 4, 2007 13:02
3-D Contaminant Dispersal Simulation Apple L S Chan Main CFD Forum 1 December 23, 1998 10:06


All times are GMT -4. The time now is 12:55.