Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Equation of State”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fluctuations in the entropy of Hawking radiation

We use the gravitational path integral (GPI) to compute the fluctuations of the Hawking radiation entropy around the Page curve in a two-dimensional model introduced by Penington Before the Page time, we find that δ S = e − S / 2 , where S is the black hole entropy. This result agrees with the Haar-averaged entropy fluctuations in a bipartite system. After the Page time, we find that δ S ∼ e − S , up to a prefactor that depends logarithmically on the width of the microcanonical energy window. This is not symmetric under exchange of subsystem sizes and so does not agree with the Haar average for a subsystem of fixed Hilbert space dimension. The discrepancy can be attributed to the fact that the black hole Hilbert space dimension is not fixed by the state preparation: even in a microcanonical ensemble with a top-hat smearing function, the GPI yields an additive fluctuation in the number of black hole states. This result, and the fact that the Page curve computed by the GPI is smooth, all point towards an ensemble interpretation of the GPI. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR reactor safety analysis. It aims to provide fast-running, modest-fidelity, whole-plant transient analyses capabilities, which are essential for fast turnaround design scoping and engineering analyses of advanced reactor concepts. While SAM is being developed as a system-level modeling and simulation tool, advanced modeling techniques being implemented include a reduced-order three-dimensional module, pseudo 3-D conjugate heat transfer modeling in reactor core, flexible and multi-scale modeling of heat transfer between fluid and structures, in addition to the advances in software environments and design, and numerical methods. SAM aims to be a generic system-level safety analysis tool for advanced non-LWRs, including Liquid-Metal-cooled fast Reactors (LMR), Molten Salt Reactors (MSR), Fluoride-salt-cooled High- temperature Reactors (FHR), and High-Temperature Gas-cooled Reactors (HTGR). SAM takes ad- vantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented computational framework (MOOSE), and its underlying meshing and finite-element library and linear and non-linear solvers, to leverage the modern advanced software environments and numerical methods. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As new code capabilities and features are added, the SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Nuclear Thermal Rocket Emulator for a Hardware-in-the-Loop Test Bed

To support NASA’s mission to use nuclear thermal rockets for future Mars missions, an instrumentation and control test bed has been built at Oak Ridge National Laboratory. The system is designed as a hardware-in-the-loop test bed for testing control elements and autonomous control algorithms for nuclear thermal propulsion rockets. The mock reactor system consists of a modular and scalable framework, using inexpensive components and open-source software. The hardware system consists of a two-phase flow loop and a mock reactor with six control drums. A single-board computer (NVIDIA Jetson) handles reactor core emulation and hosts a message queuing telemetry transport broker that allows user-deployed control algorithms to interact with the system hardware. The reactor emulator receives sensor data from the hardware and provides the simulated performance of the reactor under steady-state, transient, and fault conditions. The emulator uses a reactivity lookup table and the point kinetics equations to solve for the reactor dynamics in real time. Emulated reactor dynamics and sensor input inform the autonomous control algorithm’s decision-making in a closed-loop manner. The current system is capable of operating at 10 Hz, but faster cycle rates are an area of ongoing research. This test bed will enable NASA and other space vendors to rigorously test their autonomous control systems for NTP rockets under transient (reactor startup and shutdown), steady-state, and fault conditions to reduce development time and risk for autonomous control systems in future missions.

autonomous control↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Role of quadrupole deformation and continuum effects in the “island of inversion” nuclei 28,29,31 F

The peculiar properties of nuclei in the so-called “island of inversion” around Z = 10 and N = 20 are the focus of current nuclear physics research. Recent studies showed that 28 F has a negative-parity ground state and thus lies within the southern shore of the island of inversion, and 29 F presents a halo structure in its ground state, but it is unclear which effects, such as deformation, shell evolution due to tensor forces, or couplings to the continuum, lead to this situation. Here, we investigate the role of quadrupole deformation and continuum effects on the single-particle structure of 28,29,31 F from a relativistic mean-field approach and show how both phenomena can lead to a negative-parity ground state in 28 F and halo structures in 29,31 F. We solve the Dirac equation in the complex-momentum (Berggren) representation for a potential with quadrupole deformation at the first order obtained from relativistic mean-field calculations using the NL3 interaction and calculate the continuum level densities using the Green’s function method. We extract single-particle energies and widths from the continuum level densities to construct the Nilsson diagrams of 28,29,31 F in the continuum and analyze the evolution of both the widths and occupation probabilities of relevant Nilsson orbitals in 28 F and find that some amount of prolate deformation must be present. In addition, we calculate the density distributions for bound Nilsson orbitals near the Fermi surface in 29,31 F and reveal that, for a quadrupole deformation 0.3 ≤ β 2 ≤ 0.45 (prolate), characteristic halo tails appear at large distances. Using the relativistic mean-field approach in the complex-momentum representation with the Green’s function method, we demonstrate that in neutron-rich fluorine isotopes, while in the spherical case the pf shells are already inverted and close to the neutron emission threshold, a small amount of quadrupole deformation can dramatically reduce the gap between positive- and negative-parity states and increase the role of continuum states, ultimately leading to the negative parity in the ground state of 28 F and the halo structures in 29,31 F.

