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 775 records · Page 43

Parallel simulation via SPPARKS of on-lattice kinetic and Metropolis Monte Carlo models for materials processing

Abstract SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.

36 MATERIALS SCIENCE↗

Modeling Multiple Time Scales during Glass Formation with Phase-Field Crystals

The dynamics of glass formation in monatomic and binary liquids are studied numerically using a microscopic field theory for the evolution of the time-averaged atomic number density. A stochastic framework combining phase-field crystal free energies and dynamic density functional theory is shown to successfully describe several aspects of glass formation over multiple time scales. Agreement with mode coupling theory is demonstrated for underdamped liquids at moderate supercoolings, and a rapidly growing dynamic correlation length is found to be associated with fragile behavior.

36 MATERIALS SCIENCE↗

Phase separation versus supersolid behavior in frustrated antiferromagnets

We investigate the competition between spin supersolidity and phase separation in a frustrated spin-half model of weakly coupled dimers. We start by considering systems of hard-core bosons on the square lattice, onto which the low-energy physics of the herein investigated spin model can be mapped, and devise a criterion for gauging the interplay between supersolid order and domain-wall formation based on strong-coupling arguments. Effective bosonic models for the spin model are derived via the contractor renormalization (CORE) algorithm and we propose to combine a self-consistent cluster mean-field solution with our criterion for the occurrence of phase separation to derive the phase diagram as a function of frustration and magnetic field. In the limit of strong frustration, the model is shown to be unstable toward phase separation, in contradiction with recently published results. However, a region of stable supersolidity is identified for intermediate frustration, in a parameter range of possible experimental relevance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Finding Matrix Product State Representations of Highly Excited Eigenstates of Many-Body Localized Hamiltonians

A key property of many-body localized Hamiltonians is the area law entanglement of even highly excited eigenstates. Matrix product states (MPS) can be used to efficiently represent low entanglement (area law) wave functions in one dimension. An important application of MPS is the widely used density matrix renormalization group (DMRG) algorithm for finding ground states of one-dimensional Hamiltonians. Here, we develop two algorithms, the shift-and-invert MPS (SIMPS) and excited state DMRG which find highly excited eigenstates of many-body localized Hamiltonians. Excited state DMRG uses a modified sweeping procedure to identify eigenstates, whereas SIMPS applies the inverse of the shifted Hamiltonian to a MPS multiple times to project out the targeted eigenstate. To demonstrate the power of these methods, we verify the breakdown of the eigenstate thermalization hypothesis in the many-body localized phase of the random field Heisenberg model, show the saturation of entanglement in the many-body localized phase, and generate local excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modular approach for modeling cell motility

Modeling cell movement is a challenging task since the motility machinery is highly complex. Moreover, there is a rather broad diversity of different cell types. In order to obtain insights into generic features of the motility mechanisms of several distinct cell types, we propose a modular approach that starts with a minimal model, consisting of a phase field description of the moving cell boundary and a simplified internal dynamics. We discuss how this starting point can be extended to increase the level of detail, and how the internal dynamics "module" can be changed/adjusted to properly model various cell types. The former route allows studying specific processes involved in cell motility in the framework of a self-organized moving domain, and the latter might permit to put different cellular motility mechanisms into a unified framework.

Ziebert, F.↗

Cosmological properties of a gauged axion

We analyze the most salient cosmological features of axions in extensions of the standard model with a gauged anomalous extra U(1) symmetry. The model is built by imposing the constraint of gauge invariance in the anomalous effective action, which is extended with Wess-Zumino counterterms. These generate axionlike interactions of the axions to the gauge fields and a gauged shift symmetry. The scalar sector is assumed to acquire a nonperturbative potential after inflation, at the electroweak phase transition, which induces a mixing of the Stueckelberg field of the model with the scalars of the electroweak sector, and at the QCD phase transition. We discuss the possible mechanisms of sequential misalignments which could affect the axions of these models, and generated, in this case, at both transitions. We compute the contribution of these particles to dark matter, quantifying their relic densities as a function of the Stueckelberg mass. We also show that models with a single anomalous U(1) in general do not account for the dark energy, due to the presence of mixed U(1)-SU(3) anomalies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Running coupling and mass anomalous dimension of SU(3) gauge theory with two flavors of symmetric-representation fermions

