packedbed geometry creation
dear all
I'm trying to do some analysis on latent heat thermal energy storage using a packed bed of spheres, i need to run the simulation in 2d but the problem is that the direct treatment as a 2D model will lead to no flow passage for the fluid due to the complex packing of the spheres and the contact between them. the bed have following dimensions: 360mm diameter 460mm height 264 spheres with outer diameter of 55mm uniformly packed in eight layer anyone can help with this? 
if you are using Fluent (I just know Fluent has a porous media model) you can treat it as porous media. You'll have to figure out some of the variables, but it can be done. There are journal articles written modeling it as porous media too.

Quote:
do you have any specific article in mind to help me with this? i also wanna treat it as a porous media but i have problem with the geometry, i don't know how i should put the spheres so that have acceptable flow passage. by the way the void fraction or the bed porosity is 0.48. 
If you are using Fluent you will not have to draw every sphere etc. You will just have to draw the dimensions of the pipe/bed/whatever that is holding the spheres and then define certain variables that will let Fluent know that the section is a packed bed. The porosity will be one of them. When I did this a few years ago I had to do some measurements in the lab to get some of the other variables. Fluent's manual has a good section on porous media that will help you understand it better.
I'll try to remember to get some author's names when I get home. Jiang comes to mind, but I am not sure if it was packed beds he was doing. 
I'm using ansys fluent 13
I need the liquid fraction of the pcm (phase change material which is paraffin wax ) and the outlet temperature after the system is being totally charged... so you sure that i don't have to draw every sphere? 
get out the fluent PDF and read the porous media section.. yes, you don not have to draw all the spheres..:)
You will separate you section into an entrance length, and then the fluidized bed, and an exit section. The bed section will be defined as porous media and you will have to fill in all the necessary variables. That is what I had to do anyway. 
Rheolef NEW version 6.1: an efficient FEM C++ finite element library for solving PDE
Rheolef: an efficient FEM C++ finite element library for solving PDE
Version : 6.1 Book: http://ljk.imag.fr/membres/Pierre.Sa...ef/rheolef.pdf Home: http://ljk.imag.fr/membres/Pierre.Saramito/rheolef Distibution: sources and binaries as debian packages Keywords: finite elements, numerical simulation, partial derivative equations, C++, meshes, graphics News: * equations on a surface: implements three diferent FEM methods * improves the high order Pk Lagrange interpolation implementation * ports on intel c++ 12.0 and gnu c++ 4.7 new compiler versions Previous features: Rheolef is a programming environment that serves as a convenient laboratory for computations involving finite element methods (FEM) for solving partial differential equations (PDE). Rheolef is both a C++ library and a set of commands for unix shell programming, providing algorithms and data structures. * Algorithms refer to the most uptodate ones: preconditioned sparse solvers for linear systems, incompressible elasticity, Stokes and NavierStokes flows, characteristic method for convection dominated heat problems, etc. Also nonlinear generic algorithms such as fixed point and damped Newton methods. * Data structures fit the standard variational formulation concept: spaces, discrete fields, bilinear forms are C++ types for variables, that can be combined in any expressions, as you write it on the paper. Combined together, as a Lego game, these bricks allows the user to solve most complex nonlinear problems. The concision and readability of codes written with Rheolef is certainly a major keypoint of this environment. Main features * [NEW] Massively distributed memory finite element environment, based on MPI. * [NEW] Highorder polynomial approximation. * Poisson problems in dimension d=1,2,3. * Stokes problems (d=2,3), with TaylorHood or stabilized P1 bubbleP1 elements. * linear elasticity (d=1,2,3), including the incompressible case. * characteristic method for timedependent problems: transport, convectiondifusion, and NavierStokes equations. * input and output in various file format for meshes generators and numerical data visualization systems. Advanced features * autoadaptive mesh algorithms. * axisymetric problems. * multiregions and nonconstant coefficients. * nonlinear problems with either fixedpoint algorithms or a provided generic damped Newton solver. * 3d stereo visualization Both reference manual and users guide are available. The license is GPL. Pierre Saramito  Pierre.Saramito@imag.fr Directeur de Recherche CNRS Laboratoire Jean Kuntzmann, Grenoble, France http://wwwljk.imag.fr/membres/Pierre.Saramito 
thnx mettler it was helpfull :)

Sara:
I'm aware of 2 similar things that we've done. The first is written about on our website: http://www.pointwise.com/theconnecto...ediaCFD.shtml In another case we scripted the meshing of the flow around an automatically generated array of a variety of sphere sizes. We can't show pictures of that because it was proprietary to the client. But we could describe how it was meshed. If interested in the details, send me an email at jrc@pointwise.com and I'll forward it to the right people who can give you the answer. Best Regards 
Quote:
Could you please share the info on how you created the distribution and did the meshing for your case of randomly distributed spheres. 
All times are GMT 4. The time now is 03:51. 