Application of Finite Element Method (FEM) Instruction to Graduate Courses in Biological and Agricultural Engineering
|
|
- Kerry Spencer
- 5 years ago
- Views:
Transcription
1 Sesson: 08 Applcaton of Fnte Element Method (FEM) Instructon to Graduate Courses n Bologcal and Agrcultural Engneerng Chang S. Km, Terry H. Walker, Caye M. Drapcho Dept. of Bologcal and Agrcultural Engneerng Lousana State Unversty, Baton Rouge Abstract The applcaton of Fnte Element Methods (FEM) to a graduate level course n Bologcal Engneerng, Advanced Transport Phenomena n Bologcal Engneerng, s presented. Frst, the Galerkn Weak Statement (GWS) was ntroduced to the class to show the fundamental theory of FEM by solvng a D steady state heat transfer problem. Ths technque provdes a more accurate soluton wth the estmaton of error. The concept of error reducton through mesh refnement was also ntroduced. Each student was requred to conduct an ndependent semester project ncorporatng mathematcal modelng and smulaton of a bologcal engneerng problem. One of these projects, fxed bed on exchange modelng, s dscussed n ths paper. The outputs from these class projects llustrate that the students ganed experence n usng FEM to solve dynamc bologcal engneerng problems. Introducton Computer aded modelng of new products has allowed ndustry to quckly optmze desgn whle spendng less tme and money on physcal prototypes. Boprocess and food process engneers often deal wth complex heterogeneous system characterzed by non-newtonan behavor. Solutons to partal dfferental equatons that descrbe these complex systems are dffcult to obtan. Advantages of usng Computer-Aded Engneerng (CAE) prototypng n food and boprocess development (Datta 998; Baker et al. 999) and applcaton to mechancs of materals (Hllsman 994) have been prevously addressed They nclude: ) quck and nexpensve testng of alternatve scenaros that can result n reduced costs and ncreased profts, ) clear understandng of the nteractons between the physcal processes and ther senstvty to varous operatonal parameters, and 3) front-end engneerng before prototypng, makng the prototypes closer to optmum and reducng ther number. Development of computer models to descrbe these complex boprocessng systems s needed. Ths paper addresses the experence n ntroducng Fnte Element Methods (FEM) to a graduate level course n Bologcal Engneerng n LSU, BE 735, Advanced Transport Phenomena n Bologcal Engneerng, durng the fall semester of 00. The course ncluded Page 7..
2 two aspects of FEM development. Frst, the fundamentals of FEM usng the Galerkn Weak Statement (GWS) were ntroduced to the class wth dscusson of error analyss. Ths technque provdes a more accurate soluton wth the estmaton of error. The concept of error reducton through mesh refnement was also ntroduced. An example one-dmensonal steady state heat transfer problem was used to llustrate these concepts. Secondly, the students ganed frst-hand experence wth FEM by applcaton to ndependent semester projects n bologcal engneerng. Each project was requred to nclude computatonal modelng and smulaton of a bologcal engneerng problem, prmarly from the student s research nterest. The fundamental theory of FEM One dmensonal heat transfer was used to present the fundamental theory of FEM. Smple -D heat transfer through a plane wall s shown n Fgure. In ths problem, the equaton s gven n terms of T (temperature) on the doman S, on the boundary a < x < b wth correspondng Neumann (constant flux, q) and Drchlet (constant temperature) boundary condtons. Governng Equaton d æ dt ö - çk - s = 0 over a < x < b dx è dx ø Boundary Condtons dt - k = q at x = a dx T = Tb at x = b Fgure Conducton through a plane wall wth unform thermal conductvty Page 7..
3 The steps to obtan the fnte element soluton presented to the class are descrbed as follows (Baker and Pepper, 99):. State the governng equaton, and ntal and boundary condtons. Defne the approxmate soluton n terms of known spatal functons multpled by unknown expanson coeffcents. Determne that the sutable spatal functons exst and determne the unknown expanson coeffcents. T N N 3 3 N N å Nf N = ( x) = a ( x) + a f ( x) + a f ( x) + + a f ( x) = a ( x) f K () where T N = approxmate soluton, a s the expanson coeffcent, and N N s the tral spatal functon. The relatonshp between the exact soluton and approxmate soluton s T N N ( x) T ( x) + e ( x) = () where e N (x) s the approxmate soluton error.. Defne the Galerkn Weak Statement (GWS) to determne the expanson coeffcents that ensure an absolute mnmum approxmate error. Ths procedure mnmzes the dstance between the T(x) and T N (x) for any specfed number of terms n the approxmaton, N. To form the GWS, the weghted resdual approach s used to meet the Galerkn crteron that the weght functons must be equal to the selected tral functons to mnmze the error by makng the set of functons orthogonal. N ( x) L( T ) GWS = ò f dx º 0 (3) W where f (x) s defned below and L(T N ) s governng dfferental equaton of the system n terms of the approxmaton varable T N. The GWS for the gven -D heat transfer problem wth lnear tral space functons s: x - x - f ( x ) = for x - x x (4) x - x + - x + - x f ( x ) = for x x x + (5) x - x ( x) = 0 for x > x + or x < x - f (6) 3. Select pecewse contnuous spatal functons (lnear, quadratc, or cubc). 4. Defne a dscretzaton of the soluton spatal doman n terms of fnte elements to Page 7..3
4 evaluate the weak statement ntegrals and apply boundary condton to solve for the unknown expanson coeffcents. 5. Determne the approxmaton accuracy by determnng the error to valdate soluton qualty. The approxmate soluton of a -D heat transfer equaton obtaned through these steps can be obtaned and compared to the analytcal soluton, f the analytcal soluton exsts. Asymptotc Error Estmaton The theoretcal estmate of performance of the lnear bass Fnte Element mplementaton of the GWS s descrbed by the error sem-norm e h E, whch cannot be determned unless the exact soluton s known. Thus, t s not possble to estmate error usng a sngle soluton on mesh doman, W h. In practce, a mesh refnement process s followed usng a sequence of nested meshes, W h, W h/, W h/4, etc. The defnton of approxmaton error for two such solutons yelds T + h h h / h / + e = T = T e (7) Smlarly, the energy sem-norm relaton s gven as: T h E h h / h / + e = T = T + e. (8) E E E E The energy norm of the estmated soluton for -D heat transfer equaton can be obtaned usng followng equaton; h h h h T = E T, T º ò kt dx W. (9) M T = å ({ Q} [ DIFF]{ Q} e ) e e= The relatonshp between these two energy sem-norms s derved as h k h / e = e. (0) Thus, the asymptotc error estmate s defned as h / h / DT e = () k - D T = T - T = - e. h / h / h k h / where ( ) Page 7..4
5 The calculated error estmate for the example D heat transfer problem s tabulated n Table, whch explans the relatonshp between mesh refnement and the mprovement of approxmate soluton. Table The calculated error estmate Number of mesh T h D T h/ e h/ 4.50e e e e e e FEM software The FEM program that was selected was FEMLAB (Comsol, Inc), whch s an nteractve MATLAB-based envronment for modelng and solvng scentfc and engneerng problems based on partal dfferental equatons (PDEs). FEMLAB ntegrates computaton, vsualzaton, and programmng n an easy-to-use envronment. FEMLAB frequently uses MATLAB s syntax and data structures. One beneft of ths ntegraton s that you can save and export FEMLAB models as MATLAB programs that run drectly n that envronment, whch allows the freedom to combne FEM-based modelng, smulaton, analyss wth other engneerng algorthms. Applcaton of FEM to Bologcal Engneerng Problems Student Project An on exchange model wth lnear drvng force was developed to descrbe color removal from a bologcal mxture usng on exchange resns. Color removal n on exchange resns can be modeled as an adsorpton process of a dlute speces. Lttle nteracton between dfferent molecules n the flud s assumed, whch allows the applcaton of sngle component adsorpton models. Materal s ntroduced as a bulk flud at the top of a fxed-bed of sphercal on exchange resn beads. Local equlbrum s assumed around the adsorbent bead,.e. the adsorpton reacton s fast, allowng the amount of materal adsorbed onto the bead surface to be determned by an equlbrum sotherm relatonshp. The temperature of the columns s controlled so sothermal condtons are assumed. A concentraton boundary layer (flm) forms around the bead and wll be accounted for by the classcal lnear drvng force (LDF) approxmaton (Rce 98). The flux of materal to the bead s determned by a mass transfer coeffcent multpled by the dfference n concentraton between the bulk and flm. The beads have a bdsperse pore dstrbuton, beng an agglomeraton of many mcro-beads. It s assumed that the bead structure may be reduced nto a lumped parameter, effectve mass transfer coeffcent (Rce 98). The LDF mass transfer coeffcent now becomes an effectve parameter that s assumed to take nto account the pore dffuson effects. Page 7..5
6 Assumng that the bed porosty s constant and that there s no concentraton gradent n the radal drecton, a one-dmensonal, tme dependent axal dsperson model may be used for the fxed bed: C C C q - D + e + - e (3) ( ) 0 v 0 ax = z z t t where n 0 s velocty, C s the bulk flud concentraton, z s the axal drecton, D ax the dsperson coeffcent n the axal drecton, e s the porosty, and q s the concentraton of the adsorbed speces. Ths partal dfferental equaton (PDE) s coupled to the LDF approxmaton of the flm. q * ( - e) = k a( C - C ) c (4) t where k c s the effectve mass transfer coeffcent, a the nterfacal area, C * s the concentraton n the flm. The parameters k c and a wll be lumped together as k c a, for ths dscusson. The concentraton of the adsorbed speces, q, s related to the concentraton n the flm, C *, by an * sotherm, descrbed n general as q = f ( C ). Ths may be nverted to yeld q and substtuted nto (Equaton 4): q - ( - e) = k ca( C - f ( q) ). (5) t Two hyperbolc PDE s descrbng the adsorpton under the LDF approxmaton have been developed. These two equatons were converted to dmensonless form to reduce the number of parameters nvolved: f f f - + h Pe h q y q + ( f - x ( y) ) - ( - e) C 0 x e y q = 0 (6) = St (7) C z v 0L bulk convecton where f =, C 0 = concentraton at tme 0; h =, Pe = =, q = C 0 L D ax axal dsperson el q e - e q t =, L = total bed length, y =, x = C 0 (.e. y = ), and v 0 x - e e C 0 k c al flm mass transfer St = =. v bulk convecton 0 t, t Page 7..6
7 The model posed was solved usng FEMLAB, a plug-n to the mathematcs applcaton of MATLAB. The problem defned here s a -dmensonal, tme-dependent system wth two dependent varables (f,y) and one ndependent varable (h). The senstvty of the breakthrough curve to a parameter can be analyzed by conductng repeated smulatons. For example, the breakthrough curve s very senstve to changes to Peclet (Pe) number and nsenstve to the Stanton (St) number (Fgures 3 and 4). Fgure 3. Senstvty of adsorpton breakthrough curves to the Peclet number Page 7..7
8 Fgure 4. Senstvty of adsorpton breakthrough curves to the Stanton number Student feedback The evaluaton from the students at the end of fall semester of 00 was conducted concernng nstructon n the applcaton of computatonal fnte element methods usng FEMLAB. The evaluaton questons gven to the student are summarzed at Table. The evaluaton results are shown n Fgure 5. The majorty of the students responded that applyng computatonal methods to bologcal engneerng problems was benefcal and expect to use computatonal FEM n the future for ther other project, labs, employment, and etc. Page 7..8
9 No Somewhat Yes 0 Queston Queston Queston 3 Queston 4 Queston 5 Table Evaluaton questons, whch were asked to answer based on the scale: -no, -somewhat, and 3-yes Practcal Tranng Evaluaton USDA Hgher Educaton Challenge Grant. Dd you lke applyng computatonal fnte element methods (FEM) usng FEMLAB (and MATLAB)?. Dd the applcaton of computatonal FEM beneft you n the learnng of transport phenomena? 3. Do you understand the basc concept of numercal methods or the fnte element method? 4. Dd the nstructor (nstructors) explan the use of FEM and FEMLAB adequately? 5. Do you expect to use computatonal FEM n the future (other projects, labs, employment, etc)? Fgure 5. Student evaluatons. Summary The ntroducton of fundamental theory of FEM and the use of FEM software, FEMLAB, has been a postve step n the ntroducton of numercal technques to the Bologcal Engneerng graduate currculum. Page 7..9
10 Thus far, only lnear bass functons n fnte element analyss have been utlzed. The concept of bass functon wll be expanded to nclude hgher polynomal-degree, quadratc, and cubc fnte element bass functons. The accuracy ssue of these bass functons wll be compared usng error estmate wth energy sem-norm. Acknowledgments Ths work was supported by an USDA Hgher Educaton Challenge Grants program FY 000. The authors acknowledge the students who partcpated n the courses utlzng FEMLAB n ther projects, especally Hugh Broadhurst s contrbutons from LSU Audubon Sugar Insttute on the on exchange model. Reference Baker, A. J. and Pepper, D. W. 99. Fnte Elements --3. McGraw-Hll, New York, NY. Baker, A. J., Chambers, Z., and Taylor, M. B.999. Fnte Element Analyss for the Engneerng Scence: a Web-Based, Vdeo-Streamed Educaton Envronment at a Dstance, 999 ASEE Annual Conference Proceedngs Datta, A. K.998. Computer-Aded Engneerng n Food Processng and Product Desgn, Food Technology, 5(0):44-5. Hllsman, V. S.994. Combnng Fnte Element Analyss Software wth Mechancs of Materals, 994 ASEE Annual Conference Proceedngs Rce, RG. 98. Approxmate solutons for batch, packed tube and radal flow adsorbers comparson wth experment. Chem. Eng. Sc. 37:83-97 Bographcal nformaton Chang-Sk Km receved hs BS n Mechancal Engneerng at Kangwon Natonal Unversty, South Korea, and MS and PhD n Engneerng Scence at Lousana State Unversty. He s currently a Post Doctoral Assocate n the Department of Bologcal & Agrcultural Engneerng at Lousana State Unversty. He has research nterests n computer vson, pattern recognton, artfcal neural network, and bonformatcs. Terry H. Walker receved hs BS n Engneerng Scence and Mechancs, and MS and Ph.D. degrees n Bosystems Engneerng at the Unversty of Tennessee. He s currently an Assstant Professor n the Department of Bologcal and Agrcultural Engneerng at Lousana State Unversty. He has research nterests n boprocessng ncludng fungal fermentatons and supercrtcal flud extracton technologes. Page 7..0
11 Caye M. Drapcho receved her BS and MS degrees The Pennsylvana State Unversty. She receved her Ph. D. n Agrcultural and Bologcal Engneerng at Clemson Unversty n 993. She s currently an Assocate Professor n the Bologcal and Agrcultural Engneerng Department n the LSU AgCenter. Her areas of specalzaton are assessment of agrcultural nonpont source pollutants and bologcal reactor desgn for agrcultural waste treatment. Page 7..
Numerical 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 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 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 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 informationmodeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products
modelng of equlbrum and dynamc mult-component adsorpton n a two-layered fxed bed for purfcaton of hydrogen from methane reformng products Mohammad A. Ebrahm, Mahmood R. G. Arsalan, Shohreh Fatem * Laboratory
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 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 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 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 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 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 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 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 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 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 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 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 informationThe Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites
7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT
More information10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization
10.34 Numercal Methods Appled to Chemcal Engneerng Fall 2015 Homework #3: Systems of Nonlnear Equatons and Optmzaton Problem 1 (30 ponts). A (homogeneous) azeotrope s a composton of a multcomponent mxture
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
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 informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
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 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 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 informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
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 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 informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
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 informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
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 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 informationAnalytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures
IMS 2 Workshop Analytcal Gradent Evaluaton of Cost Functons n General Feld Solvers: A Novel Approach for Optmzaton of Mcrowave Structures P. Harscher, S. Amar* and R. Vahldeck and J. Bornemann* Swss Federal
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 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 informationNeuro-Adaptive Design - I:
Lecture 36 Neuro-Adaptve Desgn - I: A Robustfyng ool for Dynamc Inverson Desgn Dr. Radhakant Padh Asst. Professor Dept. of Aerospace Engneerng Indan Insttute of Scence - Bangalore Motvaton Perfect system
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 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 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 new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
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 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 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 informationPREDICTIVE CONTROL BY DISTRIBUTED PARAMETER SYSTEMS BLOCKSET FOR MATLAB & SIMULINK
PREDICTIVE CONTROL BY DISTRIBUTED PARAMETER SYSTEMS BLOCKSET FOR MATLAB & SIMULINK G. Hulkó, C. Belavý, P. Buček, P. Noga Insttute of automaton, measurement and appled nformatcs, Faculty of Mechancal Engneerng,
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 informationrisk and uncertainty assessment
Optmal forecastng of atmospherc qualty n ndustral regons: rsk and uncertanty assessment Vladmr Penenko Insttute of Computatonal Mathematcs and Mathematcal Geophyscs SD RAS Goal Development of theoretcal
More informationSIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD
SIMUATION OF WAVE POPAGATION IN AN HETEOGENEOUS EASTIC OD ogéro M Saldanha da Gama Unversdade do Estado do o de Janero ua Sào Francsco Xaver 54, sala 5 A 559-9, o de Janero, Brasl e-mal: rsgama@domancombr
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 informationLinear Feature Engineering 11
Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19
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 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 informationModeling of Breakage Equations Using PBE Software
Modelng of Breakage Equatons Usng PBE Software Prsclla Hll 1 Abstract Partcle breakage occurs n a varety of processng unts n chemcal processng rangng from crushers and grnders to breakage of partcles n
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 informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
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 informationPrinciples of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer
Prncples of Food and Boprocess Engneerng (FS 31) Solutons to Example Problems on Heat Transfer 1. We start wth Fourer s law of heat conducton: Q = k A ( T/ x) Rearrangng, we get: Q/A = k ( T/ x) Here,
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 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 informationParameter Estimation for Dynamic System using Unscented Kalman filter
Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,
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 informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
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 informationA LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS. Dr. Derald E. Wentzien, Wesley College, (302) ,
A LINEAR PROGRAM TO COMPARE MULTIPLE GROSS CREDIT LOSS FORECASTS Dr. Derald E. Wentzen, Wesley College, (302) 736-2574, wentzde@wesley.edu ABSTRACT A lnear programmng model s developed and used to compare
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 informationStatistics II Final Exam 26/6/18
Statstcs II Fnal Exam 26/6/18 Academc Year 2017/18 Solutons Exam duraton: 2 h 30 mn 1. (3 ponts) A town hall s conductng a study to determne the amount of leftover food produced by the restaurants n the
More informationMarkov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement
Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
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 informationCurve Fitting with the Least Square Method
WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00
ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
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 informationCinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure
nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC
COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Fnte Element Method - Jacques-Hervé SAIAC FINITE ELEMENT METHOD Jacques-Hervé SAIAC Départment de Mathématques, Conservatore Natonal des Arts et Méters, Pars,
More informationResource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis
Resource Allocaton and Decson Analss (ECON 800) Sprng 04 Foundatons of Regresson Analss Readng: Regresson Analss (ECON 800 Coursepak, Page 3) Defntons and Concepts: Regresson Analss statstcal technques
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 informationSTAT 511 FINAL EXAM NAME Spring 2001
STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte
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 informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More informationA boundary element method with analytical integration for deformation of inhomogeneous elastic materials
Journal of Physcs: Conference Seres PAPER OPEN ACCESS A boundary element method wth analytcal ntegraton for deformaton of nhomogeneous elastc materals To cte ths artcle: Moh. Ivan Azs et al 2018 J. Phys.:
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 informationFUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM
Internatonal Conference on Ceramcs, Bkaner, Inda Internatonal Journal of Modern Physcs: Conference Seres Vol. 22 (2013) 757 761 World Scentfc Publshng Company DOI: 10.1142/S2010194513010982 FUZZY GOAL
More informationMultigradient for Neural Networks for Equalizers 1
Multgradent for Neural Netorks for Equalzers 1 Chulhee ee, Jnook Go and Heeyoung Km Department of Electrcal and Electronc Engneerng Yonse Unversty 134 Shnchon-Dong, Seodaemun-Ku, Seoul 1-749, Korea ABSTRACT
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed
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 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 informationThermodynamics General
Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,
More informationTHE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE
More informationMass Transfer Processes
Mass Transfer Processes S. Majd Hassanzadeh Department of Earth Scences Faculty of Geoscences Utrecht Unversty Outlne: 1. Measures of Concentraton 2. Volatlzaton and Dssoluton 3. Adsorpton Processes 4.
More informationLecture 7: Boltzmann distribution & Thermodynamics of mixing
Prof. Tbbtt Lecture 7 etworks & Gels Lecture 7: Boltzmann dstrbuton & Thermodynamcs of mxng 1 Suggested readng Prof. Mark W. Tbbtt ETH Zürch 13 März 018 Molecular Drvng Forces Dll and Bromberg: Chapters
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 informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More informationComputation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models
Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationBayesian predictive Configural Frequency Analysis
Psychologcal Test and Assessment Modelng, Volume 54, 2012 (3), 285-292 Bayesan predctve Confgural Frequency Analyss Eduardo Gutérrez-Peña 1 Abstract Confgural Frequency Analyss s a method for cell-wse
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
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 information