Sonoluminescence Test for Equation of State in Warm Dense Matter

HIFAN 1673 SONOLUMINESCENCE TEST FOR EQUATION OF STATE IN WARM DENSE MATTER Siu-Fai Ng1, 2, J. J. Barnard3, P. T. Leung1, S. S. Yu1, 2 1. Physics Department, The Chinese University of Hong Kong, Hong Kong SAR, China 2. Lawrence Berkeley National Laboratory, California, USA 3. Lawrence Livermore National Laboratory, California, USA Accelerator Fusion Research Division Ernest Orlando Lawrence Berkeley National Laboratory University of California Berkeley, California 94720 August 2008 This work was supported by the Director, Office of Science, Office of Fusion Energy Sciences, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. SONOLUMINESCENCE TEST FOR EQUATION OF STATE IN WARM DENSE MATTER Siu-Fai Ng 1, 2, J. J. Barnard 3, P. T. Leung 1, S. S. Yu 1, 2 1. Physics Department, The Chinese University of Hong Kong, Hong Kong SAR, China 2. Lawrence Berkeley National Laboratory, California, USA 3. Lawrence Livermore National Laboratory, California, USA In experiments of Single-bubble Sonoluminescence (SBSL), the bubble is heated to temperatures of a few eV in the collapse phase of the oscillation. Our hydrodynamic simulations show that the density inside the bubble can go up to the order of 1 g/cm 3, and the electron density due to ionization is 10 21 /cm 3. So the plasma coupling constant is found to be around 1 and the gas inside the bubble is in the Warm Dense Matter (WDM) regime. We simulate the light emission of SL with an optical model for thermal radiation which takes the finite opacity of the bubble into consideration. The numerical results obtained are compared to the experimental data and found to be very sensitive to the equation of state used. As theories for the equation of state, as well as the opacity data, in the WDM regime are still very uncertain, we propose that SL may be a good low-cost experimental check for the EOS and the opacity data for matter in the WDM regime. 1. Introduction good and easy check for the EOS and opacity in WDM. Single-bubble Sonoluminescence (SBSL or SL), which was 2. Sonoluminescence Simulations discovered in 1989, is a phenomenon of periodic light In our SL simulations, we used the 1-D code DISH 8 to emission by an oscillating gas bubble trapped in the model hydrodynamics inside the bubble. The DISH code, pressure anti-node of a standing ultrasound wave in water originally written for the slab geometry, was modified for (or other fluids). Several experiments have shown that the the (assumed) perfect spherical symmetry of SL. width of the emitted light pulse is of the order of 100 ps For the outer boundary condition for the bubble, we with peak power of the order of 10 mW. 1-3 approximate the dynamics of the water outside the bubble Our simulations, as well as the results in many papers 1-7, by the commonly used Rayleigh-Plesset (RP) equation9: show that the temperature inside the bubble can go up to a 2Rɺ 3 4 R ɺ 3 ɺɺ ɺ 2 few eV and the density is around 1 g/cm in the collapse ρ 1− RR +− 1 R =− Pgas P apply c 2 3 c phase. At this temperature and density, it is believed that the (1) R Rdɺ R ɺ σ bubble has entered into the Warm Dense Matter (WDM) −+−P0 1 P gas − 4η − 2 c c dt R R regime, a regime of significant theoretical interest because of the strongly coupled effect between ions and atom. (In where ρ , c ,η and σ are density, speed of sound, viscosity general, WDM is defined to be the state of material at and surface tension of the water, R is the radius of the temperature of the order of 0.1 to 20 eV and density of 0.1 bubble as a function of time, P0 is the constant ambient to 10 solid density, with plasma coupling constant of around pressure ( P0 =1 atm ), Papply is the applied pressure 1.) Since the equation of state, opacity and other properties in the WDM regime are still very uncertain, ( Papply= − P a sin ω t ), and Pgas is the gas pressure just many experiments (such as high intensity heavy ion beam target experiments) have been proposed to probe these inside the bubble wall, which depends on the hydrodynamic properties empirically. evolution of the bubble and the equation of state used. Since the light pulse emitted in the collapse phase of SL can Our simulation includes thermal conduction with Argon be measured quite precisely, we propose that SL may be a thermal conductivity κ =0.009 + 0.00032 T 10 , where κ is in W/m K and T is in K. WDM and other high-pressure phenomena. In QEOS, the To calculate the radiation emitted (power emitted per electronic properties are obtained from a modified wavelength interval), we use the integral form of the Thomas-Fermi statistical model, while the ion thermal radiation transport equation 11 : motion is described by a multiphase equation of state R combining Debye, Gruneisen, Lindemann, and fluid-scaling Pt()= 16π2 κ tot (,) rtIrt Pl (,) λ∫0 λ λ laws. QEOS can give smooth predictions for ionization R (2) ×exp − κ tot (r′ ) dr ′ r2 dr state, pressure, energy, entropy and Helmholtz free energy ( ∫r λ ) for use in hydrodynamic simulations. tot Pl B. Saha-based EOS where κλ is the total absorption coefficient and Iλ is the We also use the well-known Saha equation to get the spectral emission by Planck blackbody radiation: ionization of the gas 14,15 : 2 3/ 2 ion Pl 2hc nnQe ζ2 ζ2π me kT EE ζ − ∆ Iλ = (3) =exp − (6) λ5 exphc / λ kT − 1 2 () nζ−1 Q ζ − 1 h kT We consider the free-free interactions between electrons This equation gives the ratio of the ion density of charge and ions and between electrons and atoms as the two state ζ and ζ −1 in terms of electron density ne , 5,11,12 dominant processes for radiation absorption . So ion temperature T and ionization energy Eζ . ∆E is the tot ff+ ff0 κλ= κ λ + κ λ and: continuum lowering term to reduce the energy needed for 1/ 2 ionization due to the interaction within ions and between 2 6 3 ff+ 4 2 π Z e nni e λ κλ = (4) 3 3 m kT 3 4 ions and electrons. Q is the partition function of the e ()4πε 0 hc m e ζ 3/ 2 individual ion with charge +ζ . Mathematically, the infinite 2 2 ff0 4e ne n 0λ 2kT d κλ =3 c + (5) 3/2 3 m3 kT series Q diverges for an isolated ion. To converge the ()4πε0 π c e ζ where n0 , ni and ne are the density of atoms, ions and series, we truncate the series when the excitation energy electrons respectively, and are dependent on the EOS ex ion Eζm ≥ E ζ +1 − ∆ E . assumed. Z is the charge of the ions. c and d are constants equal to 1.6× 10−20 m 2 /eV and We consider the ion sphere model to calculate ∆E 14 . In an −0.6 × 10−20 m 2 . ion sphere with an infinitely small +ζ ion at origin and 3. Equation of State ζ free electrons distributed uniformly in the sphere, the Results of sonoluminescence simulations depend on the net charge in the sphere is zero. Assume that the charges equation of state (EOS) used but the EOS in WDM regime outside the sphere neutralize each other and have no effect is still uncertain. Here we use two different equations of on the potential at the origin, we can calculate the extra state that are commonly used in WDM simulations: QEOS potential energy on that ion due to the free electrons in the and Saha-based EOS. sphere: A. QEOS 1/3 13 3 4 π QEOS (The quotidian equation of state) is a well-known 2 1/3 ∆E = ζ e n i (7) 2 3 equation of state model for hydrodynamic simulations of We also include the hard core effect due to interactions T0 =8 ° C and driving pressure Pa = 1.3 bar . As stable between gas particles in the Saha-based EOS model 16 . SL bubble is believed to be a noble gas bubble 1,18 , the SL Finally, we obtain the pressure and energy per unit mass bubble formed in this case is basically an Argon bubble (air (EOS): is 1 % Argon). We find by diffusive theories 19 (using gas 1 2 concentration to be the Argon concentration) that the P=+() n nkTe − n e ∆−− E () abkTn (8) 3 ambient radius r0 = 4µ m . Simulations using the two EOS described above were then performed. 1 3 e( n n ) kT an2 n E ex = +e −+ ∑ ζ ζ Simulations using Saha-based EOS show that at the instant ρ 2 ζ =0 (9) of minimum radius, the temperature at the center will reach ζ ion +∑nζ ∑ Em −ζ ∆ E its maximum of nearly 25000 K (Fig. 1A). At this high ζ =1m = 1 temperature, the Argon gas will ionize. Fig. 1B shows that ex the maximum ionization is about 22 %. Fig. 2 shows the where n is the total density of atoms and ions, Eζ is density, temperature, pressure and charge state profile at the the average excitation energy for the +ζ ions, constants stagnation point. It is obvious that the bubble is not uniform 3 −10 but has a cold and dense shell and a hot core.

Load more