We have measured the running coupling constant of SU(3) gauge theory coupled to N{sub f}=2 flavors of symmetric representation fermions, using the Schroedinger functional scheme. Our lattice action is defined with hypercubic smeared links which, along with the larger lattice sizes, bring us closer to the continuum limit than in our previous study. We observe that the coupling runs more slowly than predicted by asymptotic freedom, but we are unable to observe fixed point behavior before encountering a first order transition to a strong coupling phase. This indicates that the infrared fixed point found with the thin-link action is a lattice artifact. The slow running of the gauge coupling permits an accurate determination of the mass anomalous dimension for this theory, which we observe to be small, {gamma}{sub m} < or approx. 0.6, over the range of couplings we can reach. We also study the bulk and finite-temperature phase transitions in the strong coupling region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spontaneous symmetry breaking in a two-doublet lattice Higgs model

An SU(2) lattice gauge theory with two doublets of complex scalar fields is considered. All continuous symmetries are identified and, using the nonperturbative methods of lattice field theory, the phase diagram is mapped out by direct numerical simulation. Two-doublet models contain phase transitions that separate qualitatively distinct regions of the parameter space. In some regions global symmetries are spontaneously broken. For some special choices of the model parameters, the symmetry-breaking order parameter is calculated. The pattern of symmetry breaking is verified further through observation of Goldstone bosons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Strain phase separation: Formation of ferroelastic domain structures

Phase decomposition is a well-known process leading to the formation of two-phase mixtures. Here we show that a strain imposed on a ferroelastic crystal promotes the formation of mixed phases and domains, i.e., strain phase separation with local strains determined by a common tangent construction on the free energy versus strain curves. It is demonstrated that a domain structure can be understood using the concepts of domain/phase rule, lever rule, and coherent and incoherent strain phase separation, in a complete analogy to phase decomposition. The proposed strain phase separation model is validated using phase-field simulations and experimental observations of PbTiO 3 and BiFeO 3 thin films as examples. Finally, the proposed model provides a simple tool to guide and design domain structures of ferroelastic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational Examination of Orientation-Dependent Morphological Evolution during the Electrodeposition and Electrodissolution of Magnesium

Here a new model of electrodeposition and electrodissolution is developed and applied to the evolution of Mg deposits during anode cycling. The model captures Butler-Volmer kinetics, facet evolution, the spatially varying potential in the electrolyte, and the time-dependent electrolyte concentration. The model utilizes a diffuse interface approach, employing the phase field and smoothed boundary methods. Scanning electron microscope (SEM) images of magnesium deposited on a gold substrate show the formation of faceted deposits, often in the form of hexagonal prisms. Orientation-dependent reaction rate coefficients were parameterized using the experimental SEM images. Three-dimensional simulations of the growth of magnesium deposits yield deposit morphologies consistent with the experimental results. The simulations predict that the deposits become narrower and taller as the current density increases due to the depletion of the electrolyte concentration near the sides of the deposits. Increasing the distance between the deposits leads to increased depletion of the electrolyte surrounding the deposit. Two models relating the orientation-dependence of the deposition and dissolution reactions are presented. Finally, the morphology of the Mg deposit after one deposition-dissolution cycle is significantly different between the two orientation-dependence models, providing testable predictions that suggest the underlying physical mechanisms governing morphology evolution during deposition and dissolution.

36 MATERIALS SCIENCE↗

Less Common Topological Phenomena in Bulk Materials

