(Microsoft Word - \344\310\321\307\32327- 40) Al-Khwarizmi Engineering Journal Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, December, (2020) P. P. 27- 40 Studying the Radial and Tangential Velocity Components of the Epithelization Healing Post Photorefractive Keratectomy Surgery of the Human Eye Nebras H. Ghaeb Department of Biomedical Engineering / Al-Khwarizimi College of Engineering / University of Baghdad Email: nebras@kecbu.uobaghdad.edu.iq (Received 23 January 2019; accepted 20 September 2020) https://doi.org/10.22153/kej.2020.09.003 Abstract Photorefractive keratectomy (PRK) is the refractive technique that began with a physical scraping of the epithelial layer of cornea subsequent by laser treatment. Post this procedure to about 48 hours the removed epithelial layer regenerated to protect the eye again. The regeneration process (called re-epithelization) started from the limbus of the cornea toward the central part of it. The re-epithelization mechanism consists of a change in cell density (mitosis) and cell concentration (migration) with a velocity in two directions: radial and tangential. In the present study, an estimation for both radial (responsible for the overlapped layers toward the outward direction of the cornea) and tangential components (contour shape wave from limbus to the center) has been done for the first time, not like the previous studies that always estimate the velocity values of the re-epithelization only. Results showed that the trend shape of both components agrees with the kinematic behaviour of the mitosis and migration, where the maximum cell density fluctuated toward the central part in exponential decay shape. For a healing diameter of 2mm, the maximum redial velocity was 16.85 µ m/h, while the maximum tangential velocity was 55.48 µ m/h. These two components give a speed of re- epithelization of 58 µ m/h which agrees with the biological and practical healing speed measured of 60 µ m/h. Estimating these two components will open the way to understand the relationship between the total epithelial layer required and the total healing time to control the medication period for the patient post-surgery. Keywords: Re-epithelization, reaction – diffusion equation, PRK, speed of healing. 1. Introduction Photorefractive Keratectomy (PRK) is the correction of the visual disturbance in refraction errors for nearsightedness, farsightedness, and astigmatism using some certain laser type. It was started earlier in the eightieth of the last century and approved to be used finally by the Food and Drug Administration (FDA) in 1996 [1]. Medically this type of surgical procedure has the following steps: physical removing of the first layer of the cornea called the epithelial as it is a protective and regenerative layer, followed by applying the treatment laser that related in its energy to the refractive error shape and value, and finally protects the cornea with bandage contact lens until the epithelia build back and protect the cornea again. The epithelial building back process is called the reepithelization, which started from the limbus toward the central part of the cornea of the human eye. This process has been studied in many research works, to estimate its speed, and the encouragement factors affecting it, and Sherratt and Murray were pioneers in this field. Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 28 In 1990 Sherratt and Murray [2] begun their work in mathematical modeling of the healing of an epidermal wound. The work presented the epidermal cell density, as a function of the wound radius of curvature and time, all which examined in the reaction–diffusion equation (RDE), which state: In 1991 Sherratt and Murray [3] modified the simple RDE to a coupled equation, by including the concentration of mitosis. The differential equations are for cell density, and rate of change in chemical concentration of Epidermal Growth Factor (EGF). The mathematical model of wound healing becomes: Rate of variation in cell density, n = Migration of cell + Mitotic generation – Natural loss …(2a) Rate of variation of chemical concentration, c = Diffusion of c + Production of c – Chemical decay ...(2b) And in differential form: �� �� = �∇�� + �� . �. �2 − � ��� − �� …(3a) �� �� = �� ∇�� + ��� − �� …(3b) Where s(c), the logistic regression is a function of chemical concentration, k, Dc and λ are positive constants. ��� is a biological function. Equations (3) have been solved for different clinical applications. In ophthalmology, the corneal epithelial mapping attracted in the last five years the researchers, especially, after the presenting of the ocular coherence technology (OCT) system, which has the ability to measure and show the epithelial mapping [4]. Dale et al, [5] rewrite equations (3) considering the corneal epithelialization and solved them to find the speed of the process. The effect of concentration and density with the presence of EGF as a driving factor was considered. They proposed an analytical approximated solution verified by numerical scheme calculations. Their findings direct us to the truth that the EGF improves the healing rate or the speed of the traveling epithelial wave. The modified differential equations (based on equations (3) and customized for the epithelial corneal region) are: �� �� = ∇. ����� ∇� + �� � �� − � ��� − �� (4a) �� �� = �� ∇�� + ��� − ℎ�� � − �� (4b) Where Dn(c): is the cell diffusion coefficient, h(c): cellular degradation, �, k, Dc and λ are positive constants. Equations (4) have been utilized as a starting point for the next research studies that deal with the re-epithelization of the corneal layer (see Figure (1)). Although these studies showed many aspects of interest but still no mentioned for the cell re- epithelization velocity directions. These directions (radial and tangential) consist of the move of layers from limbus to the center in contour shape and toward the outward of the cornea to fully protect it again. Estimating the velocity of each component alone needs more mathematical manipulation to predict the required time of full protection and manage the time of medication used during this period. Rate of increase of cell density, n = Cell migration + Mitotic generation …(1) Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 29 Fig. 1. Research studies regarding the Re-epithelization of cornea. Four imperative disciplines (as shown above) were reflecting the interest of the researchers in studying the corneal re-epithelization process. These disciplines are: a. Mathematical Modelling. The modeling here started by normalizing equations (4) with its boundary and initial conditions. Then found a semi exact solution after considering sets of assumptions and simplifications [5, 6 and 7]. The final resulted manipulation may require a computer numerical scheme [5 and 11]. Some solutions utilize wave theory to solve equations (4) [9, 10, 12, and 14]. b. Biological Studies. In such studies, the biological behaviour of the cell during the reepithelization has been considered. How to control the speed through certain corticosteroids and the expected time period for the healing process [8, 11, 13, 26, and 41]. Moreover, study the cell kinetics during the whole healing process and compare between simulation program and real pathogenesis case studies [12]. The possibility of including few markers and measure the triggering signal from the stem cell to realize the optimum healing shape and rate [15, 21]. c. Refractive Outcomes. For the patients with low refraction error, the PRK will be the treatment procedure (a safe surgery procedure). In such cases, the re-epithelization process will gain back the epithelial layer of about 50 µ m which represents an extra power for the patient vision. During the healing period, the total thickness fluctuation and the visual outcomes also [13-19, 24, 26, 33-35, 38-40]. Within the last few years after the invention of the epithelial mapping which described by the OCT system many ophthalmologists begun to use it as an indication for the corneal abnormality such as dryness and keratoconus [23, 24, 27, 29, 31, 32-36, 41, 45, and 46]. d. Pharmaceutical Study. Study the normal amount of the EGF during the healing process and the possibility of utilizing some sort of pharmaceutical components to improve the supplied amount to the stem cell [5, 6, 11, and 24]. Possibility of studying the α and γ – EGF and their converse role during the re- epithelization, with the impact of including antibodies [11, 15, 24 and 26]. Assessment of effective gene and recognize the possibility of gene control to control the final healing shape and rate [43]. 2. Mathematical Model The differential equations (4) have been utilized here, applying the non-dimensionalize model recommended by Dale [5], and using the following assumptions: a. The thickness of the epithelium is littler than the required length of healing, so the model will be two-dimensional form. b. By considering linear geometry of the wound healing or strip band form with long spatial domain length the equation will be solved in the semi-infinite domain of −∞ ≤ � < Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 30 !" (Treatment Zone), and the solution will be one dimensional toward the x – direction. c. The original set of the axis is the wound boarder. d. The cellular diffusion coefficient changes linearly with EGF concentration value (c). The new differential form will be: �� �� = � �# $�%� + & �� �#' + �%(� + &( ��2 − � − � (5a) �� �� = �� �)� �#) + *�2 − 5� − ,�� ��̂.� − �� (5b) Where the dimensionless parameters have the following values: α = 0.01, β = 0.1, α1 = 0.9, β1 = 0.1, Dc = 25, σ = 4000, µ = 1.37x104, �̂ = 3.02, δ = 110 and the wound length or the TZ = 8 mm [4, 5, 6 and 7]. The initial and boundary conditions are biologically relevant to the surgical procedure or the PRK process. Figure (2) shows the center lines, wound center, epithelium boundary and the original point and axis. Fig. 2. Important notification on the human eye sketch. The boundary and initial conditions for the solution will be described graphically as shown in Figure (3): Fig. 3. Initial and Boundary Conditions selected. BEpi = Boundary of Epithelium, WC = Wound Center. Horizontal Center Line Vertical Center Line Wound Center E p it h e li u m B o u n d a ry Origin Treatment Zone (TZ) BEpi WC x O –∞ on(x,0) = n oc(x,0) = c −∞ ≤ � < 0 n(x,0) = 0 c(x,0) = 0 0 ≤ � < 0 TZ/2 = L Initial Conditions BEpi WC x O –∞ on(x,t) = n oc(x,t) = c −∞ ≤ � < 0 (x,L) = 0xn (x,L) = 0xc � = 0 TZ/2 = L Boundary Conditions Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 31 3. Proposed Solution In the present work a solution for the final equations (5) has been suggested using the Finite Difference (FD) implicit method Crank Nicolson Scheme (CNS). The CNS has the following definitions for the independent variable u, first- order in time, second-order in direction (such as x) [16]: 1234512 ∆� = ( � 789 :.( ��, <, =, �1�# , �)1 �#) � + 89: ��, <, =, �1�# , �)1 �#) �> …(6) Each part of the differential equations (5) can be define according to the followings: = = (� ?=9 :.( + =9: @ …(7a) �1 �� = 1234512 ∆� (1 st order in time for the independent variable u) … (7b) �1 �# = ( A∆# ?=9.( :.( − =95(:.( + =9.(: − =95(: @ …(7c) �)1 �#) = ( ��∆# ) �?=9.( :.( − 2=9:.( + =95(:.(@ + ?=9.(: − 2=9: + =95(: @� …(7d) The rearrangements and manipulation of the two differential equations, boundary and initial conditions will be done independently and then collect the final results as shown below: a. One Dimensional Reaction Equation From Equation (5a), rearrange the equation to be: �� �� = �%� + & �)� �#) + % ���� �#) + �%(� + &( ��2 −� − � …(8) Equation (8) will be in CNS form: �2345�2 ∆� = ( ��∆# ) B%?�9 :.( + �9: @ + &CB?�9.(:.( + �9.(: @ − 2?�9:.( + �9: @ + ?�95(:.( − �95(: @C +D (E�∆# ) B?�9.( :.( − �95(:.( + �9.(: − �95(: @?�9.(:.( − �95(:.( + �9.(: − �95(: @C + B�9:.( − �9: C 7D4� ?�9 :.( + �9: @ + &(> B1 − ?�9.(:.( + �9: @C …(9) b. One Dimensional Diffusion Equation: From Equation (5b), the CNS form is: � 2345� 2 ∆� = GH ��∆# ) B?�9.( :.( + �9:.(@ − 2?�9:.( + �9: @ + ?�95(:.( + �95(: @C + 7* �2 − 5?�9:.( + �9: @�> − ,��I234.�I2���I234.�I2� $�̂5��I234.�I2�' − *?�9:.( + �9: @..(10) Initial Conditions: �9� + �9( = �9� + �9( = 0 �JK 0 ≤ � < 0 …(11a) �9L��, 0 = �M N�O �9L��, 0 = �M J<ℎPK <ℎN� 0 ≤ � < 0 …(11b) Boundary Conditions: �9L��, 0 = �M N�O �9L��, 0 = �M J<ℎPK <ℎN� 0 ≤ � < 0 …(12a) For �# �0, < = �M N�O �# �0, < = �M �9.(:.( + �9.(: = �9Q.( + �9: N�O �9.(:.( + �9.(: = �9Q.( + �9: …(12b) c. Speed of healing Most of the studies regarding the speed of healing (U) started their modeling from the Fisher reaction-diffusion equation [2, 5, 10 and 20]. The traveling healing wave solution uses the analogy of the relative coordinate system to predict the final expected form, which may be expressed as [5]: R = �S.GH T&�� 72�� � U.V � � − �� − �&> …(13) Equation (13) has sets of constants that leads to a specific value of the healing speed of about 60 µ m/h. The modification that has been added here to equation (13) is to suggest the shape of a solution that may analogy the trend of the solution in both of the cell density and concentration. This suggested analogy predicted from the mathematical interpolation of equations (9 and 10). The new modified shape of the healing speed (U) will be: RWMX = �S.GH T&�� 72�� � U.V � � − �� − �&>. P ����.��� � sin ��MM − �MM …(14) Where noo and coo are the cell density and EGF concentration calculated from the above equations (9) and (10), as mentioned above for the state of constant speed of healing. d. Velocity components (radial and tangential): The flow velocity of the re-epithelization process is measured and estimated to be in a micro scale [2, 3 and 4]. The Reynold number, for this reason, is less than unity, and the flow could be considered as a creeping flow [43]. The description of the final equations for the two components velocity is: =\ = −R�J ] �1 − ^_�\ + _` �\`� …(15a) =a = −R b�] �1 − ^_A\ + _` A\`� …(15b) Where the R and r, are the radius of curvature of the human cornea and new build-up layer thickness Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 32 or the epithelial during the healing (see Figure (4) for more details). Fig. 4. Radial and tangential direction with the start and end of epithelization. The final algorithm for the proposed solution discussed above from point (a) to point (f), can be summarized in the flowchart as shown in Figure (5). Fig. 5. Proposed solution, equations here are discussed above from point (a) to point (f). Start Define the constant of solution: α = 0.01, β = 0.1, α1 = 0.9, β1 = 0.1, Dc = 25, σ = 4000, µ = 1.37x104, �̂ = 3.02, δ = 110 and the wound length or the TZ = 8 mm. Define the initial values: �9M ��, J , �9M ��, J Define the Boundary values: �9: �0, < , �9: �0, < , �9: �0, < , �9: �0, < , Solve Eq.(9) to find �9:.(��, < Solve Eq.(10) to find �9:.(��, < Use Eq.(14) to find the modified Healing Speed: Umod Use Eq.(15.a) to find the radial component Healing Speed: Ur Use Eq.(15.a) find the tangent component Healing Speed: Ur Change the time increment from 4h to 40h. End Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 33 4. Results The numerical solution of reaction equation (9) and diffusion equation (10), with their initial and boundary conditions (11 and 12) have been done using the open-source scientific programming language Octave (version 4.4.1). The solution of the above equations has been done after assuming the wound speed to be constant with a value of 60 µ m/h [5]. In this step, both the cell density (n) and the variation in concentration (c) are evaluated from the limbus toward the center of the cornea. In the next step, an equation for the speed of the healing will be utilized to investigate the range and the shape of the wave speed of the healing process. This time the speed of healing will be examined with respect to the change of both n and c, which will decrease the total time period of the re-epithelization process (time of healing). After all the two components were assessed also based on the overall speed of healing. Cell density and Concentration The change in cell density (n(x, t)) is a sign for cell migration, mitotic activity, and natural loss. This is what concerned in equation (2) and related specifically to the wound re-epithelization of the corneal tissue. Figure (6), shows the variation of cell density with regard to the direction of healing and time for the time period from 4 – 40 hours with ten steps changes. Fig. 6. Cell density for 10 steps of time with 4 hours per each, and through the length from the limbus to the center of healing (center of the eye). The change in chemical concentration of EGF which represents the main drive for the proliferation for epithelial cell, Figure (7) shows the change in concentration c(x, t) for 40 hours. Fig. 7. Cell concentration for EGF for 10 steps of time with 4 hours per each, and through the length from the limbus to the center of the eye. Speed of healing Figure (8) shows the healing speed variation with respect to the EGF concentration with a time of healing variation from 4 – 40 hours. Fig. 8. Speed of healing to the EGF concentration for 10 steps of time. Figure (9) shows the comparison between the first 4 hours with respect to the 12, 20, 28 and 36 hours. 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 C e ll D e n si ty n (x ,t ) X - Direction 04 hours 08 hours 12 hours 16 hours 20 hours 24 hours 28 huors 32 hours 36 hours 40 hours 0 10 20 30 40 50 60 70 80 90 100 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 E G F C O N C E N T R A T IO N C (X ,T ) X-DIRECTION 04 hours 08 hours 12 hours 16 hours 20 hours 28 hours 32 hours 36 hours 40 hours 0 20 40 60 80 100 120 140 0 200 400 S p e e d µ m /h EGF Concentration ng/ml 04 hours 08 hours 12 hours 16 hours 20 hours 24 hours 28 huors 32 hours 36 hours 40 hours Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 34 Fig. 9. Speed of healing with respect to the cell concentration of EGF for 4 steps: (a) 12 h, (b) 20h, (c) 28h and (d) 36h compared with 4 hours. The radial component of the velocity described in equation (15a) shown in figure (10) below. Fig. 10. Radial speed of healing w.r.t the EGF concentration for 10 steps of time. Figure (11) shows the comparison between the first 4 hours with respect to the 12, 20, 28 and 36 hours. 0 5 10 15 20 25 30 35 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 12 hours 0 5 10 15 20 25 30 0 100 200 300 400 500 S p e e d µ m /h EGF Concentration ng/ml 04 hours 08 hours 12 hours 16 hours 20 hours 24 hours 28 hours 32 hours 36 hours 40 hours (a) 0 20 40 60 80 100 120 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 20 hours (b) 0 20 40 60 80 100 120 140 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 36 hours (d) 0 20 40 60 80 100 120 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 28 hours (c) Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 35 0 1 2 3 4 5 6 7 8 9 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 20 hours Fig. 11. Radial speed of healing with respect the cell concentration of EGF for 4 steps: (a) 12h, (b) 20h, (c) 28h and (d) 36h compared with 4 hours. The tangential component of the velocity described in equation (15b) shown in figure (12) below. Fig. 12. Tangential speed of healing to the EGF concentration for 10 steps of time. Figure (13) shows the comparison between the first 4 hours to the 12, 20, 28 and 36 h. 0 1 2 3 4 5 6 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 12 hours 0 20 40 60 80 0 100 200 300 400 500 S p e e d µ m /h EGF Concentration ng/ml 04 hours 08 hours 12 hours 16 hours 20 hours 24 hours 28 hours 32 hours 36 hours 40 hours (b) (a) 0 0.5 1 1.5 2 2.5 3 3.5 4 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 28 hours 0 2 4 6 8 10 12 14 16 18 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 36 hours (d) (c) Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 36 Fig. 13. Tangential speed of healing with respect to the cell concentration of EGF for 4 steps: (a) 12h, (b) 20h, (c) 28h and (d) 36h compared with 4 hours. 5. Discussion Figure (6) showed the fluctuation in the cell density amount for the direction from the limbus to the center of the eye, for different time of healing started from 4h postoperatively to 40h. The behaviour here agrees with the wave theory and also satisfy the clinical evidence that the stem cells will kinetically increase the mitotic generation activity, which increases the number of cells ready for migration at the first instant after surgery. This motivated kinetic activity will be decreased exponentially reaching the new cells to the center of the cornea. For each instant of time (4h) a wave started from limbus toward the center with fluctuation period and depression rate different from one to another. Figure (7) the cell concentration behaviour started the first 4h with about flat change building up to increase the concentration near the center of the eye within time. After 40h the concentration level will reach the maximum required protection [5, 7, and 11]. It is evident from Figure (9) that the increase in the EGF concentration level will increase the healing speed. Biologically this is often due to an increase in the diffusion of the produced cell. This increase would happen after the wave of healing reaches close to the center of the corneal surface. The predicted speed of healing would alter from 14.690, 32.649, 97.208, 96.658, and 128.150 µ m/h for the 4, 12, 20, 28, 36 hours’ time of healing. The average of the ten steps maximum speed of healing is 72.872 µ m/h, whereas the average of the means of the ten steps of times is 62.516 µ m/h witch agrees with the biological and practical findings [2, 5, 10 and 20]. This leads to a practical sense that the speed of healing or diffusion rate both are compensated with respect to the position and time to recover the effect reepithelization of the surface of the cornea [17, 23, 29, and 30]. 0 2 4 6 8 10 12 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 12 hours 0 2 4 6 8 10 12 14 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 20 hours (a) (b) 0 10 20 30 40 50 60 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 28 hours 0 10 20 30 40 50 60 0 100 200 300 400 500 600 S p e e d µ m /h EGF Concentration ng/ml After 04 hours After 36 hours (c) (d) Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 37 6. Conclusions The overall healing process of epithelial post- PRK surgery has many biological and chemical effective parameters. These parameters all together affect the cells mitotic and migration, that physically represented by the speed of re- epithelization. To understand the mechanism of healing two components of velocity of re- epithelization have been estimated and compared with the biological finding. During the 40 hours of the study, the maximum radial velocity reaches about 16.85 µ m/h in the healing diameter of 2mm, whereas the maximum tangential velocity was about 55.48 µ m/h in the same healing diameter. These two components give a speed of re-epithelization of 58 µ m/h which agrees with the biological and practical healing speed of 60 µ m/h. Furthermore, we need more mathematical manipulation and prediction tools (such as artificial intelligence [47] and safety criteria [48] to check the predicted values) that finally we could estimate the required healing process per each patient after considering the other refractive error values. 7. References [1] J.L. Alio, F.A. Soria, A. Abbouda and P. Peña- García, "Fifteen years' follow-up of PRK up to 10 D of myopia: outcomes and analysis of the refractive regression," British J. Ophth. Vol. 100, iss. 5, pp. 626-32, 2016. [2] J.A. Sherratt and J.D. Murray, "Models of Epidermal Wound Healing," Proceedings of the Royal Society of London B: Biological Sciences, pp. 29-36, 1990. [3] J.A Sherratt and J.D. Murray, "Mathematical analysis of a basic model for epidermal wound healing," J. of Mathematical Biology, Vol. 29 iss. 5, pp. 389–404, 1991. [4] C.H. Karabatsas, "Studying corneal epithelium in vivo: A new tool in clinical practice," Advances in Ophth. & visual system, vol. 4, iss. 2, pp 3–5, 2016. [5] P.D. Dale, P.K. Maini and J.A. Sherratt, "Mathematical modeling of corneal epithelial wound healing," Mathematical Biosciences, Vol. 124, iss. 2, pp. 127–147, 1994. [6] P.D. Dale, J.A. Sherratt and P.K. Maini, "The Speed of Corneal Epithelial Wound Healing," Applied Mathematical Letters, Vol. 7, iss. 2, pp. 11–14, 1994. [7] P.D. Dale, J.A. Sherratt and P.K. Maini, "Corneal Epithelial Wound Healing," J. of Biological Systems, Vol. 3, iss. 4, pp. 957– 965, 1995. [8] S.W. Chang, F.R. Hu and P.K. Hou, "Corneal epithelial recovery following photorefractive keratectomy," The British J. of Ophth., Vol. 80, iss 7, pp. 663–668, 1996. [9] H. Sheardown and Y. Cheng, "Mechanisms of Corneal Epithelial Wound Healing," Chemical Engineering Science, Vol. 51, iss.19, pp. 4517–4529, 1996. [10] E.A. Gaffney, P.K. Maini, J.A. Sherratt, P.D. Dale, "Wound Healing in the Corneal Epithelial: Biologic Mechanisms and Mathematical Models," J. of Theor. Medicine, Vol. 1, iss. 1, pp. 13–23, 1996. [11] V.V. Kourenkov, O.N. Mytiagina and A.G. Pavluk, "Stimulating Re-epithelialization after Photorefractive Keratectomy," J. of Refr. Surg., Vol. iss. 15, pp. 234–237, 1999. [12] E.A. Gaffney, P.K. Maini, J.A. Sherratt and S. Tufu, The Mathematical Modelling of Cell Kinetics in Corneal Epithelial Wound Healing," J. of Theoretical Biology, Vol. 197, iss. 1, pp. 15–40, 1999. [13] R.H. Silverman, F.L. Lizzi, B.G. Ursea, L. Cozzarelli, J.A. Ketterling, C.X. Deng and J. Coleman, "Safety Levels for Exposure of Cornea and Lens to Very High-Frequency Ultrasound," J. Ultrasound Med, Vol. 20, pp. 979–986, 2001. [14] S. Serrao and M. Lombardo, "Corneal epithelial healing after photorefractive keratectomy: Analytical study," J. of Cataract and Refractive Surgery, Vol. 31, iss. 5, pp. 930–937, 2005. [15] A. Pajoohesh-Ganji and M.A. Stepp, "In search of markers for the stem cells of the corneal epithelium," Biology of the Cell / under the Auspices of the European Cell Biology Organization, Vol. 97, iss. 4, pp. 265– 76, 2005. [16] T. Callaghan, E. Khain, L.M. Sander and R.M. Ziff, "A stochastic model for wound healing," J. of Statistical Physics, Vol. 122, iss. 5, pp. 909–924, 2006. [17] D.Z. Reinstein, T.J. Archer, M. Gobbe, R.H. Silverman and D.J. Coleman, "Epithelial thickness in the normal cornea: three- dimensional display with very high frequency ultrasound," J. of Refractive Surgery, Vol. 24, iss. 6, pp. 571–581, 2008. [18] D.Z. Reinstein, S. Srivannaboon, M. Gobbe, T.J. Archer, R.H. Silverman, H. Sutton and D.J. Coleman, "Epithelial thickness profile changes induced by myopic LASIK as measured by Artemis very high-frequency Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 38 digital ultrasound," J. of Refractive Surgery, Vol. 25, iss. 5, pp. 444–450, 2009. [19] A. Ivarsen, W. Fledelius, and J.O. Hjortdal, "Three-year changes in epithelial and stromal thickness after PRK or LASIK for high myopia," Investigative Ophth. and Visual Science, Vol. 50, iss. 5, pp. 2061–2066, 2009. [20] F. Posta and T. Chou, "A mathematical model of intercellular signaling during epithelial wound healing," J. of Theoretical Biology, Vol. 266, iss. 1, pp. 70–78, 2010. [21] D.Z. Reinstein, M. Gobbe, T.J. Archer, R. H. Silverman and J. Coleman, "Epithelial, Stromal, and Total Corneal Thickness in Keratoconus: Three-dimensional Display with Artemis Very-high Frequency Digital Ultrasound, "J. of Refr. Surgery, Vol. 26, iss. 4, pp. 259–271, 2010. [22] Y. Li, O. Tan, R. Brass, J.L. Weiss, D. Huang, "Corneal Epithelial Thickness Mapping by Fourier Domain Optical Coherence Tomography in Normal and KCN Eyes," Ophth., Vol. 119, iss. 12, pp. 2425–2433, 2012. [23] A.J. Kanellopoulos, I.M. Aslanides and G. Asimellis, "Correlation between epithelial thickness in normal corneas, untreated ectatic corneas, and ectatic corneas previously treated with CXL; is overall epithelial thickness a very early ectasia prognostic factor?" Clinical Ophth., Vol. 6, iss. 1, pp. 789–800, 2012. [24] C. Du, J. Wang, L. Cui, S.A. Meixiao and Y. Yuan, "Vertical and horizontal corneal epithelial thickness profiles determined by ultra-high resolution OCT," Cornea, Vol. 1, iss. 9, pp. 1036–1043, 2012. [25] B.R. Fonslow, B.D. Stein, K.J. Webb, T. Xu, J. Choi, S. Kyu and J.R, Iii, "Wound Healing After Keratorefractive Surgery: Review of Biological and Optical Considerations," Cornea, Vol. 31, iss. 1, pp. 9–19, 2012. [26] A.J. Kanellopoulos and G. Asimellis, "Anterior segment OCT: Assisted topographic corneal epithelial thickness distribution imaging of a keratoconus patient," Case Reports in Ophth., Vol. 18, iss. 1, pp. 74-78, 2013. [27] A. Tao, Y. Shao, H. Jiang, Y. Ye, F. Lu, M. Shen and J. Wang, "Entire thickness profiles of the epithelium and contact lens in vivo imaged with high speed and high resolution OCT," Eye Contact Lens, Vol. 39, iss. 5, pp. 1–14, 2013. [28] A.J. Kanellopoulos, G. Asimellis, "In vivo 3D corneal epithelial thickness mapping as an indicator of dry eye: Preliminary clinical assessment," American J. of Ophth., Vol. 157 iss. 5, pp. 1116–1117, 2013. [29] X.J. Ma, L. Wang and D.D Koch, "Repeatability of corneal epithelial thickness measurements using Fourier-Domain OCT in normal and Post-LASIK eyes," Cornea, Vol. 32, iss. 12, pp. 1544–1548, 2013. [30] K.M. Rocha, C.P. Straziota, D. Stulting, B. Randleman, "Spectral-Domain OCT Analysis of Regional Epithelial Thickness Profiles in Keratoconus, Postoperative Corneal Ectasia, and Normal Eyes," J. of Refractive Surgery, Vol. 29, iss. 3, pp. 173–179, 2013. [31] X. Cui, J. Hong, F. Wang, S.X. Deng, Y. Yang, X. Zhu and J. Xu, "Assessment of Corneal Epithelial Thickness in Dry Eye Patients," Optom. Vis. Sci., Vol. 91, iss.12, pp. 1446–1454, 2014. [32] X. Cui, J. Hong, F. Wang, Y.J. Yang, "Assessment of corneal epithelial thickness in dry eye patients evaluated by Fourier-domain optical coherence tomography," ARVO 2014 Annual Meeting, pp. 1–2, 2014. [33] W. Zho and A. Stojanovic, "Comparison of Corneal Epithelial and Stromal Thickness Distributions between Eyes with Keratoconus and Healthy Eyes with Corneal Astigmatism ≤ 2.0 D," Plos One, Vol. 9, iss. 1, pp. 1–7, 2014. [34] R.H. Silverman, R. Urs, A. Roychoudhury, T.J. Archer, M. Gobbe, D.Z. Reinstein, "Epithelial Remodeling as Basis for Machine- Based Identification of Keratoconus" Cornea, Vol. 55, iss 3, pp. 1580–1587, 2014. [35] A.J. Kanellopoulos, G. Asimellis, "OCT corneal epithelial topographic asymmetry as a sensitive diagnostic tool for early and advancing keratoconus," Clin. Ophth., Vol. 4, iss. 8, pp. 2277–2287, 2014. [36] S. Wu, A. Tao, H. Jiang, Z. Xu, V. Perez and J. Wang, "Vertical and Horizontal Corneal Epithelial Thickness Profile Using Ultra-High Resolution and Long Scan Depth OCT," Plos One, Vol. 9, iss. 5, pp. 1–7, 2014. [37] D.Z. Reinstein, T.E. Yap, T.J. Archer, M. Gobbe, R.H. Silverman, "Comparison of Corneal Epithelial Thickness Measurement between Fourier-Domain OCT and Very High-Frequency Digital Ultrasound," J. of Refr. Surgery, Vol. 31, iss. 7, pp. 438–445, 2015. [38] M. Tang, Y. Li D. Huang, "Corneal Epithelial Remodeling after LASIK Measured by Fourier-Domain Optical Coherence Tomography," J. of Ophth, pp. 1–5, 2015. Nebras H. Ghaeb Al-Khwarizmi Engineering Journal, Vol. 16, No. 4, P.P. 27- 40 (2020) 39 [39] C. Temstet, O. Sandali, N. Bouheraoua, T. Hamiche, A. Galan, M. El Sanharawi and V. Borderie, "Corneal epithelial thickness mapping using Fourier-domain optical coherence tomography for detection of form fruste keratoconus. J. of Cataract and Refr. Surgery, Vol. 41, iss.4, pp. 812–820, 2015. [40] J. Tomás, A.M. Larra and L. Hanneken, "Corneal Regeneration after Photorefractive Keratectomy: A Review," J. of Optometry, Vol. 8, pp. 149–169, 2015. [41] X. Wang, J. Dong, Q. Wu, "Corneal thickness, epithelial thickness and axial length differences in normal and high myopia", BMC Ophthal., Vol. 15, iss. 9, pp. 1–5, 2015. [42] D.R. Koehn, K.J. Meyer and M.G. Anderson, "Genetic Evidence for Differential Regulation of Corneal Epithelial and Stromal Thickness," Investigative Ophth. & Visual Science, Vol. 56, iss. 9, pp. 5599-5607, 2015. [43] Y. Ma, X. He, X. Zhu, L. Lu, J. Zhu H. Zou, "Corneal Epithelium Thickness Profile in 614 Normal Chinese Children age 7–15 Years Old," Sc. Reports, Vol. 6, iss.10, pp. 234-242, 2016. [44] Q. Liang, H. Liang, H. Liu, Z. Pan, C. Baudouin and A. Labbé, "Ocular Surface Epithelial Thickness Evaluation in Dry Eye Patients: Clinical Correlations," J. of Ophth., Vol. 12, iss. 3, pp. 15-21, 2016. [45] C.H. Karabatsas, "Studying corneal epithelium in vivo: A New Tool in Clinical Practice," Advances in Ophth. and Visual System, Vol. 4, iss. 2, pp. 3–5, 2016. [46] D.J. Tritton, "Physical Fluid Dynamics", 2nd ed. Oxford, UK. Oxford Univ. Press, pp519, 1988. [47] A.H.A. Al Timemy, K.A. Shetha and N. H. Ghaeb, "A Proposed Artificial Intelligence Algorithm for Assessing of Risk Priority for Medical Equipment in Iraqi Hospital," Al- Khwarizmi Eng. J., Vol. 5, No. 1, pp 71-82, 2009. [48] N.H. Gheab, S.N. Saleem, "Comparison Study of EMG Using Wavelet and Neural Network," Al-Khwarizmi Eng. J., Vol. 4, No. 3, pp 108- 119, 2008. )2020( 27- 40، صفحة 4، العدد16دالمجلجلة الخوارزمي الهندسية م نبراس حسين غايب 40 دراسة مكونات السرعة الشعاعية والمماسية لشفاء النسيج الطالئي بعد جراحة اصالح قرنية العين البشرية نبراس حسين غائب جامعة بغداد /كلية الهندسة الخوارزمي/ الطب الحياتي هندسة قسم nebras@kecbu.uobaghdad.edu.iq:البريد االلكتروني الخالصة ) هي احدى تقنيات تصحيح االخطاء االنكسارية والتي تبتدأ بازالة الطبقة الطالئية االماية لقرنية العين PRKاستئصال القرنية باالنكسار الضوئي ( ساعة تقريًبا ، تتجديد الطبقة الطالئية االمامية التي إزيلت مسبقا لتعاود حماية العين ٤٨البشرية ثم يستخدم الليزر كعالج تصحيحي. بعد هذا اإلجراء بفترة تكون آلية إعادة نية باتجاه الجزء المركزي منها. والنسيج الطالئي) من طرف القربناء إعادة عملية تسمى والتي (هذه عملية التجديد تدأبتجديد. البشرية من الهجرة) بسرعة في اتجاهين: شعاعي وعرضي.وتسمى ب( هاتركيززيادة االنقسام) وب( ياتغيير في كثافة الخالخالل النسيج من هذا تكوين تجاه االسؤولة عن الطبقات المتداخلة بالموهي الشعاعية ( من المركبةكالل قيمة مركبات السرعة العادة نمو الطبقة الطالئية في هذه الدراسة ، تم تقدير م سرعة إعادة التي تقدر دائًما قي ، على عكس الدراسات السابقةالول مرةالمركز) من االطراف باتجاه ة الشكل يالعرضية (موج الخارجي للقرنية) والمركبة فقط. النمو لمركزي في شكل أظهرت النتائج أن شكل االتجاه لكال المكونين يتفق مع السلوك الحركي لالنقسام والهجرة ، حيث تتقلب كثافة الخلية القصوى نحو الجزء ا ٥٥٫٤٨ميكرومتر / ساعة، بينما كانت السرعة العرضية القصوى ١٦٫٨٥هي ، كانت أقصى سرعة إلعادة البناء ملم ٢االضمحالل األسي. لقطر عالج يبلغ ميكرومتر / ساعة والتي تتوافق مع سرعة الشفاء البيولوجية والعملية التي ٥٨ميكرومتر / ساعة. يعطي هذان المكونان سرعة إعادة تكوين النسيج الطالئي الكلية المطلوبة ووقت الشفاء الكلي للتحكم في فترة لعالقة بين الطبقة الطالئيةساعة. سيؤدي تقدير هذين المكونين إلى فتح الطريق لفهم اميكرومتر / ٦٠تبلغ العالج للمريض بعد الجراحة.