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 1,045 records · Page 58

Nonlocal PNJL model beyond mean field

A nonlocal chiral quark model is consistently extended beyond mean field using a strict 1/N{sub c} expansion scheme. It is found that the 1/N{sub c} corrections lead to a lowering of the temperature of the chiral phase transition in comparison with the mean-field result. On the other hand, near the phase transition the 1/N{sub c} expansion breaks down and a nonperturbative scheme for the inclusion of mesonic correlations is needed in order to describe the phase transition region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Chiral condensate and dressed Polyakov loop in the Nambu-Jona-Lasinio model

We investigate the chiral condensate and the dressed Polyakov loop or dual chiral condensate at finite temperature and density in the two-flavor Nambu-Jona-Lasinio model. The dressed Polyakov loop is regarded as an equivalent order parameter of deconfinement phase transition in a confining theory. We find the behavior of a dressed Polyakov loop in the absence of any confinement mechanism quite interesting. With only quark degrees of freedom present, it still shows an order parameter like behavior. It is found that in the chiral limit, the critical temperature for chiral phase transition coincides with that of the dressed Polyakov loop in the whole (T,{mu}) plane. In the case of explicit chiral symmetry breaking, it is found that the transition temperature for chiral restoration T{sub c}{sup {chi}} is smaller than that of the dressed Polyakov loop T{sub c}{sup D} in the low baryon density region where the transition is a crossover. With the increase of the current quark mass the difference between the two transition temperatures is found to be increasing. However, the two critical temperatures coincide in the high baryon density region where the phase transition is of first order. We give an explanation on the feature of T{sub c}{sup {chi}}=T{sub c}{sup D} in the case of 1st and 2nd order phase transitions, and T{sub c}{sup {chi}}<T{sub c}{sup D} in the case of crossover, and expect this feature is general and can be extended to full QCD theory. Our result might indicate that in the case of crossover, there exists a small region where chiral symmetry is restored, but the color degrees of freedom are still confined.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Particle ratios as a probe of the QCD critical temperature

We show how the measured particle ratios can be used to provide nontrivial information about the critical temperature of the QCD phase transition. This is obtained by including the effects of highly massive Hagedorn resonances on statistical models, which are used to describe hadronic yields. The inclusion of Hagedorn states creates a dependence of the thermal fits on the Hagedorn temperature, T{sub H}, which is assumed to be equal to T{sub c}, and leads to an overall improvement of thermal fits. We find that for Au+Au collisions at the Relativistic Heavy Ion Collider at {radical}(s{sub NN})=200 GeV the best square fit measure, {chi}{sup 2}, occurs at T{sub c{approx}}176 MeV and produces a chemical freeze-out temperature of 172.6 MeV and a baryon chemical potential of 39.7 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spatially resolved polarization swings in the supermassive binary black hole candidate OJ 287 with first Event Horizon Telescope observations

We present the first Event Horizon Telescope 1.3 mm observations of the supermassive binary black hole candidate OJ 287. The observations achieved an unprecedented angular resolution of 18 μas and reveal significant structural and polarization variability over just five days, marking the shortest timescale on which such changes have been directly imaged in this source. The inner jet exhibits a twisted ridgeline structure, with features displaying apparent superluminal motions up to about 22 c. The linear polarization maps reveal three main polarized features whose electric-vector position angles (EVPAs) change substantially over the time span of our observations, including a component with a radial polarization consistent with being produced by a recollimation shock. Most notably, we directly resolved two innermost jet components whose EVPAs rotate in opposite directions. The faster component, moving at 2.4 ± 0.9 μas/day (17.4 ± 6.5 c), exhibits counterclockwise EVPA swings of roughly 3.7° per day, while the slower component, with a proper motion of 1.4 ± 0.3 μas/day (10.2 ± 2.2 c), rotates clockwise at approximately 2.5° per day. Previous studies inferred helical magnetic fields in AGN jets from time-resolved or integrated polarization variability but lacked the angular resolution to directly image this effect. Our results provide spatially resolved evidence that a helical magnetic field threads the jet’s collimation and acceleration zone, ruling out models based on the superposition of unresolved components. Our analysis suggests that propagating shocks interact with a Kelvin–Helmholtz plasma instability, illuminating different phases of the helical magnetic field and producing the observed polarization spatial and temporal variability. Moreover, our model naturally accounts for the more rapid polarization rotation observed in the faster moving component. Our model predicts even more rapid swings in polarization, which could be tested with future observations featuring a more densely sampled time coverage.

OJ 287↗

Analysis and modeling of the equilibrium behaviors of U and Pu in molten LiCl-KCl/Cd system at 500 °C

