PTSOLVER: P-Theta solver for dynamic lubrication
C program, uses conservative finite difference method based
on Elrod-Adams p-theta model. Reference: "A mass-conserving algorithm for dynamical lubrication problems with cavitation", R. Ausas, M. Jai and G. Buscaglia,
to appear in ASME Journal of Tribology.
Here you can download
Both programs should compile and run without trouble. Their
coding is very simple and the name of the variables has been kept as intuitive
as possible. The squeeze program will produce a file sigma.dat containing the
location of the interface as a function of time. It will also produce files
containing the profiles of pressure and theta (fluid fraction) at different times.
The journal bearing program, on the other hand, will produce a file XY.dat
containing the X and Y coordinates of the journal center as functions of time,
together with the maximum pressure, the X and Y components of the load, and
similarly for the applied load. It will also produce the files pres.gnu and
theta.gnu, containing two-dimensional data of pressure and theta, from which
2D plots can be obtained with Gnuplot (the last 5 lines of the source code
explain how to do it). Do not hesitate in sending an e-mail in case of difficulties.
Best wishes,
RFA, MJ and GCB, 12 December 2008