SMOOTHED PARTICLE HYDRODYNAMICS METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM
|
|
- Abigail Mathews
- 5 years ago
- Views:
Transcription
1 The 5th Internatonal Symposum on Computatonal Scences (ISCS) 28 May 2012, Yogyakarta, Indonesa SMOOTHED PARTICLE HYDRODYNAMICS METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM Dede Tarwd 1,2 1 Graduate School of Natural Scence and Technology, Kanazawa Unversty Kakuma, Kanazawa , Japan 2 Faculty of Mathematcs and Natural Scences, Insttut Teknolog Bandung Jl. Ganesha 10, Bandung 40132, Indonesa emal: dede.tarwd@yahoo.com Astract. Smoothed partcle hydrodynamcs (SPH) s developed for modellng of meltng and soldfcaton. Enthalpy method s used to solve heat conducton equatons whch nvolved movng nterface etween phases. At frst, we study the meltng of floatng ce n the water for twodmensonal system. The ce ojects are assumed as sold partcles floatng n flud partcles. The flud and sold moton are governed y Naver-Stokes equaton and asc rgd dynamcs equaton, respectvely. We also propose a strategy to separate sold partcles due to meltng and soldfcaton. Numercal results are otaned and plotted for several ntal condtons. Keywords: SPH, enthalpy method, Stefan prolem, movng nterface 1. INTRODUCTION Stefan prolem s related to heat transfer prolem nvolvng phase-change such as sold to lqud (meltng) or lqud to sold (soldfcaton). Ths ssue has the attenton from many researchers snce the last few decades. Several methods are used to solve the Stefan prolem oth analytcal and numercal solutons. Due to ts complexty, the analytcal soluton s lmted to one-dmensonal and smple two-dmensonal prolem. The detal of the analytcal and numercal soluton of the one-dmensonal Stefan prolem s dscussed n [1]. As to the numercal soluton, n general, t can e dvded nto two methods, the grd/mesh-ased method and the grdless/meshless-ased method. Fnte dfference, fnte element, and fnte volume method are grd/mesh-ased methods often used n solvng Stefan prolem. Voller and Shada [2] used fnte dfference method va enthalpy formulaton for solvng two-dmensonal Stefan prolem. Moreover, the fnte element method [3-5] and the fnte volume method [6, 7] are wdely employed for solvng oth one-dmensonal and two-dmensonal Stefan prolem. The smoothed partcle hydrodynamcs (SPH) method s a meshless method. Ths method s a partcle-ased nterpolaton whch was orgnally used for modellng astrophyscal prolems. However, ths method s eventually developed n flud dynamcs prolems. In ths paper, we propose SPH method for solvng two-dmensonal Stefan prolem. The SPH method whch s appled n meltng prolem together wth flud flow could e seen n [13]. In other case, Monaghan et al. [8] used SPH method n solvng soldfcaton prolem, ut the flud flow s not consdered n the model. In the present case, we employ SPH method to smulate ce meltng prolem n twodmensonal systems. For the reverse prolem that s water soldfcaton, the mathematcal formulaton s the same ut dfferent n the ntal and oundary condton. In the meltng or soldfcaton case, two regons of the system are consdered, lqud (water) and sold (ce) separated y phase change nterface. Moreover, the nterface s always movng n each tme step. The system leads us to the two-phase Stefan prolem.
2 DEDE TARWIDI Two heat conducton equatons should e solved n sold and lqud regon, ut the oundary of the doman s not known. Our goal s to determne temperature feld T ( x, y, t ) that satsfy a heat conducton equaton nsde sold regon, a heat conducton equaton nsde lqud regon. Also, T ( x, y, t ) must satsfy nterface condton, ntal and oundary condton. 2. GOVERNING EQUATION There are three categores for governng equatons of the system: flud flow, sold moton, and heat transfer (energy equaton). The flud moton s assumed weakly compressle. The governng equatons for weakly compressle flud are the Naver- Stokes and the contnuty equatons as follow: Dv 1 = p + F (1) Dt ρ Dρ ρ Dt = v..(2) where v, ρ, and p are velocty, densty and pressure of the flud, respectvely. Whle F s external force per unt mass whch can e gravtatonal acceleraton or repulsve force per unt mass. For the present case, surface tenson s assumed not to e sgnfcant, so t can e neglected. The equaton of state n smple form s gven y 2 p = c ρ ρ + p..(3) ( ) 0 atm where c s speed of sound, ρ0s densty reference, and patm s atmospherc pressure. The choce of speed of sound s very mportant n the smulaton snce t nfluences the densty varaton. In order to keep densty fluctuatons less than 1%, c = 200gH s chosen where g s gravty acceleraton and H s heght of water level. The governng equaton of heat transfer and sold moton s dscussed further n Secton 3 and Secton 6, respectvely. 3. ENTHALPY METHOD Enthalpy method s usually used n modellng that nvolved phase-change such as meltng and soldfcaton. Ths method refers to reformulate the heat conducton equaton n dfferent phases nto the enthalpy equaton. Also, y usng enthalpy equaton the latent heat (heat of fuson) s accounted n calculatng the temperature feld. In the meltng case, ths latent heat s needed to reak up the ndng n the sold structure. Otherwse, n the soldfcaton case, the latent heat s released. One of the advantages of the enthalpy method s, the heat conducton equaton n sold and lqud regon can e solved wthout need to know the nterface poston and the heat flux n the nterface s automatcally contnuous. Moreover, f we know the relatonshp etween enthalpy and temperature then we can calculate temperature feld. If we assume heat transfer y conducton only, then the enthalpy equaton s H 1 = ( k T )...(4) t ρ
3 SPH METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM where H s enthalpy, ρ s densty and k s thermal conductvty. The relatonshp etween enthalpy and temperature s gven y [1] H Tm +, H 0 (sold) cs T = Tm, 0 < H < L (nterface) (5) H L Tm +, H L (lqud) cl T m s meltng pont, cs and c l are specfc heat of sold and lqud respectvely, and L s latent heat. 4. SPH FORMULATION The man dea of the SPH method s to represent a functon or ts dervatve nto ntegral representaton and dscretze t nto set of partcles. The detal aout SPH method s dscussed n [14]. The ntegral nterpolant of functon f ( x ) s f ( x) = f ( x ) δ ( x x ) dx...(6) Ω where δ ( x x ) s the Dract delta functon and Ω s the volume of the ntegral that contan x. If we replace the delta functon y smoothng functon or kernel functon, W, we get kernel aproxmaton of functon f ( x ) : f ( x) = f ( x ) W ( x x, h) dx (7) Ω where h s smoothng length that represents nfluence area of kernel functon. We also f x : can get the kernel aproxmaton of the dervatve of functon ( ) f ( ) = f ( ) W (, h) d x x x x x (8) Ω The followng equatons are derved y usng SPH dscretzaton: dv a p + p a (Naver-Stokes equaton ) = m aw ( ra r, h)...(9) dt ρaρ dρa (contnuty equaton) = m ( ra r ) aw ( ra r, h)... (10) dt Equaton (9) and (10) mean that the velocty and densty n a partcle can e otaned y summng up the contruton of the neghorng partcles. The contruton s governed y kernel functon W. Therefore the choce of W wll affect the accuracy of the smulaton. There are many choces of kernel functon, ut for the present smulaton, the followng cuc splne kernel s used: q + 3 q, 0 q < W ( q, h) = 2 4 ( 2 q), 1 q < 2.(11) 7πh 0, otherwse
4 DEDE TARWIDI where q = r r / h. a The heat transfer etween partcles (water, ce, and sold oundary) s represented y the temperature and enthalpy whch are carred y SPH partcles. In the present smulaton there are only two types of heat transfer to e taken nto account of the conducton and convecton. But, the heat transfer y convecton s automatcally fulflled ecause the lqud partcle s always moves accordng to the Naver-Stokes equaton [13]. The heat transfer y conducton n SPH form s gven y [9] dh ( ) ( ) (, ) a m kak ra r aw ra r h = Ta T (12) 2 2 dt ρ ρ k + k r r + η a a a where η = 0.01h. The steps n constructng the SPH heat conducton of equaton (12) s refly dscussed n [11]. Note that y usng equaton (12), n nterface etween sold and lqud phase, the heat flux s automatcally contnuous. Once we otan the enthalpy of the partcles, the temperature partcles and phase-change transton are determned y usng equaton (5). The leap-frog tme steppng s employed for equaton (9) and equaton (10) whle explct Euler tme steppng s used for equaton (12). 5. BOUNDARY TREATMENT In SPH modellng, t s very mportant to treat the partcles that approach wall oundares. In ths smulaton, the wall oundares are assumed as sold oundary partcles. We need these sold oundary partcles to prevent the water partcles penetrate the wall oundares. The water partcles approachng wall oundares wll experence repulsve force comng from sold partcles. The force s experenced y a water partcle j normal to the sold oundary partcle s gven y [10] f = n R j ( y ) P ( x ).(14) where n s normal vector of the sold oundary. Here, R y A q 1 ( ) (1 ) q ( ) 1 = and P( x) = 2 1+ cos( πx p ) where x s dstance of projecton water partcle j on tangent vector of the sold oundary partcle, y s perpendcular dstance of water partcle j from the sold oundary partcle (see Fgure 1), q = y /(2 p) and p s ntal partcle spacng. The choce of partcle spacng of sold oundary partcles depends on the ntal water partcle spacng. If p s too small, then the water partcle wll experence a repulsve force efore t approaches the wall oundary. Whereas f p s too g, some of water partcles may penetrate the oundary wall efore experencng the repulsve force. 6. FLOATING ICE MODEL There are two ways for modelng ce (sold) ojects. Frst, consder ce partcles as f they were water partcles. Hence, oth sold and lqud partcles are solved y usng Naver-Stokes equaton. In ths way the densty of ce partcles s also calculated y usng equaton (2). But, they are not updated n each tme step n order to keep the consstency of the densty of ce. However, t makes prolems ecause ce densty s 917 kg/m 3 and water densty s 1000 kg/m 3 (the dfference s too wde) whle the Naver-Stokes equaton
5 SPH METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM s desgned for weakly compressle. Consequently, ce partcles wll have low pressure, and t makes water partcles penetrate ce ojects. The second way s to consder ce partcles as sold oundary (wall) partcles. In ths way, If the water partcles approachng ce partcles, they wll experence repulsve force. The prolem s we have to fnd normal vector of ce ojects n each tme step. It s slghtly dffcult snce ce formaton s always changng. In meltng case, t may e easy to fnd normal vector of ce wall ecause at ntal we set ce partcles neatly. But, n soldfcaton case, t s rather dffcult to fnd normal vector of ce wall whch s formed from water partcles. It s ecause the ce partcles formed are unstructured. In the present smulaton, the ce ojects whch float n water are assumed as wall oundary partcles. So that for each water partcle that approaches the ce partcles, t wll experence repulsve force whch s calculated accordng to equaton (14). But, the challenge s n determnng the normal vector of ce wall snce the wall of ce s changng and movng n each tme step. The normal vector of ce partcle due to water partcle j can e determned y projectng vector r j onto the tangent vector of partcle. To fnd tangent vector, t s necessary to fnd two neghorng partcles of partcle says -1 and +1 wth the opposte poston, see Fgure 1. The tangent vector t s r + 1 r 1. Hence the unt normal vector s n = ( r proj r ) / r proj r j t j j t j where projt r j s projecton of vector r j onto tangent vector t. The advantage of ths strategy s that n s always outward normal vector. Fgure 1. Normal vector and tangent vector The repulsve force of sold partcle s calculated y summng up the contruton of surroundng water partcles [15] f = f j (15) j WPs where WPs denotes water partcles. The equatons of moton for ce ojects are derved from the asc rgd dynamcs equaton. The translatonal and rotatonal equatons of sold moton are dv dω M = mf and I = m ( r R0 ) f..(16) dt IPs dt IPs where M, I, V, Ω, and R 0 are mass, moment of nerta, velocty, rotatonal velocty, and center of mass of the ce ojects, respectvely. Here, IPs represents ce partcles. The velocty for each ce partcle s gven y
6 DEDE TARWIDI dr ( ) V Ω r R..(17) dt = + Thus the poston of ce partcles for each tme step can e calculated y ntegratng equaton (17) n tme. An ce oject may e separated nto some ce ojects durng meltng process and also new ce ojects may e formed durng soldfcaton process. To control the moton of these ojects, there should e a strategy to separate the ce partcles nto certan ce ojects, so that equaton (16) can e employed. A smple strategy to separate the sold partcles has een proposed n [12]. Frst, each ce partcle s assgned an ndex from 1 to N, where N s the numer of ce partcles. Second, each ce partcle ndex s updated teratvely y the maxmum ndex of the neghorng partcles. Here, the neghor of a partcle s defned as all partcles wth dstance less than or equal to certan numer. The last, the ce partcles ndex that elong to the same oject wll converge to the maxmum ndex of the ce partcles n that ojects. 7. RESULTS AND DISCUSSION 7.1 Physcal setup of the ce meltng The physcal setup of the system n 2D can e descred as follow. Consder a water tank and two ce cues wth the dmenson of each ce cue s m x m. The dmenson of water tank s 1 m x 1 m, and at frst flled wth water of heght 0.25 m as seen n Fgure 2. The ntal temperature of the frst and the second ce cue are C and -10 0C, respectvely, and the temperature of water s 30 0 C. The outsde oundares of the system are adaatc wth the ntal temperature s C. The heat transfer n the system s assumed to occur only y conducton and convecton. The partcles are used for ths smulaton. All parameters for ths smulaton can e seen n Tale 1. For soldfcaton case, the physcal setup s generally same ut dfferent n ntal and oundary condton. The smulaton result for soldfcaton case may e shown later n other paper. 0 Fgure 2. Physcal setup of the ce meltng Tale 1. Parameter for the smulaton No Parameter Value Unt 1 Specfc heat of ce kj/kg 0 C 2 Specfc heat of water kj/kg 0 C
7 SPH METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM 3 Conductvty of ce 2.216e-3 kj/ms 0 C 4 Conductvty of water 0.566e-3 kj/ms 0 C 5 Meltng pont 0 0C 6 Latent heat 334 kj/kg 7 Densty of ce 917 kg/m 3 8 Densty of water 1000 kg/m Smulaton results of ce meltng At frst, two ce cues fall freely from heght of 0.25 m aove water surface. These two ce cues then nteract wth the water. The forces actng on the water partcles and ce partcles are calculated y usng equaton (1) and equaton (16) respectvely. To get temperature of the partcles, heat transfer n each partcle s determned y the enthalpy equaton (4). Temperature dstruton of oth sold partcles and lqud partcles for several tme steps can e seen n Fgure 3. In that fgure, the color represents the temperature of partcles n certan poston. Each mage n Fgure 3 represents not only the phase change from sold to lqud, ut also the patterns of convecton flow. The convecton patterns are qute nterestng to e explored further. These patterns are caused y dfferences n densty etween new water partcles and old water partcles. The new water partcles mean ce partcles whch have just transformed to water partcles. In ths smulaton, the ce partcles own have fxed densty that are 917 kg/m 3, whle the water partcles have densty n the range of kg/m 3 (weakly compressle flud). However, after the ce partcles change the phase from sold to lqud, the densty of the new water partcles are changed to 1000 kg/m 3. The new water partcles caused densty dfference qute strkng wth the surroundng water partcles. There wll e movement of the partcles accordng to equaton (1) and (2). Fgure 3. Temperature dstruton for several tme steps
8 DEDE TARWIDI Fgure 4. Temperature vs tme of a sngle partcle of the frst and second ce cue The left graphc n Fgure 4 shows the temperature of a sngle partcle that ntally located on the edge of the frst ce cues. The ntal temperature of the partcle s C, then the temperature s gradually changed to 0 0 C n seconds. However, n ths condton, the state s stll n the sold phase. To change the phase of ths partcle from sold to lqud the energy s needed for amount latent heat (heat of fuson). The latent heat per partcle s m 0L = 2.99 kj. After the partcle state s changed nto lqud the temperature wll e nfluenced y the temperature of the surroundng water partcles. Total tme requred for a sngle partcle n frst ce cue to change the phase from sold to lqud s aout 2.06 seconds. Whle the rght graphc of Fgure 4 shows the temperature of a sngle partcle of the second ce cue that has the ntal temperature of C. For a sngle partcle n the second ce cue, the total tme needed to change the phase s aout 1.96 seconds. Total tme that requred y an ce partcle to melt depends on the mass and latent heat of the partcle. 6. CONCLUSION In ths paper we have shown that an SPH method can e used to smulate Stefan prolem especally n ce meltng prolem. The advantage of the SPH method s all physcal nformaton such as densty, pressure, velocty, and temperature can e stored n an SPH partcle. Wth ths way, t s easy to treat the nteracton etween sold and lqud phase whch s very hard f on the contrary, the mesh-ased method s used. By usng enthalpy method va SPH formulaton, heat conducton equaton n sold and lqud regon can e solved wthout needed to know the nterface poston etween sold and lqud. As though, ths model s new n SPH development, so t needs many mprovements to have more realstc smulaton. Acknowledgements The author would lke to thank Prof. Sero Omata for gvng very helpful comments and suggestons. The author also thanks Dr. Masak Kazama for the useful dscussons. REFERENCES [1] Alexades, V. and Solomon, A. D., 1993, Mathematcal Modelng of Meltng and Freezng Processes, Taylor & Francs, Washngton, DC.
9 SPH METHOD FOR TWO-DIMENSIONAL STEFAN PROBLEM [2] Voller, V.R. and Shada, L., 1984, Enthalpy methods for trackng a phase change oundary n two dmensons, Int. Comm. Heat and Mass Transfer, 11, 3, pp [3] Beckett, G. and Mackenze, J.A. and Roertson, M.L., 2001, A Movng Mesh Fnte Element Method for the Soluton of Two-Dmensonal Stefan Prolems, J. Comput. Phys., 168, 2, pp [4] Bonnerot, R. and Jamet, P., 1977, Numercal computaton of the free oundary for the two-dmensonal Stefan prolem y space-tme fnte elements, J. Comput. Phys., 25, 2, pp [5] Salvator, L. and Tos, N., 2009, Stefan Prolem through Extended Fnte Elements: Revew and Further Investgatons, Algorthms, 2, 3, pp [6] Lan, C.W. and Lu, C.C. and Hsu, C.M., 2002, An Adaptve Fnte Volume Method for Incompressle Heat Flow Prolems n Soldfcaton, J. Comput. Phys., 178, 2, pp [7] Voller, V. R. and Swamnathan, C. R. and Thomas, B. G., 1990, Fxed grd technques for phase change prolems: A revew, Int. J. Num. Methods n Engneerng, 30, 4, pp [8] Monaghan, J.J. and Huppert, H. E. and Worster, M.G., 2005, Soldfcaton usng smoothed partcle hydrodynamcs, J. Comput. Phys., 206, 2, pp [9] Cleary, P.W., 1998, Modellng confned mult-materal heat and mass flows usng SPH, Appled Mathematcal Modellng, 22, 12, pp [10] Monaghan, J.J. and Kos, A., 1999, Soltary waves on Cretan each, J. Waterway, Port, Coastal, Ocean Eng, 125, 3, pp [11] Cleary, P.W. and Monaghan, J.J., 1999, Conducton Modellng Usng Smoothed Partcle Hydrodynamcs, J. Comput. Phys., 148, 1, pp [12] Iwasak, K. and Uchda, H. and Doash, Y. and Nshta, T., 2010, Fast Partcle-ased Vsual Smulaton of Ice Meltng, Computer Graphcs Forum, 29, 7, pp [13] Tong, M. and Brown, D.J., Smoothed partcle hydrodynamcs modellng of the flud flow and heat transfer n the weld pool durng laser spot weldng, OP Conf. Ser.: Mater. Sc., 27. [14] Lu, M. and Lu, G., 2010, Smoothed Partcle Hydrodynamcs (SPH): an Overvew and Recent Developments, Archves of Computatonal Methods n Engneerng, 17, 1, pp [15] Monaghan, J.J., 2005, Smoothed partcle hydrodynamcs, Rep. Prog. Physc., 68, 8, pp
arxiv: v1 [physics.flu-dyn] 16 Sep 2013
Three-Dmensonal Smoothed Partcle Hydrodynamcs Method for Smulatng Free Surface Flows Rzal Dw Prayogo a,b, Chrstan Fredy Naa a a Faculty of Mathematcs and Natural Scences, Insttut Teknolog Bandung, Jl.
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationComputational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH
Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationCONTROLLED FLOW SIMULATION USING SPH METHOD
HERI COADA AIR FORCE ACADEMY ROMAIA ITERATIOAL COFERECE of SCIETIFIC PAPER AFASES 01 Brasov, 4-6 May 01 GEERAL M.R. STEFAIK ARMED FORCES ACADEMY SLOVAK REPUBLIC COTROLLED FLOW SIMULATIO USIG SPH METHOD
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More information1-Dimensional Advection-Diffusion Finite Difference Model Due to a Flow under Propagating Solitary Wave
014 4th Internatonal Conference on Future nvronment and nergy IPCB vol.61 (014) (014) IACSIT Press, Sngapore I: 10.776/IPCB. 014. V61. 6 1-mensonal Advecton-ffuson Fnte fference Model ue to a Flow under
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationFirst Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.
Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act
More informationNUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD
THERMAL SCIENCE: Year 2018, Vol. 22, No. 5, pp. 1955-1962 1955 NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD Introducton by Tomok IZUMI a* and
More informationTHE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES
THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor
More informationTotal solidification time of a liquid phase change material enclosed in cylindrical/spherical containers
Appled Thermal Engneerng 25 (2005) 488 502 www.elsever.com/locate/apthermeng Total soldfcaton tme of a lqud phase change materal enclosed n cylndrcal/sphercal contaners Levent Blr *, Zafer _ Ilken Department
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 05-03 A comparson of numercal models for one-dmensonal Stefan problems E. Javerre 1, C. Vuk, F.J. Vermolen, S. van der Zwaag ISSN 1389-6520 Reports of the Delft Insttute
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
More informationAmplification and Relaxation of Electron Spin Polarization in Semiconductor Devices
Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton
More informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationSmoothed particle hydrodynamics modelling of fluids and solids
Appled and Computatonal Mechancs 1 (2007) 521-530 Smoothed partcle hydrodynamcs modellng of fluds and solds L. Lobovský a,, J. Křen a a Department of Mechancs, Faculty of Appled Scences, UWB n Plsen, Unverztní
More informationSupplemental Material: Causal Entropic Forces
Supplemental Materal: Causal Entropc Forces A. D. Wssner-Gross 1, 2, and C. E. Freer 3 1 Insttute for Appled Computatonal Scence, Harvard Unversty, Cambrdge, Massachusetts 02138, USA 2 The Meda Laboratory,
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationA large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method
A large scale tsunam run-up smulaton and numercal evaluaton of flud force durng tsunam by usng a partcle method *Mtsuteru Asa 1), Shoch Tanabe 2) and Masaharu Isshk 3) 1), 2) Department of Cvl Engneerng,
More informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationThermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationLinear Momentum. Center of Mass.
Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html
More informationSo far: simple (planar) geometries
Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationPhysics 111: Mechanics Lecture 11
Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationComputational Astrophysics
Computatonal Astrophyscs Solvng for Gravty Alexander Knebe, Unversdad Autonoma de Madrd Computatonal Astrophyscs Solvng for Gravty the equatons full set of equatons collsonless matter (e.g. dark matter
More informationBindel, Spring 2014 Applications of Parallel Computers (CS 5220)
Bndel, Sprng 2014 Applcatons of Parallel Computers (CS 5220) 1 Dervng SPH The Naver-Stokes equatons wth gravty are ρa = p + µ 2 v + ρg. The acceleraton s the materal dervatve of velocty, and we usually
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationPhysics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1
Physcs 141. Lecture 14. Frank L. H. Wolfs Department of Physcs and Astronomy, Unversty of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Course Informaton: Lab report # 3. Exam # 2. Mult-Partcle
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationRichard Socher, Henning Peters Elements of Statistical Learning I E[X] = arg min. E[(X b) 2 ]
1 Prolem (10P) Show that f X s a random varale, then E[X] = arg mn E[(X ) 2 ] Thus a good predcton for X s E[X] f the squared dfference s used as the metrc. The followng rules are used n the proof: 1.
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationConservation of Angular Momentum = "Spin"
Page 1 of 6 Conservaton of Angular Momentum = "Spn" We can assgn a drecton to the angular velocty: drecton of = drecton of axs + rght hand rule (wth rght hand, curl fngers n drecton of rotaton, thumb ponts
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationHow Differential Equations Arise. Newton s Second Law of Motion
page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationPhysics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationConsideration of 2D Unsteady Boundary Layer Over Oscillating Flat Plate
Proceedngs of the th WSEAS Internatonal Conference on Flud Mechancs and Aerodynamcs, Elounda, Greece, August -, (pp-) Consderaton of D Unsteady Boundary Layer Over Oscllatng Flat Plate N.M. NOURI, H.R.
More informationACTM State Calculus Competition Saturday April 30, 2011
ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward
More informationChapter 8. Potential Energy and Conservation of Energy
Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationEN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics
N40: ynamcs and Vbratons Homewor 7: Rgd Body Knematcs School of ngneerng Brown Unversty 1. In the fgure below, bar AB rotates counterclocwse at 4 rad/s. What are the angular veloctes of bars BC and C?.
More informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More informationMechanics Physics 151
Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More informationNumerical Solution of Boussinesq Equations as a Model of Interfacial-wave Propagation
BULLETIN of the Malaysan Mathematcal Scences Socety http://math.usm.my/bulletn Bull. Malays. Math. Sc. Soc. (2) 28(2) (2005), 163 172 Numercal Soluton of Boussnesq Equatons as a Model of Interfacal-wave
More informationRotational Dynamics. Physics 1425 Lecture 19. Michael Fowler, UVa
Rotatonal Dynamcs Physcs 1425 Lecture 19 Mchael Fowler, UVa Rotatonal Dynamcs Newton s Frst Law: a rotatng body wll contnue to rotate at constant angular velocty as long as there s no torque actng on t.
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15
NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound
More informationAngular Momentum and Fixed Axis Rotation. 8.01t Nov 10, 2004
Angular Momentum and Fxed Axs Rotaton 8.01t Nov 10, 2004 Dynamcs: Translatonal and Rotatonal Moton Translatonal Dynamcs Total Force Torque Angular Momentum about Dynamcs of Rotaton F ext Momentum of a
More informationTurbulent Flow. Turbulent Flow
http://www.youtube.com/watch?v=xoll2kedog&feature=related http://br.youtube.com/watch?v=7kkftgx2any http://br.youtube.com/watch?v=vqhxihpvcvu 1. Caothc fluctuatons wth a wde range of frequences and
More informationAhmad Shakibaeinia Assistant Professor Department of Civil, Geological & Mining Engineering Polytechnique Montreal
Natonal Center for Atmospherc Research (NCAR) IMAGe TOY2017: Workshop on Multscale Geoscence Numercs Ahmad Shakbaena Assstant Professor Department of Cvl, Geologcal & Mnng Engneerng Polytechnque Montreal
More informationA PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS
HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationThe Quadratic Trigonometric Bézier Curve with Single Shape Parameter
J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationOptimal Control of Temperature in Fluid Flow
Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationNumerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal
More informationGrand canonical Monte Carlo simulations of bulk electrolytes and calcium channels
Grand canoncal Monte Carlo smulatons of bulk electrolytes and calcum channels Thess of Ph.D. dssertaton Prepared by: Attla Malascs M.Sc. n Chemstry Supervsor: Dr. Dezső Boda Unversty of Pannona Insttute
More information