20 ≤ A ≤ 38↗
DOE OSTI2024

Bulk Viscosity in Dense Nuclear Matter

NEUTRON star mergers are an extreme environment where, within milliseconds, the cold, dense matter in neutron stars is dramatically heated to some of the hottest temperatures encountered in astrophysics. Matter in neutron star mergers is located solidly in the interior of the QCD phase diagram, and therefore, the observation and simulation of these mergers provides an opportunity to learn about the strong interaction at a set of densities and temperatures that compliment those obtained in relativistic heavy-ion collisions. Matter during the merger of two neutron stars is far from hydrostatic equilibrium, allowing for the study of not only the equation of state of dense matter, but of transport properties as well. Perhaps the most likely transport property to play a role in neutron star mergers is the bulk viscosity, which acts to resist changes in density. Furthermore, the bulk viscosity of matter at high densities has long been anticipated to damp oscillations in isolated neutron stars, and the study of its role in neutron star mergers is just beginning.

Harris, Steven P.↗
DOE OSTI2026

A limit on the total lepton number in the Universe from BBN and the CMB

At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

General relativistic hydrodynamic simulations of binary strange star mergers

We perform fully general-relativistic simulations of binary strange star mergers considering two different approaches for thermal effects. The first uses a cold equation of state (EOS) derived from a modified version of the MIT bag model which is then supplemented by a Γ-law correction. The second approach employs a microphysical description of the finite-temperature effects. We describe results obtained with the two treatments, highlighting the influence of thermal effects. We find that the postmerger dynamics differs significantly in the two cases, leading to quantitative differences in the postmerger gravitational-wave spectrum and ejecta mass. The peak frequency of the postmerger gravitational-wave emission is consistent with the established quasi-universal relations for binary neutron star mergers and as a result, our simulations cannot distinguish between mergers of neutron stars and those of strange stars. Our models with realistic treatment of finite-temperature effects produce a significant amount of ejecta ≳0.02 M ⊙ ​. Here, the resulting flux of strangelets near the Earth, computed assuming that all neutron star mergers are in fact strange-stars mergers and that the binary considered here is representative, is in tension with experimental upper limits. As such, our results tentatively disfavor a scenario in which strange-quark matter is the lowest energy state of matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Unveiling Baryon Charge Carriers through Charge Stopping in Isobar Collisions

Utilizing a comprehensive (3 + 1)⁢D relativistic hydrodynamic framework with multiple conserved charge currents and charge-dependent lattice-QCD-based equation of state, we study the baryon and electric charge number deposition at midrapidity in isobar Ru + Ru and Zr + Zr collisions at the center of mass energy $\sqrt{s}_{NN}$ = 200 GeV. Comparing our predictions with upcoming experimental data from the Relativistic Heavy Ion Collider will shed light on the existence of baryon junctions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A hybrid porous model for full reactor core scale CFD investigation of a prismatic HTGR

Three-dimensional (3-D) Computational Fluid Dynamics (CFD) analysis of a whole nuclear reactor core is a tremendous challenge due to the large geometric volume and complex structures. Here, this research presents a hybrid porous (HP) model to simplify a prismatic High Temperature Gas-cooled Reactor (HTGR) core, so 3-D CFD investigation can be performed on a full reactor core scale. In the HP model, the prototypic small coolant channels in the nuclear fuel blocks are lumped together to form multiple equivalent large coolant channels, and then the porous medium flow model is applied to each of them. Therefore, heat transfer in fuel blocks is computed by a hybrid combination of solid energy and porous flow energy equations. The similarity between the HP model and prototypic model is achieved by deriving the porous flow permeability, inertial resistance factor, and artificial thermophysical properties. Compared with the widely used whole porous (WP) flow model, the HP model preserves more realistic geometric structures, and therefore more accurate physical processes. The General Atomics' Modular High Temperature Gas-cooled Reactor (MHTGR) design was chosen as a prototype to demonstrate the methodology. Simulations were performed using the prototypic CFD model and HP model at steady-state forced circulation, steady-state natural circulation, and transient conditions that correspond to normal operation, extended period of pressurized cool down, and short-term transients after reactor shutdown, respectively. The comparison shows good agreement between the HP model and prototypic model in the maximum fuel temperature, average solid temperature, and helium flow rate, which demonstrates the potential applicability of the HP model for a full reactor core scale simulation in the future. As a benefit, the HP model reduces the mesh quantity by a factor of 50 from a prototypic model. Correspondingly, the computation time was reduced by a factor of at least 30.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Molecular-level insights into structure and dynamics in ionic liquids and polymer gel electrolytes

