SEARCH · Engineering Papers
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.
We discuss phases of hot and dense hadronic matter using chiral Lagrangians. A two-flavored parity doublet model constrained by the nuclear matter ground state predicts chiral symmetry restoration. The model thermodynamics is shown within the mean-field approximation. A field-theoretical constraint on possible phases from the anomaly matching is also discussed.
We describe the critical temperature dependence on the pion mass and on the isospin chemical potential in remarkably good agreement with lattice data. Our framework incorporates explicit dependence on quark masses, isospin and baryonic chemical potentials for the case of two flavors through ingredients from well-known high-and low-energy theories.
Grizzly is a multi-scale multi-physics tool being developed at Idaho National Laboratory (INL) as part of the US Department of Energy's Light Water Reactor Sustainability program to provide improved safety assessments of systems, components, and structures in nuclear power plants subjected to age-related degradation. Its goal is to provide an improved scientific basis for decisions surrounding license renewal, which would permit operation of commercial nuclear power plants beyond 60 years. Grizzly is based on INL's MOOSE framework, which enables multi-physics simulations in a massively parallel computing environment. The reactor pressure vessel (RPV) has been chosen as the initial application for Grizzly. The primary concern for the RPV is that it becomes increasingly susceptible to fracture initiation at the sites of pre-existing flaws over time. This is because exposure to irradiation and high temperatures embrittles the steel. Assessing the safety of an RPV under transient conditions involves simulations at multiple scales. At the engineering scale, strongly coupled equations of heat conduction and solid mechanics are solved to simulate the global response of the RPV to accident conditions. Sub-models are used to represent regions with pre-existing flaws and the stress intensity factors are evaluated from the domain integrals. Crack propagation is then determined using a model that provides the temperature and neutron fluence dependent fracture toughness of the RPV steel. The ability of Grizzly to provide accurate engineering scale estimates of RPV safety beyond 60 years of operation strongly relies on the predictability of the embrittlement model. Multiple models for embrittlement have been proposed, and the physically based empirical model is implemented in Grizzly. However, the formation of late blooming phases greatly reduces the confidence on the extrapolated values obtained from embrittlement model, which is further exacerbated by the absence of real field data for validation. Experimentally validated lower length scale models can significantly improve the predictability of the embrittlement model beyond the calibrated range and the development of such models is another ongoing effort in Grizzly. The modeling activity in Grizzly spans across different length scales, and can be divided into: (i) Irradiation induced microstructure evolution at the atomistic and single crystal length scale using molecular dynamics, atomistic kinetic Monte-Carlo, phase-field and rate theory; (ii) Effect of the irradiation induced defects on the flow stress and fracture processes at the polycrystalline scale using crystal plasticity. Subsequently, the prediction from these models can be used to obtain irradiation induced ductile-to-brittle transition temperature (DBTT) shifts from specimen-scale simulations; (iii) Global thermo-mechanical response of the RPV under normal and accident scenarios, and probabilistic risk analyses at the engineering scale using the DBTT shifts obtained from the lower length scale models. The item (i) has been discussed in a previous report. In this extended abstract an overview of items (ii) and (iii) is provided. The extended abstract summarizes the current status of the modeling effort under Grizzly. At the engineering scale, reduced order model development is underway using the statistical tool RAVEN. This will allow probabilistic analysis considering a wide range of flaw topologies in a computationally efficient manner. Coupling the irradiation induced defect growth predictions to crystal plasticity is also in progress. (authors)
Scalar spin chirality is expected to induce a finite contribution to the Hall response at low temperatures. We observe this spin-chirality-driven Hall effect, known as the topological Hall effect, at the manganite side of the interface between La 1- x Sr x MnO 3 and SrIrO 3 . The ferromagnetic double-exchange hopping at the manganite layer, in conjunction with the Dzyaloshinskii-Moriya (DM) interaction which arises at the interface due to broken inversion symmetry and strong spin-orbit coupling from the iridate layer, could produce a skyrmion-crystal (SkX) phase in the presence of an external magnetic field. Using the Monte Carlo technique and a two-orbital spin-fermion model for manganites, supplemented by an in-plane DM interaction, we obtain phase diagrams which reveal at low temperatures a clear SkX phase and also a low-field spin-spiral phase. Increasing temperature, a skyrmion-gas phase, precursor of the SkX phase upon cooling, was identified. The topological Hall effect primarily appears in the SkX phase, as observed before in oxide heterostructures. We determine that the manganite-iridate superlattices give another useful platform to discover a plethora of unconventional magnetic and transport properties.
This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.
High–dimensional datasets are becoming ubiquitous in many applications and therefore unsupervised tensor methods to interrogate them are needed. Here, we report a new unsupervised machine learning (ML) approach (NTFk) based on nonnegative tensor factorization integrated with a custom k–means clustering. We demonstrate the ability of NTFk to extracting temporal and spatial features of phase separation of copolymers as they are modeled by self–consistent field theory. Microphase separation of block copolymers has been extensively studied both experimentally and theoretically. However, the interpretation of computer simulations and/or experimental data, representing temporal and spatial changes of molecular species concentration is still a challenging task. Thus, extracting the phase diagram from simulations or experimental data as well as the interpretation of data requires discernment of the model/experimental parameters (such as, temperature, concentrations, the number of molecular species and the interaction between species) impact on the microphase separation process. An attractive and unique aspect of the introduced ML method is that it ensures the nonnegativity of the extracted latent features. Nonnegativity is an essential constraint needed to obtain interpretable and sparse latent features that are parts–based representation of the data. The custom clustering in NTFk serves to estimate the number of latent features in the data.
We report on several efforts to search for physics beyond the standard model of particle physics at broad energy scales. The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the muon to electron conversion process, which is suppressed in the Standard Model. Mu2e will be operated at a low energy, yet can probe New Physics at very high mass scales (O(1e3 - 1e4 ) TeV). At high energies, the CMS experiment at the CERN LHC continues to deliver an impressive suite of Standard Model measurements and limits on a variety of New Physics signatures. Mu2e is under construction and slated to collect its first physics data in the coming years. This thesis describes work done during the construction phase of Mu2e and focuses on two critical areas: magnetic field modeling and statistical analysis. We describe a novel method for field modeling which we validate using a simulated dataset representing the expected magnetic field in the Detector Solenoid. This method blends a standard least-squares fitting technique that utilizes physically motivated analytical model functions with a novel physics informed network that is constructed to obey Maxwell’s equations. We show the technique can model the field with an accuracy of 10−7 despite the presence of injected noise in the pseudo-measurements at the 10−5 level. We then present preliminary results of the calibration of 3D Hall probes at the sub-10−4 level. These probes will be used to directly measure the Mu2e Detector Solenoid magnetic field on a sparse grid; these measurements serve as the input to the field model fitting. Finally, we describe the first implementation of both an unbinned shape analysis and a Bayesian interpretation applied to Mu2e pseudo-data. Up to 20% tighter limits can be set by the shape analysis compared to a standard cut & count analysis. The AlCap experiment collected data at PSI in 2015 to measure several important quantities related to nuclear muon capture on an aluminum target, which is a significant background process for Mu2e. The neutron emission from muon capture can introduce background hits in the Mu2e detectors and can increase radiation damage in various elements of the apparatus. We present measurements of the neutron group fluence and mean neutron multiplicity for muon capture on aluminum nuclei. Finally, we discuss an analysis of triboson production at CMS using an Effective Field Theory framework. Standard Model triboson production, which was first observed at CMS in 2020, has a relatively small cross section and provides direct access to both anomalous triple gauge couplings and quartic gauge couplings. These couplings, interpreted in the Standard Model Effective Field Theory, are studied in the present work. We target the boosted regime where the background rate is low and yields are enhanced when dimension-6 and dimension-8 Wilson coefficients are non-zero. We do not observe an excess in the data and therefore set bounds on the Wilson coefficients. For dimension-6 coefficients the tightest observed (expected) bounds are set on cW /Λ2 where Λ is the mass scale of new physics; the bounds are [−0.13, 0.12] TeV−2 ([−0.12, 0.12] TeV−2 ) at 95% CL. The tightest bounds in dimension-8 are set on fT,0 / Λ4 ; the observed (expected) bounds at 95% CL are [−0.63, 0.69] TeV−4 ([−0.54, 0.62] TeV−4 ). Additional results are presented which include scenarios where multiple Wilson coefficients are non-zero, the application of signal model clipping to address unitarity violation in Effective Field Theories, and a novel template fit developed for easier reinterpretation of our results.
We study p-string condensation mechanisms for fracton phases from the viewpoint of higher-form symmetry, focusing on the examples of the X-cube model and the rank-two symmetric-tensor scalar charge theory. This work is motivated by questions of the relationship between fracton phases and continuum quantum field theories, and also provides general principles to describe p-string condensation independent of specific lattice model constructions. We give a perspective on higher-form symmetry in lattice models in terms of cellular homology. Applying this perspective to the coupled-layer construction of the X-cube model, we identify a foliated 1-form symmetry that is broken in the X-cube phase, but preserved in the phase of decoupled toric code layers. Similar considerations for the scalar charge theory lead to a framed 1-form symmetry. These symmetries are distinct from standard 1-form symmetries that arise, for instance, in relativistic quantum field theory. Here, we also give a general discussion on interpreting p-string condensation, and related constructions involving gauging of symmetry, in terms of higher-form symmetry.
Motivated by the recent experiment on p-orbital-band bosons in optical lattices, we study theoretically the quantum phases of Mott insulator and superfluidity in two dimensions. The system features a superfluid phase with transversely staggered orbital current at weak interaction and a Mott insulator phase with antiferro-orbital order at strong coupling and commensurate filling. We go beyond mean-field theory and derive from a microscopic model an effective action that is capable of describing both the p-band Mott insulating and superfluid phases in strong coupling. We further calculate the excitation spectra near the quantum critical point and find two gapless modes away from the tip of the Mott lobe but four gapless modes at the tip. Our effective theory reveals how the phase coherence peak builds up in the Mott regime when approaching the critical point. We also discuss the finite-temperature phase transition of p-band superfluidity.
At large N{sub c}, cold nuclear matter is expected to form a crystal and thus spontaneously break translational symmetry. The description of chiral symmetry breaking and translational symmetry breaking can become intertwined. Here, the focus is on aspects of chiral symmetry breaking and its possible restoration that are by construction independent of the nature of translational symmetry breaking - namely spatial averages of chiral order parameters. A system will be considered to be chirally restored provided all spatially averaged chiral order parameters are zero. A critical question is whether chiral restoration in this sense is possible for phases in which chiral order parameters are locally nonzero but whose spatial averages all vanish. We show that this is not possible unless all chirally invariant observables are spatially uniform. This result is first derived for Skyrme-type models, which are based on a nonlinear sigma model and by construction break chiral symmetry on a point-by-point basis. A no-go theorem for chiral restoration (in the average sense) for all models of this type is obtained by showing that in these models there exist chirally symmetric order parameters that cannot be spatially uniform. Next, we will show that the no-go theorem applies to large N{sub c} QCD in any phase that has a nonzero but spatially varying chiral condensate. The theorem is demonstrated by showing that in a putative chirally restored phase, the field configuration can be reduced to that of a nonlinear sigma model. It is also shown that this no-go theorem is fully consistent with the vanishing of the spatial average of the chiral condensate (1/2)Tr(U) (as happens in 'half-skyrmion' configurations). This is because the chiral condensate is only one of an infinite set of chiral order parameters, some of which must be nonzero. It is also shown that while an approximation of a unit cell of a Skyrme crystal as a hypersphere does lead to a phase that is chirally restored (in the average sense), this is an artifact of the approximation.
We propose a novel approach to the solution of nematic Liquid Crystal models based on the derivation of a system of nonlinear wave equations for order parameters such that the occurrence of uniaxial and biaxial phase transitions can be interpreted as the propagation of a two-dimensional shock wave in the space of thermodynamic parameters. We obtain the exact equations of state for an integrable model of biaxial nematic liquid crystals and show that the classical transition from isotropic to uniaxial phase in absence of external fields is the result of a van der Waals type phase transition, where the jump in the order parameters is a classical shock generated from a gradient catastrophe at a non-zero isotropic field. The study of the equations of state provides the first analytical description of the rich structure of nematics phase diagrams in presence of external fields.
In F-term supergravity inflation models, scalar fields other than the inflaton generically receive a Hubble induced mass, which may restore gauge symmetries during inflation and phase transitions may occur during or after inflation as the Hubble parameter decreases. We study topological defect formation associated with such a phase transition in chaotic inflation in supergravity and obtain a severe constraint on the symmetry breaking scale which is related with the tensor-to-scalar ratio. This talk is based on the paper [1].
Turbulent flow coupled with heat transfer is investigated for a High Flux Isotope Reactor (HFIR) fuel plate. The Reynolds Averaged Navier-Stokes Models are used for fluid dynamics and the transfer of heat from a thermal nuclear fuel plate using the Multi-physics code COMSOL. Simulation outcomes are compared with experimental data from the Advanced Neutron Source Reactor Thermal Hydraulic Test Loop. The computational results for the High Flux Isotope Reactor core system provide a more physically accurate simulation of this system by modeling the turbulent flow field in conjunction with the diffusion of thermal energy within the solid and fluid phases of the model domain. Recommendations are made regarding Nusselt number correlations and material properties for future thermal hydraulic modeling efforts
Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.
We present a brief overview of axion models associated to anomalous abelian (gauge) symmetries, discussing their main phenomenological features. Among these, the mechanism of vacuum misalignment introduced at the QCD and at the electroweak phase transitions, with the appearance of periodic potentials, responsible for the generation of a mass for these types of axions.
Scalar modifications of gravity have an impact on the growth of structure. Baryon and Cold Dark Matter (CDM) perturbations grow anomalously for scales within the Compton wavelength of the scalar field. In the late time Universe when reionisation occurs, the spectrum of the 21 cm brightness temperature is thus affected. We study this effect for chameleon-f(R) models, dilatons and symmetrons. Although the f(R) models are more tightly constrained by solar system bounds, and effects on dilaton models are negligible, we find that symmetrons where the phase transition occurs before z{sub *} ∼ 12 could be detectable for a scalar field range as low as 5kpc. For all these models, the detection prospects of modified gravity effects are higher when considering modes parallel to the line of sight where very small scales can be probed. The study of the 21 cm spectrum thus offers a complementary approach to testing modified gravity with large scale structure surveys. Short scales, which would be highly non-linear in the very late time Universe when structure forms and where modified gravity effects are screened, appear in the linear spectrum of 21 cm physics, hence deviating from General Relativity in a maximal way.
Abstract Knowledge of the thermodynamic equilibria and domain structures of ferroelectrics is critical to establishing their structure–property relationships that underpin their applications from piezoelectric devices to nonlinear optics. Here, we establish the strain condition for strain phase separation and polydomain formation and analytically predict the corresponding domain volume fractions and wall orientations of, relatively low symmetry and theoretically more challenging, monoclinic ferroelectric thin films by integrating thermodynamics of ferroelectrics, strain phase equilibria theory, microelasticity, and phase‐field method. Using monoclinic K x Na 1 − x NbO 3 (0.5 < x < 1.0) thin films as a model system, we establish the polydomain strain–strain phase diagrams, from which we identify two types of monoclinic polydomain structures. The analytically predicted strain conditions of formation, domain volume fractions, and domain wall orientations for the two polydomain structures are consistent with phase‐field simulations and in good agreement with experimental results in the literature. The present study demonstrates a general, powerful analytical theoretical framework to predict the strain phase equilibria and domain wall orientations of polydomain structures applicable to both high‐ and low‐symmetry ferroelectrics and provide fundamental insights into the equilibrium domain structures of ferroelectric K x Na 1 − x NbO 3 thin films that are of technology relevance for lead‐free dielectric and piezoelectric applications.
Modern theories of everything, or theories of the grand unification of all physical interactions, try to describe the whole world starting from the first principles of quantum theory. However, the first principles operate with undetermined notions, such as the wave {psi}-function, particle, lepton and quark, de Broglie and Compton wavelengths, mass, electric charge, spin, electromagnetic field, photon, gravitation, physical vacuum, space, etc. From a logical point of view this means that such modern approach to the theory of everything is condemned to failure... Thus, what should we suggest to improve the situation? It seems quite reasonable to develop initially a theory of something, which will be able to clarify the major fundamental notions (listed above) that physics operates with every day. What would be a starting point in such approach? Of course a theory of space as such, because particles and all physical fields emerge just from space. After that, when a particle and fields (and hence the fields' carriers) are well defined and introduced in the well defined physical space, different kinds of interactions can be proposed and investigated. Moreover, we must also allow for a possible interaction of a created particle with the space that generated the appearance of the particle. The mathematical studies of Michel Bounias and the author have shown what the real physical space is, how the space is constituted, how it is arranged and what its elements are. Having constructed the real physical space we can then derive whatever we wish, in particular, such basic notions as mass, particle and charge. How are mechanics of such objects (a massive particle, a charged massive particle) organised? The appropriate theory of motion has been called a sub microscopic mechanics of particles, which is developed in the real physical space, not an abstract phase space, as conventional quantum mechanics does. A series of questions arise: can these two mechanics (submicroscopic and conventional quantum mechanics) be unified?, what can such unification bring new for us?, can such submicroscopic mechanics be a starting point for the derivation of the phenomenon of gravity?, can this new theory be a unified physical theory?, does the theory allow experimental verification? These major points have been clarified in detail. And, perhaps, the most intriguing aspect of the theory is the derivation of a new physical field associated with the notion of mass (or rather inertia of a particle, which has been called the inerton field and which represents a real sense of the particle's wave {psi}-function). This field emerges by analogy with the electromagnetic field associated with the notion of the electric charge. Yes, the postulated inerton field has being tested in a series of different experiments. Even more, the inerton field might have a number of practical applications...