Search

We have developed a gas flow model in the spent nuclear fuel during the annealing. It postulates that the gas release during an isothermal plateau at 1200{\textdegree}C corresponds to the equilibrium between overpressure gas reservoirs in the fuel sample connected to the free surface at atmospheric pressure.

We review recent advances in the design, synthesis, and modeling of active fluids. Active fluids have been at the center of many technological innovations and theoretical advances over the past two decades. Research on this new class of fluids has been inspired by the fascinating and remarkably efficient strategies that biological systems employ, leading to the development of biomimetic nano- and micro-machines and -swimmers. The review encompasses active fluids on both the nano- and micro-scale. We start with examples of biological active systems before we discuss how experimentalists leverage novel propulsion mechanisms to power nano- and micro-machines. We then examine how the study of these far-from-equilibrium systems has prompted the development of new simulation methods and theoretical models in nonquilibrium physics to account for their mechanical, thermodynamic and emergent properties. Recent advances in the field have paved the way for the design, synthesis, and modeling of autonomous systems at the nano- and micro-scale and open the door to the development of soft matter robotics.

The entropy change that occurs upon mixing two fluids has remained an intriguing topic since the dawn of statistical mechanics. In this work, we generalize the grand-isobaric ensemble to mixtures, and develop a Monte Carlo algorithm for the rapid determination of entropy in these systems. A key advantage of adiabatic ensembles is the direct connection they provide with entropy. Here, we show how the entropy of a binary mixture A-B can be readily obtained in the adiabatic grand-isobaric $(\mu_{\text{A}}$, $\mu_{\text{B}}, P, R)$ ensemble, in which $\mu_{\text{A}}$ and $\mu_{\text{B}}$ denote the chemical potential of components A and B, respectively, $P$ is the pressure, and $R$ is the heat (Ray) function, that corresponds to the total energy of the system. This, in turn, allows for the evaluation of the entropy of mixing, as well as of the Gibbs free energy of mixing. We also demonstrate that our approach performs very well both on systems modeled with simple potentials and with complex many-body force fields. Finally, this approach provides a direct route to the determination of the thermodynamic properties of mixing, and allows for the efficient detection of departures from ideal behavior in mixtures.

We extend Expanded Wang-Landau (EWL) simulations beyond classical systems and develop the EWL method for systems modeled with a tight-binding Hamiltonian. We then apply the method to determine the partition function and thus all thermodynamic properties, including the Gibbs free energy and entropy, of the fluid phases of Si. We compare the results from quantum many-body (QMB) tight binding models, which explicitly calculate the overlap between the atomic orbitals of neighboring atoms, to those obtained with classical many-body force fields (CMB), which allow to recover the tetrahedral organization in condensed phases of Si through e.g. a repulsive 3-body term that favors the ideal tetrahedral angle. Along the vapor-liquid coexistence, between 3000K and 6000K, the densities for the two coexisting phases are found to vary significantly (by $5$ orders of magnitude for the vapor and by up to 25% for the liquid) and to provide a stringent test of the models. Transitions from vapor to liquid are predicted to occur for chemical potentials that are $10-15$% higher for CMB models than for QMB models, and a ranking of the force fields is provided by comparing the predictions for the vapor pressure to the experimental data. QMB models also reveal the formation of a gap in the electronic density of states of the coexisting liquid at high temperatures. Subjecting Si to a nanoscopic confinement has a dramatic effect on the phase diagram, with e.g. at 6000K a decrease in liquid densities by about 50% for both CMB and QMB models and an increase in vapor densities between 90% (CMB) and 170% (QMB). The results presented here provide a full picture of the impact of the strategy (CMB or QMB) chosen to model many-body effects on the thermodynamic properties of the fluid phases of Si.

The determination of the critical parameters of metals has remained particularly challenging both experimentally, because of the very large temperatures involved, and theoretically, because of the many-body interactions that take place in metals. Moreover, experiments have shown that these systems exhibit an unusually strong asymmetry of their binodal. Recent theoretical work has led to new similarity laws, based on the calculation of the Zeno line and of the underlying Boyle parameters, which provided results for the critical properties of atomic and molecular systems in excellent agreement with experiments. Using the recently developed Expanded Wang-Landau (EWL) simulation method, we evaluate the grand-canonical partition function, over a wide range of conditions, for $11$ $FCC$ and $HCP$ metals ($Ag$, $Al$, $Au$, $Be$, $Cu$, $Ir$, $Ni$, $Pb$, $Pd$, $Pt$ and $Rh$), modeled with a many-body interaction potential. This allows us to calculate the binodal, Zeno line, Boyle parameters and, in turn, obtain the critical properties for these systems. We also propose two scaling laws for the enthalpy and entropy of vaporization, and identify critical exponents of $0.4$ and $1.22$ for these two laws, respectively.

