SEARCH · Engineering Papers

Results for “Phase field modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Mixeds-sourcery: Building many-body states using bubbles of nothing

We recently introduced the idea of s-sourcery [B. Swingle and J. McGreevy, Phys. Rev. B 93, 045127 (2016)], a general formalism for building many-body quantum ground states using renormalization-group-inspired quantum circuits. In this work, we define a generalized notion of s-sourcery that applies to mixed states and study its properties and applicability. For our examples we focus on thermal states of local Hamiltonians. We prove a number of theorems establishing the prevalence of mixed s-source fixed points, giving results for free fermion models, conformal field theories, holographic models, and topological phases. Thermal double states (also called thermofield double states) and the machinery of approximate conditional independence are used heavily in the constructions. For a large class of models we provide an information theoretic argument for the existence of a local Hamiltonian whose ground state is the thermal double state, and in some cases we construct such a Hamiltonian.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Comparative modeling of the disregistry and Peierls stress for dissociated edge and screw dislocations in Al

Many elementary deformation processes in metals involve the motion of dislocations. The planes of glide and specific processes dislocations prefer depend heavily on their atomic core structures. Atomistic simulations are desirable for dislocation modeling but their application to even sub-micron scale problems is in general computationally costly. Accordingly, continuum-based approaches, such as the phase-field microelasticity, phase-field dislocation dynamics (PFDD), generalized Peierls–Nabarro (GPN) models, and the concurrent atomistic–continuum (CAC) method, have attracted increasing attention in the field of dislocation modeling because they well represent both short-range cores interactions and long-range stress fields of dislocations. To better understand their similarities and differences, it is useful to compare these methods in the context of benchmark simulations and predictions. In this paper, we apply the CAC method and different PFDD variants – one of them is equivalent to a GPN model – to simulate an extended (i.e., dissociated) dislocation in Al with initially pure edge or pure screw character in terms of the disregistry. CAC and discrete forms of PFDD are also employed to calculate the Peierls stress. By conducting comprehensive convergence studies, we quantify the dependence of these measures on time/grid resolution and simulation cell size. Several important but often overlooked differences between PFDD/GPN variants are clarified. In conclusion, our work sheds light on the advantages and limitations of each method, as well as the path towards enabling them to effectively model complex dislocation processes at larger length scales.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Predicting fracture behavior in single crystal nickel using a coupled crystal plasticity phase field damage approach

Understanding the fracture behavior of single crystal metals is critical for predicting material performance under mechanical loading. Here, in this study, we investigate the fracture characteristics of single crystal nickel tensile bars using a crystal plasticity coupled phase field damage (CP-PFD) model. Experimental tensile tests were conducted on 15 specimens spanning five crystallographic orientations and three thickness variants per orientation. The results revealed two distinct fracture modes: brittle fractures with 45-degree angled surfaces and ductile fractures characterized by significant necking. The CP-PFD model successfully replicated these fracture behaviors, demonstrating strong agreement with experimental observations. The model effectively predicted the strain at which necking and fracture occurred, as well as the orientation-dependent fracture mechanisms. By comparing experimental and simulated fracture surfaces, we establish the CP-PFD model as a robust tool for predicting single crystal behavior and damage evolution. This work provides insight into the microstructural dependence of fracture behavior and establishes a predictive framework for modeling orientation-dependent damage evolution in single-crystal nickel.

Crystal plasticity↗

NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

DOE OSTI2011

Quasiparticles of strongly correlated Fermi liquids at high temperatures and in high magnetic fields

