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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2010

Phase transitions in a gas of anyons

We continue our numerical Monte Carlo simulation of a gas of closed loops on a 3 dimensional lattice, however, now in the presence of a topological term added to the action which corresponds to the total linking number between the loops. We compute the linking number using a novel approach employing certain notions from knot theory. Adding the topological term converts the particles into anyons. Interpreting the model as an effective theory that describes the 2+1-dimensional Abelian Higgs model in the asymptotic strong-coupling regime, the topological linking number simply corresponds to the addition to the action of the Chern-Simons term. The system continues to exhibit a phase transition as a function of the vortex mass as it becomes small. We find the following new results. The Chern-Simons term has no effect on the Wilson loop. On the other hand, it does effect the 't Hooft loop of a given configuration, adding the linking number of the 't Hooft loop with all of the dynamical vortex loops. We find the unexpected result that both the Wilson loop and the 't Hooft loop exhibit a perimeter law even though there are no massless particles in the theory, in both phases of the theory. It should be noted that our method suffers from numerical instabilities if the coefficient of the Chern-Simons term is too large; thus, we have restricted our results to small values of this parameter. Furthermore, interpreting the lattice loop gas as an effective theory describing the Abelian Higgs model is only known to be true in the infinite coupling limit; for strong but finite coupling this correspondence is only a conjecture, the validity of which is beyond the scope of this article.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Graphics processing unit accelerated phase field dislocation dynamics: Application to bi-metallic interfaces

We present the first high-performance computing implementation of the meso-scale phase field dislocation dynamics (PFDD) model on a graphics processing unit (GPU)-based platform. The implementation takes advantage of the portable OpenACC standard directive pragmas along with Nvidia's compute unified device architecture (CUDA) fast Fourier transform (FFT) library called CUFFT to execute the FFT computations within the PFDD formulation on the same GPU platform. The overall implementation is termed ACCPFDD-CUFFT. The package is entirely performance portable due to the use of OPENACC-CUDA inter-operability, in which calls to CUDA functions are replaced with the OPENACC data regions for a host central processing unit (CPU) and device (GPU). A comprehensive benchmark study has been conducted, which compares a number of FFT routines, the Numerical Recipes FFT (FOURN), Fastest Fourier Transform in the West (FFTW), and the CUFFT. The last one exploits the advantages of the GPU hardware for FFT calculations. The novel ACCPFDD-CUFFT implementation is verified using the analytical solutions for the stress field around an infinite edge dislocation and subsequently applied to simulate the interaction and motion of dislocations through a bi-phase copper-nickel (Cu–Ni) interface. It is demonstrated that the ACCPFDD-CUFFT implementation on a single TESLA K80 GPU offers a 27.6X speedup relative to the serial version and a 5X speedup relative to the 22-multicore Intel Xeon CPU E5-2699 v4 @ 2.20 GHz version of the code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Coupled simulation of the reactor core using CUPID/MASTER

The CUPID is a component-scale thermal hydraulics code which is aimed for the analysis of transient two-phase flows in nuclear reactor components such as the reactor vessel, steam generator, containment. This code adopts a three-dimensional, transient, two-phase and three-field model, and includes physical models and correlations of the interfacial mass, momentum, and energy transfer for the closure. In the present paper, a multi-physics simulation was performed by coupling CUPID with a three dimensional neutron kinetics code, MASTER. MASTER is merged into CUPID as a dynamic link library (DLL). The APR1400 reactor core during a control rod drop/ejection accident was simulated as an example by adopting a porous media approach to employ a fuel assembly. The following sections present the numerical modeling for the reactor core, coupling of the kinetics code, and the simulation results. And also, a preliminary study for multi-scale simulation between CUPID and system-scaled thermal hydraulics code, MARS will be introduced as well. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2010

Exploring the role of model parameters and regularization procedures in the thermodynamics of the PNJL model

The equation of state and the critical behavior around the critical end point are studied in the framework of the Polyakov-Nambu-Jona-Lasinio model. We prove that a convenient choice of the model parameters is crucial to get the correct description of isentropic trajectories. The physical relevance of the effects of the regularization procedure is insured by the agreement with general thermodynamic requirements. The results are compared with simple thermodynamic expectations and lattice data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Numerical Analysis of Molten Pool Behavior and Spatter Formation with Evaporation During Selective Laser Melting of 316L Stainless Steel

