Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power system oscillation”

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 685 records · Page 38

Performance of a High-Fidelity 4kW-Class Engineering Model PPU and Integration with HiVHAc System

The High Voltage Hall Accelerator (HiVHAc) propulsion system consists of a thruster, power processing unit (PPU), and propellant feed system. An engineering model PPU was developed by Colorado Power Electronics, Inc. funded by NASA's Small Business Innovative Research Program. This PPU uses an innovative 3-phase resonant converter to deliver 4 kW of discharge power over a wide range of input and output voltage conditions. The PPU includes a digital control interface unit that automatically controls the PPU and a xenon flow control module (XFCM). It interfaces with a control computer to receive highlevel commands and relay telemetry through a MIL-STD-1553B interface. The EM PPU was thoroughly tested at GRC for functionality and performance at temperature limits and demonstrated total efficiencies a high as 95 percent. Integrated testing of the unit was performed with the HiVHAc thruster and the XFCM to demonstrate closed-loop control of discharge current with anode flow. Initiation of the main discharge and power throttling were also successfully demonstrated and discharge oscillations were characterized.

Power Processing Unit↗

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces

The scattering-type scanning near-field optical microscope (s-SNOM) has emerged as a powerful tool for resolving nanoscale inhomogeneities in laterally heterogeneous samples. However, most analytical models used to predict the scattering near-field signals are assuming homogenous landscapes (bulk materials), resulting in inconsistencies when applied to samples with more complex configurations. In this work, we combine the point-dipole model (PDM) to the finite-element method (FEM) to account for the lateral and vertical heterogeneities while keeping the computation time manageable. Full images, spectra, or hyperspectral line profiles can be simulated by calculating the self-consistent dipole radiation demodulated at higher harmonics of the tip oscillation, mimicking real experimental procedures. Using this formalism, we clarify several important yet puzzling experimental observations in near-field images on samples with rich typography and complex material compositions, heterostructures of two-dimensional material flakes, and plasmonic antennas. The developed method serves as a basis for future investigations of nano-systems with nontrivial topography.

36 MATERIALS SCIENCE↗

Dielectric breakdown induced by picosecond laser pulses

The damage thresholds of transparent optical materials were investigated. Single picosecond pulses at 1.06 microns, 0.53 microns and 0.35 microns were obtained from a mode locked Nd-YAG oscillator-amplifier-frequency multiplier system. The pulses were Gaussian in space and time and permitted the determination of breakdown thresholds with a reproducibility of 15%. It was shown that the breakdown thresholds are characteristic of the bulk material, which included nine alkali halides, five different laser host materials, KDP, quartz, sapphire and calcium fluoride. The extension of the damage data to the ultraviolet is significant, because some indication was obtained that two- and three-photon absorption processes begin to play a role in determining the threshold. Throughout the visible region of the spectrum the threshold is still an increasing function of frequency, indicating that avalanche ionization is the dominant factor in determining the breakdown threshold. This was confirmed by a detailed study of the damage morphology with a high resolution microscope just above the threshold. The influence of self focusing is discussed, and evidence for beam distortion below the power threshold for complete self focusing is presented, confirming the theory of Marburger.

Smith, W. L.↗

Improvements in Speed and Functionality of a 670-GHz Imaging Radar

