Stability Analysis and Analytical Solution of a Nonlinear Model for Controlled Drug Release: Travelling Wave Fronts
|
|
- Alexandra Hoover
- 5 years ago
- Views:
Transcription
1 Stability Analysis an Analytical Solution of a Nonlinear Moel for Controlle Drug Release: Travelling Wave Fronts Chontita Rattanaul, an Yongwimon Lenbury Abstract In this paper, the process of rug issolution an release from a planar matri is investigate base on two couple nonlinear partial ifferential equations propose by Göran Frenning in 3. In the moelling the process rug asorption has been isregare, assuming concentration-inepenent iffusion coefficients, using perfect sin conitions, an specializing to a planar geometry, The concentration profile of the mobile, or iffusing, the resulting moel is rather comple an has been investigate only numerically an only approimate solution have been possible. In this paper it is shown that an analytical solution can be obtaine eactly in the form of a travelling wave front. We escribe the metho for fining the analytical solutions using the travelling wave coorinate when the wave is assume to be moving at constant spee. The moel system of partial ifferential equations is transforme into two couple orinary ifferential equations, which is analyse interms of the stability of its steay state. Analytical solutions are erive in three possible cases, giving travelling wave solutions. We then iscuss a comparison between the eact solutions obtaine here an the analytical short-time approimation as well as the curves obtaine from the moifie Higuchi formula reporte by Frenning in 3. Keywors Delay ifferential equations; omega limit set; persistence; signal transuction; stability. O I. INTRODUCTION NE of the most important means of rug-elivery is the utilization of the matri systems. The objectives of controlle elivery systems are to reuce the frequency an/or to increase the effectiveness of the rug through rug localization at the target site, in orer to reuce the ose require to provie uniform rug elivery. Controlle rug elivery requires the nowlege of both physical an polymer sciences so that it is possible to prouce well characterize an reproucible osage forms, which Manuscript May 7, : Revise version receive May 7,. This wor was supporte in part by the Centre of Ecellence in Mathematics, CHE, Thailan, an Mahiol University. C. Rattanaul is with the Department of Mathematics, Faculty of Science, Mahiol University, Rama 6 Roa, Bango an the Centre of Ecellence in Mathematics, CHE, 38 Si Ayutthaya Roa, Bango, Thailan ( chontita.rat@mahiol.ac.th). Y. Lenbury is with the Department of Mathematics, Faculty of Science, Mahiol University, Rama 6 Roa, Bango an the Centre of Ecellence in Mathematics, CHE, 38 Si Ayutthaya Roa, Bango, Thailan (corresponing author, phone: ; fa: ; yongwimon.len@mahiol.ac.th). control rug entry into the boy accoring the specifications of the require rug elivery profile [l]. It is well ocumente that in this type of rug elivery, the rate of rug release is principally controlle by the elivery system itself, though it may be influence by surrouning conitions, such as ph, enzymes, ions, motility an physiological conitions []. Mathematical moeling an computation has also become an important tool in the esign of pharmaceutical proucts ealing with controlle-release mechanism. As eplaine in [3], when rug release from a matri is controlle by iffusion through the polymeric matri, its release inetics obey Fic's st an n laws [3]: J = DC () c t = DC () where J represents the iffusional flu of the rug; D the is iffusional coefficient; C is the concentration of the rug; an is the istance of iffusion. When rug release is ominate by surface erosion, Hopfenberg s equation has been foun to give goo preiction of rug elivery in spherical, cylinrical, an planar geometrical forms [4]. Accoring to [3], we may istinguish between two conceptually ifferent scenarios, epening on the relative magnitues of the rug concentration an solubility in the matri. If the rug concentration is sufficiently low so that all rugs can be issolve, an the issolution process procees rapily enough, we may easily etermine the release rate [3]. In this type of release, all rugs may be assume to be completely issolve in before any release has occurre, an the rug concentration in the matri can therefore be solve from the heat conuction equation [3]. In the secon, more general, type of rug release from a matri, where the rug concentrations are higher, or the solubility is low, two forms of rugs, namely the soli an issolve forms, coeist, in the matri an the process becomes istinctly more comple. For this more general situation when rug loaing is much higher than rug solubility ( C C ), the moel propose by Higuchi [5] has s shown to perform well for planar matri uner the perfect sin assumption, although it was originally formulate for rug release from ointment bases containing rugs in suspension. Issue 3, Volume 6, 35
2 The original Higuchi moel has been subject to a great eal of generalizations an improvements [6-]. Recently, a similar moel has been stuie by Frenning an Strømme [3], who investigate the problem of rug release from spherical pellet units into the finite volume of issolution meium, uner the assumption that some of the issolve rug coul become immobilize by absorption to the pellet constituents. In [4], this moel was reformulate by isregaring rug absorption, an assuming that the iffusion coefficient is concentrationinepenent. Using perfect sin conitions, an analytical short-time approimation to the solution was erive [4]. It is the purpose of this paper to show that an eact solution can be obtaine analytically in the form of a travelling wave front. We escribe the metho for fining the analytical solutions using the travelling wave coorinate in the situation that the wave is presumably moving at a constant spee. The analytical solutions are then iscusse in comparison to the analytical short-time approimation erive in [4]. II. REFERENCED MODEL In [4], a planar matri system whose normal is in the irection was investigate. The assumptions are that the lateral imensions of the system are much larger than its thicness L, so that the process of rug release coul be effectively consiere to be one-imensional. The bounary at = is assume to be impenetrable to the rug, while the matri is in contact with the liqui at = L. The perfect sin conition assumes that the matri is in contact with a well-mie issolution meium, the volume of which is sufficiently large so that it can be assure that the rug concentration is virtually zero at all times. In orer to simplify the analysis, it was assume by Frenning [4] that liqui absorption occurs at a much faster rate than rug issolution an subsequent release. Thus, the matri which contains finely isperse soli rug is fully wette in the initial state. Also, in the initial state, entire rug is present in the soli form. Letting ( t, ) be rug concentration in the issolve phase an s( t, ) the concentration of rug in the soli phase, it is then possible to escribe rug issolution an release by the following equations [4]. t = st (3) / s = s 3 ( εc ) (4) t s where t is the time, c s the saturation constant, ε the porosity, D the rug iffusivity, the issolution rate. The initial conitions are (, ) = s(, ) = which means that all rug is assume to be present in the soli form in the initial state [4]. Similarly, the bounary conitions are = = s( t, L ) = which follow from the assumption that the interface at = is impenetrable to the rug, while the rug concentration at = L is ept at zero as a consequence of the sin conition. III. MODEL ANALYSIS 3. Traveling wave coorinate We introuce the travelling wave coorinate z = t, where is the constant velocity at which the wave is assume to be moving. By substituting z in (3) an (4), we obtain a couple system of orinary ifferential equations with respect to z: v = D + s (5) / s = s 3 ( γ ) (6) where () enotes the erivative with respect to z, stans for, an γ stans for ε cs. Integrating (5) an combining with (6) lea us to a single secon-orer ifferential equation terms of as / 3 D + v + ( v + D ) ( γ ) = (7) / 3 Now, by letting X =, an Y = X, we can write (7) as X = Y (8) / 3 Y = ( vx + DY ) ( X γ ) Y (9) / 3 D D 3. Dynamical analysis Before we erive the analytical solution, a stability analysis may be carrie out on the moel system written in the form of equations (3) an (4). The system possess only one nontrivial equilibrium state in the feasible region, namely ( X, Y ) = ( γ,). The Jacobian matri of the system (8) (9) about its steay state ( γ,) is /3 /3 J = ( DY + X) ( X γ) ( X γ) + /3 /3 /3 /3 D 3 D( DY X) 3 ( DY X) D + + At ( X, Y ) = ( γ,) the Jacobian becomes J( γ,) = / 3 γ D D whose eigenvalues are / 3 γ λ, = ± ( ) + 4 D D D which are real an opposite in signs. Therefore, the steay state ( X, Y ) = ( γ,) is a sale point. Issue 3, Volume 6, 35
3 Y = Y becoming unboune as time passes. The corresponing plots of the concentration of issolve rug an its rate of change as functions of z for the same parameter values as in Fig. - are seen in Fig. 3 an 4 respectively X = X Fig. Phase portrait of the system (8)-(9) in the case <. Here D =.8, =, v =, γ =. z Fig. 3 Plot of, the issolve rug concentration, as a function of z for the same parameter values as in Fig. -. Y Y X Fig Computer simulation of the system (8)-(9), showing the solution trajectory in the ( X, Y ) -plane in the case that D =.8, v =, =, γ =, X () =, Y () =.. In what follows, we are able to erive analytical solutions, in terms of the travelling wave coorinate, in 3 of the cases. Fig. shows the solution trajectory in the phase plane starting from the initial point at ( X, Y ) = (,.) an z Fig. 4 Plot of as a function of z for the same parameter values as in Fig. -. Issue 3, Volume 6, 353
4 IV. ANALYTICAL SOLUTION In orer to erive an analytic solution for the moel equation (5), we first let 3 / C = + D () We will see a solution of the form m n = AC + BC () so that we have 3 / C C = + D () Also, Eq. () gives 3 / D m n = C ( AC + BC ) (3) Substituting () an (3) into (7), we obtain 3 / 3/ D m n C C + C γ C ( AC BC ) / = Rearranging the above equation yiels 3 / 5 / AD m C C γ + + C C + C / 3 5 / 3 5 / 3 BD 5 / 3 C n+ + = (4) On inspection, we observe that we may fin eact solutions in the following two possible cases. Case : m =, n = For these values of m an n, Eq. (4) becomes 3 / 5 / AD C C + γ C C + C / 3 5 / 3 5 / 3 BD / + 5 / 3 C = (5) By letting AD = γ (6) equation (5) reuces to 3 / 5 / BD / C C C + C = (7) 5 / 3 5 / 3 If we let BD α = >, β = < (8) 5 / 3 BD 3 we can write Eq. (7) as C = β + α C which can be easily solve, yieling C = t an( αβ z + αk ) (9) α where K is the constant of integration. Substituting (7) into (3), rearranging an using (6), one obtains = ( tan 3/ ( αβ z+ αk) + cot / ( αβ z+ αk) 3/ ) + γ () α So that we have satisfie., z = = the following equation must be 3/ / 3 / tan αk + cot αk = γα () We can then obtain the level of rug in soli form by integrating (5) so that Fig. 5 Travelling wave solution in Case for the concentration of rug in the ilute form, given by (3), plotte as a function of for ifferent time. 3 / s = C + l or, 3 / s = tan ( αβ z + αk) + l () 3 / α where l is the constant of integration. Thus, a travelling wave solution is given as = ( tan 3/ ( αβ( t) + αk) + cot / ( αβ( t) + αk) 3/ ) α + γ (3) 3/ s = tan ( αβ ( t) + αk) + l (4) 3/ α where α, β, A, an B are given by (6) an (8), an (). Issue 3, Volume 6, 354
5 s C = (7) Jz + K where K is the constant of integration. Upon substitution (7) into (3) an using (4) we are le to AD 3 / = ( Jz + K) + γ (8) So that we have, z = = we nee the following equation to be satisfie. AD γ 3/ = (9) K Using (5) then gives s = + l (3) 3 / ( Jz + K) where l is the constant of integration. Thus, a travelling wave solution to (3)-(4) is given as AD 3/ = ( J ( t) + K) + γ (3) s = + l (3) 3 / ( J ( t) + K) where α, β, A, an B are given by (4), (6), an (9) Fig. 6 Travelling wave solution in Case for the concentration of rug s in the soli form, given by (4), plotte as a function of for ifferent time t. The analytical travelling wave solution given by (3)-(4) is shown in Fig. 5 an 6 to move in the irection of ecreasing as time elapses. Here, π α =, K =, l =, =, =, an γ = Case : m =, n = In this case, equation (4) becomes 3 / 5 / C C + γ C C / 3 5 / 3 AD 5 / BD + C + C = (3) 5 / 3 5 / 3 If we let BD = γ (4) then (3) reuces to 3 / AD 5 / C C C = 5 / 3 5 / 3 (5) or C = JC, where ( AD ) J = (6) 5 / 3 3 Solving (6), we obtain Fig. 7 Wave front solution in Case for the concentration of rug in the issolve form, given by (3), plotte as a function of for ifferent time t. Issue 3, Volume 6, 355
6 s Fig. 8 Wave front solution in Case for the concentration of rug s in the soli form, given by (3), plotte as a function of for ifferent time t. The graphs of the wave front solution given in (3)-(3) are shown in Fig. (7)-(8) to move in the irection of ecreasing as time elapses. Here, D =., K =.4, v =, =, γ =.5, an l =. 3 Case 3: m =, n = n +. For convenience, we let 3a 3b A =, B = in (), an rop the prime on the eponent n. Thus, we are now looing for a solution in the form 3 3 / n / = ( ac + bc + ) (33) such that n C = AC + B C (34) Differentiating (33) an using (34), we obtain 3 3 / n / n = ac + b( n + ) C ( A C + BC ) (35) Substituting (33) an (35) in (7), choosing A = a, an B = b, an rearranging, we are le to 3D 3 b 3D 3a a ab + + ab( n + ) C + + C D D n+ / 3/ 3D n / γ + b ( n + ) C + C / 3 3aD 5/ 3bD n+ 3/ C C (36) 5/ 3 5/ 3 + = 5/ 3 Upon observation, we see that if we let n = then Eq. (36) can be written as 9 3 γ 3bD 9 3 3/ abd + b C + + C+ a D+ a C /3 5/3 4 v v 4 3aD 5/ + C = (37) 5/3 5/3 v v Thus, if we let γ 3bD γ + =, or b = / 3 5 / 3 v v 3D (38) 3aD =, or a = 5 / 3 5 / 3 v v 3D then (37) will be satisfie. With the above choice in (38), we are in the case where v =, an (34) becomes C = ac + bc / which can be easily solve to yiel 3 / 3 az + C( z) = b e / 3 a Substituting (39) into (), we obtain (39) z 3/ 3D + D = C = γ e D which can be irectly solve, using an integrating factor, with, z = = to yiel 3 z / D+ z / D = γ ze + γ e Therefore, we have erive an analytical solution in terms of the travelling wave coorinate as follows. 3 ( ) ( ) / t D = γ t e e γ (4) an, similarly to the previous cases, 3 ( t) / D+ s = γ De + l (4) l being a constant of integration. The graphs of the wave front solution given in (4)-(4) are shown in Fig. (9)-() with =, 3 + γ =, D = 5 e, =., an l = γ + De. Issue 3, Volume 6, 356
7 V. DISCUSSION To illustrate that the analytic solution can allow us to investigate the impact of ifferent values of the system s physical properties, such as the saturation constant, the porosity, the rug iffusivity, or the issolution rate, on the waveform structures, we show in Fig. the plot of the travelling wave front of concentration of rug in issolve form for a ifferent set of parametric values, namely, γ =, D =, =., an =. Fig. 9 Wave front solution in Case 3 for the concentration of rug in the issolve form, given by (4), plotte as a function of for ifferent time t. s Fig. Wave front solution in Case 3 for the concentration of rug in the issolve form, given by (4), plotte as a function of for ifferent time t: γ =, D =, =., =, an. z = = Fig. Wave front solution in Case 3 for the concentration of rug in the issolve form, given by (4), plotte as a function of for ifferent time t. Finally, to put our analytical solution in the contet of numerical an approimate solutions reporte in earlier literatures, we refer to the results shown in [5] an [4]. In [4], fractions of release rug calculate numerically using the moel consisting of (3) an (4) is shown for ifferent values of saturation constant c s. The curves are compare with the analytical short-time approimation as well as the curves obtaine from the moifie Higuchi formula. We observe that the Higuchi solution an the analytical short-time approimation are close to the eact solution only uring the eponential phase of rug release, but offer quite a poor estimate of the release rug near the saturation perio. Issue 3, Volume 6, 357
8 Our eact solution is therefore etremely valuable as it can be use to accurately investigate effects of ifferent values of physiological or physical factors in the controlle rug release system, such as ifferent saturation constants, iffusivity, or porosity. These parameters an their impact on rug elivery are of great importance in the esign of pharmaceutical proucts which are well characterize an in reproucible osage forms, capable of controlling rug entry into the boy accoring the specifications of the require rug elivery profile. VI. CONCLUSION Intensive continuing research, both eperimentally an theoretically, has been on going, especially in the last ecae or so, on the topic of controlle-release technology, ue to its crucial role in the formulation of pharmaceutical proucts [- 4, -5]. Here, we have ae to the current nowlege by showing that an eact solution can be erive for the moel of rug release formulate by Frenning [4] where the lateral imensions of the system are assume to be much larger than its thicness L. The bounary at = is assume to be impenetrable to the rug, while the matri is in contact with the liqui at = L. The matri is assume to be in contact with a well-mie issolution meium, the volume of which is large enough to ensure that the rug concentration is virtually zero at all times. The resulting moel is relatively comple, involving a nonlinear fractional power of the epenent variable. Previous wors have been able to obtain numerical solution an, only by restriction their attention to the initial phase of rug release, an analytical approimate solution can be foun. We have relie on the wor in [7-9], using the travelling wave coorinate, an appropriate choice of combination of powers of the solution an its erivative, we have shown that an eact solution for the concentrations of rug in the soli an issolve forms can be foun as functions of the travelling coorinate, while the travelling wave velocity is foun to be negative as epecte. The methoologies for obtaining the eact solution presente in this stuy can be applie to other problems [, ] which amit travelling wave solutions. In the case that the system is escribe by a moel consisting of couple partial ifferential equations, they can be reuce to a single secon orer ifferential equation which is more easily solve. In Fig., we can see that the rug concentration in the issolve form increases (as t tens to infinity) while the wave of the issolve rug concentration wave front is moving at a constant spee of.. It might be wonere whether it woul have been better to simply solve the moel equations for numerical solutions. To answer such a question, the moel system shoul be solve numerically so that the numerical simulation coul be compare with the analytical solution, which is the subject for our future wor. It can be argue at this juncture that the ability to erive the eact solution coul be quite valuable in the investigation of impacts of various conitions on rug elivery, since the ability to accurately preict the outcome of ifferent choices of material, composition, rug properties, an others, on the release mechanism an inetics is of crucial importance in terms of meical treatment an health care purposes. ACKNOWLEDGMENT This research is supporte by Mahiol University an the Centre of Ecellence in Mathematics, CHE, Thailan, an Mahiol University. REFERENCES [] J.E. Mocel, B.C. Lippol, Zero orer rug release from hyrocolloi matrices, Pharm. Res., Vol., 993, pp [] J. Heller, Bioegraable polymers in controlle rug elivery, Crit. Rev. Ther. Drug Carrier Syst., Vol., N, 984, pp [3] J. Cran, The Mathematics of Diffusion, Ofor University Press, Lonon, 975. [4] P. Colombo, R. Bettini, G. Massimo, P.L. Catellani, P. Santi, N.A. Peppas, Drug iffusion front movement is important in rug release control from swellable matri tablets, J. Pham. Sci. Vol.87, 995, pp [5] T. Higuchi, Rate of release of meicaments from ointment bases containing rugs in suspension, J. Pharmac. Sci., 5 (96), [6] T.J. Roseman, W. Higuchi, Release of meroyprogesterone acetate from a silicone polymer, J. Pharm. Sci., Vol. 59, 97, pp [7] K. Tojo, Intrinsic release rate from matri-type rug elivery systems, J. Pharm. Sci., Vol.74, 985, pp [8] D.R. Paul, S.K. McSpaen, Diffusional release of a solute from a polymeric matri, J. Membr. Sci. Vol., 976, pp [9] P.I. Lee, Diffusional release of a solute from a polymeric matriapproimate analytical solutions, J. Membr. Sci., Vol.7, 98, pp [] A.L. Bunge, Release rates from topical formulations containing rugs in suspension, J. Control. Release, Vol.5, 998, pp [] M.J. Abehoaie, Y.-L. Cheng, Diffusional release of a isperse solute from planar an spherical matrices into finite eternal volume, J. Control. Release, Vol.43, 997, pp [] Y. Zhou, X.Y. Wu, Theoretical analyses of isperse-rug release from planar matrices with bounary layer in a finite meium, J. Control. Release, Vol.84,, pp. -3. [3] G. Frenning, M. Strømme, Drug release moele by issolution, iffusion, an immobilization, Int. J. Pharm., Vol.5, 3, pp [4] G. Frenning, Theoretical investigation of rug release from planar matri systems: effect of a finite issolution rate, J. Control. Release, Vol.9, 3, pp [5] D. Sun, V.S. Manoranjan, H.-M. Yin, Numerical solution for a couple parabolic equations arising inuction heating processes, Discrete an Continuous Dynamical System, 7, pp [6] R.L. Buren, J.D. Faires, A.C. Reynols, Numerical Analysis n Eition, Prinle, Weber an Schmit. Inc, USA, 98. [7] M.O. Gomez, C.M. Chang, V.S. Manoranjan Eact solution for contaminant transport with nonlinear sorption, Appl. Math. Lett., Vol.9, 996, pp [8] V.S. Manoranjan, T.B. Stauffer, Eact solution for contaminant transport with inetic langmuir sorption, Water Resources Research, Vol.3, No.3, 996, pp [9] C.M. Chang, V.S. Manoranjan, Travelling wave solutions of a contaminant transport moel with nonlinear sorption, Mathl. Comput. Moelling., Vol.8, No.9, 998, pp. -. [] E. Kreyszing, Avance Engineering Mathematics 8th Eition, John Wiley an Sons. Inc, USA, 3. [] P. Prasertsang, V.S. Manoranjan, Y. Lenbury, Analytical travelling wave solutions of a ental plaque moel with nonlinear sorption, Nonlinear Stuies, Vol. 8, No.,, pp Issue 3, Volume 6, 358
9 [] S.-P. Tang, An eploratory research on the nanovectors for rug elivery an for gene therapy, WSEAS Transactions on Biology an Biomeicine, Volume 5, Issue, 8, pp [3] S.-P. Tang, An eperimental comparing analysis research on the nanovectors for rug elivery an for gene therapy, WSEAS Transactions on Biology an Biomeicine, Vol. 7, No.,, pp. -. [4] T. Schroer., A. Quintilla., J. Setzier., E. Birtaian., W. Wenzel., an S. Brase., Joint eperimental an theortical Investigation of the propensity of peptois as rug carriers. WSEAS Transuctions on Biology an Biomeucine, Vol. 4, Issue, 7, pp [5] T. Hara, S. Iriyama, M. Ohya, Estimation of rug shape in rug elivery system by simulation, Proceeings of the th WSEAS International Conference on Mathematics an Computers in Biology an Chemistry, 9. Issue 3, Volume 6, 359
Some New Solutions of a Reaction Diffusion Model for Controlled Drug Release using Travelling Wave Coordinate Transformation
Some New Solutions of a Reaction Diffusion Moel for Controlle Drug Release using Travelling Wave Coorinate Transformation Chontita Rattanaul, Yongwimon Lenbury Abstract Being able to eliver rugs to a targete
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More information18 EVEN MORE CALCULUS
8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationA SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationUnit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule
Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin
More informationSummary: Differentiation
Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More information5-4 Electrostatic Boundary Value Problems
11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationP. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,
More informationMATH 205 Practice Final Exam Name:
MATH 205 Practice Final Eam Name:. (2 points) Consier the function g() = e. (a) (5 points) Ientify the zeroes, vertical asymptotes, an long-term behavior on both sies of this function. Clearly label which
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationEXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL
Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University
More informationwater adding dye partial mixing homogenization time
iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,
More informationTutorial 1 Differentiation
Tutorial 1 Differentiation What is Calculus? Calculus 微積分 Differential calculus Differentiation 微分 y lim 0 f f The relation of very small changes of ifferent quantities f f y y Integral calculus Integration
More information2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math
Lesson 2-7 Fitting a Moel to Data I BIG IDEA If you etermine from a particular set of ata that y varies irectly or inversely as, you can graph the ata to see what relationship is reasonable. Using that
More informationStudy on aero-acoustic structural interactions in fan-ducted system
Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationChapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies
hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationAsymptotics of a Small Liquid Drop on a Cone and Plate Rheometer
Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationApproaches for Predicting Collection Efficiency of Fibrous Filters
Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Moels for CFD ME 448/548 Notes Geral Recktenwal Portlan State University Department of Mechanical Engineering gerry@p.eu ME 448/548: D Convection-Di usion Equation
More informationEXACT TRAVELING WAVE SOLUTIONS FOR A NEW NON-LINEAR HEAT TRANSFER EQUATION
THERMAL SCIENCE, Year 017, Vol. 1, No. 4, pp. 1833-1838 1833 EXACT TRAVELING WAVE SOLUTIONS FOR A NEW NON-LINEAR HEAT TRANSFER EQUATION by Feng GAO a,b, Xiao-Jun YANG a,b,* c, an Yu-Feng ZHANG a School
More information5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask
5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of
More informationMath 180, Exam 2, Fall 2012 Problem 1 Solution. (a) The derivative is computed using the Chain Rule twice. 1 2 x x
. Fin erivatives of the following functions: (a) f() = tan ( 2 + ) ( ) 2 (b) f() = ln 2 + (c) f() = sin() Solution: Math 80, Eam 2, Fall 202 Problem Solution (a) The erivative is compute using the Chain
More informationBasic Differentiation Rules and Rates of Change. The Constant Rule
460_00.q //04 4:04 PM Page 07 SECTION. Basic Differentiation Rules an Rates of Change 07 Section. The slope of a horizontal line is 0. Basic Differentiation Rules an Rates of Change Fin the erivative of
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationFurther Differentiation and Applications
Avance Higher Notes (Unit ) Prerequisites: Inverse function property; prouct, quotient an chain rules; inflexion points. Maths Applications: Concavity; ifferentiability. Real-Worl Applications: Particle
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationNested Saturation with Guaranteed Real Poles 1
Neste Saturation with Guarantee Real Poles Eric N Johnson 2 an Suresh K Kannan 3 School of Aerospace Engineering Georgia Institute of Technology, Atlanta, GA 3332 Abstract The global stabilization of asymptotically
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationAdjoint Transient Sensitivity Analysis in Circuit Simulation
Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationPlanar sheath and presheath
5/11/1 Flui-Poisson System Planar sheath an presheath 1 Planar sheath an presheath A plasma between plane parallel walls evelops a positive potential which equalizes the rate of loss of electrons an ions.
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationEvaporating droplets tracking by holographic high speed video in turbulent flow
Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationThe theoretical limits to the power output of a muscle tendon complex with inertial and gravitational loads
Receive 1 November Accepte 17 March 3 Publishe online 11 June 3 The theoretical limits to the power output of a muscle tenon complex with inertial an gravitational loas Apostolos Galantis * an Roger C.
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationA General Analytical Approximation to Impulse Response of 3-D Microfluidic Channels in Molecular Communication
A General Analytical Approximation to Impulse Response of 3- Microfluiic Channels in Molecular Communication Fatih inç, Stuent Member, IEEE, Bayram Cevet Akeniz, Stuent Member, IEEE, Ali Emre Pusane, Member,
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationNON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS
ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This
More informationModelling of Three Phase Short Circuit and Measuring Parameters of a Turbo Generator for Improved Performance
Moelling of Three Phase Short Circuit an Measuring Parameters of a Turbo Generator for Improve Performance M. Olubiwe, S. O. E. Ogbogu, D. O. Dike, L. Uzoechi Dept of Electrical an Electronic Engineering,
More informationUnified kinetic model of dopant segregation during vapor-phase growth
PHYSICAL REVIEW B 72, 195419 2005 Unifie kinetic moel of opant segregation uring vapor-phase growth Craig B. Arnol 1, * an Michael J. Aziz 2 1 Department of Mechanical an Aerospace Engineering an Princeton
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationCharacterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy
JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu
More informationMULTISCALE FRICTION MODELING FOR SHEET METAL FORMING
MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of
More informationComparative Approaches of Calculation of the Back Water Curves in a Trapezoidal Channel with Weak Slope
Proceeings of the Worl Congress on Engineering Vol WCE, July 6-8,, Lonon, U.K. Comparative Approaches of Calculation of the Back Water Curves in a Trapezoial Channel with Weak Slope Fourar Ali, Chiremsel
More informationMore from Lesson 6 The Limit Definition of the Derivative and Rules for Finding Derivatives.
Math 1314 ONLINE More from Lesson 6 The Limit Definition of the Derivative an Rules for Fining Derivatives Eample 4: Use the Four-Step Process for fining the erivative of the function Then fin f (1) f(
More informationCHAPTER SEVEN. Solutions for Section x x t t4 4. ) + 4x = 7. 6( x4 3x4
CHAPTER SEVEN 7. SOLUTIONS 6 Solutions for Section 7.. 5.. 4. 5 t t + t 5 5. 5. 6. t 8 8 + t4 4. 7. 6( 4 4 ) + 4 = 4 + 4. 5q 8.. 9. We break the antierivative into two terms. Since y is an antierivative
More informationCritical Size and Particle Growth
Critical Size an article Growth rof. Sotiris E. ratsinis article Technology Laboratory Department of Mechanical an rocess Engineering, ETH Zürich, Switzerlan www.ptl.ethz.ch 1 Nucleation-Conensation A
More information