Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase behavior”

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,063 records · Page 59

Circinus X-1 revisited: Fast-timing properties in relation to spectral state

We have studied the X-ray spectral and fast-timing variations of Cir X-1 by performing a homogenous analysis of all EXOSAT ME data on this source using X-ray hardness-intensity diagrams (HIDs), color-color diagrams (CDs), and power spectra. Cir X-1 exhibits a wide range of power spectral shapes and a large variety in X-ray spectral shapes. At different epochs the power spectra variously resemble those of an atoll source, a Z source, a black-hole candidate, or are unlike any of these. At some epochs one-dimensional connected-branch patterns are seen in HID and CD, and at other times more complex structures are found. We interpret the complex behavior of Cir X-1 in terms of a model where accretion rate, orbital phase and epoch are the main determinants of the source behavior, and where the unique properties of the source are due to two special circumstances: (1) the source is the only known atoll source (accreting neutron star with a very low magnetic field) that can reach the Eddington critical accretion rate, and (2) it has a unique, highly eccentric and probably precessing orbit. Property (1) makes Cir X-1 a very important source for our understanding of the similarities in the observable properties of neutron stars and black holes as it allows to separate out black hole signatures from properties that are merely due to the presence of accretion compact with a low magnetic field.

Oosterbroek, T.↗

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗

Salt Hydrate Eutectic Thermal Energy Storage for Building Thermal Regulation (Final Technical Report)

Thermal energy storage is anticipated to play an important role in developing the power grid of the future - a power grid that meets increasing demands of users, is resistant to disruptions, but also allows for greater penetration of renewable resources. Specifically, thermal energy storage materials can be integrated into HVAC systems and building envelopes, where they can be used to shift power demands for building climate control from periods of peak demand to periods of low demand. Phase change materials (PCMs) are compelling as low-cost, high energy density thermal energy storage materials for building thermal management. However, there is a lack of high performance low-cost PCMs within the specific temperature ranges which would most effectively allow for power load shifting. Inorganic salt hydrates represent a promising class of PCMs, but their inherent limitations cause them to be currently unavailable for reliable building applications. The overarching goals of this research effort are to: 1) Discover low-cost, high volumetric density salt hydrate eutectic PCMs to store low-quality heat (10 to 40 °C); 2) Introduce a high thermal conductivity matrix to reduce the time constant for energy storage to ~0.1 to 1 hr, incorporate nucleation catalysts to decrease undercooling, and utilize microencapsulation and shape stabilization approaches, minimizing moisture loss/gain, mitigating phase separation, and maintaining stable melting behavior over the lifetime of the compounds; 3) Evaluate the impact these systems have on peak load shifting, and the potential for overall energy savings under different climatic scenarios and building configurations. These goals will be achieved by an integrated research program consisting of six cohesive research subtasks: 1) Materials discovery of eutectic salt hydrate PCMs by using computationally predicted thermodynamic equilibria, coupled with high-throughput experimental validation, 2) Rapid experimental screening of nucleation catalysts identified through robust computational databases, 3) Embedding salt hydrate PCM into a low cost and scalable high conductivity matrix, 4) Microencapsulation of salt hydrate microspheres using hybrid inorganic-polymer microencapsulation approach, 5) Shape stabilization by thermoreversible salt hydrate salogels, and 6) Analysis of end-use using thermal simulations, and characterization of mock-up energy storage finished components.

25 ENERGY STORAGE↗

Raman study of the improper ferroelectric phase transition in iron iodine boracite

The reduction in the dielectric anomaly near the phase transition temperature in improper ferroelectrics suggests that this behavior may be exploited to significantly enhance the figure-of-merit for pyroelectric infrared detection. This reduction is associated with a decoupling of the thermal and dielectric responses which, in turn, reflects the unusual nature of the lattice instability at the phase transition. In order to elucidate the characteristics of the phonon behavior connected with the lattice instability, Raman scattering of Fe3B7O13I has been examined in both the paraelectric and ferroelectric phases. The results of this study as well as results from earlier measurements on other boracites, provides insight into the nature of the mode softening.

Kim, Q.↗

Modeling void growth and movement with phase change in thermal energy storage canisters

