Editors in Chief :
Adelina Georgescu, Professor, Ilie Burdujan, Professor, Member of Academy of Romanian University of Agricultural Sciences and Veterinary Scientists, 54, Splaiul Medicine Independentei, 050094, Bucharest, “Ion Ionescu de la Brad” Iaşi ROMANIA Mihail Sadoveanu, 3, 700490 Iaşi, ROMANIA
Editorial Board:
Nuri Aksel, Professor, Dumitru Botnaru, Professor, Faculty of Applied Sciences, Chair of Superior Mathematics, Department of Applied Mechanics Technical University of Moldova, and Fluid Dynamics, Bv. Stefan cel Mare, 168, MD 2004, University of Bayreuth, D‐95440 Chişinău, REPUBLIC OF MOLDOVA Bayreuth, GERMANY
Tassos Bountis, Professor , Mitrofan Choban, Professor, Centre for Research and Applications of Nonlinear Member of Academy of Sciences of Moldova, Systems and Department of Mathematics, University of State University of Tiraspol, Iablocikin, 5, MD 2069, Patras, 26500, Patras, GREECE Chişinău, REPUBLIC OF MOLDOVA
Sanda Cleja‐Țigoiu, Professor, Constantin Fetecău, Professor, Faculty of Mathematics and Computer Science, Faculty of Mechanical Engineering, University of Bucharest, Academiei, 14, 010014, The Gh. Asachi Technical University of Iaşi, Bucharest, ROMANIA Bv. Dimitrie Mangeron, 61‐63, 700050, Iaşi, ROMANIA
Anca Veronica Ion, Senior Researcher Peter Knabner, Professor, Institute of Mathematical Statistics and Applied Chair for Applied Mathematics, Faculty of Science, Mathematics of The Romanian Academy, Friedrich Alexander University Calea 13 Septembrie, 13, 050711, Erlangen‐Nuremberg, Martensstr. 3, 91058, Erlangen, Bucharest, ROMANIA GERMANY
Boris V. Loginov, Professor, Nenad Mladenovici, Professor, Faculty of Natural Sciences, Ulyanovsk State Technical Mathematical Sciences, John Crank 210, University, Severny Venetz, 32, 432027, Brunel University, Uxbridge, UB8 3PH, Ulyanovsk, RUSSIA UNITED KINGDOM
Toader Morozan, Senior Researcher, Emilia Petrişor, Professor, Institute of Mathematics Simion Stoilow of The Romanian Department of Mathematics, Academy, Calea Griviței, 21, 010702, University Politehnica of Timişoara, Bucharest, ROMANIA Victoriei Square, 2, 300006, Timişoara, ROMANIA
Mihail Popa, Professor, Senior Researcher, Kumbakonam Rajagopal, Professor, Institute of Mathematics and Computer Science of The Department of Mathematics, Academy of Sciences of Moldova, Academiei, 5, MD 2028, Texas A&M University, Chişinău, REPUBLIC OF MOLDOVA Mailstop 3368, College Station, TX 77843‐3368, UNITED STATES OF AMERICA
Mefodie Rațiu, Professor, Senior Researcher, Mirela Ştefănescu, Professor, Institute of Mathematics and Computer Science of The Faculty of Mathematics and Computer Science, Academy of Sciences of Moldova, Academiei, 5, MD 2028, Ovidius University, Bv. Mamaia, 124, 900527, Chişinău, REPUBLIC OF MOLDOVA Constanța, ROMANIA
Nicolae Suciu, Senior Researcher, Kiyoyuki Tchizawa, Professor, Tiberiu Popoviciu Institute of Numerical Analysis, Kanrikogaku Kenkyusho, Ltd., P.O. Box 68‐1, 400110, 2‐2‐2 Sotokanda, Chiyoda‐ku, 101‐0021, Cluj‐Napoca, ROMANIA Tokyo, JAPAN
Vladilen A. Trenogin, Professor, Constantin Vârsan, Senior Researcher, Moscow Institute of Steel and Alloys, Institute of Mathematics Simion Stoilow of The Romanian B‐49, Leninsky Prospect, 4, 119049, Academy, Calea Griviței, 21, 010702, Moscow, RUSSIA Bucharest, ROMANIA
In memoriam Adelina Georgescu
April 25, 1942 – May 1, 2010
i
ii
Am cunoscut pe marea noastră matematiciană Profesor Dr. Adelina Georgescu şi am admirat-o din inimă atât pentru importanta creaţie ştiinţifică şi frumoasa activitate didactică, cât şi pentru energia inepuizabilă cu care a luptat pentru neamul românesc. Cu negrăită emoţie îmi amintesc tot ce a reuşit pentru unitatea noastră, de la primele manifestări de după 1989, Sesiuni ştiinţifice organizate în 1991-1992 şi în continuare, pentru matematicienii români de la răsărit şi apus de Prut şi lupta pentru o cât mai strânsă colaborare între ei. Spiritul de iniţiativă şi talentul deosebit de organizare stau la temelia importantei Societăţi ROMAI, pe care a creat-o în 1992 şi a seriei de Conferinţe Internaţionale CAIM iniţiată în 1993, care a ajuns în 2010 la cea de a 18-a ediţie. Doresc ca această Societate să rodească mai departe pentru înflorirea matematicii în România şi peste hotare, mai ales în Republica Moldova. La aceste minunate realizări se adaugă cei doi fii Andrei şi Sergiu Moroianu, pe care i-a crescut în dragoste pentru matematică şi care sunt astăzi la rândul lor creatori de renume internaţional în această ştiinţă şi totodată întemeietori de rodnice familii. Li se adaugă doctorii în ştiinţe matematice precum şi alţi tineri călăuziţi de Doamna Profesor pe căile acestei ştiinţe. Retrăind amintiri despre Adelina, gândul mă duce la Doamna Oltea, mama lui Ştefan cel Mare, simbol al femeii române luptând pentru dreptatea şi înflorirea României. Fie ca tradiţia făurită de Doamna Profesor Adelina Georgescu să fie continuată şi dezvoltată strălucit şi mai departe! Prof. univ. Cabiria Andreian Cazacu
Cuvant tinut in cadrul CAIM 2010 la Iasi, la lansarea cartii Matematica si viata, scrisa de Adelina Georgescu si Lucia Dragotescu.
iii
I have known our great mathematician Professor Adelina Georgescu, and I admired her from all my heart both for her important scientific creation and for the beautiful didactical activity as well as for the inexhaustible energy with which she fought for the Romanian people. With unspeakable emotion I remember all she managed to do for our unity, from the first manifestations after 1989, the sessions organized in 1991-1992, and after, mainly for the Romanian mathematicians from the East and West of Prut and the fight for their close collaboration. Her vivid spirit and organizing talent lie at the base of the important Mathematical Society ROMAI, which she founded in 1992, and of the series of Conferences on Applied and Industrial Mathematics (CAIM) that she initiated in 1993, and that reached in 2010 the 18th edition. I wish that this Society to be also fruitful from now on, for the flourishing of the mathematics from Romania and abroad, especially from Republic of Moldova. To this wonderful achievements we must add her two sons, Andrei and Sergiu Moroianu, that she grew up in the spirit of love for mathematics and that are now, at their turn, well-known mathematicians and also heads of beautiful families. To them we must add the doctors in mathematics and many other young persons that were lead by Professor Adelina Georgescu on the ways of this science. While reminding Adelina, the thought takes me to Doamna Oltea, mother of Stephen the Great, the symbol of Romanian women fighting for justice and for the flourishing of our country, Romania. Let the tradition developed by Professor Adelina Georgescu be brightly continued and developed further!
Professor Cabiria Andreian Cazacu.
Speech held by Professor Cabiria Andreian Cazacu, at CAIM 2010, in Iasi, when the book Mathematics and Life, written by Adelina Georgescu and Lucia Dragotescu was launched.
iv
A life dedicated to the whole "Knowledge" Liliana RESTUCCIA Department of Mathematics, University of Messina, Italy Considerate la vostra semenza: fatti non foste a viver come bruti, ma per seguir virtute e conoscenza. (Dante Alighieri, Divina Commedia, Inferno, Canto XXVI, versi 118-120)
This contribution contains the speech delivered by the author with the occasion of the Mini‐ symposium, organized in collaboration with D. Jou, M. S. Mongiovì and W. Muschik, as a part of the 2010 Conference of SIMAI (Italian Society of Applied and Industrial Mathematics), which was held between 21‐25 June in Cagliari, Sardegna, Italia. This Mini‐symposium, entitled "Recent Ideas in Non‐ Equilibrium Thermodynamics and Applications" was dedicated to the memory of Prof. Dr. Adelina Georgescu, co‐researcher and a great author’s friend.
We dedicate this Mini-symposium, regarding the latest ideas in Non- Equilibrium Thermodynamics and Applications to physical-mathematical models of complex materials, to Prof. Adelina Georgescu, a member of the Scientists’ Academy of Romania, a great master and an eminent scientist. Her death, on May 1st 2010, leaves behind a very great emptiness in the field of Applied Mathematics and not less in the lives of those who knew and worked with her. I am deeply grateful for her very precious advices in the organization of my research and the very enlightening discussions and remarks regarding a lot of mathematical concepts and tools, with their practical applications in different fields of science. She was my great and unforgettable friend, invaluable guide and scientific collaborator. Since 2002 I invited her as visiting professor at the University of Messina for giving lectures and seminars on: Non-linear Dynamics, Bifurcation theory, Variational problems in Mathematical-Physics, Thermodynamics of fluids, Asymptotic methods with applications to waves and shocks, Complex fluid flows, Modeling in Meteorology, Turbulence, Modeling in cancer dynamics. I always thanked her for her profound participation and for her very precious contribution in spreading her research methodologies (a scientific literature on foundations and recent progress in Mathematical-Physics due to the Romanian and Russian schools) and, especially for her publications (very useful due to their synthetic style and clarity in explaining all concepts) at the University of Messina. Some of her books and publications have been already translated from Romanian to Italian, in order to be known in Italy for their invaluable qualities. v
She also was a corresponding member of "Accademia Peloritana dei Pericolanti" of Messina; a special issue of the Academic Acts (on-line) of this institution, will be published, containing a collection of the lectures held by Prof. Adelina Georgescu during her visits at Messina. She was the founder and president of ROMAI (Romanian Society of Applied and Industrial Mathematics) since 1992. Furthermore, she founded the Institute of Applied Mathematics of Romanian Academy (IAM), which she led as director until 1995. In 1997 she became full professor at the University of Piteşti (Romania), where she founded the Chair of Applied Mathematics (the head of which she was for several years), where she created the group of research on Dynamical systems and Bifurcation. She was the editor in chief of ROMAI Journal since 2005. Since 1993, year by year she has organized the CAIM International Conferences, devoted to contributions in Applied and Industrial Mathematics, the Romanian analogous to SIMAI International Conferences in Italy. She delivered about 200 conferences at several universities and research units mostly from abroad. Moreover, she sustained short or plenary conferences at about 150 scientific meetings. She published (some are still in press) about 20 books (research monographies, three of them at Kluwer and Chapman & Hall and one at World Scientific, dictionaries, universitary texts), alone or in collaboration, and more than 200 papers in scientific journals most of them in ISI journals or refereed in MR and ZBM. She was awarded the Romanian Academy prize and many other distinctions. She was honored as a Doctor Honoris Causa of the State University of Tiraspol (Rep. of Moldova). She was a corresponding member of the Academy of Nonlinear Sciences (Moscow). The main topics of Prof. Adelina Georgescu's researches were in Hydrodynamics and applications to complex fluid flows, Turbulence, Meteorology, Perturbation theories for partial differential equations, Boundary layer theory, Hydrodynamics stability, Study of solutions of partial differential equations, Non- linear dynamics, Bifurcation theory, Variational problems in Mathematical-Physics, Modeling in cancer dynamics, Biomathematics, Synergetics. She dedicated to the progress of science her very generous and honest and highly moral and responsible life, studying very deeply mathematical and physical theories with applications and giving them new impulse. Her courageous ideas in a difficult politic framework had been well known by whole scientific world. She was a gifted and talented woman with a very great humanity and
vi
Professors Adelina Georgescu and Liliana Restuccia, at University of Messina, 2007
very broad interests. Besides her basic research items like Physics and Mathematics (Applied Mathematics, in particular), she was interested in politics, history, geography, languages, religion, literature, medicine. At the University of Messina we will be open a Study Center "Adelina Georgescu" of Applied Mathematics, dedicated to her memory, for her work be known by the present generations and those to come. This Center has as objective to spread at international level her books and publications, her research methodologies, her ideas, in general, and in particular those in Synergetics (the new synthesis of the science) and in the directions to be followed in order to find the science truths. Furthermore, this Center has also the goal of organizing international conferences, meetings, study days, lectures and seminars for discussing issues as: foundations and recent progresses in different fields of Applied Mathematics, with the participa- tion of speakers of high scientific level, researchers and students. I think that if there is a place where there is now Adelina Georgescu, well, then this is the Paradise, the Dante Alighieri’s Paradise of the whole "Knowledge", to which Adelina devoted her entire life.
vii
viii
Published scientific works of Professor ADELINA GEORGESCU
A. BOOKS
1. Adelina Georgescu, Bifurcatia, fractali si haos determinist, chapter of Enciclopedia matematica, edited by M. Iosifescu, O. Stănăşilă, D. Ştefănoiu, Editura AGIR, 2010. 2. Adelina Georgescu, L. Palese, G. Raguso, Biomatematica. Modelli, dinamica e biforcazione, Cacucci Editore, Bari, 2009. 3. Adelina Georgescu, Lidia Palese, Stability criteria of fluid flows, Series on Advances in Mathematics for applied sciences 81, World Scientific, Singapore, 2009. 4. Adelina Georgescu, George-Valentin Cârlig, Cătălin-Liviu Bichir, Ramona Radoveneanu, Matematicieni români de pretutindeni, ed. a II-a, , Seria Mat. Apl. Ind. 24, Ed. Pamantul, Pitesti, 2006. 5. Adelina Georgescu, C.-L. Bichir, G.V. Cîrlig, Matematicieni români de pretutindeni, Seria de Mat. Apl. Ind. 18, Ed. Univ. Piteşti , Piteşti, 2004. 6. C. Rocşoreanu, Adelina Georgescu, N. Giurgiţeanu, FitzHugh-Nagumo model: bifurcation and dynamics, Kluwer, Dordrecht, 2000. 7. B.-N. Nicolescu, N. Popa, Adelina Georgescu, M.Boloşteanu, Mişcări ale fluidelor cavitante.Modelare şi soluţii, Seria Mat. Apl. Ind, 3, Ed. Univ. Piteşti, Piteşti, 1999. 8. N. Popa, B.-N. Nicolescu, Adelina Georgescu, M. Boloşteanu, Modelări matematice în teoria lubrificaţiei aplicate la etanşările frontale, Seria Mat. Apl. Ind., 2, Ed. Univ. Piteşti, Piteşti, 1999. 9. Adelina Georgescu, M. Moroianu, I. Oprea, Teoria bifurcaţiei. Principii şi aplicaţii, Seria Mat. Apl. Ind. 1, Ed. Univ. Piteşti, Piteşti, 1999. 10. Adelina Georgescu, Teoria stratului limită. Turbulenţă, Ed. Gh. Asachi, Iaşi, 1997. 11. Adelina Georgescu, Asymptotic treatment of differential equations, Applied Mathematics and Mathematical Computation, 9, Chapman and Hall, London, 1995. 12. Adelina Georgescu, I. Oprea, Bifurcation theory from application viewpoint, Tipografia Univ. Timişoara, 1994, (Monografii Matematice 51). 13. Adelina Georgescu, Sinergetica. Solitoni. Fractali. Haos determinist. Turbulenţă, Tipografia Univ. Timişoara, 1992. 14. H.Dumitrescu, Adelina Georgescu, V.Ceangă, GH.Ghiţă, J.Popovici, B.Nicolescu, Al.Dumitrache, Calculul elicei, Ed.Academiei, Bucureşti, 1990. 15. Adelina Georgescu, Aproximaţii asimptotice, Ed.Tehnică, Bucureşti, 1989.
ix
16. Adelina Georgescu, Sinergetica-o nouă sinteză a ştiinţei, Ed.Tehnică, Bucureşti, 1987. 17. Adelina Georgescu, Hydrodynamic stability theory, Mechanics: Analysis, 9, Kluwer, Dordrecht, 1985. 18. St.N.Săvulescu, Adelina Georgescu, H.Dumitrescu, M.Bucur, Cercetări matematice în teoria modernă a stratului limită, Ed. Academiei, Bucureşti, 1981. 19. Adelina Georgescu, Teoria stabilităţii hidrodinamice, Ed. Ştiinţifică şi Enciclopedică, Bucureşti, 1976.
B. PAPERS
2010
1. Limit cycles by Finite Element Method for a one - parameter dynamical system associated to the Luo - Rudy I model, ROMAI J., this issue, 27-41. (with C. L. Bichir, B. Amuzescu, Gh. Nistor, et al.). 2. A linearization principle for the stability of the chemical equilibrium of a binary mixture, ROMAI J., this issue, 131-138 (with L. Palese). 3. Asymptotic waves as solutions of nonlinear PDEs deduced by double-scale method in viscoanelastic media with memory, ROMAI J., this issue, 139-154 (with L. Restuccia). 4. An application of double-scale method to the study of non-linear dissipative waves in Jeffreys media, submitted to Math. and Its Appl. (with L. Restuccia). 5. On the nonlinear stability of a binary mixture with chemical surface reactions, submitted to Math. and Its Appl. (with L. Palese). 6. Stability and self-sustained oscillations in a ventricular cardiomyocyte model, submitted to Interdisciplinary Sciences – Computational Life. (with C. L. Bichir, B. Amuzescu, M. Popescu, et al.). 7. Limit points and Hopf bifurcation points for a one - parameter dynamical system associated to the Luo - Rudy I model, to be published (with C. L. Bichir, B. Amuzescu, Gh. Nistor, et al.). 2009
8. Further results on approximate inertial manifolds for the FitzHugh-Nagumo model, Atti dell’ Accad. Peloritana dei Pericolanti, Scienze FMN, vol LXXXVII, nr. 2. (2009). (with C.-S. Nartea). 9. On the stability bounds in a problem of convection with uniform internal heat source, Atti dell’ Accad. Peloritana dei Pericolanti, Scienze FMN, acceptata. (with F.-I. Dragomirescu).
x
10. Degenerated Bogdanov-Takens bifurcations in an immuno-tumor model, Atti dell’ Accad. Peloritana dei Pericolanti, Scienze FMN, vol LXXXVII, nr. 1. (2009). (with M. Trifan). 11. Stability criteria for quasigeostrofic forced zonal flows ; I. Asymptotically vanishing linear perturbation energy, ROMAI J. 5, 1 (2009). 63-76 (with L. Palese). 12. Linear stability results in a magnetothermoconvection problem, An St. Univ. Ovidius Constanta, 17, 3 (2009), 119-129. (with F.-I. Dragomirescu). 13. Polynomial based methods for linear nonconstant coefficients eigenvalue problems, Proceedings of the Middle Volga Mathematical Society, 11, 2 (2009), 33-41. (with F.-I. Dragomirescu). 14. Thermodynamics of fluids, SIMAI e-Lecture Notes, Vol. 2 (2009), 1-26.
2008
15. A closed-form asymptotic solution of the FitzHugh-Nagumo model, Bul. Acad. St. Rep. Moldova, Seria Mat., 2 (2008), 24-34. (with Gh. Nistor, M.- N. Popescu, D. Popa). 16. A closed-form asymptotic solution of the FitzHugh-Nagumo model, Bul. Acad. St. Rep. Moldova, Seria Mat., 2 (2008), 24-34. (with Gh. Nistor, M.- N. Popescu, D. Popa). 17. Application of two spectral methods to a problem of convection with uniform internal heat source, Journal of Mathematics and Applications, 30 (2008), 43-52. (with F.-I. Dragomirescu). 18. A linear magnetic Bénard problem with Hall effect. Application of the Budiansky-DiPrima method, Trudy Srednevoljckogo Matematichescogo Obshchestva, Saransk, 10, 1 (2008), 294-306. (with L. Palese). 19. Determination of asymptotic waves in Maxwell media by double-scale method, Technische Mechanik, 28, 2 (2008), 140-151. (with L. Restuccia). 20. Approximate limit cycles for the Rayleigh model, ROMAI J. 4, 2 (2008), 73- 80. (with M. Sterpu, P. Băzăvan). 21. Lyapunov stability analysis in Taylor-Dean systems, ROMAI J. 4, 2 (2008), 81-98. (with F.-I. Dragomirescu).
2007
22. Analytical versus numerical results in a microconvection problem, Carpathian J. Math., 23, 1-2 (2007), 81-88. MR 2305839 (with F.-I. Dragomirescu) . 23. Relaxation oscillations and the “canard” phenomenon in the FitzHugh- Nagumo model, Cap. 4 in Recent Trends in Mechanics, 1, Ed. Academiei Romane, Bucuresti, 2007, 82-106 (with M.-N. Popescu, Gh. Nistor, D. Popa)
xi
24. Dynamical approach in biomathematics, ROMAI J., 2, 2 (2007), 63-76. (with L. Palese, G. Raguso).
2006
25. Some results on dynamics generated by Bazykin model, Atti Accad. Peloritana dei Pericolanti, Scienze FMN, 84, C1A0601003 (2006), 1- 10.(with R. - M. Georgescu). 26. A linear magnetic Bénard problem with tensorial electrical conductivity, Bolletttino U.M.I. (8) 9-B (2006), 197-214. (Prima data publicata ca Rapp. Int. Dipto. Mat. Univ. Bari, 17 / 03.) (with L. Palese, A. Redaelli) 27. Further study of a microconvection model using the direct method, ROMAI J., 2, 1 (2006), 77-86. (with I. Dragomirescu). 28. Bifurcations and dynamics in the lymphocytes-tumor model, The 7th Congress of SIMAI 2004, Venezia, Italy, 2004. International Conference on Mathematical Models and Methods in Biology and Medicine -MMMBM 2005, Bedlewo, Poland, May 29-June 3, 2005, trimisa la Atti Accad. Peloritana dei Pericolanti, Scienze FMN. (with M. Trifan) . 29. Asymptotic waves from the point of view of double-scale method, Atti Accad. Peloritana dei Pericolanti, Scienze FMN, 84, C1A0601005 (2006), 1- 9. (with L. Restuccia) . 30. Linear stability bounds in a convection problem for variable gravity field, Bul. Acad. St. Rep. Moldova, Mat., 3(2006), 51-56. (with I. Dragomirescu) . 31. Neutral manifolds in a penetrative convection problem. I. Expansions in Fourier series of the solutions, Sci. Annals of UASVM Iasi, 49, 2 (2006), 19-32. MR 2300509 (with A. Labianca) . 32. Mathematical models in biodynamics, Sci. Annals of UASVM Iasi, 49, 2 (2006), 361-371.
2005
33. The static bifurcation in the Gray – Scott model, Rev. Roum. Sci. Tech. – Méc. Appl., 50, 1-3 (2005), 3-13. (with R. Curtu) . 34. A nonlinear hydromagnetic stability criterion derived by a generalized energy method, Bul. Acad. St. Rep. Moldova, Seria Mat., 1(47) (2005), 85- 91. (with C.-L. Bichir, L. Palese) . 35. A direct method and its application to a linear hydromagnetic stability problem, ROMAI J., 1, 1 (2005), 67-76. (with L. Palese, A. Redaelli) . 36. Sets governing the phase portrait (approximation of the asymptotic dynamics), ROMAI J., 1, 1 (2005), 83-94. (with S.-C. Ion).
xii
37. Application of the direct method to a microconvection model, Acta Universitatis Apulensis, Alba Iulia, Mathematics - Informatics, 10 (2005), 131-142. ZB 1113.65076, MR 2240129. (with F.- I. Dragomirescu). 38. Normal forms and unfoldings: a comparative study, Sci. Annals of UASVM Iaşi, 48, 2 (2005), 15-26 (with E. Codeci) . 39. Bifurcation of planar vector fields. I. Normal forms "at the point". One zero eigenvalue, Bul. Şt. Univ. Piteşti, Ser. Mec. Apl., 1 (12) (2005), 41-47. (with D. Sârbu ) . 40. Bifurcation of planar vector fields. II. Normal forms "at the point". Hopf and cups cases, Bul. Şt. Univ. Piteşti, Ser. Mec. Apl., 1 (12) (2005), 49-55. (with D. Sârbu ). 41. Continuity of characteristics of a thin layer flow driven by a surface tension gradient, ROMAI J., 1, 2 (2005), 11-16. (with E. Borsa). 2004
42. Bifurcation in the Goodwin model I, Rev. Roum. Sci. Tech. – Méc. Appl., 49, 1-6 (2004), 13-16. (with N. Giurgiţeanu, C. Rocşoreanu). 43. Curba valorilor de bifurcaţie Hopf pentru sisteme dinamice plane, Bul. St., Seria Mec. Apl., 10 (2004), 55-62. (with E. Codeci) . 44. Dynamic bifurcation diagrams for some models in economics and biology, Acta Universitatis Apulensis, Alba Iulia, Mathematics-Informatics, 8 (2004), 156-161. 45. On instability of the magnetic Bénard problem with Hall and ion-slip effects, Intern. J. Engng. Sci., 42 (2004), 1001-1012. (with L. Palese) . 46. Liapunov method applied to the anisotropic Bénard problem, Math. Sci. Res. J., 8, 7 (2004), 196-204. (with Lidia Palese) . 2003
47. A Lie algebra of a differential generalized FitzHugh – Nagumo system, Bul. Acad. Şt., Rep. Moldova, Seria Mat. 1 (41) (2003), 18-30. (with M. Popa, C. Rocşoreanu) 48. Approximation of pressure perturbations by FEM, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 9 (2003), 31-36. (with C. - L. Bichir) . 49. Global bifurcations for FitzHugh-Nagumo model, Dynamical Systems and Applications, Trends in Mathematics: Bifurcations, Symmetry and Patterns, Birkhäuser, Basel, 2003, 197-202. (with C. Rocşoreanu, N. Giurgiţeanu) . 50. Static and dynamic bifurcation of nonlinear oscillators, Bul. Şt. Univ. Piteşti, Seria Mec. Apl., 1, 7 (2003), 133-138. 51. A Lorenz-like model for the horizontal convection flow, Int. J. Non-Linear Mech., 38 (2003), 629-644. (with E. Bucchignani, D. Mansutti) . 52. Bifurcation in biodynamics, Sci. Annals of UASVM Iasi, 46, 2 (2003), 15- 34. 53. Numerical integration of the Orr-Sommerfeld equation by wavelet methods, Bul. St. Univ. Pitesti, Seria Mat.-Inf., 9 (2003), 25-30. (with L. Bichir). xiii
2002
54. Codimension - one bifurcations for a Rayleigh model, Bul. Acad. Şt. Rep. Moldova, Seria Mat., 1 (38), (2002), 69 – 76. (with M Sterpu). 55. Static bifurcation diagram for a microeconomic model, Bul. Acad. Şt. Rep. Moldova, 3 (40) (2002), 21-26. (with L. Ungureanu, M. Popescu). 56. Improved criteria in convection problems in the presence of thermodiffusive conductivity, Analele Univ. Timişoara, 40, 2 (2002), 49-66. (with L. Palese) . 57. Existence and regularity of the solution of a problem modelling the Bénard problem, Mathematical Reports, 4 (54), 1 (2002), 87-102. (with A. – V. Ion). 58. Degenerated Bogdanov – Takens points in an advertising model, Analele Universităţii de Vest, Seria Economie; Sudiile de Economie, Timişoara, (with L. Ungureanu). 59. Global bifurcations in an advertising model, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., (with Laura Ungureanu, Liviu Ungureanu) . 60. Stability criteria for quasigeostrophic forced zonal flows. I. Asymptotically vanishing linear perturbation energy, Magnetohydrodynamics: an International J. (with L. Palese). 61. Domains of attraction for a model in enzimology, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 8 (2002), 49-58. (with R. Curtu) . 62. Heteroclinic bifurcations for the FitzHugh – Nagumo system, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 8 (2002). (with C. Rocşoreanu, N. Giurgiţeanu). 63. Topological type of some nonhyperbolic equilibria in a problem of microeconomic dynamics, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 8 (2002). (with L. Ungureanu, M. Popescu). 64. On the Misra-Progogine-Courbage theory of irreversibility, Mathematica, 44 (67), 2 (2002), 215-231. (with N. Suciu). 65. Non – Newtonian solution and viscoelastic constitutive equations, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 5 (2000). (with C. Chiujdea). 66. Normal form for the degenerated Hopf bifurcation in an economic model, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 5 (2000). (with L. Ungureanu). 67. k > 3 order degenerated Bautin bifurcation and Hopf bifurcation in a mathematical model of economical dynamics, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., (2002). (with L. Ungureanu, M. Popescu). 68. Static bifurcation diagram for a mathematical model governing the capital of a firm, Bul. Şt. Univ. Piteşti, Seria Mat.- Inf, 8 (2002), 177-181. (with L. Ungureanu). 69. Concavity of the limit cycles in the FitzHugh-Nagumo model, An. Şt. Univ. Al.I.Cuza, Iaşi, Mat. (N. S.), 47, 2 (2001), 287-298, (2002). (with C. Rocşoreanu, N. Giurgiţeanu).
xiv
2001
70. Concavity of the limit cycles in the FitzHugh-Nagumo model, Analele Univ. Iaşi, Seria I Matematica, 47, 2 (2001), 287-298. (with C. Rocşoreanu, N. Giurgiţeanu) . 71. Codimension–three bifurcations for the FitzHugh–Nagumo dynamical scheme, Mathematical Reports, 3 (53), 3 (2001), 287 – 292. (with M. Sterpu). 72. Classes of solutions for a nonlinear diffusion PDE, J. of Comput. Appl. Math., 133, 1-2 (2001), 373 - 381. (with H. Vereecken, H. Schwarze, U. Jaekel). 73. On special solutions of the Reynolds equation from lubrication, J. of Comput. Appl. Math., 133, 1-2 (2001), 367 - 372. (with B. Nicolescu, N. Popa, M. Boloşteanu). 74. Conections between saddles for the FitzHugh–Nagumo system, Int. J. Bif. Chaos, 11, 2 (2001), 533 – 540. (with C. Rocşoreanu, N. Giurgiţeanu). 75. The complete form for the Joseph extended criterion, Ann. Univ. Ferrara, Sez. VII (N. S.), Sc. Mat., 47 (2001), 9 – 22. (with L. Palese, A. Redaelli). 76. Determination of neutral stability curves for the dynamic boundary layer by splines, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 7 (2002), 15-22. (with L. Bichir). 77. Numerical integration of the Orr-Sommerfeld equation by wavelet methods, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 7 (2001), 9-14. (with L. Bichir). 78. Degenerated Bogdanov-Takens points in an advertising model, Bul. Şt. Univ. Piteşti, Seria Mat. – Inf., 7 (2001), 173-177. (with L. Ungureanu). 2000
79. Dynamics generated by the generalized Rayleigh equation. II. Periodic solutions, Mathematical Reports, 2(52), 3 (2000), 367 – 378, 2001. (with M. Sterpu, P. Băzăvan). 80. Neutral curves for the MHD Soret – Dufour driven convection, Rev. Roum. Sci. Tech. – Méc. Appl., 45, 3 (2000), 265 – 275. (with S. Mitran, L. Palese). 81. On a method in linear stability problems. Application to natural convection in a porous medium. J. of Ultrascientist of Physical Sciences, 12, 3 (2000), 324 – 336. (with L. Palese). 82. On the Misra-Prigogine-Courbage theory of irreversibility. 2. The existence of the nonunitary similarity, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 6 (2000), 213 - 222. (with N. Suciu). 83. Codimension – three bifurcation for a FitzHugh-Nagumo like system, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 6 (2000), 193 - 197. (with M. Sterpu). 84. Dynamics and bifurcation in the periodically forced FitzHugh-Nagumo system, Intern. J. of Chaos Theory and Applications, 5, 2 (2000), 63 - 79. (with M. Sterpu). xv
85. On a new method in hydrodynamic stability theory, Math. Sciences Research Hot – Line, 4, 7 (2000), 1 – 16. (with L. Palese, A. Redaelli). 86. Hopf and homoclinic bifurcations in a biodynamical system, Bul. Şt. Univ. Baia Mare, Seria Mat–Inf., 16, 1(2000), 131–142. (with C. Rocşoreanu, N. Giurgiţeanu). 87. Dynamics of the cavitation spherical bubble. II. Linear and affine approximation, Rev. Roum. Sci. Tech. – Méc. Appl., 45, 2 (2000), 163 – 175. (with B.-N. Nicolescu). 88. Degenerated Hopf bifurcation in the FitzHugh – Nagumo system. II. Bautin bifurcation, Mathematica- Anal. Numér. Theor. Approx., 29, 1 (2000), 97 – 109. (with C. Rocşoreanu, N. Giurgiţeanu). 1999
89. Set of attraction of certain initial data in a nonlinear diffusion problem, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 3 (1999), 235-261 (with H. Schwarze, H. Vereecken, U. Jaekel). 90. Coincidence of the linear and nonlinear stability bounds in a horizontal thermal convection problem, Intern. J. Nonlin. Mech., 34, 4 (1999), 603-613. (with D. Mansutti). 91. New types of codimension-one and-two bifurcations in the plane, Inst. Matem. Acad. Rom, Preprint No. 12/1999. (with C. Rocşoreanu, N. Giurgiţeanu). 92. Regimes with two or three limit cycles in the FitzHugh-Nagumo system, ZAMM 79, Supplement 2 (1999). (with C. Rocşoreanu, N. Giurgiţeanu). 93. Asymptotic analysis of solute transport with linear nonequilibrium sorption in porous media, Transp. Porous Media, 36, 2 (1999), 189-210 (with H. Vereecken, U. Jaekel). 94. Symmetry of the solution of the nonlinear Reynolds equation describing mechanical face seals, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 3 (1999), 333-343. (with B. Nicolescu, N. Popa). 95. Hopf bifurcation and canard phenomenon in the FitzHugh-Nagumo model, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 3 (1999), 217-233. (with C. Rocşoreanu, N. Giurgiţeanu). 96. Convecţia termică cu efect Marangoni. Condiţii de echilibru, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 3 (1999), 345-350. (with Gh. Nistor). 97. Applications of coarse-grained and stochastic averages to transport processes in porous media, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 3 (1999), 435-445. (with N. Suciu, C. Vamos, U. Jaekel, H. Vereeken). 98. Modelul continuu multiplicator-accelerator. II. Cazul liniar pentru anumite valori negative ale parametrilor şi cazul neliniar, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 3 (1999), 263-266. (with C. Georgescu).
xvi
99. Investigation of the normalized Gierer-Meinhardt system by center manifold method, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 3 (1999), 277-283. (with A. Ionescu). 100. Dynamics and bifurcations in a biological model, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 4 (1999), 137 – 153. (with N. Giurgiţeanu, C. Rocşoreanu). 101. On an inertial manifold in the dynamics of gas bubbles, Rev. Roum. Sci. Tech. – Méc. Appl., 44, 6 (1999), 629 – 631. (with B. Nicolescu). 1998
102. Neutral thermal hydrodynamic and hydromagnetic stability hypersurfaces for a micropolar fluid layer, Indian J. Pure and Appl. Math.,29, 6 (1998), 575-582. (with M. Gavrilescu, L. Palese). 103. On the Misra-Prigogine-Courbage theory of irreversibility, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 2 (1998), 169-188. (with N. Suciu). 104. On the mechanism of drag reduction in Maxwell fluids, Bul. Şt. Univ. Piteşti, Seria Matematică şi Informatică, 2 (1998), 107-114. (with C. Chiujdea). 105. Transport processes in porous media. 1. Continuous modelling, Romanian J. Hydrology Water Resources, 5, 1-2 (1998) 39-56. (with N. Suciu, C. Vamoş, U. Jaekel, H. Vereecken). 106. Equilibria and relaxation oscillations of the nodal system of the heart. 2. Hopf bifurcation, Rev. Roum. Sci. Tech.-Méc. Appl., 43, 3, (1998), 403 – 414 (with C. Rocşoreanu, N. Giurgiţeanu). 107. Neutral surfaces for Soret – Dufour – driven convective instability, Rev. Roum. Sci. Tech. – Méc. Appl., 43, 2 (1998), 251 – 260. (with L. Palese, L. Pascu). 1997
108. On the existence and on the fractal and Hausdorff dimensions of some global attractor, Nonlinear Anal., Theory, Methods & Applications, 30, 8, (1997), 5527-5532. (with A. V. Ion). 109. Stability spectrum estimates for confined fluids, Rev. Roum. Math. Pures Appl., 42,1-2 (1997), 37-51. (with L. Palese). 110. Thermosolutal instability of a compresible Soret-Dufour mixture with Hall and ion-slip currents through a porous medium, Rev. Roum. Sci. Tech.-Méc. Appl., 42, 3-4 (1997), 279-296 (with L. Palese, D. Paşca, D. Bonea). 111. Studiul portretului de fază. II. Punctele de inflexiune ale traiectoriilor de fază ale sistemului dinamic Van der Pol, St. Cerc. Mec. Apl., 56, 1-2 (1997), 15-31. (with N. Giurgiţeanu). 112. Studiul portretului de fază. III. Influenţa liniarizării asupra sistemului dinamic neliniar, St. Cerc. Mec. Apl., 56, 3 - 4 (1997), 141-153. (with N. Giurgiţeanu).
xvii
113. Hydrodynamical equations for one-dimensional systems of inelastic particles, Phys. Rev., E (3), 55, 5 (1997), 6277-6280. (with C. Vamoş, N. Suciu). 114. Modelul continuu multiplicator-accelerator. Cazul liniar, Bul. Şt. Seria Mat.-Inform. Univ. Piteşti, 1 (1997), 95-104. (with C. Georgescu). 115. Bifurcation in the Goodwin model from economics. II, Bul. Şt. Seria Mat.- Inform. Univ. Piteşti, 1 (1997), 105-112. (with N. Giurgiţeanu, C. Rocşoreanu). 116. Studiul portretului de fază IV. Absenţa bifurcaţiei canard, St. Cerc. Mec. Apl., 56, 5-6 (1997), 297-305. (with N. Giurgiţeanu, C. Rocşoreanu). 117. Degenerated Hopf bifurcation in the FitzHugh-Nagumo system. 1. Bogdanov-Takens bifurcation, Analele Univ. din Timişoara, 35, 2 (1997), 285-298. (with C. Rocşoreanu, N. Giurgiţeanu). 1996
118. Neutral stability hypersurfaces for an anisotropic MHD thermodiffusive mixture. III. Detection of false secular manifolds among the bifurcation characteristic manifolds, Rev. Roum. Math. Pures Appl.,41,1-2 (1996), 35- 49. (with L. Palese). 119. Balance equations for physical systems with corpuscular structure, Physica A, 227 (1996), 81-92. (with C. Vamoş, N. Suciu, I. Turcu). 120. Balance equations for a finite number of material points, Stud. Cerc. Mat., 48, 1-2 (1996), 115-127. (with C. Vamoş, N. Suciu). 121. A nonlinear stability criterion for a layer of a binary mixture, ZAMM Supplement 2, 76 (1996), 529-530. (with L. Palese). 122. Extension of the Joseph's criterion on the nonlinear stability of mechanical equilibria in the presence of thermodiffusive conductivity, Theoret. Comput. Fluid Dyn., 8, 6 (1996), 403-413.( with L. Palese). 123. Coarse grained averages in porous media, KFA / ICG - 4 Internal Report No. 501296/1996. (with N. Suciu, C. Vamoş, U. Jaekel, H. Vereecken). 124. Asymptotic analysis of nonlinear equilibrium solute transport in porous media, Water Ressources Research, 32, 10, (1996), 3093-3098. (with U. Jaekel, H. Vereecken). 1995
125. Amélioration des estimations de Prodi pour le spectre, C. R. Acad. Sci. Paris Sér. I Math., 320, 7 (1995), 891-896. (with L. Palese). 126. Sulla stabilitá globale del equilibrio meccanico per una miscela binaria in presenza di effetti Soret e Dufour, Rapp. Dipto. Mat., Univ. Bari, 8/1995. (with L. Palese, A. Redaelli).
xviii
1994
127. Nonlinear stability criteria for l MHD flows. I. IsothermaI isotropic case, Rev. Roum. Math. Pures Appl., 39, 2 (1994), 131-146, (cu M. Maiellaro, L. Palese). 1993
128. Synergetics and synergetic method to study processes in hierarchical systems, Noesis, 18 (1993), 121-127. 129. The application of the shooting method to the hydrodynamic stability of the Poiseuille flow in channels and pipes, Computing, 4 (1993), 3-6. (with R. Florea). 130. Stability of a binary mixture in a porous medium with Hall and ion-slip effect and Soret-Dufour currents, Analele Univ. din Oradea, 3 (1993), 92-96. (with L. Palese, D. Paşca). 131. Metode de determinare a curbei neutrale în stabilitatea Bénard, St. Cerc. Mec. Apl., 52, 4 (1993), 267-276. (with I. Oprea, D. Paşca). 132. Critical hydromagnetic stability of a thermodiffusive state, Rev. Roum. Math. Pures Appl., 38, 10 (1993), 831-840. (with L. Palese, D. Paşca, M. Buican). 133. Direcţii de cercetare principale în teoria sistemelor dinamice, Stud. Cerc. Mec .Apl., 52, 2 (1993), 153-171. 134. Bifurcation manifolds in multiparametric linear stability of continua, ZAMM, 73, 7-8 (1993), T831-T833. (with D. Paşca, S. Grădinaru, M. Gavrilescu). 135. Balance equations for the vector fields defined on orientable manifolds, Tensor (N. S.), 54 (1993), 88-90. (with C. Vamoş, N. Suciu).
1992
136. Aspecte ale modelării stratului limită al atmosferei. I, Stud. Cerc. Mec. Apl., 51, 1 (1992) ), 25-41. 137. Studiul calitativ al ecuaţiilor diferenţiale, St. Cerc. Mec. Apl., 51, 3 (1992), 317-326. 138. Models of asymptotic approximation governing the atmospheric motion over a low obstacle, Stud. Cerc. Mat., 44, 3 (1992), 237-252. (with G. Marinoschi). 139. Linear stability of a turbulent flow of Maxwell fluids in pipes, Rev. Roum. Math. Pures Appl., 37, 7 (1992), 579-586 (with C. Chiujdea, R. Florea). 140. Bifurcation problems in linear stability of continua, Quaderni. Dipto. di Mat. Univ. Bari, 1 (1992). (with I. Oprea). 141. Aspecte ale modelării stratului limită al atmosferei. II, Stud. Cerc. Mec. Apl., 51, 2 (1992). xix
1991
142. Modelarea matematică în mecanica fluidelor, St. Cerc. Mec. Apl., 50, 3-4 (1991), 295-298. 143. Efectul Toms, St. Cerc. Mec. Apl., 50, 5-6 (1991), 305-321. (with C. Chiujdea). 144. Linear stability of a turbulent flow of Maxwell fluids in pipes, Université de Metz, 21/1991, (with C. Chiujdea, R. Florea). 145. Evolution of the concept of asymptotic approximation, Noesis, 17 (1991), 45-50. 1990
146. Models of asymptotic approximation for synoptical flows, Zeitschrift für Meteorologie, 40, 1 (1990), 14-20. (with C. Vamoş). 147. Neutral stability curves for a thermal convection problem. II. The case of multiple solutions of the characteristic equation, Acta Mechanica, 81 (1990), 115-119. (with I. Oprea). 148. Metode numerice în teoria bifurcaţiei. III. Soluţii periodice, Stud. Cerc. Fiz., 42, 1 (1990), 117-125. (with I. Oprea). 149. Scenarii de turbulenţă în cadrul haosului determinist, Stud. Cerc. Mec. Apl., 49, 4 (1990), 413-417. (with C. Vamoş, N. Suciu). 150. Asimptote oblice din punctul de vedere al aproximaţiei asimptotice, Gaz. Mat. M, 2 (1990), 58-60. 151. On a hydrodynamic-social analogy, Rev. Roum. Philos. Logique, 34, 1-2 (1990), 100-102. 1989
152. Bifurcation manifolds in a multiparametric eigenvalue problem for linear hydromagnetic stability theory, Mathematica - Anal. Numér. Théor. Approx., 18, 2 (1989), 123-138 (with I. Oprea, C. Oprea). 153. Model de aproximaţie asimptotică de ordinul patru pentru ecuaţiile meteorologice primitive când numărul Rossby tinde la zero, St. Cerc. Meteorolgie, 3 (1989), 13-21 (with C. Vamoş). 154. Filtred equations as an asymptotic approximation model, Meteorology and Hydrology, 19, 2 (1989), 21-22 (with C. Vamoş). 155. Boundary layer separation I. Bubbles on leading edges, Rev. Roum. Sci. Tech.-Méc. Appl., 34, 5 (1989), 509-525, (with H. Dumitrescu, Al. Dumitrache). 156. Stabilitatea mişcării lichidelor pe un plan înclinat, Stud. Cerc. Mec. Apl., 48, 5 (1989). 157. Comparative study of the analytic methods used to solve problems in hydromagnetic stability theory, An. Univ. Bucureşti, Mat., 38, 1 (1989), 15- 20 (with A. Setelecan). 158. Lagrange and the calculus of variations, Noesis, 15 (1989), 29-35. xx
159. Fractalii şi unele aplicaţii ale lor, Stud. Cerc. Fiz., 41, 3 (1989), 269-288. 160. Suprafeţe neutrale bifurcate într-o problemă de inhibiţie a convecţiei termice datorită unui câmp magnetic, Stud. Cerc. Mec. Apl., 48, 3 (1989), 263-278. (with I. Oprea). 1988
161. Metode numerice în teoria bifurcaţiei. II. Soluţii staţionare în cazul infinit dimensional, Stud. Cerc. Fiz., 40, 1 (1988), 7-18. (with I. Oprea). 162. Bifurcation (catastrophe) surfaces in multiparametric eigenvalue problems in hydromagnetic stability theory, Bull. Inst. Politehn. Bucureşti, Ser. Construc. Maş, 50 (1988), 9-12 (with I. Oprea). 163. The bifurcation curve of characteristic equation provides the bifurcation point of the neutral curve of some elastic stability, Mathematica - Anal. Numér. Théor. Approx., 17, 2 (1988), 141-145 (with I. Oprea). 1987
164. Metode de rezolvare a unor probleme de valori proprii care apar în stabilitatea hidrodinamică liniară, St.Cerc.Fiz., 39, 1 (1987), 3-25 (with A. Setelecan). 165. Notă asupra unor probleme izoperimetrice în calculul elicei de randament maxim, St.Cerc.Mec.Apl., 46, 5 (1987), 478-482. 166. Exact solutions for some instability of Bénard type, Rev. Roum. Phys., 32, 4 (1987), 391-397. 1986
167. Proiectarea aerodinamică a elicei de randament maxim, St.Cerc. Mec.Apl., 45, 2 (1986), 129-141 (with H. Dumitrecsu, Al. Dumitrache). 168. Bifurcaţia stratului limită, BITNAV, 3 (1986), 148-149. 169. Metode numerice în teoria bifurcaţiei, Stud. Cerc. Fiz., 38, 10 (1986), 912- 924 (with I. Oprea). 1984
170. Echilibrul plasmei în sisteme toroidale şi stabilitatea sa macroscopică, St.Cerc.Fiz., 36, 1 (1984), 86-110. 1983
171. Neustanovivseesia ploscoe dvijenie tipa Puazeilia dlia jidkostei Rivlina- Eriksena, PMM, 47, 2 (1983), 342-344. (with S.S. Chetti). 172. Stabilitatea şi ramificarea în contextul sinergeticii, St.Cerc.Mec.Apl., 42, 2 (1983), 174-180. 1982
173. Bifurcation (catastrophe) surfaces for a problem in hydromagnetic stability, Rev. Roum. Math. Pures Appl., 27, 3 (1982), 335-337. xxi
174. Neutral stability curves for a thermal convection problem, Analele Univ. din Craiova, Secţia Mat. Fiz.-Chim., X (1982), 51-53 (with I. Oprea). 175. Characteristic equations for some eigenvalue problems in hydromagnetic stability theory, Mathematica, Cluj, 24 (47), 1-2 (1982), 31-41. 1981
176. On the nonexistence of regular solutions of a Blasius-like equation in the theory of the boundary layer of finite depth, Rev. Roum. Math. Pures Appl., 26, 6 (1981), 849-854 (with M. Moroianu). 177. . Recent results in fluid mechanics, Preprint 2, Univ. "Babeş - Bolyai", Fac. Mat., Cluj-Napoca, 1981. 178. On a Bénard convection in the presence of dielectrophoretic forces, J. Appl. Mech., 48 , 4 (1981), 980-981 (with O. Polotzka). 1980
179. Neutral stability curves for a thermal convection problem, Acta Mechanica, 37 (1980), 165-168 (with V. Cardoş). 1978
180. Bounds for linear characteristics of Couette and Poiseuille flows, Rev. Roum. Math. Pures et Appl., 23, 5 (1978), 707-720 (with Tr. Bădoiu). 1977
181. Stability of the Couette flow of a viscoelastic fluid. II, Rev. Roum. Math. Pures et Appl., 22, 9 (1977), 1223-1233. (with O. Polotzka). 182. Metode analitice în studiul fenomenologic al stabilităţii mişcării fluidelor vâscoase incompresibile descrise de soluţii generalizate ale ecuaţiilor Navier-Stokes, St.Cerc.Mat., 29, 6 (1977), 603-619. 1976
183. Universal criteria of hydrodynamic stability, Rev. Roum. Math. Pures et Appl., 21, 3 (1976), 287-302. 1973
184. Stability of the Couette flow of a viscoelastic fluid, Rev. Roum. Math. Pures et Appl., 18, 9 (1973), 1371-1374. 185. Teorema lui Squire pentru o mişcare într-un mediu poros, Petrol şi Gaze, 24, 11 (1973). 1972
186. Linear Couette flow stability for arbitrary gap between two rotating cylinders, Rev. Roum. Math. Pures et Appl., 17, 4 (1972). 187. Stability of spiral flow and of the flow in a curved channel, Rev. Roum. Math. Pures et Appl., 17, 3 (1972), 353-357. xxii
1971
188. Theorems of Joseph's type in hydrodynamic stability theory, Rev.Roum.Math.Pures et Appl., 16, 3 (1971), 355-362. 189. On the Kelvin-Helmholtz instability in presence of porous media, Rev. Roum. Math. Pures et Appl., 16, 1 (1971), 27-39 (with Şt. I. Gheorghiţă). 190. On the neutral stability of the Couette flow between two rotating cylinders, Rev.Roum.Math.Pures et Appl. 16, 4 (1971), 499-502. 191. Instability of two superposed liquids in a circular tube in the presence of a porous medium, Rev.Roum.Math.Pures et Appl., 16, 5 (1971), 677-680, (with Şt. I. Gheorghiţă). 1970
192. Note on Joseph’s inequalities in stability theory, ZAMP, 21, 2 (1970), 258- 260. 193. Sur la stabilité linéaire des mouvements plans des fluides, Comptes Rendus, Paris, Série A, 271 (1970), 559-561. 194. . Sufficient conditions for linear stability of two Ladyzhenskaya type fluids, Rev.Roum.Math.Pures et Appl., 15, 6 (1970), 819-823. 195. Contribuţii la studiul stabilităţii liniare a mişcării fluidelor, St. Cerc. Mat., 22, 9 (1970), 1247-1333. (The PhD Thesis of A. G.) 1969
196. Asupra soluţiilor asimptotice ale lui Heisenberg, St. Cerc. Mat., 21, 5 (1969), 747-750. 197. On a relationship between Heisenberg and Tollmien solutions, Rev.Roum.Math.Pures et Appl., 14, 7 (1969), 991-998. 198. Improvement of one of Joseph's theorems and one of its applications, Rev. Roum. Math. Pures et Appl., 14, 8 (1969), 1089-1092. 199. Criterii de stabilitate liniară a mişcării plan paralele a unui fluid nenewtonian, St.Cerc.Mat., 21, 7 (1969), 1027-1036. 200. On the generalized Tollmien solutions of the Rayleigh equation for a general velocity profile, Bull.Math. de la Soc. Sci.Math. de la R.S. de Roumanie, 13 (61), 2 (1969), 147-158. 1966
201. Corecţii de compresibilitate pentru profile von Mises, St.Cerc.Mat., 18, 2(1966), 301-308.
xxiii
ROMAI J., 6, 2(2010), 1–14
ON SPECIAL DIFFERENTIAL SUBORDINATIONS USING MULTIPLIER TRANSFORMATION AND RUSCHEWEYH DERIVATIVE Alina Alb Lupa¸s Department of Mathematics and Computer Science, University of Oradea, Romania [email protected] Abstract In the present paper we define a new operator using the multiplier transformation and α α Ruscheweyh derivative. Denote by RIn,λ,l the operator given by RIn,λ,l : A → A, α n n RIn,λ,l f (z) = (1−α)R f (z)+αI (n, λ, l) f (z), for z ∈ U, where R f (z) denote the Ruscheweyh derivative, I (n, λ, l) f (z) is the multiplier transformation and An = { f ∈ H(U): f (z) = n+1 z + an+1z + ..., z ∈ U} is the class of normalized analytic functions with A1 = A. A certain subclass, denoted by RIn (δ, λ, l, α) , of analytic functions in the open unit disc is introduced by means of the new operator. By making use of the concept of differ- ential subordination we will derive various properties and characteristics of the class RIn (δ, λ, l, α) . Also, several differential subordinations are established regardind the α operator RIn,λ,l. Keywords: differential subordination, convex function, best dominant, differential operator, Ruscheweyh derivative, multiplier transformations. 2000 MSC: 30C45, 30A20, 34A40. .
1. INTRODUCTION Denote by U the unit disc of the complex plane, U = {z ∈ C : |z| < 1} and H(U) the space of holomorphic functions in U. Let X∞ p j A (p, n) = { f ∈ H(U): f (z) = z + a jz , z ∈ U}, j=p+n with A (1, n) = An, A (1, 1) = A1 = A and
n n+1 H[a, n] = { f ∈ H(U): f (z) = a + anz + an+1z + ..., z ∈ U}, where p, n ∈ N, a ∈n C. o z f 00(z) Denote by K = f ∈ A : Re f 0(z) + 1 > 0, z ∈ U , the class of normalized convex functions in U. If f and g are analytic functions in U, we say that f is subordinate to g, written f ≺ g, if there is a function w analytic in U, with w(0) = 0, |w(z)| < 1, for all z ∈ U
1 2 Alina Alb Lupa¸s such that f (z) = g(w(z)) for all z ∈ U. If g is univalent, then f ≺ g if and only if f (0) = g(0) and f (U) ⊆ g(U). Let ψ : C3 × U → C and h an univalent function in U. If p is analytic in U and satisfies the (second-order) differential subordination
ψ(p(z), zp0(z), z2 p00(z); z) ≺ h(z), for z ∈ U, (1) then p is called a solution of the differential subordination. The univalent function q is called a dominant of the solutions of the differential subordination, or more simply a dominant, if p ≺ q for all p satisfying (1). A dominant eq that satisfies eq ≺ q for all dominants q of (1) is said to be the best dominant of (1). The best dominant is unique up to a rotation of U. Definition 1.1. (Ruscheweyh [3]) For f ∈ A, n ∈ N, the operator Rn is defined by Rn : A → A,
R0 f (z) = f (z) R1 f (z) = z f 0 (z) ... (n + 1) Rn+1 f (z) = z Rn f (z) 0 + nRn f (z) , for z ∈ U.
P∞ j n P∞ n j Remark 1.1. If f ∈ A, f (z) = z + j=2 a jz , then R f (z) = z + j=2 Cn+ j−1a jz , for z ∈ U. Definition 1.2. [6] For f ∈ A(p, n), p, n ∈ N, m ∈ N∪ {0}, λ, l ≥ 0, the operator Ip (m, λ, l) f (z) is defined by the following infinite series ! X∞ p + λ ( j − 1) + l m I (m, λ, l) f (z):= zp + a z j. p p + l j j=p+n
Remark 1.2. It follows from the above definition that
Ip (0, λ, l) f (z) = f (z),
0 (p + l) Ip (m + 1, λ, l) f (z) = p(1 − λ) + l Ip (m, λ, l) f (z) + λz Ip (m, λ, l) f (z) , for z ∈ U.
Remark 1.3. If p = 1, n = 1, we have A(1, 1) = A1 = A, I1 (m, λ, l) f (z) = I (m, λ, l) and
(l + 1) I (m + 1, λ, l) f (z) = [l + 1 − λ] I (m, λ, l) f (z) + λz (I (m, λ, l) f (z))0 ,
for z ∈ U. On special differential subordinations using multiplier transformation... 3
P∞ j Remark 1.4. If f ∈ A, f (z) = z + j=2 a jz , then ! X∞ 1 + λ ( j − 1) + l m I (m, λ, l) f (z) = z + a z j, l + 1 j j=2
for z ∈ U.
m Remark 1.5. For l = 0, λ ≥ 0, the operator Dλ = I (m, λ, 0) was introduced and studied by Al-Oboudi [5], which reduced to the S˘al˘agean differential operator S m = m I (m, 1, 0) [4] for λ = 1. The operator Il = I (m, 1, l) was studied recently by Cho and Srivastava [7] and Cho and Kim [8]. The operator Im = I (m, 1, 1) was studied δ by Uralegaddi and Somanatha [12], the operator Dλ = I (δ, λ, 0), with δ ∈ R, δ ≥ 0, was introduced by Acu and Owa [1]. Lemma 1.1. (Hallenbeck and Ruscheweyh [1, Th. 3.1.6, p. 71]) Let h be a convex function with h(0) = a, and let γ ∈ C\{0} be a complex number with Re γ ≥ 0. If p ∈ H[a, n] and 1 p(z) + zp0(z) ≺ h(z), for z ∈ U, γ then p(z) ≺ g(z) ≺ h(z), for z ∈ U, R γ z γ/n−1 where g(z) = nzγ/n 0 h(t)t dt, for z ∈ U. Lemma 1.2. (Miller and Mocanu [1]) Let g be a convex function in U and let h(z) = g(z) + nαzg0(z), for z ∈ U, where α > 0 and n is a positive integer. n n+1 If p(z) = g(0) + pnz + pn+1z + ..., for z ∈ U is holomorphic in U and
p(z) + αzp0(z) ≺ h(z), for z ∈ U, then p(z) ≺ g(z), for z ∈ U, and this result is sharp.
2. MAIN RESULTS α Definition 2.1. Let α, λ, l ≥ 0, n ∈ N. Denote by RIn,λ,l the operator given by α RIn,λ,l : A → A,
α n RIn,λ,l f (z) = (1 − α)R f (z) + αI (n, λ, l) f (z), for z ∈ U. P Remark 2.1. If f ∈ A, f (z) = z + ∞ a z j, then n j=2 j o α P∞ 1+λ( j−1)+l n n j RIn,λ,l f (z) = z + j=2 α l+1 + (1 − α) Cn+ j−1 a jz , for z ∈ U. 4 Alina Alb Lupa¸s
0 n Remark 2.2. For α = 0, RIn,λ,l f (z) = R f (z), where z ∈ U and for α = 1, 1 RIn,λ,l f (z) = I (n, λ, l) f (z), where z ∈ U, which was studied in [3]. α n For l = 0, we obtain RIn,λ,0 f (z) = RD1,α f (z) which was studied in [4] and for α n l = 0 and λ = 1, we obtain RIn,1,0 f (z) = Lα f (z) which was studied in [2]. α 0 0 For n = 0, RI0,λ,l f (z) = (1 − α) R f (z) + αI (0, λ, l) f (z) = f (z) = R f (z) = I (0, λ, l) f (z), where z ∈ U. Definition 2.2. Let δ ∈ [0, 1), α, λ, l ≥ 0 and n ∈ N. A function f ∈ A is said to be in the class RIn (δ, λ, l, α) if it satisfies the inequality α 0 Re RIn,λ,l f (z) > δ, for z ∈ U. (2)
Theorem 2.1. The set RIn (δ, λ, l, α) is convex. Proof. Let the functions X∞ j f j (z) = z + a jkz , for k = 1, 2, z ∈ U, j=2 be in the class RIn (δ, λ, l, α). It is sufficient to show that the function
h (z) = η1 f1 (z) + η2 f2 (z) is in the class RIn (δ, λ, l, α) with η1 andη2 nonnegative such that η1 + η2 = 1. P∞ j Since h (z) = z + j=2 η1a j1 + η2a j2 z , for z ∈ U, then ( ! ) X∞ 1 + λ ( j − 1) + l n RIα h (z) = z+ α + (1 − α) Cn η a + η a z j, for z ∈ U. n,λ,l l + 1 n+ j−1 1 j1 2 j2 j=2 (3) Differentiating (3) we obtain
∞ ( !n ) 0 X 1 + λ ( j − 1) + l RIα h (z) = 1+ α + (1 − α) Cn η a + η a jz j−1, n,λ,l l + 1 n+ j−1 1 j1 2 j2 j=2 for z ∈ U. Hence n o α 0 P∞ 1+λ( j−1)+l n n j−1 Re RIn,λ,lh (z) = 1 + Re η1 j=2 j α l+1 + (1 − α) Cn+ j−1 a j1z + n o P∞ 1+λ( j−1)+l n n j−1 +Re η2 j=2 j α l+1 + (1 − α) Cn+ j−1 a j2z . (4) Taking into account that f1, f2 ∈ RIn (δ, λ, l, α) we deduce ( ! ) X∞ 1 + λ ( j − 1) + l n Re η j α + (1 − α) Cn a z j−1 > η (δ − 1) , for k = 1, 2. k l + 1 n+ j−1 jk k j=2 (5) On special differential subordinations using multiplier transformation... 5
Using (5) we get from (4) α 0 Re RIn,λ,lh (z) > 1 + η1 (δ − 1) + η2 (δ − 1) = δ, for z ∈ U, which is equivalent that RIn (δ, λ, l, α) is convex. 1 0 Theorem 2.2. Let g be a convex function in U and let h (z) = g (z)+ c+2 zg (z) , where z ∈ U, c > 0. R z ( ) ( ) ( )( ) c+2 c ( ) If f ∈ RIn δ, λ, l, α and F z = Ic f z = zc+1 0 t f t dt, for z ∈ U, then α 0 RIn,λ,l f (z) ≺ h (z) , for z ∈ U, (6) implies α 0 RIn,λ,lF (z) ≺ g (z) , for z ∈ U, and this result is sharp. Proof. We obtain that Z z zc+1F (z) = (c + 2) tc f (t) dt. (7) 0 Differentiating (7), with respect to z, we have (c + 1) F (z) + zF0 (z) = (c + 2) f (z) and α α 0 α (c + 1) RIn,λ,lF (z) + z RIn,λ,lF (z) = (c + 2) RIn,λ,l f (z) , for z ∈ U. (8) Differentiating (8) we have
0 1 00 0 RIα F (z) + z RIα F (z) = RIα f (z) , for z ∈ U. (9) n,λ,l c + 2 n,λ,l n,λ,l Using (9), the differential subordination (6) becomes
0 1 00 1 RIα F (z) + z RIα F (z) ≺ g (z) + zg0 (z) . (10) n,λ,l c + 2 n,λ,l c + 2 If we denote α 0 p (z) = RIn,λ,lF (z) , for z ∈ U, (11) then p ∈ H [1, 1] . Replacing (11) in (10) we obtain 1 1 p (z) + zp0 (z) ≺ g (z) + zg0 (z) , for z ∈ U. c + 2 c + 2 Using Lemma 1.2 we have α 0 p (z) ≺ g (z) , for z ∈ U, i.e. RIn,λ,lF (z) ≺ g (z) , for z ∈ U, 6 Alina Alb Lupa¸s and g is the best dominant. Example 2.1. If f ∈ RD 1, 1, 0, 1 , then f 0 (z) + z f 00 (z) ≺ 3−2z implies F0 (z) + 1 2 3(1−z)2 R z 00 ( ) 1 ( ) 3 ( ) zF z ≺ 1−z , where F z = z2 0 t f t dt. 1+(2δ−1)z Theorem 2.3. Let h (z) = 1+z , where δ ∈ [0, 1) and c > 0. R z ( )( ) c+2 c ( ) If α, λ, l ≥ 0, n ∈ N and Ic f z = zc+1 0 t f t dt, for z ∈ U, then ∗ Ic [RIn (δ, λ, l, α)] ⊂ RIn δ , λ, l, α , (12)
R 1 c+1 ∗ ( )( ) t where δ = 2δ − 1 + c + 2 2 − 2δ 0 t+1 dt. Proof. The function h is convex and using the same steps as in the proof of Theorem 2.2 we get from the hypothesis of Theorem 2.3 that 1 p (z) + zp0 (z) ≺ h (z) , c + 2 where p (z) is defined in (11). Using Lemma 1.1 we deduce that α 0 p (z) ≺ g (z) ≺ h (z) , i.e. RIn,λ,lF (z) ≺ g (z) ≺ h (z) , where Z Z c + 2 z 1 + (2δ − 1) t (c + 2)(2 − 2δ) z tc+1 g (z) = tc+1 dt = 2δ − 1 + dt. c+2 c+2 z 0 1 + t z 0 t + 1 Since g is convex and g (U) is symmetric with respect to the real axis, we deduce
Z 1 c+1 α 0 ∗ t Re RIn,λ,lF (z) ≥ minRe g (z) = Re g (1) = δ = 2δ−1+(c + 2)(2 − 2δ) dt. |z|=1 0 t + 1 (13) From (13) we deduce inclusion (12). Theorem 2.4. Let g be a convex function, g(0) = 1 and let h be the function h(z) = g(z) + zg0(z), for z ∈ U. If α, λ, l ≥ 0, n ∈ N, f ∈ A and satisfies the differential subordination α 0 RIn,λ,l f (z) ≺ h(z), for z ∈ U, (14) then RIα f (z) n,λ,l ≺ g(z), for z ∈ U, z and this result is sharp. On special differential subordinations using multiplier transformation... 7
α Proof. By using the properties of operator RIn,λ,l, we have ( ! ) X∞ 1 + λ ( j − 1) + l n RIα f (z) = z + α + (1 − α) Cn a z j, for z ∈ U. n,λ,l l + 1 n+ j−1 j j=2 Consider X∞ ( !n ) 1 + λ ( j − 1) + l n j z + α + (1 − α) C a jz α l + 1 n+ j−1 RI f (z) j=2 p(z) = n,λ,l = = z z 2 = 1 + p1z + p2z + ..., for z ∈ U. We deduce that p ∈ H[1, 1]. α Let RIn,λ,l f (z) = zp(z), for z ∈ U. By differentiating we obtain α 0 0 RIn,λ,l f (z) = p(z) + zp (z), for z ∈ U. Then (14) becomes p(z) + zp0(z) ≺ h(z) = g(z) + zg0(z), for z ∈ U. By using Lemma 1.2, we have RIα f (z) p(z) ≺ g(z), for z ∈ U, i.e. n,λ,l ≺ g(z), for z ∈ U. z
Theorem 2.5. Let h be an holomorphic function which satisfies the inequality zh00(z) 1 Re 1 + h0(z) > − 2 , for z ∈ U, and h(0) = 1. If α, λ, l ≥ 0, n ∈ N, f ∈ A and satisfies the differential subordination α 0 RIn,λ,l f (z) ≺ h(z), for z ∈ U, (15) then RIα f (z) n,λ,l ≺ q(z), for z ∈ U, z R 1 z where q(z) = z 0 h(t)dt. The function q is convex and it is the best dominant. Proof. Let n n o α P∞ 1+λ( j−1)+l n j RI f (z) z + j=2 α l+1 + (1 − α) Cn+ j−1 a jz p(z) = n,λ,l = z z ( ! ) X∞ 1 + λ ( j − 1) + l n X∞ = 1 + α + (1 − α) Cn a z j−1 = 1 + p z j−1, l + 1 n+ j−1 j j j=2 j=2 8 Alina Alb Lupa¸s for z ∈ U, p ∈ H[1, 1]. α 0 0 Differentiating, we obtain RIn,λ,l f (z) = p(z)+zp (z), for z ∈ U and (15) becomes p(z) + zp0(z) ≺ h(z), for z ∈ U.
Using Lemma 1.1, we have
α Z RI f (z) 1 z p(z) ≺ q(z), for z ∈ U, i.e. n,λ,l ≺ q(z) = h(t)dt, for z ∈ U, z z 0 and q is the best dominant.
Theorem 2.6. Let g be a convex function such that g (0) = 1 and let h be the function h (z) = g (z) + zg0 (z), for z ∈ U. If α, λ, l ≥ 0, n ∈ N, f ∈ A and the differential subordination zRIα f (z)0 n+1,λ,l α ≺ h (z) , for z ∈ U (16) RIn,λ,l f (z) holds, then α RIn+1,λ,l f (z) α ≺ g (z) , for z ∈ U, RIn,λ,l f (z) and this result is sharp. P Proof. For f ∈ A, f (z) = z + ∞ a z j we have n j=2 j o α P∞ 1+λ( j−1)+l n n j RIn,λ,l f (z) = z + j=2 α l+1 + (1 − α) Cn+ j−1 a jz , for z ∈ U. Consider P∞ 1+λ( j−1)+l n+1 n+1 j α z + α + (1 − α) C a jz RIn+1,λ,l f (z) j=2 l+1 n+ j p(z) = α = P n n o . RI f (z) ∞ 1+λ( j−1)+l n j n,λ,l z + j=2 α l+1 + (1 − α) Cn+ j−1 a jz We have 0 0 RIα f (z) RIα f (z) 0 n+1,λ,l n,λ,l p (z) = α − p (z) · α RIn,λ,l f (z) RIn,λ,l f (z) and we obtain zRIα f (z)0 0 n+1,λ,l p (z) + z · p (z) = α RIn,λ,l f (z) . Relation (16) becomes
p(z) + zp0(z) ≺ h(z) = g(z) + zg0(z), for z ∈ U. On special differential subordinations using multiplier transformation... 9
By using Lemma 1.2, we have
α RIn+1,λ,l f (z) p(z) ≺ g(z), for z ∈ U, i.e. α ≺ g(z), for z ∈ U. RIn,λ,l f (z)
Theorem 2.7. Let g be a convex function such that g(0) = 0 and let h be the function h(z) = g(z) + zg0(z), for z ∈ U. If α, λ, l ≥ 0, n ∈ N, f ∈ A and the differential subordination
α α (n + 1) RIn+1,λ,l f (z) − (n − 1) RIn,λ,l f (z) − ! l + 1 α n + 1 − I (n + 1, λ, l) f (z) − I (n, λ, l) f (z) ≺ h(z), for z ∈ U (17) λ holds, then α RIn,λ,l f (z) ≺ g(z), for z ∈ U. This result is sharp. Proof. Let
α n p(z) = RIn,λ,α f (z) = (1 − α)R f (z) + αI (n, λ, l) f (z) = (18) ( ! ) X∞ 1 + λ ( j − 1) + l n = z + α + (1 − α) Cn a z j = p z + p z2 + .... l + 1 n+ j−1 j 1 2 j=2 We deduce that p ∈ H[0, 1]. α n By using the properties of operators RIn,λ,l, R and I (n, λ, l), after a short calcula- tion, we obtain p (z) + zp0 (z) = (n + 1) RIα f (z) − (n − 1) RIα f (z) − n+1,λ,l n,λ,l l+1 −α n + 1 − λ [I (n + 1, λ, l) f (z) − I (n, λ, l) f (z)]. Using the notation in (18), the differential subordination becomes
p(z) + zp0(z) ≺ h(z) = g(z) + zg0(z).
By using Lemma 1.2, we have
α p(z) ≺ g(z), for z ∈ U, i.e. RIn,λ,l f (z) ≺ g(z), for z ∈ U, and this result is sharp.
1+(2β−1)z Theorem 2.8. Let h(z) = 1+z be a convex function in U, where 0 ≤ β < 1. 10 Alina Alb Lupa¸s
If α, λ, l ≥ 0, n ∈ N, f ∈ A and satisfies the differential subordination
α α (n + 1) RIn+1,λ,l f (z) − (n − 1) RIn,λ,l f (z) − ! l + 1 −α n + 1 − I (n + 1, λ, l) f (z) − I (n, λ, l) f (z) ≺ h(z), for z ∈ U, (19) λ then α RIn,λ,l f (z) ≺ q(z), for z ∈ U, ln(1+z) where q is given by q(z) = 2β − 1 + 2(1 − β) z , for z ∈ U. The function q is convex and it is the best dominant. Proof. Following the same steps as in the proof of Theorem 2.7 and considering α p(z) = RIn,λ,l f (z), the differential subordination (19) becomes 1 + (2β − 1)z p(z) + zp0(z) ≺ h(z) = , for z ∈ U. 1 + z By using Lemma 1.1 for γ = 1 and n = 1, we have p(z) ≺ q(z), i.e.,
Z z Z z α 1 1 1 + (2β − 1)t 1 RIn,λ,l f (z) ≺ q(z) = h(t)dt = dt = 2β−1+2(1−β) ln(z+1), z 0 z 0 1 + t z for z ∈ U.
Theoremh 2.9.i Let h be an holomorphic function which satisfies the inequality zh00(z) 1 Re 1 + h0(z) > − 2 , for z ∈ U, and h (0) = 0. If α, λ, l ≥ 0, n ∈ N, f ∈ A and satisfies the differential subordination
α α (n + 1) RIn+1,λ,l f (z) − (n − 1) RIn,λ,l f (z) − ! l + 1 α n + 1 − I (n + 1, λ, l) f (z) − I (n, λ, l) f (z) ≺ h(z), for z ∈ U, (20) λ then α RIn,λ,l f (z) ≺ q(z), for z ∈ U, R 1 z where q is given by q(z) = z 0 h(t)dt. The function q is convex and it is the best dominant.
α α Proof. Using the properties of operator RIn,λ,l and considering p (z) = RIn,λ,l f (z), we obtain 0 α α p(z) + zp (z) = (n + 1) RIn+1,λ,l f (z) − (n − 1) RIn,λ,l f (z) − ! l + 1 α n + 1 − I (n + 1, λ, l) f (z) − I (n, λ, l) f (z) , for z ∈ U. λ On special differential subordinations using multiplier transformation... 11
Then (20) becomes p(z) + zp0(z) ≺ h(z), for z ∈ U. Since p ∈ H[0, 1], using Lemma 1.1, we deduce Z z α 1 p(z) ≺ q(z), for z ∈ U, i.e. RIn,λ,l f (z) ≺ q(z) = h(t)dt, for z ∈ U, z 0 and q is the best dominant. Theorem 2.10. Let g be a convex function such that g(0) = 1 and let h be the function h(z) = g(z) + zg0(z), for z ∈ U. If α, λ, l ≥ 0, n ∈ N, f ∈ A and the differential subordination (n + 1)(n + 2) (n + 1)(2n + 1) n2 RIα f (z) − RIα f (z) + RIα f (z) + z n+2,λ,l z n+1,λ,l z n,λ,l " # α (l + 1)2 − (n + 1)(n + 2) I (n + 2, λ, l) f (z) − z λ2 ¯ α 2 (l + 1 − λ) l + 1 − (n + 1)(2n + 1) I (n + 1, λ, l) f (z) + z λ2 " # α (l + 1 − λ)2 − n2 I (n, λ, l) f (z) ≺ h(z), for z ∈ U (21) z λ2 holds, then α 0 [RIn,λ,l f (z)] ≺ g(z), for z ∈ U. This result is sharp. Proof. Let α 0 n 0 0 p(z) = RIn,λ,l f (z) = (1 − α) R f (z) + α (I (n, λ, l) f (z)) (22) ( ! ) X∞ 1 + λ ( j − 1) + l n = 1 + α + (1 − α) Cn ja z j−1 = 1 + p z + p z2 + .... l + 1 n+ j−1 j 1 2 j=2 We deduce that p ∈ H[1, 1]. α n By using the properties of operators RIn,λ,l, R and I (n, λ, l), after a short calcula- tion, we obtain
p (z) + zp0 (z) = (n+1)(n+2) RIα f (z) − (n+1)(2n+1) RIα f (z) + z n+2,λ,l z n+1,λ,l 2 2 + n RIα f (z) + α (l+1) − (n + 1)(n + 2) I (n + 2, λ, l) f (z) − z n,λ,l z λ2 2(l+1−λ)(l¯+1) − α − (n + 1)(2n + 1) I (n + 1, λ, l) f (z) + z λ2 2 + α (l+1−λ) − n2 I (n, λ, l) f (z) . z λ2 12 Alina Alb Lupa¸s
Using the notation in (22), the differential subordination becomes
p(z) + zp0(z) ≺ h(z) = g(z) + zg0(z).
By using Lemma 1.2, we have α 0 p(z) ≺ g(z), for z ∈ U, i.e. RIn,λ,l f (z) ≺ g(z), for z ∈ U, and this result is sharp.
Example 2.2. If n = 1, α = 1, λ = 1, l = 0, f ∈ A, we deduce that f 0(z) + 3z f 00(z) + z2 f 000(z) ≺ g(z) + zg0(z), which yields that f 0(z) + z f 00(z) ≺ g(z), for z ∈ U. 1+(2β−1)z Theorem 2.11. Let h(z) = 1+z be a convex function in U, where 0 ≤ β < 1. If α, λ, l ≥ 0, n ∈ N, f ∈ A and satisfies the differential subordination
(n + 1)(n + 2) (n + 1)(2n + 1) n2 RIα f (z) − RIα f (z) + RIα f (z) + z n+2,λ,l z n+1,λ,l z n,λ,l " # α (l + 1)2 − (n + 1)(n + 2) I (n + 2, λ, l) f (z) − z λ2 ¯ α 2 (l + 1 − λ) l + 1 − (n + 1)(2n + 1) I (n + 1, λ, l) f (z) + z λ2 " # α (l + 1 − λ)2 − n2 I (n, λ, l) f (z) ≺ h(z), for z ∈ U, (23) z λ2 then α 0 RIn,λ,l f (z) ≺ q(z), for z ∈ U, ln(1+z) where q is given by q(z) = 2β − 1 + 2(1 − β) z , for z ∈ U. The function q is convex and it is the best dominant.
Proof. Following the same steps as in the proof of Theorem 2.10 and considering α 0 p(z) = RIn,λ,l f (z) , the differential subordination (23) becomes 1 + (2β − 1)z p(z) + zp0(z) ≺ h(z) = , for z ∈ U. 1 + z By using Lemma 1.1 for γ = 1 and n = 1, we have p(z) ≺ q(z), i.e., Z z Z z α 0 1 1 1 + (2β − 1)t 1 RIn,λ,l f (z) ≺ q(z) = h(t)dt = dt = 2β−1+2(1−β) ln(z+1), z 0 z 0 1 + t z for z ∈ U. On special differential subordinations using multiplier transformation... 13
Theoremh 2.12.i Let h be an holomorphic function which satisfies the inequality zh00(z) 1 Re 1 + h0(z) > − 2 , for z ∈ U, and h (0) = 1. If α, λ, l ≥ 0, n ∈ N, f ∈ A and satisfies the differential subordination (n + 1)(n + 2) (n + 1)(2n + 1) n2 RIα f (z) − RIα f (z) + RIα f (z) + z n+2,λ,l z n+1,λ,l z n,λ,l " # α (l + 1)2 − (n + 1)(n + 2) I (n + 2, λ, l) f (z) − z λ2 ¯ α 2 (l + 1 − λ) l + 1 − (n + 1)(2n + 1) I (n + 1, λ, l) f (z) + z λ2 " # α (l + 1 − λ)2 − n2 I (n, λ, l) f (z) ≺ h(z), for z ∈ U, (24) z λ2 then α 0 RIn,λ,l f (z) ≺ q(z), for z ∈ U, R 1 z where q is given by q(z) = z 0 h(t)dt. The function q is convex and it is the best dominant. α α 0 Proof. Using the properties of operator RIn,λ,l and considering p (z) = RIn,λ,l f (z) , we obtain (n + 1)(n + 2) (n + 1)(2n + 1) n2 p(z)+zp0(z) = RIα f (z)− RIα f (z)+ RIα f (z) + z n+2,λ,l z n+1,λ,l z n,λ,l " # α (l + 1)2 − (n + 1)(n + 2) I (n + 2, λ, l) f (z) − z λ2 ¯ α 2 (l + 1 − λ) l + 1 − (n + 1)(2n + 1) I (n + 1, λ, l) f (z) + z λ2 " # α (l + 1 − λ)2 − n2 I (n, λ, l) f (z) , for z ∈ U. z λ2 Then (24) becomes
p(z) + zp0(z) ≺ h(z), for z ∈ U.
Since p ∈ H[1, 1], using Lemma 1.1, we deduce Z z α 0 1 p(z) ≺ q(z), for z ∈ U, i.e. RIn,λ,l f (z) ≺ q(z) = h(t)dt, for z ∈ U, z 0 and q is the best dominant. 14 Alina Alb Lupa¸s References
[1] M. Acu, S. Owa, Note on a class of starlike functions, RIMS, Kyoto, 2006. [2] A. Alb Lupas¸, On special differential subordinations using S˘al˘agean and Ruscheweyh operators, Math. Inequ. Appl., 12, 4(2009), 781-790. [3] A. Alb Lupas¸, A special comprehensive class of analytic functions defined by multiplier transfor- mation, J. Comput. Anal. Appl., 12, 2(2010), 387-395. [4] A. Alb Lupas¸, On special differential subordinations using a generalized S˘al˘agean operator and Ruscheweyh derivative, J. Comput. Anal. Appl., 2011 (accepted). [5] F.M. Al-Oboudi, On univalent functions defined by a generalized S˘al˘agean operator, Ind. J. Math. Math. Sci., 27(2004), 1429-1436. [6] A. Catas¸,˘ On certain class of p-valent functions defined by new multiplier transformations, Pro- ceedings Book of the International Symposium on Geometric Function Theory and Applications, August 20-24, 2007, TC Istanbul Kultur University, Turkey, 241-250. [7] N.E. Cho, H.M. Srivastava, Argument estimates of certain analytic functions defined by a class of multiplier transformations, Math. Comput. Modelling, 37, 1-2(2003), 39-49. [8] N.E. Cho, T.H. Kim, Multiplier transformations and strongly close-to-close functions, Bull. Ko- rean Math. Soc., 40, 3(2003), 399-410. [9] S.S. Miller, P.T. Mocanu, Differential Subordinations. Theory and Applications, Marcel Dekker Inc., New York, 2000. [10] St. Ruscheweyh, New criteria for univalent functions, Proc. Amet. Math. Soc., 49(1975), 109- 115. [11] G. St. Sal˘ agean,˘ Subclasses of univalent functions, Lecture Notes in Math., Springer, Berlin, 1013(1983), 362-372. [12] B.A. Uralegaddi, C. Somanatha, Certain classes of univalent functions, Current topics in analytic function theory, World. Sci. Publishing, River Edge, N.Y., (1992), 371-374. ROMAI J., 6, 2(2010), 15–26
A PRODUCT FORMULA APPROACH TO AN INVERSE PROBLEM GOVERNED BY NONLINEAR PHASE-FIELD TRANSITION SYSTEM. CASE 1D Tommaso Benincasa1, Costic˘aMoro¸sanu2 1University of Bologna, Italy 2University ”Alexandru Ioan Cuza” of Ia¸si,Romania Abstract In this paper we study an inverse problem, in one space dimension case, connected with the industrial solidification process called casting wire, as an optimal control prob- lem governed by nonlinear phase-field system with nonhomogeneous Cauchy-Neumann boundary conditions. We prove the convergence of an iterative scheme of fractional steps type for the optimal control problem. Moreover, necessary optimality conditions are established for the approximating process. The advantage of such approach leads to the idea of a numerical algorithm in order to approximate the original optimal control problem.
Keywords: nonlinear parabolic systems, phase-field models, optimality conditions, applications, time- dependent initial-boundary value problems, fractional steps method, inverse problems. 2000 MSC: 35K55, 49N15, 62P30, 65M12, 65M32.
1. INTRODUCTION; SETTING OF THE PROBLEM Phase-field models, strongly studied in recent years, describe the phase transitions between two different phases in a pure material by a system of nonlinear parabolic equations. These models can be viewed as extensions of the classical Stefan problem in two phases. Consequently, the interface boundary between the phases can be con- structed from the so-called phase function and, phenomena associated with surface tension and supercooling are incorporated into the model. The mathematical literature concerning the optimal control problems associated with such models is in a deep process of development, as the models are suitable for many modern applications. One of them is the subject of this paper.
1.1. PROBLEM FORMULATION
Denote by Ω = (0, b1) ⊂ R, 0 < b1 < +∞. Let T > 0 and we set: Q = (0, T) × Ω, Σ0 = {(t, x) ∈ Q, t = f (x)}, Σ = (0, T) × {b1},
15 16 Tommaso Benincasa, Costic˘aMoro¸sanu where t = f (x) is considered to be the equation of the moving boundary separating the liquid and solid phases, 0 = f (b0), 0 < b0 < b1 (see Figure 1.1). Consider the following nonlinear parabolic system in one space dimension: ` ρc ut + ϕt = kuxx, in Q, 2 (1.1) 1 τϕ = ξ2ϕ + (ϕ − ϕ3) + 2u, in Q, t xx 2a subject to the non-homogeneous Cauchy-Neumann boundary conditions
ux + hu = w(t), ϕx = 0 on Σ, (1.2)
ux = 0, ϕx = 0 on Σ0, (1.3) and to the initial conditions
u(0, x) = u0(x), ϕ(0, x) = ϕ0(x) on Ω0 = [b0, b1], (1.4) where u is the reduced temperature distribution, ϕ is the phase function used to distin- guish between the phase of Ω, u0, ϕ0 : Ω −→ R are given functions, w : [0, T] −→ R is the boundary control (the temperature surrounding at x = b1), w ∈ U where U = {v ∈ L∞([0, T]), −R ≤ v(t) ≤ 0 a.e. t ∈ [0, T]}; the positive pa- rameters ρ, c, τ, ξ, `, k, h, a, have the following physical meaning: ρ - is the density, c - is the heat capacity, τ - is the relaxation time, ξ - is the length scale of the interface, ` - denotes the latent heat, k - the heat conductivity, h - the heat transfer coefficient and a is an probabilistic measure on the individual atoms (a depends on ξ).
Figure 1.1.Geometrical image of the elements in inverse problem (Pinv). The mathematical model (1.1), introduced by Caginalp [3], has been established in the literature as an extension of the classical two phase Stefan problem to capture the effects of surface tension, supercooling, and superheating. A product formula approach to an inverse problem... 17
As regards the existence, it is known that under appropriate conditions on u0, ϕ0 2,1 ∞ and w, the state system (1.1)-(1.4) has a unique solution u, ϕ ∈ W = Wp (Q)∩L (Q), p > 3/2 (see Proposition 2.1 in [6]). Given the positive numbers d1, d2, we define: the pure liquid region : {(t, x) ∈ Q, u(t, x) > d2 and ϕ(t, x) ≥ 1 + d1}, the pure solid region : {(t, x) ∈ Q, u(t, x) < −d2 and ϕ(t, x) ≤ −1 − d1}, the separating region : {(t, x) ∈ Q, |u(t, x)| < d2, |ϕ(t, x)| ≤ 1 + d1}. and we set Q0 = (t, x) ∈ Q, f (x) ≤ t ≤ T . Consider the following inverse problem: ∞ (Pinv) Given Σ0 f ind the boundary control w ∈ L ([0, T]) such that Q0 is in the liquid region, Q1 = Q \ Q0 is in the solid region and a neighbourhood of Σ0 is the separating region between the liquid and the solid region.
(Pinv) is in general ”ill posed” and a common way to treat this inverse problem is to reformulate it as an optimal control problem with an appropriate cost functional. Consequently, we will concern in the present paper with an optimal control problem associated to the inverse problemZ (Pinv), namely: β (P) Minimize L (w) = (u(t, x) − δ )+ 2 · χ dtdx+ 0 2 2 Q0 Q Z ZT 1 1 + (ϕ(t, x) − 1 − δ )2 · χ dtdx + w2(t) dt, 2 1 Q0 2 Q 0 on all (u, ϕ) solution of the system (1.1)-(1.4) and for all w ∈ U β > 0 is a given constant. In the above statement we denoted by u+ the positive part of u, i.e. ( u, if u > 0, u+ = 0, if u ≤ 0. We point out that problem (P) is an optimal problem with boundary control w(t) depending on time variable t ∈ [0, T], being dictated by the industrial solidification process like casting wire.
1.2. APPROXIMATING PROCESS We associate to the nonlinear system (1.1)-(1.4) the following approximanting scheme (ε > 0): ` ρcuε + ϕε = kuε t 2 t xx ε in Q0 = (t, x) ∈ Q, ε ≤ t ≤ T , (1.5) 1 τϕε = ξ2ϕε + ϕε + 2uε, t xx 2a 18 Tommaso Benincasa, Costic˘aMoro¸sanu
ε ε ε ε ux + hu = w(t), ϕx = 0 on Σ = [ε, T] × {b1}, (1.6) ε ε ε ux = 0, ϕx = 0 on Σ0 = {(t, x) ∈ Q, ε ≤ t ≤ T}, (1.7) ε ε ε u (ε, x) = u0(x) ϕ+(ε, x) = z(ε, ϕ−(ε, x)) on Ωε. (1.8) ε where z(ε, ϕ−(ε, x)) is the solution of the Cauchy problem: 1 z0(s) + z3(s) = 0, s ∈ (0, ε), 2a (1.9) ε ε z(0) = ϕ−(ε, x), ϕ−(0, x) = ϕ0(x), ε ε ε ε and ϕ+(ε, x) = limt↓ε ϕ (t, x), ϕ−(ε, x) = limt↑ε ϕ (t, x). The convergence and weak stability of the approximating scheme (1.5)-(1.9), in a more general case (w(t, x) in place of w(t)), was studied in the paper [2]. Corresponding to the approximating scheme (1.5)-(1.9), we will consider the ap- proximating optimal control problem: Z β (Pε) Minimize Lε(w) = (uε(t, x) − δ )+ 2 · χ dtdx+ 0 2 2 Q0 Q Z ZT 1 1 + (ϕε(t, x) − 1 − δ )2 · χ dtdx + w2(t) dt, 2 1 Q0 2 Q 0 on all (uε, ϕε) solution of (1.5)-(1.9) corresponding to w ∈ U. The main result of the present paper (Theorem 2.1) says that problem (P) can be approximated for ε → 0 by the sequence of optimal control problems (Pε) and so the computation of the approximate boundary control w(t) can be substituted by computation of an approximate control of (Pε). In Section 2 we prove the convergence results regarding the sequence of optimal control problem (Pε). Such a convergence scheme was studied (for an optimal control problem governed by nonlinear parabolic variational inequalities) by Barbu [1]. For other works in this context see [6] and references therein. Necessary optimality con- ditions for the approximating process (Pε) (Theorem 3.1) and, a conceptual algorithm of gradient type are established in the last Section.
2. THE CONVERGENCE OF PROBLEM (Pε) The main result of this paper is ∗ ε Theorem 2.1. Let {wε} be a sequence of optimal controllers for problem (P ). Then ε lim inf L (w) = inf {L0(w); w ∈ U} (2.1) ε→0 0 and ∗ lim L0(wε) = inf {L0(w); w ∈ U}. (2.2) ε→0 A product formula approach to an inverse problem... 19
∗ Moreover, every weak limit point of {wε} is an optimal controller for problem (P). Remark 2.1. Theorem 2.1 amounts to saying that (Pε) approximates problem (P) ∗ ε and, an optimal controller {wε} of (P ) is a suboptimal controller for problem (P). The main ingredient in the proof of the Theorem 2.1 is the following Lemma.
∗ ε Lemma 2.1. If {wε} is a sequence of optimal controllers for problems (P ) then there exists {εn} −→ 0 such that
∗ ∗ ∞ wεn −→ w weakly star in L (Σ), (2.3)
∗ ∗ 2 1 uεn −→ u strongly in L ((0, T); H (Ω)), (2.4) ∗ ∗ 2 1 ϕεn −→ ϕ strongly in L ((0, T); H (Ω)), (2.5) ∗ ∗ ∗ ∗ ∗ wεn wεn ∗ where (uεn , ϕεn , wεn ) = (uεn , ϕεn , wεn ) is the solution to (1.5)-(1.8) corresponding to ∗ ∗ ∗ ∗ w∗ w∗ ∗ w = wεn and (u , ϕ , w ) = (u , ϕ , w ) is the solution to (1.1)-(1.4) corresponding to w = w∗. Proof. Details on the demonstration of this Lemma can be found in the work [6, Lemma 3.1, pp. 11]. We omit them.
We can now give the proof of Theorem 2.1.
∗ ε ∗ ∗ ∗ Proof. Let {wε} be an optimal controller for problem (P ) and let (uε, ϕε, wε) be the ∗ corresponding solution of (1.5)-(1.8) with w = wε. Lemma 2.1 above allows us to ∗ ∞ conclude that there exist w ∈ L ([0, T]) and {εn} such that relations (2.3)-(2.5) are valid. Since: Z β u −→ (u(t, x) − δ )+ 2 · χ dtdx, 2 2 Q0 ZQ 1 ϕ −→ (ϕ(t, x) − 1 − δ )2 · χ dtdx, 2 1 Q0 Q ZT 1 w −→ w2(t) dt 2 0 are convex continuous functions, it follows that these are weakly lower semicontinu- ous functions. Hence ∗ εn ∗ L0(w ) ≤ lim inf L (w ). (2.6) n−→∞ 0 εn ∗ ∗ Letw ¯ be an optimal controller for problem (P). Since wεn is an optimal controller for problem (Pεn ) it follows that
Lεn (w∗ ) ≤ L (w ¯ ε). 0 εn 0 20 Tommaso Benincasa, Costic˘aMoro¸sanu
w¯ ∗ w¯ ∗ w¯ ∗ w¯ ∗ 2 1 But (see (2.4) and (2.5)) uεn −→ u ϕεn −→ ϕ strongly in L ((0, T); H (Ω)) and so, the latter inequalities implies
εn ε ε lim L (w ¯ ) ≤ L0(w ¯ ). (2.7) n−→∞ 0 From (2.6)-(2.7) we get
L (w∗) ≤ lim inf Lεn (w∗ ) ≤ lim sup Lεn (w∗ ) ≤ L (w ¯ ∗). 0 0 εn 0 εn 0 n−→∞ n−→∞ Hence lim inf Lεn (w∗ ) = L (w ¯ ∗) = inf{L (w), w ∈ U} 0 εn 0 0 εn−→0 and then (2.1) holds. ∗ ∗ wε wε ∗ To prove (2.2) we set:u ¯ε = u ,ϕ ¯ ε = ϕ (we recall that wε is chosen to be optimal ε in (P )). On a subsequence {εn} we have
∗ 0 ∞ wε −→ w weakly star in L ([0, T]), 2 1 u¯εn −→ u strongly in L ((0, T), H (Ω)), 2 1 ϕ¯ εn −→ ϕ strongly in L ((0, T), H (Ω)), where (u, ϕ, w0) satisfy (1.1)-(1.4), i.e., (u, ϕ) = (uw0 , ϕw0 ). Therefore, we derive
0 L0(w ) ≤ inf P and, because {εn} was chosen arbitrarily, (2.2) follows. ∗ ε Now, taking into account that wε is an optimal controller for problem (P ), it fol- lows that ε ∗ ε L0(wε) ≤ L0(w) ∀w ∈ U. On the other part, on the basis of relation (2.6), we can put
∗ ε ∗ L0(w ) ≤ lim inf L (wε) ε→0 0 and thus, along with previous inequality, we may conclude that
∗ ε L0(w ) ≤ lim L (w) ∀w ∈ U. ε→0 0 Consequently ∗ L0(w ) ≤ L0(w) ∀w ∈ U i.e., the weak limit point w∗ is a suboptimal controller for problem (P). This com- pletes the proof of Theorem 2.1. A product formula approach to an inverse problem... 21 3. NECESSARY OPTIMALITY CONDITIONS IN (Pε) Let (uε, ϕε, w) be the solution of (1.5)-(1.8) and letw ˜ ∈ L∞[0, T]) be arbitrary but fixed and λ > 0. Set wλ = w + λw˜ and let (uλ,ε, ϕλ,ε) be the solution of (1.5)-(1.8) corresponding to wλ, that is: ` ρcuλ,ε + ϕλ,ε = kuλ,ε, t 2 t xx ε in Q0, (3.1) 1 τϕλ,ε = ξ2ϕλ,ε + ϕλ,ε + 2uλ,ε, t xx 2a subject to non-homogeneous Cauchy-Neumann boundary conditions:
λ,ε λ,ε λ λ,ε ε ux + hu = w , ϕx = 0 on Σ , (3.2)
λ,ε λ,ε ε ux = 0, ϕx = 0 on Σ0, (3.3) and initial conditions:
λ,ε λ,ε λ u (ε, x) = u0(x), ϕ+ (ε, x) = z (ε, ϕ0(x)) on Ωε, (3.4)
λ where z (ε, ϕ0(x)) is the solution of the Cauchy problem: 0 1 3 zλ(s) + zλ(s) = 0, s ∈ (0, ε), 2a (3.5) λ λ,ε λ,ε z (0) = ϕ˜ − (ε, x), ϕ˜ − (0, x) = ϕ0(x). Subtracting (1.5)-(1.8) from (3.1)-(3.4) and dividing by λ > 0, we get
λ,ε ε λ,ε ε λ,ε ε u − u ` ϕ − ϕ u − u ρc + = k , λ t 2 λ t λ xx in Qε, (3.6) 0 ϕλ,ε − ϕε ϕλ,ε − ϕε 1 ϕλ,ε − ϕε uλ,ε − uε τ = ξ2 + + 2 , λ t λ xx 2a λ λ uλ,ε − uε uλ,ε − uε wλ − w ϕλ,ε − ϕε + h = , = 0 on Σε, (3.7) λ x λ λ λ x uλ,ε − uε ϕλ,ε − ϕε = 0, = 0 on Σε, (3.8) λ x λ x 0
uλ,ε(ε, x) − uε(ε, x) ϕλ,ε(ε, x) − ϕε (ε, x) zλ(ε, ϕ (x)) − z(ε, ϕ (x)) = 0, + + = 0 0 on Ω . λ λ λ ε (3.9) 22 Tommaso Benincasa, Costic˘aMoro¸sanu
Letting λ tend to zero in (3.6)-(3.9) we get the system in variation (3.10)-(3.13) below ` ρcu˜ε + ϕ˜ ε = ku˜ε , t 2 t xx ε in Q0, (3.10) 1 τϕ˜ ε = ξ2ϕ˜ ε + ϕ˜ ε + 2˜uε, t xx 2a ε ε ε ε u˜ x + hu˜ = w˜ , ϕ˜ x = 0 on Σ , (3.11) ε ε ε u˜ x = 0, ϕ˜ x = 0 on Σ0, (3.12) ε ε u˜ (ε, x) = 0, ϕ˜ +(ε, x) = η(ε, x) on Ωε. (3.13) uλ,ε − uε whereu ˜ε = lim , etc., and λ−→0 λ zλ(ε, ϕ (x)) − z(ε, ϕ (x)) η(ε, x) = lim 0 0 = λ−→0 λ
λ ε = z (ε, ϕ0(x)) · ϕ˜ −(ε, x) + z˜(ε, ϕ0(x)) = z˜(ε, ϕ0(x)) with η(·, x) the solution of the Cauchy problem 3 η0(s, x) + z2(s, x)η(s, x) = 0, s ∈ (0, ε), 2a (3.14) ε η(0, x) = ϕ˜ −(ε, x), that is Zε 3 η(s, x) = exp − z(t, ·)2dt ϕ˜ ε (ε, x). (3.15) 2a − 0 We now introduce the adjoint state system. For this, the system (3.10) can be written in the abstract form: ε! ε! u˜ u˜ ε ε = A ε in Q0 ϕ˜ t ϕ˜ where (here ∆ϕ = ϕxx) ` ` 1 1 k∆− − ξ2∆+ ρc τ 2τρc 2a A = , 2 1 1 ξ2∆+ τ τ 2a ( ) ∂ψ ∂γ D(A) = (ψ, γ) ∈ H2(Ω) × H2(Ω); + hψ ∈ L2(∂Ω), = 0 . ∂ν ∂ν A product formula approach to an inverse problem... 23
Then ` 2 1 k∆− ρc τ τ ∗ A = , ` 1 1 1 − ξ2∆+ ξ2∆+ 2τρc 2a τ 2a
( ) ∂ψ ∂γ 2ρc ∂ψ D(A∗) = (ψ, γ) ∈ H2(Ω) × H2(Ω); + hψ = 0, = . ∂ν ∂ν ` ∂ν Thus, the adjoint state system is
ε! ε! ∂ ε ! p ∗ p ∂uε L0(w) ε = −A ε + ∂ ε q t q ∂ϕε L0(w) i.e. k ` 2 pε + pε − pε + qε = β(uε − d )+ · χ , Qε, t xx 2 Q0 0 ρc τρc τ (3.16) pε + hpε = 0, on Σε, x ε ε p−(ε, x) = 0, p−(T, x) = 0 x ∈ ΩT , 2 2 ε `ξ ε ` ε ξ ε 1 ε ε ε qt − pxx − p + qxx + q =(ϕ −1−d1) · χQ , in Q , 2τρc 4aτρc τ 2aτ 0 0 ` qε = pε, on Σε, x 2ρc x (3.17) ε ε qx = 0, on Σ0, ε Z ε 3 2 ε ε q−(ε, x) = exp z (t, ·)dt q+(ε, x), q−(T, x) = 0, x ∈ ΩT . 2a 0 Let us introduce the cost functional 1 Lε(w) = Lε(w) + I (w) 1 0 2 U where, as usually, IU(w) is the indicator function of the set U. If w∗ is an optimal controller of problem (Pε), then Lε(w∗ + λw˜ ) − Lε(w∗) 1 1 ≥ 0 ∀λ > 0. λ that leads to (letting λ tend to zero) 24 Tommaso Benincasa, Costic˘aMoro¸sanu R Z ε ε + ε ε β u˜ (u − d2) · χQ0 dtdx + ϕ˜ (ϕ − 1 − d1) · χQ0 dtdx+ (3.18) Q Q ZT ∗ 0 ∗ ∗ + w w˜ dt + IU(w , w˜ ) ≥ 0 ∀w˜ ∈ TU(w ). 0 ε ε Multiplying (3.16)1 byu ˜ and (3.17)1 byϕ ˜ , using integration by parts and Green’s formula, we get R R R ε ε k ε ε ` ε ε pt u˜ dt dx + ρc p u˜ xxdt dx − τρc p u˜ dt dx+ Qε Qε Qε 0 R 0 R 0 + 2 qεu˜εdt dx + k pεu˜ε − pεu˜ε dt dγ = τ ρc x x (3.19) Qε Σε R0 ε + ε = β (u − d2) · χQ0 u˜ dtdx, ε Q0
R R R ε ε `ξ2 ε ε ` ε ε qt ϕ˜ dt dx − 2τρc p ϕ˜ xxdt dx − 4aτρc p ϕ˜ dt dx+ Qε Qε Qε 0 R 0 R 0 `ξ2 ε ε ε ε ξ2 ε ε ε ε + 2τρc pxϕ˜ − p ϕ˜ x dt dγ + τ q ϕ˜ x − ϕ˜ qx dtdγ+ (3.20) Σε Σε 2 R R R + ξ qεϕ˜ ε dt dx + 1 qεϕ˜ εdt dx = (ϕε − 1 − d ) · χ ϕ˜ ε dtdx. τ xx 2aτ 1 Q0 ε ε ε Q0 Q0 Q0 ε ε Now we multiply (3.11) by p , (3.16)2 byu ˜ , by subtraction we get
ε ε ε ε ε pxu˜ − u˜ x p = −p w˜ . (3.21)
Adding (3.19)-(3.20) and taking into account (3.11)2, (3.17)2, (3.21), we obtain R R R ε ε ε ε k ε pt u˜ dt dx + qt ϕ˜ dt dx + ρc p wdt˜ dγ + ε ε ε Q0 Q0 Σ R h 2 i ε k ε `ξ ε ` ε ` ε + p u˜ xx − ϕ˜ xx − ϕ˜ − u˜ dt dx + Qε ρc 2τρc 4aτρc τρc R0 h i ε ξ2 ε 1 ε 2 ε + q τ ϕ˜ xx + 2aτ ϕ˜ + τ u˜ dt dx = Qε 0 R R ε + ε ε ε = β (u − d2) · χQ0 u˜ dtdx + (ϕ − 1 − d1) · χQ0 ϕ˜ dtdx, ε ε Q0 Q0 i.e., making use of equations in (3.10), the last relation leads to R R ε ε ε ε ε ε ε ε k ε pt u˜ + p u˜t + qt ϕ˜ + q ϕ˜ t dt dx + ρc p wdt˜ dγ = Qε Σε 0 R R ε + ε ε ε = β (u − d2) · χQ0 u˜ dtdx + (ϕ − 1 − d1) · χQ0 ϕ˜ dtdx, ε ε Q0 Q0 A product formula approach to an inverse problem... 25
By Fubini’s theorem and definition of distributional derivative, the latter relation give us Z Z Z k pεwdt˜ dγ = β (uε − d )+ · χ u˜ε dtdx + (ϕε − 1 − d ) · χ ϕ˜ ε dtdx, ρc 2 Q0 1 Q0 ε ε ε Σ Q0 Q0 and then (3.18) becomes
Z ZT k pεw˜ dtdγ + w∗w˜ dt + I0 (w∗, w˜ ) ≥ 0 ∀w˜ ∈ T (w∗) ρc U U Σε 0 or ZT k pε(s, b ) + w∗(s) w˜ (s) ds + I0 (w∗, w˜ ) ≥ 0 ∀w˜ ∈ T (w∗). ρc 1 U U 0 The last inequality is equivalent with
∗ −r(t) ∈ ∂IU(w ) a.p.t. (t, x) ∈ [0, T],
k ε where r(t) = ρc p (t, b1) + w(t), and thus we can conclude that ( 0, if r(t) > 0, w∗(t) = (3.22) −R, i f r(t) < 0.
Summing up, we have proved the following maximum principle for problem (Pε) Theorem 3.1. Let (u∗,ε, ϕ∗,ε, w∗) be optimal in problem (Pε). Then the optimal con- trol is given by (3.22) where (pε, qε) satisfy along with u∗,ε, ϕ∗,ε the dual system (3.16) − (3.17). Now we will present a numerical algorithm of gradient type in order to compute the approximating optimal control stated by Theorem 3.1. Algorithm InvPHT1D (Inverse PHase Transition case 1D)
P0. Choose w(0) ∈ U and set iter= 0; Choose ε > 0; P1. Compute z(ε, ·) from (1.9); P2. Compute (uε,iter, ϕε,iter) from (1.5)-(1.8); P3. Compute (pε,iter, qε,iter) from (3.16)-(3.17); P4. For t ∈ [0, T], compute k riter(t) = · pε,iter(t, b ) + witer; ρc 1 26 Tommaso Benincasa, Costic˘aMoro¸sanu
P5. Set ( 0, if riter(t) > 0, w˜ iter(t) = −R, if riter(t) < 0.
P6. Compute λiter ∈ [0, 1] solution of the minimization process:
ε iter iter min {L0(λw + (1 − λ)w ˜ , λ ∈ [0, 1]; iter+1 iter iter Set w = λiterw + (1 − λiter)w ˜ ; P7. If k witer+1 − witer k ≤ η /* the stopping criterion */ then STOP else iter:= iter+1; Go to P1. In the above, the variable iter represents the number of iterations after which the ε ε algorithm InvPHT1D found the optimal value of the cost functional L0(w) in (P ). ε iter+1 ε iter Remark 3.1. The stopping criterion in P7 could be k L0(w ) − L0(w ) k ≤ η, where η is a prescribed precision.
4. CONCLUSIONS The main novelty brought by this work is that the computation of the approximate solution corresponding to the nonlinear system (1.1) is replaced with calculation of the approximate solution for an ordinary equation and a linear system (compare step P1 in [5] with the steps P1-P2 in present paper). Numerical implementation of the conceptual algorithm InvPHT1D remain an open problem. We wish only to draw attention to the type of boundary condition considered here (see (1.2)) namely that they fully cover industrial application pro- posed by us for numerical simulations - a matter for further investigation.
References [1] V. Barbu, A product formula approach to nonlinear optimal control problems, SIAM J. Control Optim., 26(1988), 496-520. [2] T. Benincasa, C. Moros¸anu, Fractional steps scheme to approximate the phase-field transition system with nonhomogeneous Cauchy-Neumann boundary conditions, Numer. Funct. Anal. and Optimiz., 30, 3-4(2009), 199-213. [3] G. Caginalp, An analysis of a phase field model of a free boundary, Arch. Rat. Mech. Anal., 92(1986), 205-245. [4] 4. M. Heinkenschloss, F. Troltzsch,¨ Analysis of the Lagrange-SQP- Newton Method for the Con- trol of a Phase Field Equation, Control & Cybernetics, 28, 2(1999), 177–211. [5] C. Moros¸anu, Numerical approach of an inverse problem in the phase field equations, An. S¸t. Univ. ”Al.I. Cuza” Ias¸i, T XXXIX, s. I-a, 4(1993), 419-436. [6] C. Moros¸anu, Boundary optimal control problem for the phase-field transition system using frac- tional steps method, Control & Cybernetics, 32, 1(2003), 05-32. ROMAI J., 6, 2(2010), 27–39
LIMIT CYCLES BY FINITE ELEMENT METHOD FOR A ONE - PARAMETER DYNAMICAL SYSTEM ASSOCIATED TO THE LUO - RUDY I MODEL C˘at˘alinLiviu Bichir1, Adelina Georgescu2, Bogdan Amuzescu 3, Gheorghe Nistor 4, Marin Popescu 5, Maria-Luiza Flonta 6, Alexandru Dan Corlan 7, Istvan Svab 8 1 Rostirea Maths Research, Galat¸i, Romania, 2 Academy of Romanian Scientists, Bucharest, Romania, 3, 6, 8 Faculty of Biology, University of Bucharest, Romania, 4, 5 University of Pite¸sti,Romania, 7 Bucharest University Emergency Hospital, Romania, [email protected], [email protected], [email protected], [email protected], fl[email protected], [email protected], [email protected] Abstract A one - parameter dynamical system is associated to the mathematical problem gov- erning the membrane excitability of a ventricular cardiomyocyte, according to the Luo- Rudy I model. Limit cycles are described by the solutions of an extended system. A finite element method time approximation (FEM) is used in order to formulate the ap- proximate problem. Starting from a Hopf bifurcation point, approximate limit cycles are obtained, step by step, using an arc-length-continuation method and Newton’s method. Some numerical results are presented.
Keywords: limit cycle, finite element method time approximation, Luo-Rudy I model, arc-length- continuation method, Newton’s method. 2000 MSC: 37N25 37G15 37M20 65L60 90C53 37J25.
1. INTRODUCTION The well-known Hodgkin-Huxley model of the squid giant axon [16] represented a huge leap forward in comparison with the earlier models of excitable systems built from abstract sets of equations or from electrical circuits including non-linear com- ponents, e.g. [33]. The pioneering work of Denis Noble’s group made the transition from neuronal excitability models, characterized by Na+ and K+ conductances with fast gating kinetics, to cardiomyocyte electrophysiology models, a field expanding steadily for over five decades [23]. Nowadays, complex models accurately reproduc- ing transmembrane voltage changes as well as ion concentration dynamics between various subcellular compartments and buffering systems are incorporated into de- tailed anatomical models of the entire heart [24]. The Luo-Rudy I model of isolated guinea pig ventricular cardiomyocyte [21] was developed in the early 1990s starting
27 28 C. L. Bichir, A. Georgescu, B. Amuzescu, Gh. Nistor, M. Popescu, et al from the Beeler-Reuter model [1]. It includes more recent experimental data related to gating and permeation properties of several types of ion channels, obtained in the late 1980s with the advent of the patch-clamp technique [22]. The model com- prises only three time and voltage-dependent ion currents (fast sodium current, slow inward current, time-dependent potassium current) plus three background currents (time-independent and plateau potassium current, background current), their dynam- ics being described by Hodgkin-Huxley type equations. This apparent simplicity, compared to more recent multicompartment models, renders it adequate for mathe- matical analysis using methods of linear stability and bifurcation theory. Nowadays, there exist numerous software packages for the numerical study of finite - dimensional dynamical systems, for example MATCONT, CL−MATCONT, CL−MATCONTM [7], [15], AUTO [8]. In [19], [8], [7], [15]; the periodic boundary value problems used to locate limit cycles are approximated using orthogonal col- location method. Finite differences method is also considered. In this paper, limit cycles are obtained for the dynamical system associated to the Luo-Rudy I model by using finite element method time approximation (FEM).
2. LUO-RUDY I MODEL The mathematical problem governing the membrane excitability of a ventricular cardiomyocyte, according to the Luo-Rudy I model [21], is a Cauchy problem
u(0) = u0 , (1) for the system of first order ordinary differential equations du = F(η, u) , (2) dt where u = (u1,..., u8) = (V,[Ca]i, h, j, m, d, f , X), η = (η1, . . . , η13) = (Ist, Cm, gNa, 8 13 gsi, gKp, gb,[Na]0,[Na]i,[K]0,[K]i, PRNaK, Eb, T), M = R , F : R × M → M, F = (F1,..., F8),
1 3 F1(η, u) = − [Ist + η3u3u4u5(u1 − ENa(η7, η8, η13)) η2 +η4u6u7(u1 − c1 + c2 ln u2) +gK(η10)Xi(u1)(u1 − EK(η7, η8, η9, η10, η11, η13))u8
+gK1(η10)K1∞(η9, η10, η13, u1)(u1 − EK1(η9, η10, η13)) +η5Kp(u1)(u1 − EK p(η9, η10, η13)) + η6(u1 − η12)] , F2(η, u) = −c3η4u6u7(u1 − c1 + c2 ln u2) + c4(c5 − u2) ,
F`(η, u) = α`(u1) − (α`(u1) + β`(u1))u` , ` = 3,..., 8 .
For the definition of variables V,[Ca]i, h, j, m, d, f , X, parameters Ist, Cm, gNa, gsi, gK p, gb,[Na]0,[Na]i,[K]0,[K]i, PRNaK, Eb, T, constants c1,..., c5, functions Limit cycles by Finite Element Method for a one - parameter dynamical system... 29 gK, ENa, EK, EK1, EKp, K1∞, Xi, K p, α`, β`, default values of parameters and initial values of variables in the Luo-Rudy I model, the reader is referred to [21]. The reader is also referred to [20] for the continuity of the model, and to [4] for the treatment of the vector field F singularities. F is of class C2 on the domain of interest.
3. THE ONE - PARAMETER DYNAMICAL SYSTEM ASSOCIATED TO THE LUO - RUDY I MODEL We performed the study of the dynamical system associated with the Cauchy prob- lem (1), (2) by considering only the parameter η1 = Ist and fixing the rest of param- eters. Denote λ = η1 = Ist and η∗ the vector of the fixed values of η2, . . . , η13. Let F : R × M → M, F(λ, u) = F(λ, η∗, u), F = (F1,..., F8). Consider the dynamical system associated with the Cauchy problem (1), (3), where du = F(λ, u) . (3) dt The equilibrium points of this problem are solutions of the equation
F(λ, u) = 0 . (4)
The existence of the solutions and the number were established by graphical repre- sentation in [4], for the domain of interest. The equilibrium curve (the bifurcation di- agram) was obtained in [4], via an arc-length-continuation method [13] and Newton’s method [12], starting from a solution obtained by solving a nonlinear least-squares problem [13] for a value of λ for which the system has one solution. In [4], the re- sults are obtained by reducing (4) to a system of two equations in (u1, u2) = (V, [Ca]i). Here, we used directly (4).
4. EXTENDED SYSTEM METHOD FOR LIMIT CYCLES The extended system in (λ, T, u) du − TF(λ, u) = 0, dτ u(0) − u(1) = 0, (5) R1 dw(t) < u(t), > dt = 0, 0 dt was introduced, in [19], [8], [7], in order to locate limit cycles of a general problem (1), (3). Here, T is the unknown period of the cycle, w is a component of a known 30 C. L. Bichir, A. Georgescu, B. Amuzescu, Gh. Nistor, M. Popescu, et al reference solution (λ,ˆ Tˆ, w) of (5), and, in our case, tuv X8 X8 2 8 < u, v >= uivi and kuk = ui for u, v ∈ R . i=1 i=1
The system (5) becomes determined in a continuation process. In order to approx- imate and solve it by FEM time approximation, let us obtain the weak form of (5) in the sequel. For this, we consider the function spaces: dx X = {x ∈ L2(0, 1; R8); ∈ L2(0, 1; R8), dt x = (x1,..., x8), xi(0) = xi(1), i = 1,..., 8} . dv V = {v ∈ L2(0, 1; R); ∈ L2(0, 1; R), v(0) = v(1)} . dt
The weak form of (5) is the problem in (λ, T, u) ∈ R × R × X Z1 Z1 dv(τ) u (τ) dτ + T F (λ, u(τ))v(τ)dτ = 0, i dτ i 0 0 ∀v ∈ V, i = 1,..., 8, (6) Z1 dw(t) < u(t), > dt = 0 . dt 0 5. ARC-LENGTH-CONTINUATION METHOD FOR (6) Following the usual practice ([17], [18], [7], [8], [12], [13], [14], [15], [19], [25], [27], [28], [29]), we also use an arc-length-continuation method in order to formulate an algorithm to solve (6) approximatively. Glowinski ([13], following Keller [17], [18]) and Doedel ([8], where Keller’s name is also cited) chose a continuation equation written in our case as
Z1 du(t) dT dλ k k2dt + ( )2 + ( )2 = 1 , (7) ds ds ds 0 where s is the curvilinear abscissa. Let (λ0, u0) be a Hopf bifurcation point, ±β0i a pair of purely imaginary eigenval- 0 0 0 0 0 ues of of the Jacobian matrix DuF(λ , u ), and a nonzero complex vector g = gr +igi . (λ0, u0) is located on the equilibrium curve during a continuation procedure using Limit cycles by Finite Element Method for a one - parameter dynamical system... 31
0 0 0 0 0 8 8 8 some test functions ([19], [14], [7]). (λ , β , u , gr , gi ) ∈ R × R × R × R × R is the solution of the extended system ([27], [28], [29])
F(λ, u) D F(λ, u)g + βg u r i D F(λ, u)g − βg u i r = 0 , (8) g − 1 r,k gi,k where k is a fixed index of gr and of gi, 1 ≤ k ≤ 8. To solve (6), the extended system formed by (6) and (7), parametrized by s, was considered. Let 4s be an arc-length step and λn u(λ4s), T n T(n4s), un u(n4s). We have the algorithm (following the cases from [13], [8], [28], [29]): Algorithm 1. 1. take the Hopf bifurcation point (λ0, u0) and T 0 = 2π/β0; 0 0 retain gr , gi ; 2. for n = 0, (λ1, T 1, u1) ∈ R × R × X is obtained ([8], [29]) by (13) below
1 Z X8 dφ (t) u1(t) i dt = 0 , (9) i dt 0 i=1 and Z1 8 X 1 0 (ui (t) − ui (t))φi(t) dt = 4s , (10) 0 i=1 where 0 0 φ(t) = sin(2πt)gr + cos(2πt)gi , (11) by using Newton’s method with the initial iteration
(u1)0(t) = u0 + 4s φ(t) , (T 1)0 = T 0 , (λ1)0 = λ0 . (12)
3. for n ≥ 1, assuming that (λn−1,T n−1, un−1), (λn,T n, un) are known, (λn+1, T n+1, un+1) ∈ R × R × X is obtained from
1 1 Z dv(τ) Z un+1(τ) dτ + T n+1 F (λn+1, un+1(τ))v(τ)dτ = 0 , (13) i dτ i 0 0 ∀v ∈ V, i = 1,..., 8,
1 Z X8 dun(t) un+1(t) i dt = 0, (14) i dt 0 i=1 32 C. L. Bichir, A. Georgescu, B. Amuzescu, Gh. Nistor, M. Popescu, et al
Z1 X8 un(t) − un−1(t) (un+1(t) − un(t)) i i dt+ i i 4s 0 i=1 T n − T n−1 λn − λn−1 +(T n+1 − T n) + (λn+1 − λn) = 4s,, (15) 4s 4s by using Newton’s method with the initial iteration
((λn+1)0, (T n+1)0, (un+1)0) = (λn, T n, un). (16)
6. NEWTON’S METHOD FOR THE STEPS OF ALGORITHM 1 In (15) (n ≥ 1), let us denote λ∗ = λn, T ∗ = T n, u∗ = un, λ∗∗ = (λn − λn−1)/4s, T ∗∗ = (T n − T n−1)/4s, u∗∗ = (un − un−1)/4s. We write (13), (9), (10) (the iteration n = 0) in the same general form as (13), (14), (15). So denote u∗ = u0, u∗∗ = φ and consider λ∗ = λ0, T ∗ = T 0, λ∗∗ = 0, T ∗∗ = 0 in (15) and consider u∗ = u0 = φ in (14). Each step of Algorithm 1, given (λ∗, T ∗, u∗), (λ∗∗, T ∗∗, u∗∗), calculates (λn+1, T n+1, un+1) ∈ R × R × X, n ≥ 0, by (13),
Z1 X8 du∗(t) un+1(t) i dt = 0 , (17) i dt 0 i=1 and
Z1 X8 n+1 ∗ ∗∗ n+1 ∗ ∗∗ n+1 ∗ ∗∗ (ui (t) − ui (t))ui (t) dt + +(T − T )T + (λ − λ )λ = 4s . (18) 0 i=1 Newton’s method applied (13), (17) and (18), for n ≥ 0, leads to:
let ((λ1)0,(T 1)0,(u1)0), given by (12), be an initial iteration (m = 0), if n = 0;
let ((λn+1)0,(T n+1)0,(un+1)0), given by (16), be an initial iteration (m = 0), if n ≥ 1;
calculate (λn+1, T n+1, un+1) as the solution of the algorithm: for m ≥ 0, ((λn+1)m+1,(T n+1)m+1,(un+1)m+1) = (λm+1, T m+1, um+1) ∈ R × R × X is obtained by Limit cycles by Finite Element Method for a one - parameter dynamical system... 33
Z1 Z1 dv(τ) um+1(τ) dτ + T m+1 F (λm, um(τ))v(τ)dτ+ i dτ i 0 0 Z1 m m m m+1 m+1 +T DFi(λ , u (τ))(λ , u (τ))v(τ)dτ = (19) 0 Z1 m m m m m = T DFi(λ , u (τ))(λ , u (τ))v(τ)dτ , ∀v ∈ V, i = 1,..., 8. 0 Z1 X8 du∗(t) um+1(t) i dt = 0, (20) i dt 0 i=1
Z1 X8 m+1 ∗∗ m+1 ∗∗ m+1 ∗∗ ui (t)ui (t) dt + T T + λ λ = (21) 0 i=1 Z1 X8 ∗ ∗∗ ∗ ∗∗ ∗ ∗∗ = ui (t)ui (t) dt + T T + λ λ + 4s . 0 i=1 7. APPROXIMATION OF PROBLEM (19), (20), (21) BY FINITE ELEMENT METHOD TIME APPROXIMATION In order to perform this approximation, let us divide the interval [0, 1] in N + 1 subintervals K = K j = [t j, t j+1], 0 ≤ j ≤ N, where 0 = t0 < t1 < . . . < tN+1 = 1. The sets K represent a triangulation Th of [0, 1]. Let us approximate the spaces V and X by the spaces
Vh = {v : [0, 1] → R; v ∈ C[0, 1], v(0) = v(1), v|K ∈ Pk(K), ∀K ∈ Th},
8 Xh = {x : [0, 1] → R ; x = (x1,..., x8), xi ∈ Vh, i = 1,..., 8} , respectively, where Pk(K) is the space of polynomials in t of degree less than or equal to k defined on K, k ≥ 2. Let k = 2. An element K ∈ Th has three nodal points. To obtain a function uh ∈ Xh reduces to obtain a function vh ∈ Vh. In order to obtain a function vh ∈ Vh, we use a basis of functions of Vh. Let JK = {1, 2, 3} be the local numeration for the nodes of K, where 1, 3 correspond to t j, t j+1 respectively and 2 corresponds to a node between t j and t j+1. Let {ψi, i ∈ JK} be the local quadratic basis of functions on K corresponding to the local nodes. Let J = {1,..., 2N + 1} be the global numeration 34 C. L. Bichir, A. Georgescu, B. Amuzescu, Gh. Nistor, M. Popescu, et al for the nodes of [0, 1]. The two numerations are related by a matrix L whose elements are the elements j ∈ J. Its rows are indexed by the elements K ∈ Th (by the number of the element K in a certain fixed numeration with elements from the set {1,..., N}) and its columns, by the local numeration i ∈ JK, that is j = L(K, i). A function vh ∈ Vh is defined by its values v j from the nodes j ∈ J, X X vh(t) = v j ψi(t) , (22) K∈ Th i∈ JK , j=L(K,i) and a function uh ∈ Xh is defined by its values u j from the nodes j ∈ J, X X uh(t) = u j ψi(t) . (23) K∈ Th i∈ JK , j=L(K,i)
So, an unknown function uh = ((uh)1, ...,(uh)8) is reduced to the unknowns u j, u j = ((u j)1, ...,(u j)8), j ∈ J. In (19), (20), (21), approximate (λm+1, T m+1, um+1) ∈ R × R × X by (λm+1, T m+1, m+1 m+1 uh ) ∈ R × R × Xh. Taking uh = uh, uh given by (23), and v = ψ`, for all ` ∈ JK, for all K ∈ Th, we obtain the discrete variant of problem (19), (20), (21) as the 8·(2N+1) following problem in (λ, T, u1, ..., u2N+1) ∈ R × R × R , written suitable for the assembly process, X X Z Z m+1 ψ` (τ) m+1 m m (u )n ψ (τ) dτ + T Fi(λ , u (τ)) ψ dτ + j i dτ h ` K∈ Th i∈ JK , j=L(K,i) K K X Z m m+1 m m + T < u j , (DuFn(λ , uh (τ))ψi(τ)) ψ` dτ > + i∈ JK , j=L(K,i) K Z m+1 m m +λ (DλFn(λ , u (τ))) ψ dτ = (24) h ` K X Z Z m m m m m m m = T (DuFn(λ , u (τ))u (τ)) ψ dτ + (DλFn(λ , u (τ))λ ) ψ dτ , h h ` h ` K∈ Th K K n = 1,..., 8, m+1 m+1 u0 = u2N+1 , (25) Z X X du∗(τ) < um+1, ψ (τ) h dτ > = 0 , (26) j i dτ K∈ Th i∈ JK , j=L(K,i) K Limit cycles by Finite Element Method for a one - parameter dynamical system... 35
X X Z m+1 ∗∗ m+1 ∗∗ m+1 ∗∗ < u j , ψi(τ) uh (τ) dτ > +T T + λ λ = (27) K∈ Th i∈ JK , j=L(K,i) K
X Z ∗ ∗∗ ∗ ∗∗ ∗ ∗∗ = < uh(τ), uh (τ) > dτ + T T + λ λ + 4s , K∈ Th K
for all ` ∈ JK, for all K ∈ Th. 8. NUMERICAL RESULTS Based on [30] and on the computer programs for [2] and [3], relations (24), (25), (26), (27) and the algorithm at the end of section 5 furnished the numerical results of this section. Let (λ0, u0) be the Hopf bifurcation point located during the construction of the equilibrium curve by a continuation procedure in [5]. 0 0 0 0 0 The solution (λ , β , u , gr , gi ) of (8), calculated in [5], is λ0 = −1.0140472901 , β0 = 0.0162886062, u0 = (−24.3132508542, 0.0034641214, 0.0, 0.0, 0.9176777444, 0.5025242162, 4920204612, 0.5071561613), 0 gr = (1.0, 0.0000468233, 0.0, 0.0, 0.0093195354, 0.0198748652, −0.0072420216, 0.0001706577), 0 gi = (0.0, 0.0000029062, 0.0, 0.0, −0.0000171311, −0.0118192370, 0.0136415789, −0.0017802907). 0 0 The eigenvalues of the Jacobian matrix DuF(λ , u ), calculated by the QR al- gorithm, are ± 0.0162886062 i, −8.8611865338, −0.1026761869, −0.0647560667, −0.0024565181, −1.7398266947, −0.2049715178. These data are considered in the step 1 of the algorithm at the end of section 5. We took Cm = 1, gNa = 23, gsi = 0.09, gK = 0.282, gK1 = 0.6047, gK p = 0.0183, Gb = 0.03921, [Na]0 = 140, [Na]i = 18, [K]0 = 5.4, [K]i = 145, PRNaK = 0.01833, Eb = −59.87, T = 310. In order to solve (6) numerically by the algorithm at the end of section 5 and by (24), (25), (26), (27), we performedR calculations using 4s = 1.0 and 500 iterations in the continuation process. Integrals f (τ) dτ were calculated using Gauss integration K formula with three integration points. Figure 1 and 2 present some results obtained using 20 elements K (41 nodes) (N = 20, J = {1,..., 41} in section 7). The curves of the projections of the limit cycles, on the planes indicates in figure, are plots generated from values calculated in the nodes, corresponding to a fixed value of the parameter. Two projections of some limit cycles and of a part of the equilibrium curve (marked by ”N”) are presented in Fig. 1. The Hopf bifurcation point is marked by ”•”. 36 C. L. Bichir, A. Georgescu, B. Amuzescu, Gh. Nistor, M. Popescu, et al
1
0.9
0.8
0.7 f 0.6
0.5
0.4
0 −20 −40 0 −60 −1 −0.5 −80 −2 −1.5 V I st
0.6
0.5
0.4
X 0.3
0.2
0.1
0 −20 −40 −60 0 −80 −2 −1.5 −1 −0.5 V I st
Fig. 1. Two projections of limit cycles and of a part of the equilibrium curve (marked by ”N”). The Hopf bifurcation point is marked by ”•”. Limit cycles by Finite Element Method for a one - parameter dynamical system... 37
−3 −3 x 10 x 10 5 5 i
4 h 0 [Ca] 3 −5 −30 −20 −10 0 −30 −20 −10 0 V V −3 x 10 1 5 0.9 m j 0 0.8
0.7 −5 −30 −20 −10 0 −30 −20 −10 0 V V
0.8 0.8 f d 0.6 0.6 0.4 0.4 0.2 −30 −20 −10 0 −30 −20 −10 0 V V
0.55
X 0.5
0.45 −30 −20 −10 0 V
Fig. 2. Projections of two limit cycles calculated for Ist = −1.2000465026 and for Ist = −1.2000183729 (marked by ”x”) (20 elements, 41 nodes). 38 C. L. Bichir, A. Georgescu, B. Amuzescu, Gh. Nistor, M. Popescu, et al
In Fig. 2, there are represented the projections of the plots of two limit cycles calculated for Ist = −1.2000465026 (iteration 148) and for Ist = −1.2000183729 (iteration 248, marked by ”x” in figure). The results obtained are relevant from a biological point of view, pointing to un- stable electrical behavior of the modeled system in certain conditions, translated into oscillatory regimes such as early afterdepolarizations [32] or self-sustained oscilla- tions [4], which may in turn synchronize, resulting in life-threatening arrhythmias: premature ventricular complexes or torsades-de-pointes, degenerating in rapid poly- morphic ventricular tacycardia or fibrillation [26]. Acknowlegdements: This research was partially supported from grant PNCDI2 61- 010 to M-LF by the Romanian Ministry of Education, Research, and Innovation.
References
[1] G. W. Beeler , H. Reuter, Reconstruction of the action potential of ventricular myocardial fibres, J. Physiol. 268(1977), 177-210. [2] C. L. Bichir, A.Georgescu, Approximation of pressure perturbations by FEM, Scientific Bulletin of the Pites¸ti University, the Mathematics-Informatics Series, 9 (2003), 31-36. [3] C. L. Bichir, A numerical study by FEM and FVM of a problem which presents a simple limit point, ROMAI J., 4, 2(2008), 45-56, http://www.romai.ro, http://rj.romai.ro. [4] C. L. Bichir, B. Amuzescu, A. Georgescu, M. Popescu, Ghe. Nistor, I. Svab, M. L. Flonta, A. D. Corlan, Stability and self-sustained oscillations in a ventricular cardiomyocyte model, submitted to Interdisciplinary Sciences - Computational Life Sciences, Springer. [5] C. L. Bichir, A. Georgescu, B. Amuzescu, Ghe. Nistor, M. Popescu, M. L. Flonta, A. D. Corlan, I. Svab, Limit points and Hopf bifurcation points for a one - parameter dynamical system associated to the Luo - Rudy I model, to be published. [6] C. Cuvelier, A.Segal, A.A.van Steenhoven, Finite Element Methods and Navier-Stokes Equa- tions, Reidel, Amsterdam, 1986. [7] A. Dhooge, W. Govaerts, Yu.A. Kuznetsov, W. Mestrom, A.M. Riet, B. Sautois, MATCONT and CL−MATCONT: Continuation toolboxes in MATLAB, 2006, http://www.matcont.ugent.be/manual.pdf [8] E. Doedel, Lecture Notes on Numerical Analysis of Nonlinear Equations, 2007, http://cmvl.cs.concordia.ca/publications/notes.ps.gz, from the Home Page of the AUTO Web Site, http://indy.cs.concordia.ca/auto/. [9] A.Georgescu, M.Moroianu, I.Oprea, Bifurcation Theory. Principles and Applications, Applied and Industrial Mathematics Series, 1, University of Pites¸ti, 1999. [10] W. J. Gibb , M. B. Wagner, M. D. Lesh, Effects of simulated potassium blockade on the dynamics of triggered cardiac activity, J. theor. Biol 168(1994), 245-257. [11] V.Girault, P.-A.Raviart, Finite Element Approximation of the Navier-Stokes Equations, Springer, Berlin, 1979. [12] V.Girault, P.-A.Raviart, Finite Element Methods for Navier-Stokes Equations.Theory and Algo- rithms, Springer, Berlin, 1986. [13] R.Glowinski, Numerical Methods for Nonlinear Variational Problems, Springer, New York, 1984. Limit cycles by Finite Element Method for a one - parameter dynamical system... 39
[14] W.J.F. Govaerts, Numerical methods for Bifurcations of Dynamical Equilibria, SIAM, Philadel- phia, 2000. [15] W. Govaerts, Yu. A. Kuznetsov, R. Khoshsiar Ghaziani, H.G.E. Meijer, Cl MatContM: A toolbox for continuation and bifurcation of cycles of maps, 2008, http://www.matcont.ugent.be/doc−cl−matcontM.pdf [16] A. L. Hodgkin, A. F. Huxley, A quantitative description of membrane current and its application to conduction and excitation in nerve, J. Physiol., 117 (1952), 500-544. [17] H.B.Keller, Numerical Solution of Bifurcation Eigenvalue Problems, in Applications in Bifurca- tion Theory, ed. by P.Rabinowitz, Academic, New York, 1977. [18] H.B.Keller, Global Homotopies and Newton Methods, in Recent Advances in Numerical Meth- ods, ed. by C. de Boor, G.H.Golub, Academic, New York, 1978. [19] Yu. A. Kuznetsov, Elements of Applied Bifurcation Theory, Springer, New York, 1998. [20] L. Livshitz, Y. Rudy, Uniqueness and stability of action potential models during rest, pacing, and conduction using problem - solving environment, Biophysical J., 97 (2009), 1265-1276. [21] C.H. Luo, Y. Rudy, A model of the ventricular cardiac action potential. Depolarization, repolar- ization, and their interaction, Circ. Res., 68 (1991), 1501-1526. [22] E. Neher, B. Sakmann, Single-channel currents recorded from membrane of denervated frog muscle fibres, Nature, 260 (1976), 799-802. [23] D. Noble, Modelling the heart: insights, failures and progress, Bioessays, 24 (2002), 1155-1163. [24] D. Noble, From the Hodgkin-Huxley axon to the virtual hear, J. Physiol., 580 (2007), 15-22. Epub 2006 Oct 2005. [25] T.S.Parker, L.O.Chua, Practical Numerical Algorithms for Chaotic Systems, Springer, New York, 1989. [26] D. Sato, L. H. Xie, A. A. Sovari, D. X. Tran, N. Morita, F. Xie, H. Karagueuzian, A. Garfinkel, J. N. Weiss, Z. Qu , Synchronization of chaotic early afterdepolarizations in the genesis of cardiac arrhythmias, Proc. Natl. Acad. Sci. USA, 106 (2009), 2983-2988. Epub 2009 Feb 2913. [27] R. Seydel, Numerical computation of branch points in nonlinear equations, Numer. Math., 33 (1979), 339-352. [28] R. Seydel, Nonlinear Computation, invited lecture and paper presented at the Distinguished Ple- nary Lecture session on Nonlinear Science in the 21st Century, 4th IEEE International Workshops on Cellular Neural Networks and Applications, and Nonlinear Dynamics of Electronic Systems, Sevilla, June, 26, 1996. [29] R. Seydel, Practical Bifurcation and Stability Analysis, Springer, New York, 2010. [30] C. Taylor, T.G. Hughes, Finite Element Programming of the Navier-Stokes Equations, Pineridge Press, Swansea, U.K., 1981. [31] R.Temam, Navier-Stokes equations. Theory and numerical analysis, North-Holland, Amsterdam, 1979. [32] D. X, Tran, D. Sato, A. Yochelis, J. N. Weiss, A. Garfinkel, Z. Qu, Bifurcation and chaos in a model of cardiac early afterdepolarizations, Phys. Rev. Lett., 102:258103 (2009). Epub 252009 Jun 258125. [33] B. Van der Pol, J. Van der Mark, The heartbeat considered as a relaxation oscillation and an electrical model of the heart, Phil. Mag. (suppl.), 6 (1928), 763-775.
ROMAI J., 6, 2(2010), 41–53
THE FACTORIZATION OF THE RIGHT PRODUCT OF TWO SUBCATEGORIES Dumitru Botnaru, Alina T¸urcanu Technical University of Moldova, Chi¸sin˘au,Republic of Moldova [email protected], alina−[email protected] Abstract In the category of locally convex spaces the right product of two subcategories is a reflective subcategory. In the topological completely regular spaces a similar property is not always true. The factorization of this product according to a structure of factorization leads always to a reflective subcategory. Thus, some well known compactifications in the topology appear as this type of factorization.
Keywords: reflective and coreflective subcategories, the right product of two subcategories, τ-complete spaces. 2000 MSC: 18A20, 18 B30. INTRODUCTION
In the sequel we use the following notations: Eu (resp. Mu ) denotes the class of universal epi (resp. mono); e b Ep, (resp. Mp) the class of precise epi (resp. mono) : Ep= Mu, Mp= Eu; C2V - the category of Hausdorff locally convex topological vector space; Th - the category of Tikhonov spaces (the completely regular Hausdorff spaces); if K (resp. R) is coreflective (resp. reflective) subcategory, then k: C →R (r: C→R) is the coreflector (resp. reflector) functor; (Epi, M f ) - (the class of epimorphisms, the class of strict monomorphisms) = (the class of mappings with dense image, the class of topological inclusions with closed image); (E f , Mono) - (the class of strict epimorphisms, the class of monomorphisms) (Eu, Mp)-(the class of universal epimorphisms, the class of precise monomor- phisms)=(the class of surjective mapping, the class of topological inclusions). For concepts from general topology see [7], from topology of locally convex spaces see [8], and for those related to factorization structures see also [8]. The right product of a coreflective and reflective subcategory was introduced and examined in the paper [5]. Necessary and sufficient conditions for the product to be a reflective subcategory were identified. In the category C2V of Hausdorff locally convex topological vector spaces many cases when this product is a reflective subcat- egory were found. The examination of the right product of two subcategories is requested by the following situations:
41 42 Dumitru Botnaru, Alina T¸urcanu
1. The right product appears in natural way when studying the semi-reflexive subcategories [2]. 2. The relative torsions theories, that are so frequent in the category C2V, can be performed as right product theories. 3. Many reflective subcategories can be explained as the right product of two subcategories of certain type [3]. In the category C2V, as well as in the category of Tikhonov spaces, there are ex- amples when this product is not a reflective subcategory (Theorem 1.1). The properties of right product factorization are examined in this paper (Lemma 1.2). There are stipulated the conditions when this factorization defines a reflective subcategory (Theorem 1.3). In Section 2 it is proved that τ-complete spaces [6], [11] could be constructed as such factorizations (Theorem 2.3). In Section 3 it is shown how the subcategory of τ-complete spaces could be per- formed in two ways: either varying in the product the coreflective subcategory, or varying the factorization structure. In Section 4 there are formulated some issues for the category of Tikhonov spaces.
1. THE RIGHT PRODUCT OF TWO SUBCATEGORIES Let K be a coreflective subcategory, and R - a reflective subcategory of the cate- gory C. For any object X of the category C assume that rX : X → rX is R-replique of object X, and kX : kX → X and kkrX : krX → rX are the K-corepliques of respective objects. Then rXkX = krXg (1) for some morphism g. Since g = k(rX) we can write the preceding equality
rXkX = krXk(rX). (2)
We assume that in the category C pushout squares exist and we construct the pushout square on the morphisms kX and k(rX):
vXkX = gXk(rX). (3)
Then there exists a morphism uX so that
rX = uXvX, (4)
krX = uXgX. (5) The factorization of the right product of two subcategories 43
Figure 1.1
We denote by K ∗d R the full subcategory of all objects of category C, isomorphic with the objects of vX, X ∈| C | form.
Definition 1.1. The subcategory K ∗d R is called the d-product or the right product of subcategories K and R.
Theorem 1.1. Let C be a category with pull-back and pushout squares, R a monore- flective subcategory, K - a epireflective subcategory, and V = K ∗d R. Then, the following affirnations are equivalent: 1. V is a reflective subcategory of category C. 2. For any object X of category C the morphism vX is V-replique of object X. 3. For any object X of category C the morphism vX is an epi. 4. For any object X of category C the morphism uX is R-replique of object vX. 5. X ∈| V |⇔ rXkX is K-coreplique of object rX.
The proving is performed as in the case of category C2V ([5], Theorem 2.5). We mention that there are known various cases when subcategory V is reflective (see [5], Theorems 3.2.-3.4. and [3] Theorem 2.6.). Let us examine the right product of two subcategories in the category Th of Tikhonov spaces. Let K be a coreflective subcategory. Then it contains the subcategory D of the spaces with discreteT topology. It follows that it is a monoreflective subcategory, and then it is (Eu Mono) - coreflective, where Eu is class of universal epimorphisms (continued and surjective maps). It is obvious that in the category Th the reflective subcategory R is monoreflective iff when R includes the subcategory of compact spaces: Comp ⊂ R. Let Comp ⊂ R. We reffer to Figure 1.1. Let’s presume that V is a reflective sub- category. Then uX is Comp-replique of object vX. So, uX is a topological inclusion. And from the equality (5) it results that uX is a surjective application. Therefore, uX is an isomorphism. Thus, we proved:
Theorem 1.2. Let K be a coreflective subcategory, R - a reflective subcategory in the category Th and Comp ⊂ R. Then, two cases are possible: 1. K ∗d R = R. 2. K ∗d R is not a reflective subcategory of Th category. 44 Dumitru Botnaru, Alina T¸urcanu
But some examples of reflective subcategory of Th category show they can be obtained as a simple modification of the right product of two subcategories. In what follows we will describe this modification. Definition 1.2. The class of morphisms A of category C is called right-stable if, because g0 f = f 0g, it is a pushout square, and from f ∈ A it follows that f 0 ∈ A.
In the categories Th and C2V the pair (Eu, Mp) = (the class of surjective mor- phisms, the class of topological inclusions) is a factorization structure. Therefore the class Eu is right-stable. Lemma 1.1. Let C be a category with pull-back and pushout squares in which the class Eu is right-stable, K - a monoreflective subcategory, and R - a reflective sub- category. Then for any abject X of the category C, uX is a monomorphism. T Proof. Any monoreflective subcategory is (Eu Mono)-coreflective. We use the no- X tations from the beginning of the section. In the pushout square (3) k ∈ Eu and X rX according to the hypothesis that g ∈ Eu. In the equality (5) we have k ∈ Mono X X and g ∈ Eu. According to the Lemma 1.3 [5], we deduce that u ∈ Mono.
Corollary 1.1. In the categories Th and C2V for any two subcategories one coreflec- tive and other reflective, uX is always a monomorphism. Lemma 1.2. Assume in the category C and in the subcategory K and R, for any object X, uX is a monomorphism. Then: 1. For any object X of category C the morphism gX is K-coreplique of object vX : gX = kvX. 2. The correspondence X 7−→ vX defines a covariant functor v : C → C. 3. R is a monoreflective subcategory of the category K ∗d R. X 4. For any object X of the subcategory K ∗d R the morphism v is sectionable. Proof. 1. Consider f : A → vX, where A ∈| K |. Then uX f = krXg for some morphism g.
Figure 1.2 The factorization of the right product of two subcategories 45
We have uX f = krXg = uXgXg i.e. uX f = uXgXg, and since uX is a mono, we deduce that f = gXg. The uniqueness of the morphism g which satisfies the preceding equality results from the fact that krX is K-coreplique of object rX. So, we proved that gX = kvX. 2. Let us define the functor v on morphisms. Let f : X → Y ∈ C.
Figure 1.3 Then Y X r f = f1r , (6) for some morphism f1, where f1 = r( f );
rX rY f1k = k f2, (7) for some morphism f2, where f2 = k( f1). We have Y Y X Y X X X u v f k = (from (4), for object Y) = r f k =(from (6))= f1r k =(from (2)) rX X rY X Y vY X = f1k k(r ) = (from (7)) = k f2k(r )(= from (5) for object Y) = u k f2k(r ), i.e.
Y Y X Y vY X u v f k = u k f2k(r ). (8) 46 Dumitru Botnaru, Alina T¸urcanu
Since uY is a mono, from equality (8) we obtain
Y X vY X (v f )k = (k f2)k(r ). (9) From equality (9) and the fact that (3) is an pushout square, we deduce that
Y X v f = f3v , (10)
vY vX k f2 = f3k . (11)
We define v( f ) = f3. 3. Let us get look again to Fig.1.1. Let X ∈| R |. Then rX, and with him and k(rX) X and v , are isomorphisms. Therefore R ⊂ K ∗d R. Further on, let be rvX : vX → rvX R replique of object vX. Then uX = t · rvX X vX for some morphism t. Since u is a mono of category C , we deduceT that r is the same. The only thing we can add is: the subcategory R is (Epi Mu)-reflective in the category K ∗d R. 4. Let’s complete the diagram from Fig. 1.1 with an analogous diagram built for the object vX.
Figure 1.4 We have uX = tXrvX (12) for some morphism tX. Then tXkrvX = krXhX (13) for some morphism hX = k(tX). We have The factorization of the right product of two subcategories 47 krXhXk(rvX)= (from (13)) =tXkrvXk(rvX)=(from (2) for object vX) = tXrvXkvX =(from (12))= uXkvX = (from (10) since gX = kvx) = krX, i.e. krXhXk(rvX) = krX. (14) Therefore hXk(rvX) = 1. (15) Then in the pushout square vvXkvX = kvvXk(rvX) (16) the morphism k(rvX) is sectional. Therefore the morphism vvX is the same. We can assert that rvXvX = lXrX (17) for some morphism lX. We have tXlXrX = (from (17)) = tXrvXvX = (from (1)) = uXvX = (from (4)) = rX i.e. tXlXrX = rX (18) or tXlX = 1. (19)
We assume that the conditions of preceding lemma are satisfied, and (P, I) is a factorization structure in the category C. Let L be the full subcategory of category C comprising I-subobjects of objects of subcategories K ∗d R. For any object X of category C let vX = iXlX (20) be (P, I)-factorization of morphism vX. Theorem 1.3. The correspondence X 7−→ (lX, lX) defines the category L as a P-reflective subcategory of the category C. Proof. Let be A ∈| L |, and f : X −→ A ∈ C. We show that the morphism f is extended through morphism lX. According to the hypothesis there exists an object B ∈| K ∗d R | and a morphism i : A −→ B ∈ I. Then
vB(i f ) = v(i f )vX, (21) or
(vBi) f = (v(i f )iX)lX. (22) Since vB is sectional we deduce that vBi ∈ I, and lX ∈ P. Therefore lX ⊥ vBi, i.e. f = glX, (23) 48 Dumitru Botnaru, Alina T¸urcanu
v(i f )iX = vBig, (24) for some morphism g. The uniqueness of morphism g, that satisfies equality (23), results from the fact that lX is an epi.
Figure 1.5
2. THE SUBCATEGORY OF τ-COMPLETE SPACES We examine some coreflective subcategories of the category Th of Tikhonov spaces. Definition 2.1. Let τ be an cardinal. 1. In a topological space the intersection of τ open sets is called Gτ -set. − − 2. Pτ -change of the space (X, t) is called the space (X, tτ ), where the basis of − topology tτ is formed by Gα-set, α < τ. 3. Pτ-change of space (X, t) is called the space (X, tτ), where the basis of topology tτ is formed by Gτ-sets. We note with K−(τ) (respectively K(τ)) the full subcategory of the category Th − comprising all spaces (X, t) for which t = tτ (respectively, t = tτ). We observe that K(τ) = K−(τ+), where τ+ is the first cardinal which follows τ. If τ is limiting cardinal, then K−(τ) , K−(λ), for any cardinal λ. Therefore, it is enough to examine only the subcategories K−(τ). It is easily checked that K−(τ) (similar by and K(τ)) are the coreflective subcate- gories of category Th with coreflective functors.
− − − − Pτ : Th −→ K (τ), Pτ (X, t) = (X, tτ ),
Pτ : Th −→ K(τ), Pτ(X, t) = (X, tτ).
We mention the following properties of the subcategories K−(τ): 1. K−(τ) = Th for τ ≤ ω. 2. Let α < β be. Then K−(α) ⊇ K−(β). 3. Let Disc be the subcategory of discrete spaces. Then ∩{K−(τ)/τ} = ∩{K(τ)/τ} = Disc. Therefore, we can conclude that Disc = K−(∞) = K(∞), considering that τ < ∞ for every cardinal τ. The factorization of the right product of two subcategories 49
Theorem 2.1. ([4], Theorem1.2). Consider ω ≤ α < β. Then:
K−(β) ⊂ K−(α) and K−(β) , K−(α).
Definition 2.2. Let τ be an cardinal. The Tikhonov space X is called Q(τ)-space (respectively Q−(τ)-space), if X is closed in K(τ)-coreplique (respectively K−(τ)- coreplique) of space βX, where βX is Comp-replique of spaces X. We note with Q(τ) (respectively, Q−(τ) ) - the full subcategory of all Q(τ)-spaces (respectively, Q−(τ) -spaces). The categories Q(τ) have been studied by A. Cigoghidze [6], and the Q−(τ) by H. Herrlich [11]. − Let (X, t) be a Tikhonov space, (Y, u) = β(X, t), and (Y, uτ) − K (τ) -coreplique of space (Y, u). Let X be the closure of the set X in the space (Y, uτ). On set X we − induce the topology u0 out of space (Y, u). The topology space (X, u0) we note vt X, − or vt (X, t).
Figure 2.1
Evidently, Q(τ) = Q−(τ+) for any limiting cardinals τ. But for limiting cardinals τ, the subcategories Q−(τ) are of other form that Q(τ). For a subcategory A of category Th we note with PrA the subcategory that con- tains the products of objects of category A, and with M f (A) the subcategory that contains M f -subobjects of objects of subcategory A. Theorem 2.2. ([4], Theorem 2.6). 1. The subcategory Q−(τ) is a monoreflective subcategory (therefore also epireflective) of category Th with reflector functor − − vt : Th −→ Q (τ). 2. Q−(ω) = Q(n) = Comp, n ∈ N. − 3. Q (ωI) = Q(ω) = Q - the subcategory of Hewitt spaces. − τ 4.Q (τ) = M f Pr(R(τ)), where R(τ) = [−1; 1] \ {−1; 1}. − 5. Q (τ) = M f Pr(E(τ)), where E(τ) = Πα<τR(α). 6. α < β and ω < β. Then Q−(α) ⊂ Q−(β) and Q−(α) , Q−(β) . Corollary 2.1. Let τ be a limiting cardinal. Then
Q(τ) = M f Pr(∪Q(λ): λ > τ). 50 Dumitru Botnaru, Alina T¸urcanu
Remark 2.1. In [11] is defined the problem of existence of generators for subcate- − gories Q (τ), i.e. if there is a space Aτ, so that − Q (τ) = M f Pr(Aτ). Ignoring the case in [6], the problem is solved for subcategories Q(τ) and fully in the precedent theorem.
− Theorem 2.3. The reflective subcategory Q (τ) is (Epi, M f )-factorization of the − right product K (τ) ∗d Comp. Proof. For Tikhonov space (X, t) let (Y, u) = β(X, t) the Stone-Cech˘ compactification, X Y − and pτ :(X, tτ) −→ (X, t) and pτ :(Y, uτ) −→ (Y, u) the Kτ -coreplique of respective objects.
Figure 2.2 On the set Y we examine the inductive topology m, which is not mandatory being X Y the Tikhonov topology, defined by applications β and pτ : X −1 Y −1 G ∈ m ⇐⇒ [(β ) (G) ∈ t and (pτt) (G) ∈ uτ]. The square X X gpτ(β ) = hpτ , (25) X X is the pushout square constructed on morphisms pτ and pτ(β( )) in the category of topological spaces. In the sequel we construct the pushout square on these morphisms in the category Th. Let F be the set of continued defined functions (Y, m) with values in the field of real numbers R: F = Hom((Y, m), R), and m(F) - the topology defined on the set Y of this set F: The factorization of the right product of two subcategories 51
m(F) = { f −1(G)| f ∈ F and G is open in R}. Let l be the canonic application. Then
X X X (lh)pτ = g pτ(β ) (26) X X is the pushout square constructed on morphisms pτ and pτ(β ) in the category Th. Let Z be the closure of set X in the space (Y, m(F)), and let the topology v be the one induced from space (Y, m(F)). Since lX is an epi in the category Th it follows that uXiX is the Stone-Cech compactification of space (Z, v). Thus we can consider that the topology v is the one induced from application uXiX = βZ on set Z from space β(X, t). The closure of set X in the spaces (Y, m). A set A is closed in the space (Y, m) iff −1 −1 the set h (A) is closed in the space (X, t) and g (A) is closed in the space (Y, uτ). If −1 X ⊂ A, then the set A is closed iff the set g (A) is closed in the space (Y, uτ).
The closure of set X in the spaces (Y, m) and (Y, uτ) coincides: clmX = cluτ X. The last set we denote as T : T = cluτ X.
Let us prove that the closure of set X in the space (Y, uτ) coincides with Z: cluτ X = clm(F)X. Firstly we mention that T ⊂ Z and will prove the reverse inclusion. Let be y ∈ Y\T. Then exists a set H ⊂ Y \ T, that comprise point y so that: 1. H is closed in the topology u. 2. H is a Gτ-set in the topology u. 3. H is closed and open in the topology uτ. 4. H remains open and closed in the topology m. 5. H remains open and closed in the topology m(F). The theorem is proved.
3. THE CASE OF THE SUBCATEGORY Q−(τ) Let be R and L two reflective subcategories of the category C and L ⊂ R. If C is local and colocal small with projective limits, then there exists a class L(R) of factorization structure in the category C with the following property. For any object X of category C let be rX : X −→ rX and lX : X −→ lX the respective replique. Then lX = vXrX (27) for some morphism vX. If the subcategory L is monoreflective, then in the written equality all morphisms are bimorphisms. We note U = {rX|X ∈ |C| }, V = {vX|X ∈ |C| }. Always U ⊥ V. 52 Dumitru Botnaru, Alina T¸urcanu
Assume P00(R) = P00 = Vq; I00(R) = I00 = Vqx; P0(R) = P0 = Uxq; I0(R) = Ux. Theorem 3.1. [1] Let C be local and colocal small with projective limits and the subcategory L is monoreflective. Then: 1. (P0, I0) and (P00, I00) are structures of factorization in the subcategory C. 2. Let (P, I) be a structure of factorization in the category C. The following affir- mations are equivalent: a) for any object X of category C the equality (1) is (P, I)-factorization of mor- phism lX; b) P0 ⊂ P ⊂ P00. Relying on this theorem we assert
Theorem 3.2. Let τ be a cardinal, τ > ℵ0. 1. The reflective subcategory Q−(τ) could be obtained factorizing the right product X − (the morphism v ) K (τ) ∗d Comp following the structure of factorization (Epi, M f ). 2. In the category Th there ar the structures of factorization (P, I), P0 (Q−(τ)) ⊂ P ⊂ P00(Q−(τ)), so that the reflective subcategory Q−(τ) could be obtained doing the − (P, I)- factorization of right product K (ω) ∗d Comp. 4. PROBLEMS In paper [4] some classes of coreflective subcategories in the category Th are ex- amined.
− Definition 4.1. Let τ be an cardinal. The topological space (X, t) is called kτ -space, if any function definite on set X and continue on every compact K ⊂ X, | K |< τ is continue on space (X, t).
− − Let C (τ) a full subcategory of all kτ -spaces. Theorem 4.1. The subcategory C−(τ) is a coreflective subcategory of category Th. Definition 4.2. The weight w(X, t) of the topological space (X, t) is called minimal of cardinals |B|, where B is basis of spaces (X, t).
− Definition 4.3. Let τ be an cardinal. The Tikhonov space is called bτ -space if every function definite on set X and continue on every compact K ⊂ X, w(K) < τ is continue on space (X, t).
− − Let B (τ) is a subcategory full of all bτ -spaces. Theorem 4.2. 1. The subcategory B−(τ) is a coreflective subcategory of category Th. 2. Let be τ ≤ ω. Then B−(τ) = Disc. 3. ∩B−(τ) = C - subcategory of k-functionals space. The factorization of the right product of two subcategories 53
4. Let ω ≤ α < β and α a regular cardinal and β ≥ α+. Then B−(α) ⊂ B−(β) and B−(α) , B−(β).
Problem 4.1. 1. Let us describe the subcategory R which is (Epi, M f )-factorization − of the right product C (τ) ∗d Comp. 2. R-replique of an arbitrary object X of category Th is also the closing of set X in the C−(τ)-coreplique of space βX (with induce topology from βX)? − 3. The same thing is valid for the right product B (τ) ∗d Comp. References
[1] Botnaru D., Structure bicategorielles complementaires, ROMAI J., 5, 2(2009), 5-27. [2] Botnaru D., Cerbu O., Semireflexive product of two subcategories, Proc.of the Sixth Congress of Romanian Mathematicians, Bucharest, 2007, V.1, p. 5-19. [3] Botnaru D., Cerbu O., Some properties of semireflexivity subcategories (submitted). [4] Botnaru D., Robu R., Some categorical aspects of Tikhonov spaces, Scientific Annals. of MSU, The series ”Physical-Mathematical Sciences”, Chis¸inau,˘ 2000, 87-90 (in Romanian). [5] Botnaru D., Turcanu A., Les produits de gauche et de droit de deux sous-categories, Acta et coment, Chis¸inau,˘ III(2003), 57-73. [6] Cigoghidze A.C., About properties close of compactness, Uspehi mat.nauk, 35, 6 (216)(1980), 177-178. [7] Engelking R., General topology, Warszawa, 1985. [8] Grothendieck A., Topological vector spaces, Gordon and Breach, New York, London, Paris, 1973. [9] Herrlich H., Fortsetzbarkeit stetiger abbildungen und kompactheitsgrad topologischer Raume, Math. Z., 96, 1(1967), 64-72.
ROMAI J., 6, 2(2010), 55–67
CLASSIFICATION OF A CLASS OF QUADRATIC DIFFERENTIAL SYSTEMS WITH DERIVATIONS Ilie Burdujan University of Agricultural Sciences and Veterinary Medicine ”Ion Ionescu de la Brad” Ia¸si, Romania burdujan [email protected] Abstract The classification up to a center-affinity of a class of homogeneous quadratic differential systems defined on R3, which have a semisimple derivation with a 2-dimensional kernel is achieved. It is proved that there exist eighteen classes of affinely nonequivalent such systems.
Keywords: homogeneous quadratic dynamical systems, nilpotent derivation. 2000 MSC: Primary 34G20, Secondary 34L30, 15A69.
1. INTRODUCTION The study of a homogeneous quadratic differential system (S ) (briefly, HQDS) on a B space A can be achieved by using suitable algebraic tools. This is possible because a commutative algebra A(·) is naturally associated with any such a system (S ). In general, algebra A(·) is not associative. There exists (see [11], [1]) a 1-to- 1 correspondence between the classes of affinely equivalent HQDSs on A and the classes of isomorphic commutative algebras on A. A strong connection between the qualitative properties of HQDSs and the properties - invariant up to an isomorphism - of the corresponding commutative algebras emerges from this correspondence. In particular, any derivation/automorphism of (S ) is a derivation/automorphism of its associated algebra A(·) (see [5], [9]). As it is well known (see, for example, [1]), if a derivation of A(·) having only real eigenvalues exists, then its semisimple and nilpo- tent parts are necessarily derivations, too. The commutative algebras on R3 having a nilpotent derivation were already classified (see [4], [6]). The problem of clas- sification of 3-dimensional real commutative algebras having at least a semisimple derivation is an intricate task and it is still unsolved. In this paper, in order to save the space, we shall classify - up to a center-affinity - only a particular class of HQDSs on R3 which have a semisimple derivation. It is proved that this class of HQDSs is naturally partitioned in eighteen mutually disjoint subclasses consisting of nonequivalent systems, i.e., every system in one subclass is not equivalent to any system in another subclass; moreover, two systems in the same subclass that correspond to different values of parameters are nonequivalent, too.
55 56 Ilie Burdujan
Our results are based on the remark: the problem of classification up to an affine equivalence of HQDSs on Rn is equivalent with the problem of classification up to an isomorphism of all real n-dimensional commutative algebras. That is why we shall classify, up to an isomorphism, the class of real 3-dimensional commutative algebras that correspond to the analyzed class of HQDSs. The lattices of the subalgebras of these algebras are used as a main tool for solving the isomorphism problems. In a forthcoming paper we shall classify the set of all HQDSs on R3 which have a semisimple derivation.
2. REAL COMMUTATIVE ALGEBRAS ON R3 HAVING A SEMISIMPLE DERIVATION Let us consider a nontrivial real 3-dimensional commutative algebra A(·). Then, any basis in A allows to identify A with R3. Now, let us suppose that A(·) has a nonzero semisimple derivation. Then Der A necessarily contains a semisimple derivation D having its spectrum SpecD of the form SpecD = {1, λ, µ} (λ, µ ∈ R). Accordingly, there exists a basis B = (e1, e2, e3) such that D(e1) = e1, D(e2) = λe2, D(e3) = µe3. 2 2 2 Obviously, the case λ = µ = 1 is not possible. Since D(e1) = 2e1, then either e1 = 0 2 (if 2 < SpecD) or e1 , 0 (if 2 ∈ SpecD). More exactly, it follows 2 e1 = 0 i f 2 < SpecD 2 e1 = κe2 i f λ = 2 , µ 2 e1 = ωe3 i f µ = 2 , λ 2 e1 = κe2 + ωe3 i f λ = µ = 2. Each of the following equations can be exploited in a similar way: 2 2 D(e1e2) = (1 + λ)e1e2, D(e1e3) = (1 + µ)e1e3, D(e2) = 2λe2, 2 2 D(e2e3) = (λ + µ)e2e3, D(e3) = 2µe3. In fact, at least one of the elements of the set {2, 1 + λ, 1 + µ, 2λ, 2µ, λ + µ} belongs to SpecD. A simple inspection of papers [4], [6], [10] provides the following cases of the spectrum of D that have to be considered: {1, 0, 0}, {1, 1, 0}, {1, 1, 2}, {1, 2, 0}, {1, 2, 2}, {1, 2, 3}, {1, 2, 4}, {1, −1, 0}. In what follows we deal with a class of algebras which have a semisimple derivation D with SpecD = {1, 0, 0}, only. This is the case when dimR Ker D = 2.
3. CASE SPEC D = {1, 0, 0}
There exists a basis B = {e1, e2, e3} such that
D(e1) = e1, D(e2) = 0, D(e3) = 0. Classification of a class of quadratic differential systems with derivations 57
Then, from the equations
2 2 D(e1) = 2e1, D(e1e2) = e1e2, D(e1e3) = e1e3, 2 2 D(e2) = D(e2e3) = D(e3) = 0 emerges the following multiplication table of algebra A(·)
2 2 e1 = 0 e2 = αe2 + βe3 e1e2 = κe1 e2e3 = γe2 + δe3 2 e1e3 = ωe1 e3 = εe2 + θe3. It results that A = Im D ⊕ Ker D where Re1 is an ideal and Ker D is a subalgebra of A(·). The algebra Ker D can have (see [11]) either a basis of nilpotent elements of in- dex two or a single 1-subspace of such nilpotent elements, only. Among the algebras of the former kind, there exists one (see [11]) which has the following multiplication table: 2 2 e2 = 0, e2e3 = e3, e3 = γe2 + e3 (γ , 0). It suggests us to deal with algebras defined by the following multiplication table:
2 2 Table T e1 = 0 e2 = 0 e1e2 = αe1 e2e3 = e3 2 e1e3 = βe1 e3 = γe2 + e3 with α, β, γ ∈ R, γ , 0. We shall denote by A(α, β, γ) any algebra on R3 with the multiplication table T in basis B.
Algebra A(α, β, γ) Any algebra A = A(α, β, γ)(γ , 0) is not simple because it contains the ideal Im D. Distinct triples of real numbers (α, β, γ)(γ , 0) determine non-isomorphic algebras. Indeed, the following result is proved by a straightforward checking. Theorem 3.1. The algebras A(α, β, γ)(γ , 0) and A(p, q, r)(r , 0) are isomorphic if and only if α = p, β = q, γ = r. In order to find a partition for class A(α, β, γ)(γ , 0) into subclasses of noniso- morphic subalgebras we look for the family of its 1-dimensional subalgebras. To this end we need to find its annihilator elements, nilpotent (of index two) elements and idempotent (of index two) elements. This problem is solved in the next three propositions. Proposition 3.1. The set Ann (A) of all annihilator elements of A = A(α, β, γ) is: ( {0} i f α2 + β2 , 0 Ann (A) = Im D i f α = β = 0. 58 Ilie Burdujan
Proposition 3.2. The set N(A) of all nilpotent elements (of index two) of A = A(α, β, γ) is: ( Re ∪ Re i f α , 0 N(A) = 1 2 SpanR{e1, e2} i f α = 0. Proposition 3.3. The set of idempotent elements (of index two) of any algebra A(α, β, γ) (γ , 0) is: ◦ 1 1 I(A) = ∅ if γ < − 8 , ◦ 1 1 2 I(A) = {− 2 e2 + 2e3} if γ = − 8 and α − 4β + 1 , 0, ◦ 1 1 3 I(A) = {xe1 − 2 e2 + 2e3 | x ∈ R} if γ = − 8 and α − 4β + 1 = 0, ◦ 1 4 I(A) = {y1e2 + z1e3, y2e2 + z2e3 | x ∈ R} if γ > − 8 , γ , 0, 2(αy1 + βz1) − 1 , 0 and 2(αy2 + βz2) − 1 , 0, ◦ 1 5 I(A) = {y1e2 + z1e3} ∪ {xe1 + y2e2 + z2e3 | x ∈ R} if γ > − 8 , γ , 0, 2(αy1 + βz1) − 1 , 0 and 2(αy2 + βz2) − 1 = 0, ◦ 1 6 I(A) = {xe1 + y1e2 + z1e3, | x ∈ R} ∪ {y2e2 + z2e3} if γ > − 8 , γ , 0, 2(αy1 + βz1) − 1 = 0 and 2(αy2 + βz2) − 1 , 0, ◦ 7 I(A) = {xe1 + y1e2 + z1e3 | x ∈ R} ∪ {xe1 + y2e2 + z2e3 | x ∈ R} if 1 γ > − 8 , γ = 0, 2(αy1 + βz1) − 1 = 0 and 2(αy2 + βz2) − 1 , 0. 1 Corollary 3.1. Every algebra A(α, β, γ)(γ < − 8 ) is not isomorphic to any algebra 1 1 A(p, q, 8 ) or A(p, q, r)(r > − 8 , r , 0). Proof. It is enough to notice that any two such classes of algebras have distinct fam- ilies of idempotent elements.
Properties of algebra A(α, β, γ) In what follows, we list the properties of these eighteen subclasses of algebras that help us to solve easily the isomorphism problem. Moreover, note that the lattices of subalgebras of each subclass of algebras were obtained.
◦ 1 2 2 Properties of algebra A1 =A(α, β, γ)(α , 0, γ < − 8 , (α − 1) + (2β − 1) , 0) • AnnA = {0} • N(A) = Re1 ∪ Re2 • I(A) = ∅ • Subalgebras: Re1, Re2, SpanR{e1, e2}, SpanR{e2, e3} • Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1 Classification of a class of quadratic differential systems with derivations 59
◦ 1 1 Properties of algebra A2 =A(1, 2 , γ)(γ < − 8 ) • AnnA = {0} • N(A) = Re1 ∪ Re2 • I(A) = ∅ 2 2 • Subalgebras: Re1, Re2, SpanR{e2, ae1 + be3} (a + b , 0) • Ideals : Re1 2 • A = A 0 0 1 0 0 0 • Der A = RD ⊕ RD1 where [D1] = and [D, D1] = D1 0 0 0 x 0 y 0 1 0 • Aut (A) = {T(x, y) | x, y ∈ R, x , 0} where T(x, y) = . 0 0 1
◦ 1 Properties of algebra A3 =A(0, β, γ)(β , 0, γ < − 8 ) • AnnA = {0} • N(A) = SpanR{e1, e2} • I(A) = ∅ 2 2 • Subalgebras: R(pe1 + qe2)(p + q , 0), SpanR{e1, e2}, SpanR{e2, e3} • Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 Properties of algebra A4 =A(α, β, − 8 )(α , 0, α − 4β + 1 , 0) • AnnA = {0} • N(A) = Re1 ∪ Re2 1 • I(A) = {− 2 e2 + 2e3} 1 • Subalgebras: Re1, Re2, R(− 2 e2 + 2e3), SpanR{e1, e2}, S panR{e2, e3}, 1 SpanR{e1, − 2 e2 + 2e3} • Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1 60 Ilie Burdujan
◦ 1 1 Properties of algebra A5 =A(0, β, − 8 )(β , 0, β , 4 ) • AnnA = {0} • N(A) = SpanR{e1, e2} 1 • I(A) = {− 2 e2 + 2e3} 2 2 1 • Subalgebras: Re1, Re2, R(pe1+qe2)(p +q , 0), R(− 2 e2+2e3), SpanR{e1, e2}, 1 SpanR{e2, e3}, S panR{e1, − 2 e2 + 2e3} • Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 1 1 Properties of algebra A6 =A(4β − 1, β, − 8 )(β < { 4 , 2 }) • AnnA = {0} • N(A) = Re1 ∪ Re2 1 • I(A) = {xe1 − 2 e2 + 2e3 | x ∈ R} 1 • Subalgebras: Re1, Re2, R(xe1 − 2 e2 + 2e3)(x ∈ R), SpanR{e1, e2}, 1 SpanR{e2, e3}, SpanR{e1, − 2 e2 + 2e3} • Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 1 Properties of algebra A7 =A(1, 2 , − 8 ) • AnnA = {0} • N(A) = Re1 ∪ Re2 1 • I(A) = {xe1 − 2 e2 + 2e3 | x ∈ R} 1 • Subalgebras: Re1, Re2, R(xe1 − 2 e2 + 2e3)(x ∈ R), SpanR{e1, e2}, 1 SpanR{e2, e3}, S panR{e1, − 2 e2 + 2e3}, S panR{e2, ae1 + be3} (a, b ∈ R)}, 1 cx Span {xe − e + 2e , e + be + ce } (b, c ∈ R, 4b + c , 0) R 1 2 2 3 2 1 2 3 • Ideals : Re1 2 • A = A 0 0 1 0 0 0 • Der A = RD ⊕ RD1 where [D1] = and [D, D1] = D1 0 0 0 Classification of a class of quadratic differential systems with derivations 61 x 0 y 0 1 0 • Aut (A) = {T(x, y) | x, y ∈ R, x , 0} where T(x, y) = . 0 0 1
◦ 1 1 Properties of algebra A8 =A(0, 4 , − 8 ) • AnnA = {0} • N(A) = SpanR{e1, e2} 1 • I(A) = {xe1 − 2 e2 + 2e3 | x ∈ R} 2 2 • Subalgebras: R(pe1 + qe2)(p + q , 0), SpanR{e1, e2}, SpanR{e2, e3}, 1 ∗ SpanR{e1, − 2 e2 + 2e3}, SpanR{pe1 + qe2, e3}(p, q ∈ R ), 1 SpanR{− 2 e2 + 2e3, ae1 + be2 + ce3} (a, b, c ∈ R, 4b + c , 0), 1 SpanR{xe1 − 2 e2 + 2e3, (1 − 2xq)e1 + qe2} (x, q ∈ R, xq , 0) • Ideals : Re1 2 • A = A 0 4 1 0 0 0 • Der A = RD ⊕ RD1 where [D1] = and [D, D1] = D1 0 0 0 x 4y y 0 1 0 • Aut (A) = {T(x, y) | x, y ∈ R, x , 0} where T(x, y) = . 0 0 1
◦ 1 Properties of algebra A9 =A(α, β, γ)(α , 0, γ , 0, γ > − 8 , 2(αy1 + βz1) − 1 , 0, 2(αy2 + βz2) − 1 , 0) • AnnA = {0} • N(A) = Re1 ∪ Re2 • I(A) = {y1e2 + z1e3, y2e2 + z2e3} • Subalgebras: Re1, Re2, R(y1e2 + z1e3), R(y2e2 + z2e3), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{y1e2 + z1e3, e1}, SpanR{y1e2 + z1e3, e2}, SpanR{y2e2 + z2e3, e1}, SpanR{y2e2 + z2e3, e2}, SpanR{y1e2 +z1e3, ae1 +b[γ(α−z1 −αz1 +2βz1 +1)e2 +2(β+γz1 +αγz1)e3]} (a, b ∈ R), SpanR{y2e2 +z2e3, ae1 +b[γ(α−z2 −αz2 +2βz2 +1)e2 +(β+γz2 +αγz2)e3]} (a, b ∈ R) • Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 Properties of algebra A10 =A(0, β, γ)(γ , 0, γ > − 8 , 2βz1 − 1 , 0, 2βz2 − 1 , 0, β , 0) • AnnA = {0} 62 Ilie Burdujan
• N(A) = SpanR{e1, e2} • I(A) = {y1e2 + z1e3, y2e2 + z2e3} 2 2 • Subalgebras: R(pe1 + qe2)(p + q , 0), R(y1e2 + z1e3), R(y2e2 + z2e3), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{y1e2 + z1e3, e1}, SpanR{y2e2 + z2e3, e1}, SpanR{y1e2 + z1e3, ae1 + b[γ(1 − z1 + 2βz1)e2 + 2(β + γz1)e3]} (a, b ∈ R, γ , 1), SpanR{y2e2 + z2e3, ae1 + b[γ(1 − z2 + 2βz2)e2 + 2(β + γz2)e3]} (a, b ∈ R, γ , 1) There exists other subalgebras if α, β, γ are connected by some appropriate equa- tions, namely: p 2 1 SpanR{pe1 +qe2, q (b+βc)e1 +be2 +ce3} i f γ = 2β −β, β , 4 (b, c, p, q ∈ R, q , 0) • Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 Properties of algebra A11 =A(α, β, γ)(α , 0, γ , 0, γ > − 8 , 2(αy1 + βz1) − 1 , 0, 2(αy2 + βz2) − 1 = 0) • AnnA = {0} • N(A) = Re1 ∪ Re2 • I(A) = {y1e2 + z1e3} ∪ {xe1 + y2e2 + z2e3 | x ∈ R} • Subalgebras: Re1, Re2, R(y1e2+z1e3), R(xe1+y2e2+z2e3)(x ∈ R), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{y1e2 + z1e3, e1}, SpanR{y2e2 + z2e3, e1}, If x , 0 a biparametric family of 2-dimensional algebras of type SpanR{xe1 + y2e2 + z2e3, ae1 + be2 + ce3} is obtained; they correspond to the following conditions: 2 2 4x(1 − α)b + 4γx(βz2 − 1)c + 2(z2 − 2)ab + (4γz2 + z2 − 1)ac+ +4x(γz2 + αγz2 − β + 1)bc = 0, 3 2 2 2 −γxc + 2xb c + xbc − 2(y2 − αy2 + βz2)abc + (−y2 + 2βy2 + γz2)ac − 2 −2αz2ab = 0.
• Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 Properties of algebra A12 = A(0, β, γ)(γ , 0, γ > − 8 , 2βz1 −1 , 0, 2βz2 −1 = 2 0, i.e. γ = 2β − β, β = −γz1, β , 0) • AnnA = {0} • N(A) =SpanR{e1, e2} • I(A) = {y1e2 + z1e3} ∪ {xe1 + y2e2 + z2e3 | x ∈ R} Classification of a class of quadratic differential systems with derivations 63
• Subalgebras: R(pe1 + qe2), R(y1e2 + z1e3), R(xe1 + y2e2 + z2e3)(x ∈ R), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{y1e2 + z1e3, e1}, SpanR{y2e2 + z2e3, e1}, p Span{pe1 + qe2, q (b + βc)e1 + be2 + ce3} (b, c, p, q ∈ R, pq , 0), SpanR{y1e2 + z1e3, ae1 + b[γ(2βz1 − z1 + 1)e2 + 2(β + γz1)e3]} (a, b ∈ R, β , −1), 2β−1 SpanR{xe1 + y2e2 + z2e3, ae1 + b( 2 e2 + e3)} (x, a, b ∈ R), 4βa − a − 4β2 xb SpanR{xe1 + y2e2 + z2e3, ae1 + 4βx e2 + be3} (x, a, b ∈ R, x , 0) • Ideals : Re1 2 • A = A 0 1 −γz 1 0 0 0 • Der A = RD ⊕ RD1 where [D1] = and [D, D1] = D1 0 0 0 x y −γz y 1 ∗ 0 1 0 • Aut (A) = {T(x, y) | x ∈ R } where T(x) = . 0 0 1 ◦ 1 Properties of algebra A13 =A(α, β, γ)(α , 0, γ , 0, γ > − 8 , 2(αy1 +βz1)−1 = 2 2 2 0 , 2(αy1 +βz2)−1 , 0, (α−1) +(2β−1) , 0, i.e. 2γ(1−α) +(2β−1)(α−2β) = 0) • AnnA = {0} • N(A) = Re1 ∪ Re2 • I(A) = {xe1 + y1e2 + z1e3 | x ∈ R} ∪ {y2e2 + z2e3} • Subalgebras: Re1, Re2, R(y1e2+z1e3), R(xe1+y2e2+z2e3)(x ∈ R), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{y1e2 + z1e3, e1}, SpanR{y2e2 + z2e3, e1}, If x , 0 a 2-parametric family of 2-dimensional algebras of type S panR{xe1 + y1e2 + z1e3, ae1 + be2 + ce3} is obtained; they correspond to the following conditions:
2 2 4x(1 − α)b + 4γx(βz1 − 1)c + 2(z1 − 2)ab + (4γz1 + z1 − 1)ac+ +4x(γz1 + αγz2 − β + 1)bc = 0, 2 2 2 2 −γxc + 2xb c + xbc − 2(y1 − αy1 + βz1)abc + (−y1 + 2βy1 + γz1)ac − 2 −2αz1ab = 0.
• Ideals : Re1 • A2 = A • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 1 Properties of algebra A14 =A(1, 2 , γ)(γ , 0, γ > − 8 ) • AnnA = {0} • N(A) = Re1 ∪ Re2 • I(A) = {xe1 + y1e2 + z1e3 | x ∈ R} ∪ {y2e2 + z2e3} 64 Ilie Burdujan
• Subalgebras: Re1, Re2, R(xe1+y1e2+z1e3)(x ∈ R), R(y2e2+z2e3), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{e1, y1e2 + z1e3}, SpanR{e1, y2e2 + z2e3}, SpanR{e2, ae1 + be3} (a, b ∈ R), SpanR{xe1 + y1e2 + z1e3, −4xe1 + z1(e2 − 4e3)}, ∗ SpanR{y2e2 + z2e3, ae1 + b(y1e2 + z1e3)(a, b ∈ R ) • Ideals : Re1 2 • A = A 0 0 1 0 0 0 • Der A = RD ⊕ RD1 where [D1] = and [D, D1] = D1 0 0 0 x 0 y 0 1 0 • Aut (A) = {T(x, y) | x, y ∈ R, x , 0} where T(x, y) = . 0 0 1 ◦ 1 Properties of algebra A15 =A(0, β, γ)(γ , 0, γ > − 8 , 2βz1 − 1 = 0, 2βz2 − 1 , 2 0, i.e. γ = 2β − β, β = −γz2) • AnnA = {0} • N(A) = S panR{e1, e2} • I(A) = {xe1 + y1e2 + z1e3 | x ∈ R} ∪ {y2e2 + z2e3} • Subalgebras: R(pe1 + qe2)(p, q ∈ R), R(xe1 + y1e2 + z1e3)(x ∈ R), R(y2e2 + z2e3), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{y1e2 + z1e3, e1}, p SpanR{y2e2 +z2e3, e1}, Span{pe1 +qe2, q (b+βc)e1+be2+ce3} (b, c, p, q ∈ R, pq , 0), SpanR{y2e2 + z2e3, ae1 + b[γ(2βz2 − z2 + 1)e2 + 2(β + γz2)e3]} (a, b ∈ R, β , −1), 2β−1 SpanR{xe1 + y2e2 + z2e3, ae1 + b( 2 e2 + e3)} (x, a, b ∈ R), 4βa − a − 4β2 xb SpanR{xe1 + y1e2 + z1e3, ae1 + 4βx e2 + be3} (x, a, b ∈ R, x , 0) • Ideals : Re1 2 • A = A 0 1 −γz 2 0 0 0 • Der A = RD ⊕ RD1 where [D1] = and [D, D1] = D1 0 0 0 x y −γz y 2 ∗ 0 1 0 • Aut (A) = {T(x, y) | x ∈ R } where T(x) = . 0 0 1
◦ 1 Properties of algebra A16 =A(0, 0, γ)(γ < − 8 )
• AnnA = Re1 • N(A) = S panR{e1, e2} • I(A) = ∅ 2 2 • Subalgebras: R(pe1 + qe2)(p + q , 0), SpanR{e1, e2}, SpanR{e2, e3} • Ideals : Re1, SpanR{e1, e2} • A2 = A Classification of a class of quadratic differential systems with derivations 65
• Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 Properties of algebra A17 =A(0, 0, − 8 )
• AnnA = Re1 • N(A) = SpanR{e1, e2} 1 • I(A) = {− 2 e2 + 2e3} 2 2 1 • Subalgebras: R(pe1 + qe2)(p + q , 0), R(− 2 e2 + 2e3), SpanR{e1, e2}, 1 1 SpanR{e2, e3}, SpanR{e1, − 2 e2 + 2e3}, SpanR{e2, − 2 e2 + 2e3} • Ideals : Re1, SpanR{e1, e2} 2 1 • A = SpanR{− 8 e2 + e3, e3} • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1
◦ 1 Properties of algebra A18 =A(0, 0, γ)(γ , 0, γ > − 8 )
• AnnA = Re1 • N(A) = SpanR{e1, e2} • I(A) = {y1e2 + z1e3} ∪ {y2e2 + z2e3} 2 2 • Subalgebras: R(pe1 + qe2)(p + q , 0), R(y1e2 + z1e3), R(y2e2 + z2e3), SpanR{e1, e2}, SpanR{e2, e3}, SpanR{e1, y1e2 + z1e3}, SpanR{e1, y2e2 + z2e3} • Ideals : Re1, SpanR{e1, e2} 2 • A = SpanR{γe2 + e3, e3} • Der A = RD x 0 0 ∗ 0 1 0 • Aut (A) = {T(x) | x ∈ R } where T(x) = . 0 0 1 These lists of properties of classes Ai◦ allow us to assert the following result.
Theorem 3.2. Every algebra of class Ai◦ is not isomorphic to any algebra of class A j◦ (i , j, i, j ∈ {1, 2, ..., 18}).
Accordingly, the class of homogeneous quadratic differential systems x0 = 2αxy + 2βxz (S ) y0 = γz2 z0 = 2yz + z2 66 Ilie Burdujan is naturally divided into eighteen classes of affinely nonequivalent systems. The last two equations of (S) supply the prime integral F(y, z), namely 2y2 + yz − 2γz2 2γz − y ln + p 1 arctg p i f γ < − 1 4 8 y γ |8γ + 1| 2γy |8γ + 1| 2 (4y + z) 8y 1 F(y, z) = ln 4 − i f γ = − 8 4y z + 4y ! 1 p 2 2 √ |2y + yz − 2γz | 2 8γ+1 2γz − (1 + 8γ + 1)y 1 4 p i f γ > − . y 2γz − (1 − 8γ + 1)y 8 Consequently, any trajectory of (S) lies on a cylinder F(y, z) = const.
4. THE LATTICE OF SUBALGEBRAS OF A(α, β, γ) It is well known that isomorphic algebras have necessarily isomorphic lattices of their subalgebras. Consequently, the lattices of all subalgebras of two algebras be- come a natural tool in proving whether these algebras are or are not isomorphic. Moreover, the lattice of subalgebras of any algebra A delivers a partition of its vector space in nonempty mutually disjoint cells. In its turn, this partition provides a partition of the set of all solutions of the associated homogeneous quadratic differ- ential system. In the previous section we have obtained the lattice of all subalgebras of any algebra A(α, β, γ), even if, this finding is a customary task. For particular classes of algebras this problem can be easily solved, namely for algebras A with Der A = RD (i.e. algebras with fewer symmetries). For example, the lattice of subalgebras for algebra A1◦ is
{{0}, Re1, Re2, S panR{e1, e2}, S panR{e2, e3}, A}. These subalgebras define a partition of R3 by means of their ground subspaces Ox, Oy, xOy and yOz. Accordingly, any solution of the corresponding HQDS (S) has its orbit in one of the cells bounded by these subalgebras or it is contained in one of these subalgebras. Thus, some orbits lay either on the axis Ox or on Oy, other orbits are contained in the plane xOy or yOz (but have no intersection with Ox and Oy) and the other ones are contained inside of one of the four quadrants defined in R3 by planes xOy and yOz. The lattice of subalgebras for algebra A2◦ is 2 2 {{0}, Re1, Re2, S panR{e2, ae1 + be3} (a + b , 0), A}. A comparison of the two lattices allows us to decide that A1◦ and A2◦ are not isomor- phic. On the other hand, each solution of the corresponding HQDS (S) is contained in a subalgebra, i.e. it is necessarily a null torsion curve. Indeed, a straightfor- ward computation confirms us that the triple scalar product (r’, r”, r”’) = 0 where 2 2 r’ = (2xy + xz)e1 + (γz )e2 + (2yz + z )e3. Moreover, each trajectory of a HQDS associated to one of the algebras A16◦, A17◦, A18◦ has null torsion, too. Classification of a class of quadratic differential systems with derivations 67 5. CONCLUSIONS The analyzed class of commutative algebras is divided into eighteen classes of non-isomorphic algebras. Accordingly, the corresponding class of homogeneous quadratic differential systems is partitioned into eighteen classes of affinely nonequiv- alent systems. Since A is a vector direct sum of an ideal and a subalgebra, each such a system is decoupled into subsystems that can be solved. The most part of classes of algebras have Der A = RD, so that the correspond- ing orbits of solutions have just a few symmetries. The algebras in the classes A2◦, A7◦, A8◦, A12◦, A14◦, A15◦ are the only algebras that have two-dimensional derivation algebras.
References
[1] I. B, Quadratic differential systems, Publ. House PIM, 2008.(Romanian) [2] I. B, A classification of a class of homogeneous quadratic dynamical systems on R3 with derivations, Bull. I.P. Ias¸i, Sect. Matematica, Mecanica teoretica, Fizica, T.LIV (2008), 37–47. [3] I. B, Homogeneous quadratic dynamical systems on R3 having derivations with complex eigenvalues, Libertas Mathematica, XXVIII (2008), 69–92. [4] I. B, Classification of Quadratic Differential Systems on R3 having a nilpotent of Order 3 Derivation, Libertas Mathematica, XXIX (2009), 47-64. [5] I. B, Automorphisms and derivations of homogeneous quadratic differential systems, Ro- mai J., 6, (12010), 15-28. [6] I. B, Classification of Quadratic Differential Systems on R3 having a nilpotent of Order 2 Derivation (in press). [7] T. D, Classifications and Analysis of Two-dimensional Real Homogeneous Quadratic Differ- ential Equation Systems, J. Diff. Eqs., 32(1979), 311–334. [8] I. K, Algebras with many derivations, ”Aspects of Mathematics and its Applications” (ed. J.A. Barroso), Elsevier Science Publishers B.V., 1986, 431. [9] K. M. K , A. A. S, Quadratic Dynamical Systems and Algebras, J. of Diff. Eqs, 117 (1995), 67–127 . [10] K. M. K , A. A. S, Automorphisms ans Derivations of Differential Equations and Al- gebras, Rocky Mountain J. of Mathematics, 24, 1,(1994), 135–153. [11] L. M, Quadratic Differential Equations and Non-associative Algebras, in ”Contributions to the Theory of Nonlinear Oscillations” Annals of Mathematics Studies, no. 45, Princeton Uni- versity Press, Princeton, N. Y., 1960. [12] A.A. S, R. W, Introduction to Lie groups and Lie algebras, Academic Press, New York, 1973. [13] R. D. S, An Introduction to Nonassociative Algebras, Academic Press, New York, London, 1966. [14] N. I. V, K. S. S˘ı, Geometrical Classification of Quadratic differential systems (Rus- sian), Differentialnye Uravnenje, 13, 5(1977), 803–814.
ROMAI J., 6, 2(2010), 69–81
ON THE 3D-FLOW OF THE HEAVY AND VISCOUS LIQUID IN A SUCTION PUMP CHAMBER Mircea Dimitrie Cazacu1, Cabiria Andreian Cazacu2 1Polytechnic University of Bucharest, Romania 2”Simion Stoilow” Institute of Mathematics of the Romanian Academy, Bucharest, Romania [email protected], [email protected] Abstract We present the essence of the theoretical and experimental researches, developed in the field of the three-dimensional flow with free surface of a heavy and viscous liquid in the suction chamber of a pump, with or without air carrying away or appearance of the dangerous cavitational phenomenon, problem of a special importance for a faultless working of a pump plant, especially of these relevant to the nuclear electric power plants, which constituted the object of our former preoccupations. From the mathematical point of view we present the solving particularities of the two-dimensional boundary problems, specific to the three-dimensional domain limits of the pump suction chamber flow, as well as the possibility to obtain a stable numerical solution. From the practical point of view, we show the possibilities of laboratory modelling of this complex phenomenon, by simultaneously use of many similitude criterions, as well as the manner to determine the characteristic curves, specific to each designed pump suction chamber. Because the drawings of the pump station of Romanian nuclear power station Cer- navoda on Danube are classified, we shall present the theoretical and experimental re- searches performed on the model existent in the Laboratory of New Technologies of Energy Conversion and Magneto-Hydrodynamics, founded after the world energetic cri- sis from 1973 [1] in the POLITEHNICA University of Bucharest, Power Engineering Faculty, Hydraulics, Hydraulic Machines and Environment Engineering Department.
Keywords: liquid three-dimensional flow with free surface, three-dimensional flow in the pump suction chamber, complex modelling of liquid three-dimensional flow with free surface. 2000 MSC: 35 G 30, 35 J 25, 35 J 40, 35 Q 30, 65 N 12, 76 D 17.
1. INTRODUCTION The theoretical research concerning the three-dimensional flow of the real liquid in a suction chamber of a pump, [2], [3], is justified by the flow effect with or with- out appearance of whirls, with or without air training [4] or by the appearance of the cavitation destructive phenomenon at the liquid entry in a pump, as well as by their negative influence on the hydraulic and energetic pump efficiency, the quick cavita- tional erosions, the machine vibrations followed by the fast wear, eventually by the
69 70 Mircea Dimitrie Cazacu, Cabiria Andreian Cazacu pump suction less and by their pumping cease, extremely dangerous for the pumping station, for instance of a nuclear electric power plant [1]. These situations may come in the sight in special work conditions linked with the low level of the water in the river bed or in the pump suction chamber, or by an inefficient design of this installation for different cases of pumps working [5], [6].
2. NUMERICAL INTEGRATION OF THE HEAVY AND VISCOUS LIQUID THREE-DIMENSIONAL FLOW IN A PUMP SUCTION CHAMBER We present the mathematical specific boundary conditions [2], [3], which inter- vene in the problem of the three-dimensional viscous liquid steady flow into the suc- tion chamber of a pump plant, represented in the Figure 1. Thus, for the uniform flow of the viscous and heavy liquid in the entrance and go- ing out sections of the domain, we must consider the two-dimensional special bound- ary problems: - for the fluid entrance section we must solve a mixed Dirichlet-Neumann prob- lem for a Poisson equation with an unknown a priori constant, depending of the slope angle α of the suction chamber bottom, which can be determined by a successive cal- culus cycle of the computer program, satisfying a normalization condition, which, for the given value of mean input velocity and entrance section dimensions, ensures the desired flow-rate, - for the fluid going out pipe section we must solve a Dirichlet problem for another 0 Poisson equation, having the pressure drop gradient pz in the suction pipe also as unknown a priori constant, which can be determined also by a successive calculus cycle satisfying a normalization condition and which establishes the pressure drop gradient on the pipe length for the same flow-rate, satisfying another normalization condition.
2.1. PARTIAL DIFFERENTIAL EQUATIONS OF A VISCOUS AND HEAVY LIQUID STEADY FLOW The three-dimensional motion equations in Cartesian trihedron (Fig.1) are
0 0 0 0 1 0 U + U U + U V + U W + P = ν(U” 2 + U” 2 + U” 2 ) + g sin α, (1) T X Y Z ρ X X Y Z
0 0 0 0 1 0 V + V U + V V = V W + P = ν(V” 2 + V” 2 + V” 2 ), (2) T X Y Z ρ Y X Y Z
0 0 0 0 1 W + W U + W V + W W + = ν(W” 2 + W” 2 + W” 2 ) + g cos α, (3) T X Y Z ρ X Y Z On the 3D-flow of the heavy and viscous liquid in a suction pump chamber 71 the mass conservation equation for an incompressible fluid being 0 0 0 UX + VY + WZ = 0, (4) where, as usual, (U, V, W) is the velocity of the fluid, P is its the pressure, ρ - the fluid density, and ν - its kinematical viscosity. T and X, Y, Z represent time and spatial independent variables, respectively.
2.2. THE DIMENSIONLESS FORM OF THE EQUATION SYSTEM For more generality of the numerical solving [3], [5] - [8] we shall work with dimensionless equations, by using the characteristic physical magnitudes of flow: - the suction chamber width B H, approximately equal with its height, - the liquid mean velocity of the entrance in suction chamber Um = Q/BH, - the air atmospheric pressure on its free surface P0, - and the period T0 of whirl appearance in the three dimensional flow. With the new dimensionless variables and functions: X Y Z T U V W P x = , y = , z = , t = , u = , v = , w = , p = , (5) H H H T0 Um Um Um P0 the partial differential equations (1) to (4) become in the dimensionless form:
0 0 0 0 0 1 sin α 0 Sh ut + uxu + uyv + uzw + Eu px = ∆ u + , (1 ) Re x,y,z Fr 72 Mircea Dimitrie Cazacu, Cabiria Andreian Cazacu
0 0 0 0 0 1 0 Sh vt + vxu + vyv + vzw + Eu py = ∆ v, (2 ) Re x,y,z
0 0 0 0 0 1 cos α 0 Sh wt + wxu + wyv + wzw + Eu pz = ∆ w + , (3 ) Re x,y,z Fr in which one puts into evidence the following criteria of hydrodynamic flow simi- 2 larity, numbers of: Strouhal Sh = H/T0Um, Euler Eu = P0/ρUm, Reynolds Re = 2 HUm/ν and Froude Fr = Um/gH; the mass conservation equation is an invariant 0 0 0 0 ux + vy + wz = 0. (4 ) 2.3. THE NUMERICAL SOLVING METHOD The numerical integration of partial differential equation system was performed by an iterative calculus of unknown functions u, v and w, given by the algebraic rela- tions associated to the partial differential equations, by introducing their expressions deduced from the finite Taylor’s series developments [8] in a cubical grid, having the same step δx = δy = δz = χ:
0 f1,2,5 − f3,4,6 f1,2,5 − 2 f0 + f3,4,6 f = and f ”2 2 2 = . (6) x,y,z 2χ x ,y ,z χ2 2.4. THE BOUNDARY CONDITIONS FOR THE THREE-DIMENSIONAL STEADY FLOW - on the solid walls, due to the molecular adhesion condition, we consider
uSW = vSW = wSW = 0 (7) and consequently, from equation (3’) for α ≈ 0 the hydrostatic pressure distribution becomes γH pSW(z) = 1 + z = 1 + Ar z, (8) P0 in which we introduced the hydrostatic similarity number Ar = γH/P0, dedicated to the Archimedean lift discoverer, - on the free surface, considered to be a plane due to the liquid important weight and small flow velocities, we have the conditions
pFS(x, y, 0) = 1 and wFS(x, y, 0) = 0, (9) excepting a null measure set of points constituted by the whirl centers on the free surface, in whose neighbourhood the fluid has a descent motion (Fig. 1). Neglecting the liquid friction with the air, we shall cancel both shearing stress components, obtaining: 0 0 0 τzx |FS= uz + wx = uz |FS= 0, → u6 = u5 (10) 0 0 0 τxy |FS= vx + wy = vz |FS= 0, → v6 = v5. On the 3D-flow of the heavy and viscous liquid in a suction pump chamber 73
The value w6 is calculated from the mass conservation equation (4’), written in finite differences (that is not unstable in this local appliance [5] [6]) u − u v − v w − w 1 3 + 2 4 + 5 6 = 0 → w | = w + u − u + v − v , (11) 2x 2x 2x 6 FS 5 1 3 2 4 - in the entrance section, we consider the uniformly and steady flow at the normal depth of the current parallel with the channel bottom slope i = sin α, which leads us to the condition v |ES= w |ES= 0 and to the hydrostatic pressure repartition, the dis- tribution of u(y, z) velocity component being obtained from the first motion equation (1’), which in these boundary conditions become a Poisson type equation with the constant α a priori unknown, written also in finite differences
Re 1 X4 χ2 Re ∆ u + sin α = 0 → u0 = ui + sin α. (12) x,z Fr 4 4 Fr i=1 The numerical solving is possible by an iterative cycle on the computer, until to the obtaining of slope value i, necessary to the given Re and Fr numbers corresponding to the flow velocity Um, when the heavy force component g sin α is balanced by the interior friction forces, acting by the liquid adhesion to the channel bottom and its solid walls. The boundary conditions for the Poisson equation (12) for solving in the frame of plane the mixed problem Dirichlet-Neumann
USW = 0 on the solid walls, 0 =zx|FS= µUz = 0 on the free surface, lead us to solve a mixed problem Dirichlet-Neumann, whose normalization con- dition is ZB ZH Ui(Y, Z)dYdZ = Q = BHUm, (13) 0 0 or, in the iterative numerical calculus for the dimensionless case,
j −1 j −1 j −1 kXM−1 XN χ2 XN χ2 XN χ2 kXM−1 χ2 u + u + u + u = q = 1, (130) k,j 2 kM,j 8 k=2,j 4 k,j=2 k=2 j=2 j=2 j=2 k=2 representing the flow-rate which is brought by this velocity repartition in the entrance section, for the given values of Um, therefore Re and Fr numbers, as well the consid- ered grid step χ. For the starting of the calculus it is necessary to introduce only a single velocity value in a point as a seed, for instance uS = 2 in the point j = 20 and k = 40 in the middle of the free surface, in the rest of the domain the initial arbitrary values being 74 Mircea Dimitrie Cazacu, Cabiria Andreian Cazacu equal with zero, or, by admittance of a paraboloidal repartition a the initial given arbitrary values concerning the velocity distribution in the domain, in the shape
2 1 1 2 1 2 uva(y, z) = 8(1 − z ) · ( − y) · ( + y) = 8(1 − z ) · ( − y ), 2 2 4 (14) 2 uva( j, k) = 4us · a ( j − 1) · (k − 1) · [2(k − 1) · a][2( j − 1) · a], because we have 1 y = ( j − 1)a − and y = 1 − (k − 1)a, (15) 2 - in the exit section, we shall consider that the streamlines are parallel with the vertical pipeline walls, uES = vES = 0, the velocity component w(x, y) distribution in the uniformly and steady flow verifying also a Poisson type equation with an a priori 0 unknown constant pz, but this time due to the ignorance of the pressure drop on the 0 vertical pipeline pz, deduced from the third motion equation (3’) in approximation cos α ≈ 1, determined above by solving the entry section problem
X4 0 Re 1 2 1 0 Re Eu pz = cos α + ∆ w, → w0 = [ wi + χ Re( − Eu pz)]. (16) Fr x,y 4 Fr i=1 The boundary condition to solve the Poisson equation (16) in the frame of Dirichlet problem, is given by the liquid adhesion wES = 0 on the interior solid wall of suction pipe. The Poisson type equation solving, in which one does not know a priori the pres- 0 sure drop pz along the pipe, can be made by numerical way, also by an iterative 0 cycle on computer, yielding up to the arbitrary values pz & 0.1, which exceed a few the pressure distribution deduced from the hydrostatic repartition (8), until the ve- locity integration wES in the exit section (Fig. 2) shall equalize the unitary value of dimensionless flow-rate, calculated in the case of non-symmetrical or symmetrical supposed velocity repartitions, using the forms of [6] X20 X10 1 X X9 1 X20 X q ≈ χ2 · w + · w + · w = 1, nonsym ij 8 ij 8 ij i=17 j=7 i∈{17, 2} j=8 i=18 j∈{7, 1} (17)
2 qsym ≈ χ · [9w17, 8 + 4(w18, 8 + w17, 7)] = 1. 2.5. CONCLUSIONS ON THE USED BOUNDARY SPECIFIC PROBLEMS IN THE 3D-FLOWS The Dirichlet problem, proposed by the German mathematician Peter Gustav Lejeune Dirichlet, being generally specific to the equations of elliptic type, corre- sponding to the repartition of scalar physic magnitudes given by harmonic functions On the 3D-flow of the heavy and viscous liquid in a suction pump chamber 75
(for example of the temperature distribution in a homogeneous medium, without con- vection i.e. U = V = 0, and without interior heat sources in virtue of Laplace equa- tion) has here a specific character. The Neumann problem formulated by the German physicist Carl Gottfried Neumann, consists in giving, on the boundary C of the domain D, the 0 values of the normal derivative of the function, fn.
To specify the unknown function values in the domain, it is necessary also in this case to fulfil the so-called normalization condition, giving the flux through the frontier line
Z 0 φ = fndl. (18) C
In the three-dimensional case of mixed boundary problems, by combination of these two classical problems, to solve the Poisson equations with unknown constants, intervene the well-known normalization conditions, even for the Dirichlet problem only. Initially, the computational program has been made in Cobol [4]], due to the great number of the data, working successively to relax the three velocity components values in a plane, considering their values also in the two neighbouring planes, the graphical representation of the three-dimension flow being shown by the three veloc- ity components as in Figure 2. 76 Mircea Dimitrie Cazacu, Cabiria Andreian Cazacu 2.6. THE IMPORTANCE OF THE COMPLEX PHENOMENA FROM A PUMP SUCTION CHAMBER FLOW The faultless working of the pump plant for the nuclear electric power stations, which must work even in the case of nuclear reaction stopping, yet 10 to 12 hours till the complete cooling of the reactor, is the main aim to ensure the nuclear security. In the second part of this paper we present the experimental research performed, regarding the complex modelling of the spatial flow and the specific phenomena, which can be produced in the pump suction chamber. The research has been started in the year 1977 together with my former student and then collaborator Dr. Eng. Costel Iancu [2] for the pump station design of the nuclear power station Cernavoda in Romania [1] and continued by the construction of the Model installation endowed with the measuring apparatus for a good working research of a pump suction (Fig. 1 to 3) in the Laboratory of new technologies of energy conversion and environment protection, realized in the year 1979 with the help of the former Eng. PhD student Dumitru Ispas [3]-[6].
3. THE COMPLEX MODELLING OF THE HEAVY AND VISCOUS LIQUID 3D- FLOW IN A PUMP SUCTION CHAMBER After doctoral Thesis defending in 1957, M.D.Cazacu published a work on the lab- oratory complex modelling of the complicated water-hammer phenomenon [9], con- sidering not only the dimensionless motion equation, but also the mass conservation and state equation of the compressible liquid, as well as the pipe elastic deformation with influence upon the sound propagation velocity. He decided at that time to pay attention to this very important phenomenon, taking into consideration the high security level of a pump plant for a nuclear electric power plant, since the research on the laboratory model was not conclusive in the consulted technical literature in this field.
3.1. THE EXPERIMENTAL INSTALLATION DESCRIPTION AND MEASURING APPARATUS The experimental installation is presented in Fig. 3, and it consists in 1-stained glass chamber, 2-supporting feet, 3-tightening cap, 4-air tap, 5-flow reassure, 6- cen- trifugal pump, 7-direct current electric motor, 8-pump pressing pipeline, 9-discharge control valve, 10-flow meter, 11-pump suction pipeline, 12-suction basket, 13-emptying tap, 14-air carrying away whirl, 15-whirl on the bottom, 16-blades against the whirls. On the 3D-flow of the heavy and viscous liquid in a suction pump chamber 77
The water level in the pump suction chamber may be changed through the tap 13, the water leaving the chamber through the pump suction pipe 12 producing two whirls without or with air carrying away 14, which due to the rotational flow ω caused by the pump suction flow rates, smaller then the nominal flow rate, change their positions: one remaining on the free surface, but at the other chamber side the second removing on the chamber bottom, where with the solid particles uncoupled from the concrete and concentrated in the whirl axis, due to the bigger flow velocities, leads at a strong erosion of the pump suction chamber.
For the measuring of different physical magnitudes we utilized the apparatus: - a mechanic tachometer contacted with the electric motor shaft end, to measure the pump number of rotations n, - a diaphragm to measure the pump flow rate Q, - a graduated rule from transparent plastic to measure the water level Z in the suction chamber, fixed at the suction chamber glass, - a manual chronometer to measure the frequency of whirl appearance, - a precision metallic mano-vacuum-meter or a mercury manometer tube to mea- sure the air pressure at the suction chamber free surface.
3.2. THE FLOW COMPLEX MODELLING USING MORE SIMILARITY CRITERIONS The short history concerning these inefficient foregoing modelling attempts began any decades before, when the researchers tried to utilize the very known Froude’s criterion for the modelling of this gravitational flow in a pump suction chamber, but which gave on the Model, smaller then the Reality, built at the length scale 78 Mircea Dimitrie Cazacu, Cabiria Andreian Cazacu
λ = H/H0 = D/D0 > 1, lesser velocities at which the whirls could not to appear
V2 V02 V Q Fr = = Fr0 = → β = = λ0,5 ≥ 1 and k = = λ2β = λ2,5, (19) gH g0H0 V0 Q0 in which we noted with β the velocity ratio and with k the flow rate ratio. For this reason, in the last decade of that period, the researchers, without no justi- fication, have enlarged the velocities on the Model till the same values as the veloc- ities of the Reality, introducing the so called equal velocity criterion, for which M. D. Cazacu gave a proof considering Weber’s criterion, thinking at the free surface whirls with air carrying away similarity σ σ0 We = = We0 = → β = 1 = λ0 and then k = λ2, (20) ρV2 ρ0V02 where σ is the air-water superficial tension. As a good expert in the recommendations of the pump manufacturers and pump station designers, that the water access to the pump suction baskets must be as possi- ble direct and because the estimation of the separation of viscous fluid flow is given by Reynold’s number, he considered also this kind of modelling, supposing that the water cinematic viscosity is the same VD V0D0 V υ Re = = Re0 = → β = = λ−1 ≈ λ−1 and k = λ, (21) υ υ0 V0 υ0 in which the velocities on the Model are greater than those from the Reality. Another very important problem was the modelling of the influence of pump work- ing, especially at lesser flow rates than that nominal. He introduced the pump mod- elling by its rapidity criterion nq
nQ1/2 n0Q01/2 n = = n0 = , (22) q H3/4 q H03/4 from which we can choose the pump number of rotations on the laboratory Model
k1/2 n0 = n = nλ1/4β1/2 (23) λ3/4 and determine the convenient geometric scale in the case of composed similarity 0 Froude and pump rapidity nq, in which case he gave in Table 1 the pressure P0 nec- essary at the water free surface in suction chamber, to benefit also of the cavitation phenomenon modelling at pump or at the strong erosive whirl, which ends on a side wall or on the bottom of the suction chamber built from the concrete, in the aim to eliminate his appearance in the Reality [5]. On the 3D-flow of the heavy and viscous liquid in a suction pump chamber 79
Table 1
n’[rot/min] 3.000 1.500 1.000 750 600 n = 500 λ = (n0/n)2 36 9 4 2, 25 1, 44 1 P0 − P0 (T 0) 0 v = P0 − Pv(T) = 1/β2 = 1/λ = σ/σ0 0, 0277 0, 111 0, 25 0, 444 0, 6944 1 0 P0(bar) 0, 261 0, 342 0, 43 0, 5778 0, 7678 1
In this table one sees the advantage of the length scale λ = 4, corresponding to the standardized number of rotations 1,000 rot/min of the pump model, inclusively under pressure that we must assure at the model free surface, as well as that which assures only the cavitational phenomenon similarity, for which we must have also accomplished the equality relation of the cavitational criterion
P − P (T) P0 − P0 (T 0) σ = 0 v = σ0 = 0 v . (24) ρν2/2 ρ0ν02/2
Finally, in the case of Strouhal’s similarity, concerning the whirl appearance fre- quency in time, we shall have the relation, from which we can deduce the frequency of whirl appearance in reality as function of that obtained on the model
D f D0 f 0 D0V Sh = = Sh0 = → f = f 0 . (25) V V0 V0D
The symmetrical position (Fig. 1) of the two whirls from the right part can be modified by the pump working at littler flow rates (Fig. 3), by ω rotation impressed to the chamber flow as in the left part, when one whirl arrive to us and the other being disposed vertically on the bottom chamber (Fig. 4), causing its accelerated wear.
3.3. CONCLUSIONS AND IMPORTANT REMARKS With this occasion we shall observe that the greater models become too expensive and the smaller models introduce forces of other nature, which modify completely the development of the hydrodynamic process. On the basis of this modelling method one sees, that we can study a complex hydro-aerodynamic similarity using 6 criterions, as: Froude, Weber, Reynolds, pump rapidity, cavitation criterion and Strouhal, not only individually but also simultane- ously. 80 Mircea Dimitrie Cazacu, Cabiria Andreian Cazacu
In this way we can also determine the convenient technical-economically length scale of the Model, to have also a good modelling similarity. In the same time we can determine the whirl apparition frequency in reality, that is very important because the erosion of the pump suction chamber walls is however a problem of time and we can also determine the zones that can be covered with metallic plates and one can determine for any pump plant, working at different exploitation situations its characteristic curves (Fig. 4), which determine the region of a good operation.
References
[1] M.D.Cazacu, C.Iancu, Studiu privind aspirat¸ia ¸siaduct¸iunea pompelor mari de circulat¸ie din termocentrale CNE ¸sia criteriilor de similitudine pentru verificarea pe model (Study regarding the suction and delivery of the great circulation pumps from the Electric Nuclear Power Sta- tions and the similarity criterions for the check on the model). Contr. nr.201/1977, Beneficiary − ICEMENERG − Bucharest. [2] M.D.Cazacu, Gh.Baran, A.Ciocanea, High education and research development concerning the recoverable, inexhaustible and new energy sources in the POLITEHNICA University of Bucharest. Internat. Conf. on Energy and Environment, 23-25 Oct.2003, Polytechnic Univ. of Bucharest, Sect. 5 - Energetic and Educational Politics, Ses.2 - Educational and Environmental protection Politics, 5, 120-125. [3] M.D.Cazacu, On the boundary conditions in three-dimensional viscous flows. The 5th Congress of Romanian Mathematicians, June 22-28, 2003, University of Pitesti, Romania, p. 24. [4] M.D.Cazacu, Boundary mathematical problems in viscous liquid three-dimensional flows. Inter- nat. Conf. on Theory and Applications of Mathematics and Informatics - ICTAMI 2007, Alba Iulia, Acta Universitatis Apulensis, Mathematics-Informatics, No. 15 / 2008, 281 - 289. [5] C.Iancu, Contribut¸ii teoretice ¸siexperimentale la studiul curgerii fluidelor cu mi¸scare de rotat¸ie ¸sivˆartejlegat de suprafet¸e libere (Theoretical and experimental contributions at the study of fluid flows with rotation motion and whirl bound of free surfaces). Doctoral Thesis defended in 1970 at the Polytechnic Institute of Bucharest, Romania. On the 3D-flow of the heavy and viscous liquid in a suction pump chamber 81
[6] M.D.Cazacu, D.Ispas, Cercet˘ariteoretice ¸siexperimentale asupra curgerii lichidului vˆıscosˆın camera de aspirat¸ie a unei pompe (Theoretical and experimental research on the viscous liquid flow in a pump suction chamber), Conf. of Hydraulic Machines and Hydrodynamics, Timisoara, Romania, 18 - 19 Oct. 1985, vol. III, 67-74. [7] M.D.Cazacu, D.A.Ispas, M.Gh.Petcu, Pipe-channel suction flow researches. Proc. of XXth IAHR Congress (Subject B.a ), Moscow SSRU, September 5 - 9, 1983, VI, 214 - 221. [8] M.D.Cazacu, On the solution stability in the numerical integration of non-linear with partial differential equations. Proceedings of the Internat. Soc. of Analysis its Appl. and Computation, 17-21 Sept. 2002, Yerevan, Armenia. Complex analysis, Differential equations and Related top- ics. Publishing House Gitityun, 2004, Vol. III, 99-112. [9] D.Dumitrescu, M.D.Cazacu, Theoretische und experimentelle Betrachtungen ¨uberdie Str¨omung z¨aherFl¨ussigkeiten um eine Platte bei kleinen und mittleren Reynoldszahlen. Zeitschrift fur¨ Angewandte Mathematik und Mechanik, 1, 50, 1970, 257- 280. [10] M.D.Cazacu, La similitude hydrodynamique concernant le ph´enom`enedu coup de b´elier. Rev. Roum. de Sci. Techn. Mec.´ Appl., Rom. Acad. Publ. House, Tome 12, No. 2, 1967, Bucharest, 317- 327.
ROMAI J., 6, 2(2010), 83–91
FIXED POINT THEOREMS IN E - METRIC SPACES Mitrofan M. Choban, Laurent¸iu I. Calmut¸chi Department of Mathematics, Tiraspol State University, Chi¸sin˘au,Republic of Moldova [email protected] Abstract The aim of the present article is to give some general methods in the fixed point theory for mappings of general topological spaces. Using the notions of the multi-metric space and of E-metric space, we prove the analogous of several theorems about common fixed points of commutative semigroups of mappings.
Keywords: fixed point, m-scale, semifield, multi-metric space, E-metric space, pseudo-metric. 2000 MSC: 54H25, 54E15, 54H13, 12J17, 54E40.
1. INTRODUCTION In the sequel, any space is considered to be Tychonoff and non-empty. We use the terminology from [6, 7]. Let R be the space of real numbers, ω = {0, 1, 2, ...} and N = {1, 2, ...}. In [1, 2] the metric spaces over topological semifields were introduced. The gen- eral concept of the metrizability of spaces is contained in [14]. The notion of a topological semifield may be generalized in the following way. We say that E is a metric scale or, briefly, an m-scale if: 1. E is a topological algebra over the field of reals, R; 2. E is a commutative ring with the unit 1 , 0; 3. E is a vector lattice and 0 ≤ xy provided 0 ≤ x and 0 ≤ y; 4. For any neighborhood U of 0 in E there exists a neighborhood V of 0 in E such that if x ∈ V and 0 ≤ y ≤ x, then y ∈ U. 5. For any non-empty upper bounded set A of E there exists the suppremum ∨A. (The set A is upper bounded if there exists b ∈ E such that x ≤ b for any x ∈ A.) If the set A is not upper bounded, then we put ∨A = ∞. From the condition 4 it follows: 6. If 0 ≤ yn ≤ xn and limn→∞ xn = 0, then limn→∞yn = 0 too. From the condition 5 it follows: 7. If the non-empty set A ⊆ E is lower bounded, then there exists the infimum ∧A. Any topological semifield is an m-scale. A topological product of m-scales is an m-scale.
83 84 Mitrofan M. Choban, Laurent¸iu I. Calmut¸chi
Let E be an m-scale. Denote by E−1 = {x : x · y = 1 for some y ∈ E} the set of all invertible elements of E. By N(0, E) we denote some base of the space E at the point 0. n We consider that 0 ≤ x 1 if 0 ≤ x < 1, 1 − x is invertible and limn→∞ x = 0. We put E(+,1) = {x ∈ E : 0 ≤ x 1}. If t ∈ R and 0 ≤ t < 1, then t · 1 ∈ E(+,1). We identify t with t · 1 ∈ E for each t ∈ R. A mapping d : X × X −→ E is called a pseudo-metric over m-scale E or a pseudo- E-metric if it satisfies the following axioms: d(x, x) = 0, (∀)x ∈ E; d(x, y) = d(y, x), (∀)x, y ∈ E; d(x, y) ≤ d(x, z) + d(z, y)(∀)x, y, z ∈ E.
Every pseudo-E-metric is non-negative, i.e. d(x, y) ≥ 0 for all x, y ∈ X. An E-pseudo-metric d is called an E-metric if it satisfies the following axiom:
d(x, y) = 0 if and only if x = y. The ordered triple (X, d, E) is called a metric space over m-scale E or an E-metric space if d is an E-metric on X. Let d be a pseudo-E-metric on X. If A ⊆ B, B ⊆ X, x ∈ X, r ∈ E and r ≥ 0, then d(A, B) = in f {d(x, y): x ∈ A, y ∈ B} and B(x, d, r) = {y ∈ X : d(x, y) < r}. If U ∈ N(0, E), then we put B(x, d, U) = {y ∈ X : d(x, y) ∈ U} for any x ∈ X. The family {B(x, d, U): x ∈ X, U ∈ N(0, E)} is the base of the topology T(d) of the pseudo-E-metric space (X, d, E). If d is an E-metric, then the space (X, T(d)) is a Tychonoff (completely regular and Hausdorff) space. The pseudo-E-metric d generates on X the canonical equivalence relation: x ∼ y iff d(x, y) = 0. Let X/d be the quotient set with the canonical projection πd : X −→ X/ρ ¯ −1 −1 and the metric d(u, v) = ρ(πd (u), πd (v)).
Definition 1.1. Let A be a non-empty set and {Eα : α ∈ A} be a family of m-scales. A multi-Eα-metric space is a pair (X, P), where X is a set and P = {dα : X × X −→ Eα : α ∈ A} is a non-empty family of pseudo-Eα-metrics on X satisfying the condition: x = y if and only if d(x, y) = 0 for each d ∈ P.
Remark 1.1. Let A be a non-empty set and {Eα : α ∈ A} be a family of m-scales. Then E = Π{Eα : α ∈ A} is an m-scale. If (X, P) is a multi-Eα-metric space, then: 1. d(x, y) = (dα(x, y): α ∈) is an E-metric on X and T(d) is the topology of the multi-Eα-metric space X. 2. Let ρα(x, y) = {bβ ∈ Eβ : β ∈ A}, where bα = dα(x, y) and bβ = 0 for all β , α. Then Q = {ρα : α ∈ A} is a family of pseudo-E-metrics on X and T(d) is the topology of the multi-E-metric space (X, Q). Fixed point theorems in E - metric spaces 85
Thus we may consider that E = Eα for any α ∈ A in the Definition 1.1. In particular, any multi-E-metric space is an E0-metric space. If E = R, then the pseudo-E-metric is called a pseudo-metric and the pseudo-E- metric space is called a pseudo-metric space. Fix a multi-E-metric space (X, P). A subset V ⊆ X is called P-open if for any x ∈ V there exists a set U ∈ N(0, E) and a finite set A = A(x, U) ⊆ P such that B(x, A, U) = ∩{B(x, ρ, U): ρ ∈ A} ⊆ V. The family T(P) of all P-open subsets is a completely regular Hausdorff topology on X. If T is a completely regular Hausdorff topology on X, then T = T(P) for some family P of pseudo-E-metrics on X (see [6]). If d(x, y) = (ρ(x, y): ρP), then d is an EP-metric on X. Let X be a space and F be a non-empty set of mappings f : X −→ X. If Fix( f ) , ∅ for each f ∈ F, then X is called a fixed point space relative to F. By Fix(F) = {x ∈ X : f (x) = x for any f ∈ F} we denote the set of all common fixed points of the mappings f ∈ F. If ϕ : X −→ X is a mapping, then by Coin(F, ϕ) = {x ∈ X : f (x) = ϕ(x) for any f ∈ F} we denote the set of all coincidence points of the mappings f ∈ F and mapping ϕ. The excellent book [7] contains the fixed point theory for metric spaces with the important applications in distinct domains. Several results for general topological spaces with interesting applications are contained in the surveys [7, 12]. Distinct generalizations of the Banach fixed point principle were proposed in [3, 7, 12]. Dis- tinct classes of metrizable spaces were examined in [8, 9, 14]. In [16], the following general problem arose it was arisen the following general problem: to find topolog- ical analogies of the Banach fixed point principle. Some solutions to this problem were proposed in [5, 17]. This article is a continuation of the works [5, 16, 17, 11]. For each family γ of subsets of a set X and any x ∈ X the set S t(x, γ) = ∪{U ∈ γ : x ∈ U} is the star of the point x relative to γ. We put S t(A, γ) = ∪{S t(x, γ): x ∈ A}.
2. COMPLETE MULTI-METRIC SPACES Fix an m-scale E and an E-metric space (X, d, E). A set {xµ : µ ∈ M} is a net or a generalized sequence if M is a directed set. A point x ∈ X is a limit of the net {xµ : µ ∈ M} and we put limµ∈M xµ = x if for any U ∈ N(0, E) there exists λ ∈ M such that d(x, xµ) ∈ U for any µ ≥ λ. A net {xµ : µ ∈ M} is called fundamental if for any U ∈ N(0, E) there exists λ ∈ M such that d(xµ, xη) ∈ U for any µ, η ≥ λ. Any convergent net is fundamental. The limit of a fundamental sequence is unique (if the limit exists). The space (X, d, E) is called complete if any fundamental net is convergent (see [6, 1, 2]).
Example 2.1. Let A be an uncountable set. We put Iα = I = [0, 1] for any α ∈ A. A Let Y = I = Π{Iα : α ∈ A}, 0A = (0α : α ∈ A) ∈ Y and X = Y \{0A}. The space E = RA is a topological field and an m-scale. By construction, X ⊆ IA ⊆ E. We put d(x, y) = (|xα − yα| : α ∈ A) and dα(x, y) = |xα − yα| for any pair of points 86 Mitrofan M. Choban, Laurent¸iu I. Calmut¸chi
x = (xα : α ∈ A) ∈ Y and y = (yα : α ∈ A) ∈ Y. Let P = {dα : α ∈ A}. Then (X, d, E) is an E-metric space and (X, P) is a multi-metric space. Obviously, T(d) = T(P). The −n space Y is compact and the space X is pseudocompact. We put xn = (2α : α ∈ A). Then {xn : n ∈ N} is a fundamental sequence. We have limP xn = 0A. Thus the spaces (X, d, E) and (X, P) are not complete.
Remark 2.1. Let (X, d, E) be an E-metric space and the space E be first countable, i.e. it is metrizable in the usual sense. Then the space X is metrizable in the usual sense too.
3. CLOSURE OPERATORS ON PSEUDO-E-METRIC SPACES Fix an m-scale E, a topological space X and a pseudo-E-metric d on a space X. Denote by exp(X) the set of all subsets of X. A family γ ⊆ exp(X) is centered if ∩η , ∅ for any finite subfamily η ⊆ γ. If A is a non-empty subset of X, then diam(A) = ∨{d(x, y): x, y ∈ A}. Consider that diam(∅) = 0. A closure operator on the pseudo-E-metric space (X, d, E) is a mapping C : exp(X) −→ exp(X) satisfying the conditions: 1. if A ⊆ B, then C(A) ⊆ C(B); 2. A ⊆ C(A) and diam(C(A)) = diam(A); 3. C(C(A)) = C(A). Let C be a closure operator on X. The set A ⊆ X is C-closed if C(A) = A. The space X is called C-compact if ∩γ , ∅ for any centred family γ of non-empty C- closed sets. If the space X is compact and C(A) = clX A is the closure of the set A, then C is a closure operator on X and X is a C-compact space.
4. NIEMYTZKI FIXED POINT THEOREM FOR MULTI-METRIC SPACES The next theorem for compact metric spaces due to V. Niemytzki (see [15, 12]) and L.P.Belluce and W.A.Kirk [4].
Theorem 4.1. Let f, ϕ : X −→ X be mappings of a countably compact multi-metric space (X, P, R) into itself, ϕ(X) = X and d(ϕ(x), ϕ(y)) ≥ d(x, y) for all points x, y ∈ X and d ∈ P. Assume that d( f (x), f (y)) < d(x, y) for any pair of distinct points x, y ∈ X and any d ∈ P. Then: 1. The set Fix( f ) = {x ∈ X : f (x) = x} is non-empty and it is a singleton. 2. The set Coin( f, ϕ) = {x ∈ X : f (x) = ϕ(x)} is non-empty and it is a singleton.
Theorem 4.2. Let f, ϕ : X −→ X be mappings of a countably compact multi-metric space (X, P, R) into itself, ϕ(X) = X and d(ϕ(x), ϕ(y)) > d(x, y) for any pair of distinct Fixed point theorems in E - metric spaces 87 points x, y ∈ X and any d ∈ P. Assume that d( f (x), f (y)) ≤ d(x, y) for all points x, y ∈ X and d ∈ P. Then: 1. the set Fix(ϕ) = {x ∈ X : ϕ(x) = x} is non-empty and it is a singleton. 2. the set Coin( f, ϕ) = {x ∈ X : f (x) = ϕ(x)} is non-empty and it is a singleton.
Proof. The mapping ϕ is one-to-one and the mappings f and ϕ−1 are uniformly con- tinuous. Moreover, if g = ϕ−1 ◦ f , i.e. g(x) = ϕ−1( f (x)) for any x ∈ X, then d(g(x), g(y)) < d(x, y) for any pair of distinct points x, y ∈ X and any d ∈ P. Fix ρ ∈ P. Consider the projection πρ : X −→ X/ρ of (X, ρ) onto the metric space (X/ρ, ρ¯) and the mapping gρ : X/ρ −→ X/ρ, where g(πρ(x)) = πρ(gρ(x)) and ρ(x, y) = ρ¯(g(x), g(y)) for all x, y ∈ X. By construction, X/ρ is a compact metric space andρ ¯(g(x), g(y)) < ρ¯(x, y) for all x, y ∈ X/ρ. Suppose that for any n ∈ ω there exists xn, yn ∈ X such that −n ρ(g(xn), g(yn) ≥ (1 − 2 )ρ(x.y). Since X/ρ is a metrizable compact space, there exist an infinite subset L ⊆ ω, an infinite subset M ⊆ L and two points x, y ∈ X such that: −n - πρ(x) is the limit of the subsequence {πρ(xn): n ∈ L} and ρ(x, xn) < 2 for any n ∈ L; −n - πρ(y) is the limit of the subsequence {πρ(yn): n ∈ M} and ρ(y, yn) < 2 for any n ∈ M. There exists m ∈ ω such that 2−m+2 < ρ(x, y) − ρ(g(x), g(y)). By construction, ρ(g(x), g(y)) = lim ρ(g(xn), g(yn)) = lim ρ(xn, yn) = ρ(x, y), a contradiction. Thus there exists a number kρ < 1 such that ρ(g(x), g(y)) ≤ k · ρ(x, y) for all x, y ∈ X. Fix x = x0. We put xn = g(xn−1) for any n ∈ N. There exists b ∈ E such that b · (1 − k) = 1. Then d(xn, xn+1) = d(g(xn−1), g(xn)) ≤ k · d(xn−1, xn) for any n ∈ N. Hence n d(xn, xn+1) ≤ k · d(x0, x1) for any n ∈ N. Obviously, d(xn, xm) ≤ d(xn, xn+1 + n n+1 m−1 m ... + d(xm−1, xm) ≤ (k + k + ... + k + k ) · d(x0, x1). Therefore d(xn, xm) ≤ n n+m n (k − k ) · b · d(x0, x1) provided n, m ∈ N and n ≤ m. Since limn→∞k = 0, the sequence {xn : n ∈ ω} is fundamental in (X, ρ). Since the space X is countably compact and ρ is an arbitrary pseudo-metric, there exists a unique point b ∈ X such that b = lim xn. Then g(b) = lim g(xn) = lim xn+1 = b. Thus f (b) = ϕ(b). The assertion 2 is proved. The assertion 1 follows from the assertion 2. The proof is complete.
The assertion 1 for compact metric spaces was proved by V.Niemytzki (see [15, 12]) in 1936.
5. MAPPINGS OF E-METRIC SPACES WITH DIMINISHING ORBITAL DIAMETERS Fix an m-scale E. Let (X, d, E) be a pseudo-E-metric space and let X be endowed with the topology T(d) generated by the pseudo-E-metric d. Obviously, clA is a closure operator on X. 88 Mitrofan M. Choban, Laurent¸iu I. Calmut¸chi
Consider a commutative semigroup F of mappings f : X −→ X. We can assume that the identity mapping eX : X −→ X is an element of F. For any point x ∈ X the set F(x) = {(g(x)) : g ∈ F} is the F-orbit of x. Let r(x, F, d) = ∧{diam(F(g(x)) : g ∈ F}. Since F is a semigroup, the family F2(x) = {F(g(x)) : g ∈ F} is centered. The set L ⊆ X is F-invariant if g(L) ⊆ L for any g ∈ F. Let In(x, F, d) = ∩{H ⊆ X : x ∈ X, clH = H, H is F-invariant}. Obviously, F(x) ⊆ In(x, F, d). A semigroup F is said to have diminishing orbital diameters (d.o.d.) (see [4, 10, 13]) if for every point x ∈ X we have r(x, F) < diam(F(x)) provided diam(F(x)) > 0. Our definition is more general then that in [13]. The set VCoin(F) = {x ∈ X : d(x, f (x)) = 0 for any f ∈ F} is the set of all d-fixed points of the family F. If d is a E-metric, then VCoin(F) = Coin(F).
Theorem 5.1. Let F : X −→ X be a commutative semigroup of continuous mappings with diminishing orbital diameters of the pseudo-E-metric space (X, d, E) and the space (X, T(d)) is compact. Then VCoin(F) , ∅. Moreover, clF(x) ∩ VCoin(F) , ∅ for any point x ∈ X.
Proof. We can assume that the identity mapping eX ∈ F. Since the mappings F are continuous, the closure of an invariant set is invariant. Fix a point b ∈ X. If diam(F(b)) = 0, then b ∈ F(b) ⊆ VCoin(F). Assume that diam(F(b) > 0. Since F has d.o.d. and the space (X, T(d)) is compact, there exists a minimal invariant d-closed subset H of X such that diam(H) < ∞ and H ⊆ clF(b). We affirm that diam(H) = 0. Assume that diam(H) > 0. Fix c ∈ H. If diam(F(c)) = 0), then P = clF(c) ⊆ H, diam(P) = 0 and P is a d-closed invariant set, in contradiction with the minimality of the set H. Thus 0 < diam(F(c)) ≤ diam(H). Since F has d.o.d., there exists a point y ∈ F(c) such that diam(F(y)) < diam((F(c)). Then B = clF(y) ⊆ H, diam(B) < diam(H) and B is a d-closed invariant set, a contradiction with the minimality of the set H. Therefore diam(H) = 0. Then H ⊆ clF(x) ∩ VCoin(F).
Corollary 5.1. Let F : X −→ X be a commutative semigroup of mappings with di- minishing orbital diameters of the compact E-metric space (X, d, E). Then Coin(F) , ∅. Moreover, clF(x) ∩ Coin(F) , ∅ for any point x ∈ X.
The Corollary 4.2 for the compact metric spaces was proved in [4, 10, 13].
Remark 5.1. If E = R and F : X −→ X is a commutative semigroup of continuous mappings with diminishing orbital diameters of the pseudo-E-metric space (X, d, E), then d( f (x), f (y)) ≤ d(x, y) for all x, y ∈ X and f ∈ F. In particular, the mappings F are continuous. Fixed point theorems in E - metric spaces 89
6. MAPPINGS OF MULTI-E-METRIC SPACES WITH DIMINISHING ORBITAL DIAMETERS Fix an m-scale E. Let (X, P) be a multi-E-metric space and let X be endowed with the topology T(P) generated by the pseudo-E-metrics P. Obviously, clA is a closure operator on X. For any A ⊆ X and d ∈ P the element diamd(A) = ∧{d(x, y): x, y ∈ A} is the d-diameter of the set A. We put diam(A) = {diamd(A): d ∈ P}. The element b = (bd : d ∈ P) is called an E-P-number if bd ∈ {∞} ∪ E for any d ∈ P. If b = (bd : d ∈ P) and c = (cd : d ∈ P) are two E-P-numbers, then b < c if and only if bd ≤ cd for all d ∈ P and bρ < cρ for some ρ ∈ P. By construction, diam(A) is an E-P-number. Consider a commutative semigroup F of mappings f : X −→ X. We can assume that the identity mapping eX : X −→ X is an element of F. For any point x ∈ X the set F(x) = {(g(x): g ∈ F} is the F-orbit of x. Let r(x, F, d) = ∧{diamd(F(g(x)) : g ∈ F} and r(x, F) = (r(x, F, d): d ∈ P). By construction, r(x, F) is an E-P-number. Let In(x, F) = ∩{H ⊆ X : x ∈ X, clH = H, H is F-invariant}. Obviously, F(x) ⊆ In(x, F). A semigroup F is said to have diminishing orbital diameters (d.o.d.) if for every point x ∈ X we have r(x, F) < diam(F(x)) provided diam(F(x)) > 0.
Theorem 6.1. Let F : X −→ X be a commutative semigroup of continuous mappings with diminishing orbital diameters of the multi-E-metric space (X, P) and the space (X, T(P)) is compact. Then Coin(F) , ∅. Moreover, clF(x) ∩ Coin(F) , ∅ for any point x ∈ X.
Proof. The proof is similar to the proof of Therem 4.1. We can assume that the identity mapping eX ∈ F. Since the mappings F are continuous, the closure of an invariant set is invariant. Fix a point b ∈ X. If diam(F(b)) = 0, then b ∈ F(b) ⊆ Coin(F). Assume that diam(F(b) > 0. Since F has d.o.d. and the space (X, T(d)) is compact, there exists a minimal invariant closed subset H of X such that diam(H) < diam(F(b)) and H ⊆ clF(b). We affirm that diam(H) = 0. Assume that diam(H) > 0. Fix c ∈ H. If diam(F(c)) = 0), then P = clF(c) ⊆ H, diam(P) = 0 and P is a closed invariant set, a contradiction with the minimality of the set H. Thus 0 < diam(F(c)) ≤ diam(H). Since F has d.o.d., there exists a point y ∈ F(c) such that diam(F(y)) < diam((F(c)). Then B = clF(y) ⊆ H, diam(B) < diam(H) and B is a closed invariant set, in contradiction with the minimality of the set H. Therefore diam(H) = 0. Then H ⊆ clF(x) ∩ Coin(F).
By the virtue of the following example, the method of multi-metrics is more effec- tive than the method of metrics. 90 Mitrofan M. Choban, Laurent¸iu I. Calmut¸chi
Example 5.2. Let Xn = [0, 1] and ρn(x, y) = (n + 1) · |x − y| for all x, y ∈ Xn and n ∈ ω. We put X = Π{Xn : n ∈ ω} and d(x, y) = (n + 1) · |xn − yn| for all x = (xn : n ∈ ω) ∈ X, y = (yn : n ∈ ω) ∈ X and n ∈ ω. If P = {dn : n ∈ ω}, −1 then (X, P) is a metrizable multi-metric space. Let f (x) = (2 xn : n ∈ ω) and −n−1 g(x) = (2 xn : n ∈ ω) for any x = (xn : n ∈ ω) ∈ X. Denote by F the semigroup with identity generated by the mappings f and g. Since f g = g f , the semigroup F is ω commutative. If t ∈ R, then t¯ = (xn = t : n ∈ ω) ∈ R . If 0¯ < x = (xn : n ∈ ω) ∈ X, then diam(F(g(x)) < diam(F(x) and diam(F(g(x)) < 1.¯ Moreover, r(x, F) = 0 for any x ∈ X. Thus the semigroup F has diminishing orbital diameters. By construction, −1 diam(F(t¯)) = (2 ·(n+1)·t : n ∈ ω) for any t ∈ [0, 1]. If t > 0, then sup{diamdn (F(t¯): n ∈ ω} = ∞.
7. ASYMPTOTICAL PROPERTIES OF MAPPINGS Fix an m-scale E. Let E(−,1) = E(+,1) ∩ E−1. If ∈ E(−,1) and 0 < t < 1, then t · , t · 1 ∈ E(−,1). Consider an E-metric space (X, d, E) and a non-empty commutative semigroup F : X −→ X of continuous mappings of X into itself. A semigroup F : X −→ X is called asymptotically non-expansive (see [10, 13]) if for any x, y ∈ X there exists g ∈ F such that d( f g(x), f g(y)) ≤ d(x, y) for all f ∈ F. If A ⊆ X, then the set AF = {z ∈ X : there exists x ∈ A such that for every f ∈ F and ∈ E(−,1) there exists g ∈ F with d( f g(x), z) < } is called the F-closure of the set A in X.
Proposition 7.1. Let F : X −→ X be a non-empty commutative semigroup of con- tinuous asymptotically non-expansive mappings of an E-metric space (X, d, E) into itself. If z ∈ XF, then: 1. For all f ∈ F and ∈ E(−,1) there exists g ∈ F such that d( f g(z), z) < }. 2. For any f ∈ F the mapping f |F(z) is an isometrical embedding of F(z) into F(z). 3. diam(F(z)) = diam(F(g(z))) for all g ∈ F.
Proof. Is similar to the proof of Propositions 1 and 2 from [10].
Corollary 7.1. ([13], Proposition 1, for metric space). Let F : X −→ X be a non- empty commutative semigroup of continuous asymptotically non-expansive mappings of an E-metric space (X, d, E) into itself. If the semigroup F has diminishing orbital diameters and z ∈ XF, then F(z) = {z} and z ∈ Coin(F). References
[1] M. Ya. Antonovskii, V. G. Boltyanskii, T. A. Sarymsakov, An outline of the theory of topological semifields, Russian Math. Surveys, 21:4 (1966), 163-192. Fixed point theorems in E - metric spaces 91
[2] M. Ya. Antonovskii, V. G. Boltyanskii, T. A. Sarymsakov, Topological semifields and their ap- plications to general topology, American Math. Society. Translations, Series 2, 106, 1977. [3] V. Berinde, M. Choban, Remarks on some completeness conditions involved in several common fixed point theorems, Creative Mathematics and Informatics, 19:1(2010), 3-12. [4] L. P. Belluce, W. A. Kirk, Fixed point theorems for certain classes of nonexpansive mappings, Proceed. Amer. Math. Soc., 20 (1969), 141-146. [5] L. Calmut¸chi, M. Choban, On mappings with fixed points, Buletin S¸tiint¸ific. Universitatea din Pites¸ti, Matematica s¸i Informatica,˘ 3 (1999), 91-96. [6] R. Engelking, General Topology, PWN. Warszawa, 1977. [7] A. Granas, J. Dugundji, Fixed point theory, Springer, New York, 2003. [8] G. Gruenhage, Generalized Metric spaces, In. Handbook of Set-Theoretic Topology, K. Kunen and J. E. Vaughan, eds., Elsevier, Amsterdam, 1984. [9] G. Gruenhage, Metrizable spaces and generalizations, In: Recent Progress in General Topology II, M. Husekˇ and J. van Mill, eds., Elsevier, Amsterdam, 2002, 201-225. [10] R.D.Molmes, P.P.Narayanaswami, On asymptotically nonexpansive semigroups of mappings, Canadian Math. Bull., 13 (1970), 209-214. [11] K. Iseki, On Banach theorem on contractive mappings, Proceed. Japan Academy, 41 :2 (1965), 145-146. [12] A. A. Ivanov, Fixed points of mappings of metric spaces, Journal of Soviet Mathematics 12:1 (1979), 1-64. [13] Mo Tak Kiang, Semigroups with diminishing orbital diameters, Pacific J. Math. 41:1 (1972), 143-151. [14] S. I. Nedev, M.M.Choban, General conception of the metrizability of topological spaces, Annales of the Sophia University, Mathematics, 65 (1973), 111-165. [15] V. Niemytzki, The method of fixed points in analysis, Uspekhi Matem. Nauk, 1 (1936), 141-174. [16] J. Stepfans, S. Watson, W. Just, A topological fixed point theorem for compact Hausdorff spaces, York University, Preprint, 1991. [17] I. V. Yashchenko, Fixed points of mappings and convergence of the iterations of the dual map- pings, In: Obshchaya Topologia: Otobrajenia, Proizvedenia i Razmernosti Prostranstv, Moskva: Moskov. Gosud. Uni-t, 1995, 131-142.
ROMAI J., 6, 2(2010), 93–105
DISLOCATIONS AND DISCLINATIONS IN FINITE ELASTO-PLASTICITY Sanda Cleja-T¸igoiu Faculty of Mathematics and Computer Science, University of Bucharest, Romania [email protected] Abstract The paper deals with analytical description of the dislocation and disclination, which are lattice defects of the crystalline materials. The evolution equations for disclinations, having as source the screw dislocations are derived within the constitutive framework of second order plasticity developed by the author in the papers appeared in Int. J. Fracture (2007), (2010).
Keywords: dislocation, disclination, elastic and plastic distortions, micro and macro forces, balance equations, connections, evolution and constitutive equations. 2000 MSC: 74C99, 74A20.
1. INTRODUCTION The plastic deformability of metals, which are crystalline materials, is produced because of the existence of lattice defects inside the microstructure. The lattice de- fects, among which the dislocations, disclinations and point defects are mathemati- cally modeled by the differential geometry concepts, as torsion, curvature and metric property of certain connection, see Kroner¨ [10], [11], de Wit [7], but without any elasto-plastic constitutive equations. The elastic models for crystal defects can be found in Teodosiu [14]. We are not dealing with curved space but with curved geom- etry in flat space as stated de Wit [7]. 1. The nature of the geometry is determined by the linear connection Γ, fixed by its coefficients, the curvature tensor R, the Cartan torsion or torsion tensor S; 2. the metric tensor C, to measure the distance; 3. the non-metricity measure Q, of the connection relative to the metric tensor, i.e. in terms of Γ and C. Geometry for which R, S, Q are non-vanishing is non-metric, non-Riemannian. We restrict ourself to the case Q = 0, which means that the geometry is metric. If R = 0 the geometry is called flat. If S = 0 the geometry is called symmetric. If S = 0 the geometry is called Riemannian. If R = 0, S = 0 geometry is called Euclidean. In this paper the dislocations and disclinations are lattice defects of interest. The dislocations are characterized by the Cartan torsion, or by the non-zero curl of plastic connection, which means that the plastic distorsion can not be derived from a certain potential. The disclinations are characterized by a non-vanishing curvature R.
93 94 Sanda Cleja-T¸igoiu
The mathematical description of the continuously distributed dislocations is given by Noll [13], and the differential geometry support within the context of continuum mechanics can be found also in [12]. In this paper a peculiar mathematical problem is analyzed: find the disclinations, which are solutions of the appropriate evolution equation in such a way that the micro balance equation are satisfied, when the distribution of the dislocations is given. To give the mathematical description of the problem, we precise the general con- stitutive framework which is able to capture the dislocations and disclinations. We mention here the direction developed by Clayton et al. [3] within a micropolar elasto- plastic model, in order to emphasize the translational (dislocation) and rotational (disclinaton) defects. The behaviour of elasto-plastic body is described within the constitutive frame- work of second order plasticity, see Cleja-T¸igoiu [5], [6], based on the decomposition rule of the motion connection into the elastic and plastic second order deformations, see Cleja-T¸igoiu [4], and on the existence of configurations with torsion. The so called configuration with torsion is denoted by K, and it is described through a pair, composed by Fp, an invertible second order tensor which is called plastic distorsion (p) and Γ, a third order tensor field which is called plastic connection. (p) The pair (Fp, Γ) defines the second order plastic deformation with respect to the reference configuration of the body B. The actual configuration of the body is asso- ciated with t he motion function, which is defined at every time t by χ(·, t): B −→ E, E being the Euclidean (flat space) with a three dimensional vector space V. The sec- ond order elastic deformation, as a consequence of the decomposition rule, is defined with respect to the so called configuration with torsion K, which is time dependent. Two type of forces, macro and micro forces, are considered in the model. They are viewed as pairs of stress (a second order tensorial field) and stress momentum (a third order tensorial field), and they obey their own balance laws. The micro forces satisfy the viscoplastic type constitutive equations, in Kt. The evolution equations for Kt, which means the plastic distortion and plastic connection, have to be given. They have been derived to be compatible with an appropriate dissipation postulate. The energetic arguments, like a virtual power principle, macro and micro balance equations, and especially a dissipation postulate (see energy imbalance principle, for instance in Gurtin [9]) in order to prove thermomechanics restrictions, see Gurtin [9], Cleja-T¸igoiu [5]. In our adopted formalism the measure of the dislocations is characterized by the non-vanishing plastic curl (see Bilby [1], Noll [13]), while the disclinations has been related to a certain second order tensor Λ which enters the expression of the plas- tic connection and generates a non-zero curvature, apart from de Wit [7], where a measure of disclination is considered to be a second order curvature tensor. The following notations, definitions and relationships will be used in the paper: Dislocations and disclinations in finite elasto-plasticity 95
u · v, u × v, u ⊗ v denote scalar, cross and tensorial products of vectors; a ⊗ b and a ⊗ b ⊗ c are a second order tensor and a third order tensor defined by (a ⊗ b)u = a(b · u), (a ⊗ b ⊗ c)u = (a ⊗ b)(c · u), for all vectors u. For A ∈ Lin (Lin- the space of second order tensors), we introduce: the tensorial product A ⊗ a for a ∈ V, is a third order tensor, with the property (A ⊗ a)v = A(a · v), ∀v ∈ V. I is the identity tensor in Lin, AT denotes the transpose of A ∈ Lin, ∇A is the derivative (or the gradient) of the field A in a coordinate system {xa} ∂Ai j (with respect to the reference configuration), ∇A = ei ⊗ e j ⊗ ek. ∂xk Definition 1.1. The curl operator is defined for any smooth second order field, say A, through (curlA)(u × v) = (∇A)u)v − ((∇A)v)u, ∀ u, v, z ∈ V. (1) Lin(V, Lin) = {N : V −→ Lin, linear}− defines the space of all third order i j k tensors and it is given by N = Ni jki ⊗ i ⊗ i . T The scalar product of two second order tensor A, B is A · B := tr(AB ) = Ai jBi j, and the scalar product of third order tensors is given by N · M = Ni jk Mi jk, in a Carte- sian coordinate system. AB denotes the composition of A, B ∈ Lin. The compositions of A ∈ Lin and N, a third order tensor, defines the appropriate third order tensorial i j i p fields, via the formulae NAu = Ni jkAkpupi ⊗ i , and ANu = Ai jN jpkuki ⊗ i , which are written in a Cartesian basis, for any vector u. The affine connection is defined in a coordinate system by its coordinate represen- tation i m k Γ = Γmkei ⊗ e ⊗ e . (2) We introduce the third order tensor field Γ[F1, F2], which is generated by a third order field Γ together with the second order tensors F1, F2 through the formula
(Γ[F1, F2]u)v = (Γ(F1u)) F2v, ∀u, v ∈ V. (3)
For any Λ1, Λ2 ∈ Lin we define a third order tensor associated with them, denoted Λ1 × Λ2, by ((Λ1 × Λ2)u)v = (Λ1u) × (Λ2v), ∀u, v ∈ V. (4)
For A, a third order tensor, we define the vector field tr(2)A through the relationship written for all vectors (tr(2)A) · u = tr(Au). (5)
Three types of second order tensors, A B, A r B and A l B will be associated with any pair A, B of third order tensors, following the rules written for all L ∈ Lin
(A B) · L = A[I, L] · B = AiskLsnBink (A r B) · L = A · (LB) = Ai jkLinBn jk (6) (A l B) · L = A · (BL) = Ai jkBi jnLkn. 96 Sanda Cleja-T¸igoiu 2. GEOMETRIC RELATIONSHIPS Let F(X, t) = ∇χ(X, t) be the deformation gradient at time t, X ∈ B, and Γ = F−1∇F be the motion connection or the material connection. ∇χF is a gradient in the actual configuration, while the gradient in the configuration with torsion K, ∇KF, is calculated by p −1 ∇KF := (∇F)(F ) . (7) Ax.1 The decomposition of the second order deformation, (F, Γ), associated with (e) (p) e p the motion of the body B, into elastic, (F , Γ K), and plastic, (F , Γ), second order deformations is given by F = FeFp, (p) (e) (8) p −1 p p −1 Γ =Γ +(F ) Γ K [F , F ], Γ = F ∇F.
(p) Here the plastic connection with respect to the configuration with torsion, Γ K, is (p) related with the plastic connection,Γ, previously defined with respect to the reference configuration, by the following relationships
(p) (p) p p −1 p −1 (9) Γ K= −F Γ [(F ) , (F ) ]. The plastic metric tensor Cp and strain gradient C are defined with respect to the reference configuration, while the elastic metric tensor Ce is defined in configuration with torsion by Cp := (Fp)T Fp, Ce = (Fp)−T C(Fp)−1, C = FT F. (10) Definition 2.1. The Bilby’s type plastic connection is defined in a coordinate system through (p) A:= (Fp)−1∇Fp. (11) Remark 2.1. The Cartan torsion that belongs to the Bilby’s connection is given by (Su)v = (Fp)−1 ((∇Fp)u)v − ((∇Fp)u)v , while the curvature tensor is vanishing. Moreover, if the second order torsion tensor N is defined by N(u × v) = (Su)v, then it has the representation N = (Fp)−1curlFp. Consequently, we can say that the torsion tensor S is equivalent to curlFp. Let us introduce the expression for the plastic connection with respect to the ref- erence configuration built by Cleja-Tigoiu in [5], which has metric property with respect to Cp, and that allows a represented under the form
(p) (p) Γ=A +(Cp)−1(Λ × I), (12) Dislocations and disclinations in finite elasto-plasticity 97 where the third order tensor Λ × I, generated by the second order (covariant) tensor Λ is defined by (4), namely ((Λ × I)u)v = Λu × v, ∀ u, v ∈ V. (13) Λ is called the disclination tensor.
3. SCREW DISLOCATION First the skew dislocation will be defined in connection with the definition of the Burgers vector. The Burgers vector can be defined in terms of the plastic distortion p F , by considering a closed curve (circuit) C0 in the reference configuration of the body and A0 a surface with normal N bounded by C0 Z Z Z p p b = F dX = (curl(F ))NdA = αKnKdAK, (14) C0 A0 AK where αK is Noll’s dislocation density in [13] 1 α ≡ (curl(Fp))(Fp)T . (15) K detFp The expression of the Burgers vector can be aproximate by the formula p b ' curl(F )N area(A0). (16) In crystal plasticity the presence of the defects inside the crystals is measured by non- vanishing Burgers vector. The integral representation (14) shows that non-zero curl of plastic distortion, supposed to be continuum and non-zero in a certain material neighborhood, leads to a non-vanishing Burgers vector. Definition 3.1. We say that the plastic distortion Fp characterizes a screw dislocation if the generated Burgers vector through a circuit with the appropriate normal N is collinear with the normal, i.e. b k N, in contrast with the edge dislocation when b ⊥ N. p Let us introduce a Cartesian basis (e1, e2, e3) and a plastic distortion F which defines a screw dislocation that correspond to a Burgers vector directed to e3, given by p F = I + e3 ⊗ τ, τ ⊥ e3, with (17) 2 p p p p τ : D ⊂ R −→ V, τ := F31e1 + F32e2, J := det(F ) = 1. The mathematical description of the problem related to the screw dislocation within the finite elasticity is performed by Cermelli and Gurtin [2]. If we consider that τ = τ(x1, x2), the curl of the plastic distortion, see (1), is given by p p ∂F ∂F curl Fp = 31 − 32 e ⊗ e , (18) ∂x2 ∂x1 3 3 98 Sanda Cleja-T¸igoiu and Bilby’s type plastic connection
(p) ∂τβ A:= (Fp)−1∇Fp = e ⊗ ∇τ, ∇τ = e ⊗ e , β, δ ∈ {1, 2}. (19) 3 ∂xδ β δ The representation given in (18) justifies the name attributed to the plastic distor- tion introduced in (17), taking into account the expression for the Burgers vector, for a plane curve with the normal e3. (p) Let us remark its appropriate trace tr(2) A, defined by the formula (5)
(p) (p) ∂τ3 tr A ·u := tr(A u) = tr((e ⊗ ∇τ)u) = e ∇τ(u) = uk = 0, (20) (2) 3 3 ∂xk is zero, since τ3 : τ · e3 = 0, τ being a vector function orthogonal on e3. The plastic metric tensor Cp has the representation
p C = I + τ ⊗ e3 + e3 ⊗ τ + τ ⊗ τ. (21) Let us introduce the disclination tensor Λ represented in terms of Frank vector ω, like in Cleja-T¸igoiu [5] Λ := ηω ⊗ ζ, (22) ζ the tangent vector line for the disclination field, and the scalar valued function η have to be defined. It follows that ∇Λ := ω ⊗ ζ ⊗ ∇η, Λ˙ := η˙ω ⊗ ζ, (23) if we take constant value for ζ. Hypothesis. We assume that Frank and Burgers vectors are orthogonal, ω · b = 0, here b = e3. This hypothesis corresponds to the physical meaning attributed to these types of lattice defects. See for instance [3]. Remark 3.1. In Fig.1 the Burgers vector produced by the screw dislocation and the rotation produced by the disclinations have been plotted, following the comments from [3]. 4. MACRO AND MICRO FORCES Within the constitutive framework developed in [6], we consider that the free en- ergy in the reference configuration can be introduce in terms of the deformation fields listed below (p) (p) ψ = ψ(Ce, Fp, A, Λ, ∇Λ) ≡ ψev(C, Fp, A) + ψd(Λ, ∇Λ), (24) but in contrast with [6] we have no influence of the elastic connection. Concerning the free energy function we assume that the defect energy is given by κ ψd(∇Λ) = 2 β2∇Λ · ∇Λ, (25) 2 2 Dislocations and disclinations in finite elasto-plasticity 99
Fig. 1. Lattice defects: Screw dislocaton with Burgers vector b, disclination with Frank vector ω.
with β2 a length scale parameter. Ax.2 The macro forces are T− the non-symmetric Cauchy stress and macro stress momentum, µ, which is described by a third order tensor field, satisfy the macro balance equation div T + ρb = 0, −2Ta = div µ, (26)
Ta is the skew symmetric part of the stress tensor. The equations (26) are similar to those proposed in Fleck et al. [8], see also Cleja-Tigoiu [4], [6] We denote by ρ, ρ,˜ ρ0 the mass densities in the actual configuration, in the config- uration with torsion and in the reference one. When we pull back the macro stress µ to the configuration with torsion, the appro- priate expression for the macro stress momentum µK is derived
1 e T 1 e −T e −T e p −1 µK := (F ) µ[(F ) , (F ) ], where F = F(F ) . (27) ρ˜0 eρ
Ax.3 The micro stress denoted by Υλ and micro momentum respectively, denoted by µλ, which are associated with the disclinations satisfy in K their own balance equation (see [6]):
λ λ λ p λ p λ p −T λ Υ = divK µ + ρ˜B ⇐⇒ J Υ = div J µ (F ) + ρ˜B . (28)
Hereρ ˜Bλ is mass density of the couple body force, J p =| detFp | . Ax.4 The micro stress Υp and micro stress momentum µp are associated with the plastic mechanism, and they satisfy their own balance equation in the configuration 100 Sanda Cleja-T¸igoiu with torsion K : p p λ Υ = divK (µ − µK) + ρ˜B ⇐⇒ (29) p p p p p −T p J Υ = div J (µ − µK)(F ) + ρ˜B . Hereρ ˜Bp is appropriate mass density of the couple body force.The equation (29) could be found in a modified form in [5]. The constitutive restrictions imposed by the free energy imbalance, that have been obtained by Cleja-T¸igoiu in [6], are adapted to the proposed here model as it follows. The elastic type constitutive functions are derived from the free energy function, viewed like a potential for Cauchy stress 1 T = 2F(∂ ψ)FT , µ = ∂ ψ, (30) ρ C Γ but here µ = 0, as a consequence of the supposition made in (24) that the free energy function is not dependent on Γ. We postulate here an energetic (non-dissipative) constitutive equation for the micro stress momentum associated with plastic behaviour and disclinations, respectively
1 p µ0 = ∂(p)ψ ρ0 A (31) 1 λ d 2 µ0 = ∂∇Λψ (∇Λ) = κ2β2 ω ⊗ ζ ⊗ ∇η. ρ0 The last equality is a consequence of (25) together with (23). The expressions for the micro stress plastic momentum µp and the macro stress momentum µK both of them being considered with respect to K, pulled back to the reference configuration cab be expressed under the form
1 p p T 1 p p −T p −T µ0 := (F ) µ [(F ) , (F ) ], ρ0 eρ (32) 1 p T 1 p −T p −T µ0 := (F ) µK[(F ) , (F ) ], ρ0 eρ as well as the micro stress momentum associated with the disclination in K can be λ expressed in terms of µ0 by
1 λ p −T 1 λ p T p T µ := (F ) µ0[(F ) , (F ) ]. (33) ρ˜ eρ0 The viscoplastic type, dissipative evolution equations for plastic distorsion Fp and Λ have been postulated in [6] to be compatible with the appropriate dissipation inequal- ity. The viscoplastic evolution equation for plastic distortion, giving rise to the rate Dislocations and disclinations in finite elasto-plasticity 101 of plastic distortion with respect to the reference configuration is given by lp = −(Fp)−1F˙ p 1 p p T p (34) (Σ0 − Σ0 ) + (F ) ∂Fp ψ = ξ1 l . ρ0 Definition 4.1. The Mandel type stress measures in the reference configuration are associated with the appropriate stresses as it follows
1 p 1 p T p p −T 1 1 T −T Σ0 = (F ) Υ (F ) , Σ0 = F TF (35) ρ0 ρ˜ ρ0 ρ and 1 λ 1 p T λ p −T Σ0 = (F ) Υ (F ) . (36) ρ0 ρ˜ The viscoplastic type, dissipative evolution equations are postulated for disclina- tion Λ (p) 1 λ 1 λ Σ0 − ∂Λ ψ + A µ0 − ρ0 ρ0 (37) (p) (p) 1 λ 1 λ ˙ − µ0 r A − µ0 tr(2)(A) = ξ3 Λ. ρ0 ρ0 Ax.5 The scalar constitutive functions ξ1, ξ3 are defined in such a way to be compat- ible with the dissipation inequality p p ˙ ˙ ξ1 l · l + ξ3 Λ · Λ ≥ 0. (38) Note that the dissipation inequality is reduced to the expression written in the left hand side of (39) if the viscoplastic type equations (34) and (38) have been accepted. λ The Mandel type stress measure associated with the disclination mechanism, Σ0 , is λ λ related to the micro stress Υ via (37), while µ0 is expressed in (31). 5. DISCLINATIONS GENERATED BY A PLASTIC DISTORTION We suppose that the second order disclination tensor Λ is described in terms of Frank vector ω, with the scalar intensity function η and the disclination line ζ, say constant during the deformation process, have to be found. We are now able to solve the problem: Find the disclination tensor Λ, having the expression (22), to be solution of the evolution equation (38) with the micro stresses related by (37), and which is compatible with the micro balance equation associated with the disclination, (28)2. We take into account the physically motivated hypothesis that the Frank and Burg- ers vectors are (fixed) orthogonal, as well as that the plastic distorsion (17) character- izes the screw dislocation, then
ω · e3 = 0, τ · e3 = 0, ∇(τ)u · e3 = 0, ∀ u ∈ V. (39) 102 Sanda Cleja-T¸igoiu
First of all we eliminate the micro stress from (38), via (28) together with (37). By definitions and the hypotheses, from (33) together with (31)2 we get λ p −T p −T p T µ (F ) = ρ0(F ) ∂∇Λψ[I, (F ) ]. (40) As a direct consequence of (24) together with (25) and (23) we can prove the follow- ing relationship
1 λ p −T 2 µ (F ) = κ2β2 {ω ⊗ ζ ⊗ ∇η + (ζ · e3)(ω ⊗ ζ ⊗ ∇η)}. (41) ρ0 Consequently, we apply the divergence operator to (42) and then