Using molecular simulation, we analyze the capillary condensation and evaporation processes for Argon confined in a cylindrical nanopore. For this purpose, we define the entropy of the adsorbed fluid as a reaction coordinate and determine the free energy associated with both processes along entropic pathways. For capillary condensation, we identify a complex free energy profile resulting from the multi-stage nature of this phenomenon. We find capillary condensation to proceed through the nucleation of a liquid bridge across the nanopore, followed by its expansion throughout the pore to give rise to the stable phase of high density. In the case of capillary evaporation, the free energy profile along the entropy pathway also exhibits different regimes, corresponding to the initial destabilization of the layered structure of the fluid followed by the formation, and subsequent expansion, of a bubble across the nanopore.

Using molecular simulation, we study the nucleation of liquid droplets from binary mixtures and determine the free energy of nucleation along entropic pathways. To this aim, we develop the $\mu_1 \mu_2 VT-S$ method, based on the grand-canonical ensemble modeling the binary mixture, and use the entropy of the system $S$ as the reaction coordinate to drive the formation of the liquid droplet. This approach builds on the advantages of the grand-canonical ensemble, which allows for the direct calculation of the entropy of the system and lets the composition of the system free to vary throughout the nucleation process. Starting from a metastable supersaturated vapor, we are able to form a liquid droplet by gradually decreasing the value of $S$, through a series of umbrella sampling simulations, until a liquid droplet of a critical size has formed. The $\mu_1 \mu_2 VT-S$ method also allows us to calculate the free energy barrier associated with the nucleation process, to shed light on the relation between supersaturation and free energy of nucleation, and to analyze the interplay between the size of the droplet and its composition during the nucleation process.

Using the entropy $S$ as a reaction coordinate, we determine the free energy barrier associated with the formation of a liquid droplet from a supersaturated vapor for atomic and molecular fluids. For this purpose, we develop the $\mu VT-S$ simulation method that combines the advantages of the grand-canonical ensemble, that allows for a direct evaluation of the entropy, and of the umbrella sampling method, that is well suited to the study of an activated process like nucleation. Applying this approach to an atomic system such as $Ar$ allows us to test the method. The results show that the $\mu VT-S$ method gives the correct dependence on supersaturation of the height of the free energy barrier and of the size of the critical droplet, when compared to predictions from classical nucleation theory and to previous simulation results. In addition, it provides insight into the relation between entropy and droplet formation throughout this process. An additional advantage of the $\mu VT-S$ approach is its direct transferability to molecular systems, since it uses the entropy of the system as the reaction coordinate. Applications of the $\mu VT-S$ simulation method to $N_2$ and $CO_2$ are presented and discussed in this work, showing the versatility of the $\mu VT-S$ approach.

Using molecular simulations, we determine the locus of ideal compressibility, or Zeno line, for a series of ionic compounds. We find that the shape of this thermodynamic contour follows a linear law, leading to the determination of the Boyle parameters. We also show that a similarity law, based on the Boyle parameters, yields accurate critical data when compared to the experiment. Furthermore, we show that the Boyle density scales linearly with the size-asymmetry, providing a direct route to establish a correspondence between the thermodynamic properties of different ionic compounds.

Using molecular simulation, we study the nucleation process from supercooled liquid alloys of Cu and Pd. The simulations reveal a complex interplay between the size of the crystal nucleus and its composition that greatly impacts the crystallization process on the nanoscale. In particular, we find that the free energy of nucleation strongly depends on the composition of the alloy, with a free energy barrier for the equimolar alloy that is more than two times larger than for the pure metals. We attribute this increase in free energy to the variations in composition occurring both at the surface and in the core of the nucleus. The local changes in composition are then analyzed by comparing the surface energies of the two metals and by taking into account the competition between crystallization and demixing that takes place at the interface between the nucleus and the surrounding liquid.