A scheme was developed to model the thermal hydrodynamic behavior of thermal energy storage salts. The model included buoyancy, surface tension, viscosity, phases change with density difference, and void growth and movement. The energy, momentum, and continuity equations were solved using a finite volume formulation. The momentum equation was divided into two pieces. The void growth and void movement are modeled between the two pieces of the momentum equations. Results showed this scheme was able to predict the behavior of thermal energy storage salts.

Darling, Douglas↗

Modeling void growth and movement with phase change in thermal energy storage canisters

A scheme was developed to model the thermal hydrodynamic behavior of thermal energy storage salts. The model included buoyancy, surface tension, viscosity, phases change with density difference, and void growth and movement. The energy, momentum, and continuity equations were solved using a finite volume formulation. The momentum equation was divided into two pieces. The void growth and void movement are modeled between the two pieces of the momentum equations. Results showed this scheme was able to predict the behavior of thermal energy storage salts.

Darling, Douglas↗

Condensed matter physics of planets - Puzzles, progress and predictions

Attention is given to some of the major unresolved issues concerned with the physics of planetary interiors. The important advances in observations, and experimental and theoretical investigations are briefly reviewed, and some areas for further study are identified, including: the characteristics of atomic and electronic degrees of freedom at the high pressures and temperatures typical of a condensed planetary core; the behavior of water at megabar pressures; and the nature of the core-alloy in the earth and in the core mantle phase boundary. Consideration is also given to the behavior of carbon at high pressures and temperatures in the presence of oxygen and hydrogen; the behavior of the volatile ice assemblage in Titan at pressures of 2-40 kbar; and the electrical conductivities of matter under planetary core conditions.

Stevenson, D. J.↗

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microwave Oscillators Based on Nonlinear WGM Resonators

Optical oscillators that exploit resonantly enhanced four-wave mixing in nonlinear whispering-gallery-mode (WGM) resonators are under investigation for potential utility as low-power, ultra-miniature sources of stable, spectrally pure microwave signals. There are numerous potential uses for such oscillators in radar systems, communication systems, and scientific instrumentation. The resonator in an oscillator of this type is made of a crystalline material that exhibits cubic Kerr nonlinearity, which supports the four-photon parametric process also known as four-wave mixing. The oscillator can be characterized as all-optical in the sense that the entire process of generation of the microwave signal takes place within the WGM resonator. The resonantly enhanced four-wave mixing yields coherent, phase-modulated optical signals at frequencies governed by the resonator structure. The frequency of the phase-modulation signal, which is in the microwave range, equals the difference between the frequencies of the optical signals; hence, this frequency is also governed by the resonator structure. Hence, further, the microwave signal is stable and can be used as a reference signal. The figure schematically depicts the apparatus used in a proof-of-principle experiment. Linearly polarized pump light was generated by an yttrium aluminum garnet laser at a wavelength of 1.32 microns. By use of a 90:10 fiber-optic splitter and optical fibers, some of the laser light was sent into a delay line and some was transmitted to one face of glass coupling prism, that, in turn, coupled the laser light into a crystalline CaF2 WGM disk resonator that had a resonance quality factor (Q) of 6x10(exp 9). The output light of the resonator was collected via another face of the coupling prism and a single-mode optical fiber, which transmitted the light to a 50:50 fiber-optic splitter. One output of this splitter was sent to a slow photodiode to obtain a DC signal for locking the laser to a particular resonator mode. The other output of this splitter was combined with the delayed laser signal in another 50:50 fiber-optic splitter used as a combiner. The output.of the combiner was fed to a fast photodiode that demodulated light and generated microwave signal. In this optical configuration, the resonator was incorporated into one arm of a Mach-Zehnder interferometer, which was necessary for the following reasons: It was found that when the output of the resonator was sent directly to a fast photodiode, the output of the photodiode did not include a measurable microwave signal. However, when the resonator was placed in an arm of the interferometer and the delay in the other arm was set at the correct value, the microwave signal appeared. Such behavior is distinctly characteristic of phase-modulated light. The phase-modulation signal had a frequency of about 8 GHz, corresponding to the free spectral range of the resonator. The spectral width of this microwave signal was less than 200 Hz. The threshold pump power for generating the microwave signal was about 1 mW. It would be possible to reduce the threshold power by several orders of magnitude if resonators could be made from crystalline materials in dimensions comparable to those of micro-resonators heretofore made from fused silica.

Maleki, Lute↗

Towards understanding the development and characteristics of under-expanded flash boiling jets

