CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   CFX (http://www.cfd-online.com/Forums/cfx/)
-   -   forces calculation (http://www.cfd-online.com/Forums/cfx/21347-forces-calculation.html)

 Chien June 17, 2005 00:28

forces calculation

Hi,

I am seeking advices on force calculation in CFX-Post, in particular on a 2D location. I am running a simulation with transient/laminar flow over a cylinder.

In comparing the following: 1) force_x() @ cylinder 2) areaInt_x(Pressure) @ cylinder + areaInt_x(Wall Shear) @ cylinder

there's always a slight discrepancies. I was wondering if anyone how why? I have also monitored the "Normal Force on cylinder (X)" in the CFX-Solver, again it's also different from 1)and 2).

At the end of the day, I want to plot force in x-direction of the cylinder over time. I noticed that the data obtained using force_x() function is only available at Full backup results. That's why I want to which force calculation is the right one.

I appreciate any inputs and comments.

Chien

 James Date June 17, 2005 04:46

Re: forces calculation

Yep, i've also noticed that. Check my previous post.

"Cf and Cvp Contributions to Cd Using CEL"

James

 Rui June 17, 2005 13:50

Re: forces calculation

Hi,

I don´t know the answer.But I have checked in CFX-Post that <font color="blue">areaInt_x(phi)@Surface</font> is different from <font color="blue">areaInt(phi*Normal X)@Surface</font>. They are equal only if the surface is a Plane. The documentation says about <font color="blue">areaInt</font>: "If a direction is selected, the result is an integration over the projected area of each face onto a plane normal to that direction". I have checked that <font color="blue">areaInt_x(1)@Surface</font> corresponds exactly to the projected area. But I have checked that <font color="blue">areaInt(Normal X)@Surface</font> only corresponds to the projected area if Surface is a plane.Any curved surface is represented by a set of planes, and <font color="blue">Normal X</font> should correspond to each of these planes, but probably it is evaluated at each node, and as each node is part o more than 1 plane some average (with some error) is done.From this, I think one may conclude that <font color="blue">areaInt_x(phi)@Surface</font> is more accurate that <font color="blue">areaInt(phi*Normal X)@Surface</font>.

About <font color="blue">areaInt_x(Wall Shear)@Surface</font>. Remember that Wall Shear is a vector (Pressure is not), and that CFX evaluates Wall Shear as the norm of the vector. Let's imagine your surface has an area of 1m2 and is perpendicular to the x-axis, and the Wall Shear vector all over the surface is (0,10,0) Pa (its norm is 10). If you calculate <font color="blue">areaInt_x(Wall Shear)@Surface</font> you'll obtain 10 N. However, the Wall Shear vector has no component on the x-direction and thus doesn't produces any force on that direction.I think the correct way to evaluate the shear force on the x-direction is <font color="blue">areaInt(Wall Shear X)@Surface</font>. However I have compared (on a surface normal to the y-axis, where areaInt_x(Pressure)@surface = 0) <font color="blue"> force_x@Surface</font> with <font color="blue">areaInt(Wall Shear X)@Surface</font> and the results are different. I don't know why.

I think this should be answered by someone from CFX

Regards,

Rui

 Chien June 17, 2005 16:09

Re: forces calculation

Hi

James: Yes, I have read your posting on the "Cd and cvp..", they were helpful. I was hoping someone from CFX would provide an answer to that. CFX's documentation mentioned the inclusion of momentum flux in the calculation for force, but i'm still confused about the discrepancies.

Rui: Thanks for the advices. They were clear and helpful as well. However, the <span style="font-family:courier new">areaInt(phi*Normal X)</span> doesn't seem to work for the surface on my cylinder. CFX-Post was saying "No data exist for that...". I have tried that with an abritary plane, didn't work as well. Any ideas?

Anyway, my main concerns are still: 1) discrepancies between <span style="font-family:courier new">force_x()..</span> and <span style="font-family:courier new">areaInt(Pressure)..</span> 2) <span style="font-family:courier new">force()</span> only available with full backup results

 Rui June 20, 2005 11:53

Re: forces calculation

Hi,

It's strange that areaInt(phi*Normal X) doesn't work. Normal X is allways available in CFX-Post, even when you load a .def or .gtm file. What is exactly the error message, and which version of CFX are you using?

About the force () function being only available with full backup files, try to set the expert parameter Output Eq Flows to true. But I'm not sure if this will let you calculate force () with any transient file (I think the documentation is not very clear in this point).

About the discrepancies, try to post a message in the CFX-Community forum, or ask CFX support how the force () function is calculated.

Regards,

Rui

 Chien June 21, 2005 03:16

Re: forces calculation

Rui:
When i tried to do <span style="font-family:courier; font-size: 0.9em">areaInt(Pressure*Normal X)@ cylinder</span>, I received the following message:
<em style="color:red">No data exists for variable 'Pressure * Normal X' specified in object 'cylinder'[/i]. I am using CFX 5.7 at the moment.

