Engineering Papers⌕ Search

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.

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

At least 721 records · Page 40

Baryogenesis in the two doublet and inert singlet extension of the Standard Model

We investigate an extension of the Standard Model containing two Higgs doublets and a singlet scalar field (2HDSM). We show that the model can have a strongly first-order phase transition and give rise to the observed baryon asymmetry of the Universe, consistent with all experimental constraints. In particular, the constraints from the electron and neutron electric dipole moments are less constraining here than in pure two-Higgs-doublet model (2HDM). The two-step, first-order transition in 2HDSM, induced by the singlet field, may lead to strong supercooling and low nucleation temperatures in comparison with the critical temperature, T{sub n}≪T{sub c}, which can significantly alter the usual phase-transition pattern in 2HD models with T{sub n}≈T{sub c}. Furthermore, the singlet field can be the dark matter particle. However, in models with a strong first-order transition its abundance is typically but a thousandth of the observed dark matter abundance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

How slowly can the early Universe expand?

When the expansion of the Universe is dominated by a perfect fluid with equation of state parameter w and a sound speed c s satisfying w = $c^{2}_{s}$ ≤ 1 , the Hubble parameter H and time t satisfy the bound Ht ≥ 1/3 . There has been recent interest in “ultraslow” expansion laws with Ht <1/3 (sometimes described as “fast expanding” models). Here, we examine various models that can produce ultraslow expansion: scalar fields with negative potentials, barotropic fluids, braneworld models, or a loitering phase in the early Universe. Scalar field models and barotropic models for ultraslow expansion are unstable to evolution toward w = 1 or w → ∞ in the former case and w → ∞ in the latter case. Braneworld models can yield ultraslow expansion but require an expansion law beyond the standard Friedman equation. Loitering early universe models can produce a quasistatic expansion phase in the early Universe but require an exotic negative-density component. These results suggest that appeals to an ultraslow expansion phase in the early Universe should be approached with some caution, although the loitering early universe may be worthy of further investigation. These results do not apply to ultraslow contracting models.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding Microstructure Variability in Vapor-Deposited Energetic Materials by Using Phase-Field Methods

Critical components, such as detonators, in Sandia's stockpile contain heterogeneous materials whose performance and reliability depend on accurate, predictive models of coupled, complex phenomena to predict their synthesis, processing, and operation. Ongoing research in energetic materials has shown that microstructural properties, such as density, pore-size, morphology, and specific surface area are correlated to their initiation threshold and detonation behavior. However, experiments to study these specific characteristics of energetic materials are challenging and time consuming. Therefore, in this work, we turn to mesoscale modeling methods that may be capable of reproducing some observed phenomena to refine and predict outcomes beforehand. Even so, we have no physics-based modeling capability to predict how the microstructure of an energetic material will evolve over near- and long-term time scales. Thus, the goal of this work is to (i) identify any knowledge gaps in how the underlying microstructure forms and evolves during the synthesis process, and (ii) develop and test a mesoscale phase-field model for vapor deposition to capture critical mechanisms of microstructure formation, evolution, and variability in vapor-deposited energetic materials, such as processing conditions, material properties, and substrate interactions. Based on this work, the phase-field method is shown to be a valuable tool for developing the necessary models containing coupled, complex phenomena to investigate and understand the synthesis and processing of energetic materials.

36 MATERIALS SCIENCE↗

Exploration of the Phase Field Framework Marmot to Include Anisotropic Grain Boundaries with Molecular Dynamics [Thesis]

Using molecular dynamic (MD) calculations, meso-scale phase field simulations behave differently by including anisotropic grain boundary energy in the model for transient energy minimization in uranium-dioxide. As a preliminary step, MD will be used to find face centered cubic (FCC) copper grain boundary energies and to explore its anisotropic behavior in phase field simulations. MARMOT is a C++ code that uses object-oriented programming to implement the Allen-Cahn and Cahn-Hilliard equations for phase field modeling. MARMOT uses constants for mobility and grain boundary energy for the entire domain of the simulation. For MARMOT to reach its full capability, these constants need to vary within the domain. This paper explores the expansion of MARMOT to include an anisotropic view of grain boundary energy in transient phase field simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

