Impurities in inelastic Maxwell models

Size: px
Start display at page:

Download "Impurities in inelastic Maxwell models"

Transcription

1 Impurities in inelastic Maxwell moels Vicente Garzó Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain Abstract. Transport properties of impurities immerse in a granular gas unergoing homogenous cooling state are stuie. The results are obtaine from the Boltzmann-Lorentz kinetic equation for inelastic Maxwell moels in imensions. The kinetic equation is solve by means of the Chapman-Enskog metho up to first orer in the ensity graient of impurities. The mass an heat fluxes associate with impurities are etermine an the corresponing transport coefficients are ientifie. Both transport coefficients (iffusion an Dufour coefficients) are exactly obtaine in terms of the coefficients of restitution for the impurity-gas an gas-gas collisions as well as the ratios of mass an iameters. The results are compare with those obtaine for inelastic har spheres in the leaing Sonine approximation an by means of Monte Carlo simulations. The comparison shows goo agreement between both interaction moels, especially in the case of the iffusion coefficient. 1. INTRODUCTION Granular fluis are usually moelle by an iealize system of smooth har spheres with inelastic collisions. In the low-ensity limit, the Boltzmann equation has been conveniently moifie to account for the inelasticity of binary collisions an the Navier-Stokes transport coefficients have been compute by solving the corresponing kinetic equation by means of the Chapman-Enskog metho [1]. This task has been carrie out in the past few years in the case of a monocomponent gas [] as well as for multicomponent systems [3]. In both cases, as happens for elastic collisions, all the transport coefficients are given in terms of the solutions to linear integral equations which are solve approximately by taking the leaing terms in a Sonine polynomial expansion. In spite of this approximation, the kinetic theory results compare quite well with those obtaine from numerical solutions of the Boltzmann equation by using the irect simulation Monte Carlo (DSMC) metho [4]. The fact that the agreement between kinetic theory an DSMC results even occurs for strong issipation is a further testimony of the utility of the kinetic theory an the hyroynamic escriptions for granular fluis beyon the weak issipation limit. The main mathematical ifficulty in solving the Boltzmann equation for inelastic har spheres (IHS) comes from the form of the collision rate, which is proportional to the magnitue of the relative velocity of the two colliing particles. In the case of elastic collisions, an alternative to overcome this problem is to consier the repulsive Maxwell potential, for which the collision rate is inepenent of the relative velocity. This property simplifies the mathematical structure of the Boltzmann collision operator an allows one to get exact solutions of the Boltzmann equation for far from equilibrium situations [5]. In the case of inelastic gases, the so-calle inelastic Maxwell moels (IMM) [6] have been wiely use in the past few years as a toy moel to characterize the effect of the inelasticity of collisions on the properties of granular gases. The IMM share with elastic Maxwell molecules the property that the collision rate is velocity inepenent but their scattering rules are the same as for IHS. Most of the stuies performe in the context of IMM have been evote to homogeneous states, mainly to the analysis of the overpopulate high energy tails [7]. Much less is known for inhomogeneous situations, especially for multicomponent systems. One of the few stuies consiering mixtures of IMM has been carrie out by the author [8] in the uniform shear flow problem. The rheological properties (shear an normal stresses) of the mixture have been exactly evaluate in terms of issipation an parameters of the mixture an surprisingly, the results of IMM show an excellent agreement with those obtaine for IHS by means of Gra s approximation an Monte Carlo simulations. This fact stimulates the search for exact solutions for IMM which can be confronte with the results obtaine for IHS by using approximate analytical methos an computer simulations. The goal of this paper is to etermine the iffusion coefficient of impurities immerse in a gas of inelastic Maxwell molecules. The gas is in the homogeneous cooling state (HCS), namely, a homogeneous state where the temperature uniformly ecreases in time ue to collisional cooling. Diffusion is generate in the system by the presence of a 83

2 Report Documentation Page Form Approve OMB No Public reporting buren for the collection of information is estimate to average 1 hour per response, incluing the time for reviewing instructions, searching existing ata sources, gathering an maintaining the ata neee, an completing an reviewing the collection of information. Sen comments regaring this buren estimate or any other aspect of this collection of information, incluing suggestions for reucing this buren, to Washington Heaquarters Services, Directorate for Information Operations an Reports, 115 Jefferson Davis Highway, Suite 14, Arlington VA -43. Responents shoul be aware that notwithstaning any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it oes not isplay a currently vali OMB control number. 1. REPORT DATE 13 JUL 5 4. TITLE AND SUBTITLE Impurities in inelastic Maxwell moels. REPORT TYPE N/A 3. DATES COVERED - 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 1. SPONSOR/MONITOR S ACRONYM(S) 1. DISTRIBUTION/AVAILABILITY STATEMENT Approve for public release, istribution unlimite 11. SPONSOR/MONITOR S REPORT NUMBER(S) 13. SUPPLEMENTARY NOTES See also ADM179, International Symposium on Rarefie Gas Dynamics (4th) Hel in Monopoli (Bari), Italy on 1-16 July ABSTRACT 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT UU a. REPORT unclassifie b. ABSTRACT unclassifie c. THIS PAGE unclassifie 18. NUMBER OF PAGES 6 19a. NAME OF RESPONSIBLE PERSON Stanar Form 98 (Rev. 8-98) Prescribe by ANSI St Z39-18

