New tilting-pad hidrodynamic bearing models

[rodrigomoliveira1]

PR #690 is about implementing new types of bearing models to ROSS.
Certainly it'll bring a lot of discussion along, so I'm opening this issue for this purpose. Any suggestions, enhancements, new bearing models, references, examples (how to use), etc, please, report here to keep track of what's happening on ROSS.

Moreover, always add a brief summary in each commit, it helps a lot to understand the changes in the code.

(@verg1lio, @hiagopinacio, @bela-bs, @marcusfilipesr, @JuliaMota, @davidjgm)

[verg1lio]

Great to know @rodrigomoliveira1, and it's a pleasure to meet you @JuliaMota and @davidjgm. Our team will be working on it, mostly @elyqg, @marcusfilipesr and myself (@verg1lio). @bela-bs should be busy finishing the defect models for a while, and @hiagopinacio is our Python/GIT guy when we have critical decisions to make.

Other pull requests will be made for the other bearing types, in total those are: 1) Cylindrical (still to be opened); 2) Tilting-pad (open!) and 3) Thrust (being validated in Matlab), all being thermo-hydro-dynamic bearing models. We will be going for a object-oriented and modular code, using the same iterative procedures to converge the temperature, pressure and velocity fields whenever possible, so it should prove to be a very interesting code to work on. Also, an abstract bearing class with shared properties for the others to inherit is the current go-to idea, but we are open to your inputs as always. We will have a more clear view of the problem as things progress further.

Anyways, we are currently doing the hard translation of the Matlab code to Python, and will follow that with its validation and conference against the working Matlab version. After this procedure is done, we will modularize and make it object-oriented instead of linear.

It is a pleasure to be working on a challenging problem such as this, and with a capable team to top it off! Let's get to it guys.

[JuliaMota]

Good Morning!

The initiative to include a hydrodynamic tilting bearing model is very good. Currently, we have a functional code validated with the literature that calculates the pressure field, forces, equilibrium position and dynamic coefficients for a cylindrical hydrodynamic bearing. Recently, the possibility of two new geometries has been added: the elliptical bearing (or lemon bearing) and the worn bearing. Both validated and with the same returns as the cylindrical bearing.

Including the hydrodynamic tilting bearing is one of the objectives that we have not yet achieved. I believe that the code already implemented in FluidFlow can be used in some way also for this case. Anyway, I think your model should be in the FluidFlow directory, not in a separate location.

I am available for any clarification.

(@verg1lio, @hiagopinacio, @bela-bs, @marcusfilipesr )

[rodrigomoliveira1]

I agree with @JuliaMota. These new files should be located in fluidflow folder, because:

1. it belongs to the same group of classes and functions. Separating them could lead to repeated code.
2. you could take advantage of several existing methods in FluidFlow package.
3. adding new folders require adding __init__.py files to install and start the package.

[davidjgm]

Hi @verg1lio, @hiagopinacio, @bela-bs and @marcusfilipesr,
I would be pleased to help.
As your code is being committed, we can see how the currently code can be integrated with yours, since we have already done some work related to cylindrical bearings. Indeed, a model of tilting-pad bearings is a must-have in the repository; existing methods in the FluidFlow package should help.

[verg1lio]

Oh that would be all the better @JuliaMota. Our code should be compatible, or at least adaptable, to a standardized formulation without too much trouble. It is our focus at the moment, and the translation should progress quickly too. @davidjgm assistance would also be greatly appreciated, since we do our formulations ourselves, and despite knowing the physics behind it we do not have much, if any, experience with the FluidFlow package.

I believe that we should dedicate some time to familiarize ourselves with what's available at FluidFlow, before commiting to a full conversion of the more extensive convergence methods. That could save us some time and redundant work at the very least.

What is you take on the matter @marcusfilipesr, @elyqg and @hiagopinacio?

[marcusfilipesr]

I am delighted to know that you guys are going to be around us through this process, I'm sure that your help will be essential.

And I agree with @verg1lio when saying that we should dedicate some time to understand your FluidFlow package.

[JuliaMota]

We published an article at an event of the National Laboratory for Scientific Computing. It contains all the initial formulation for the pressure field and checks with the models used for short and infinitely long bearings. It is still an initial study, but it can help them better understand the formulation.