Here in this research, the n-hexane jets discharged from a straight channel were studied numerically to reveal the development and characteristics of under-expanded flashing jets. Firstly, the development of flash boiling jets was investigated and compared to that of two-phase jets without phase transition. It was shown that the two-phase jets without phase transition can be under-expanded only if the flow is choked. In contrast, flash boiling jets can be under-expanded whether or not the choked conditions are verified. Secondly, the factors affecting the flashing jet behaviors were discussed. The relaxation time in phase transition, the expansion ratio, and the injection pressure have been identified as three key factors determining the jet characteristics. With high enough expansion ratio and low enough relaxation time, the expansion was so strong that the flow became transonic and shock structures similar to those of under-expanded gaseous jets appeared near the orifice exit. Finally, two types of choking have been identified for flash boiling jets, i.e., vaporization-induced choking and expansion-induced choking. The former one is mainly ascribed to the vapor formation in the orifice which reduces the local speed of sound, while the latter one is due to the expansion of under-expanded two-phase flow and the resultant flow acceleration.

42 ENGINEERING↗

Recovery of the Disrupted Quasi-Biennial Oscillation

There occurred a unique and significant disruption of the Quasi-Biennial Oscillation (QBO) during the Northern Hemisphere winter of 2015-16. Here we document the return of the QBO to its normal downward phase speed and period based on Singapore soundings, MERRA-2 re-analysis (Modern Era Retrospective Reanalysis for Research and Applications), and a simple QBO model. Daily averaged zonal winds from 100-10 hPa are used to characterized the behavior of the QBO's amplitude and phase as seen in the first two EOFs (Empirical Orthogonal Functions). These EOFs capture the QBO structure and evolution during the pre-disruption, disruption, and post-disruption times. Results show that the amplitude and phase returned to normal by June 2016, however the post-disruption QBO phase was delayed relative to the pre-disruption phase by four tenths of a QBO cycle (~11 months). A rapid, seasonal, phase shift of this magnitude is shown to be unique in the QBO observational record.

Coy, Lawrence↗

Cross phases of temperature-gradient-driven turbulence as a model basis for I -mode particle transport

As a model for understanding the type of transport behavior characteristic of the tokamak I mode, cross-phase physics for particle-transport is studied analytically for turbulence dominated by either ion-temperature-gradient (ITG) or electron-temperature-gradient (ETG) instability. I mode is a transport-barrier regime of reduced thermal transport but essentially unaffected particle transport. It is assumed that ITG turbulence applies to the baseline L mode, ETG to I mode, and that E × B flow shear is stronger in I mode, lowering all fluxes. In ITG turbulence, particle transport is governed by trapped electrons. Sensitivity to collisions produces the well-known temperature-gradient-driven pinch that offsets density-gradient-driven outward diffusion, weakening particle transport in L mode. In ETG turbulence, nonadiabatic ions are collisionless. Nonzero transport requires an ion spectrum feature whose magnetic-drift resonance supplies the necessary cross phase. If frequencies of order the ion diamagnetic drift frequency dominate the ion part of the spectrum, as would occur with weakly unstable ITG turbulence, all components of the particle transport are outward and can offset flow-shear-induced flux reductions to produce a flux that is similar to the ITG L-mode particle flux. Nonlinear frequencies are potentially relevant and discussed in relation to I mode.

Physics↗

Ultrafast spin-nematic and ferroelectric phase transitions induced by femtosecond light pulses

We report, optically induced phase transitions of the manganite Pr 1/3 Ca 2/3 MnO 3 have been simulated by using a model Hamiltonian that captures the dynamics of strongly correlated charge, orbital, lattice, and spin degrees of freedom. Its parameters have been extracted from first-principles calculations. Beyond a critical intensity of a femtosecond light pulse, the material undergoes an ultrafast and nonthermal magnetic phase transition from a noncollinear to collinear antiferromagnetic phase. The light-pulse excites selectively either a spin-nematic or a ferroelectric phase, depending on the light polarization. The behavior can be traced to an optically induced ferromagnetic coupling between Mn trimers, i.e., polarons which are delocalized over three Mn sites. The polarization guides the polymerization of the polaronic crystal into distinct patterns of ferromagnetic chains determining the target phase.

36 MATERIALS SCIENCE↗

Metamagnetism of few-layer topological antiferromagnets