"Prediction and subsequent discovery of topological insulators is considered to be one of the main results in condensed matter physics in the last decades. Not surprisingly, it has received major attention of both researchers and funding agencies. This attention is well-deserved; yet, one cannot but note that nearly all this research revolves about essentially the same concept: electronic excitations with linear dispersion, covering, of course, such diverse and intriguing phenomena as topological insulator, bulk Dirac states (or Weyl, if not spin degenerate), Mayorana fermions. In this project, we will address, mainly, other topological phenomena, such as topologically nontrivial magnetic patterns (as, for instance, topological Hall and related phenomena). Specifically, we propose three interrelated trusts: (1) Time-reversal symmetry breaking nonrelativistic antiferromagnets, called altermagnets. These are materials that break Kramers degeneracy of electronic bands, despite having zero net magnetization by symmetry and being fully collinear, and not necessarily non-centrosymmetric. The corresponding band structure is very similar to the band structure in non-centrosymmetric spin-orbital materials, but materials that we propose to study are distinctly different, first and foremost in the sense that despite sharing many aspects of their electronic properties with the latter, they break the time-reversal symmetry without either spin-orbit coupling or lack of inversion symmetry. (2) Topologically nontrivial magnetic spirals. The PI has been engaged with the experimental group of Dr. Ghimire at GMU investigating Dirac materials with helical magnetism, based on stacked magnetic Kagome layers, with a generic formula of RMn6Sn6. In particular, Y Mn6Sn6 demonstrates a component of the Hall effect that is naturally interpreted in terms of a topological spin texture, as well as linear magnetoresistance. Our calculation identify Dirac states that are robust with respect to the spiral formation, and let us derive an advance mean-field model explaining the observed phase diagram. This model predicted four distinct phases, with very distinct properties, which have now been seen in neutron experiments. The same compound is known to demonstrate topological Hall effect in a particular magnetic phase, and only at elevated temperature. Based on our understanding of the phase diagram, we have worked out a phenomenological theory of a chiral (skyrmionic) response to an external magnetic field, similar to the nematic response to external strain in Fe-based superconductors, which is possible in a centrosymmetric lattice and without interplanar Dzyaloshinskii-Moriya interaction. This phenomenological theory agrees quantitatively with the experiment. It is in our plans to research other similar materials for this effect. (3) Search for 3D analogues of Fe-based superconductors. We want to investigate materials that can be viewed as 3D analogues of FeSe. Specifically, we want materials that are good metals and host antiferromagnetism, which can be suppressed by pressure and generate an s-wave superconductivity, as in Fe-based superconductors. We have in mind some candidates already. This work will proceed in close collaboration with the experimental group of Prof. Nirmal Ghimire in the same department, whose expertise lies in sample making, magnetometry and transport measurements of materials with complex magnetic structures."

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of HF-CRPA CCQE Model on the Latest NOvA Results

NOvA is a long-baseline neutrino experiment based at Fermilab, dedicated to measuring various neutrino oscillation parameters. Recently, the NOvA collaboration presented new results at NEUTRINO 2024. A significant enhancement in the modeling of charged-current quasi elastic (CCQE) interactions has been achieved with the implementation of the Hartree-Fock (HF) mean-field model, incorporating continuum random phase approximation (CRPA) corrections. This HF-CRPA model offers substantial improvements in the low-energy region. In this presentation, I will discuss the impact of the HF-CRPA model on the latest three-flavor oscillation results of the NOvA experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

How to approach continuum physics in the lattice Weinberg-Salam model

We investigate the lattice Weinberg-Salam model without fermions numerically for the realistic choice of coupling constants correspondent to the value of the Weinberg angle {theta}{sub W{approx}}30 deg., and bare fine structure constant around {alpha}{approx}(1/150). We consider the values of the scalar self-coupling corresponding to Higgs mass M{sub H{approx}}100, 150, 270 GeV. It has been found that nonperturbative effects become important while approaching continuum physics within the lattice model. When the ultraviolet cutoff {Lambda}=({pi}/a) (where a is the lattice spacing) is increased and achieves the value around 1 TeV, one encounters the fluctuational region (on the phase diagram of the lattice model), where the fluctuations of the scalar field become strong. The classical Nambu monopole can be considered as an embryo of the unphysical symmetric phase within the physical phase. In the fluctuational region quantum Nambu monopoles are dense, and therefore, the use of the perturbation expansion around the trivial vacuum in this region is limited. Further increase of the cutoff is accompanied by a transition to the region of the phase diagram, where the scalar field is not condensed (this happens at the value of {Lambda} around 1.4 TeV for the considered lattice sizes). Within this region further increase of the cutoff is possible, although we do not observe this in detail due to the strong fluctuations of the gauge boson correlator. Both above mentioned regions look unphysical. Therefore we come to the conclusion that the maximal value of the cutoff admitted within lattice electroweak theory cannot exceed the value of the order of 1 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Three-dimensional simulation of unstable gravity-driven infiltration of water into a porous medium