Significant improvements have been made in the instrument originally described in a prior NASA Tech Briefs article: Improved Speed and Functionality of a 580-GHz Imaging Radar (NPO-45156), Vol. 34, No. 7 (July 2010), p. 51. First, the wideband YIG oscillator has been replaced with a JPL-designed and built phase-locked, low-noise chirp source. Second, further refinements to the data acquisition and signal processing software have been performed by moving critical code sections to C code, and compiling those sections to Windows DLLs, which are then invoked from the main LabVIEW executive. This system is an active, single-pixel scanned imager operating at 670 GHz. The actual chirp signals for the RF and LO chains were generated by a pair of MITEQ 2.5 3.3 GHz chirp sources. Agilent benchtop synthesizers operating at fixed frequencies around 13 GHz were then used to up-convert the chirp sources to 15.5 16.3 GHz. The resulting signals were then multiplied 36 times by a combination of off-the-shelf millimeter- wave components, and JPL-built 200- GHz doublers and 300- and 600-GHz triplers. The power required to drive the submillimeter-wave multipliers was provided by JPL-built W-band amplifiers. The receive and transmit signal paths were combined using a thin, high-resistivity silicon wafer as a beam splitter. While the results at present are encouraging, the system still lacks sufficient speed to be usable for practical applications in a contraband detection. Ideally, an image acquisition speed of ten seconds, or a factor of 30 improvement, is desired. However, the system improvements to date have resulted in a factor of five increase in signal acquisition speed, as well as enhanced signal processing algorithms, permitting clearer imaging of contraband objects hidden underneath clothing. In particular, advances in three distinct areas have enabled these performance enhancements: base source phase noise reduction, chirp rate, and signal processing. Additionally, a second pixel was added, automatically reducing the imaging time by a factor of two. Although adding a second pixel to the system doubles the amount of submillimeter components required, some savings in microwave hardware can be realized by using a common low-noise source.

Dengler, Robert J.↗

Light Echos in Kerr Geometry: A Source of High Frequency QPOs from Random X-ray Bursts

We propose that high frequency quasi-periodic oscillations (HFQPOs) can be produced from randomly-formed X-ray bursts (flashes) by plasma interior to the ergosphere of a rapidly-rotating black hole. We show by direct computation of their orbits that the photons comprising the observed X-ray light curves, if due to a multitude of such flashes, are affected significantly by the black hole's dragging of inertial frames; the photons of each such burst arrive to an observer at infinity in multiple (double or triple), distinct 'bunches' separated by a roughly constant time lag of t/M approximately equal to 14, regardless of the bursts' azimuthal position. We argue that every other such 'bunch' represents photons that follow trajectories with an additional orbit around the black hole at the photon circular orbit radius (a photon 'echo'). The presence of this constant lag in the response function of the system leads to a QPO feature in its power density spectra, even though the corresponding light curve consists of a totally stochastic signal. This effect is by and large due to the black hole spin and is shown to gradually diminish as the spin parameter a decreases or the radial position of the burst moves outside the static limit surface (ergosphere). Our calculations indicate that for a black hole with Kerr parameter of a/M=0.99 and mass of M=10*Msun the QPO is expected at a frequency of approximately 1.3-1.4 kHz. We discuss the plausibility and observational implications of our model/results as well as its limitations.

Fukumura, K.↗

Wind Energy Forecasting with the Weather Research and Forecasting Model

This was a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Siemens Energy, Inc. (Siemens) to develop a wind resource forecasting tool. LLNL was to develop an independent high-resolution mesoscale modeling capability forecasting tool that could be implemented in conjunction with existing wind farm control and monitoring software to provide forecasting of wind resources using local observations of winds and temperature. Research with LLNL’s state-of-the-art large-eddy simulation meteorological prediction model, based on the community WRF model and innovative turbulence parameterizations, would improve that model’s applicability to large wind farms offshore and in complex terrain. The modeling capability would include uncertainty quantification. Finally, the application of the modeling tool and existing global climate change predictions would enable the delineation of the likely effects of climate change on wind resources. Siemens was to provide high time resolution hub-height wind speed and other meteorological data streams, including temperature profiles from wind farms, for LLNL to incorporate into the modeling system, to validate and tune this forecasting model for their locations of interest. These data streams would also be used for longer-term studies of correlations of wind resources to climate oscillations to indicate how long-term climate change trends may affect the available wind resource. Siemens would also provide information and observations of turbine wakes for incorporation into the modeling tool. By implementing state-of-the-art turbulence parameterizations into a simulation model and/or ensembles of simulation models, and by integrating real-time hub height wind speed and other meteorological datastreams from wind farms into that model or ensemble of models, LLNL would develop a forecasting tool that could be implemented by Siemens as an add-on to existing wind farm control and monitoring software to provide owners with useful resource forecasting. The desired outcome was that the accuracy level of the output would be sufficient to substantiate power output commitments. The final deliverable for this work would consist of a document outlining the algorithms and software tools that could be integrated into Siemens Wind Park Supervisor.

17 WIND ENERGY↗

LACE flight dynamics experiment