3 weak concentration graient, which also inuces a heat flux. Both mass an heat fluxes efine the relevant transport coefficients of the problem: the iffusion coefficient an the Dufour coefficient. Here, I get explicit expressions for these coefficients by solving the corresponing Boltzmann-Lorentz equation for impurities by means of the Chapman- Enskog metho through the first orer in the concentration graient. These transport coefficients are given in terms of the coefficients of restitution for the impurity-gas an gas-gas collisions as well as the masses an iameters of impurities an gas particles. The epenence of these coefficients on the parameters of the system is illustrate an compare with known results erive for IHS [3, 9].. DESCRIPTION OF THE PROBLEM Let us consier a ilute gas of inelastic Maxwell molecules of mass m, iameter σ, an interparticle coefficient of restitution α. In the low-ensity regime, its velocity istribution function f (v) obeys the nonlinear Boltzmann equation (BE). We assume that the gas is in the HCS, namely, a homogeneous solution of the BE in which all the time epenence of f (v,t) is through the granular temperature T (t). In this case, the BE can be written as 1 ζ (v f ) = J[v f, f ], (1) v where the Boltzmann collision operator J[v f, f ] for IMM is J [v 1 f, f ] = ω v σ [ α 1 f (v nω 1,t) f (v,t) f (v 1,t) f (v,t) ]. () Here, n is the number ensity of the gas particles, ω is an effective collision frequency (to be chosen later), Ω = π / /Γ(/) is the total soli angle in imensions, an σ is a unit vector irecte along the centers of the two colliing spheres. Upon writing Eq. (1), I have taken into account the balance equation for the temperature T 1 t T = ζ, where the cooling rate ζ is given by [1] ζ (α) = 1 α ω. (3) The explicit form of the velocity istribution function f (v, t) is not explicitly known, although its velocity moments can be explicitly compute [1]. We assume now that some impurities of mass m an iameter σ are ae to the system. Given that its molar fraction is negligible, the state of the granular gas is not isturbe by the presence of impurities an so the velocity istribution function f (v, t) obeys the BE (1). Furthermore, collisions among impurities themselves can neglecte versus the impurity-gas collisions, which are characterize by a coefficient of restitution α. As sai in the Introuction, I want to etermine the transport properties of impurities generate in the system by the presence of a weak concentration graient x, x = n /n being the mole fraction of impurities. Uner these conitions, the velocity istribution function f (r,v,t) of impurities verifies the Boltzmann-Lorentz equation t f + v f = J[v f, f ], (4) where J [v 1 f, f ] = ω v σ [ α 1 f (r,v n Ω 1;t) f (r,v ;t) f (r,v 1 ;t) f (r,v ;t) ], (5) ω is an effective collision frequency (to be chosen also later) for impurity-gas collisions, an v 1 = v 1 m ( ) 1 + α 1 ( σ g 1 ) σ, v = v + m ( ) 1 + α 1 ( σ g 1 ) σ, (6) m + m m + m where g 1 = v 1 v is the relative velocity of the colliing pair. The partial temperature of impurities T is efine as n T = v m v f (v), (7) where n is the number ensity of impurities. In contrast to what happens for normal fluis, if impurities an gas particles are mechanically ifferent, the partial temperature T can iffer from the granular temperature T when 84

4 issipation is present. The failure of energy equipartition in granular fluis has even observe in real experiments [11]. Momentum an energy are not collisional invariants of the Boltzmann-Lorentz collision operator J[ f, f ]. Only the number ensity of impurities is conserve: where the mass flux j is efine as t n + j m =, (8) j = vm v f (v). (9) In the absence of iffusion, impurities are also in HCS an Eq. (4) becomes 1 ζ v (v f ) = J[v f, f ], (1) where use has been mae of the relation T 1 t T = ζ, where ζ is the cooling rate associate with the partial temperature T of impurities. Its expression is [8] ζ = ω [ µ(1 + α ) 1 µ ] (1 + α )(1 + θ), (11) where µ = m/(m + m ) an θ = m T /mt is the mean square velocity of gas particles relative to that of impurities. The balance equations for the temperatures T an T lea to the following time evolution equation for the temperature ratio γ T /T : γ 1 t γ = ζ ζ. (1) The fact that f epens on time only through T (t) necessarily implies that the temperature ratio γ must be inepenent of time, an so Eq. (1) gives the HCS conition ζ (t) = ζ (t). The impurity equilibrates to a common HCS with ifferent temperatures for the impurity an gas particles. This implies a breakown of the energy equipartition. In orer to get the explicit epenence of γ on the parameter space, one still nees to fix the parameters ω an ω. To optimize the agreement with the IHS results, the cooling rates ζ an ζ of IMM given by Eqs. (3) an (11), respectively, are ajuste to be the same as the ones given for IHS [1]. Since ζ an ζ are not exactly known for IHS, one can estimate them by taking their local equilibrium approximations. In this case, the collision frequencies ω an ω are given by ω = π Ω ν, ω = Ω π ( σ σ ) 1 ( ) 1 + θ 1/ ν, (13) θ where σ = (σ +σ )/ an ν = nσ 1 T /m is an effective collision frequency. With these choices, the epenence of the temperature ratio T /T on the parameters of the system given for IMM presents an excellent agreement with the one foun for IHS, even for strong issipation [13]. 3. TRANSPORT PROPERTIES OF IMPURITIES The aim of this Section is to get the transport properties (mass an heat fluxes) of impurities in the limit of small concentration graient x. These fluxes are etermine by solving the Boltzmann-Lorentz equation (4) from the Chapman-Enskog metho [1] aapte to issipative systems. Thus, I look for a normal solution in which all the space an time epenence of f is through the hyroynamic fiels. In our problem, this means that f (r,v,t) = f [v x (t),t (t)]. The normal solution is generate by expaning f in powers of x as f = f () + ε f (1) +, (14) where each factor ε correspons to the implicit factor x. The time erivative is also expane as t = () t +ε (1), where the action of the ifferent operators (k) t can be obtaine from the balance equation (8). They are given by () t x =, () t T = T ζ, (k) t t + x = j(k 1), (k) t T =, k 1, (15) m n 85