General model of phospholipid bilayers in fluid phase within the single chain mean field theory

Coarse-grained model for saturated phospholipids: 1,2-didecanoyl-sn-glycero-3-phosphocholine (DCPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) and unsaturated phospholipids: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2- dioleoyl-sn-glycero-3-phosphocholine (DOPC) is introduced within the single chain mean field theory. A single set of parameters adjusted for DMPC bilayers gives an adequate description of equilibrium and mechanical properties of a range of saturated lipid molecules that differ only in length of their hydrophobic tails and unsaturated (POPC, DOPC) phospholipids which have double bonds in the tails. A double bond is modeled with a fixed angle of 120°, while the rest of the parameters are kept the same as saturated lipids. The thickness of the bilayer and its hydrophobic core, the compressibility, and the equilibrium area per lipid correspond to experimentally measured values for each lipid, changing linearly with the length of the tail. The model for unsaturated phospholipids also fetches main thermodynamical properties of the bilayers. This model is used for an accurate estimation of the free energies of the compressed or stretched bilayers in stacks or multilayers and gives reasonable estimates for free energies. The proposed model may further be used for studies of mixtures of lipids, small molecule inclusions, interactions of bilayers with embedded proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Matching time and spatial scales of rapid solidification: Dynamic TEM experiments coupled to CALPHAD-informed phase-field simulations

A combination of dynamic transmission electron microscopy (DTEM) experiments and CALPHAD-informed phase-field simulations was used to study rapid solidification in Cu–Ni thin-film alloys. Experiments—conducted in the DTEM—consisted of in situ laser melting and determination of the solidification kinetics by monitoring the solid–liquid interface and the overall microstructure evolution (time-resolved measurements) during the solidification process. Modelling of the Cu–Ni alloy microstructure evolution was based on a phase-field model that included realistic Gibbs energies and diffusion coefficients from the CALPHAD framework (thermodynamic and mobility databases). DTEM and post mortem experiments highlighted the formation of microsegregation-free columnar grains with interface velocities varying from ~0.1 to ~0.6 m s –1 . After an 'incubation' time, the velocity of the planar solid–liquid interface accelerated until solidification was complete. In addition, a decrease of the temperature gradient induced a decrease in the interface velocity. The modelling strategy permitted the simulation (in 1D and 2D) of the solidification process from the initially diffusion-controlled to the nearly partitionless regimes. Lastly, results of DTEM experiments and phase-field simulations (grain morphology, solute distribution, and solid–liquid interface velocity) were consistent at similar time (μs) and spatial scales (μm).

36 MATERIALS SCIENCE↗

Development of a Grain Boundary Pinning Model that Considers Particle Size Distribution Using the Phase Field Method

Grain boundary (GB) migration significantly impacts material behavior. However, GB migration is slowed or even halted by resistive pressure applied by pores or particles. Zener’s original investigation of particle pinning, and subsequent modifications by other researchers, describe the resistive pressure for various spatial distributions of particles with respect to GBs. In this work, we develop a pinning model that considers the impact of the particle size distribution and we verify it by comparing to mesoscale phase field and Monte Carlo simulations. Resistive pressure expressions are developed that are functions of the percentage of GB area covered by particles and of the particle volume fraction for any spatial distribution of particles. In both expressions, the mean value of the resistive pressure decreases with increasing standard deviation of the particle radius.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Calibration of Multi-Parameter Models of Avascular Tumor Growth Using Time Resolved Microscopy Data