Strongly correlated Fermi systems are among the most intriguing, best experimentally studied and fundamental systems in physics. There is, however, lack of theoretical understanding in this field of physics. The ideas based on the concepts like Kondo lattice and involving quantum and thermal fluctuations at a quantum critical point have been used to explain the unusual physics. Alas, being suggested to describe one property, these approaches fail to explain the others. This means a real crisis in theory suggesting that there is a hidden fundamental law of nature. It turns out that the hidden fundamental law is well forgotten old one directly related to the Landau-Migdal quasiparticles, while the basic properties and the scaling behavior of the strongly correlated systems can be described within the framework of the fermion condensation quantum phase transition (FCQPT). The phase transition comprises the extended quasiparticle paradigm that allows us to explain the non-Fermi liquid (NFL) behavior observed in these systems. In contrast to the Landau paradigm stating that the quasiparticle effective mass is a constant, the effective mass of new quasiparticles strongly depends on temperature, magnetic field, pressure, and other parameters. Our observations are in good agreement with experimental facts and show that FCQPT is responsible for the observed NFL behavior and quasiparticles survive both high temperatures and high magnetic fields.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Implementation of a Pressure Correction Term in the CTF Subchannel Code

CTF is a modernized and improved version of the legacy subchannel code COBRA-TF, being jointly developed and maintained by Oak Ridge National Laboratory (ORNL) and North Carolina State University. The code was adopted by ORNL for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) in 2012 to help address CASL challenge problems. Since that time, activities related to CTF have included implementing software quality assurance measures, implementing new features and models, performing validation and verification testing, establishing and supporting a CTF user group, and developing the code for use in coupled applications. CTF is applicable to single- and two-phase flows in light water reactor geometries at normal and accident operating conditions. The two-phase flow model consists of the six-equation two-phase flow model augmented by equations for the droplet field, and it also includes appropriate source terms. This work focuses on the numerical stability of the six-equation two-phase flow model and is part of a recent effort to improve CTF capabilities for simulating flows in boiling water reactors (BWRs). The six-equation two-phase flow model is not strongly hyperbolic, meaning that some of its eigenvalues are either not real or not distinct. In a Cauchy sense, and for an initial-value problem, strong hyperbolicity was found to be equivalent to well-posedness, which ensures the uniqueness of the solution in the space-time used. This concept is also of importance when designing high-order numerical methods and also when implementing boundary conditions that rely on real eigenvalues such as Riemann solvers. A tremendous amount of research is available in the literature on how to hyperbolize the six-equation two-phase flow model by adding ad hoc terms to the system of equations. Such terms are often referred to as virtual mass or interfacial pressure correction terms. From a mathematical perspective, this approach adds terms to the phasic momentum equation containing first-order partial derivatives of the unknowns to modify the Jacobian matrix, thus the imaginary part of the eigenvalues. For this work, an interfacial pressure correction (IPC) term on the same model as what was used for the CATHARE system code, is implemented in CTF. This allows for recovery of a two-pressure formulation. In the following, the six-equation two-phase flow model is recalled, along with the IPC term and its influence on the eigenvalues. Then, numerical results and convergence studies are presented for a water-faucet problem at a pressure of 10 bar. Finally, conclusions and future work are addressed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗

DOE OSTI2013

Phase transitions in quark matter under strong magnetic fields

In this work we use de SU(2) Nambu-Jona-Lasinio model to study the chiral transition at finite temperature, chemical potential and magnetic field. We show how the magnetic field affects the location of the critical end-point in the phase diagram, the constituent quark masses and the spinodal lines concerning the first order transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗

DOE OSTI2011

A Micro-Scale Model for Oxygen Reduction on LSM-YSZ Cathode

In this study, a micro-scale model is developed to simulate the oxygen reduction on LSM-YSZ composite cathode. The model incorporates the effects of cathode microstructural properties on the local transport phenomena and electrochemistry inside the cathode. A detailed reaction mechanism is used in the model which has two parallel routes for oxygen conversion into oxide ions, namely two-phase boundary and three-phase boundary pathways. The model predicts field distributions of local thermodynamic values, over-potential, Faradaic current and other parameters relevant to cathode performance. Electrochemical impedance simulations are performed using the current model to analyze the contribution of various processes to the overall impedance.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2021

C-FOG: Life of Coastal Fog