In the field of pyroprocessing of used nuclear fuels, liquid cadmium is applied in different ways to recover transuranics from molten salt, while effecting some degree of separations between transuranics and lanthanides. In this work, a model of the separation performance between uranium and plutonium across the phase boundary between molten salt and liquid cadmium is presented. The salt is the LiCl-KCl eutectic with solute concentrations of UCl 3 and PuCl 3 . The cadmium is the liquid cadmium phase field from the Cd-U-Pu ternary alloy system. Lastly, the model used to describe the behavior of uranium and plutonium in the salt/cadmium system is applied to describe the behavior of plutonium and lanthanide in the salt/cadmium system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Microemulsions in the driven Widom-Rowlinson lattice gas

An investigation of the two-dimensional Widom-Rowlinson lattice gas under an applied drive uncovered a remarkable nonequilibrium steady state in which uniform stripes (reminiscent of an equilibrium lamellar phase) form perpendicular to the drive direction [R. Dickman and R. K. P. Zia, Phys. Rev. E 97, 062126 (2018)]. Here we study this model at low particle densities in two and three dimensions, where we find a disordered phase with a characteristic length scale (a “microemulsion”) along the drive direction. Further we develop a continuum theory of this disordered phase to derive a coarse-grained field-theoretic action for the nonequilibrium dynamics. The action has the form of two coupled driven diffusive systems with different characteristic velocities, generated by an interplay between the particle repulsion and the drive. We then show how fluctuation corrections in the field theory may generate the characteristic features of the microemulsion phase, including a peak in the static structure factor corresponding to the characteristic length scale. This work lays the foundation for understanding the stripe phenomenon more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Coverage Self-Consistent Microkinetic Model for Vapor-Phase Formic Acid Decomposition over Pd/C Catalysts

An iterative approach utilizing density functional theory (DFT, PW91-GGA)-informed mean-field microkinetic models and reaction kinetics experiments is used to determine the reaction mechanism and the active site for formic acid (HCOOH, FA) decomposition over a Pd/C catalyst. Models parametrized using DFT energetics on clean Pd(100) and Pd(111) required large corrections to the DFT energetics for capturing our experimental data. Further, both Pd(111) and Pd(100) models predicted a high coverage of adsorbed CO (CO*), inconsistent with the assumption of a clean surface at which the rate parameters for these models were calculated. To better represent the active site under reaction conditions and explicitly account for the presence of CO*, subsequent microkinetic models were formulated using DFT energetics that were calculated on partially (5/9 ML) CO*-covered Pd (111) and (100) facets. Upon parameter adjustment, the resultant 5/9 ML CO*-covered Pd(100) model, although consistent in terms of CO* coverage, was unable to capture the dehydration path measured in the experiments and was, therefore, deemed not to offer an accurate representation of the active site for FA decomposition over Pd/C. In contrast, a partially CO*-covered Pd(111) model was better at representing the catalytic active site, as in addition to being consistent in terms of CO* coverages, it required small adjustments of the DFT parameters to accurately capture the experimental data set (both dehydrogenation and dehydration). Our results suggest that the reaction occurs via the spectroscopically elusive carboxyl (COOH*) intermediate and that spectator CO*-assisted decomposition pathways play an important role under typical experimental conditions. In addition, our study highlights the importance of striving for coverage self-consistent microkinetic models and for including spectator-assisted mechanisms in order to develop an improved picture of the active site under reaction conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adaptive variational quantum minimally entangled typical thermal states for finite temperature simulations

Scalable quantum algorithms for the simulation of quantum many-body systems in thermal equilibrium are important for predicting properties of quantum matter at finite temperatures. Here we describe and benchmark a quantum computing version of the minimally entangled typical thermal states (METTS) algorithm for which we adopt an adaptive variational approach to perform the required quantum imaginary time evolution. The algorithm, which we name AVQMETTS, dynamically generates compact and problem-specific quantum circuits, which are suitable for noisy intermediate-scale quantum (NISQ) hardware. We benchmark AVQMETTS on statevector simulators and perform thermal energy calculations of integrable and nonintegrable quantum spin models in one and two dimensions and demonstrate an approximately linear system-size scaling of the circuit complexity. We further map out the finite-temperature phase transition line of the two-dimensional transverse field Ising model. Finally, we study the impact of noise on AVQMETTS calculations using a phenomenological noise model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chiral restoration at finite temperature with meson loop corrections

