Pressure limit of hydrogen spontaneous ignition in a T-shaped channel
|
|
- Ralf Mitchell
- 5 years ago
- Views:
Transcription
1 4 th International Conference on Hydrogen Safety, September 2011, San Francisco, USA Pressre limit of hydrogen spontaneos ignition in a T-shaped channel Maim Bragin, Dmitriy Makarov, Vladimir Molkov Hydrogen Safety Engineering and Research (HySAFER) Centre, University of Ulster
2 Backgrond Aim and obectives Eperiment in T-shape channel (Golb et al., 2010) LES EDC model Modelling of membrane opening Instant verss finite membrane opening Simlation reslts Otline Dynamics of spontaneos ignition at different storage pressre Conclsions
3 Backgrond Pressre relief device (PRD) activation Sdden release of hydrogen into air What is the lower pressre limit for spontaneos ignition in T-shape channel (PRD mock-p)?
4 Aim and obectives Aim: Develop contemporary model for predictive simlation of the phenomenon of hydrogen spontaneos ignition dring sdden release from pressre relief devices (to be applied as a tool for hydrogen safety engineering) Obectives: Develop Large Eddy Simlation (LES) model based on the Eddy Dissipation Concept (EDC), fll chemistry, and modelling of a brst disk opening in a finite time Validate the LES EDC model and simlations against eperimental data by Golb et al. (2010) on spontaneos ignition in T-shape channel (mock-p PRD) Understand dynamics of the ignition process and eplain scattering in reported eperimental data (Golb et al., 2010) Establish the lower pressre limit for spontaneos ignition in T-shape channel
5 Eperiment (Golb et al., 2010) Brst disk rptre at storage pressres from 13.5 to 29 bar. Reported: no ignition at 13.5 bar, and spontaneos ignition at 29 bar (private commnication ignition at 24.3 bar) Brst disk ID 4 ID 16 ID 10 ID High-pressre section (hydrogen) PRD (air)
6 Comptational domain High pressre tbes PRD Eternal domain Brst disk
7 LES eqations ( ) 0 = t ρ ρ ( ) i i k k i i i i i g p t ρ δ μ ρ ρ = 3 2 ( ) ( ) ( )= p E E t ρ ρ m m m p Y Sc h T c Pr [ μ μ ( ) ( ) m m m H R Y Sc Y Y t = μ ρ ρ Continity eqation: Momentm conservation eqation: Energy conservation eqation: Species conservation eqation: m C m i k k i i i H R ] 3 2 δ μ
8 RNG trblence model The ective viscosity in RNG SGS model is calclated as μ 2 μs μ = μ 1 H μ , where μ = ( 0.157VCV ) 2Si Si In highly trblent flows, μ μ s and the RNG model redces to Smagorinsky model in laminar flow regions the argment of Heaviside fnction becomes negative and the model recovers moleclar viscosity, = μ. The ective Prandtl and Schmidt nmbers are calclated following the RNG theory by Yakhot prely from the theoretical eqation 1 1 N N μ N N S = s i = ρ μ μ 1 2 i i
9 EDC combstion model 1/2 Eddy Dissipation Concept (EDC) gives an epression for a combstion rate based on an assmption that chemical reactions occr in so-called fine strctres of Kolmogorov s scale, where the dissipation of trblence energy takes place. In EDC model a sorce term in species ( ) 2 transport eqation is modelled as: ρ ξ [ ] ( R ) m = Ym Y 3 m τ 1 ξ R m is a net rate of prodction of specie m by chemical reactions, * ξ ( ) is the length fraction of the fine scale trblent strctres where the reaction occrs, Y m is the fine scale species mass fraction (specie m) after reacting over the time, τ Y m is a species mass fraction for specie m in the srronding fine scales state.
10 The length fraction of fine strctres: ξ = η / SGS 1 3 where = μ ρ L ) and L SGS = 0.157V SGS t EDC combstion model 2/2 /( SGS 3 μ SGS and the Kolmogorov s velocity η = ρ LSGS Characteristic sb-grid eddy and Kolmogorov timescales are: 1 2 τ SGS = L SGS / SGS and τ μ L SGS = η 3 ρ SGS * The volme fraction of the fine scales is calclated as ξ 3 and species are assmed to react in the fine strctres over a time scale τ = τ η. It is assmes that all the fine scales in the cell are perfectly * stirred reactors with a residence time τ 1 4
11 Chemistry mechanism Mechanism for H2 combstion (with NO), 13 species and 37 reactions Units are cm3, mol, s, KJ and K HO2=OHO 2.E H2O=OHH 1.8E H2OO=OHOH 5.9E H2OH=H2OH 1.17E HO2M=HO2M 2.3E H2/1./ H2O/6.5/ O2/0.4/ N2/0.4/ HHO2=OHOH 1.5E HHO2=H2O2 2.5E OHHO2=H2OO2 2.E HHM=H2M 1.8E H2/1./ H2O/6.5/ O2/0.4/ N2/0.4/ HOHM=H2OM 2.2E H2/1./ H2O/6.5/ O2/0.4/ N2/0.4/ HO2HO2=H2O2O2 2.E H2O2M=OHOHM 1.3E H2O2OH=H2OHO2 1.E OHO2=OHO2 2.E HHO2=OH2O 5.E HOM=OHM 6.2E OOM=O2M 6.17E H2O2H=H2OOH 1.E H2O2H=HO2H2 4.79E OOHM=HO2M 1.E H2O2=OHOH 1.7E ON2=NNO 1.82E ONO=NO2 3.8E HNO=NOH 2.63E NOM=NOM 3.98E N2M=NNM 3.72E N2OO=NONO 6.92E N2OO=N2O2 1.E N2ON=N2NO 1.E NHO2=NOOH 1.E N2OH=N2OH 7.6E HNOO=NOOH 5.01E HNOOH=NOH2O 1.26E NOHO2=HNOO2 2.E HNOHO2=NOH2O2 3.16E HNOH=NOH2 1.26E HNOM=HNOM 1.78E GUTHEIL, E., BALAKRISHNAN, G. & WILLIAMS, F. A. (1993) Strctre and etinction of hydrogen-air diffsion flames. IN PETERS, N. & ROGG, B. (Eds.) Redced kinetic mechanisms for applications in combstion systems. New York, Springer-Verlag
12 Brst disk opening sections Brst disk opening process was approimated by a 10 step process. Sections were open from 1 to 10 at times calclated according to the table in the net slide.
13 Brst disk opening time Following Spence & Woods (1964), the following formla was sed to obtain the rptre time of a brst disk: t = ρbd k p 1 2 where ρ the density of the diaphragm material (annealed copper), b and d are the thickness and the diameter of the diaphragm, k constant, 0.92 p driving pressre Opening times of sections for each of the simlations are listed below: Section Opening time, μs [1.35 MPa case] Opening time, μs [1.65 MPa case] Opening time, μs [2.43 MPa case] Opening time, μs [2.9 MPa case]
14 Instant vs finite opening time
15 Simlation reslts
16 P=29 bar Ignition initiated in the radial vent channel and etingished. Combstion is reinitiated in a nmber of spots otside the PRD. Concentration of hydrogen in these spots st before the ignition is in the range 29-36% by vol.
17 P=24.3 bar Ignition is initiated in the radial channel (closer to the pstream wall second shock reflection) similar to 29 bar case Reignition otside is not observed. Process took 10μs and was eperimentally registered ( ignition in private commnication)
18 P=16.5 bar Temperatre Ignition spot is located in the aial channel (not in a symmetry plane). The combstion propagates pstream (ntil oygen is consmed)! Iso-srfaces are: Ble 550 K, Green 1500 K, Red 2400 K
19 P=16.5 bar OH mole fraction Ignition is confirmed by a hydroyl mole fraction distribtion. Light sensor in eperiments was located otside and cold not have registered ignition. Iso-srfaces are: Ble , Green 0.002, Red 0.02
20 P=13.5 bar no ignition Hydroyl mole fraction is (not shown). Note: temperatre scale is different in this pictre.
21 Conclsions The LES model based on EDC with detailed kinetics and inertial brst disk was sccessflly applied to reprodce and eplain eperimental observations in T-shape channel (PRD mock-p). The lower pressre limit for spontaneos ignition in a T- shaped PRD is established: there is no ignition at storage pressre 13.5 bar; there is ignition followed by selfetinction at pressres 16.5 and 24.3 bar; and there is sstainable ignition (reignition otside the channel) for storage pressre 29 bar. The ignition is located in the are of the second reflection of initial shock and where the concentration of hydrogen in air is in the range 29-36% vol. The model can be sed as a contemporary tool for design of innovative pressre relief devices.
22 MSc in Hydrogen Safety Engineering (distance learning corse): 22
Lewis number and curvature effects on sound generation by premixed flame annihilation
Center for Trblence Research Proceedings of the Smmer Program 2 28 Lewis nmber and crvatre effects on sond generation by premixed flame annihilation By M. Talei, M. J. Brear AND E. R. Hawkes A nmerical
More informationOn the importance of horizontal turbulent transport in high resolution mesoscale simulations over cities. A. Martilli (CIEMAT, Spain), B. R.
On the importance of horizontal trblent transport in high resoltion mesoscale simlations over cities. A. Martilli (CIEMAT, Spain), B. R. Rotnno, P. Sllivan, E. G. Patton, M. LeMone (NCAR, USA) In an rban
More informationWall treatment in Large Eddy Simulation
Wall treatment in arge Edd Simlation David Monfort Sofiane Benhamadoche (ED R&D) Pierre Sagat (Université Pierre et Marie Crie) 9 novembre 007 Code_Satrne User Meeting Wall treatment in arge Edd Simlation
More informationCourse Outline. Boundary Layer Flashback Core Flow Flashback and Combustion Induced Vortex Breakdown
Corse Otline A) Introdction and Otlook B) Flame Aerodynamics and Flashback C) Flame Stretch, Edge Flames, and Flame Stabilization Concepts D) Distrbance Propagation and Generation in Reacting Flows E)
More informationSpring Semester 2011 April 5, 2011
METR 130: Lectre 4 - Reynolds Averaged Conservation Eqations - Trblent Flxes (Definition and typical ABL profiles, CBL and SBL) - Trblence Closre Problem & Parameterization Spring Semester 011 April 5,
More informationModeling Effort on Chamber Clearing for IFE Liquid Chambers at UCLA
Modeling Effort on Chamber Clearing for IFE Liqid Chambers at UCLA Presented by: P. Calderoni own Meeting on IFE Liqid Wall Chamber Dynamics Livermore CA May 5-6 3 Otline his presentation will address
More information5.1 Heat removal by coolant flow
5. Convective Heat Transfer 5.1 Heat removal by coolant flow Fel pellet Bond layer Cladding tbe Heat is transferred from the srfaces of the fel rods to the coolant. T Temperatre at center of fc fel pellet
More informationNATURAL CONVECTION No mechanical force to push the fluid pump, fan etc. No predefined fluid flowrate and velocity can t prescribe Reynolds
NATURA CONVECTION No mechanical force to psh the flid pmp, fan etc. No predefined flid flowrate and velocit can t prescribe Renolds nmber Flid moves as a reslt of densit difference Flid velocit established
More informationPROBLEMS
PROBLEMS------------------------------------------------ - 7- Thermodynamic Variables and the Eqation of State 1. Compter (a) the nmber of moles and (b) the nmber of molecles in 1.00 cm of an ideal gas
More informationBurning Rate, Kinetic Coupling, and Mechanism Reduction
Brning Rate, Kinetic Copling, and Mechanism Redction Yigang J Mechanical and Aerospace Engineering Princeton University 27 AFOSR MURI Kick-Off Meeting Generation of Comprehensive Srrogate Kinetic Models
More informationCFD-Simulation thermoakustischer Resonanzeffekte zur Bestimmung der Flammentransferfunktion
CFD-Simlation thermoakstischer Resonanzeffekte zr Bestimmng der Flammentransferfnktion Ator: Dennis Paschke Technische Universität Berlin Institt für Strömngsmechanik nd Technische Akstik FG Experimentelle
More informationReducing the Carbon Footprint of Transportation Fuels through Computational Modelling
Reducing the Carbon Footprint of Transportation Fuels through Computational Modelling http://www.conserve-energy-future.com Computational Modelling Group Founded in 1999 Head Prof. Currently 20 members
More informationEDEXCEL NATIONAL CERTIFICATE/DIPLOMA. PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 13 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS
EDEXCEL NATIONAL CERTIFICATE/DIPLOMA PRINCIPLES AND APPLICATIONS of FLUID MECHANICS UNIT 3 NQF LEVEL 3 OUTCOME 3 - HYDRODYNAMICS TUTORIAL - PIPE FLOW CONTENT Be able to determine the parameters of pipeline
More informationTurbulence and boundary layers
Trblence and bondary layers Weather and trblence Big whorls hae little whorls which feed on the elocity; and little whorls hae lesser whorls and so on to iscosity Lewis Fry Richardson Momentm eqations
More informationApplying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB
IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile
More informationTHE ACOUSTIC RESPONSE OF BURNER-STABILIZED PREMIXED FLAT FLAMES
Proceedings of the Combstion Institte, Volme 29, 2002/pp. 115 122 THE ACOUSTIC RESPONSE OF BURNER-STABILIZED PREMIXED FLAT FLAMES K. R. A. M. SCHREEL, R. ROOK and L. P. H. DE GOEY Eindhoven University
More informationMicroscale physics of fluid flows
Microscale physics of flid flows By Nishanth Dongari Senior Undergradate Department of Mechanical Engineering Indian Institte of Technology, Bombay Spervised by Dr. Sman Chakraborty Ot line What is microflidics
More informationMicroscopic Properties of Gases
icroscopic Properties of Gases So far we he seen the gas laws. These came from observations. In this section we want to look at a theory that explains the gas laws: The kinetic theory of gases or The kinetic
More informationNumerical Simulation of Three Dimensional Flow in Water Tank of Marine Fish Larvae
Copyright c 27 ICCES ICCES, vol.4, no.1, pp.19-24, 27 Nmerical Simlation of Three Dimensional Flo in Water Tank of Marine Fish Larvae Shigeaki Shiotani 1, Atsshi Hagiara 2 and Yoshitaka Sakakra 3 Smmary
More informationTwo identical, flat, square plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHADED areas.
Two identical flat sqare plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHAE areas. F > F A. A B F > F B. B A C. FA = FB. It depends on whether the bondary
More informationMODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION
MODELLING OF TURBULENT ENERGY FLUX IN CANONICAL SHOCK-TURBULENCE INTERACTION Rssell Qadros, Krishnend Sinha Department of Aerospace Engineering Indian Institte of Technology Bombay Mmbai, India 476 Johan
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationToward a multilevel QM/MM methodology for performing molecular dynamic simulations of complex reactive processes
Toward a multilevel QM/MM methodology for performing molecular dynamic simulations of complex reactive processes Michael R. Salazar Department of Chemistry Union University, Jackson, TN 38305 msalazar@uu.edu
More informationUncertainties of measurement
Uncertainties of measrement Laboratory tas A temperatre sensor is connected as a voltage divider according to the schematic diagram on Fig.. The temperatre sensor is a thermistor type B5764K [] with nominal
More informationρ u = u. (1) w z will become certain time, and at a certain point in space, the value of
THE CONDITIONS NECESSARY FOR DISCONTINUOUS MOTION IN GASES G I Taylor Proceedings of the Royal Society A vol LXXXIV (90) pp 37-377 The possibility of the propagation of a srface of discontinity in a gas
More informationRadiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction
Radiation Effects on Heat and Mass Transfer over an Exponentially Accelerated Infinite Vertical Plate with Chemical Reaction A. Ahmed, M. N.Sarki, M. Ahmad Abstract In this paper the stdy of nsteady flow
More informationMEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY
MEASUREMENT OF TURBULENCE STATISTICS USING HOT WIRE ANEMOMETRY Mrgan Thangadrai +, Atl Kmar Son *, Mritynjay Singh +, Sbhendra *, Vinoth Kmar ++, Ram Pyare Singh +, Pradip K Chatterjee + + Thermal Engineering,
More informationNumerical Model for Studying Cloud Formation Processes in the Tropics
Astralian Jornal of Basic and Applied Sciences, 5(2): 189-193, 211 ISSN 1991-8178 Nmerical Model for Stdying Clod Formation Processes in the Tropics Chantawan Noisri, Dsadee Skawat Department of Mathematics
More informationCHAPTER 8 CONVECTION IN EXTERNAL TURBULENT FLOW
CHAPTER 8 CONVECTION IN EXTERNAL TURBULENT FLOW 8.1 Introdction Common phsical phenomenon, bt comple Still relies on empirical data and rdimentar conceptal drawings Tremendos growth in research over last
More informationChapter 1: Differential Form of Basic Equations
MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)
More informationDILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS
Forth International Conference on CFD in the Oil and Gas, Metallrgical & Process Indstries SINTEF / NTNU Trondheim, Noray 6-8 Jne 005 DILUTE GAS-LIQUID FLOWS WITH LIQUID FILMS ON WALLS John MORUD 1 1 SINTEF
More information5. The Bernoulli Equation
5. The Bernolli Eqation [This material relates predominantly to modles ELP034, ELP035] 5. Work and Energy 5. Bernolli s Eqation 5.3 An example of the se of Bernolli s eqation 5.4 Pressre head, velocity
More informationTheoretical study on impingement heat transfer with single-phase free-surface slot jets
International Jornal of Heat and Mass Transfer 4 (25) 331 336 www.elsevier.com/locate/ijhmt Theoretical stdy on impingement heat transfer with single-phase free-srface slot jets Y.C. Chen, C.F. Ma *, M.
More informationUNIT V BOUNDARY LAYER INTRODUCTION
UNIT V BOUNDARY LAYER INTRODUCTION The variation of velocity from zero to free-stream velocity in the direction normal to the bondary takes place in a narrow region in the vicinity of solid bondary. This
More informationSIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS
THERMAL SCIENCE, Year 011, Vol. 15, No. 1, pp. 45-55 45 SIMULATION OF TURBULENT FLOW AND HEAT TRANSFER OVER A BACKWARD-FACING STEP WITH RIBS TURBULATORS b Khdheer S. MUSHATET Mechanical Engineering Department,
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationHybrid modelling and model reduction for control & optimisation
Hybrid modelling and model redction for control & optimisation based on research done by RWTH-Aachen and TU Delft presented by Johan Grievink Models for control and optimiation market and environmental
More informationAn Experimental Facility for Studying Flame Propagation at Elevated Pressures
Paper # 7LT-21 Topic: Laminar and Trblent Flames 8 th U. S. National Combstion Meeting Organized by the Western States Section of the Combstion Institte and hosted by the University of Utah May 19-22,
More informationMaterial Transport with Air Jet
Material Transport with Air Jet Dr. István Patkó Bdapest Tech Doberdó út 6, H-1034 Bdapest, Hngary patko@bmf.h Abstract: In the field of indstry, there are only a very few examples of material transport
More informationQuantitative Characterization of Pressure-related Turbulence Transport Terms using Simultaneous Nonintrusive Pressure and Velocity Measurement*
Qantitative Characterization of Pressre-related Trblence Transport Terms sing Simltaneos Nonintrsive Pressre and Velocity Measrement* Xiaofeng Li and Joseph Katz Department of Aerospace Engineering San
More informationA Model-Free Adaptive Control of Pulsed GTAW
A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,
More informationComputational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model
Comptational Flid Dynamics Simlation and Wind Tnnel Testing on Microlight Model Iskandar Shah Bin Ishak Department of Aeronatics and Atomotive, Universiti Teknologi Malaysia T.M. Kit Universiti Teknologi
More informationA Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane
Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani
More informationHOW TO USE CHEMKIN 경원테크
HOW TO USE CHEMKIN 경원테크 Agenda CHEMKIN-PRO vs. Old Chemkin II CHEMKIN-PRO Overview Advanced Feature: Particle Tracking Advanced Feature: Reaction Path Analysis Example : Turbulent jet with kinetics & mixing
More informationLift-Off and Blow-Out of Under-Expanded Hydrogen Jets: Experiments versus Simulations
Lift-Off and Blow-Out of Under-Expanded Hydrogen Jets: Experiments versus Simulations Shentsov V. 1, *, Sakatsume R. 2, Makarov D. 1, Takeno K. 2, Molkov V. 1 1 University of Ulster, Hydrogen Safety Engineering
More informationThe Numerical Simulation of Enhanced Heat Transfer Tubes
Aailable online at.sciencedirect.com Phsics Procedia 4 (01 70 79 01 International Conference on Applied Phsics and Indstrial Engineering The Nmerical Simlation of Enhanced Heat Transfer Tbes Li Xiaoan,
More informationPrandl established a universal velocity profile for flow parallel to the bed given by
EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient
More informationParametrization of turbulent fluxes in the outer layer
Parametriation of trblent flxes in the oter layer Irina Sand Oerie of models Blk models Local -closre -profile closre 0 order 1 st order ED/F closre non-local TE closre 1.5 th order Crrent closre in the
More informationProcess Modeling of Wellbore Leakage for GCS Risk Assessment
Process Modeling of Wellbore Leakage for CS Risk Assessment Crtis M. Oldenbrg Earth Sciences Division Lawrence Berkeley National Laboratory Contribtors: Leha Pan (LBNL) Bill Carey (LANL) IEAH- Modeling
More informationIJAET International Journal of Application of Engineering and Technology ISSN: Vol.1 No.1
IJAET International Jornal of Application of Engineering and Technology ISSN: 395-3594 Vol1 No1 ANALYSIS OF SUPERSONIC FLOWS IN THE E -LAVAL NOZZLE AT 1 INTO A SUENLY EXPANE UCT AT L/=WITH CAVITY ASPECT
More informationLarge Eddy Simulation Of Flow Past A Two-dimensional Hill
Large Eddy Simlation Of Flow Past A Two-dimensional Hill Sankara N.Vengadesan ) and Akihiko Nakayama ) ) Research Associate, Email: vengades@kobe-.ac.jp, ) Professor, Email: nakayama@kobe-.ac.jp Gradate
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationToward a multilevel QM/MM methodology for performing molecular dynamic simulations of complex reactive processes
Toward a multilevel QM/MM methodology for performing molecular dynamic simulations of complex reactive processes Michael R. Salazar Department of Chemistry Union University, Jackson, TN 38305 msalazar@uu.edu
More informationInternational Journal of Scientific & Engineering Research, Volume 5, Issue 3, March ISSN
International Jornal of Scientific & Engineering Research, Volme 5, Isse 3, March-4 83 ISSN 9-558 Doble Dispersion effects on free convection along a vertical Wavy Srface in Poros Media with Variable Properties
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 02-01 The Krylov accelerated SIMPLER method for incompressible flow C. Vik and A. Saghir ISSN 1389-6520 Reports of the Department of Applied Mathematical Analysis
More informationarxiv: v1 [physics.flu-dyn] 4 Sep 2013
THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv:1309.1728v1 physics.fl-dyn]
More informationFlood flow at the confluence of compound river channels
Rier Basin Management VIII 37 Flood flow at the conflence of compond rier channels T. Ishikawa 1, R. Akoh 1 & N. Arai 2 1 Department of Enironmental Science and Technology, Tokyo Institte of Technology,
More informationStudy of the diffusion operator by the SPH method
IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait
More information1 JAXA Special Pblication JAXA-SP-1-E Small-scale trblence affects flow fields arond a blff body and therefore it governs characteristics of cross-sec
First International Symposim on Fltter and its Application, 1 11 IEXPERIMENTAL STUDY ON TURBULENCE PARTIAL SIMULATION FOR BLUFF BODY Hiroshi Katschi +1 and Hitoshi Yamada + +1 Yokohama National University,
More informationResidual Currents Induced by Asymmetric Tidal Mixing in Weakly Stratified Narrow Estuaries
SEPTEMBER 21 C H E N G E T A L. 2135 Residal Crrents Indced by Asymmetric Tidal Mixing in Weakly Stratified Narrow Estaries PENG CHENG* AND ARNOLDO VALLE-LEVINSON Department of Civil Coastal Engineering,
More informationFundamentals of magnetohydrodynamics
Fndamentals of magnetohydrodynamics Part II Daniel Gómez 1, Email: dgomez@df.ba.ar Webpage: astro.df.ba.ar (1 Institto de Astronomía y Física del Espacio, CONICET, Argentina ( Departamento de Física, Universidad
More information2 nd Joint Summer School on Fuel Cell and Hydrogen Technology September 2012, Crete, Greece. Hydrogen fires
2 nd Joint Summer School on Fuel Cell and Hydrogen Technology 17 28 September 2012, Crete, Greece Hydrogen fires Sile Brennan (on behalf of the HySAFER group) Hydrogen Safety Engineering and Research Centre
More informationMomentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary
Momentm Eqation Interest in the momentm eqation: Qantification of proplsion rates esign strctres for power generation esign of pipeline systems to withstand forces at bends and other places where the flow
More informationCHEMICAL REACTION EFFECTS ON FLOW PAST AN EXPONENTIALLY ACCELERATED VERTICAL PLATE WITH VARIABLE TEMPERATURE. R. Muthucumaraswamy and V.
International Jornal of Atomotive and Mechanical Engineering (IJAME) ISSN: 9-8649 (int); ISSN: 18-166 (Online); Volme pp. 31-38 Jly-December 1 niversiti Malaysia Pahang DOI: http://dx.doi.org/1.158/ijame..11.11.19
More informationHigh speed analysis of high pressure combustion in a constant volume cell
High speed analysis of high pressre combstion in a constant volme cell.j.m. Frijters *, R.J.H. Klein-Dowel, S.S. Manski, L.M.T. Somers and R.S.G. Baert Section Combstion Technology Eindhoven University
More informationSources of Non Stationarity in the Semivariogram
Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationNumerical Simulation of Heat and Mass Transfer in the Evaporation of Glycols Liquid Film along an Insulated Vertical Channel
Colmbia nternational Pblishing American Jornal of Heat and Mass Transfer doi:1.7726/ajhmt.215.15 Research Article Nmerical Simlation of Heat and Mass Transfer in the Evaporation of lcols iqid Film along
More informationFLUID FLOW FOR CHEMICAL ENGINEERING
EKC FLUID FLOW FOR CHEMICL ENGINEERING CHTER 8 (SOLUTION WI EXERCISE): TRNSORTTION SYSTEM & FLUID METERING Dr Mohd zmier hmad Tel: +60 (4) 5996459 Email: chazmier@eng.sm.my . Benzene at 7.8 o C is pmped
More informationProceedings of the ASME nd International Conference on Ocean, Offshore and Arctic Engineering OMAE2013 June 9-14, 2013, Nantes, France
Proceedings of the ASME 203 32nd International Conference on Ocean, Offshore and Arctic Engineering OMAE203 Jne 9-4, 203, Nantes, France OMAE203-0207 THE EFFICIENCY OF TIDAL FENCES: A BRIEF REVIEW AND
More informationFLUID MECHANICS D203 SAE SOLUTIONS TUTORIAL 2 APPLICATIONS OF BERNOULLI SELF ASSESSMENT EXERCISE 3
FLUI MEHNIS 0 SE SOLUTIONS TUTORIL PPLITIONS OF ERNOULLI SELF SSESSMENT EXERISE Take the density of water to be 997 kg/m throghot nless otherwise stated.. Ventri meter is 50 mm bore diameter at inlet and
More informationThe Bow Shock and the Magnetosheath
Chapter 6 The Bow Shock and the Magnetosheath The solar wind plasma travels sally at speeds which are faster than any flid plasma wave relative to the magnetosphere. Therefore a standing shock wave forms
More informationM. Moeini 1 and M. R. Chamani 2. Department of Civil Engineering, Amirkabir University of Technology, Tehran, Iran 2
Jornal of Applied Flid Mechanics, Vol. 1, No. 4, pp. 171-177, 17. Available online at www.jafmonline.net, ISSN 1735-357, EISSN 1735-3645. DOI: 1.18869/acadpb.jafm.73.41.7544 New Perspectives on the Laminar
More informationNumerical Study on Bouncing and Separation Collision Between Two Droplets Considering the Collision-Induced Breakup
Jornal of Mechanical Science and Technology (007) 585~59 Jornal of Mechanical Science and Technology Nmerical Stdy on Boncing and Separation Collision Between Two Droplets Considering the Collision-Indced
More informationLift-off and blow-out of under-expanded hydrogen jets: experiments versus simulations
Lift-off and blow-out of under-expanded hydrogen jets: experiments versus simulations Shentsov, V. 1, Sakatsume, R. 2, Makarov, D. 1, Takeno, K. 2, Molkov, V. 1 1 Hydrogen Safety Engineering and Research
More informationProceedings of the 51st Anniversary Conference of KSME PHYSICAL MODELING OF ATMOSPHERIC FLOW AND ENVIRONMENTAL APPLICATIONS B. R.
Proceedings of the 51st Anniversary Conference of KSME PHYSICAL MODELING OF ATMOSPHERIC FLOW AND ENVIRONMENTAL APPLICATIONS B. R. WHITE Faclty of Mechanical and Aeronatical Engineering University of California,
More informationScienceDirect. River morphology modeling at the downstream of Progo River post eruption 2010 of Mount Merapi
Available online at www.sciencedirect.com ScienceDirect Procedia Environmental Sciences 8 (015 ) 148 157 The 5th Sstainable Ftre for Hman Secrity (SstaiN 014) River morphology modeling at the downstream
More informationThermal balance of a wall with PCM-enhanced thermal insulation
Thermal balance of a wall with PCM-enhanced thermal inslation E. Kossecka Institte of Fndamental Technological esearch of the Polish Academy of Sciences, Warsaw, Poland J. Kośny Oak idge National aboratory;
More informationSteady State and Transient Thermal Analysis of Switched Reluctance Machine
Steady State and Transient Thermal Analysis of Switched Relctance Machine E. Annie Elisabeth Jebaseeli and S. Paramasivam Abstract This paper presents the two dimensional (-D) steady state and transient
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationActive Flux Schemes for Advection Diffusion
AIAA Aviation - Jne, Dallas, TX nd AIAA Comptational Flid Dynamics Conference AIAA - Active Fl Schemes for Advection Diffsion Hiroaki Nishikawa National Institte of Aerospace, Hampton, VA 3, USA Downloaded
More informationHomotopy Perturbation Method for Solving Linear Boundary Value Problems
International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop
More informationBy Dr. Salah Salman. Problem (1)
Chemical Eng. De. Problem ( Solved Problems Samles in Flid Flow 0 A late of size 60 cm x 60 cm slides over a lane inclined to the horizontal at an angle of 0. It is searated from the lane with a film of
More informationAppendix A: The Fully Developed Velocity Profile for Turbulent Duct Flows
Appendix A: The lly Developed Velocity Profile for Trblent Dct lows This appendix discsses the hydrodynamically flly developed velocity profile for pipe and channel flows. The geometry nder consideration
More informationVIBRATION MEASUREMENT UNCERTAINTY AND RELIABILITY DIAGNOSTICS RESULTS IN ROTATING SYSTEMS
VIBRATIO MEASUREMET UCERTAITY AD RELIABILITY DIAGOSTICS RESULTS I ROTATIG SYSTEMS. Introdction M. Eidkevicite, V. Volkovas anas University of Technology, Lithania The rotating machinery technical state
More informationExercises in Combustion Technology
Exercises in Combustion Technology Exercise 4: Turbulent Premixed Flames Turbulent Flow: Task 1: Estimation of Turbulence Quantities Borghi-Peters diagram for premixed combustion Task 2: Derivation of
More informationA GNF framework for turbulent flow models of drag reducing fluids and proposal for a k ε type closure
J. Non-Newtonian Flid Mech. 114 2003 149 184 A GNF framewor for trblent flow models of drag redcing flids and proposal for a ε type closre F. T. Pinho Centro de Estdos de Fenómenos de Transporte, DEMEGI,
More informationPhysicsAndMathsTutor.com
. Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the
More informationThe SISTEM method. LOS ascending
The SISTEM method Simltaneos and Integrated Strain Tensor Estimation from geodetic and satellite deformation Measrements A new global approach to obtain three-dimensional displacement maps by integrating
More informationTransient Approach to Radiative Heat Transfer Free Convection Flow with Ramped Wall Temperature
Jornal of Applied Flid Mechanics, Vol. 5, No., pp. 9-1, 1. Available online at www.jafmonline.net, ISSN 175-57, EISSN 175-645. Transient Approach to Radiative Heat Transfer Free Convection Flow with Ramped
More informationReduction of over-determined systems of differential equations
Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical
More informationChapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS
Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system
More informationEntropy ISSN
Entrop 3, 5, 56-518 56 Entrop ISSN 199-43 www.mdpi.org/entrop/ Entrop Generation Dring Flid Flow Between wo Parallel Plates With Moving Bottom Plate Latife Berrin Erba 1, Mehmet Ş. Ercan, Birsen Sülüş
More informationAnalysis of Crack Propagation in Asphalt Concrete Using Cohesive Crack Model
90 TRANSPORTATION RESEARCH RECORD 1317 Analysis of Crack Propagation in Asphalt Concrete Using Cohesive Crack Model YEOU-SHANG }ENQ AND }IA-DER PERNG A cohesive crack model, which is similar to the Dgdale
More informationTHE REDUCTION IN FINESTRUCTURE CONTAMINATION OF INTERNAL WAVE ESTIMATES FROM A TOWED THERMISTOR CHAIN
DANIEL C. DUBBEL THE REDUCTION IN FINESTRUCTURE CONTAMINATION OF INTERNAL WAVE ESTIMATES FROM A TOWED THERMISTOR CHAIN Estimates of internal wave displacements based on towed thermistor array data have
More informationV. Hadron quantum numbers
V. Hadron qantm nmbers Characteristics of a hadron: 1) Mass 2) Qantm nmbers arising from space-time symmetries : total spin J, parity P, charge conjgation C. Common notation: 1 -- + 2 J P (e.g. for proton:
More informationUncertainty Evaluation of Toluene Determination in Room Air by Thermal Desorption Gas Chromatography
International Conference on Civil, Transportation and Environment (ICCTE 06) ncertainty Evalation of Tolene Determination in Room Air by Thermal Desorption Gas Chromatography Xiaoyan Wen, a,yanhi Gao,
More informationAvailable online at ScienceDirect. Procedia Engineering 150 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 150 (016 ) 635 640 International Conference on Indstrial Engineering, ICIE 016 Research of Dynamics of the Rotor with Three-film
More informationChapter 6 Momentum Transfer in an External Laminar Boundary Layer
6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned
More information