During the selective laser-melting process, material evaporation and resultant spatter are common phenomena that bring about many defects. However, the underlying physical phenomena such as molten pool behavior and spatter formation, as evaporation occurred, are sparsely understood and difficult to observe during the process. Thus, a three–dimensional powder-scale model was established to investigate the thermal and flow behavior of the molten pool, the morphology evolution of the molten pool and keyhole, and the spatter formation with evaporation in the selective laser-melting processing of 316L stainless steel. Phase transitions and variations in interfacial force were taken into account in this model. The modified phase-field method was applied to trace the melt–gas interface. The results showed that keyhole formed in molten pool under recoil pressure, and that there were some differences in the temperature distribution and flow behavior inside and outside the keyhole. In addition, the dimension and surface morphology of the molten pool and keyhole depth altered and gradually stabilized during the process. Moreover, droplet spatter formation at rear of molten pool proceeded roughly as follows: the backward ejection of melt vapor would induce a depression in molten pool and a bump on the molten pool surface, the bump moved backward and subsequently a liquid column formed, then the liquid column formed a neck and gradually pinched off, resulting in the droplet spatter. Furthermore, the characteristic morphology of the scan track and dimension of the molten pool were obtained through experiments, which showed good agreement with the simulation results.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Quantum phase transitions and spontaneous symmetry-breaking in Dicke model

