Transactions on Engineering Sciences vol 5, 1994 WIT Press, ISSN
|
|
- Samuel Dorsey
- 5 years ago
- Views:
Transcription
1 Generalize heat conuction transfer functions in two imensional linear systems of assigne geometry F. Marcotullio & A. Ponticiello Dipartimento i Energetica- Universit i L 'Aquila Monteluco i Roio L'Aquila, Italy ABSTRACT The present paper eals with the search for transfer functions vali for the whole class of object of an assigne two imensional geometry. The basic criterion was to turn the governing conuction equation an the transient bounary conitions as well into imensionless variables. The conclusion was that the same the TFs hol for all systems provie the Nusselt number (Nu) an the Fourier moulus (Fo) are equal. An appropriate TF's coefficients table has been set-up for the central an the corner points of a homogeneous soli square profile. The imensionless time interval AFo was selecte accoring to the sampling theorem, i.e. to the cut-off frequency of the system which epens on Nu. In aition, the Nu incremental steps were set in orer to provie "equal" percent eviation (10%) an thus avoi the necessity of temperature interpolations. The generalize TF's coefficients enables a wier use of the TF approach which is very useful for thermal transient analysis. INTRODUCTION The Z-transfer function approach represents a very powerful tool for preicting temperature histories in multiimensional thermal fiels, e.g. Mascheroni [1]. It is particularly useful, for instance, in simulating the ynamic behaviour of linear systems uner ifferent transient bounary conitions, e.g. Stephenson an Mitalas [2]. The search for the transfer functions (TF) in multiimensional heat conuction by a finite-element technique was envisage by the authors [3]. The basic criterion was that of rearranging an conensing the matrix finite-element equation to reuce the egrees of freeom to only one. The conensing proceure will yiel, in general, aitional terms (higher orer time-enva-
2 36 Heat Transfer tives) in the resulting "lumpe" matrix equation. At this stage, the numerical solution in the time-omain iscretization was obtaine explicitly using an appropriate recurrence formula. The coefficients of the Z-transfer functions were thereby compute from the cite recurrence scheme by matching the system free response in orer to set the correct recurrence parameters. In this stuy, the matrix conuction equation is turne into imensionless variables reaching the conclusion that the same TFs hol for the whole class of objects having the same geometric shape (plates, ros, bars, cyliners...) provie the Nusselt number (Nu) an the Fourier moulus (Fo) are equal. This expecte but relevant result enables a wier use of the TF approach, otherwise relegate to rare applications since the TF's ientification presupposes a certain amount of skilfulness which is infrequent in common practice. So far, it will suffice to set-up, at preselecte points an for iffering geometric shapes of practical interest, appropriate tables listing the series of the TF's coefficients for pairs of the inepenent imensionless variables Nu an Fo BASIC CONSIDERATIONS The governing conuction equation in a uniform thermal conuctivity soli ii (x,y,z) free of sources is given by: x* Let our task be the search for the general temperature history within the soli unergoing the following initial an transient bounary conitions: To turn Equation (1) into imensionless variables let: ' - - where L<. is a characteristic length an Fo is the Fourier moulus. By substituting in Equation (1) we obtain:
3 Heat Transfer 37 with the bounary conition expresse as: T (A\ ul ~. XT,, /nr T \ = Q (5) vv < where Nu=hLJk is the Nusselt number. It follows that the analytical solution T will be a function of the imensionless variables %,T,^, Fo, Nu an will not epen explicitly upon the soli thermophysical properties which, instea, are inclue in Fo (iffusivity a) an Nu (conuctivity *). The same conclusions can be rawn by performing a numerical analysis base, for example, on the finite-element metho. In this case, the iscretization of the whole omain Q (%,T,Q in a certain number M of finite elements leas to the well-known matrix equation Ca+Ka-f = 0 (6) here a is the vector of the unknown noal temperature of the system an a the first time-erivative with respect to the Fourier moulus Fo. All the matrices are assemble from the element submatnces in the stanar manner with submatnces K, CT an f given by: K = " j\$ $ T) ri t c^q ; (7a) (7b) (7c) where a is the ambient flui temperature varying with time (forcing function) an Ni, N^ the interpolation an weighting functions. Now, since the
4 38 Heat Transfer transfer functions (TF) at any selecte noe (,T,Q of the fiel are obtaine from the matrix Equation (6) by means of the proceure escribe in [3], it follows that exactly the same TF's hol for any assigne geometry provie the interpolation functions N(,r,Q, the Nusselt number Nu an the imensionless time interval AFo=Aia/L/ are equal. The transfer functions are given as sets of coefficients each relate to its own forcing function represente in general by m inepenent surrouning temperatures. The time-series of the preselecte noal temperature # +, are linke to each timeseries of the input # +, by means of a convolution equation such as: m ^ Nu AND AFo SELECTION To set-up a table listing the series of TF's coefficients for pairs of the inepenent imensionless variables Nu an AFo to be use without interpolation, some appropriate criteria shoul be followe. Firstly, the Fourier moulus iscretization AFo is selecte at least equal to X/15 accoring to the sampling theorem, X being the ominant eigenvalue of the matrix K C obtaine from Equation (6). The above conition implies that the system frequency response to the unitary step-function of the ambient flui temperature is ampe out at approximately 2% of its maximum value. In fact, the temperature history at any arbitrary point is given by: which, as time progresses, can be aequately represente by: T o) The Laplace transfer function is: an by setting the Laplace parameter s=27wi the moulus IW(s)l takes the form:
5 Heat Transfer 39 Table 1 - TF's coefficients for a soli homogeneous square profile 2^ i Nu AFo xlo* L;= t L ^~ I 4-/Vi-y+ tf-/v,-/ = <> 2=0 /= * b* b~ b> b^ b^ b- n n n n n n
6 40 Heat Transfer 1 (12) If a cut-off frequency v, is chosen so as to obtain for Iwl a 2% attenuation, then v,«50/27txi an, corresponingly, AFo=l/2v«.=A,i/15. As it clearly appears from the set of Equations (7), since A,, epens upon Nu then the above time interval selection criterion allows us to set-up the TF's coefficients table with only one entry (Nu). Another problem that arose was the search for the appropriate Nusselt incremental steps ensuring "equal" percent temperature variations. For general engineering purposes, a 10% sprea will suffice to hol the transient response within 5% maximum eviation an thus avoi teious temperature interpolations. As an application example, the TF's coefficients for the central an the corner points of an infinite homogenous plate 0.6m by 0.6m of a soli square profile subject to transient convective flow are reporte in Table 1. Assuming the characteristic length L<=L=0.6m, the unit surface conuctance h=5.8wm^k"'an the thermal conuctivity k=0.65wm"*k% then Nu=hLJ( ^=5.35 against 4.8, which is the nearest Nu value in Table 1. Temperature histories for both noes 1 an 2 were plotte by using the series of coefficients corresponing to Nu=4.8 an AFo=6.510^ an assuming (Figure 1) an ambient flui linear rise temperature (3 C/h) an (Figure 2) a triangular pulse 10 C amplitue an 2Ai uration (AT=AFo-L//a=1.42h). The plots were checke in terms of actual time intervals against the results eriving from the use of the exact sets of TF's coefficients (Table 2). It is confirme that all eviations are within the chosen 5% maximum error range. Table 2 - Exact FT's Coefficients ( Ai=1.42h) Nu n+1 n n n-2 n-3 n b* CONCLUSIONS The following comments can be emphasize from the above stuy: the availability of the transfer function coefficients for the whole class of objects of assigne geometry in terms of imensionless variables such as the Nusselt number Nu an the time interval AFo, enables a wier use of the
7 Heat Transfer =0.65 Wm'iK'1 = 1700kgm" =0.60m = " =5.8 Wm'^K Noe! - Nu=4.8 Noel-Nu=5.35 Noe2 - Nu=4.8 Noe2-Nu= ,81 15,66 18,5 21,36 5,69 ' 8,54 11,39 14,23 17,08 19,93 Time (h) Figure 1 - Temperature histories at noes 1 an 2 with an ambient flui linear rise temperature (3 C/h) '"t" k=0.65 Wnr'K = 1700kgm'* c=836jkg-'k-' L =0.60 m a = ,85 5,69 8,54 11,39 Time (h) 14,23 19,93 21,35 Figure 2 - Temperature histories at noes 1 an 2 with an ambientfluitriangular pulse 2Ar=2.85 h uration an 10 C amplitue
8 42 Heat Transfer transfer function approach. appropriate tables listing the TF's coefficients have been set-up by choosing for the Nusselt number incremental steps of "equal" percent eviation (10%). This criterion permits the irect use of the sets of coefficients by referring to the closest Nu value, without the necessity of teious temperature interpolations. Moreover, by choosing the time interval AFo=l/2Vc=X/15, which epens on Nu, the TF's coefficients table turns into a more suitable one-entry table. REFERENCES 1. Mascheroni, R.H., Sanz, P.D. an Dominguez, M. 'A New Way to Preict Thermal Histories in Multiimensional Heat Conuction: the Z- Transfer Function Metho' Int. Comm. Heat & Mass Transfer, Vol.14, pp , Stephenson, D.G. an Mitalas, G.P. 'Calculation of Heat Conuction Transfer Function for Multilayer Slabs' ASHRAE Transactions, Marcotullio, F an Ponticiello, A. "Determination of Transfer Functions in Multiimensional Heat Conuction by Means of a Finite Element Technique' (Eite by R.W. Lewis, University of Swansea U.K.), Vol.VIII, part 1, sect. 2, pp , Proceeings of the Eighth International Conference on Numerical Methos in Thermal Problems, Swansea, Wales, Salvaori, U.O., Reynoso, R.O. an Mascheroni, R.H. 'The Use of Z- Transfer Functions to Evaluate Temperature an Quality Changes in Refrigerate Storage' Proceeings of meetings of Comm. B2, C2, Dl, D2/3, Int. Institute of Refrigeration, Paris, pp , ASHRAE Hanbook of Founamentals, 1977, Chap Zienkiewicz, O C The Finite Element Metho McGraw-Hill, Lonon, Huebner, K.H. The Finite Element Metho for Engineers J Wiley&Sons, New York, Rao, S.S. The Finite Element Metho in Engineering Pergamon Press, Oxfor, Schneier, P.J. Conuction Heat Transfer Aison-Wesley Publishing Company, Inc, 1974
APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationSemianalytical method of lines for solving elliptic partial dierential equations
Chemical Engineering Science 59 (2004) 781 788 wwwelseviercom/locate/ces Semianalytical metho of lines for solving elliptic partial ierential equations Venkat R Subramanian a;, Ralph E White b a Department
More informationPredictive control of synchronous generator: a multiciterial optimization approach
Preictive control of synchronous generator: a multiciterial optimization approach Marián Mrosko, Eva Miklovičová, Ján Murgaš Abstract The paper eals with the preictive control esign for nonlinear systems.
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationOptimal LQR Control of Structures using Linear Modal Model
Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationMULTI-SCALE METHODS FOR THE SOLUTION OF THE RADIATIVE TRANSFER EQUATION
MULTI-SCALE METHODS FOR THE SOLUTION OF THE RADIATIVE TRANSFER EQUATION Pero J. Coelho *,, Nicolas Crouseilles **, Pero Pereira * an Maxime Roger *** * LAETA, IDMEC, Dept. of Mechanical Engineering, Instituto
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationDAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1
11 th National Congress on Theoretical an Applie Mechanics, 2-5 Sept. 2009, Borovets, Bulgaria DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 E. MANOACH Institute of Mechanics,
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationRecommendations: Part 7: Transient Creep for service and accident conditions
Materials an Structures/Matériaux et Constructions, Vol. 31, June 1998, pp 290-295 RILEM TECHNICAL COMMITTEES RILEM TC 129-MHT: TEST METHODS FOR MECHANICAL PROPERTIES OF CONCRETE AT HIGH TEMPERATURES Recommenations:
More informationModeling time-varying storage components in PSpice
Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC,
More informationSIMULATION OF DIRECT CONTACT CONDENSATION OF STEAM JETS SUBMERGED IN SUBCOOLED WATER BY MEANS OF A ONE-DIMENSIONAL TWO-FLUID MODEL
HEFAT014 10 th International Conference on Heat Transfer, Flui Mechanics an Thermoynamics 14 16 July 014 Orlano, Floria SIMULATION OF DIRECT CONTACT CONDENSATION OF STEAM JETS SUBMERGED IN SUBCOOLED WATER
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationTIME SAMPLING OF DIFFUSION SYSTEMS USING SEMIGROUP DECOMPOSITION METHODS
TIME SAMPLING OF DIFFUSION SYSTEMS USING SEMIGROUP DECOMPOSITION METHODS J. M. Lemos J. M. F. Moura INESC-ID/IST, R. Alves Reol 9, -29 Lisboa, Portugal, jlml@inesc.pt Carnegie Mellon University, Pittsburgh
More informationA SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationNON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS
ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationA Short Note on Self-Similar Solution to Unconfined Flow in an Aquifer with Accretion
Open Journal o Flui Dynamics, 5, 5, 5-57 Publishe Online March 5 in SciRes. http://www.scirp.org/journal/oj http://x.oi.org/.46/oj.5.57 A Short Note on Sel-Similar Solution to Unconine Flow in an Aquier
More informationEntropy Generation Study of MHD Thermosolutal Convection in a Square Cavity for Different Prandtl Numbers
International Journal of Mechanics an Applications. ; (): -9 DOI:. 593/j.mechanics..3 Entropy Generation Stuy of MHD Thermosolutal Convection in a Square Cavity for Different Prantl Numbers Mounir Bouabi
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationarxiv:quant-ph/ v1 29 Jun 2001
Atomic wave packet basis for quantum information Ashok Muthukrishnan an C. R. Strou, Jr. The Institute of Optics, University of Rochester, Rochester, New York 14627 (March 15, 2001) arxiv:quant-ph/0106165
More informationImproving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers
International Journal of Statistics an Probability; Vol 6, No 5; September 207 ISSN 927-7032 E-ISSN 927-7040 Publishe by Canaian Center of Science an Eucation Improving Estimation Accuracy in Nonranomize
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationAn analytical investigation into filmwise condensation on a horizontal tube in a porous medium with suction at the tube surface
Heat Mass Transfer (29) 45:355 361 DOI 1.17/s231-8-436-y ORIGINAL An analytical investigation into filmwise conensation on a horizontal tube in a porous meium with suction at the tube surface Tong Bou
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More information[Kumar*, 5(2): February, 2016] ISSN: (I2OR), Publication Impact Factor: 3.785
[Kumar*, 5(): February, 6] ISSN: 77-9655 (IOR), Publication Impact Factor: 3.785 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY THERMOSOLUTAL CONVECTION IN A HETEROGENEOUS VISCO-ELASTIC
More informationSome New Thoughts on the Multipoint Method for Reactor Physics Applications. Sandra Dulla, Piero Ravetto, Paolo Saracco,
Jeju, Korea, April 16-20, 2017, on USB 2017 Some New Thoughts on the Multipoint Metho for Reactor Physics Applications Sanra Dulla, Piero Ravetto, Paolo Saracco, Politecnico i Torino, Dipartimento Energia,
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationA variational approach to Bogoliubov excitations and dynamics of dipolar Bose Einstein condensates
IOP PUBLISHING JOURNAL OF PHYSICS B: ATOMIC, MOLECULAR AND OPTICAL PHYSICS J. Phys. B: At. Mol. Opt. Phys. () 5 (pp) oi:./5-5///5 A variational approach to Bogoliubov excitations an ynamics of ipolar Bose
More informationDesigning Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations
EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through
More informationInfluence of Radiation on Product Yields in a Film Boiling Reactor
R&D NOTES Influence of Raiation on Prouct Yiels in a Film Boiling Reactor C. Thomas Aveisian, Wing Tsang, Terence Daviovits, an Jonah R. Allaben Sibley School of Mechanical an Aerospace Engineering, Cornell
More informationApproaches for Predicting Collection Efficiency of Fibrous Filters
Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationShape Effect on Blind Frequency for Depth Inversion in Pulsed Thermography
Shape Effect on Blin Frequency for Depth Inversion in Pulse Thermography M. Genest 1, E. Grinzato 2, P. Bison 2, S. Marinetti 2 C. Ibarra-Castaneo 1, X. Malague 1 1 Electrical an Computing Eng. Dept.,
More informationStochastic Eutectic Growth
VOLUME 72, NUMBER 5 PHYSCAL REVEW LETTERS Stochastic Eutectic Growth 31 JANUARY 1994 K. R. Eler, Frangois Drolet, J. M. Kosterlitz 2 an Martin Granti Department of Physics, McGill University, Rutherfor
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationarxiv: v1 [hep-lat] 19 Nov 2013
HU-EP-13/69 SFB/CPP-13-98 DESY 13-225 Applicability of Quasi-Monte Carlo for lattice systems arxiv:1311.4726v1 [hep-lat] 19 ov 2013, a,b Tobias Hartung, c Karl Jansen, b Hernan Leovey, Anreas Griewank
More informationEssential Considerations for Buckling Analysis
Worlwie Aerospace Conference an Technology Showcase, Toulouse, France, Sept. 24-26, 2001 Essential Consierations for Buckling Analysis 2001-120 Sang H. Lee MSC.Software Corporation, 2 MacArthur Place,
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationInvestigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model
ORIGINAL ARTICLE Investigation of local loa effect on amping characteristics of synchronous generator using transfer-function block-iagram moel Pichai Aree Abstract of synchronous generator using transfer-function
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationEstimation of the Maximum Domination Value in Multi-Dimensional Data Sets
Proceeings of the 4th East-European Conference on Avances in Databases an Information Systems ADBIS) 200 Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Eleftherios Tiakas, Apostolos.
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationAverage value of position for the anharmonic oscillator: Classical versus quantum results
verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationFINITE POINTSET METHOD (FPM): A MESHFREE APPROACH FOR INCOMPRESSIBLE FLOW SIMULATIONS APPLIED TO COMPOSITE MATERIALS
FINITE POINTSET METHOD (FPM): A MESHFREE APPROACH FOR INCOMPRESSIBLE FLOW SIMULATIONS APPLIED TO COMPOSITE MATERIALS Alain Trameçon 1, Patrick e Luca 1, Christophe Binetruy 2 an Jörg Kuhnert 3 1 ESI-GROUP
More informationThe Second Order Contribution to Wave Crest Amplitude Random Simulations and NewWave
The Secon Orer Contribution to Wave Crest Amplitue Ranom Simulations an NewWave Thomas A.A. Acock Department of Engineering Science University of Oxfor Parks Roa Oxfor Unite Kingom Scott Draper School
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E.
More informationStudy on aero-acoustic structural interactions in fan-ducted system
Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,
More informationEvaporating droplets tracking by holographic high speed video in turbulent flow
Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)
Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution
More informationNumerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number
Journal of Applie Flui Mechanics, Vol. 9, No. 2, pp. 643-651, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Investigation of Non-Stationary Parameters on Effective
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationADIT DEBRIS PROJECTION DUE TO AN EXPLOSION IN AN UNDERGROUND AMMUNITION STORAGE MAGAZINE
ADIT DEBRIS PROJECTION DUE TO AN EXPLOSION IN AN UNDERGROUND AMMUNITION STORAGE MAGAZINE Froe Opsvik, Knut Bråtveit Holm an Svein Rollvik Forsvarets forskningsinstitutt, FFI Norwegian Defence Research
More informationSharp Fickian Behavior of Electrogenerated Chemiluminescence Waves in Conjugated Polymer Films
pubs.acs.org/jpcc Sharp Fickian Behavior of Electrogenerate Chemiluminescence Waves in Conjugate Polymer Films Davi A. Ewars* Department of Mathematical Sciences, University of Delaware, Newark, Delaware
More informationTopic 7: Convergence of Random Variables
Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Moels for CFD ME 448/548 Notes Geral Recktenwal Portlan State University Department of Mechanical Engineering gerry@p.eu ME 448/548: D Convection-Di usion Equation
More informationAdjoint Transient Sensitivity Analysis in Circuit Simulation
Ajoint Transient Sensitivity Analysis in Circuit Simulation Z. Ilievski 1, H. Xu 1, A. Verhoeven 1, E.J.W. ter Maten 1,2, W.H.A. Schilers 1,2 an R.M.M. Mattheij 1 1 Technische Universiteit Einhoven; e-mail:
More information(3-3) = (Gauss s law) (3-6)
tatic Electric Fiels Electrostatics is the stuy of the effects of electric charges at rest, an the static electric fiels, which are cause by stationary electric charges. In the euctive approach, few funamental
More informationSlide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)
Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More information