5 D(α) /D(1) m /m=8 σ /σ= α FIGURE 1. Plot of the reuce iffusion coefficient D(α)/D(1) as a function of the coefficient of restitution α = α in the three-imensional case for σ /σ = an m /m = 8. The soli line correspons to the exact result obtaine here for IMM while the ashe line is the result erive for IHS in the secon Sonine approximation.the symbols refer to Monte Carlo simulations for IHS. where j (k) = vm v f (k). (16) The zeroth-orer approximation f () is the solution of Eq. (1) but taking into account now the local epenence on the mole fraction x. Since f () is isotropic, it follows that the flux of impurities vanishes at this orer, i.e., j () an so (1) t x =. To first orer in ε, one has the kinetic equation () t f (1) J[ f (1) (, f ] = = (1) t ( f () x ) + v 1 f () =, ) v 1 x. (17) The secon equality follows from the balance equations (15) an the space epenence of f () through x. The mass flux j (1) an the contribution to the heat flux coming from impurities q (1) can be easily obtaine from Eq. (17). The heat flux is efine as q (1) = v m v v f (1). (18) To get these fluxes, one multiplies both sies of Eq. (17) by m v an 1 m v v an integrates over v. The result for the mass an heat fluxes is, respectively ( ) () t + ν j (1) = γ p x, (19) where p = nt is the hyrostatic pressure an ( ) () t + ν q q (1) = βj (1) pt ( + m c = 8 [ m ( + ) 4n T vv 4 f (v) 1 + c ) γ x, () ] ( + ). (1) 4 Upon writing Eqs. (19) an (), use has been mae of the collisional moments [13] vm vj[ f (1), f ] = ν j (1), () v m v vj[ f (1), f ] = ν q q (1) βj (1), (3) 86

6 L q (α) m /m=8 σ =σ 6 4 m /m= FIGURE. Plot of the reuce Dufour coefficient L q (α) as a function of the coefficient of restitution α = α in the threeimensional case for σ = σ, an two ifferent values of the mass ratio: m /m = an m /m = 8. The soli lines correspon to the exact results obtaine here for IMM while the ashe lines are the results erive for IHS in the first Sonine approximation. α where β = ν The mass an heat fluxes have the forms ν = Ω ( ) σ 1 ( ) 1 + θ 1/ µ ν, (4) π σ θ ν q = ν + {µ(1 + α )[ + 8 3µ(1 + α )] 3( + )}, (5) T m µ(1 + α )[ 3µ(1 + α ) + ]. (6) j (1) = m D x, (7) q (1) = T D q x, (8) where D is the iffusion coefficient an D q is a cross transport coefficient measuring energy transport ue to iffusion of impurities. Dimensional analysis requires that D T 1/ an D q T 1/. Consequently, () t j (1) = ζ T T j (1) = m ζ D x, () t q (1) = ζ T T q (1) = 3 ζ T D q x. (9) Inserting Eqs. (9) into Eqs. (19) an (), respectively, the transport coefficients D an D q can be easily ientifie from (7) an (8). The result is D = p γ m ν ν 1 ζ, (3) D q = n (m /m)β D + (1 + c )γ m ν 3 ζ. (31) νq Here, ν = ν /ν, β = βm/(t ν), ν q = ν q /ν, an ζ = ζ /ν. A useful way of characterizing the heat transport in a mixture ue to a concentration graient is through the Dufour coefficient. In our problem, the (imensionless) Dufour coefficient L q is efine by the relation J q q (1) + T j(1) = nt m m ν L q x. (3) The coefficient L q can be easily written in terms of the coefficients D an D q as L q = m ( ν T D q + ) p D. (33) 87