The Low Power Atmospheric Compensation Experiment (LACE) is scheduled for launch in late 1989 into a 556 km altitude circular orbit of 43 deg inclination. The LACE flight dynamics experiment is an experiment secondary to the primary LACE mission. The purpose of the experiment is to provide on-orbit systems identification of the LACE spacecraft. The structure of the LACE spacecraft is of special interest to the CSI community. It incorporates 3 deployable/retractable booms of maximum length 45.72 m (150 ft) mounted on a rectangular parallelepiped bus of mass 1,200 kg. The zenith directed gravity gradient boom is mounted on the top of the bus; the retroreflector boom is mounted forward and deployed along the velocity vector, the balance boom is mounted and pointed aft. Attitude stabilization is accomplished by means of gravity gradient torques and by a momentum wheel. The LACE flight dynamics experiment is designed to measure modal frequencies, damping ratios, and oscillation amplitudes of the LACE spacecraft, as well as the vibration intensity generated by boom deployments and retractions. It is anticipated that this experiment will provide an opportunity for improvements in the accuracy of computer simulations of flexible structures and multibody dynamics.

Fisher, Shalom↗

SALSA - A Synthesis Array for Lunar Submillimeter Astronomy

A conceptual design is described for a lunar submillimeter wavelength interferometer called SALSA, a Synthesis Array for Lunar Submillimeter Astronomy. Its design is greatly simplified over conventional submillimeter wavelength arrays because it takes advantage of a beam waveguide to relay signals to a centrally-located receiver system. The array is optimized to synthesize images at wavelengths between 60 and 300 microns, and it has an angular resolution of 10 milliarcsecs at 60 microns. Observations in this region provide unique insights into fundamental issues in astrophysics such as star formation, but are impossible from the earth's surface because of atmospheric absorption. The baseline design for SALSA consists of twelve, 3.5-meter diameter antennas arranged in a Y-shaped configuration consisting of three 0.5-km long arms, each with four antennas distributed according to a power-law function of distance. Such a beam-waveguide approach has significant advantages, since all the high-power-consumption items (active cryogenics, correlator) are at a central location and could share a single power source. In addition, the antennas can be quite simple and maintenance is greatly reduced. The primary technology challenges for developing such an array are the coherent detectors (mixers and local oscillator sources), and ensuring the dimensional stability of the optical elements under the extreme lunar thermal environment (85-385 K).

Mahoney, M. J.↗

Influence of modulation frequency in rubidium cell frequency standards

The error signal which is used to control the frequency of the quartz crystal oscillator of a passive rubidium cell frequency standard is considered. The value of the slope of this signal, for an interrogation frequency close to the atomic transition frequency is calculated and measured for various phase (or frequency) modulation waveforms, and for several values of the modulation frequency. A theoretical analysis is made using a model which applies to a system in which the optical pumping rate, the relaxation rates and the RF field are homogeneous. Results are given for sine-wave phase modulation, square-wave frequency modulation and square-wave phase modulation. The influence of the modulation frequency on the slope of the error signal is specified. It is shown that the modulation frequency can be chosen as large as twice the non-saturated full-width at half-maximum without a drastic loss of the sensitivity to an offset of the interrogation frequency from center line, provided that the power saturation factor and the amplitude of modulation are properly adjusted.

Audoin, C.↗

RF Device for Acquiring Images of the Human Body

