Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase Field Model”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phase Field Sintering Simulations of Tagged UO 2

Isotopic taggants are being considered as an additive for UO 2 manufacturing to improve nuclear security. Before taggants can be adopted, it must be determined whether taggants have any effect on the fuel microstructure or performance. As part of this effort, a grand potential formulation of the phase field method was used to model untagged, Cr-tagged, and Fe-tagged UO 2 . The simulations were compared to determine whether the taggants had any effect on fuel density or grain size. The results suggest that the taggants will not affect the density or grain size, but these findings may not be reliable because of high uncertainty. Causes of the uncertainty were discussed, and future work is proposed to reduce the model uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

High-fidelity multiphysics simulation of BWR assembly with coupled TORT-TD/CTF

This paper describes the application of the coupled codes TORT-TD and CTF to the pin-by-pin modeling of a BWR fuel assembly with thermal-hydraulic feedback. TORT-TD, developed at GRS, is a time-dependent three dimensional discrete ordinates code. CTF is the PSU's improved version of the subchannel code COBRA-TF, which uses a two-fluid, three-field model to represent two-phase flow with entrained droplets, and is commonly applied to evaluate LWR safety margins. The coupled codes system TORT-TD/CTF, already applied to several PWR cases involving MOX, was adapted from PWR to BWR applications. The purpose of the research described in this paper is to verify the coupling for modeling two-phase flow at the pin cell level. Using an ATRIUM-10 assembly, the system's steady-state capabilities were tested on two cases: one without control blade insertion and another with partially inserted blades. The influence of the neutron absorber on local axial and radial parameters is presented. The description of an inlet flow reduction transient is an example for the time-dependent capability of the coupled system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Quantum simulation of discrete curved spacetime by the Bose–Hubbard model: From analog acoustic black hole to quantum phase transition

We present a theoretical scheme to simulate quantum field theory in a discrete curved spacetime based on the Bose–Hubbard model describing a Bose–Einstein condensate trapped inside an optical lattice. Using the Bose–Hubbard Hamiltonian, we first introduce a hydrodynamic presentation of the system evolution in discrete space. We then show that the phase (density) fluctuations of the trapped bosons inside an optical lattice in the superfluid (Mott insulator) state obey the Klein–Gordon equation for a massless scalar field propagating in a discrete curved spacetime. We derive the effective metrics associated with the superfluid and Mott-insulator phases and, in particular, we find that in the superfluid phase the metric exhibits a singularity which can be considered as the manifestation of an analog acoustic black hole. The proposed approach is found to provide a suitable platform for quantum simulation of various spacetime metrics through adjusting the system parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2017

Multicanonical sampling of the space of states of ℋ(2, n)-vector models

Problems of temperature behavior of specific heat are solved by the entropy simulation method for Ising models on a simple square lattice and a square spin ice (SSI) lattice with nearest neighbor interaction, models of hexagonal lattices with short-range (SR) dipole interaction, as well as with long-range (LR) dipole interaction and free boundary conditions, and models of spin quasilattices with finite interaction radius. It is established that systems of a finite number of Ising spins with LR dipole interaction can have unusual thermodynamic properties characterized by several specific-heat peaks in the absence of an external magnetic field. For a parallel multicanonical sampling method, optimal schemes are found empirically for partitioning the space of states into energy bands for Ising and SSI models, methods of concatenation and renormalization of histograms are discussed, and a flatness criterion of histograms is proposed. It is established that there is no phase transition in a model with nearest neighbor interaction on a hexagonal lattice, while the temperature behavior of specific heat exhibits singularity in the same model, in case of LR interaction. A spin quasilattice is found that exhibits a nonzero value of residual entropy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Engineering infinite-range SU($\textit{n}$) interactions with spin-orbit-coupled fermions in an optical lattice

Here, we study multilevel fermions in an optical lattice described by the Hubbard model with on-site SU($\textit{n}$) -symmetric interactions. We show that in an appropriate parameter regime this system can be mapped onto a spin model with all-to-all SU($\textit{n}$)-symmetric couplings. Raman pulses that address internal spin states modify the atomic dispersion relation and induce spin-orbit coupling, which can act as a synthetic inhomogeneous magnetic field that competes with the SU($\textit{n}$) exchange interactions. We investigate the mean-field dynamical phase diagram of the resulting model as a function of $\textit{n}$ and different initial configurations that are accessible with Raman pulses. Consistent with previous studies for $\textit{n}$ = 2 , we find that for some initial states the spin model exhibits two distinct dynamical phases that obey simple scaling relations with $\textit{n}$. Moreover, for $\textit{n}$ > 2 we find that dynamical behavior can be highly sensitive to initial intraspin coherences. Our predictions are readily testable in current experiments with ultracold alkaline-earth-metal(-like) atoms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Competition between Abelian and Zeeman magnetic field effects in a two dimensional ultracold gas of fermions