We investigate the pattern of chiral-symmetry restoration of QCD for N{sub c}=3 and N{sub f}=2 at finite temperature (T) beyond the chiral limit. To this end, we employ the instanton-vacuum configuration for the flavor SU(2) sector and the Harrington-Shepard caloron for modifying relevant instanton parameters as functions of T. The meson loop corrections (MLC), which correspond to 1/N{sub c} corrections, are also taken into account to reproduce appropriate m{sub q} dependences of chiral order parameters. We compute the chiral condensate as a function of T and/or m{sub q}. We observe that MLC play an important role to have a correct universality-class behavior of chiral-restoration patterns in this framework, depending on m{sub q}: Second-order phase transition in the chiral limit m{sub q}=0 and cross-over for m{sub q{ne}}0. Without MLC, all the restoration patterns are crossover, due to simple saddle-point approximations. It turns out that T{sub c}{sup {chi}=}159 MeV in the chiral limit and T{sub c}{sup {chi}=}(177,186,196) MeV for m{sub q}=(5,10,15) MeV, using the phenomenological choices for the instanton parameters at T=0.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Phase Equilibria in the ZnO-“FeO”-SiO{sub 2} System in Reducing Atmosphere and in the ZnO-“FeO”-SiO{sub 2}-“Cu{sub 2}O” System in Equilibrium with Liquid Copper Metal at 1250 °C (1523 K)

Recent experimental studies in the ZnO-“FeO”-SiO{sub 2} system in reducing atmosphere demonstrated significant discrepancies with the current FactSage thermodynamic model developed using previous experimental data in this system in equilibrium with metallic iron and air. The present experimental study on phase equilibria in the ZnO-“FeO”-SiO{sub 2}-“Cu{sub 2}O” system in equilibrium with liquid copper at 1250 °C (1523 K) at low copper oxide concentrations in slag was initiated and undertaken to resolve these discrepancies. A high-temperature equilibration–rapid quenching–electron-probe X-ray microanalysis (EPMA) technique using a primary phase substrate support and closed system approach with Cu metal introduced to determine effective equilibrium oxygen partial pressure from the Cu{sub metal}/Cu{sub 2}O{sub slag} equilibria was applied to provide accurate information on the liquidus and corresponding solid compositions in the spinel, willemite, and tridymite primary phase fields. The present results confirmed the accuracy of the FactSage model, resolved discrepancies, and demonstrated significant uncertainties in the recent studies by other authors on the system in the open reducing atmosphere. The present study shows how this closed system approach can be used to obtain key thermodynamic data on phase equilibria in systems containing volatile metal species, overcoming the limitations and uncertainties encountered in conventional open gas/condensed phase equilibration with these systems. The study highlights the importance of the focus on obtaining accurate experimental data and the risks of misleading information from inadequate experimental control and analysis. The study also demonstrates that continuing in-depth critical review and analysis of the elemental reactions taking place in complex systems is an essential step in phase equilibrium research.

36 MATERIALS SCIENCE↗

Superradiant phase in field-matter interactions

We show that semiclassical states adapted to the symmetry of the Hamiltonian are an excellent approximation to the exact quantum solution of the ground and first excited states of the Dicke model. Their overlap with the exact quantum states is very close to 1 except in a close vicinity of the quantum phase transition. Furthermore, they have analytic forms in terms of the model parameters and allow us to calculate analytically the expectation values of field and matter observables. Some of these differ considerably from results obtained via the standard coherent states and by means of Holstein-Primakoff series expansion of the Dicke Hamiltonian. Comparison with exact solutions obtained numerically supports our results. In particular, it is shown that the expectation values of the number of photons and of the number of excited atoms have no singularities at the phase transition. We comment on why other authors have previously found otherwise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Influence of pions on the hadron-quark phase transition

In this work we present the features of the hadron-quark phase transition diagrams in which the pions are included in the system. To construct such diagrams we use two different models in the description of the hadronic and quark sectors. At the quark level, we consider two distinct parametrizations of the Polyakov-Nambu-Jona-Lasinio (PNJL) models. In the hadronic side, we use a well known relativistic mean-field (RMF) nonlinear Walecka model. We show that the effect of the pions on the hadron-quark phase diagrams is to move the critical end point (CEP) of the transitions lines. Such an effect also depends on the value of the critical temperature (T{sub 0}) in the pure gauge sector used to parametrize the PNJL models. Here we treat the phase transitions using two values for T{sub 0}, namely, T{sub 0}= 270 MeV and T{sub 0}= 190 MeV. The last value is used to reproduce lattice QCD data for the transition temperature at zero chemical potential.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗

Toward High Fidelity LES Simulations of 5 x 5 PWR Fuel Bundle with Mixing Vane