A safe, non-invasive method for forming images through clothing of large groups of people, in order to search for concealed weapons either made of metal or not, has been developed. A millimeter wavelength scanner designed in a unique, ring-shaped configuration can obtain a full 360 image of the body with a resolution of less than a millimeter in only a few seconds. Millimeter waves readily penetrate normal clothing, but are highly reflected by the human body and concealed objects. Millimeter wave signals are nonionizing and are harmless to human tissues when used at low power levels. The imager (see figure) consists of a thin base that supports a small-diameter vertical post about 7 ft (=2.13 m) tall. Attached to the post is a square-shaped ring 2 in. (=5 cm) wide and 3 ft (=91 cm) on a side. The ring is oriented horizontally, and is supported halfway along one side by a connection to a linear bearing on the vertical post. A planar RF circuit board is mounted to the inside of each side of the ring. Each circuit board contains an array of 30 receivers, one transmitter, and digitization electronics. Each array element has a printed-circuit patch antenna coupled to a pair of mixers by a 90 coupler. The mixers receive a reference local oscillator signal to a subharmonic of the transmitter frequency. A single local oscillator line feeds all 30 receivers on the board. The resulting MHz IF signals are amplified and carried to the edge of the board where they are demodulated and digitized. The transmitted signal is derived from the local oscillator at a frequency offset determined by a crystal oscillator. One antenna centrally located on each side of the square ring provides the source illumination power. The total transmitted power is less than 100 mW, resulting in an exposure level that is completely safe to humans. The output signals from all four circuit boards are fed via serial connection to a data processing computer. The computer processes the approximately 1-MB data set into a three-dimensional image in a matter of seconds. The innovation is to configure the receiver array in a ring topology surrounding the scanned object. The ring is then scanned vertically to cover the necessary two-dimensional surface. This fabrication of the ring is made possible by using planar antenna and circuit technology. A planar circuit board serves as a medium for both antennas and signal processing components. Using this technique, parts counts are kept low, and the cost per element is a small fraction of a waveguide-based system.

Gaier, Todd C.↗

Synoptic Meteorology Explains Temperate Forest Carbon Uptake

Abstract While substantial attention has been paid to the effects of both global climate oscillations and local meteorological conditions on the interannual variability of ecosystem carbon exchange, the relationship between the interannual variability of synoptic meteorology and ecosystem carbon exchange has not been well studied. Here we use a clustering algorithm to identify a summertime cyclonic precipitation system northwest of the Great Lakes to determine (a) the association at a daily scale between the occurrence of this system and the local meteorology and net ecosystem exchange at three Great Lakes region forested eddy covariance sites and (b) the association between the seasonal prevalence of this system and the summertime net ecosystem exchange of these sites. We find that temperature, in addition to precipitation and cloud cover, is an important explanatory factor for the suppression of net ecosystem productivity that occurs during these cyclonic events in this region. In addition, the prevalence of this cyclonic system can explain a significant proportion of the interannual variability in summertime forest ecosystem exchange in this region. This explanatory power is not due to a simple accumulation of low‐productivity days that cooccur with this meteorological event, but rather a broader association between the frequency of these events and several aspects of prevailing seasonal conditions. This work demonstrates the usefulness of conceptualizing meteorology in terms of synoptic systems for explaining the interannual variability of regional carbon fluxes.

Randazzo, Nina A.↗

Backscatter-Immune Injection-Locked Brillouin Laser in Silicon

As self-sustained oscillators, lasers possess the unusual ability to spontaneously synchronize. These nonlinear dynamics are the basis for a simple yet powerful stabilization technique known as injection locking, in which a laser’s frequency and phase can be controlled by an injected signal. Because of its inherent simplicity and favorable noise characteristics, injection locking has become a workhorse for coherent amplification and high-fidelity signal synthesis in applications ranging from precision atomic spectroscopy to distributed sensing. Within integrated photonics, however, these injection-locking dynamics remain relatively untapped—despite significant potential for technological and scientific impact. In this work, we demonstrate injection locking in a silicon photonic Brillouin laser. Injection locking of this monolithic device is remarkably robust, allowing us to tune the laser emission by a significant fraction of the Brillouin gain bandwidth. Harnessing these dynamics, we demonstrate amplification of small signals by more than 23 dB. Moreover, we demonstrate that the injection-locking dynamics of this system are inherently nonreciprocal, yielding unidirectional control and backscatter immunity in an all-silicon system. This device physics opens the door to strategies for phase-noise reduction, low-noise amplification, and backscatter immunity in silicon photonics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Implementation and Initial Validation of a 100-Kilowatt Class Nested-Channel Hall Thruster

