Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite temperature”

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 469 records · Page 26

Hands-on with OWL: the Oak–Ridge Wang–Landau Monte Carlo software suite

The Oak–Ridge Wang–Landau (OWL) package is an open-source scientific software specialized for large-scale, Monte Carlo simulations for the study of materials properties at finite temperature. In this paper, we discuss the main features and capabilities of OWL, followed by detailed descriptions of building and running the code. The readers will be guided through the usage and functionality of the code with a few hands-on examples. This paper is based on a tutorial on OWL given at the 32nd Center for Simulational Physics Workshop on Recent Developments in Computer Simulation Studies in Condensed Matter Physics.

Li, Ying Wai↗

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗

Quantum Joule expansion of one-dimensional systems

Here, we investigate the Joule expansion of nonintegrable quantum systems that contain bosons or spinless fermions in one-dimensional lattices. A barrier initially confines the particles to be in half of the system in a thermal state described by the canonical ensemble and is removed at time t = 0. We investigate the properties of the time-evolved density matrix, the diagonal ensemble density matrix, and the corresponding canonical ensemble density matrix with an effective temperature determined by the total energy conservation using exact diagonalization. The weights for the diagonal ensemble and the canonical ensemble match well for high initial temperatures that correspond to negative effective final temperatures after the expansion. At long times after the barrier is removed, the time-evolved Rényi entropy of subsystems bigger than half can equilibrate to the thermal entropy with exponentially small fluctuations. The time-evolved reduced density matrix at long times can be approximated by a thermal density matrix for small subsystems. Few-body observables, like the momentum distribution function, can be approximated by a thermal expectation of the canonical ensemble with strongly suppressed fluctuations. The negative effective temperatures for finite systems go to non-negative temperatures in the thermodynamic limit for bosons, but there is a true thermodynamic effect for fermions, which is confirmed by finite-temperature density matrix renormalization group calculations. We propose the Joule expansion as a way to dynamically create negative temperature states for fermion systems with repulsive interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the instability and energy flux of lower hybrid waves in the Venus plasma mantle

Waves generated near the lower hybrid resonance frequency by the modified two stream instability have been invoked as a possible source of energy flux into the topside ionosphere of Venus. These waves are observed above the ionopause in a region known as the plasma mantle. The plasma within the mantle appears to be a mixture of magnetosheath and ionospheric plasmas. Since the magnetosheath electrons and ions have temperatures of several tens of eV, any instability analysis of the modified two stream instability requires the inclusion of finite electron and ion temperatures. Finite temperature effects are likely to reduce the growth rate of the instability. Furthermore, the lower hybrid waves are only quasi-electrostatic, and the energy flux of the waves is mainly carried by parallel Poynting flux. The magnetic field in the mantle is draped over the ionopause. Lower hybrid waves therefore cannot transport any significant wave energy to lower altitudes, and so do not act as a source of additional heat to the topside ionosphere.

Strangeway, R. J.↗

Determining microwave brightness temperatures from precipitating horizontally finite and vertically structured clouds

Microwave radiances that may be measured from satellite-borne radiometers operating at 37 GHz have been computed as a function of rainfall rates from horizontally finite precipitating clouds that contain both ice and liquid hydrometeors. It is found that precipitating ice at the top of the cloud depresses brightness temperatures significantly. Detailed comparison between finite clouds and the equivalent sections of plane-parallel clouds are made. Footprint averaged brightness temperatures from finite clouds are found to deviate considerably from the sectioned plane-parallel approximation, especially for large rain rates. Better agreement is achieved if a plane-parallel source function is used in the finite cloud model. The effect of shape and orientation of the precipitating cells with respect to the satellite has also been considered. It is found that the largest errors introduced by plane-parallel theory besides the footprint-filling errors are encountered when precipitating cells have large fractions of their surface below the freezing level or contain large amounts of ice.

Kummerow, C.↗

Real-time dynamics of the Schwinger model as an open quantum system with Neural Density Operators

Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1 + 1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NuScale Pressure and Temperature Limits Methodology Using Finite Element Analysis

Per 10 CFR 50 Appendix G, the pressure-temperature (P-T) limits curves and minimum temperature must be established to provide adequate margins for ferritic pressure-retaining components of the reactor coolant pressure boundary; this is to protect against brittle failure during any normal operating conditions, including anticipated operational occurrences and system hydrostatic tests, to which the pressure boundary may be subjected over its service lifetime. Specifically, ASME Code Section XI Appendix G procedures must be used for P-T limits calculation considering the pressure and temperature at various operating transient conditions. However, the elastic fracture mechanics solutions in Section XI are only suitable for cylindrical reactor pressure vessel (RPV) beltline without geometric discontinuities. Hence, these solutions are not suitable for postulated flaws near the core support blocks attached to the NuScale RPV inside surface, which is part of the beltline. As a result, NuScale has used finite element analysis to calculate thermal stress and stress intensity factor for the postulated flaws. The solutions using finite-element analysis have been validated using the formulations for straight cylinders for both axial and circumferential flaws. In addition, special considerations are given to the RPV beltline nil-ductility transition temperature (RTNDT) due to neutron irradiation occurring at lower temperatures than conventional plants. The paper summarizes the methodology and finite-element models used to develop P-T limits curves for NuScale RPV at the end of its 60-year design life.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal energy of a charm-meson molecule in a pion gas