7 In the elastic case (α = α = 1), ζ =, T = T, an as expecte L q =. Figures 1 an show the epenence of the reuce coefficients D(α)/D(1) an L q (α), respectively, on the (common) coefficient of restitution α = α for ifferent systems. Here, D(1) is the elastic value of the iffusion coefficient. In the case of the iffusion coefficient, it is apparent that the results for IMM an IHS are practically inistinguishable over the whole range of values of α stuie. Moreover, both theories present an excellent agreement with Monte Carlo simulations. More iscrepancies between IMM an IHS are foun in the case of the Dufour coefficient L q, especially with increasing issipation. However, at a qualitative level, the IMM preictions compare well with the ones erive for IHS. 4. CONCLUDING REMARKS Diffusion of impurities in a ilute granular gas of IMM unergoing HCS has been analyze. The relevant transport coefficients of the problem (iffusion an Dufour coefficients) have been etermine from a Chapman-Enskog solution [1] of the Boltzmann-Lorentz equation up to first orer in the concentration graient. The expressions (3), (31), an (33) erive here for the iffusion D an Dufour L q coefficients are exact (within the context of IMM) an constitute the main goal of this paper. This contrasts with the previous results erive for IHS [3, 9], where the coefficients D an L q were obtaine by consiering the leaing terms in a Sonine polynomial expansion of the velocity istribution function of impurities. To make contact with the results reporte for IHS, the collision frequencies ω an ω (which can be seen as free parameters of the moel) nee to be chosen. Here, I have taken these frequencies to reprouce the cooling rates ζ an ζ of IHS in the local equilibrium approximation. With this choice, the comparison with known results for IHS shows that in general, the IMM preictions are reasonably goo for not too large issipation, especially for the iffusion coefficient which is relate to a first-egree velocity moment of the velocity istribution function. In the case of the Dufour coefficient, although the IMM capture qualitatively well the trens obtaine for IHS, the isagreement between both interaction potentials is more significant than the one foun for the iffusion coefficient D. In this sense, more iscrepancies between IMM an IHS are expecte when one consiers higher egree velocity moments. ACKNOWLEDGMENTS Partial support from the Ministerio e Eucación y Ciencia (Spain) through Grant No. FIS is acknowlege. REFERENCES 1. Chapman, S., an Cowling, T. G., The Mathematical Theory of Nonuniform Gases (Cambrige University, Cambrige, 197).. Brey, J. J., Dufty, J. W., Kim, C. S., an Santos, A., Phys. Rev. E 58, (1998). 3. Garzó, V., an Dufty, J. W., Phys. Fluis 14, (). 4. Bir G. A., Molecular Gas Dynamics an the Direct Simulation Monte Carlo of Gas Flows (Clarenon, Oxfor, 1994). 5. Garzó, V., an Santos, A., Kinetic Theory of Gases in Shear Flows. Nonlinear Transport (Kluwer Acaemic, Dorrecht, 3). 6. Bobylev, A. V., Carrillo, J. A., an Gamba, I., J. Stat. Phys. 98, (); Ernst, M. H., an Brito, R., J. Stat. Phys. 19, (). 7. See for instance, Ernst, M. H., an Brito, R., Phys. Rev. E 65, (); Krapivsky, P. L., an Ben-Naim, E., Phys. Rev. E 66, (); Bobylev, A. V., an Cercignani, C., J. Stat. Phys. 11, (3). 8. Garzó, V., J. Stat. Phys. 11, (3). 9. Garzó, V., an Montanero, J. M., Phys. Rev. E 69, (4). 1. Santos, A., Physica A 31, (3). 11. Wilman, R. D., an Parker, D. J., Phys. Rev. Lett. 88, (); Feitosa, K., an Menon, N., Phys. Rev. Lett. 88, (). 1. Garzó, V., an Dufty, J. W., Phys. Rev. E 6, (1999). 13. Garzó, V., an Astillero, A., e-print con-mat/

BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES

BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES Vicente Garzó Departamento de Física, Universidad de Extremadura Badajoz, SPAIN Collaborations Antonio Astillero, Universidad de Extremadura José

More information

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 21 Nov 2003

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 21 Nov 2003 Diffusion of impurities in a granular gas Vicente Garzó Departamento de Física, Universidad de Extremadura, E-0607 Badajoz, Spain arxiv:cond-mat/0309639v [cond-mat.stat-mech] Nov 003 José María Montanero

More information

A wall-function approach to incorporating Knudsen-layer effects in gas micro flow simulations

A wall-function approach to incorporating Knudsen-layer effects in gas micro flow simulations A wall-function approach to incorporating Knudsen-layer effects in gas micro flow simulations D. A. Lockerby 1, J. M. Reese 2 and M. A. Gallis 3 1 Department of Mechanical Engineering, King s College London,

More information

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications Joseph R. Gabriel Naval Undersea Warfare Center Newport, Rl 02841 Steven M. Kay University of Rhode

More information

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE A wavelet-based generalization of the multifractal formalism from scalar to vector valued d- dimensional random fields : from theoretical concepts to experimental applications P. Kestener and A. Arneodo

More information

Report Documentation Page