The X3 is a 100-kilowatt class nested-channel Hall thruster developed by the Plasmadynamics and Electric Propulsion Laboratory at the University of Michigan in collaboration with the Air Force Research Laboratory and NASA. The cathode, magnetic circuit, boron nitride channel rings, and anodes all required specific design considerations during thruster development, and thermal modeling was used to properly account for thermal growth in material selection and component design. A number of facility upgrades were required at the University of Michigan to facilitate operation of the X3. These upgrades included a re-worked propellant feed system, a completely redesigned power and telemetry break-out box, and numerous updates to thruster handling equipment. The X3 was tested on xenon propellant at two current densities, 37% and 73% of the nominal design value. It was operated to a maximum steady-state discharge power of 60.8 kilowatts. The tests presented here served as an initial validation of thruster operation. Thruster behavior was monitored with telemetry, photography and high-speed current probes. The photography showed a uniform plume throughout testing. At constant current density, reductions in mass flow rate of 18% and 26% were observed in the three-channel operating configuration as compared to the superposition of each channel running individually. The high-speed current probes showed that the thruster was stable at all operating points and that the channels influence each other when more than one is operating simultaneously. Additionally, the ratio of peak-to-peak AC-coupled discharge current oscillations to mean discharge current did not exceed 51% for any operating points reported here, and did not exceed 17% at the higher current density.

High Power↗

A phenomenological wobbling model for isolated pulsars and the braking index

ABSTRACT An isolated pulsar is a rotating neutron star possessing a very high magnetic dipole moment, thus providing a powerful radiating mechanism. These stars loose rotational energy E through various processes, including a plasma wind originating from a highly magnetized magnetosphere and the emission of magnetic dipole radiation (MDR). Such phenomena produce a time decreasing angular velocity Ω(t) of the pulsar that is usually quantified in terms of its braking index. Although these mechanisms are widely acknowledged as the primary drivers of the spin evolution of isolated pulsars, it is plausible that other contributing factors influencing this effect have yet to be comprehensively investigated. Most of young isolated pulsars present a braking index different from that given by the MDR and plasma wind processes. Working in the weak field (Newtonian) limit, we take in this work a step forward in describing the evolution of such a system by allowing the star’s shape to wobble around an ellipsoidal configuration as a backreaction effect produced by its rotational deceleration. It is assumed that an internal damping of the oscillations occurs, thus introducing another form of energy loss in the system, and this phenomenon may be related to the deviation of the braking index from the models based on $\dot{E} \sim -\Omega ^4$ predictions. Numerical calculations suggest that the average braking index for typical isolated pulsars can be thus explained by a simple phenomenological model.

Araujo, E. C. A. (ORCID:0000000213520879)↗

Effect of Adding a Regenerator to Kornhauser's MIT "Two-Space" (Gas-Spring+Heat Exchanger) Test Rig

This study employed entropy-based second law post-processing analysis to characterize the various thermodynamic losses inside a 3-space solution domain (gas spring+heat exchanger+regenerator) operating under conditions of oscillating pressure and oscillating flow. The 3- space solution domain is adapted from the 2-space solution domain (gas spring+heat exchanger) in Kornhauser's MIT test rig by modifying the heat exchanger space to include a porous regenerator system. A thermal nonequilibrium model which assumes that the regenerator porous matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle is employed. An important and primary objective of this study is the development and application of a thermodynamic loss post-processor to characterize the major thermodynamic losses inside the 3-space model. It is anticipated that the experience gained from thermodynamic loss analysis of the simple 3-space model can be extrapolated to more complex systems like the Stirling engine. It is hoped that successful development of loss post-processors will facilitate the improvement of the optimization capability of Stirling engine analysis codes through better understanding of the heat transfer and power losses. It is also anticipated that the incorporation of a successful thermal nonequilibrium model of the regenerator in Stirling engine CFD analysis codes, will improve our ability to accurately model Stirling regenerators relative to current multidimensional thermal-equilibrium porous media models.

Ebiana, Asuquo B.↗

Slosh-Induced Ullage Collapse Validation by Computational Fluid Dynamics and Application to In-space Propellant Dynamics