The thermal corrections to the propagator of a loosely bound charm-meson molecule in a pion gas are calculated to next-to-leading order in the heavy-meson expansion using a zero-range effective field theory. Ultraviolet divergences in the charm-meson-pair self energy are canceled by corrections to the charm-meson-pair contact vertex. Terms that are singular at the charm-meson-pair threshold can be absorbed into thermal corrections to the rest energies and kinetic masses of the charm-meson constituents. The remaining terms reduce to a thermal correction to the binding momentum that is proportional to the pion number density and suppressed by the pion/charm-meson mass ratio. The correction gives a tiny decrease in the binding energy of the charm-meson molecule relative to the charm-meson-pair threshold in the pion gas and a change in its thermal width that is small compared to the thermal widths of the charm-meson constituents. These results are encouraging for the prospects of observing X(3872) and $T$$^{+}_{cc}$ (3875) in the expanding hadron gas produced by heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The LPM effect in sequential bremsstrahlung: 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ corrections

An important question concerning in-medium high-energy parton showers in a quark-gluon plasma or other QCD medium is whether consecutive splittings of the partons in a given shower can be treated as quantum mechanically independent, or whether the formation times for two consecutive splittings instead have significant overlap. Various previous calculations of the effect of overlapping formation times have either (i) restricted attention to a soft bremsstrahlung limit, or else (ii) used the large-N c limit (where N c = 3 is the number of quark colors). In this paper, we make a first study of the accuracy of the large-N c limit used by those calculations of overlap effects that avoid a soft bremsstrahlung approximation. Specifically, we calculate the 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ correction to previous N c = ∞ results for overlap g → gg → ggg of two consecutive gluon splittings g → gg. At order 1/$ {\mathrm{N}}_{\mathrm{c}}^2 $ , there is interesting and non-trivial color dynamics that must be accounted for during the overlap of the formation times.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The LPM effect in sequential bremsstrahlung: incorporation of “instantaneous” interactions for QCD

The splitting processes of bremsstrahlung and pair production in a medium are coherent over large distances in the very high energy limit, which leads to a suppression known as the Landau-Pomeranchuk-Migdal (LPM) effect. We continue study of the case when the coherence lengths (formation lengths) of two consecutive splitting processes overlap, avoiding soft-emission approximations. Previous work made a “nearly-complete” calculation of the effect of overlapping formation times on gluonic splittings such as g → gg → ggg (with simplifying assumptions such as an infinite QCD medium and the large-N c limit). In this paper, we extend those previous rate calculations from nearly-complete to complete by including processes involving the exchange of longitudinally-polarized gluons. In the context of Lightcone Pertubation Theory, used earlier for the “nearly-complete” calculation, such exchanges are instantaneous in lightcone time and have their own diagrammatic representation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An analysis of temperature effect in a finite journal bearing with spatial tilt and viscous dissipation

The analysis presented herein deals with the evaluation of the pressure, velocity, and temperature profiles in a finite-length plane journal bearing. The geometry of the case under study consists of a spatially tilted shaft. The two-dimensional Reynolds equation accounts for the variation of the clearance gap h with x and z and is used to model the pressure field. The latter is solved for a variety of shaft tilt angles and then used to calculate the two-dimensional flow field. Finally, the flow field is used in the energy equation to solve for the film temperature profile, when the effect of viscous dissipation is taken into account.

Braun, M. J.↗

Finite element simulation of temperature dependent free surface flows

The method of Engelman and Sani (1984) for a finite-element simulation of incompressible surface flows with a free and/or moving fluid interface, such as encountered in crystal growth and coating and polymer technology, is extended to temperature-dependent flows, including the effect of temperature-dependent surface tension. The basic algorithm of Saito and Scriven (1981) and Ruschak (1980) has been generalized and implemented in a robust and versatile finite-element code that can be employed with relative ease for the simulation of free-surface problems in complex geometries. As a result, the costly dependence on the Newton-Raphson algorithm has been eliminated by replacing it with a quasi-Newton iterative method, which nearly retains the superior convergence properties of the Newton-Raphson method.

Engelman, M. S.↗

Spontaneous breaking of finite group symmetries at all temperatures

We study conformal field theories with finite group symmetries with spontaneous symmetry breaking (SSB) phases which persist at all temperatures. We work with two \lambda \phi^4 λ ϕ 4 theories coupled through their mass terms. The two \lambda \phi^4 λ ϕ 4 theories are chosen to preserve either the cubic symmetry group, or the tetrahedral symmetry group. A one-loop calculation in the 4-\epsilon 4 − ϵ expansion shows that there exist infinitely many fixed points that can host an all temperature SSB phase. An analysis of the renormalization group (RG) stability matrix of these fixed points, reveals that their spectrum contains at least three relevant operators. In other words, these fixed points are tetracritical points.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The temperature distribution in an infinite medium resulting from a plane source of finite duration.

The temperature distribution in an infinite slab during and following heat generation by a plane source of finite duration is studied. An expression for the instantaneous plane source is obtained. By integration of this expression, the temperature distribution following the end of heating can be obtained. Further, the treatment is extended to time-varying heat generation functions, and an example is presented.

Montealegre, A.↗