Report Documentation Page Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems Aerodynamic Issues of Unmanned Air Vehicles Fluid-Structure Interaction High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems Raymond E. Gordnier Computational Sciences

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

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 information

Diagonal Representation of Certain Matrices

Diagonal Representation of Certain Matrices Diagonal Representation of Certain Matrices Mark Tygert Research Report YALEU/DCS/RR-33 December 2, 2004 Abstract An explicit expression is provided for the characteristic polynomial of a matrix M of the

More information

Report Documentation Page

Report Documentation Page Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

ADIT 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 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 information

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO William L. Dimpfl Space Science Applications Laboratory The Aerospace Corporation

More information

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA A report (dated September 20, 2011) on scientific research carried out under Grant: FA2386-10-1-4150 First-principles determination of thermal properties in nano-structured hexagonal solids with doping

More information

System Reliability Simulation and Optimization by Component Reliability Allocation

System Reliability Simulation and Optimization by Component Reliability Allocation System Reliability Simulation and Optimization by Component Reliability Allocation Zissimos P. Mourelatos Professor and Head Mechanical Engineering Department Oakland University Rochester MI 48309 Report

More information

Thermo-Kinetic Model of Burning for Polymeric Materials

Thermo-Kinetic Model of Burning for Polymeric Materials Thermo-Kinetic Model of Burning for Polymeric Materials Stanislav I. Stoliarov a, Sean Crowley b, Richard Lyon b a University of Maryland, Fire Protection Engineering, College Park, MD 20742 b FAA W. J.

More information

Hypersonic Flow over a Wedge with a Particle Flux Method

Hypersonic Flow over a Wedge with a Particle Flux Method Hypersonic Flow over a Wedge with a Particle Flux Method M. N. Macrossan, M. V. Metchnik and P. A. Pinto Centre for Hypersonics, School of Engineering, University of Queensland, St Lucia, 47, Australia

More information

Kinetic Theory for Binary Granular Mixtures at Low Density

Kinetic Theory for Binary Granular Mixtures at Low Density 10 Kinetic Theory for Binary Granular Mixtures at Low Density V. Garzó Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain vicenteg@unex.es Many features of granular media can be

More information

Dynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University

Dynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University Dynamics of Droplet-Droplet and Droplet-Film Collision C. K. Law Princeton University The physical phenomena of droplet-droplet and droplet-film collision in the head-on orientation were studied experimentally

More information

Introduction to the Vlasov-Poisson system

Introduction 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 information

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E. Concept V- Hull Floor Configuration in Two Dimensions Dr. Bijan Khatib-Shahidi & Rob E. Smith 10 November 2010 : Dist A. Approved for public release Report Documentation Page Form Approved OMB No. 0704-0188

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE 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 information

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT) AL/EQ-TP-1993-0307 REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT) John T. Mourand, John C. Crittenden, David W. Hand, David L. Perram, Sawang Notthakun Department of Chemical Engineering

More information

DSMC Calculation of Vortex Shedding behind a Flat Plate with a New Intermolecular Collision Scheme

DSMC Calculation of Vortex Shedding behind a Flat Plate with a New Intermolecular Collision Scheme DSMC Calculation of Vortex Shedding behind a Flat Plate with a New Intermolecular Collision Scheme M. Usami and K. Mizuguchi Dept. of Mechanical Eng., Mie University, 1515 Kamihama-cho, Tsu 514-8507, Japan

More information

Attribution Concepts for Sub-meter Resolution Ground Physics Models

Attribution Concepts for Sub-meter Resolution Ground Physics Models Attribution Concepts for Sub-meter Resolution Ground Physics Models 76 th MORS Symposium US Coast Guard Academy Approved for public release distribution. 2 Report Documentation Page Form Approved OMB No.

More information

Flow of a Rarefied Gas between Parallel and Almost Parallel Plates

Flow of a Rarefied Gas between Parallel and Almost Parallel Plates Flow of a Rarefied Gas between Parallel and Almost Parallel Plates Carlo Cercignani, Maria Lampis and Silvia Lorenzani Dipartimento di Matematica, Politecnico di Milano, Milano, Italy 033 Abstract. Rarefied

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal 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 information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data.

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data. Coherence of Low-Frequency Sound Signals Propagating through a Fluctuating Ocean: Analysis and Theoretical Interpretation of 2004 NPAL Experimental Data Alexander G. Voronovich NOAA/ESRL, PSD4, 325 Broadway,

More information

Report Documentation Page

Report Documentation Page Inhibition of blood cholinesterase activity is a poor predictor of acetylcholinesterase inhibition in brain regions of guinea pigs exposed to repeated doses of low levels of soman. Sally M. Anderson Report

More information

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects Chris Harrison Centre for Marine

More information

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model for Dopant and Impurity Segregation During Vapor Phase Growth Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016 Noname manuscript No. (will be inserte by the eitor) Scaling properties of the number of ranom sequential asorption iterations neee to generate saturate ranom packing arxiv:607.06668v2 [con-mat.stat-mech]

More information

Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method

Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method Manohar D. Deshpande, NASA Langley Research Center, Hmapton, VA., USA Abstract