A method to find the Excited-States Quantum Phase Transitions (ESQPT's) from paritysymmetry in the Dicke model is studied and presented. This method allows us to stablish a critical energy where ESQPT's take places, and divides the whole energy spectrum in two regions with different properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Assimilation of Wave Imaging Radar Observations for Real-time Wave-by-Wave Forecasting

This project addressed Topic 3: “Wave Measurement Instrumentation for Feed Forward Controls” under the FOA number DE-FOA-0000971. The overall goal of the program was to develop a phase-resolving wave forecasting technique for application to the active control of Wave Energy Conversion (WEC) devices. We have developed an approach that couples a wave imaging marine radar with a phase-resolving linear wave model for real-time wave field reconstruction and forward propagation of the wave field in space and time. The scope of the project was to develop and assess the performance of this novel forecasting system. Specific project goals were as follows: Develop and verify a fast, GPU-based (Graphical Processing Unit) wave propagation model suitable for phase-resolved computation of nearshore wave transformation over variable bathymetry; Compare the accuracy and speed of performance of the wave model against a deep water model in their ability to predict wave field transformation in the intermediate water depths (50 to 70 m) typical of planned WEC sites; Develop and implement a variational assimilation algorithm that can ingest wave imaging radar observations and estimate the time-varying wave conditions offshore of the domain of interest such that the observed wave field is best reconstructed throughout the domain and then use this to produce model forecasts for a given WEC location; Collect wave-resolving marine radar data, along with relevant in situ wave data, at a suitable wave energy test site, apply the algorithm to the field data, assess performance, and identify any necessary improvements; and Develop a production cost estimate that addresses the affordability of the wave forecasting technology and include in the Final Report. The developed forecasting algorithm (“Wavecast”) was evaluated for both speed and accuracy against a substantial synthetic dataset. Early in the project, performance tests definitively demonstrated that the system was capable of forecasting in real-time, as the GPU-based wave model backbone was very computationally efficient. The data assimilation algorithm was developed on a polar grid domain in order to match the sampling characteristics of the observation system (wave imaging marine radar). For verification purposes, a substantial set of synthetic wave data (i.e. forward runs of the wave model) were generated to be used as ground truth for comparison to the reconstructions and forecasts produced by Wavecast. For these synthetic cases, Wavecast demonstrated very good accuracy, for example, typical forecast correlation coefficients were between 0.84-0.95 when compared to the input data. Dependencies on shadowing, observational noise, and forecast horizon were also identified. During the second year of the project, a short field deployment was conducted in order to assess forecast accuracy under field conditions. For this, a radar was installed on a fishing vessel and observations were collected at the South Energy Test Site (SETS) off the coast of Newport, OR. At the SETS site, simultaneous in situ wave observations were also available owing to an ongoing field project funded separately. Unfortunately, the position and heading information that was available for the fishing vessel were not of sufficient accuracy in order to validate the forecast in a phase-resolving sense. Instead, a spectral comparison was made between the Wavecast forecast and the data from the in situ wave buoy. Although the wave and wind conditions during the field test were complex, the comparison showed a promising reconstruction of the wave spectral shape, where both peaks in the bimodal spectrum were represented. However, the total reconstructed spectral energy (across all directions and frequencies) was limited to 44% of the observed spectrum. Overall, wave-by-wave forecasting using a data assimilation approach based on wave imaging radar observations and a physics-based wave model shows promise for short-term phase-resolved predictions. Two recommendations for future work are as follows: first, we would recommend additional focused field campaigns for algorithm validation. The field campaign should be long enough to capture a range of wave conditions relevant to the target application and WEC site. In addition, it will be crucial to make sure the vessel of choice has high accuracy position and heading instrumentation (this instrumentation is commercially available but not standard on commercial fishing vessels). The second recommendation is to expand the model physics in the wave model backbone to include some nonlinear effects. Specifically, the third-order correction to the wave speed due to amplitude dispersion would be the next step in order to more accurately represent the phase speeds of large amplitude waves.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Quantification of Dead Lithium on Graphite Anode under Fast Charging Conditions

A series of computational and experimental studies were conducted to understand the onset of lithium plating and subsequent quantification of dead lithium on graphite electrodes in the design of fast charging batteries. The experiments include titration and relaxation studies for detecting initiation of lithium metal plating for various SOC and C-rates, which are compared against the thermodynamically consistent phase field computational results. The collaborative study on "model graphite electrode" with 2.18 mAh cm-2 nominal capacity at 25 degrees C demonstrates: (1) the macroscopic voltage response during relaxation studies indicate the reintercalation of plated lithium into the graphite anode; (2) for SOC below 60% and low C-Rates, there is no dead lithium; (3) for SOC between 60% to 80%, and C-Rates in the range of 4C-6C show dead lithium both in experiments and simulations; (4) at 100% SOC and 4C-6C rates, large amounts of dead lithium are observed. The study presented here allows us to evaluate the effects of the physical properties of the electrochemical system on plating and stripping kinetics and the amount of dead lithium on graphite electrodes, which determines the cell capacity loss under fast charge.

batteries↗
DOE OSTI2025

Extending Interfaces in 3D to Achieve Superior Nanoscale Strength in Ti/Nb Nanolaminates

Tuning the atomic-level structure of nanolaminates enables high strength, increased deformability, and the ability to absorb and mitigate damage due to varied crystalline defects, including dislocations. Here, we present the enhanced strength of Ti/Nb nanolaminates containing thick 3D interfaces (3DIs), relative to their chemically abrupt 2D counterparts. We examine the effects of crystallographic alignment and compositional gradients on mechanical behavior via experiments and phase-field-dislocation dynamics (PFDD) modeling. Mechanical testing reveals that nanolaminates containing thicker 3DIs demonstrate a 28% hardness enhancement compared with sharp-interface nanolaminates. PFDD modeling shows that the critical resolved shear stress (CRSS) increases with the 3DI thickness. Gradual compositional transitions in 3DIs were confirmed via scanning transmission electron microscopy and high-resolution transmission electron microscopy, showing sharp crystallographic transitions and a heightened interface topography. The findings establish a positive function between the 3DI thickness and mechanical robustness for hexagonal-closest-packed-containing composites, emphasizing the role of defect–interface interactions in tailoring the mechanical performance and providing a foundation for future interfacial engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultrasonic oscillatory two-phase flow in microchannels

Experimental and numerical investigations are performed to provide an assessment of the transport behavior of an ultrasonic oscillatory two-phase flow in a microchannel. The work is inspired by the flow observed in an innovative ultrasonic fabric drying device using a piezoelectric bimorph transducer with microchannels, where a water-air two-phase flow is transported by harmonically oscillating microchannels. The flow exhibits highly unsteady behavior as the water and air interact with each other during the vibration cycles, making it significantly different from the well-studied steady flow in microchannels. Here, the computational fluid dynamics (CFD) modeling is realized by combing the turbulence Reynolds-averaged Navier-Stokes (RANS) k – ω model with the phase-field method to resolve the dynamics of the two-phase flow. The numerical results are qualitatively validated by the experiment. Through parametric studies, we specifically examined the effects of vibration conditions (i.e., frequency and amplitude), microchannel taper angle, and wall surface contact angle (i.e., wettability) on the flow rate through the microchannel. The results will advance the potential applications where oscillatory or general unsteady microchannel two-phase flows may be present.

42 ENGINEERING↗
DOE OSTI2020

ICME for Creep of Ni-Base Superalloys in Advanced Ultra-Supercritical Steam Turbines

The goal of this project contains three major aspects: (1) critical assessment of existing experimental data, (2) critical assessment of existing modeling capabilities and (3) development of new modeling capabilities that are critical but currently missing in predicting long-term creep behavior of Ni-base superalloys used in advanced ultra- supercritical (A-USC) steam turbines for coal-fired power plants. The project uses advanced materials informatics to datamine the experimental data, uncertainty quantification tools to rigorously assess existing creep-life models, and integration of phase-field and FFT-based crystal plasticity to develop a multi-scale, physics-based, microstructure-sensitive creep model.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2015

Strangeness driven phase transitions in compressed baryonic matter and their relevance for neutron stars and core collapsing supernovae

We discuss the thermodynamics of compressed baryonic matter with strangeness within non-relativistic mean-field models with effective interactions. The phase diagram of the full baryonic octet under strangeness equilibrium is built and discussed in connection with its relevance for core-collapse supernovae and neutron stars. A simplified framework corresponding to (n, p, Λ)(+e)-mixtures is employed in order to test the sensitivity of the existence of a phase transition on the (poorely constrained) interaction coupling constants and the compatibility between important hyperonic abundances and 2M{sub ⊙} neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