1
|
Clérouin J, Blanchet A, Blancard C, Faussurier G, Soubiran F, Bethkenhagen M. Equivalence between pressure- and structure-defined ionization in hot dense carbon. Phys Rev E 2022; 106:045204. [PMID: 36397512 DOI: 10.1103/physreve.106.045204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Accepted: 09/23/2022] [Indexed: 06/16/2023]
Abstract
The determination of the ionization of a system in the hot dense regime is a long standing issue. Recent studies have shown inconsistencies between standard predictions using average atom models and evaluations deduced from electronic transport properties computed with quantum molecular dynamics simulations [Bethkenhagen et al., Phys. Rev. Res. 2, 023260 (2020)]2643-156410.1103/PhysRevResearch.2.023260. Here, we propose a definition of the ionization based on its effect on the plasma structure as given by the pair distribution function (PDF), and on the concept of effective one-component plasma (eOCP). We also introduce a definition based on the total pressure and on a modelization of the electronic pressure. We show the equivalence of these definitions on two studies of carbon along the 100 eV isotherm and the 10 g/cm^{3} isochor. Simulations along the 100 eV isotherm are obtained with the newly implemented Ext. First principles molecular dynamics (Fpmd) method in Abinit for densities ranging from 1 to 500 g/cm^{3}and along the 10 g/cm^{3} isochor with the recently published Spectral quadrature DFT (Sqdft) simulations, between 8 and 860 eV. The resulting ionizations are compared to the predictions of the average-atom code Qaam which is based on the muffin-tin approximation. A disagreement between the eOCP and the actual PDFs (non-OCP behavior) is interpreted as the onset of bonding in the system.
Collapse
Affiliation(s)
- Jean Clérouin
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Augustin Blanchet
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Christophe Blancard
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Gérald Faussurier
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - François Soubiran
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Mandy Bethkenhagen
- CNRS, École Normale Supérieure de Lyon, Laboratoire de Géologie de Lyon LGLTPE UMR 5276, Centre Blaise Pascal, 46 allée d'Italie Lyon 69364, France
| |
Collapse
|
2
|
Wu J, González-Cataldo F, Soubiran F, Militzer B. The phase diagrams of beryllium and magnesium oxide at megabar pressures. J Phys Condens Matter 2022; 34:144003. [PMID: 35026747 DOI: 10.1088/1361-648x/ac4b2a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Accepted: 01/13/2022] [Indexed: 06/14/2023]
Abstract
We performab initiosimulations of beryllium (Be) and magnesium oxide (MgO) at megabar pressures and compare their structural and thermodynamic properties. We make a detailed comparison of our two recently derived phase diagrams of Be (Wuet al2021Phys. Rev.B104014103) and MgO (Soubiran and Militzer 2020Phys. Rev. Lett.125175701) using the thermodynamic integration technique, as they exhibit striking similarities regarding their shape. We explore whether the Lindemann criterion can explain the melting temperatures of these materials through the calculation of the Debye temperature at high pressure. From our free energy calculations, we find that the melting line of both materials is well represented by the Simon-Glazel fitTm(P) =T0(1 +P/a)1/c, whereT0= 1564 K,a= 15.8037 GPa andc= 2.4154 for Be, whileT0= 3010 K,a= 10.5797 GPa andc= 2.8683 for the MgO in the B1. For the B2 phase, we use the valuesa= 26.1163 GPa andc= 2.2426. Both materials exhibit negative Clapeyron slopes on the boundaries between the two solid phases that are strongly affected by anharmonic effects, which also influence the location of the solid-solid-liquid triple point. We find that the quasi-harmonic approximation underestimates the stability range of the low-pressure phases, namely hcp for Be and B1 for MgO. We also compute the phonon dispersion relations at low and high pressure for each of the phases of these materials, and also explore how the phonon density of states is modified by temperature. Finally, we derive secondary shock Hugoniot curves in addition to the principal Hugoniot curve for both materials, and study their offsets in pressure between solid and liquid branches.
Collapse
Affiliation(s)
- Jizhou Wu
- Department of Earth and Planetary Science, University of California, Berkeley, CA 94720, United States of America
| | - Felipe González-Cataldo
- Department of Earth and Planetary Science, University of California, Berkeley, CA 94720, United States of America
| | | | - Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, CA 94720, United States of America
- Department of Astronomy, University of California, Berkeley, CA 94720, United States of America
| |
Collapse
|
3
|
Guarguaglini M, Soubiran F, Hernandez JA, Benuzzi-Mounaix A, Bolis R, Brambrink E, Vinci T, Ravasio A. Electrical conductivity of warm dense silica from double-shock experiments. Nat Commun 2021; 12:840. [PMID: 33547308 PMCID: PMC7865001 DOI: 10.1038/s41467-021-21046-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Accepted: 01/04/2021] [Indexed: 11/09/2022] Open
Abstract
Understanding materials behaviour under extreme thermodynamic conditions is fundamental in many branches of science, including High-Energy-Density physics, fusion research, material and planetary science. Silica (SiO2) is of primary importance as a key component of rocky planets' mantles. Dynamic compression is the most promising approach to explore molten silicates under extreme conditions. Although most experimental studies are restricted to the Hugoniot curve, a wider range of conditions must be reached to distill temperature and pressure effects. Here we present direct measurements of equation of state and two-colour reflectivity of double-shocked α-quartz on a large ensemble of thermodynamic conditions, which were until now unexplored. Combining experimental reflectivity data with numerical simulations we determine the electrical conductivity. The latter is almost constant with pressure while highly dependent on temperature, which is consistent with simulations results. Based on our findings, we conclude that dynamo processes are likely in Super-Earths' mantles.
Collapse
Affiliation(s)
- M Guarguaglini
- LULI, CNRS, CEA, École Polytechnique - Institut Polytechnique de Paris, Palaiseau cedex, France.,Sorbonne Université, Faculté des Sciences et Ingénierie, Laboratoire d'utilisation des lasers intenses (LULI), CNRS, Paris, France
| | - F Soubiran
- École Normale Supérieure de Lyon, Université Lyon 1, Laboratoire de Géologie de Lyon, Lyon, France.,CEA DAM-DIF, Arpajon, France
| | - J-A Hernandez
- LULI, CNRS, CEA, École Polytechnique - Institut Polytechnique de Paris, Palaiseau cedex, France.,Sorbonne Université, Faculté des Sciences et Ingénierie, Laboratoire d'utilisation des lasers intenses (LULI), CNRS, Paris, France.,Centre for Earth Evolution and Dynamics, University of Oslo, Oslo, Norway
| | - A Benuzzi-Mounaix
- LULI, CNRS, CEA, École Polytechnique - Institut Polytechnique de Paris, Palaiseau cedex, France.,Sorbonne Université, Faculté des Sciences et Ingénierie, Laboratoire d'utilisation des lasers intenses (LULI), CNRS, Paris, France
| | - R Bolis
- LULI, CNRS, CEA, École Polytechnique - Institut Polytechnique de Paris, Palaiseau cedex, France.,Sorbonne Université, Faculté des Sciences et Ingénierie, Laboratoire d'utilisation des lasers intenses (LULI), CNRS, Paris, France
| | - E Brambrink
- LULI, CNRS, CEA, École Polytechnique - Institut Polytechnique de Paris, Palaiseau cedex, France.,Sorbonne Université, Faculté des Sciences et Ingénierie, Laboratoire d'utilisation des lasers intenses (LULI), CNRS, Paris, France
| | - T Vinci
- LULI, CNRS, CEA, École Polytechnique - Institut Polytechnique de Paris, Palaiseau cedex, France.,Sorbonne Université, Faculté des Sciences et Ingénierie, Laboratoire d'utilisation des lasers intenses (LULI), CNRS, Paris, France
| | - A Ravasio
- LULI, CNRS, CEA, École Polytechnique - Institut Polytechnique de Paris, Palaiseau cedex, France. .,Sorbonne Université, Faculté des Sciences et Ingénierie, Laboratoire d'utilisation des lasers intenses (LULI), CNRS, Paris, France.
| |
Collapse
|
4
|
Li Z, Winisdoerffer C, Soubiran F, Caracas R. Ab initio Gibbs ensemble Monte Carlo simulations of the liquid-vapor equilibrium and the critical point of sodium. Phys Chem Chem Phys 2021; 23:311-319. [PMID: 33347522 DOI: 10.1039/d0cp04158k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The ab initio (ai) Gibbs ensemble (GE) Monte Carlo (MC) method coupled with Kohn-Sham density functional theory is successful in predicting the liquid-vapour equilibrium of insulating systems. Here we show that the aiGEMC method can be used to study also metallic systems, where the excited electronic states play an important role and cannot be neglected. For this we include the electronic free energy in the formulation of the effective energy of the system to be used in the acceptance criteria for the MC moves. The application of this aiGEMC method to sodium yields a good agreement with available experimental data on the liquid-vapour equilibrium densities. We predict a critical point for sodium at 2338 ± 108 K and 0.24 ± 0.03 g cm-3. The liquid structure stemming from aiGEMC simulations is very similar to the one from ab initio molecular dynamics. Since this method can determine phase transition without computing the Gibbs free energy, it may offer a new possibility to study other materials with a reasonable computational cost.
Collapse
Affiliation(s)
- Zhi Li
- CNRS, École Normale Supérieure de Lyon, Laboratoire de Géologie de Lyon LGLTPE UMR 5276, 46 allée d'Italie, Lyon 69364, France.
| | | | | | | |
Collapse
|
5
|
Militzer B, González-Cataldo F, Zhang S, Driver KP, Soubiran F. First-principles equation of state database for warm dense matter computation. Phys Rev E 2021; 103:013203. [PMID: 33601631 DOI: 10.1103/physreve.103.013203] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Accepted: 12/11/2020] [Indexed: 06/12/2023]
Abstract
We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF, B_{4}C, BN, CH_{4}, CH_{2}, C_{2}H_{3}, CH, C_{2}H, MgO, and MgSiO_{3}. For all these materials, we provide the pressure and internal energy over a density-temperature range from ∼0.5 to 50 g cm^{-3} and from ∼10^{4} to 10^{9} K, which are based on ∼5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L- and K-shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. We derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H_{2}O, H_{2}O_{2}, Al_{2}O_{3}, CO, and CO_{2} to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.
Collapse
Affiliation(s)
- Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Department of Astronomy, University of California, Berkeley, California 94720, USA
| | - Felipe González-Cataldo
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Shuai Zhang
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| | - Kevin P Driver
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - François Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- CEA DAM-DIF, 91297 Arpajon, France
| |
Collapse
|
6
|
Soubiran F, Militzer B. Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime. Phys Rev Lett 2020; 125:175701. [PMID: 33156661 DOI: 10.1103/physrevlett.125.175701] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Revised: 06/15/2020] [Accepted: 09/17/2020] [Indexed: 06/11/2023]
Abstract
With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20 000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T=10000 K and P=370 GPa, which is higher in pressure than was inferred with quasiharmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. This may render it difficult to observe this phase with shock experiments because of kinetic effects.
Collapse
Affiliation(s)
- François Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- École Normale Supérieure de Lyon, Université Lyon 1, Laboratoire de Géologie de Lyon, CNRS UMR 5276, 69364 Lyon Cedex 07, France
- CEA DAM-DIF, 91297 Arpajon, France
| | - Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Department of Astronomy, University of California, Berkeley, California 94720, USA
| |
Collapse
|
7
|
Soubiran F, González-Cataldo F, Driver KP, Zhang S, Militzer B. Magnesium oxide at extreme temperatures and pressures studied with first-principles simulations. J Chem Phys 2019; 151:214104. [DOI: 10.1063/1.5126624] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- François Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- École Normale Supérieure de Lyon, Université Lyon 1, Laboratoire de Géologie de Lyon, CNRS UMR 5276, 69364 Lyon Cedex 07, France
| | - Felipe González-Cataldo
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Kevin P. Driver
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Shuai Zhang
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14623, USA
| | - Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Department of Astronomy, University of California, Berkeley, California 94720, USA
| |
Collapse
|
8
|
Abstract
We perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of 0.1-32.4gcm^{-3} and 10^{4}-10^{8} K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multishock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multishock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and 〈Z〉 to reveal the evolution of the plasma structure and ionization behavior.
Collapse
Affiliation(s)
- K P Driver
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - F Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - B Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Department of Astronomy, University of California, Berkeley, California 94720, USA
| |
Collapse
|
9
|
Zhang S, Militzer B, Benedict LX, Soubiran F, Sterne PA, Driver KP. Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas. J Chem Phys 2018; 148:102318. [DOI: 10.1063/1.5001208] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shuai Zhang
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Department of Astronomy, University of California, Berkeley, California 94720, USA
| | - Lorin X. Benedict
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - François Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Philip A. Sterne
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Kevin P. Driver
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
10
|
Zhang S, Driver KP, Soubiran F, Militzer B. First-principles equation of state and shock compression predictions of warm dense hydrocarbons. Phys Rev E 2017; 96:013204. [PMID: 29347225 DOI: 10.1103/physreve.96.013204] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2017] [Indexed: 06/07/2023]
Abstract
We use path integral Monte Carlo and density functional molecular dynamics to construct a coherent set of equations of state (EOS) for a series of hydrocarbon materials with various C:H ratios (2:1, 1:1, 2:3, 1:2, and 1:4) over the range of 0.07-22.4gcm^{-3} and 6.7×10^{3}-1.29×10^{8}K. The shock Hugoniot curve derived for each material displays a single compression maximum corresponding to K-shell ionization. For C:H = 1:1, the compression maximum occurs at 4.7-fold of the initial density and we show radiation effects significantly increase the shock compression ratio above 2 Gbar, surpassing relativistic effects. The single-peaked structure of the Hugoniot curves contrasts with previous work on higher-Z plasmas, which exhibit a two-peak structure corresponding to both K- and L-shell ionization. Analysis of the electronic density of states reveals that the change in Hugoniot structure is due to merging of the L-shell eigenstates in carbon, while they remain distinct for higher-Z elements. Finally, we show that the isobaric-isothermal linear mixing rule for carbon and hydrogen EOS is a reasonable approximation with errors better than 1% for stellar-core conditions.
Collapse
Affiliation(s)
- Shuai Zhang
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Kevin P Driver
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - François Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA and Department of Astronomy, University of California, Berkeley, California 94720, USA
| |
Collapse
|
11
|
Zhang S, Driver KP, Soubiran F, Militzer B. Equation of state and shock compression of warm dense sodium—A first-principles study. J Chem Phys 2017; 146:074505. [DOI: 10.1063/1.4976559] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shuai Zhang
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Kevin P. Driver
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - François Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Department of Astronomy, University of California, Berkeley, California 94720, USA
| |
Collapse
|
12
|
Driver KP, Soubiran F, Zhang S, Militzer B. First-principles equation of state and electronic properties of warm dense oxygen. J Chem Phys 2015; 143:164507. [DOI: 10.1063/1.4934348] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- K. P. Driver
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - F. Soubiran
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - Shuai Zhang
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
| | - B. Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA
- Department of Astronomy, University of California, Berkeley, California 94720, USA
| |
Collapse
|
13
|
|
14
|
Vidal V, Soubiran F, Divoux T, Géminard JC. Degassing cascades in a shear-thinning viscoelastic fluid. Phys Rev E Stat Nonlin Soft Matter Phys 2011; 84:066302. [PMID: 22304185 DOI: 10.1103/physreve.84.066302] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/28/2011] [Indexed: 05/31/2023]
Abstract
We report the experimental study of the degassing dynamics through a thin layer of shear-thinning viscoelastic fluid when a constant air flow is imposed at its bottom. The fluid is an aqueous solution of cetyltrimethylammonium bromide (CTAB) and sodium salicylate (NaSal). Over a large range of parameters, the air is periodically released through a series of successive bubbles, hereafter named cascades. Each cascade is followed by a continuous degassing, lasting for several seconds, corresponding to an open channel crossing the fluid layer. The periodicity between two cascades does not depend on the injected flow rate. Inside one cascade, the properties of the overpressure signal associated with the successive bubbles vary continuously. The pressure threshold above which the fluid starts flowing, fluid deformation and pressure drop due to degassing through the thin fluid layer can be simply described by a Maxwell model.
Collapse
Affiliation(s)
- Valérie Vidal
- Université de Lyon, Laboratoire de Physique, École Normale Supérieure de Lyon, CNRS, 46 Allée d'Italie, FR-69364 Lyon cedex 07, France
| | | | | | | |
Collapse
|