More information

A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS

A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS Michael Thiel*, Truong Vu Bang Giang and Achim Dreher Institute of

More information

Testing Turbulence Closure Models Against Oceanic Turbulence Measurements

Testing Turbulence Closure Models Against Oceanic Turbulence Measurements Testing Turbulence Closure Models Against Oceanic Turbulence Measurements J. H. Trowbridge Woods Hole Oceanographic Institution Woods Hole, MA 02543 phone: 508-289-2296 fax: 508-457-2194 e-mail: jtrowbridge@whoi.edu

More information

Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight

Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight Dale P. Winebrenner Polar Science Center/Applied Physics Laboratory University of Washington

More information

SW06 Shallow Water Acoustics Experiment Data Analysis

SW06 Shallow Water Acoustics Experiment Data Analysis DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. SW06 Shallow Water Acoustics Experiment Data Analysis James F. Lynch MS #12, Woods Hole Oceanographic Institution Woods

More information

Direct Numerical Simulation of Aeolian Tones

Direct Numerical Simulation of Aeolian Tones THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN, KOREA, OCTOBER 27-30, 2003 Direct Numerical Simulation of Aeolian Tones Osamu Inoue 1 1. Institute of Fluid Science, Tohoku University,2-1-1 Katahira,

More information

Extension of the BLT Equation to Incorporate Electromagnetic Field Propagation

Extension of the BLT Equation to Incorporate Electromagnetic Field Propagation Extension of the BLT Equation to Incorporate Electromagnetic Field Propagation Fredrick M. Tesche Chalmers M. Butler Holcombe Department of Electrical and Computer Engineering 336 Fluor Daniel EIB Clemson

More information

Broadband matched-field source localization in the East China Sea*

Broadband matched-field source localization in the East China Sea* Broadband matched-field source localization in the East China Sea* Renhe Zhang Zhenglin Li Jin Yan Zhaohui Peng Fenghua Li National Laboratory of Acoustics, Institute of Acoustics, Chinese Academy of Sciences,

More information

HIGH-POWER SOLID-STATE LASER: LETHALITY TESTING AND MODELING

HIGH-POWER SOLID-STATE LASER: LETHALITY TESTING AND MODELING HIGH-POWER SOLID-STATE LASER: LETHALITY TESTING AND MODELING R. P. Abbott, C. D. Boley, S. N. Fochs, L. A. Nattrass, J. M. Parker, A. M. Rubenchik, J. A. Smith, and R. M. Yamamoto* University of California

More information

Parameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction

Parameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction Parameterizing the Effects of Upper-Ocean Large Eddies on Air-Sea Interaction Ming Li Horn Point Laboratory University of Maryland Center for Environmental Science 2020 Horn Point Road, Cambridge, MD 21613

More information

Crowd Behavior Modeling in COMBAT XXI

Crowd Behavior Modeling in COMBAT XXI Crowd Behavior Modeling in COMBAT XXI Imre Balogh MOVES Research Associate Professor ilbalogh@nps.edu July 2010 831-656-7582 http://movesinstitute.org Report Documentation Page Form Approved OMB No. 0704-0188

More information

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures AFOSR grant #FA9550-10-1-0061 Program manager: Dr. Les Lee PI: C.T. Sun Purdue University West Lafayette, Indiana

More information

CRS Report for Congress

CRS Report for Congress CRS Report for Congress Received through the CRS Web Order Code RS21396 Updated May 26, 2006 Summary Iraq: Map Sources Hannah Fischer Information Research Specialist Knowledge Services Group This report

More information

High Resolution Surface Characterization from Marine Radar Measurements

High Resolution Surface Characterization from Marine Radar Measurements DISTRIBUTION STATEMENT A: Distribution approved for public release; distribution is unlimited High Resolution Surface Characterization from Marine Radar Measurements Hans C. Graber CSTARS - University

More information

Modulation Instability of Spatially-Incoherent Light Beams and Pattern Formation in Incoherent Wave Systems

Modulation Instability of Spatially-Incoherent Light Beams and Pattern Formation in Incoherent Wave Systems Modulation Instability of Spatially-Incoherent Light Beams and Pattern Formation in Incoherent Wave Systems Detlef Kip, (1,2) Marin Soljacic, (1,3) Mordechai Segev, (1,4) Evgenia Eugenieva, (5) and Demetrios

More information

Design for Lifecycle Cost using Time-Dependent Reliability

Design for Lifecycle Cost using Time-Dependent Reliability Design for Lifecycle Cost using Time-Dependent Reliability Amandeep Singh Zissimos P. Mourelatos Jing Li Mechanical Engineering Department Oakland University Rochester, MI 48309, USA : Distribution Statement

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

uniform distribution theory

uniform distribution theory Uniform Distribution Theory 9 (2014), no.1, 21 25 uniform distribution theory ON THE FIRST DIGITS OF THE FIBONACCI NUMBERS AND THEIR EULER FUNCTION Florian Luca Pantelimon Stănică ABSTRACT. Here, we show

More information

THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS

THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS Annales Univ. Sci. Budapest., Sect. Comp. 40 (2013) nn nnn THE EULER FUNCTION OF FIBONACCI AND LUCAS NUMBERS AND FACTORIALS Florian Luca (Morelia, Mexico) Pantelimon Stănică (Monterey, USA) Dedicated to

More information

arxiv: v3 [cond-mat.stat-mech] 2 Dec 2015

arxiv: v3 [cond-mat.stat-mech] 2 Dec 2015 Inelastic Maxwell models for monodisperse gas-solid flows Aleksander Kubicki Departamento de Física, Universidad de Extremadura, Universidad de Extremadura, E-0607 Badajoz, Spain arxiv:4.345v3 [cond-mat.stat-mech]

More information

Analysis of Mixing and Dynamics Associated with the Dissolution of Hurricane-Induced Cold Wakes

Analysis of Mixing and Dynamics Associated with the Dissolution of Hurricane-Induced Cold Wakes DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. Analysis of Mixing and Dynamics Associated with the Dissolution of Hurricane-Induced Cold Wakes Carol Anne Clayson Dept.

More information

Babylonian resistor networks

Babylonian resistor networks IOP PUBLISHING Eur. J. Phys. 33 (2012) 531 537 EUROPEAN JOURNAL OF PHYSICS doi:10.1088/0143-0807/33/3/531 Babylonian resistor networks Carl E Mungan 1 and Trevor C Lipscombe 2 1 Physics Department, US

More information

Modeling the Impact of Extreme Events on Margin Sedimentation

Modeling the Impact of Extreme Events on Margin Sedimentation Modeling the Impact of Extreme Events on Margin Sedimentation Jasim Imran Department of Civil and Environmental Engineering, University of South Carolina, 3 Main Street, Columbia, SC 2928. phone: (83)

More information

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography

More information

Ocean Acoustics Turbulence Study

Ocean Acoustics Turbulence Study Ocean Acoustics Turbulence Study PI John Oeschger Coastal Systems Station 6703 West Highway 98 Panama City, FL 32407 phone: (850) 230-7054 fax: (850) 234-4886 email: OeschgerJW@ncsc.navy.mil CO-PI Louis

More information

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS) Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS) D.N. Kurk, T.E. Beegle, S.C. Spence and R.J. Swatski Report Documentation

More information

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,

More information

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10 Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The 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 information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles 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 information

Fleet Maintenance Simulation With Insufficient Data

Fleet Maintenance Simulation With Insufficient Data Fleet Maintenance Simulation With Insufficient Data Zissimos P. Mourelatos Mechanical Engineering Department Oakland University mourelat@oakland.edu Ground Robotics Reliability Center (GRRC) Seminar 17

More information

Weak Turbulence in Ocean Waves

Weak Turbulence in Ocean Waves Weak Turbulence in Ocean Waves Yuri V. Lvov Rensselaer Polytechnic Institute Troy, New York 12180-3590 Phone: (518) 276-6893 fax: (518) 276-4824 email: lvov@rpi.edu Award Number N000140210528 http://www.rpi.edu/~lvovy

More information

Impact of Typhoons on the Western Pacific Ocean DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds

Impact of Typhoons on the Western Pacific Ocean DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds Impact of Typhoons on the Western Pacific Ocean DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds Ramsey R. Harcourt Applied Physics Laboratory, University of Washington,

More information

Scattering of Internal Gravity Waves at Finite Topography

Scattering of Internal Gravity Waves at Finite Topography Scattering of Internal Gravity Waves at Finite Topography Peter Muller University of Hawaii Department of Oceanography 1000 Pope Road, MSB 429 Honolulu, HI 96822 phone: (808)956-8081 fax: (808)956-9164

More information

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models PI Isaac Ginis Graduate School of Oceanography

More information

Grant Number: N IP

Grant Number: N IP Coherence of Low-frequency Sound Signals Propagating through a Fluctuating Ocean: Analysis and Theoretical Interpretation of 2004 NPAL Experimental Data Alexander G. Voronovich NOAA/ESRL, PSD3, 325 Broadway,

More information

Real-Time Environmental Information Network and Analysis System (REINAS)

Real-Time Environmental Information Network and Analysis System (REINAS) Calhoun: The NPS Institutional Archive Faculty and Researcher Publications Faculty and Researcher Publications 1998-09 Real-Time Environmental Information Network and Analysis System (REINAS) Nuss, Wendell

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Range-Dependent Acoustic Propagation in Shallow Water with Elastic Bottom Effects

Range-Dependent Acoustic Propagation in Shallow Water with Elastic Bottom Effects DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Range-Dependent Acoustic Propagation in Shallow Water with Elastic Bottom Effects Robert I. Odom Applied Physics Laboratory

More information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012 arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,

More information

Wavelets and Affine Distributions A Time-Frequency Perspective

Wavelets and Affine Distributions A Time-Frequency Perspective Wavelets and Affine Distributions A Time-Frequency Perspective Franz Hlawatsch Institute of Communications and Radio-Frequency Engineering Vienna University of Technology INSTITUT FÜR NACHRICHTENTECHNIK

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 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 information