You can access the article in the proceedings of the XIII EAMC - 2020 event, available at: [http://www.eamc.lncc.br/PastEditions.php]

The article is called "Simplifications of Navier-Stokes in Journal Bearing Simulation", it was written by me, by professor @juvianna and by @flaviorangel , who contributed a lot with all the initial part of the code.

[verg1lio]

We published an article at an event of the National Laboratory for Scientific Computing. It contains all the initial formulation for the pressure field and checks with the models used for short and infinitely long bearings. It is still an initial study, but it can help them better understand the formulation.

You can access the article in the proceedings of the XIII EAMC - 2020 event, available at: [http://www.eamc.lncc.br/PastEditions.php]

The article is called "Simplifications of Navier-Stokes in Journal Bearing Simulation", it was written by me, by professor @juvianna and by @flaviorangel , who contributed a lot with all the initial part of the code.

Alright, and thanks for the reference @JuliaMota! I will check it out now and see if it is compatible to what we have implemented here. Out material is based mostly on a Thesis and a Dissertation, and on some smaller works, which we can provide here should you wish it.

[verg1lio]

Ok, so apparently your model solves the Navier-Stokes equation for the pressure and fluid velocity fields correct? Starting from the geometry, then making the dimensionless analysis with pressure field and finally moving to the Poisson formulation to get the fluid velocity field.

Our model have some differences, but we will evaluate them better tomorrow morning in an internal meeting (@marcusfilipesr @hiagopinacio @verg1lio @elyqg), and get back to you guys after we have a consensus on how to proceed. Again, thanks a lot @JuliaMota @rodrigomoliveira1 and @davidjgm for the assistance, and we will keep updating our course of action and progress here.

Have a good evening!

[JuliaMota]

We published an article at an event of the National Laboratory for Scientific Computing. It contains all the initial formulation for the pressure field and checks with the models used for short and infinitely long bearings. It is still an initial study, but it can help them better understand the formulation.
You can access the article in the proceedings of the XIII EAMC - 2020 event, available at: [http://www.eamc.lncc.br/PastEditions.php]
The article is called "Simplifications of Navier-Stokes in Journal Bearing Simulation", it was written by me, by professor @juvianna and by @flaviorangel , who contributed a lot with all the initial part of the code.

Alright, and thanks for the reference @JuliaMota! I will check it out now and see if it is compatible to what we have implemented here. Out material is based mostly on a Thesis and a Dissertation, and on some smaller works, which we can provide here should you wish it.

Of course! You can give me the references. I see it as soon as possible. As I said, this is an initial study. We've done a lot of things after that. But I believe that by synthesizing FluidFlow modeling well, it will be of great help in this beginning. Whenever I want to pass other references, including some that justify the use of this model. The parameterization of the FluidFlow is a little different from the others in order to facilitate the implementation of other types of bearings and, in the future, to enable an adaptation also for seals. However, it is different only in the way it is described, as it follows the same premises as the classic lubrication theory.

[JuliaMota]

Ok, so apparently your model solves the Navier-Stokes equation for the pressure and fluid velocity fields correct? Starting from the geometry, then making the dimensionless analysis with pressure field and finally moving to the Poisson formulation to get the fluid velocity field.

Our model have some differences, but we will evaluate them better tomorrow morning in an internal meeting (@marcusfilipesr @hiagopinacio @verg1lio @elyqg), and get back to you guys after we have a consensus on how to proceed. Again, thanks a lot @JuliaMota @rodrigomoliveira1 and @davidjgm for the assistance, and we will keep updating our course of action and progress here.

Have a good evening!

That's right! This is basically the way. Any questions, I'm here.

[verg1lio]

Hello there, and good afternoon!

I will keep working on the hard translation of the matlab code into python for the moment, since there's still a lot to go. But we will soon be having the meeting to define our path from this point on. We hope to clarify some doubts on what and how should we implement the bearing models, and their integration and functioning together with the main ROSS code.

Should you (@rodrigomoliveira1 , @JuliaMota , @davidjgm or @raphaeltimbo) need any info, just ask away anytime. We also had a general updating of our LMEst fork to the latest ROSS version recently, and i just pushed it here.

So cheers, and have a great week!
Vergilio

[verg1lio]

@JuliaMota , good afternoon!

You marked me in the Petrobrás MS Teams, but unfortunately I do not have access to that. I saw the mention, but cannot see the message. If possible, could you either give me access or resend the message here?

Thanks, and have a good week!