Abstract Two of the central challenges of using mathematical models for predicting the spatiotemporal development of tumors is the lack of appropriate data to calibrate the parameters of the model, and quantitative characterization of the uncertainties in both the experimental data and the modeling process itself. We present a sequence of experiments, with increasing complexity, designed to systematically calibrate the rates of apoptosis, proliferation, and necrosis, as well as mobility, within a phase-field tumor growth model. The in vitro experiments characterize the proliferation and death of human liver carcinoma cells under different initial cell concentrations, nutrient availabilities, and treatment conditions. A Bayesian framework is employed to quantify the uncertainties in model parameters. The average difference between the calibration and the data, across all time points is between 11.54% and 14.04% for the apoptosis experiments, 7.33% and 23.30% for the proliferation experiments, and 8.12% and 31.55% for the necrosis experiments. The results indicate the proposed experiment-computational approach is generalizable and appropriate for step-by-step calibration of multi-parameter models, yielding accurate estimations of model parameters related to rates of proliferation, apoptosis, and necrosis.

60 APPLIED LIFE SCIENCES↗

Modelling Accident Tolerant Fuel Concepts

The catastrophic events that occurred at the Fukushima-Daiichi nuclear power plant in 2011 have led to widespread interest in research of alternative fuels and claddings that are proposed to be accident tolerant. The United States Department of Energy (DOE) through its Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded an Accident Tolerant Fuel (ATF) High Impact Problem (HIP). The ATF HIP is a three-year project to perform research on two accident tolerant concepts. The final outcome of the ATF HIP will be an in-depth report to the DOE Advanced Fuels Campaign (AFC) giving a recommendation on whether either of the two concepts should be included in their lead test assembly scheduled for placement into a commercial reactor in 2022. The two ATF concepts under investigation in the HIP are uranium silicide fuel and iron-chromium-aluminum (FeCrAl) alloy cladding. Utilizing the expertise of three national laboratory participants (Idaho National Laboratory, Los Alamos National Laboratory, and Argonne National Laboratory), a comprehensive multiscale approach to modeling is being used that includes atomistic modeling, molecular dynamics, rate theory, phase-field, and fuel performance simulations. Model development and fuel performance analysis are critical since a full suite of experimental studies will not be complete before AFC must prioritize concepts for focused development. In this paper, we present simulations of the two proposed accident tolerance fuel systems: U3Si2 fuel with Zircaloy-4 cladding, and UO2 fuel with FeCrAl cladding. Sensitivity analyses are completed using Sandia National Laboratories’ Dakota software to determine which input parameters (e.g., fuel specific heat) have the greatest influence on the output metrics of interest (e.g., fuel centerline temperature). We also outline the multiscale modelling approach being employed. Considerable additional work is required prior to preparing the recommendation report for the Advanced Fuels Campaign.

Accident Tolerant Fuel↗

Quantum phase transition dynamics in the two-dimensional transverse-field Ising model

The quantum Kibble-Zurek mechanism (QKZM) predicts universal dynamical behavior near the quantum phase transitions (QPTs). It is now well understood for the one-dimensional quantum matter. Higher-dimensional systems, however, remain a challenge, complicated by the fundamentally different character of the associated QPTs and their underlying conformal field theories. In this work, we take the first steps toward theoretical exploration of the QKZM in two dimensions for interacting quantum matter. We study the dynamical crossing of the QPT in the paradigmatic Ising model by a joint effort of modern state-of-the-art numerical methods, including artificial neural networks and tensor networks. As a central result, we quantify universal QKZM behavior close to the QPT. We also note that, upon traversing further into the ferromagnetic regime, deviations from the QKZM prediction appear. We explain the observed behavior by proposing an extended QKZM taking into account spectral information as well as phase ordering. Our work provides a testing platform for higher-dimensional quantum simulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Domain patterns and super-elasticity of freestanding BiFeO 3 membranes via phase-field simulations