Using molecular simulation, we determine Ginzburg-Landau free energy functions for molecular fluids. To this aim, we extend the Expanded Wang-Landau method to calculate the partition functions, number distributions and Landau free energies for $Ar$, $CO_2$ and $H_2O$. We then parametrize a coarse-grained free energy function of the density order parameter and assess the performance of this free energy function on its ability to model the onset of criticality in these systems. The resulting parameters can be readily used in hybrid atomistic/continuum simulations that connect the microscopic and mesoscopic length scales.

We use molecular simulation to compute the thermodynamic properties of 7 rigid models for water (SPC/E, TIP3P, TIP4P, TIP4P/2005, TIP4P/Ew, TIP5P, OPC) over a wide range of temperature and pressure. Carrying out Expanded Wang-Landau simulations, we obtain a high accuracy estimate for the grand-canonical partition function which, in turn, provides access to all properties, including the free energy and entropy, both at the vapor-liquid coexistence and for compressed water. The results at coexistence highlight the close connection between the behavior of the statistical functions and the thermodynamic properties. They show that the subgroup (SPC/E,TIP4P/2005,TIP4P/Ew) provides the most accurate account of the vapor-liquid equilibrium properties. For compressed water, the comparison of the simulation results to the experimental data establishes that the TIP4P/Ew model performs best among the 7 models considered here, and captures the experimental trends for the dependence of entropy and molar Gibbs free energy on pressure.

Using molecular simulation, we assess the impact of an electric field on the properties of water, modeled with the SPC/E potential, over a wide range of states and conditions. Electric fields of the order of $0.1V/$\AA and beyond are found to have a significant impact on the grand-canonical partition function of water, resulting in shifts in the chemical potential at the vapor-liquid coexistence of up to $20$%. This, in turn, leads to increases in the critical temperatures by close to $7$% for a field of $0.2V/$\AA, to lower vapor pressures, and to much larger entropies of vaporization (by up to $35$%). We interpret these results in terms of the greater density change at the transition and of the increased structural order resulting from the applied field. The thermodynamics of compressed liquids and of supercritical water are also analyzed over a wide range of pressures, leading to the determination of the Zeno line and of the curve of ideal enthalpy that span the supercritical region of the phase diagram. Rescaling the phase diagrams obtained for the different field strength by their respective critical properties allows us to draw a correspondence between these systems for fields of up to $0.2V/$\AA.