C-FOG is a comprehensive bi-national project dealing with the formation, persistence, and dissipation (life cycle) of fog in coastal areas (coastal fog) controlled by land, marine, and atmospheric processes. Given its inherent complexity, coastal-fog literature has mainly focused on case studies, and there is a continuing need for research that integrates across processes (e.g., air–sea–land interactions, environmental flow, aerosol transport, and chemistry), dynamics (two-phase flow and turbulence), microphysics (nucleation, droplet characterization), and thermodynamics (heat transfer and phase changes) through field observations and modeling. Central to C-FOG was a field campaign in eastern Canada from 1 September to 8 October 2018, covering four land sites in Newfoundland and Nova Scotia and an adjacent coastal strip transected by the Research Vessel Hugh R. Sharp . An array of in situ, path-integrating, and remote sensing instruments gathered data across a swath of space–time scales relevant to fog life cycle. Satellite and reanalysis products, routine meteorological observations, numerical weather prediction model (WRF and COAMPS) outputs, large-eddy simulations, and phenomenological modeling underpin the interpretation of field observations in a multiscale and multiplatform framework that helps identify and remedy numerical model deficiencies. An overview of the C-FOG field campaign and some preliminary analysis/findings are presented in this paper.

58 GEOSCIENCES↗

DOE OSTI2019

Modeling Sintering Processes of Nanoparticle Inks

The In-Pile Instrumentation program conducts research on advanced manufacturing of sensors for use within nuclear reactor experiments. The process for making these sensors involves the sintering of nanoparticles within solvent. The nanoparticles must be sintered in a way that minimizes pores and maximizes continuity. It has also been observed that the nanoparticle films may fall off the substrate, so maximizing adhesion is also important. To model the sintering process, the phase field method was employed for one and two nanoparticle - substrate systems and the contact angle was measured and compared to existing experimental and theoretical studies.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Electrochemically-Driven Phase Transformations in Battery Compounds (Final Technical Report)

Ion-insertion compounds for advanced batteries frequently exhibit phase transformations as the concentration of the working ion varies. Under the large electrochemical driving forces inherent to practical use, systems are often driven far from equilibrium and exhibit phase transition behavior not seen in other materials. An improved understanding of phase transformation pathways in electrode materials upon cycling will lead to new materials design concepts and electrochemical duty cycle management strategies that improve capacity utilization at high charge/discharge rates, reduce voltage and capacity hysteresis, and extend battery life. The main goal of this project is to develop a predictive understanding of the phase transition behavior of battery compounds when electrochemically driven far from equilibrium through a combined experimental- theoretical approach. Research focuses on elucidating how the phase transition behaviors in a diverse group of selected model systems (olivines, lithium metal, etc.) are regulated by various factors including transformation strains, plasticity, metastable transformation pathways, electrode microstructure, surface reaction and ion diffusion kinetics, and what are the unique features of phase transitions in battery compounds. To this end, operando X-ray-based characterization techniques are applied to interrogate the structure and composition evolution of electrode materials and interfaces at nano- and meso-scales. Mesoscale modeling techniques such as phase-field simulation are employed to shed light on the experimental observations and establish a unified framework for describing phase transition behaviors in battery compounds.

25 ENERGY STORAGE↗

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↗

NASA NTRS1999

The Dynamics of Vapor Bubbles in Acoustic Pressure Fields

In spite of a superficial similarity with gas bubbles, the intimate coupling between dynamical and thermal processes confers to oscillating vapor bubbles some unique characteristics. This paper examines numerically the validity of some asymptotic-theory predictions such as the existence of two resonant radii and a limit size for a given sound amplitude and frequency. It is found that a small vapor bubble in a sound field of sufficient amplitude grows quickly through resonance and continues to grow thereafter at a very slow rate, seemingly indefinitely. Resonance phenomena therefore play a role for a few cycles at most, and reaching a limit size-if one exists at all-is found to require far more than several tens of thousands of cycles. It is also found that some small bubbles may grow or collapse depending on the phase of the sound field. The model accounts in detail for the thermo-fluid-mechanic processes in the vapor. In the second part of the paper, an approximate formulation valid for bubbles small with respect to the thermal penetration length in the vapor is derived and its accuracy examined, The present findings have implications for acoustically enhanced boiling heat transfer and other special applications such as boiling in microgravity.

Hao, Y.↗

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↗