We report designing new electrolytes requires a better understanding of the correlation between their transport properties and their molecular structure. In this work, we present a detailed study of ionic liquids and polymer gel electrolytes probing their structure and dynamics by nuclear magnetic resonance (NMR) spectroscopy. In particular, ammonium- and phosphonium-based ionic liquids (ILs) are combined with different LiTFSI concentrations, and then with different ratios of poly(methylmethacrylate) polymer. The temperature dependence of self-diffusion coefficients of mobile species, D Li+ , D TFSI- and D P4441+ or D N4441+ , measured by pulsed field gradient (PFG) NMR spectroscopy obeys the Arrhenius equation. Diffusivity of [P4441][TFSI] ILs is found to be greater than those of the [N4441][TFSI] ILs in both liquid and gel electrolytes. Solid state NMR experiments including 13 C and 19 F MAS, 13 C{ 19 F}, 13 C{ 1 H} CPMAS probe the local structure and molecular-level interactions between ions and polymer in gel electrolytes, particularly for samples with higher PMMA content (≥25 wt%). Finally, the fast field cycling relaxometry has been used to unveil the rotational and translational dynamics of P4441 + or N4441 + and TFSI - by measuring 1 H and 19 F R 1 relaxation rate profiles at different temperatures. A comprehensive NMR analysis including relaxation studies at low magnetic field provides decisive new insights regarding the formation of ionic clusters and the interaction of ions with the polymer chain in the case of the gel electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A Parameterized Neutrino Emission Model to Study Mass Ejection in Failed Core-collapse Supernovae

Some massive stars end their lives as failed core-collapse supernovae (CCSNe) and become black holes (BHs). Although in this class of phenomena the stalled supernova (SN) shock is not revived, the outer stellar envelope can still be partially ejected. This occurs because the hydrodynamic equilibrium of the star is disrupted by the gravitational mass loss of the protoneutron star (PNS) due to neutrino emission. We develop a simple parameterized model that emulates PNS evolution and its neutrino emission and use it to simulate failed CCSNe in spherical symmetry for a wide range of progenitor stars. Our model allows us to study mass ejection of failed CCSNe where the PNS collapses into a BH within ~100 ms and up to ~10 6 s. We perform failed CCSNe simulations for 262 different pre-SN progenitors and determine how the energy and mass of the ejecta depend on progenitor properties and the equation of state (EOS) of dense matter. In the case of a future failed CCSN observation, the trends obtained in our simulations can be used to place constraints on the pre-SN progenitor characteristics, the EOS, and on PNS properties at BH formation time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Extension of the high-resolution thermal-hydraulics code ESCOT to hexagonal core geometries for multi-physics calculations

The extension of the capabilities of the pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT to analyze hexagonal fueled cores and its performance are presented. ESCOT is an accurate yet fast core thermal-hydraulics solution aiming at high-fidelity and high-resolution multi-physics core analysis in the framework of massively parallel computing platforms. Its algorithm solution is based on the four-equation drift-flux model for two-phase calculations, these are numerically solved by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models such as turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analysis. ESCOT is parallelized by a double (radial and axial) domain decomposition that enables its highly parallelized execution. The coupling of the code with the neutronics whole core solver for hexagonal geometries nTRACER is described. The newly implemented ESCOT features are validated by comparing single assembly and full core steady state nTRACER-ESCOT solutions with nTRACER standalone internal one-dimensional T/H solver results. The validation problems are based on the VVER 440 and VVER 1000 cores. ESCOT results show differences within an acceptable range with respect to the simple 1D nTRACER built-in solver. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Bottomonium suppression in an open quantum system using the quantum trajectories method

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Critical Unresolved Region Integral Experiment Execution

Integral Experiments are a key aspect of validating nuclear data behavior and simulation capabilities. Further, integral experiments are the earliest form of validation of equations and codes. They have been used extensively since the 1940s. Today, simulation and predictive capabilities have improved greatly from 75 years ago, but nonetheless, integral experiments are still needed. As simulation capabilities improve, their uncertainties get much smaller. The current focus of many new experiments is understanding the intermediate energy range. The intermediate energy range does not have a wide application space and it can be challenging to design experiments that are sensitive to neutron energies in this region. The Zeus experiments evaluated highly enriched uranium (HEU) in the intermediate energy region. Building on the success and knowledge gained from the Zeus experiments, the Critical Unresolved Region Integral Experiment (CURIE) experiments were designed to evaluate HEU in the narrower unresolved resonance region (URR). These experiments were executed during June and July 2020 at the National Criticality Experiments Research Center (NCERC). The National Criticality Experiments Research Center (NCERC), operated by Los Alamos National Laboratory, is the only general purpose critical experiments facility in the United States of America. NCERC regularly designs and executes critical experiments and other measurements useful to a wide variety of fields including nuclear criticality safety, commercial nuclear energy, nonproliferation, and nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Bose-Einstein condensation phenomenology in systems with repulsive interactions

The role of repulsive interactions in statistical systems of Bose particles is investigated. Three different phenomenological frameworks are considered: A mean-field model, an excluded volume model, and a model with a medium-dependent effective mass. All three models are tuned to yield similar equations of state, with only minor deviations from the ideal Bose gas at small chemical potentials. Our analysis indicates, however, that these models lead to qualitatively different results for the Bose-Einstein condensation phenomenon. We discuss the different aspects of this phenomenon, namely, an onset of the Bose-Einstein condensation, particle number fluctuations, and a behavior of the Bose condensate. The obtained results can be helpful for interpreting the lattice QCD data at small temperature and large isospin chemical potential and the data on multiple pion production in high-energy nuclear collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