Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

BR-Wamit #26

Open
treppe8k opened this issue May 29, 2023 · 20 comments
Open

BR-Wamit #26

treppe8k opened this issue May 29, 2023 · 20 comments

Comments

@treppe8k
Copy link

treppe8k commented May 29, 2023 •
edited

Dear Iñaki,

  1. Recently, i was testing BR with wamit outputs .1 and .2 files and exporting them to csv files
    My geometry is a hemisphere R=2.5 with gravity center =1.25 and ULEN=2.5.
    I was testing with ULEN = {1 , 2.5} and something is wrong with the scale

To fix it i had that ...
A= A * (rho/1000) * ulen^PowRad;
B= B * 100
Fe= Fe * (rho/10000) * grav * Aw * ulen^PowDiff;

As usual in the pitch case, PowDiff=3, PowRad=5, Aw=1 , grav=9.81 , rho=1025

Can you review these options in BEMRosetta? The idea is that at least it shows the exported csv as the original data output of wamit
I can transform them using code.

thanks
@izabala123
Copy link
Collaborator

Dear Treppe

When saving in .1, .2, ... files, the LEN of the calculation is not saved, so BR imagines itself that it is 1. However, .out files save the LEN dimension.Internally, BR follows the Wamit criteria in chapter 3. For example, for A and B.

image

Are you asking to change the LEN in BR?

Best regards
Iñaki

@treppe8k
Copy link
Author

treppe8k commented May 29, 2023 •
edited

Dear Iñaki,
"My mistake" was to use BR with outputs wamit runned with ULEN = 2.5. (although should to work !! using the formulas)
With ULEN = 1 everything is ok!

It is possible export wamit outputs (.1,..2) to the CSV files, for each DOF, using command line?
It is possible export nemoh TEC outputs to the CSV files, for each DOF, using command line?
(the nemoh dat outputs are not read.)
In the picture it is done manually.

csv

thanks

@izabala123
Copy link
Collaborator

izabala123 commented May 30, 2023 via email

@treppe8k
Copy link
Author

Dear Iñaki,
It is enough if each file exported has 2 cols , i.e {w,A} or {w,B} or {w,Fe}

bemrosetta_cl -bem "inputfilepath" -exportToCSV -"Fe-heave" "outputpath"
bemrosetta_cl -bem "inputfilepath" -exportToCSV -"A-pitch" "outputpath"

The utility is that csv is universal format and it can be used by any programming language for complex postprocessing

@izabala123
Copy link
Collaborator

Hi Treppe

I have added two new formats to save as .csv ('csv mat' and 'csv table').
Please test them to see if that is what you need.
I have added the second to the command line version, if you set .csv as the extension of the output file.

image

Best regards
Iñaki

@treppe8k
Copy link
Author

treppe8k commented Jun 2, 2023

Dear Inaki
yes, the CSV implementation is perfect... thanks a lot
i have a last point to solve...

in wamit, to run the solver, the gdf used has ulen=2.5. To display the outputs wamit, BR is configured with ULEN=2.5 but the plots of A, B, Fe are not displayed correctly. Could you verify it?
with ulen=1 everything is correct.
ulen

@izabala123
Copy link
Collaborator

izabala123 commented Jun 2, 2023 via email

@treppe8k
Copy link
Author

treppe8k commented Jun 2, 2023

i used bem data handling and this happened with .1 .2 files

ulen1
hemisphere R=2.5 , ULEN =2.5 CG 1.25
.1 generated with HO GDF1 (with LID mesh)
buoy01.txt

thanks

@izabala123
Copy link
Collaborator

izabala123 commented Jun 3, 2023 via email

@treppe8k
Copy link
Author

treppe8k commented Jun 5, 2023

dear Iñaki
When i use the *.out file everything right
thanks so much by your feedback

@izabala123
Copy link
Collaborator

izabala123 commented Jun 5, 2023 via email

@treppe8k
Copy link
Author

Dear Inaki

I detected the same issue but this time with pitch RAO results (using .4 and out files)
Case: truncated cylinder R=10 T=10 with CG=-7.5, GDF with ULEN =10

RAO Pitch
wec.4.txt
The issue is in BR because when i use Excel and i apply the formula everything is ok
form

thanks

@izabala123
Copy link
Collaborator

Hi Treppe

You are right. This is not implemented in BEMRosetta as I supposed that it was not necessary to scale nondimensional parameters. However, Wamit manual, just before that equation, says 'The non-dimensional definitions of the body motions are'.
So I will include it.

Thank you!
Iñaki

@izabala123
Copy link
Collaborator

It can be found here; https://www.wamit.com/manualv7.4/wamit_v74manualch3.html#x5-210003.3
A screenshot here:
image

@izabala123
Copy link
Collaborator

Hi Treppe

I have changed how RAO is handled. I hope now it works well with scale. Please check it for me.
Binaries are here.

Iñaki

@treppe8k
Copy link
Author

treppe8k commented Jun 19, 2023 •
edited

Hi Iñaki
the issue persist
RAO .4
raopich
RAO .out
rao pitch out

.4 files
wec.ulen1.4.txt
wec.ulen10.4.txt

ps
it would be excellent whether add the option for to see the rotation motions in deg/m and the phase in deg
in future versions ...

thanks

@izabala123
Copy link
Collaborator

Hi Treppe

.4 files does not have the length scale L so BR assigns an 1. Without the .gdf, BR does not know it.

Iñaki

@treppe8k
Copy link
Author

HI Inaki

The results shown up above were tested with .out files also
ulen10
ulen1

[treppe]

@izabala123
Copy link
Collaborator

Sorry Treppe, I am a bit lost.

If it suits you, we could have a 30 min chat to solve this. Teams or other system is good for me.
If you agree, you can send me your meeting invite through my personal email.

Best regards
Iñaki

@treppe8k
Copy link
Author

Dear Iñaki

My apologies for long time offline, i am back.

I did testings in Excel with WAMIT, HAMS and NEMOH (to compare pitch cylinder truncated coeff. A, B, Fe
and after compare them with the reading of BemRosetta )
Wamit with ULEN !=1 is perfect
Nemoh is perfect
Hams , there is a fail in the B Damping results (the conversion formula must be wrong) . A, Fe are OK

A meeting would be good, I will send e-mail

thanks

treppe

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@izabala123 @treppe8k