The ground state of ultracold fermions in the presence of effects of orbital and Zeeman magnetic fields is analyzed. Five different states are found: unpolarized superconducting state, partially and fully polarized normal states and phase separated regions, partially or fully polarized. The system, in the presence of orbital synthetic magnetic field effects, shows non-monotonous changes of the phase boundaries when electron concentration is varied. We observe not only reentrant phenomena, but also density dependent oscillations of different areas of the phase diagram. Moreover the chemical potential shows oscillatory behavior and discontinuities with respect to changes in the number of fermions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

Primordial black holes and magnetic fields in conformal neutrino mass models

Sufficiently strong and long-lasting first-order phase transitions can produce primordial black holes (PBHs) that contribute substantially to the dark matter abundance of the Universe, and can produce large-scale primordial magnetic fields. Here, we study these mechanisms in a generic class of conformal U(1)' models that also explain active neutrino oscillation data via the type-I seesaw mechanism. We find that phase transitions that occur at seesaw scales between 10 4 GeV and 10 11 GeV produce gravitational wave signals (from the dynamics of the phase transition and from the decay of cosmic string loops) at LISA/ET that can be correlated with microlensing signals of PBHs at the Roman Space Telescope, while scales near 10 11 GeV can be correlated with Hawking evaporation signals at future gamma-ray telescopes. LISA can probe the entire range of PBH masses between 1 × 10 -16 M ⊙ and 8 × 10 -11 M ⊙ if PBHs fully account for the dark matter abundance. For Z' masses between 40 TeV and 10 4 TeV, and 10 TeV right-handed neutrinos, helical magnetic fields can be produced with magnitudes ≳ 0.5 pG and coherence lengths ≳ 0.008 Mpc, above current blazar lower bounds.

cosmological phase transitions↗
DOE OSTI2011

Modeling PCM-Enhanced Insulation System and Benchmarking EnergyPlus against Controlled Field Data

Phase-change materials (PCM) used in building envelopes appear to be a promising technology to reduce energy consumption and reduce/shift peak load. However, due to complexity in modeling the dynamic behavior of PCMs, current modeling tools either lack an accurate way of predicting the performance and impact of PCMs in buildings or validation of predicted or measured performance is not available. This paper presents a model of a PCM-enhanced dynamic-insulation system in EnergyPlus (E+) and compares the simulation results against field-measured data. Laboratory tests to evaluate thermal properties and to characterize the PCM and PCM-enhanced cellulose insulation system are also presented in this paper. Results indicate that the predicted daily average heat flux through walls from the E+ simulation was within 9% of field measured data. Future analysis will allow us to predict annual energy savings from the use of PCM in buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Integrated Lift- and Wall-Force Closures in Computational Fluid Dynamics (CMFD)

In this report, we provide the necessary background to put this work in context and for a more detailed background, the reader is referred to CASL 2013 report. Central constraints on the design of LWRs are thermal-structural limits, e.g., maximum fuel pin temperature and maximum allowable working pressure in a pressure vessel. In both of these cases, the knowledge of the spatially and temporally varying two-phase flow is required to make accurate predictions as the distribution of both the phases have an impact on heat extraction as well as neutron moderation. Engineered two-phase flow systems have complex geometry and often high Reynolds numbers which increase the computational cost of an analysis, and generally make direct numerical simulations of the flow field impossible even on the fastest computers today. Consequently, a number of models which seek to reduce this computational cost have been created over the years. One of these called the two-fluid model has become the standard in three dimensional CFD (Computational Fluid Dynamics) codes and is and is also referred to as CMFD (Computational Multiphase Fluid Dynamics). The two-fluid model splits each conserved quantity into two fields, one for each phase that is treated as interspersed continua. Splitting one field into two requires additional interfacial closure relationships. In the two-fluid model, the momentum flux between the phases is governed by a number of these closure relationships. Generally, the net force on the dispersed phase is decomposed into inertial, added mass, buoyancy, drag, lift, and wall forces. In addition, there is a time-dependent Basset force, as well as turbulent dispersion effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