The present work illustrates one example of a computational fluid dynamics (CFD) application in the nuclear engineering field. Since CFD is one of the most powerful tools available at the actual state of the art that could in principle help and assist the design process, a compromise should be established between the accuracy of the CFD numerical model (RANS - LES - DNS) and the computational cost of the simulations. In this spirit, this work will present the computational results from two different CFD model categories, Large Eddy Simulations (LES) and Reynolds-averaged Navier-Stokes (RANS) models. Both treat the flow turbulence issue in two separate ways, which differ mostly in the resolution of the turbulence statistics and the computational cost. The single-phase hydraulic problem presented involves the steady-state turbulent flow field modelling and resolution around a scaled PWR mixing vane grid, designed by Westinghouse and tested by Texas A and M University. The role of the mixing vane/spacer grid is to enhance the turbulent mixing and hence heat transfer of the core coolant around fuel bundles, thus increasing the departure of nucleate boiling-critical heat flux (DNB-CHF) point, as well as important structural functions for the rods bundles. The authors numerically studied the flow field characteristics at the region downstream of the spacer grids and validated their results to particle imagery velocimetry (PIV) measurements. Dominguez and Hassan performed two-dimensional time resolved PIV (TR-PIV) measurements within a 5 x 5 PWR rod bundle with spacer-grids that utilized a matching index of refraction (MIR) approach. Their averaged data showed a rather uniform direction velocity profile after the spacer grid while the downstream illustrated a no slip wall profile. The available PIV flow field experimental results from Texas A and M show complex structures of the flow field downstream the mixing vanes with the formation of secondary flows that tend to enhance the turbulent heat transfer, due to the mixing vanes geometry itself. Computational results will be validated with the available experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effects of natural zeolites on field-scale geologic noble gas transport

Improving predictive models for noble gas transport through natural materials at the field-scale is an essential component of improving US nuclear monitoring capabilities. Several field-scale experiments with a gas transport component have been conducted at the Nevada National Security Site (Non-Proliferation Experiment, Underground Nuclear Explosion Signatures Experiment). However, the models associated with these experiments have not treated zeolite minerals as gas adsorbing phases. This is significant as zeolites are a common alteration mineral with a high abundance at these field sites and are shown here to significantly fractionate noble gases during field-scale transport. This fractionation and associated retardation can complicate gas transport predictions by reducing the signal-to-noise ratio to the detector (e.g. mass spectrometers or radiation detectors) enough to mask the signal or make the data difficult to interpret. Omitting adsorption-related retardation data of noble gases in predictive gas transport models therefore results in systematic errors in model predictions where zeolites are present.Herein is presented noble gas adsorption data collected on zeolitized and non-zeolitized tuff. Experimental results were obtained using a unique piezometric adsorption system designed and built for this study. Data collected were then related to pure-phase mineral analyses conducted on clinoptilolite, mordenite, and quartz. These results quantify the adsorption capacity of materials present in field-scale systems, enabling the modeling of low-permeability rocks as significant sorption reservoirs vital to bulk transport predictions.

58 GEOSCIENCES↗

Topological and magnetic phase transitions in the bilayer Kitaev-Ising model

We investigate the phase diagram of a bilayer Kitaev honeycomb model with Ising interlayer interactions, deriving effective models via perturbation theory and performing Majorana mean-field theory calculations. We show that a diverse array of magnetic and topological phase transitions occur, depending on the direction of the interlayer Ising interaction and the relative sign of Kitaev interactions. When two layers have the same sign of the Kitaev interaction, a first-order transition from a Kitaev spin liquid to a magnetically ordered state takes place. The magnetic order points along the Ising axis and it is (anti)ferromagnetic for (anti)ferromagnetic Kitaev interactions. However, when two layers have opposite signs of the Kitaev interaction, we observe a notable weakening of magnetic ordering tendencies and the Kitaev spin liquid survives up to a remarkably larger interlayer exchange. Our mean-field analysis suggests the emergence of an intermediate gapped ℤ 2 spin-liquid state, which eventually becomes unstable upon vison condensation. The confined phase is described by a highly frustrated 120° compass model. We furthermore use perturbation theory to study the model with the Ising axis pointing along the $\hat{z}$ axis or lying in the 𝑥⁢𝑦 plane. In both cases, our analysis reveals the formation of one-dimensional Ising chains, which remain decoupled in perturbation theory, resulting in a subextensive ground-state degeneracy. Our results highlight the interplay between topological order and magnetic ordering tendencies in bilayer quantum spin liquids.

Kitaev model↗

Dynamics of electrified liquid metal surface using shallow water model

A shallow water model that incorporates surface tension and electric field effects is developed to investigate the dynamics of an electrified liquid surface. The computational model is verified against the Zakharov–Kuznetsov equation and is applied to study the growth and damping of the electrified liquid surface. A linear wave analysis is performed under a shallow water theory assuming an analytic solution of the electric field, similar to the Tonks–Frenkel instability. The electrified liquid surface grows or dampens based on the balance of the electric field, surface tension, and gravitational forces. As a result, the numerical results obtained from the electrified shallow water solver are in good agreement with the theoretical analysis.

Computational models↗