Super-elasticity of functional ferroelectric oxides offers promises for integrating ferroelectric films into flexible electronics. However, super-elastic deformation is a complex phenomenon related to possibly multiple concurrent mechanisms. Fundamentally understanding how multiple mechanisms contribute to the super-elasticity of ferroelectric oxides is crucial to realizing their potential flexible electronic applications. In this study, we employ phase-field simulations to model the dynamics of ferroelectric domain patterns of freestanding BiFeO 3 membranes to understand the origin of their super-elasticity under substantial bending deformation (5% strain). It is demonstrated that both a reversible Rhombohedral-Tetragonal (R-T) phase transition and a nearly reversible domain evolution of BiFeO 3 membranes contribute to accommodating the large deformation and thus their super-elasticity. The dynamics of domain evolution also reveal the formation of an exotic ferroelectric vortex and polarization rotation before the phase transition. We constructed a diagram of phases and domain patterns as a function of the membrane thickness and bending angle, which allows one to readily predict the emergence of T phase and ferroelectric vortex in bent BFO membranes. These results not only provide fundamental understanding of mesoscale super-elastic mechanisms but also reveal exotic domain states of ferroelectric membranes.

36 MATERIALS SCIENCE↗

Helical Phase Inflation and Monodromy in Supergravity Theory

We study helical phase inflation which realizes “monodromy inflation” in supergravity theory. In the model, inflation is driven by the phase component of a complex field whose potential possesses helicoid structure. We construct phase monodromy based on explicitly breaking global U ( 1 ) symmetry in the superpotential. By integrating out heavy fields, the phase monodromy from single complex scalar field is realized and the model fulfills natural inflation. The phase-axion alignment is achieved from explicitly symmetry breaking and gives super-Planckian phase decay constant. The F -term scalar potential provides strong field stabilization for all the scalars except inflaton, which is protected by the approximate global U ( 1 ) symmetry. Besides, we show that helical phase inflation can be naturally realized in no-scale supergravity with S U ( 2,1 ) / S U ( 2 ) × U ( 1 ) symmetry since the supergravity setup needed for phase monodromy is automatically provided in the no-scale Kähler potential. We also demonstrate that helical phase inflation can be reduced to another well-known supergravity inflation model with shift symmetry. Helical phase inflation is free from the UV-sensitivity problem although there is super-Planckian field excursion, and it suggests that inflation can be effectively studied based on supersymmetric field theory while a UV-completed framework is not prerequisite.

Li, Tianjun↗

Dendrite-resolved, full-melt-pool phase-field simulations to reveal non-steady-state effects and to test an approximate model

In this report we study the epitaxial, columnar growth of (multiply oriented) dendrites/cells for a spot melt in a polycrystalline Al–Cu substrate using two-dimensional, phase-field, direct numerical simulations (DNS) at the full-melt-pool scale. Our main objective is to compare the expensive DNS model to a much cheaper but approximate “line” model in which a single-crystal phase-field simulation is confined to a narrow rectangular geometry. To perform this comparison, we develop algorithms that automatically extract quantities of interest (QoIs) from both DNS and line models. These QoIs allow us to quantitatively assess the assumptions in the line model and help us analyze its discrepancy with the DNS model. We consider four sets of heat source parameters, mimicking welding and additive manufacturing conditions, that create a combination shallow and deep melt pools. Our largest DNS simulation used 16K x 14K grid points in space. Our main findings can be summarized as follows. Under AM conditions, the QoIs of line models are in excellent agreement with the full DNS results for both shallow and deep melt pools. Under welding conditions, the primary spacing of the DNS model is smaller than the prediction of line model. We identify a geometric crowding effect that accounts for the discrepancies between the DNS and line models. We propose two potential mechanisms that determine the response of the microstructure to geometric crowding.

36 MATERIALS SCIENCE↗

Roadmap on multiscale materials modeling