A multi-species exchange model for fully fluctuating polymer field theory simulations

Field-theoretic models have been used extensively to study the phase behavior of inhomogeneous polymer melts and solutions, both in self-consistent mean-field calculations and in numerical simulations of the full theory capturing composition fluctuations. The models commonly used can be grouped into two categories, namely, species models and exchange models. Species models involve integrations of functionals that explicitly depend on fields originating both from species density operators and their conjugate chemical potential fields. In contrast, exchange models retain only linear combinations of the chemical potential fields. In the two-component case, development of exchange models has been instrumental in enabling stable complex Langevin (CL) simulations of the full complex-valued theory. No comparable stable CL approach has yet been established for field theories of the species type. Here, we introduce an extension of the exchange model to an arbitrary number of components, namely, the multi-species exchange (MSE) model, which greatly expands the classes of soft material systems that can be accessed by the complex Langevin simulation technique. We demonstrate the stability and accuracy of the MSE-CL sampling approach using numerical simulations of triblock and tetrablock terpolymer melts, and tetrablock quaterpolymer melts. This method should enable studies of a wide range of fluctuation phenomena in multiblock/multi-species polymer blends and composites.

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

General monogamy property of global quantum discord and the application

We provide a family of general monogamy inequalities for global quantum discord (GQD), which can be considered as an extension of the usual discord monogamy inequality. It can be shown that those inequalities are satisfied under the similar condition for the holding of usual monogamy relation. We find that there is an intrinsic connection among them. Furthermore, we present a different type of monogamy inequality and prove that it holds under the condition that the bipartite GQDs do not increase when tracing out some subsystems. We also study the residual GQD based on the second type of monogamy inequality. As applications of those quantities, we investigate the GQDs and residual GQD in characterizing the quantum phase transition in the transverse field Ising model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Imprints of the QCD phase transition on the spectrum of gravitational waves

We have investigated effects of the QCD phase transition on the relic GW spectrum applying several equations of state for the strongly interacting matter: Besides the bag model, which describes a first-order transition, we use recent data from lattice calculations featuring a crossover. Finally, we include a short period of inflation during the transition which allows for a first-order phase transition at finite baryon density. Our results show that the QCD transition imprints a step into the spectrum of GWs. Within the first two scenarios, entropy conservation leads to a step-size determined by the relativistic degrees of freedom before and after the transition. The inflation of the third scenario much stronger attenuates the high frequency modes: An inflationary model being consistent with observation entails suppression of the spectral energy density by a factor of < or approx. 10{sup -12}.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Topological Mott transition in a two band model of spinless fermions with on-site Coulomb repulsion

Highlights: • The effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling lead to the 2D topological semimetal. • The ground state phase diagram is calculated. • The topological semimetal phase is characterized by the zero energy Majorana states. In the framework of mean field approach, we study topological Mott transition in a two band model of spinless fermions on a square lattice at half filling. We consider the combined effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling. The ground state phase diagram is calculated as a function of the strength of the spin–orbit Rashba coupling and the Coulomb repulsion. The spin–orbit Rashba coupling leads to a distinct phase of matter, the topological semimetal. We study a new type of phase transition between the non-topological insulator and topological semimetal states. Topological phase state is characterized by the zero energy Majorana states, which there are in defined region of the wave vectors and are localized at the boundaries of the sample. The region of existence of the zero energy Majorana states tends to zero at the point of the Mott phase transition. The zero energy Majorana states are dispersionless (they can be considered as flat bands), the Chern number and Hall conductance are equal to zero (note in two dimensional model).

36 MATERIALS SCIENCE↗
DOE OSTI2012

Lattice Modeling for SPEAR3

We use measured or simulated magnetic fields for dipoles and quadrupoles to build a lattice model for SPEAR3. In a non-symplectic approach the phase space coordinate mapping on the fields is based on numerical integration of the equation of motion. In a symplectic approach we approximate the fields with proper fringe field models. Complication of the use of rectangular gradient dipoles in SPEAR3 is considered. Results of the model is compared to measurements on the real machine.

43 PARTICLE ACCELERATORS↗