Acceleration of Gas Bubble-Free Surface Interaction Computation Using Basis Preconditioners
|
|
- William Stokes
- 5 years ago
- Views:
Transcription
1 Acceleraton of Gas Bubble-Free Surface Interacton Computaton Usng Bass Precondtoners K. L. Tan, B. C. Khoo and J. K. Whte ABSTRACT The computaton of gas bubble-free surface nteracton entals a tme-steppng algorthm whereby a lnear system s solved at each tme-teraton. In our nvestgaton, the lnear systems are derved from a desngularzed boundary ntegral formulaton and are poorly condtoned. Ths leads to poor convergence rates when Krylov subspace methods are used to solve these systems. The convergence rates may however be mproved wth proper precondtonng. We lmt our nvestgaton to gas bubbles ntated at depths suffcently small such that a spke forms on the free surface durng the later stages of evoluton. Bubble dynamcs dctate that for gas bubbles ntated at such depths, the stages through whch the gas bubble and free surface evolve are smlar. Based on ths fact, we propose to perform one computaton run for a gas bubble ntated at one partcular depth, obtan a judcous set of a pror bass precondtoners from ths run and thereafter, use ths set of precondtoners on computaton runs for gas bubble ntated at dfferent depths. The computaton tme taken by the proposed method s, n general, 50% and 0% of the tme taken by the present method (wthout precondtonng) wth termnatng crtera of 1.0e-5 and 1.0e-7 n the nfnty-norm respectvely usng the B-conjugate Gradent Stablzed solver. The present method further enables computaton to an nfnty-norm termnatng crteron of 1.0e-10 n a shorter tme compared to the present method wth a crteron of 1.0e-5. K. L. Tan s a student wth the Sngapore-MIT Allance (SMA), Natonal Unversty of Sngapore, 10 Kent Rdge Cresent, Sngapore E-mal: smap9055@nus.edu.sg. B. C. Khoo and J. K. Whte are SMA fellows. INTRODUCTION Bubbles surrounded by a lqud can be found n many applcatons. Examples nclude bubble jet prnters, bolng water and underwater explosons. The chemcal ndustry uses bubble columns to enhance chemcal reactons and mxng. The study of bubbles can be dvded nto two man areas. The frst area, bubbles that do not undergo much change n volume, s studed extensvely. A seemngly smple problem of ths knd, the termnal rse velocty of a sngle bubble n water, s stll not yet fully understood. The second area entals bubble dynamcs pertanng to oscllatng bubbles. Examples of these nclude cavtaton bubbles and underwater exploson bubbles. Underwater exploson bubbles s an actve area of defense research. An exploson or gas bubble s formed when an underwater exploson takes place. Due to nerta, ths bubble wll overexpand and thereafter collapse. If ths collapse takes place near a sold surface, a hgh-speed jet wll be formed, drected towards the sold surface. On the other hand, f t takes place near a free water surface, the hghspeed jet wll then be drected away from the free surface and a water plume or spke wll be formed on the surface. All these problems have been studed usng the boundary ntegral method. Ths method necesstates the meshng of the surface only and not the entre doman thereby greatly reducng the sze of matrces to be solved. In ths paper, we seek to reduce the computaton tme taken to evolve a sngle gas bubble nteractng wth a free surface through use of precondtonng. Lmtng our nvestgaton to gas bubbles ntated at depths such that a spke forms on the free surface, the precondtoners are
2 chosen and computed based on some a pror knowledge of evoluton of such bubbles ntated at dfferent depths. MATHEMATICAL FORMULATION Mathematcal Formulaton Assumng the flud doman bounded by the gas bubble and free surface to be nvscd, ncompressble and rrotatonal, the flow wthn the doman Ω may then be descrbed by a potental feld φ ( x,y,z), where the flud velocty s gven by φ, that satsfes Laplace s equaton φ = 0. The knematc and dynamc boundary condtons governng the moton of the gas bubble wth no buoyancy force, and free surface are Dp = φ, Dt Tme Dscretzaton The Euler Forward tme dscretzaton scheme s used. The tme dscretzed equatons would thus be of the form 1 k k p = p + φ t for the update of the poston of the feld ponts, and and λ 1 k 1 k V0 φ = φ φ ε t k V 1 k 1 k φ = φ φ t for Ω S, for Ω Σ for the update of potental at the feld ponts. smplfcaton, we shall wrte the above two equatons as ( ) 1 k k k φ = φ + F φ,v t For and λ V0 = 1+ φ ε Dφ 1 Dt V for Ω S, Dφ 1 = 1 + φ for Ω Σ Dt where = ( x,y,z) p are the spatal coordnates of the feld pont n 3-D, t s the dmensonless tme varable, γ the dmensonless depth of ntaton, ε the dmensonless strength parameter, V the volume of the gas bubble, V 0 the ntal volume of the gas bubble, λ the rato of specfc heats, Ω the doman boundary, S the porton of Ω defned by the gas bubble and Σ the porton of Ω defned by the free surface (Best [1] and Wang et al. []). The far feld boundary condton s prescrbed as ( ) 0 φ p as p = x + y + z. The value of parameters used n our nvestgaton s ε = 100, for exploson bubble, and γ = 1. 4, for datomc gas (Best [1]). wth F representng the approprate functon for Ω S or Ω Σ. k φ s updated by solvng Laplace s equaton wth boundary condtons gven by φ k. In our code, φ s computed drectly through the ntegral equaton by dfferentatng the Greens functon n the three prncpal drectons. An alternatve method for computng φ would be to obtan the normal component of φ usng the ntegral equaton and the tangental components usng fnte dfferences (Wang et al. []). Soluton of Laplace s Equaton The Desngularsed Indrect Boundary Integral Method (DIBIM) s used to solve Laplace s equaton (Zhang et al. [3]) n order to derve the updates to φ. Based on the above method, we have ( p) φ = σ ( q ) p q where σ s the unknown source strength, = ( x,y,z) q are the spatal coordnates of the source pont, and p q s the
3 Eucldean dstance between the feld pont and the source pont. Intal Condtons The ntal bubble radus for the chosen value of ε s , obtaned by solvng the Raylegh equaton modfed for gas bubble (Best [1]). Assumng an ntal quescent free surface and notng the ntal velocty of the gas bubble surface to be neglgble, the ntal condton of φ = 0 may be prescrbed for the boundary. DISCRETIZATION The ntal geometry of the gas bubble s a dscretzed sphere of 49 nodes wth radus located at ( x,y,z ) = ( 00,,γ ). The ntal quescent free surface s dscretzed as a flat rectangular surface comprsng 79 nodes spannng ( x,y,z ) = ( 5:5, 5:5, 0). The applcaton of a truncated dscretzed free surface mesh s an approxmaton made possble by the far feld boundary condton prescrbed above. A pont collaton scheme s used wth the DIBIM descrbed above. The collocaton s performed on the nodes, 49 for the gas bubble and 79 for the free surface, defned by the dscretzaton. In order to have the number of unknowns equal to the number of known equatons, we set the number of source ponts to be equal to the number of dscretzed nodes, wth 49 ponts assocated wth the gas bubble and 79 ponts assocated wth the free surface. Ths results n a system wth N = 11 degrees-of-freedom. The source ponts are postoned n accordance wth Zhang et al. [3]. As mentoned above, the Euler Forward tme-steppng scheme s used to evolve the gas bubble and free surface. The tme-step s chosen such that the bound on the maxmum change n potental φ for any node s restrcted to 0.04 at each tme-teraton, φ 1 V0 max 1 + φ + ε V t = mn φ 1 max 1 + φ j λ, The value of 0.04 s typcal of the range of φ chosen for such computatons (Best [1] and Wang et al. []). The expected nodal potental values would range from 1 to 3. PRESENT IMPLEMENTATION φ Let = ( x, y, z) p be the poston of feld pont j, j N! be the vector of potentals at the N feld ponts, σ ponts and K be the number of tme-teratons to be taken. N! be the vector of source strengths at the N source The present mplementaton s as follows: 0 Set ntal condtons: p = p j = 1: N, φ = 0 where p j j s the ntal poston of the feld pont j. Do k = 1: K 1. Form matrx G where G j = p q. Solve Gσ = φ for σ usng teratve solver wth set termnatng crteron 3. Compute materal dervatve φ =Hσ where H = j p q j 4. Compute tme-step t 5. Update poston of feld ponts ( ) = + t j = 1: N j j j 1 k p p φ 6. Update potental of feld ponts ( ) φ = φ + F φ, z, V t j = 1: N 1 k j j j where F s the knematc boundary condton functon assocated wth node j 7. Compute new volume V End The teratve solver used s the B-conjugate Gradent Stablzed (BICGSTAB) solver. The G matrces generated j 0 j
4 based on the DIBIM above have been found to be poorly condtoned. A termnatng crteron of r < 1.0e 5 N where r! s the vector of resduals, s used. Ths lberal termnatng crteron s chosen because of the poor convergence characterstcs assocated wth usng a Krylov Subspace solver on ll-condtoned matrces. The DIBIM s also lmted n that t starts to fal as opposte sectons of the bubble surface approach each other. In our nvestgaton, we set K = 150 so as to end our computaton before ths scenaro arses. (c) The output based on the above mplementaton for the evoluton of free surface nteracton wth a sngle gas bubble ntated at γ = 1.0 s shown n Fgure 1. (d) (a) (e) (b) Fgure 1: Evoluton of Free Surface Interacton wth Gas Bubble Intated at 1.0 γ = PROPOSED METHOD The bound on convergence rate of a Krylov Subspace solver for the soluton of a general lnear system Ax = b
5 s a functon of the condton number of the matrx A, κ ( A ). In general, the better the condton of A, the hgher s the convergence rate and vce-versa. Block Dagonal Precondtonng The use of precondtonng to acheve better convergence propertes for teratve solvers s not new. There are also many forms of precondtoners (Saad [5]). In our nvestgaton, we choose to use block dagonal precondtoners. The ratonale for such a choce s explaned below. As stated n Step 1 of the mplementaton shown before, the entres n the G matrx are G j = p q. Ths means that the magntude of entry G s nversely j proportonal to the dstance between feld pont j and source pont. Let the 49 nodes assocated wth the gas bubble be ndexed j = 1: 49 and the 79 nodes assocated wth the free surface be ndexed j = 493 :11. The resultng G matrx generated would be of the form G1 G3 G = G4 G where G1! s the sub-matrx assocated wth the gas bubble nodes only, G! s the sub-matrx assocated wth the free surface nodes only, and G3! and G4! are the cross-coupled terms. j ponts n accordance wth Zhang et al. [3], t s expected that the value of all entres n matrces G3 and G4 to be much smaller n magntude as compared to the value of the entres n matrces G1 and G. As the gas bubble evolves, t s observed that only the top secton of the bubble surface comes nto close proxmty wth the free surface. Ths mples that most of the couplng occurs wthn matrces G1 and G, and not wthn G3 and G4. Based on the above observaton, f we are able to construct a precondtoner M for the G generated at every tme-teraton such that G1 0 M = 0 G then (rght-)precondtonng the matrx G wth M would yeld I G3 G GM =, G4 G1 I resultng n a better condtoned matrx. It would however be computatonally expensve to nvert sub-matrces G1 and G at every tme-teraton for use as precondtoners. It s observed from Fgure 1a that durng the ntal stages of evoluton, the dstance between two bubble nodes s much smaller compared to the dstance between a bubble node and a free surface node. Because of the postonng of the source
6 The need to compute precondtoners at every tmeteraton can be elmnated f we relax the requrement that the block precondtoner sub-matrces be exact nverses of sub-matrces G1 and G. In fact, our am s just to precondton the G matrx such that t has better convergence propertes and not to convert the dagonal submatrces nto dentty sub-matrces. Fgures a, b and c show the gas bubble and free surface for γ = 1.0 at tmeteratons 60, 65 and 70 respectvely. It s observed from Fgure that the change n geometry over several tmeteratons s lmted. Ths mples that the change n value of entres n the G matrx from one tme-teraton to the next s also lmted. In ths case, we can compute the precondtoner for the G matrx generated at a partcular tme-teraton and make use of ths precondtoner to precondton the G matrces generated at tme-teratons about that partcular tme-teraton from whch the precondtoner s obtaned. In ths case, what we obtan s (a) G1 0 M = 0 G and after (rght-)precondtonng (b) (c) Fgure : Gas Bubble and Free Surface Geometry at Tmeteratons 60 (a), 65 (b) and 70 (c) for 1.0 γ =. I G3 G GM = G4 G1 I whch would stll result n a relatvely well-condtoned matrx. By judcously choosng a handful of G matrces throughout the entre computaton run to generate the precondtoners from, effectve precondtonng may be acheved for the entre computaton run. The above method may mprove the condton of the matrx and therefore reduce the tme taken to solve the lnear system at each tme-teraton. However, ths method would requre a pror knowledge of the G matrces so as to compute the precondtoners. Ths would mply that, n order to accelerate the computaton, we would have to perform the sad computaton frst, whch s somewhat neffcent. An alternatve method s to decde a pror at whch tme-
7 teraton to compute the precondtoners. Ths would however necesstate the nverson of matrces durng the computaton run thus makng the computaton run expensve. In fact, there may not even be any savngs n computaton tme n usng ths method over the present method. A Pror Computed Bass Precondtoners Gven n Fgure 3 s the evoluton of free surface nteracton wth a sngle gas bubble ntated at γ = 1.1. Comparng Fgures 3 and 1, t s observed that the gas bubble n both cases, though ntated at dfferent depths, goes through the same sequence of evoluton wth small dfferences n geometry. (d) (a) (b) (e) Fgure 3: Evoluton of Free Surface Interacton wth Gas Bubble Intated at 1.1 γ = As mentoned before, our am of precondtonng s not to generate dentty sub-matrces but smply to mprove the condton of the G matrces generated such that they have better convergence propertes. It s therefore plausble to make use of precondtoners generated from the computaton run of γ = 1.0 for the computaton run of γ = 1.1. Ths clam s substantated by the observaton above that for a sngle gas bubble ntated at varous depths such that a spke forms on the free surface, the gas bubble and free surface would go through the same sequence of evoluton wth margnal dfferences n geometry In ths case, the need to have a pror knowledge of the G matrces for each new computaton run of a dfferent γ s elmnated. In fact, we need only perform one computaton run of one value of γ usng the present method, compute the precondtoners from the G matrces obtaned at judcously chosen tme-teratons and then make use of these precondtoners for other computaton runs of any value of γ that results n spke formaton on the free surface. (c)
8 Proposed Implementaton The proposed mplementaton s gven below. Perform only once to obtan the bass precondtoners: Judcously select the tme-teratons to compute the bass precondtoners from. Let the number of bass precondtoners selected be n#. For a selected value of γ = γ#, perform computaton run usng the present method. If tme-teraton k = one of the n# pre-selected tme-teratons above, 1. Extract sub-matrces G1 and G, and compute and G. Form M 3. Wrte M to fle G1 At the end of the computaton run, there should be n# a pror computed bass precondtoners. To obtan the soluton to any value of γ : Set ntal condtons. Do k = 1: K 1. Form matrx G. Precondton G wth the most approprate M from the set of n# bass precondtoners 3. Solve Gσ = φ for σ usng teratve solver wth set termnatng crteron 4. Compute materal dervatve φ =Hσ 5. Compute tme-step t 6. Update poston of feld pont ( ) = + t j = 1: N j j j 1 k p p φ 7. Update potental of feld ponts 1 k φ φ φ ( ) = + F, V t j = 1: N j j where F s the knematc boundary condton functon assocated wth node j 8. Compute new volume V End In Step above, the most approprate M would be the precondtoner whose geometry whch t s computed from, most closely matches the geometry at tme-teraton k, whch matrx G s computed from. Rght-precondtonng s used snce ths preserves the resdual of the orgnal lnear system. Selecton of Precondtoners In the selecton of precondtoners, we would not want to select too large a number snce ths would degrade the performance of the proposed method through excessve readng of fles. On the other hand, we would not want to select too low a number of precondtoners snce ths would not effcently precondton the G matrces and also degrade the performance of the proposed method. For the computaton of the bass precondtoners, we use γ = 1.0. Ths value of γ s used because t s approxmately the mean of the range of γ values that we are nterested n nvestgatng. The selecton crteron we choose s based on the amount of change n geometry wthn a range of tme-teratons, rather than smply a unform dstrbuton based on the number of tme-teratons or geometry of bubble. The ratonale for the above crteron s that f the geometry changes mnmally over a wde range of tme-teratons, there s no need for many precondtoners over that range of tme-teratons on the bass of the lmted change n geometry. On the other hand, f the change n geometry s tremendous over a lmted range of tme-teratons, there s no need for many precondtoners over that range of geometry on the bass of the lmted number of tme-teratons taken to evolve over that range of geometry. Snce the volume of the gas bubble s solely a functon of the geometry of the bubble, we can replace geometry n the above crteron wth volume of the gas bubble. We can also use the volume to establsh the most approprate precondtoner to use. The alternatve vew s that we can determne the range of volumes that we use a partcular precondtoner for. Based on the above crteron, we look to the volumeteraton plot gven n Fgure 4 to select the volume, and thus the tme-teraton to compute the precondtoners from.
9 It s observed from Fgure 4 that the ncrease n volume durng the frst 40 tme-teratons s mnmal. The rate of ncrease pcks up dramatcally at about tme-teraton 60, taperng off upon nearng maxmum volume. The rate of change n volume durng the collapse phase s slghtly less than that durng the rapd expanson phase. Based on the selecton crteron detaled above, the volumes and the correspondng tme-teratons we choose to compute the precondtoners for varous ranges of volume are gven n Table 1. volume tme-teraton Fgure 4: Gas Bubble Volume aganst Tme-teraton Plot for 1.0 γ =. Desred Volume Selected Tmeteraton Range of Volume Selected For Intal Table 1: Volumes Chosen for Computaton of Precondtoners RESULTS AND DISCUSSION The cpu-tme (CPU) and wall-clock (WC) tme taken to complete a computaton run of 150 tme-teratons for varous γ s s gven below. The teratve soluton method used s the BICGSTAB method wth the termnatng crteron beng the nfnty norm of the resdual as s wth the present method. Tables, 3 and 4 show the tmes taken wth termnatng crtera of 1.0e-5, 1.0e-7 and 1.0e-10 respectvely. Table 5 shows a comparson of the total number of teratons taken by the teratve solver durng the entre computaton run of 150 tme-teratons. It s observed from the results above that wth a termnatng crteron of 1.0e-5, the proposed method reduces the tme taken by almost 50% over the present method. Upon tghtenng the termnatng crteron to 1.0e-7, the proposed method reduces the tme taken by almost 80% over the present method. There s no tmng recorded for the present method wth a termnatng crteron of 1e-10. However, an alternatve way of nterpretng the results would be that wth the proposed method, we can acheve an mprovement of 5 orders n the nfnty norm of resdual over the present method whle takng less tme to complete the computaton run. Further evdence of the postve effect of the proposed precondtonng method on the convergence characterstcs of the G matrces formed may be observed n Table 5, whch shows the total number of teratons taken by the teratve solver for an entre computaton run of 150 tme-teratons at varous γ s. Fgure 5 shows the plots of number of teratons taken by the BICGSTAB solver at each tme-teraton for 0.8 γ = usng the present method and the proposed method wth a termnatng crteron of 1.0e-7. It s observed from Fgure 5 that the convergence rate s mproved throughout the entre computaton run. Havng determned the range of volume assocated wth each precondtoner, a new precondtoner s read nto memory from fle as and when necessary durng the computaton.
10 Proposed Method Present Method Rato CPU (s) WC (mn:s) CPU (s) WC (mn:s) CPU WC γ = : : γ = : : γ = : : γ = : : γ = : : Table : Tme taken for 150 Tme-teratons wth Iteratve Solver Termnatng Crteron of 1.0e-5 Proposed Method Present Method Rato CPU (s) WC (mn:s) CPU (s) WC (mn:s) CPU WC γ = : : γ = : : γ = : : γ = : : γ = : : Table 3: Tme taken for 150 Tme-teratons wth Iteratve Solver Termnatng Crteron of 1.0e-7 Proposed Method CPU (s) WC (mn:s) γ = :08 γ = : γ = :55 γ = :56 γ = :09 Table 4: Tme taken for 150 Tme-teratons wth Iteratve Solver Termnatng Crteron of 1.0e-10 Tolerance 1.0e-5 1.0e-7 1.0e-10 Method Proposed Present Rato Proposed Present Rato Proposed γ = γ = γ = γ = γ = Table 5: Total Number of Iteratons taken by Iteratve Solver for 150 Tme-teratons
11 number of BICGSTAB teratons present method proposed method 3. Zhang, Y. L., Yeo, K. S., Khoo, B. C. and Chong, W. K., 1999, Smulaton of Three-dmensonal Bubbles Usng Desngluarsed Boundary Integral Method, Internatonal Journal for Numercal Methods n Fluds Vol. 31, No. 8, pp Zhang, Y. L., Yeo, K. S., Khoo, B. C. and Chong, W. K., 1998, Three-dmensonal Computaton of Bubbles Near a Free Surface, Journal of Computatonal Physcs 146, pp Saad, Y., 1996, Iteratve Methods for Sparse Lnear Systems, PWS Publshng Company tme-teraton Fgure 5: Comparson of Number of Iteratons taken by Iteratve Solver at each Tme-teraton between Present Method and Proposed Method CONCLUSIONS Based on some a pror knowledge of the evoluton of free surface nteractng wth a sngle gas bubble ntated at the range of depths under nvestgaton, a set of bass precondtoners s computed from the computaton run made for one partcular depth. Ths set of bass precondtoners s then used to precondton the matrces generated at each tmeteraton for computaton runs made for dfferent depths. Usng ths proposed method, the tme taken to complete 150 tme-teratons of evoluton of free surface nteractng wth a sngle bubble s 50% and 0% of the tme taken by the present method wth termnatng crtera of 1.0e-5 and 1.0e-7 respectvely for the teratve solver used. Furthermore, the proposed method allows computaton to a tolerance of 1.0e- 10 n a tme less than that taken by the present method to compute to a tolerance of 1.0e-5. REFERENCES 1. Best, J. P., 1991, The Dynamcs of Underwater Explosons, Ph.D. Thess, The Unversty of Wollongong, Australa.. Wang, Q. X., Yeo, K. S., Khoo, B. C. and Lam, K. Y., 1996, Nonlnear Interacton Between Gas Bubble and Free Surface, Computers & Fluds Vol. 5, No. 7, pp
Chapter 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 informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
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 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 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 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 informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
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 informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationVector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.
Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm
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 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 informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationRelaxation Methods for Iterative Solution to Linear Systems of Equations
Relaxaton Methods for Iteratve Soluton to Lnear Systems of Equatons Gerald Recktenwald Portland State Unversty Mechancal Engneerng Department gerry@pdx.edu Overvew Techncal topcs Basc Concepts Statonary
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 informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
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 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 informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationYong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )
Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often
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 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 informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
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 informationx = , so that calculated
Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to
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 informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationDFT with Planewaves pseudopotential accuracy (LDA, PBE) Fast time to solution 1 step in minutes (not hours!!!) to be useful for MD
LLNL-PRES-673679 Ths work was performed under the auspces of the U.S. Department of Energy by under contract DE-AC52-07NA27344. Lawrence Lvermore Natonal Securty, LLC Sequoa, IBM BGQ, 1,572,864 cores O(N)
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 informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
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 informationPsychology 282 Lecture #24 Outline Regression Diagnostics: Outliers
Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.
More information8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS
SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationInvestigation of a New Monte Carlo Method for the Transitional Gas Flow
Investgaton of a New Monte Carlo Method for the Transtonal Gas Flow X. Luo and Chr. Day Karlsruhe Insttute of Technology(KIT) Insttute for Techncal Physcs 7602 Karlsruhe Germany Abstract. The Drect Smulaton
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
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 informationAppendix B: Resampling Algorithms
407 Appendx B: Resamplng Algorthms A common problem of all partcle flters s the degeneracy of weghts, whch conssts of the unbounded ncrease of the varance of the mportance weghts ω [ ] of the partcles
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 informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
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 informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
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 informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
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 information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
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 informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationarxiv: 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 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 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 informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
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 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 informationGrover s Algorithm + Quantum Zeno Effect + Vaidman
Grover s Algorthm + Quantum Zeno Effect + Vadman CS 294-2 Bomb 10/12/04 Fall 2004 Lecture 11 Grover s algorthm Recall that Grover s algorthm for searchng over a space of sze wors as follows: consder the
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
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 informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationThe Study of Teaching-learning-based Optimization Algorithm
Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute
More informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
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 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 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 informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationFeb 14: Spatial analysis of data fields
Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationP R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering /
Theory and Applcatons of Pattern Recognton 003, Rob Polkar, Rowan Unversty, Glassboro, NJ Lecture 4 Bayes Classfcaton Rule Dept. of Electrcal and Computer Engneerng 0909.40.0 / 0909.504.04 Theory & Applcatons
More informationProblem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)
Problem adapted reduced models based on Reacton-Dffuson Manfolds (REDIMs) V Bykov, U Maas Thrty-Second Internatonal Symposum on ombuston, Montreal, anada, 3-8 August, 8 Problem Statement: Smulaton of reactng
More informationLimited Dependent Variables
Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages
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 informationNorms, Condition Numbers, Eigenvalues and Eigenvectors
Norms, Condton Numbers, Egenvalues and Egenvectors 1 Norms A norm s a measure of the sze of a matrx or a vector For vectors the common norms are: N a 2 = ( x 2 1/2 the Eucldean Norm (1a b 1 = =1 N x (1b
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
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 information[7] R.S. Varga, Matrix Iterative Analysis, Prentice-Hall, Englewood Clis, New Jersey, (1962).
[7] R.S. Varga, Matrx Iteratve Analyss, Prentce-Hall, Englewood ls, New Jersey, (962). [8] J. Zhang, Multgrd soluton of the convecton-duson equaton wth large Reynolds number, n Prelmnary Proceedngs of
More informationLOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin
Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
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 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 informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
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 informationSuppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl
RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com
More informationAdjoint Methods of Sensitivity Analysis for Lyapunov Equation. Boping Wang 1, Kun Yan 2. University of Technology, Dalian , P. R.
th World Congress on Structural and Multdscplnary Optmsaton 7 th - th, June 5, Sydney Australa Adjont Methods of Senstvty Analyss for Lyapunov Equaton Bopng Wang, Kun Yan Department of Mechancal and Aerospace
More informationρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to
THE INVERSE POWER METHOD (or INVERSE ITERATION) -- applcaton of the Power method to A some fxed constant ρ (whch s called a shft), x λ ρ If the egenpars of A are { ( λ, x ) } ( ), or (more usually) to,
More information