Calculation of Power Density with MCNP in TRIGA Reactor
|
|
- Beverly O’Connor’
- 6 years ago
- Views:
Transcription
1 Internatonal Conerence Nuclear Energy or New Europe 26 Portorož, Slovena, September 18-21, 26 Calculaton o Power Densty wth MCNP n TRIGA Reactor ABSTRACT Luka Snoj, Matjaž Ravnk Jože Stean Insttute Jamova 39, SI-1 Ljubljana, Slovena Luka.Snoj@js.s, Matjaz.Ravnk@js.s Modern Monte Carlo codes (e.g. MCNP allow calculaton o power densty dstrbuton n 3-D geometry assumng detaled geometry wthout unt-cell homogenzaton. To normalze MCNP calculaton by the steady-state thermal power o a reactor, one must use approprate scalng actors. The descrpton o the scalng actors s not adequately descrbed n the MCNP manual and requres detaled knowledge o the code model. As the applcaton o MCNP or power densty calculaton n TRIGA reactors has not been reported n open lterature, the procedure o calculatng power densty wth MCNP and ts normalzaton to the power level o a reactor s descrbed n the paper. 1 INTRODUCTION Monte Carlo computer code, MCNP, s a very powerul and versatle tool or partcle transport calculatons. It can be used or transport o neutrons, photons and electrons. Transport o neutrons s o specal nterest or a reactor physcst. MCNP code can be used or calculatons o multplcaton actor, reacton rates, saturated actvtes, neutron luxes and spectra, power peakng actors, reacton rate dstrbutons, sheldng etc. Its man advantage s the ablty to handle complcated geometres. MCNP also provdes seven standard tally types [1]. All talles are normalzed to one "startng" partcle except n KCODE crtcalty problems, whch are normalzed to one sson neutron. In order to normalze the result by the thermal power o a system, one must use approprate scalng actors. As the procedure o calculatng the scalng actors or KCODE calculaton s not adequately descrbed n the MCNP manual, t wll be descrbed n the paper. The present paper wll ocus on standard cell lux (track length estmate o cell lux F4 tally [1] whch s o greatest nterest or reactor physcs calculatons. However the results can be appled also to other neutron talles. 2 NORMALISATION OF F4 TALLIES IN A CRITICALITY CALCULATION The multplcaton actor s one o the most mportant propertes o a reactor or other system made o ssle materal. In MCNP the most common way to calculate the multplcaton actor s through the use o KCODE card. It s mportant to note that all the standard MCNP talles can be made durng a crtcalty calculaton. Snce the MCNP results are normalsed to one source neutron, the result has to be properly scaled n order to get absolute comparson to the measured quanttes (lux, reacton rate, sson densty, etc.. The F4 tally results can be scaled to a desred sson neutron 19.1
2 19.2 source (power level or total neutron pulse strength. The scalng actor can be entered on the FM (tally multpler card or can be appled later n data processng. The neutron brth rate n a ssle system can be calculated rom the released energy per P unt tme.e. the power o the system. The system producng power P needs ssons per w second, where w denotes eectve energy released per sson event. Although the value o w wll vary somewhat wth the type o reactor and the detaled core composton, t s typcally Pν o the order o 198 MeV or steady state condton. Ths sson rate produces neutrons w per second, where ν denotes the average number o neutrons released per sson. (The value o ν s lsted n the MCNP output le n the box contanng the nal k e result and represents the value averaged over ssle sotopes and neutron energes. Thereore to normalze an F4 tally by the steady-state thermal power o a crtcal system, the ollowng scalng actor n unts o sson neutrons per unt tme should be used neutron P[ W] ν sson S = 13 J MeV w MeV sson (1. The upper scalng actor s approprate or crtcal.e. steady-state power level systems only (k = 1. KCODE talles or subcrtcal and supercrtcal systems do not nclude any multplcaton eects because sson s treated as absorpton. Thereore one must multply the equaton (1 by 1 on the rght-hand sde or subcrtcal and supercrtcal systems, ke respectvely. It s mportant to note that the scalng actor or subcrtcal and supercrtcal systems s vald only when one uses neutron source dstrbuton dentcal n space and energy to the source dstrbuton obtaned rom the soluton o an egenvalue problem wth k e 1. To conclude, when one wants to scale the calculated F4 tally "lux", Φ F4, one must use the ollowng equaton neutron P[ W] ν neutron sson 1 1 Φ = Φ cm s J MeV MeV sson 2 F ke cm w where Φ denotes the actual total neutron lux n the system., (2 3 CALCULATION OF POWER DENSITY DISTRIBUTION 3.1 Calculaton o power densty The energy released n a nuclear sson reacton s dstrbuted among a varety o reacton products. The majorty o the sson energy appears as the knetc energy o the sson ragments and s deposted essentally at the pont o sson. About 97 % o the Proceedngs o the Internatonal Conerence Nuclear Energy or New Europe, 26
3 19.3 recoverable sson energy s deposted drectly n the ssle materal [2]. In our calculaton o the power densty dstrbuton we wll assume that power densty s proportonal to sson densty. In other word that means that we assume that all o the recoverable sson energy s deposted at the pont o sson. We wll also assume that there are no temperature eedback eects and that there s only one ssle sotope n the system. Power densty, dened as the energy deposted n the ssle materal per unt volume per unt tme, can be wrtten as p( r = w ( E, r Σ ( E, r ϕ( E, r de, (3 where Σ and φ denote macroscopc sson cross secton and neutron spectrum, respectvely. Neutron spectrum φ(e,r s normalzed such as ϕ( E, r de =Φ( r, (4 where Φ(r s total neutron lux n cm -2 s -1. The same s vald also or MCNP F4 tally ϕ F ( E, r de =Φ ( r. (5 4 F4 Equaton (3 represents thermal power densty at poston r n the sson system. Hence the total power generated by the sson system s just the ntegral o the power densty over the total volume where Σ. 3 P = d r w ( E, r Σ ( E, r ϕ( E, r de V, Σ (6 Assumng that w and number densty o the ssle materal do not depend sgncantly on the energy and the poston n the ssle system, equaton (3 can be wrtten as p( r = wn σ (, (, r E, (7 E r ϕ E d where N denotes ssle materal atom densty and σ denotes mcroscopc total sson cross secton. In MCNP the F4 cell lux tally, Φ F4 (r, s averaged over the volume o the samplng cell. Thereore, when we calculate sson densty n a certan cell, Δ V ( runs over all spatal cells o the reactor, we can omt the spatal dependence rom the ntegral n equaton (7 and obtan p( r = w N σ ( (, E ϕ E de r V (8 Δ The ntegral can be qute easly calculated by usng the tally multpler (FM card [1], whch s used to calculate, or Δ V, any quantty o the orm Proceedngs o the Internatonal Conerence Nuclear Energy or New Europe, 26
4 Emax 19.4 FR = C σr( E ϕf4( E de, (9 Emn where the constant C s any arbtrary scalar quantty that can be used or normalsaton, σ R (E s mcroscopc cross secton or reacton R taken rom MCNP cross-secton lbrares and E mn and E max denote mnmum and maxmum neutron energy n the system (usually and 2 MeV, respectvely. F R s the quantty or reacton R that s calculated by MCNP and ts value can be ound n the output le. Any number o ENDF or specal reactons can be used as a multpler as long as they are present n the MCNP cross-secton lbrares. I the C entry s negatve (or type 4 tally only, C s replaced by C tmes the atom densty o the cell where the tally s made. Thus or calculatng power densty we use F4 tally wth tally multpler -6, whch s the mcroscopc total sson cross secton, and normalse the result to a desred sson neutron source (power level or total neutron pulse strength. Usng equatons (8, (9 and the scalng actor S we obtan p( r = Pν M w w Emax σ ( ϕf4(, r E Δ V (1 mn N E E de It s nterestng that the power densty s ndependent on w. Fnally we can wrte power densty as p( r = Pν NMF, r V. (11 Note that p(r s a uncton o r, as t s constant nsde each Δ V. I we want to obtan "smooth" p(r dstrbuton, all Δ V must be small. It s mportant to note that σ n equaton (1 s actually macroscopc total sson cross secton normalzed to one atom o the sampled materal. Thereore n case the sampled materal contans ssle and non-ssle components (e.g. UO 2, UCO, UZrH, UO 2 + PuO 2, etc. the ntegral n equaton (1 can be smply multpled by the sampled materal atom densty,.e. a mxture o several nucldes, and not the atom densty o the ssle materal only. I the sampled materal contans several ssle sotopes, mcroscopc total sson cross secton, σ, s calculated as the weghted average over all ssle sotopes, wth regard to ther atom racton. The same s true orν. Δ 3.2 Alternatve opton or calculaton o power densty The alternatve opton or calculaton o power densty s much smpler, but t has some lmtatons. Ths opton s partcularly useul when we want to calculate the power densty n ndvdual cells or uel elements. Frst we calculate sson densty (normalsed "per atom".e. calculate F n each cell contanng ssle materal and obtan E max F = Φ( E σ ( E de E mn (, (12 Proceedngs o the Internatonal Conerence Nuclear Energy or New Europe, 26
5 19.5 where denotes the cell ndex. Snce the power produced n one cell s lnearly proportonal to the number o ssons n that cell, we can calculate the power produced n cell, P, by multplyng the ssle system thermal power, P, by the relatve number o ssons n cell. As the relatve number o ssons n cell s proportonal to the product o sson densty n cell and the volume o cell, we obtan P ( F ( F V = P. (13 V Usng the denton o power densty as beng the power produced per unt volume we obtan the power densty n cell as p ( F ( F = P. (14 V The man lmtatons o ths procedure o sson densty calculaton are that we have to know the exact volume o each cell and that we have to sample all cells contanng ssle materal, whch s not necessary or the prevously descrbed opton. Note that V can be the volume o the uel nsde the uel element or very small volume o arbtrary materal composton and needs not to contan only ssle materal. In both cases p corresponds to the average power densty n that volume. 3.3 Calculaton o power densty dstrbuton The easest way or calculatng the power densty dstrbuton can be by usng the supermposed mesh tally card, FMESHn (or the present only type 4 talles are permtted. FMESH card allows the user to dene a mesh tally supermposed over the problem geometry. By deault, the mesh tally calculates the track length estmate o partcle lux, averaged over a mesh cell. I we use the tally multpler card or sson together wth the mesh tally and scale the results to approprate power level, we can obtan power densty dstrbuton. When usng tally multpler card together wth supermposed mesh tally card t s recommended to set materal number to, whch causes that the reacton cross sectons or the materal n whch the partcle s travellng are used. Thus we do not have to worry n whch materal partcle s travellng. The FMESH card s extremely powerul and useul method also or calculatng lux dstrbutons, power peakng actors, etc.. By usng the very ne mesh we can calculate local power peakngs and power densty dstrbutons near the water channels, whch s very dcult not mpossble to do wth determnstc methods. 4 POWER DENSITY DISTRIBUTION IN TRIGA MARK II CORE The supermposed mesh tally eature o the MCNP code was used to calculate detaled power densty dstrbuton and power peakng actors o varous hypothetcal mxed cores o TRIGA Mark II reactor located at Joze Stean Insttute. The mesh o the mesh tally was so ne that power densty dstrbuton was calculated wth a resoluton o 1mm 1mm. The dsadvantage o such hgh resoluton s the very large number o samplng cells (5 5 Proceedngs o the Internatonal Conerence Nuclear Energy or New Europe, 26
6 19.6 and consequently relatvely hgh relatve error o the tally. That s the man reason why we calculated the power densty dstrbuton averaged over the uel heght. The results are presented n Fgure 1. Fgure 1: Power densty dstrbuton (rel. unts n mxed core wth 8.5 w/o standard uel and 6 LEU uel elements n the D rng. The numbers on x n y axs are dstances rom the centre o the core n cm. 5 CONCLUSIONS We have shown how to scale MCNP tally results to a desred power level o the system. The scalng actors can be derved rom the system power, multplcaton actor, average number o neutrons released per sson and average energy produced per sson event. When one s amlar wth normalzng the results to a certan power level, the MCNP becomes a versatle calculaton tool or varous reactor physcs calculatons. One can calculate reacton rates, saturated actvtes, neutron luxes and spectra, power peakng actors etc.. The MCNP ablty to handle complcated geometres and use o the FMESH card eature o the MCNP code enables one to calculate detaled neutron lux and varous reacton rate dstrbutons wth a resoluton o 1mm 1mm or less (depends on the sze o the system, the computer hardware and sotware. The later eature s especally useul or calculaton o power densty dstrbuton and local power peakng actors calculaton snce t can reproduce properly the eects o local power gradents due to small heterogenetes n the core. REFERENCES [1] J.J. Bresmester, "MCNP5- A General Monte Carlo N-Partcle Transport Code, Verson 5 Los Alamos Natonal Laboratory", March, 25 [2] J. J. Duderstadt, L. J. Hamlton, Nuclear Reactor Analyss, John Wley & Sons, 1976 Proceedngs o the Internatonal Conerence Nuclear Energy or New Europe, 26
Physics 2A Chapters 6 - Work & Energy Fall 2017
Physcs A Chapters 6 - Work & Energy Fall 017 These notes are eght pages. A quck summary: The work-energy theorem s a combnaton o Chap and Chap 4 equatons. Work s dened as the product o the orce actng on
More informationMeasurement of Radiation: Exposure. Purpose. Quantitative description of radiation
Measurement of Radaton: Exposure George Starkschall, Ph.D. Department of Radaton Physcs U.T. M.D. Anderson Cancer Center Purpose To ntroduce the concept of radaton exposure and to descrbe and evaluate
More informationUniversity of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014
Lecture 16 8/4/14 Unversty o Washngton Department o Chemstry Chemstry 452/456 Summer Quarter 214. Real Vapors and Fugacty Henry s Law accounts or the propertes o extremely dlute soluton. s shown n Fgure
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 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 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 informationPhysics for Scientists and Engineers. Chapter 9 Impulse and Momentum
Physcs or Scentsts and Engneers Chapter 9 Impulse and Momentum Sprng, 008 Ho Jung Pak Lnear Momentum Lnear momentum o an object o mass m movng wth a velocty v s dened to be p mv Momentum and lnear momentum
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 informationSPANC -- SPlitpole ANalysis Code User Manual
Functonal Descrpton of Code SPANC -- SPltpole ANalyss Code User Manual Author: Dale Vsser Date: 14 January 00 Spanc s a code created by Dale Vsser for easer calbratons of poston spectra from magnetc spectrometer
More information5.15 MICRO GAMMA SCANNING ON THE HIGH BURNUP PWR FUEL SAMPLES
5.15 MICRO GAMMA SCANNING ON THE HIGH BURNUP PWR FUEL SAMPLES Hyoung-Mun Kwon, Hang-Seog Seo, Hyung-Kwon Lee, Duck-Kee Mn, Yong-Bum Chun Post Irradaton Examnaton Faclty, Korea Atomc Energy Research Insttute
More informationWeek 6, Chapter 7 Sect 1-5
Week 6, Chapter 7 Sect 1-5 Work and Knetc Energy Lecture Quz The frctonal force of the floor on a large sutcase s least when the sutcase s A.pushed by a force parallel to the floor. B.dragged by a force
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 informationExperiment 1 Mass, volume and density
Experment 1 Mass, volume and densty Purpose 1. Famlarze wth basc measurement tools such as verner calper, mcrometer, and laboratory balance. 2. Learn how to use the concepts of sgnfcant fgures, expermental
More informationNEUTRON FLUX IN THE EXPOSURE ROOM OF THE TRIGA MARK II REACTOR IN LJUBLJANA
Internatonal Conference Nuclear Energy n Central Europe 2000 Golf Hotel, Bled, Slovena, September 11-14, 2000 NEUTRON FLUX IN THE EXPOSURE ROOM OF THE TRIGA MARK II REACTOR IN LJUBLJANA Edvard S. Krštof
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationEntropy generation in a chemical reaction
Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationFinite Difference Method
7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use
More informationGeneral Tips on How to Do Well in Physics Exams. 1. Establish a good habit in keeping track of your steps. For example, when you use the equation
General Tps on How to Do Well n Physcs Exams 1. Establsh a good habt n keepng track o your steps. For example when you use the equaton 1 1 1 + = d d to solve or d o you should rst rewrte t as 1 1 1 = d
More informationMarkov Chain Monte Carlo Lecture 6
where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways
More informationDesigning of Combined Continuous Lot By Lot Acceptance Sampling Plan
Internatonal Journal o Scentc Research Engneerng & Technology (IJSRET), ISSN 78 02 709 Desgnng o Combned Contnuous Lot By Lot Acceptance Samplng Plan S. Subhalakshm 1 Dr. S. Muthulakshm 2 1 Research Scholar,
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 informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
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 information2016 Wiley. Study Session 2: Ethical and Professional Standards Application
6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton
More informationThe exponential map of GL(N)
The exponental map of GLN arxv:hep-th/9604049v 9 Apr 996 Alexander Laufer Department of physcs Unversty of Konstanz P.O. 5560 M 678 78434 KONSTANZ Aprl 9, 996 Abstract A fnte expanson of the exponental
More informationONE-DIMENSIONAL COLLISIONS
Purpose Theory ONE-DIMENSIONAL COLLISIONS a. To very the law o conservaton o lnear momentum n one-dmensonal collsons. b. To study conservaton o energy and lnear momentum n both elastc and nelastc onedmensonal
More informationComputational Biology Lecture 8: Substitution matrices Saad Mneimneh
Computatonal Bology Lecture 8: Substtuton matrces Saad Mnemneh As we have ntroduced last tme, smple scorng schemes lke + or a match, - or a msmatch and -2 or a gap are not justable bologcally, especally
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationPhysics 2A Chapter 3 HW Solutions
Phscs A Chapter 3 HW Solutons Chapter 3 Conceptual Queston: 4, 6, 8, Problems: 5,, 8, 7, 3, 44, 46, 69, 70, 73 Q3.4. Reason: (a) C = A+ B onl A and B are n the same drecton. Sze does not matter. (b) C
More informationarxiv:cs.cv/ Jun 2000
Correlaton over Decomposed Sgnals: A Non-Lnear Approach to Fast and Effectve Sequences Comparson Lucano da Fontoura Costa arxv:cs.cv/0006040 28 Jun 2000 Cybernetc Vson Research Group IFSC Unversty of São
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 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 informationA Simple Research of Divisor Graphs
The 29th Workshop on Combnatoral Mathematcs and Computaton Theory A Smple Research o Dvsor Graphs Yu-png Tsao General Educaton Center Chna Unversty o Technology Tape Tawan yp-tsao@cuteedutw Tape Tawan
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 informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More 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 informationChapter 6. Operational Amplifier. inputs can be defined as the average of the sum of the two signals.
6 Operatonal mpler Chapter 6 Operatonal mpler CC Symbol: nput nput Output EE () Non-nvertng termnal, () nvertng termnal nput mpedance : Few mega (ery hgh), Output mpedance : Less than (ery low) Derental
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 informationEndogenous timing in a mixed oligopoly consisting of a single public firm and foreign competitors. Abstract
Endogenous tmng n a mxed olgopoly consstng o a sngle publc rm and oregn compettors Yuanzhu Lu Chna Economcs and Management Academy, Central Unversty o Fnance and Economcs Abstract We nvestgate endogenous
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
TESTING THE SPECTRAL DECONVOLUTION ALGORITHM TOOL (SDAT) WITH XE SPECTRA Steven R. Begalsk, Kendra M. Foltz Begalsk, and Derek A. Haas The Unversty of Texas at Austn Sponsored by Army Space and Mssle Defense
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015
Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want
More informationPHYS 1441 Section 002 Lecture #15
PHYS 1441 Secton 00 Lecture #15 Monday, March 18, 013 Work wth rcton Potental Energy Gravtatonal Potental Energy Elastc Potental Energy Mechancal Energy Conservaton Announcements Mdterm comprehensve exam
More informationTHEOREMS OF QUANTUM MECHANICS
THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn
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 informationModule 11 Design of Joints for Special Loading. Version 2 ME, IIT Kharagpur
Module 11 Desgn o Jonts or Specal Loadng Verson ME, IIT Kharagpur Lesson 1 Desgn o Eccentrcally Loaded Bolted/Rveted Jonts Verson ME, IIT Kharagpur Instructonal Objectves: At the end o ths lesson, the
More informationSpring Force and Power
Lecture 13 Chapter 9 Sprng Force and Power Yeah, energy s better than orces. What s net? Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi IN THIS CHAPTER, you wll learn how to solve problems
More informationPeriod & Frequency. Work and Energy. Methods of Energy Transfer: Energy. Work-KE Theorem 3/4/16. Ranking: Which has the greatest kinetic energy?
Perod & Frequency Perod (T): Tme to complete one ull rotaton Frequency (): Number o rotatons completed per second. = 1/T, T = 1/ v = πr/t Work and Energy Work: W = F!d (pcks out parallel components) F
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 information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
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 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 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 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 informationSupplemental Material: Causal Entropic Forces
Supplemental Materal: Causal Entropc Forces A. D. Wssner-Gross 1, 2, and C. E. Freer 3 1 Insttute for Appled Computatonal Scence, Harvard Unversty, Cambrdge, Massachusetts 02138, USA 2 The Meda Laboratory,
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationSOLITARY BURN-UP WAVE SOLUTION IN A MULTI-GROUP DIFFUSION-BURNUP COUPLED SYSTEM
3th Internatonal Conerence on Emergng uclear Energy Systems (ICEES 27) Istanbul Turkey June 3-8 27 on CD-ROM Gaz Unversty Ankara Turkey (27) SOLITARY BUR-UP WAVE SOLUTIO I A MULTI-GROUP DIFFUSIO-BURUP
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationEMU Physics Department
Physcs 0 Lecture 8 Potental Energy and Conservaton Assst. Pro. Dr. Al ÖVGÜN EMU Physcs Department www.aovgun.com Denton o Work W q The work, W, done by a constant orce on an object s dened as the product
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationNeutral-Current Neutrino-Nucleus Inelastic Reactions for Core Collapse Supernovae
Neutral-Current Neutrno-Nucleus Inelastc Reactons for Core Collapse Supernovae A. Juodagalvs Teornės Fzkos r Astronomjos Insttutas, Lthuana E-mal: andrusj@tpa.lt J. M. Sampao Centro de Físca Nuclear da
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 informationModeling and Simulation NETW 707
Modelng and Smulaton NETW 707 Lecture 5 Tests for Random Numbers Course Instructor: Dr.-Ing. Magge Mashaly magge.ezzat@guc.edu.eg C3.220 1 Propertes of Random Numbers Random Number Generators (RNGs) must
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 information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2019 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons Our goal to derve the form of the abstract quanttes n rate equatons, such as synaptc
More informationConvergence of random processes
DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More information), it produces a response (output function g (x)
Lnear Systems Revew Notes adapted from notes by Mchael Braun Typcally n electrcal engneerng, one s concerned wth functons of tme, such as a voltage waveform System descrpton s therefore defned n the domans
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 informationOn the correction of the h-index for career length
1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat
More informationMAE140 - Linear Circuits - Winter 16 Midterm, February 5
Instructons ME140 - Lnear Crcuts - Wnter 16 Mdterm, February 5 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator
More information: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:
764: Numercal Analyss Topc : Soluton o Nonlnear Equatons Lectures 5-: UIN Malang Read Chapters 5 and 6 o the tetbook 764_Topc Lecture 5 Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationHomework Notes Week 7
Homework Notes Week 7 Math 4 Sprng 4 #4 (a Complete the proof n example 5 that s an nner product (the Frobenus nner product on M n n (F In the example propertes (a and (d have already been verfed so we
More informationJoint Statistical Meetings - Biopharmaceutical Section
Iteratve Ch-Square Test for Equvalence of Multple Treatment Groups Te-Hua Ng*, U.S. Food and Drug Admnstraton 1401 Rockvlle Pke, #200S, HFM-217, Rockvlle, MD 20852-1448 Key Words: Equvalence Testng; Actve
More informationElastic Collisions. Definition: two point masses on which no external forces act collide without losing any energy.
Elastc Collsons Defnton: to pont asses on hch no external forces act collde thout losng any energy v Prerequstes: θ θ collsons n one denson conservaton of oentu and energy occurs frequently n everyday
More informationDepartment of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution
Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable
More informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More informationLecture 3 Stat102, Spring 2007
Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture
More informationExercises of Chapter 2
Exercses of Chapter Chuang-Cheh Ln Department of Computer Scence and Informaton Engneerng, Natonal Chung Cheng Unversty, Mng-Hsung, Chay 61, Tawan. Exercse.6. Suppose that we ndependently roll two standard
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationModerator & Moderator System
NPTL Chemcal ngneerng Nuclear Reactor Technology Moderator & Moderator System K.S. Rajan Professor, School of Chemcal & Botechnology SASTRA Unversty Jont Intatve of IITs and IISc Funded by MHRD Page of
More informationLecture 2 Solution of Nonlinear Equations ( Root Finding Problems )
Lecture Soluton o Nonlnear Equatons Root Fndng Problems Dentons Classcaton o Methods Analytcal Solutons Graphcal Methods Numercal Methods Bracketng Methods Open Methods Convergence Notatons Root Fndng
More informationThe Expectation-Maximization Algorithm
The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.
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 informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationWork is the change in energy of a system (neglecting heat transfer). To examine what could
Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes
More informationPHYS 1441 Section 002 Lecture #16
PHYS 1441 Secton 00 Lecture #16 Monday, Mar. 4, 008 Potental Energy Conservatve and Non-conservatve Forces Conservaton o Mechancal Energy Power Today s homework s homework #8, due 9pm, Monday, Mar. 31!!
More informationOld Dominion University Physics 420 Spring 2010
Projects Structure o Project Reports: 1 Introducton. Brely summarze the nature o the physcal system. Theory. Descrbe equatons selected or the project. Dscuss relevance and lmtatons o the equatons. 3 Method.
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 informationThis chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.
1 Unform Electron Gas Ths chapter llustrates the dea that all propertes of the homogeneous electron gas (HEG) can be calculated from electron densty. Intutve Representaton of Densty Electron densty n s
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 informationSociété de Calcul Mathématique SA
Socété de Calcul Mathématque SA Outls d'ade à la décson Tools for decson help Probablstc Studes: Normalzng the Hstograms Bernard Beauzamy December, 202 I. General constructon of the hstogram Any probablstc
More information36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to
ChE Lecture Notes - D. Keer, 5/9/98 Lecture 6,7,8 - Rootndng n systems o equatons (A) Theory (B) Problems (C) MATLAB Applcatons Tet: Supplementary notes rom Instructor 6. Why s t mportant to be able to
More information