Infiltration of water in dry porous media is subject to a powerful gravity-driven instability. Although the phenomenon of unstable infiltration is understood, its description using continuum mathematical models has posed a significant challenge for several decades. The classical model of water flow in the unsaturated flow, the Richards equation, is unable to reproduce the instability. In this work, we present a computational study of a model of unsaturated flow in porous media that extends the Richards equation and is capable of predicting the instability and captures the key features of gravity fingering quantitatively. The extended model is based on a phase-field formulation and is fourth-order in space. The new model poses a set of challenges for numerical discretizations, such as resolution of evolving interfaces, stiffness in space and time, treatment of singularly perturbed equations, and discretization of higher-order spatial partial–differential operators. We develop a numerical algorithm based on Isogeometric Analysis, a generalization of the finite element method that permits the use of globally-smooth basis functions, leading to a simple and efficient discretization of higher-order spatial operators in variational form. We highlight the accuracy, efficiency and robustness of our method with several examples in two and three dimensions in both homogeneous and strongly heterogeneous media. We simulate, for the first time, unstable gravity-driven infiltration in three dimensions, and confirm that the new theory reproduces the fundamental features of water infiltration into a porous medium. Our findings are consistent with classical experimental observations that demonstrate a transition from stable to unstable fronts depending on the infiltration flux.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An extended/generalized phase-field finite element method for crack growth with global-local enrichment

In this paper, an extended/generalized finite element method (XFEM/GFEM) for simulating quasistatic crack growth based on a phase-field method is presented. The method relies on approximations to solutions associated with two different scales: a global scale, that is, structural and discretized with a coarse mesh, and a local scale encapsulating the fractured region, that is, discretized with a fine mesh. A stable XFEM/GFEM is employed to embed the displacement and damage fields at the global scale. The proposed method accommodates approximation spaces that evolve between load steps, while preserving a fixed background mesh for the structural problem. In addition, a prediction-correction algorithm is employed to facilitate the dynamic evolution of the confined crack regions within a load step. Several numerical examples of benchmark problems in two- and three-dimensional quasistatic fracture are provided to demonstrate the approach.

multiscale↗

Investigation of LSM-YSZ Composite Cathode Performance Degradation with a Multistep Charge Transfer Model

The long-term performance degradation of LSM/YSZ composite cathodes is investigated via calibrated multiphysics simulation with a multistep oxygen reduction reaction (ORR) model and structural coarsening data from a phase field study. Multiphysics simulations, including a multistep ORR mechanism, are first developed for better understanding of the limiting processes in solid oxide fuel cells. The multistep ORR mechanism includes two parallel pathways, surface pathway and bulk pathway, which consist of elementary reaction steps. The multiphysics model is simultaneously calibrated with experimental polarization curves and impedance behavior for various air/fuel supply conditions. To our knowledge, this is the first time that a multistep ORR model is calibrated with such datasets. Next, the calibrated simulations are utilized to simulate a 2D half-cell constructed with measured microstructural data and random heterogeneity. Finally, the long-term performance degradation of the half-cell is predicted by the calibrated multiphysics model coupled with structural coarsening trends simulated using a phase field-based coarsening model. Degradation of both polarization curves and impedance behavior is investigated. Thorough analyses, including changes of contributions from different pathways, the resistance components, and overall reaction order, are performed to provide more insights into cathode performance degradation due to grain coarsening phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