We determine the impact of the Friedel oscillations on the phase behavior, critical properties and thermodynamic contours in films ($2D$) and bulk phases ($3D$). Using Expanded Wang-Landau simulations, we calculate the grand-canonical partition function and, in turn, the thermodynamic properties of systems modeled with a linear combination of the Lennard-Jones and Dzugutov potentials, weighted by a parameter $X$ ($0

Crystallization often proceeds through successive stages that lead to a gradual increase in organization. Using molecular simulation, we determine the nucleation pathway for solid solutions of copper and gold. We identify a new nucleation mechanism (liquid$\to$$L1_2$~precursor$\to$solid solution), involving a chemically ordered intermediate that is more organized than the end product. This nucleation pathway arises from the low formation energy of $L1_2$ clusters which, in turn, promote crystal nucleation. We also show that this mechanism is composition-dependent since the high formation energy of other ordered phases precludes them from acting as precursors.

Using molecular simulations, we shed light on the mechanism underlying crystal nucleation in metal alloys and unravel the interplay between crystal nucleation and glass transition, as the conditions of crystallization lie close to this transition. While decreasing the temperature of crystallization usually results in a lower free energy barrier, we find an unexpected reversal of behavior for glass-forming alloys as the temperature of crystallization approaches the glass transition. For this purpose, we simulate the crystallization process in two glass-forming Copper alloys, $Ag_6Cu_4$, which has a positive heat of mixing, and in $CuZr$, characterized by a large negative heat of mixing. Our results allow us to identify that this unusual behavior is directly correlated with a non-monotonic temperature dependence for the formation energy of connected icosahedral structures, which are incompatible with crystalline order and impede the development of the crystal nucleus, leading to an unexpectedly larger free energy barrier at low temperature. This, in turn, promotes the formation of a predominantly closed-packed critical nucleus, with fewer defects, thereby suggesting a new way to control the structure of the crystal nucleus, which is of key importance in catalysis.

Using nonequilibrium molecular dynamics (NEMD) simulations, we study the properties of supercooled liquids of Si under shear at T=1060K over a range of densities encompassing the low-density liquid (LDL) and high-density liquid (HDL) forms. This enables us to generate nonequilibrium steady-states of the LDL and HDL polymorphs, that remain stabilized in their liquid forms for as long as the shear is applied. This is unlike the LDL and HDL forms at rest, which are metastable under those conditions and, when at rest, rapidly undergo a transition towards the crystal, i.e. the thermodynamically stable equilibrium phase. In particular, through a detailed analysis of the structural and energetic features of the liquids under shear, we identify the range of densities, as well as the range of shear rates, that give rise to the two forms. We also show how the competition between shear and tetrahedral order impacts the two-body entropy in steady-states of Si under shear. These results open the door to new ways of utilizing shear to stabilize forms that are metastable at rest and can exhibit unique properties, since, for instance, experiments on Si have shown that HDL is metallic, with no band gap, while LDL is semimetallic, with a pseudogap.

Entropy has become increasingly central to characterize, understand and even guide assembly, self-organization and phase transition processes. In this work, we build on the analogous role of partition functions (or free energies) in isothermal ensembles and that of entropy in adiabatic ensembles. In particular, we show that the grand-isobaric adiabatic $(\mu,P,R)$ ensemble, or Ray ensemble, provides a direct route to determine the entropy. This allows us to follow the variations of entropy with the thermodynamic conditions and thus to explore phase transitions. We test this approach by carrying out Monte Carlo simulations on Argon and Copper in bulk phases and at phase boundaries and assess the reliability and accuracy of the method through comparisons with the results from flat-histogram simulations in isothermal ensembles and with the experimental data. Advantages of the approach are multifold and include the direct determination of the $\mu-P$ relation, without any evaluation of pressure via the virial expression, the precise control of the system size and of the number of atoms via the input value of $R$, and the straightforward computation of enthalpy differences for isentropic processes, which are key quantities to determine the efficiency of thermodynamic cycles. A new insight brought by these simulations is the highly symmetric pattern exhibited by both systems along the transition, as shown by scaled temperature-entropy and pressure-entropy plots.

It has been recently shown that nanobubbles exhibit a remarkable and unexpected stability. The lifetime of nanobubbles, formed either within liquids or on hydrophobic surfaces, can exceed by more than 10 orders of magnitude the theoretical expectation, which predicts an almost immediate dissolution due to the very high Laplace internal pressure in such small bubbles. This unexpected property of nanobubbles has made them leading candidates for energy applications, e.g. as high-pressure nanoreactors in fuel cells, and for biological systems, as transport systems for gas delivery to membranes and cells. Here we use molecular simulation to shed light on the molecular mechanisms accounting for the formation and stabilization of nanobubbles under an hydrophobic nanoconfinement. Using an entropic reaction coordinate, we elucidate the nucleation pathway and determine the formation free energy of nanobubbles in water confined in carbon nanotubes. We identify a critical volume for which the existence of nanobubbles is thermodynamically favored, as the free energy profile flattens around this critical volume, and mechanically favored, since the nanoconfined fluid pressure, along the nanotube axis, is positive at this juncture. We also show that the stabilization process is assisted by the hydrophobic nature of the carbon nanotube and by the formation of strong hydrogen bonds at the interface. Caroline Desgranges and Jerome Delhommelle

We study in this paper the time evolution of stock markets using a statistical physics approach. Each agent is represented by a spin having a number of discrete states $q$ or continuous states, describing the tendency of the agent for buying or selling. The market ambiance is represented by a parameter $T$ which plays the role of the temperature in physics. We show that there is a critical value of $T$, say $T_c$, where strong fluctuations between individual states lead to a disordered situation in which there is no majority: the numbers of sellers and buyers are equal, namely the market clearing. We have considered three models: $q=3$ ( sell, buy, wait), $q=5$ (5 states between absolutely buy and absolutely sell), and $q=\infty$. The specific measure, by the government or by economic organisms, is parameterized by $H$ applied on the market at the time $t_1$ and removed at the time $t_2$. We have used Monte Carlo simulations to study the time evolution of the price as functions of those parameters. Many striking results are obtained. In particular we show that the price strongly fluctuates near $T_c$ and there exists a critical value $H_c$ above which the boosting effect remains after $H$ is removed. This happens only if $H$ is applied in the critical region. Otherwise, the effect of $H$ lasts only during the time of the application of $H$. The second party of the paper deals with the price variation using a time-dependent mean-field theory. By supposing that the sellers and the buyers belong to two distinct communities with their characteristics different in both intra-group and inter-group interactions, we find the price oscillation with time., Comment: to appear in Physica A (2021)

We combine machine learning (ML) with Monte Carlo (MC) simulations to study the crystal nucleation process. Using ML, we evaluate the canonical partition function of the system over the range of densities and temperatures spanned during crystallization. We achieve this on the example of the Lennard-Jones system by training an artificial neural network using, as a reference dataset, equations of state for the Helmholtz free energy for the liquid and solid phases. The accuracy of the ML predictions is tested over a wide range of thermodynamic conditions, and results are shown to provide an accurate estimate for the canonical partition function, when compared to the results from flat-histogram simulations. Then, the ML predictions are used to calculate the entropy of the system during MC simulations in the isothermal-isobaric ensemble. This approach is shown to yield results in very good agreement with the experimental data for both the liquid and solid phases of Argon. Finally, taking entropy as a reaction coordinate and using the umbrella sampling technique, we are able to determine the Gibbs free energy profile for the crystal nucleation process. In particular, we obtain a free energy barrier in very good agreement with the results from previous simulation studies. The approach developed here can be readily extended to molecular systems and complex fluids, and is especially promising for the study of entropy-driven processes.

The thermodynamic Bethe ansatz equations of the Coqblin-Schrieffer model have been solved numerically. The full N=6 (J=5/2) degeneracy of the Hund's rule ionic ground state of Ce is taken into account. Results for the temperature dependent magnetic susceptibility parallel and perpendicular to the crystal axis are presented. The deviations, due to the Kondo effect, to the non-interacting ion results are pointed out., Comment: 4 pages, 2 figures, accepted for publication in Journal of Magnetism and Magnetic Materials

We describe a new method to automatically discriminate between patients with Alzheimer's disease (AD) or mild cognitive impairment (MCI) and elderly controls, based on multidimensional classification of hippocampal shape features. This approach uses spherical harmonics (SPHARM) coefficients to model the shape of the hippocampi, which are segmented from magnetic resonance images (MRI) using a fully automatic method that we previously developed. SPHARM coefficients are used as features in a classification procedure based on support vector machines (SVM). The most relevant features for classification are selected using a bagging strategy. We evaluate the accuracy of our method in a group of 23 patients with AD (10 males, 13 females, age $\pm$ standard-deviation (SD) = 73 $\pm$ 6 years, mini-mental score (MMS) = 24.4 $\pm$ 2.8), 23 patients with amnestic MCI (10 males, 13 females, age $\pm$ SD = 74 $\pm$ 8 years, MMS = 27.3 $\pm$ 1.4) and 25 elderly healthy controls (13 males, 12 females, age $\pm$ SD = 64 $\pm$ 8 years), using leave-one-out cross-validation. For AD vs controls, we obtain a correct classification rate of 94%, a sensitivity of 96%, and a specificity of 92%. For MCI vs controls, we obtain a classification rate of 83%, a sensitivity of 83%, and a specificity of 84%. This accuracy is superior to that of hippocampal volumetry and is comparable to recently published SVM-based whole-brain classification methods, which relied on a different strategy. This new method may become a useful tool to assist in the diagnosis of Alzheimer's disease., Comment: Data used in the preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database

In this note we determine the set of expansions such that a partial cube is planar if and only if it arises by a sequence of such expansions from a single vertex. This corrects a result of Peterin., Comment: 4 pages, 3 figures

The thermodynamic Bethe ansatz equations for the Coqblin-Schrieffer model have been solved for the first time to obtain the magnetic susceptibility in the presence of crystal fields for non-zero temperatures. For the case of N = 4 effective ionic states an analytic expression for the limiting values of the pseudo-energies has been found facilitating the numerical solution for various crystal and magnetic field configurations. The single-impurity model applies to a wide range of dense Kondo systems and has been used before to explain apparent non-Fermi-liquid behavior. The flattening off of the susceptibility curves at a substantially higher temperature than the specific heat is shown to be a general feature of the Coqblin-Schrieffer thermodynamics., Comment: 13 pages, 6 figures

Because of permanent use-dependent brain plasticity, all lifelong individuals' experiences are believed to influence the cognitive aging quality. In older individuals, both former and current musical practices have been associated with better verbal skills, visual memory, processing speed, and planning function. This work sought for an interaction between musical practice and cognitive aging by comparing musician and non-musician individuals for two lifetime periods (middle and late adulthood). Long-term memory, auditory-verbal short-term memory, processing speed, non-verbal reasoning, and verbal fluencies were assessed. In Study 1, measures of processing speed and auditory-verbal short-term memory were significantly better performed by musicians compared with controls, but both groups displayed the same age-related differences. For verbal fluencies, musicians scored higher than controls and displayed different age effects. In Study 2, we found that lifetime period at training onset (childhood vs. adulthood) was associated with phonemic, but not semantic, fluency performances (musicians who had started to practice in adulthood did not perform better on phonemic fluency than non-musicians). Current frequency of training did not account for musicians' scores on either of these two measures. These patterns of results are discussed by setting the hypothesis of a transformative effect of musical practice against a non-causal explanation.

The thermodynamic Bethe ansatz equations for the N = 6 Coqblin-Schrieffer model with crystal fields have been solved numerically. The realistic case of three Kramers doublets with arbitrary splittings has been studied for the first time. The specific heat has been calculated for representative combinations of the ionic energy splittings providing ample material for comparison with experimental results for Cerium impurities and compounds., Comment: 9 pages, 5 figures

Cracking is observed when a UO_{2} single crystal is oxidised in air. Previous studies led to the hypothesis that cracking occurs once a critical depth of U_{3}O_{7} oxidised layer is reached. We present some \mu-Laue X-ray diffraction results, which evidence that the U_{3}O_{7} layer, grown by topotaxy on UO_{2}, is made of domains with different crystalline orientations. This observation was used to perform a modelling of oxidation coupling chemical and mechanical parameters, which showed that the domain patterning induces stress localisation. This result is discussed in comparison with stress localisation observed in thin layer deposited on a substrate and used to propose an interpretation of UO_{2} oxidation and cracking.

In PWR, the Zircaloy based clad is the first safety barrier of the fuel rod, it must prevent the dispersion of the radioactive elements, which are formed by fission inside the UO2 pellets filling the clad. We focus here on internal corrosion that occurs when the clad is in tight contact with the UO2 pellet. In this situation, with temperature of 400^{\circ}C on the internal surface of the clad, a layer of oxidised Zircaloy is formed with a thickness ranging from 5 to 15 $\mu$m. In this paper, we will underline the specific behaviour of this internal corrosion layer compared to wet corrosion of Zircaloy. Simulations will underline the differences of stress field and their influences on corresponding dissolved oxygen profiles. The reasons for these differences will be discussed as function of the mechanical state at inner surface of the clad which is highly compressed. Differences between mechanical conditions generated by an inner or outer corrosion of the clad are studied and their influences on the diffusion phenomena are highlighted.

The growth of a U3O7 oxide layer during the anionic oxidation of UO2 pellets induced very important mechanical stresses due to the crystallographic lattice parameters differences between UO2 and its oxide. These stresses, combined with the chemical processes of oxidation, can lead to the cracking of the system, called chemical fragmentation. We study the crystallographic orientation of the oxide lattice growing at the surface of UO2, pointing the fact that epitaxy relations at interface govern the coexistence of UO2 and U3O7. In this work, several results are given: - Determination of the epitaxy relations between the substrate and its oxide thanks to the Bollmann's method; epitaxy strains are deduced. - Study of the coexistence of different domains in the U3O7 (crystallographic compatibility conditions at the interface between two phases: Hadamard conditions). - FEM simulations of a multi-domain U3O7 connected to a UO2 substrate explain the existence of a critical thickness of U3O7 layer.

The aim of the present work is to introduce a thermodynamic model to describe the growth of an oxide layer on a metallic substrate. More precisely, this paper offers a study of oxygen dissolution into a solid, and its consequences on the apparition of mechanical stresses. They strongly influence the oxidation processes and may be, in some materials, responsible for cracking. To realize this study, mechanical considerations are introduced into the classical diffusion laws. Simulations were made for the particular case of uranium dioxide, which undergoes the chemical fragmentation. According to our simulations, the hypothesis of a compression stress field into the oxidised UO2 compound near the internal interface is consistent with the interpretation of the mechanisms of oxidation observed experimentally.