θ x = f ( x,t) could be written as

θ 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 information

Modelling and numerical methods for the diffusion of impurities in a gas

Modelling and numerical methods for the diffusion of impurities in a gas INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi

More information

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance Predictive Model for Archaeological Resources Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the

More information

Catalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids

Catalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids Catalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids William M. Nelson,, PhD IL Waste Management & Research Center 2003 Joint Service Scientific Conference on Chemical & Biological Defense Research

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization 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 information

Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG'

Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG' Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG' S.A. Burke 1, J.W. Wright 2, M.K. Robinson, B. Bronk 3, and R.L. Warren 1 1 Battelle Dugway Operations 2 U.S. Army Soldier

More information

KINETIC THEORY OF GRANULAR SOME APPLICATIONS

KINETIC THEORY OF GRANULAR SOME APPLICATIONS KINETIC THEORY OF GRANULAR BINARY MIXTURES: SOME APPLICATIONS Vicente Garzó Departamento t de Física, Universidad id dde Extremadura, Badajoz, Spain OUTLINE 1. Introduction 2. Granular binary mixtures.

More information

Grupo demecanica Celeste, Facultad de Ciencias, Valladolid, Spain. Dept. Applied Mathematics, University of Alicante, Alicante, Spain

Grupo demecanica Celeste, Facultad de Ciencias, Valladolid, Spain. Dept. Applied Mathematics, University of Alicante, Alicante, Spain Advances in the Unied Theory of the Rotation of the Nonrigid Earth Juan Getino Grupo demecanica Celeste, Facultad de Ciencias, Valladolid, Spain Jose M.Ferrandiz Dept. Applied Mathematics, University of

More information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A

More information

On Applying Point-Interval Logic to Criminal Forensics

On Applying Point-Interval Logic to Criminal Forensics On Applying Point-Interval Logic to Criminal Forensics (Student Paper) Mashhood Ishaque Abbas K. Zaidi Alexander H. Levis Command and Control Research and Technology Symposium 1 Report Documentation Page

More information

Chapter 2 Lagrangian Modeling

Chapter 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 information

Playing Abstract games with Hidden States (Spatial and Non-Spatial).

Playing Abstract games with Hidden States (Spatial and Non-Spatial). Playing Abstract games with Hidden States (Spatial and Non-Spatial). Gregory Calbert, Hing-Wah Kwok Peter Smet, Jason Scholz, Michael Webb VE Group, C2D, DSTO. Report Documentation Page Form Approved OMB

More information

STUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR

STUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR STUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR William R. Creasy EAI Corporation, 1308 Continental Drive, Suite J, Abingdon MD 21009 David J. McGarvey, Jeffrey S. Rice, Richard O'Connor,

More information

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION CETN-I-6 3185 DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION PURPOSE : To present a parameterized model of a directional spectrum of the sea surface using an energy spectrum and a value for the

More information

Chapter 6: Energy-Momentum Tensors

Chapter 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 information

Marginal Sea - Open Ocean Exchange

Marginal Sea - Open Ocean Exchange Marginal Sea - Open Ocean Exchange Michael A. Spall Mail Stop 21 Department of Physical Oceanography Woods Hole Oceanographic Institution Woods Hole, MA 02543-1541 phone: (508) 289-3342 fax: (508) 457-2181

More information

Contract No. N C0123

Contract No. N C0123 Particle Size Distribution and Optical Volume Scattering Function in the Mid and Upper Water Column of Optically Deep Coastal Regions: Transport from the Bottom Boundary Layer Y. C. Agrawal Sequoia Scientific,

More information

MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS

MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS Mingfu Li, Hsin-Min Peng, and Chia-Shu Liao Telecommunication Labs., Chunghwa

More information

Shallow Water Fluctuations and Communications

Shallow Water Fluctuations and Communications Shallow Water Fluctuations and Communications H.C. Song Marine Physical Laboratory Scripps Institution of oceanography La Jolla, CA 92093-0238 phone: (858) 534-0954 fax: (858) 534-7641 email: hcsong@mpl.ucsd.edu

More information

Spectral Flow, the Magnus Force, and the. Josephson-Anderson Relation

Spectral Flow, the Magnus Force, and the. Josephson-Anderson Relation Spectral Flow, the Magnus Force, an the arxiv:con-mat/9602094v1 16 Feb 1996 Josephson-Anerson Relation P. Ao Department of Theoretical Physics Umeå University, S-901 87, Umeå, SWEDEN October 18, 2018 Abstract

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

INFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS)

INFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS) AFRL-MN-EG-TP-2006-7403 INFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS) C. Michael Lindsay, National Research Council, Post Doctoral Research Associate Mario E. Fajardo

More information

Exact Solution of a Constrained. Optimization Problem in Thermoelectric Cooling

Exact Solution of a Constrained. Optimization Problem in Thermoelectric Cooling Applied Mathematical Sciences, Vol., 8, no. 4, 77-86 Exact Solution of a Constrained Optimization Problem in Thermoelectric Cooling Hongyun Wang Department of Applied Mathematics and Statistics University

More information