# Implementation of CEL Expressions

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

 July 2, 2014, 08:43 Implementation of CEL Expressions #1 New Member   rainer arenskrieger Join Date: Jun 2014 Posts: 12 Rep Power: 4 Hi, I am currently running a simulation on the fluid flow in porous media. The main purpose of the simulation is the FSI of the fluid with a deforming wall and vice versa. To make a porous media deformable I have to use the Fluid Domain and use the Source Terms. I am doing a multiphase simulation. So I can only add Bulk Sources to a subdomain. In the first step: a assign a Permability under the momentum source/loss model/permeability tab this gives me a solution. In a second step I want to implement anisotropic permeability. This should be possible by assigning the Momentum Souces to Darcys law. The permeability loss model equals to S = - viscosity*Ux/Permeability 2. try: In CEL syntax this is: Sx = - viscosity * u /Permeability Sy = - viscosity * v /Permeability Sy = - viscosity * w /Permeability These Expressions are put into Bulk Sources/ Momentum Source/ General Momentum Source tab. Viscosity and u,v,w are variables and Permeability is an Constant defined in an other CEL expression. this gives error while running the simulation although it should give the exactly same results as 1. 3.try as I am not sure how to deal with the viscosity, I tried to calculate the viscosity on every element by myself. I did this by using the Volume Fraction of each element and the calculating an average viscosity: Sx = - (Fluid 1.viscosity+(Air.viscosity - Fluid 1.viscosity)*Air.vf) * u /Permeability Sy=... Sz=... I also get errors when trying this. 4.try The Momentum Sources are only relevant to the fluid, as the viscosity of Air is very small. So in the TAB Fluid Souces/Fluid 1 ... I carried out the same steps mentioned in the 1st to 3rd try. Again I get errors. This time also when assinging a constant permeability. Does someone know what my mistake is. I guess the velocites might cause a problem. If I these Problems are solved there are more of them coming up. I do not only want to simulate anisotropic permeability but also variable permeability: What I tried so far is: Instead of assigning a constant permeability to the Sources/Bulk Sources/ Momentum Source/Permeability Tab is assigned a variable permeability of the following kind: Permeability = f(Total Mesh Displacement Z) So as the mesh moves I want the permeability to change. The values of this function are reasonable but still I get errors when running the simulation. I am now wondering whether Total Mesh Displacement Z is a Field variable which has a values for each element? Way to many questions... Can someone give me a hint on how to solve some of my problems? Thank you in advance.

 July 2, 2014, 18:14 #2 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,638 Rep Power: 98 I am no expert on porous flow so will not comment on those questions. Is Total Mesh Displacement a field variable? Yes. You can see what variables are available and where they exist in the CFX reference manual.

 July 7, 2014, 05:08 #3 New Member   rainer arenskrieger Join Date: Jun 2014 Posts: 12 Rep Power: 4 When running the simulation with a variable Permeability= f(Total Mesh Displacement) I get the following Error: | ERROR #001100279 has occurred in subroutine ErrAction. | | Message: | | Failed to calculate sources in GETSOU for LOSS_F1 from directory | | "LOSSM/RESIST" on location "Porous". Error return code is: BADL According to the manual this means: Variable can not be calculated atthe current LOCALE. As my experience using Expressions is still limited I wonder whether it is possible to use an expression of the Form Permeability= a*(Total Mesh Displacement Z)^b I would like this Expression to calculate the Permeability at every element as a function of Total Mesh Displacement Z. So as a result the Permeability should be constant at every point that is not deformed during the simulation where as the the Permeability is recalculated at every timestep for deformed elements. Is this somehow possible? Or do I have to try using User Functions and Subroutines? Thank you for your advice

 July 7, 2014, 05:57 #4 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,638 Rep Power: 98 You should be able to implement the CEL expression you describe, but you will need to properly account for units. Any value raised to a real power needs to be unitless, and a will need to have appropriate units to make it permeability. But your expression will make the permeability zero at any location not moving. That does not seem to agree with your description of what you are trying to do.

 July 7, 2014, 06:58 #5 New Member   rainer arenskrieger Join Date: Jun 2014 Posts: 12 Rep Power: 4 I did pay attention to units. the Expression I really use is different to the one I gave in my last post. It is Permeability = A*B^(c*Fibre Volume Fraction) while Fibre Volume Fraction = D/(initial cavity height + Total Mesh Displacement Z) But I think I figured out my fault: When running the FSI simulation the variable Total Mesh Displacement is only assigned to the Interface, so a limited number of elements and not all of them. But in my expression I use Total Mesh Displacement as a Variable in the whole fluid Domain. Is there a variable in the fluid domain that gives me a ratio of the initial Mesh-Element Volume / Current-Mesh-Element Volume ? So the change in Element Volume due to the Mesh Deformation? If the is something similar to this I could a current Fibre Volume Fraction which would lead to Permeability. Thank you for still answering my questions.

 Tags cel, fsi, momentum sources, permeability, total mesh displacement

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Niru CFX 0 May 10, 2013 13:13 foo7 CFX 4 September 21, 2011 09:15 rbarrett CFX 3 July 6, 2011 10:13 RAMAN6089 CFX 13 December 1, 2009 00:17 Jervds CFX 0 March 4, 2008 11:03

All times are GMT -4. The time now is 10:07.