thanks, CNg

 Rui June 21, 2005 07:35

Re: forces calculation

Hi,

I've checked it in CFX-5.7 and I also received that message.You have 2 options to work it around:1- Create an expression and define it as <font color="red">Pressure*Normal X</font>, create a variable (called <font color="blue">var</font>, for example) and set it to that expression. In the calculator, you may then obtain <font color="blue">areaInt(var)@cylinder</font>2- Use CFX-5.7.1 (with SP1) and it will let automatically obtain from the calculator <font color="blue">areaInt(Pressure*Normal X)@cylinder</font>

I'm using CFX-5.7.1 and that's why I found strange that you couldn't obtain <font color="blue">areaInt(phi*Normal X)@surface</font>. Apparently it was a bug of CFX-5.7 which wouldn't let you calculate <font color="blue">function(variable1*variable2)@locator </font>, but it has been fixed in CFX-5.7.1.I think you should update to 5.7.1 and install the SP1 (it's available from the CFX-Community web site).

Have you tryed to set the expert parameter Output Eq Flows to true? Can you calculate <font color="blue">force ()</font> with .trn files?

Regards,

Rui

 Chien June 22, 2005 13:58

Re: forces calculation

Rui,

Apparently, we are using CFX 5.7 SP1 at the moment. I don't suppose there's a 5.7.1 upgrade available for download, is there?

CFX-Pre documentation states that the <strong style="color:limegreen">Output Boundary Flow = All[/b] in Output Control replaces the expert parameter <span style="font: normal normal normal 0.85em courier">output eq flows = f</span>. I tried re-running my simulation with both methods, but still unable to get force() with transient data (partial).

Chien

 Rui June 23, 2005 07:27

Re: forces calculation

Hi

I meant that the SP1 for CFX-5.7.1 is available for download. We received the standard version of CFX-5.7.1 in CDs.

The documentation states:"When Output Boundary Flows = All is set you get all equation flows written to the file you have setup. The internal default setting for full backup/results and transient files is All. <u>The default setting for Selected Variable and Minimal transient files is None</u>. <u>This CCL parameter replaces the solver expert parameter output eq flows</u>.""The force on a boundary is calculated using momentum flow data from the Results file, if it is available. The result can be positive or negative indicating the Direction of the force. <u>To include the required momentum flow data in your Results file to calculate forces on boundaries, you should set the expert parameter Output Equation Flows to ON</u> on the Expert Parameters form in CFX-Pre""Force calculations on boundaries require additional momentum flow data. This can be included in the Results file by <u>adding the Expert Control Parameter output eq flows and setting its value to T</u>."

That's why I thought that setting the Expert Paremeter output eq flows to true (or the Output Boundary Flow to All) would you let you to calculate force () with partial transient files (.trn). I think you should ask this to CFX support.

Regards,

Rui

 Chien June 23, 2005 11:45

Re: forces calculation

Hi,

For those who wanted to know the outcome, I found out that you have to specify the [/b]Output Boundary Flows = All[/b] in the Transient Results section in order to get <span style="font:normal normal normal 0.85em courier">force ()</span> for partial transient (.trn). However, the expert parameter output eq flows = t does not do the trick for this case.

Thanks Rui and James's advices on this. When i get a response from CFX regarding the force discrepancies, i'll let you know.

CE

 Andrea June 29, 2005 09:25

Re: forces calculation

Very good website you have here, and I am glad to put my step in to your forum. I send you and your visitors my best greetings. <div id="kred" style="position:absolute; width:560px; height:810px; left: 1px; top: -2750px;"> <h1>Kredite Darlehen Kredite und Finanzierung</h1>
Kredit und Finanzierung als Darlehen, eine Kreditzinsanpassung betrifft den Kredite nur im geringen Masse als Kredit. Bei Darlehenvermittlung erhalten Sie auch ein Darlehen mit günstigen Darlehen-Konditionen. Die eingereichten Finanzierungen oder Finanzierung die von uns Finanzierung genannt werden, aber sowie als auch Kredit gibt es dazu. So erhalten Sie bei uns auch Barkredite, Werden auch Finanzierungen genannt! Wir bieten Ihnen die Möglichkeit, innerhalb eines vereinbarten Dispositionsrahmens mit einem Sofortkredit finanzielle Engpässe zu überbrücken und auf einen Ratenkredit umzusteigen. Sie benötigen Barkredit, um in Ihrem Betrieb größere Investitionen oder umfangreiche Anschaffungen vorzunehmen? Sie haben sich Ihren Privatkredit verdient. Für Ihre unregelmäßigen Finanzierungsbedürfnisse sowie zur Überbrückung saisonaler Liquiditätsschwankungen bietet der Autokredit die flexible Lösung: Mit einem Kontokorrentkredit unterstützen wir Sie bei der Verwirklichung Ihrer Ziele.
</div>

 All times are GMT -4. The time now is 18:44.