CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > SU2

Center of Pressure Output

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 2 Post By pcg

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 13, 2021, 00:50
Default Center of Pressure Output
  #1
New Member
 
Juan Salazar
Join Date: Jun 2019
Posts: 19
Rep Power: 7
saladbowl is on a distinguished road
Hello,

I am a recent user of SU2 and I noticed that the center of pressure result in the forces_breakdown.dat file is rather awkward. I find the following line in the file for an airfoil simulation (2D):

Code:
Center of Pressure: X=0.300838Y+0.239842 m.
The relevant source code I found is in CFlowOutput.cpp (version 7.1.1 "Blackbird"). In earlier versions it is in output_structure.cpp

Code:
    /*--- Begin forces breakdown info. ---*/

    Breakdown_file << fixed;
    Breakdown_file << "\n" << "\n" <<"Forces breakdown:" << "\n" << "\n";

    if (nDim == 3) {
      su2double m = solver_container[FLOW_SOL]->GetTotal_CFz()/solver_container[FLOW_SOL]->GetTotal_CFx();
      su2double term = (Total_CoPz/m)-Total_CoPx;

      if (term > 0) Breakdown_file << "Center of Pressure: X="  << 1/m <<"Z-"<< term << "." << "\n\n";
      else Breakdown_file << "Center of Pressure: X="  << 1/m <<"Z+"<< fabs(term);
      if (config->GetSystemMeasurements() == SI) Breakdown_file << " m." << "\n\n";
      else Breakdown_file << " in." << "\n\n";
    }
    else {
      su2double m = solver_container[FLOW_SOL]->GetTotal_CFy()/solver_container[FLOW_SOL]->GetTotal_CFx();
      su2double term = (Total_CoPy/m)-Total_CoPx;
      if (term > 0) Breakdown_file << "Center of Pressure: X="  << 1/m <<"Y-"<< term << "." << "\n\n";
      else Breakdown_file << "Center of Pressure: X="  << 1/m <<"Y+"<< fabs(term);
      if (config->GetSystemMeasurements() == SI) Breakdown_file << " m." << "\n\n";
      else Breakdown_file << " in." << "\n\n";
    }

The way it is being written, it looks like an equation. What am I missing here?
saladbowl is offline   Reply With Quote

Old   April 13, 2021, 05:41
Default
  #2
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 466
Rep Power: 13
pcg is on a distinguished road
You are not missing anything, it is the equation for a line.
It is discussed here https://github.com/su2code/SU2/discussions/1227
(you'll need to scroll down a bit)
saladbowl and ugurtan666 like this.
pcg is offline   Reply With Quote

Reply

Tags
center of pressure, forces_breakdown, output, su2


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Simple Box - Gravity with Pressure Outlet - Unrealistic Reverse Flow pyccknn FLUENT 2 December 1, 2021 18:31
How is total pressure calculated? looee CFX 15 October 30, 2019 09:14
Reference point for calculating center of pressure krihamm CFX 1 January 31, 2018 05:36
Pressure Outlet Guage pressure Mohsin FLUENT 36 April 29, 2016 18:16
static vs. total pressure auf dem feld FLUENT 17 February 26, 2016 14:04


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