Ullage collapse predictions made using Computational Fluid Dynamics (CFD) were validated then made for an in-space cryogenic propellant unsettling application. Ullage collapse is the decrease in propellant tank pressure due to cooling of the ullage gas, which can be significant for cryogenic propulsion systems. Vehicles may be developed that are robust to temperature and pressure changes caused by heat leaks, e.g., propellant tanks with high maximum operation pressures, technologies yielding active thermal conditioning, and systems with ample pressurant for conditioning. These methods are robust because additional mass, cost, and complexity are used to address cryogenic fluid management (CFM). As vehicles are optimized for increased payload capacity, tracking of events such as ullage collapse and self-pressurization in tanks will become increasingly important for propellant inventory calculations. Some ullage collapse events can be studied through ground experimentation, but propellant unsettling in and around in-space coast mission phases must be studied in the absence of an acceleration field. Computational Fluid Dynamics is a powerful tool to address this need, enabling prediction of complex physics prior to flight. Validation of the chosen CFD tool and modeling methodology was done using data from a previously conducted subscale slosh-induced ullage collapse ground test performed in part by the Japan Aerospace Exploration Agency (JAXA). A small tank was partially filled with liquid nitrogen and pressurized with gaseous nitrogen at various temperatures. The tank was oscillated laterally to induce high amplitude slosh dynamics resulting in wave breakup. Successful qualitative validation of slosh wave dynamics and quantitative validation of ullage collapse promoted confidence in extending ullage collapse predictions to a flight application. Mechanisms of ullage collapse caused by main engine cutoff while in-space were identified. Both slosh dynamics leading to engine shutdown and body dynamics following shutdown were shown to contribute to propellant unsettling and ullage collapse.

J M Brodnick↗

A New Interferometer for Monitoring Atmospheric Phase Fluctuations

Water vapor in the Earth's troposphere introduces an extra electrical path in the propagation of radio signals through the atmosphere. The distribution of water vapor is irregular and distorts the wavefronts of incoming radio waves, limiting the angular resolution that can be achieved with ground-based telescopes. The level of fluctuations depends both on the location of the site ,and on the prevailing atmospheric conditions. The ability to measure the fluctuations is therefore important when choosing a site for a new instrument, and for scheduling observations of existing telescopes. Existing phase monitors are radio interferometers that monitor monochromatic beacon tones from geostationary communications satellites at a frequency of about 12 GHz. They have a classical heterodyne design based on two satellite receiving antennas; each has a front-end for amplifying and down-converting the incoming signals using a local oscillator that is phase-locked to a common reference frequency. In addition to multiple phase-locked loops these instruments require expensive phase-stable cabling to reduce the effects of thermal drift. The new system uses two consumer 18" digital satellite TV dishes to monitor satellite TV broadcast signals over a bandwidth of 500 MHz (12.2 to 12.7 GHz). The novel design eliminates the need for phase-locked loops and thermally stable components, and uses a pair of Gilbert Cell multipliers to perform the broadband correlation. A phase monitor has been been built and deployed at the site of the Berkeley-Illinois-Maryland Association Millimeter Array in Northern California, and has been operating successfully since June 1998, measuring the difference in electrical path length for parallel lines of sight to the satellite separated by a baseline of 100 m. With a hardware cost of approximately $4000, it is much cheaper than previous instruments, and the low power requirements and high reliability make the system suitable for site testing in remote locations.

Lay, Oliver↗

Classical discrete time crystals

The spontaneous breaking of time-translation symmetry in periodically driven quantum systems leads to a new phase of matter: the discrete time crystal (DTC). This phase exhibits collective subharmonic oscillations that depend upon an interplay of non-equilibrium driving, many-body interactions and the breakdown of ergodicity. However, subharmonic responses are also a well-known feature of classical dynamical systems ranging from predator–prey models to Faraday waves and a.c.-driven charge density waves. This raises the question of whether these classical phenomena display the same rigidity characteristic of a quantum DTC. In this work, we explore this question in the context of periodically driven Hamiltonian dynamics coupled to a finite-temperature bath, which provides both friction and, crucially, noise. Focusing on one-dimensional chains, where in equilibrium any transition would be forbidden at finite temperature, we provide evidence that the combination of noise and interactions drives a sharp, first-order dynamical phase transition between a discrete time-translation invariant phase and an activated classical discrete time crystal (CDTC) in which time-translation symmetry is broken out to exponentially long timescales. Power-law correlations are present along a first-order line, which terminates at a critical point. We analyse the transition by mapping it to the locked-to-sliding transition of a d.c.-driven charge density wave. To conclude, building upon results from the field of probabilistic cellular automata, we conjecture the existence of classical time crystals with true long-range order, where time-translation symmetry is broken out to infinite times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