Modeling and simulation is transforming modern materials science, becoming an important tool for the discovery of new materials and material phenomena, for gaining insight into the processes that govern materials behavior, and, increasingly, for quantitative predictions that can be used as part of a design tool in full partnership with experimental synthesis and characterization. Modeling and simulation is the essential bridge from good science to good engineering, spanning from fundamental understanding of materials behavior to deliberate design of new materials technologies leveraging new properties and processes. This Roadmap presents a broad overview of the extensive impact computational modeling has had in materials science in the past few decades, and offers focused perspectives on where the path forward lies as this rapidly expanding field evolves to meet the challenges of the next few decades. The Roadmap offers perspectives on advances within disciplines as diverse as phase field methods to model mesoscale behavior and molecular dynamics methods to deduce the fundamental atomic-scale dynamical processes governing materials response, to the challenges involved in the interdisciplinary research that tackles complex materials problems where the governing phenomena span different scales of materials behavior requiring multiscale approaches. The shift from understanding fundamental materials behavior to development of quantitative approaches to explain and predict experimental observations requires advances in the methods and practice in simulations for reproducibility and reliability, and interacting with a computational ecosystem that integrates new theory development, innovative applications, and an increasingly integrated software and computational infrastructure that takes advantage of the increasingly powerful computational methods and computing hardware.

36 MATERIALS SCIENCE↗

Nanoscale Origins of Ferroelastic Domain Wall Mobility in Ferroelectric Multilayers

Here we investigate the nanoscale origins of ferroelastic domain wall motion in ferroelectric multilayer thin films that lead to giant electromechanical responses. We present direct evidence for complex underpinning factors that result in ferroelastic domain wall mobility using a combination of atomic-level aberration corrected scanning transmission electron microscopy and phase-field simulations in model epitaxial (001) tetragonal (T) PbZr x Ti 1-x O 3 (PZT)/rhombohedral (R) PbZr x Ti 1-x O 3 (PZT) bilayer heterostructures. The local electric dipole distribution is imaged on an atomic scale for a ferroelastic domain wall that nucleates in the R-layer and cuts through the composition breaking the T/R interface. Our studies reveal a highly complex polarization rotation domain structure that is nearly on the knife-edge at the vicinity of this wall. Induced phases, namely tetragonal-like and rhombohedral-like monoclinic were observed close to the interface, and exotic domain arrangements, such as a half-four-fold closure structure, are observed. Phase field simulations show this is due to the minimization of the excessive elastic and electrostatic energies driven by the enormous strain gradient present at the location of the ferroelastic domain walls. Thus, in response to an applied stimulus, such as an electric field, any polarization reorientation must minimize the elastic and electrostatic discontinuities due to this strain gradient, which would induce a dramatic rearrangement of the domain structure. This insight into the origins of ferroelastic domain wall motion will allow researchers to better “craft” such multilayered ferroelectric systems with precisely tailored domain wall functionality and enhanced sensitivity, which can be exploited for the next generation of integrated piezoelectric technologies.

36 MATERIALS SCIENCE↗

Seasonal and Interannual Patterns and Controls of Hydrological Fluxes in an Amazon Floodplain Lake With a Surface-Subsurface Process Model

Floodplain lakes represent important aquatic ecosystems, and field-based estimates of their water budgets are difficult to obtain, especially over multiple years. We examine the hydrological fluxes for an Amazon floodplain lake connected to the Solimões River using a process-based hydrologic model. Water exchanges between the river and lake agree well with field estimates, including the timing of different hydrological phases. However, beyond available field data, modeling results show that the seven simulated years all differed from each other. These interannual differences were caused by the interplay between phases when water levels were rising with river-water flowing into the lake (RWRI), versus rising with lake-water flowing out to the river (RWLO). This exchange determines the river-water content in the lake ( C L ). Maximum C L occurred before river levels peaked because local catchment contributions can be sufficient to push lake-water out to the river, even as river levels rise. Numerical experiments show that the seasonal distribution of local rainfall, local catchment size, and interannual variability in both climate and river stage can contribute to differing dynamics of C L in a floodplain lake. Their impacts vary among phases: river-rise dominates the RWRI, whereas local hydrological processes dominate the RWLO and receding-water phases. Intermediate-to-long-term rainfall accumulation controls C L during the RWLO phase, whereas annual precipitation accumulation is important for C L during low water. Finally, our model generalizes beyond limited available field studies and offers potential to better understand floodplain lakes in other areas and how regional versus local changes in climate may affect their hydrological dynamics.

54 ENVIRONMENTAL SCIENCES↗