MnBi 2 Te 4 (MBT) materials are a promising class of antiferromagnetic topological insulators whose films provide access to novel and technologically important topological phases, including quantum anomalous Hall states and axion insulators. MBT device behavior is expected to be sensitive to the various collinear and noncollinear magnetic phases that are accessible in applied magnetic fields. Here, we use classical Monte Carlo simulations and electronic structure models to calculate the ground state magnetic phase diagram as well as topological and optical properties for few-layer films with up to six septuple layers. Using magnetic interaction parameters appropriate for MBT, we find that it is possible to prepare a variety of different magnetic stacking sequences, some of which have sufficient symmetry to disallow nonreciprocal optical response and Hall transport coefficients. Other stacking arrangements do yield large Faraday and Kerr signals, even when the ground state Chern number vanishes.

36 MATERIALS SCIENCE↗

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Observing the advection of sea ice in the Weddell Sea using buoy and satellite passive microwave data

Data from four buoys tracked by Nimbus 6 and concurrent ice concentrations retrieved from Nimbus 7 scanning multichannel microwave radiometer data are used to investigate the progress and behavior of an area of sea ice as it drifts from the southwestern Weddell Sea. The overall drift characteristics and their relationship to ice edge displacement are examined within the framework of four zones. Three phases are identified in the large-scale behavior of the Weddell Sea ice cover, namely, a rapid equatorward and eastward advance, a quasi-equilibrium phase, and a period of rapid recession. Outbreaks of cold continental air alternate with incursions of relatively warm air from the north; warm conditions are recorded as far as 1200 km in from the ice edge in winter. Closed loops in the buoy trajectories, which are clockwise to the south of 63 deg S, reverse to become anticlockwise to the north. A coherence is observed in the response of the buoys to the passage of storms, even though the buoys separated by a distance of over 100 km.

Massom, Robert A.↗

Phase transition phenomenology with nonparametric representations of the neutron star equation of state

Astrophysical observations of neutron stars probe the structure of dense nuclear matter and have the potential to reveal phase transitions at high densities. Most recent analyses are based on parametrized models of the equation of state with a finite number of parameters and occasionally include extra parameters intended to capture phase-transition phenomenology. However, such models restrict the types of behavior allowed and may not match the true equation of state. We introduce a complementary approach that extracts phase transitions directly from the equation of state without relying on, and thus being restricted by, an underlying parametrization. We then constrain the presence of phase transitions in neutron stars with astrophysical data. Current pulsar mass, tidal deformability, and mass-radius measurements disfavor only the strongest of possible phase transitions (latent energy per particle ≳ 100 MeV). Weaker phase transitions are consistent with observations. Here we further investigate the prospects for measuring phase transitions with future gravitational-wave observations and find that catalogs of O (100) events will (at best) yield Bayes factors of ~ 10 : 1 in favor of phase transitions even when the true equation of state contains very strong phase transitions. Our results reinforce the idea that neutron star observations will primarily constrain trends in macroscopic properties rather than detailed microscopic behavior. Fine-tuned equation of state models will likely remain unconstrained in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗

Holographic Kolmogorov-Sinai entropy and the quantum Lyapunov spectrum

In classical chaotic systems the entropy, averaged over initial phase space distributions, follows a universal behavior. While approaching thermal equilibrium it passes through a stage where it grows linearly, while the growth rate, the Kolmogorov-Sinai entropy (rate), is given by the sum over all positive Lyapunov exponents. A natural question is whether a similar relation is valid for quantum systems. We argue that the Maldacena-Shenker-Stanford bound on quantum Lyapunov exponents implies that the upper bound on the growth rate of the entropy, averaged over states in Hilbert space that evolve towards a thermal state with temperature T, should be given by πT times the thermal state’s von Neumann entropy. Strongly coupled, large N theories with black hole duals should saturate the bound. To test this we study a large number of isotropization processes of random, spatially homogeneous, far from equilibrium initial states in large N, $\mathcal{N}$ = 4 Super Yang Mills theory at strong coupling and compute the ensemble averaged growth rate of the dual black hole’s apparent horizon area. We find both an analogous behavior as in classical chaotic systems and numerical evidence that the conjectured bound on averaged entropy growth is saturated granted that the Lyapunov exponents are degenerate and given by λ i = ±2πT. This fits to the behavior of classical systems with plus/minus symmetric Lyapunov spectra, a symmetry which implies the validity of Liouville’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