Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theoretical description and modelling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Self-Consistent Ring Current Modeling with Propagating Electromagnetic Ion Cyclotron Waves in the Presence of Heavy Ions

The self-consistent treatment of the RC ion dynamics and EMlC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. Under certain conditions, relativistic electrons, with energies greater than or equal to 1 MeV, can be removed from the outer radiation belt by EMlC wave scattering during a magnetic storm (Summers and Thorne, 2003; Albert, 2003). That is why the modeling of EMlC waves is critical and timely issue in magnetospheric physics. This study will generalize the self-consistent theoretical description of RC ions and EMlC waves that has been developed by Khazanov et al. [2002, 2003] and include the heavy ions and propagation effects of EMlC waves in the global dynamic of self-consistent RC - EMlC waves coupling. The results of our newly developed model that will be presented at Huntsville 2006 meeting, focusing mainly on the dynamic of EMlC waves and comparison of these results with the previous global RC modeling studies devoted to EMlC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

Khazanov, G. V.↗
NASA NTRS2006

Self-Consistent Ring Current/Electromagnetic Ion Cyclotron Waves Modeling

The self-consistent treatment of the RC ion dynamics and EMIC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. For example, the EMlC waves cause the RC decay on a time scale of about one hour or less during the main phase of storms. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt by EMIC wave scattering during a magnetic storm. That is why the modeling of EMIC waves is critical and timely issue in magnetospheric physics. This study will generalize the self-consistent theoretical description of RC ions and EMIC waves that has been developed by Khazanov et al. [2002, 2003] and include the heavy ions and propagation effects of EMIC waves in the global dynamic of self-consistent RC - EMIC waves coupling. The results of our newly developed model that will be presented at the meeting, focusing mainly on the dynamic of EMIC waves and comparison of these results with the previous global RC modeling studies devoted to EMIC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

Khazanov, G. V.↗
NASA NTRS1982

Statistical mechanics of dilute liquid mixtures of He3 in He4

The Woo-Tan-Massey quantum-mechanical theory of dilute mixtures of He-3 in liquid He-4 is extended to incorporate superfluid flow. Then statistical-mechanical principles are employed to find explicit formulas for thermodyanmic and hydrodynamic properties of mixtures, such as internal energy, free energy, entropy, pressure, chemical potential, and normal-fluid and superfluid densities. This paper is the first in a planned series that applies methods originally developed for pure liquid He-4 to dilute mixtures as a means of establishing a comprehensive theory that includes a theoretical bridge between a microscopic description of the mixtures and a hydrodynamic two-fluid model. Ultimately effects of applied electric and magnetic fields on dilute solutions are also planned to be treated.

Jackson, H. W.↗
NASA NTRS1999

Ultrasonic Nondestructive Characterization of Adhesive Bonds

Adhesives and adhesive joints are widely used in various industrial applications to reduce weight and costs, and to increase reliability. For example, advances in aerospace technology have been made possible, in part, through the use of lightweight materials and weight-saving structural designs. Joints, in particular, have been and continue to be areas in which weight can be trimmed from an airframe through the use of novel attachment techniques. In order to save weight over traditional riveted designs, to avoid the introduction of stress concentrations associated with rivet holes, and to take full advantage of advanced composite materials, engineers and designers have been specifying an ever-increasing number of adhesively bonded joints for use on airframes. Nondestructive characterization for quality control and remaining life prediction has been a key enabling technology for the effective use of adhesive joints. Conventional linear ultrasonic techniques generally can only detect flaws (delamination, cracks, voids, etc) in the joint assembly. However, more important to structural reliability is the bond strength. Although strength, in principle, cannot be measured nondestructively, a slight change in material nonlinearity may indicate the onset of failure. Furthermore, microstructural variations due to aging or under-curing may also cause changes in the third order elastic constants, which are related to the ultrasonic nonlinear parameter of the polymer adhesive. It is therefore reasonable to anticipate a correlation between changes in the ultrasonic nonlinear acoustic parameter and the remaining bond strength. It has been observed that higher harmonics of the fundamental frequency are generated when an ultrasonic wave passes through a nonlinear material. It seems that such nonlinearity can be effectively used to characterize bond strength. Several theories have been developed to model this nonlinear effect. Based on a microscopic description of the nonlinear interface binding force, a quantitative method was presented. Recently, a comparison between the experimental and simulated results based on a similar theoretical model was presented. A through-transmission setup for water immersion mode-converted shear waves was used to analyze the ultrasonic nonlinear parameter of an adhesive bond. In addition, ultrasonic guided waves have been used to analyze adhesive or diffusion bonded joints. In this paper, the ultrasonic nonlinear parameter is used to characterize the curing state of a polymer/aluminum adhesive joint. Ultrasonic through-transmission tests were conducted on samples cured under various conditions. The magnitude of the second order harmonic was measured and the corresponding ultrasonic nonlinear parameter was evaluated. A fairly good correlation between the curing condition and the nonlinear parameter is observed. The results show that the nonlinear parameter might be used as a good indicator of the cure state for adhesive joints.

Qu, Jianmin↗
NASA NTRS2015

Thermodynamic Derivation of the Activation Energy for Ice Nucleation

Cirrus clouds play a key role in the radiative and hydrological balance of the upper troposphere. Their correct representation in atmospheric models requires an understanding of the microscopic processes leading to ice nucleation. A key parameter in the theoretical description of ice nucleation is the activation energy, which controls the flux of water molecules from the bulk of the liquid to the solid during the early stages of ice formation. In most studies it is estimated by direct association with the bulk properties of water, typically viscosity and self-diffusivity. As the environment in the ice-liquid interface may differ from that of the bulk, this approach may introduce bias in calculated nucleation rates. In this work a theoretical model is proposed to describe the transfer of water molecules across the ice-liquid interface. Within this framework the activation energy naturally emerges from the combination of the energy required to break hydrogen bonds in the liquid, i.e., the bulk diffusion process, and the work dissipated from the molecular rearrangement of water molecules within the ice-liquid interface. The new expression is introduced into a generalized form of classical nucleation theory. Even though no nucleation rate measurements are used to fit any of the parameters of the theory the predicted nucleation rate is in good agreement with experimental results, even at temperature as low as 190 K, where it tends to be underestimated by most models. It is shown that the activation energy has a strong dependency on temperature and a weak dependency on water activity. Such dependencies are masked by thermodynamic effects at temperatures typical of homogeneous freezing of cloud droplets; however, they may affect the formation of ice in haze aerosol particles. The new model provides an independent estimation of the activation energy and the homogeneous ice nucleation rate, and it may help to improve the interpretation of experimental results and the development of parameterizations for cloud formation.

Cirrus clouds↗
NASA NTRS1999

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗
NASA NTRS1993

Mathematical description of complex chemical kinetics and application to CFD modeling codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗
NASA NTRS1993

Mathematical Description of Complex Chemical Kinetics and Application to CFD Modeling Codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗
DOE OSTI2025

Measurement of the 28 Si ⁢(𝑛,𝑛′⁢𝛾) cross section with 𝑛, 𝛾, and correlated 𝑛−𝛾 angular distributions

Silicon has become an unavoidable element in the circuitry central to everyday life. In turn, the interactions of silicon isotopes with neutrons for nuclear physics applications, among other motivations, have become increasingly important to understand. The dominant isotope of silicon, 28 Si, is thus of primary interest for enhanced understanding for neutron transport calculations and related investigations. Unfortunately, the existing measurement database for neutron scattering reactions on 28 Si is minimal, and nuclear data evaluations on this topic have not been updated for decades. This article details new measurements of the 28 Si ⁢(𝑛,𝑛′⁢𝛾) reaction utilizing multiple analysis methods available within the correlated gamma neutron array for scattering (CoGNAC). Specifically, high-precision near-threshold results and high-incident-energy results were obtained using the 𝛾-only and correlated 𝑛−𝛾 techniques. First-ever measurements of the correlated 𝑛−𝛾 angular distribution for particles emitted following population of the first excited state in 28 Si were obtained as well, which provide unique insight into theoretical descriptions of the inelastic neutron scattering reaction mechanism itself and detailed guidance for nuclear reaction models. The results agree well with literature data where they exist, and substantially expand on the current database for neutron reactions on 28 Si .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

Magnetospheric plasma modeling (0-100 keV)

Spacecraft surface charging, which is primarily a current balance phenomenon, is in general a function of the dominant currents to and from the vehicle's surface. Within the near-earth magnetosphere the dominant currents to the surface are the ambient space plasma fluxes between approximately 0 and 100 keV. A major effort to understand the near-earth environment was initiated when spacecraft charging became a major issue. The present paper has the objective to summarize the basic features of the models which have resulted from this effort. A description is given of four categories of models, based primarily on the degree of empirical and theoretical input. Types of quantitative models are discussed, taking into account definitions, statistical models, analytic models, static models, and time-dependent models. Engineering models are also considered, giving attention to baseline models and 'worst-case' models.

Garrett, H. B.↗
NASA NTRS1989

Planetary rings - Theory

Theoretical models of planetary-ring dynamics are examined in a brief analytical review. The mathematical description of streamlines and streamline interactions is outlined; the redistribution of angular momentum due to collisions between particles is explained; and problems in the modeling of broad, narrow, and arc rings are discussed.

Borderies, Nicole↗
DOE OSTI2023

𝛽 decay of 36 Mg and 36 Al: Identification of a 𝛽-decaying isomer in 36 Al

The level structure of 36 Al has been studied via 𝛽 decay of 36 Mg at the Facility for Rare Isotope Beams (FRIB) and the National Superconducting Cyclotron Laboratory (NSCL). A long-lived isomer in 36 Al was identified which decays by 𝛽 to an excited state of 36 Si. The ground state and the isomeric state of 36 Al were found to populate different energy levels of 36 Si. Furthermore, the results from the two data sets in the present work complement each other. Configuration interaction calculations performed with the FSU shell-model Hamiltonians provide reasonable descriptions to the experimental observations and offer insight into future improvements of the theoretical interpretation.

20 ≤ A ≤ 38↗
NASA NTRS2011

Effective Albedo of Vegetated Terrain at L-Band

This paper derives an explicit expression for an effective albedo of vegetated terrain from the zero- and multiple- order radiative transfer (RT) model comparison. The formulation establishes a direct physical link between the effective vegetation parameterization and the theoretical description of absorption and scattering within the canopy. The paper will present an evaluation of the derived albedo for corn canopies with data taken during an experiment at Alabama A&M Winfield A. Thomas Agricultural Research Station near Huntsville, Alabama in June, 1998. The test site consisted of two 50-m x 60-m plots - one with a bare surface and the other with grass cover - and four 30-m x 50-m plots of corn at different planting densities. One corn field was planted at a full density of 9.5 plants/sq m while the others were planted at 1/3, 1/2 and 2/3 of the full density. The fields were observed with a truck-mounted L-band radiometer at incident angle of 15 degree for the period of two weeks. Soil moisture (SM) changed daily due to irrigation and natural rainfall. Variations in gravimetric SM from 18 % to 34 % were seen during this period. Ground truth data, including careful characterization of the corn size and orientation statistics, and its dielectric, was also collected and used to simulate the effective albedo for the vegetation. The single-scattering albedo is defined as the fractional power scattered from individual vegetation constituents with respect to canopy extinction. It represents single-scattering properties of vegetation elements only, and is independent of ground properties. The values of the albedo get higher when there is dense vegetation (i.e. forest, mature corn, etc.) with scatterers, such as branches and trunks (or stalks in the case of corn), which are large with respect to the wavelength. This large albedo leads to a reduction in brightness temperature in the zero-order RT solution (known as tau-omega model). Higher-order multiple-scattering RT solutions are required for proper representation of scattering within vegetation. In this paper, an expression for an effective albedo for the whole canopy including the ground is derived for use in the zero-order RT model-based SM retrieval. This effective albedo takes into account of all the processes taking place within the canopy, including multiple-scattering. This new formulation will be presented and its importance for microwave SM retrieval will be evaluated for corn canopies in conjunction with the detailed ground truth data obtained during the experiment at Alabama in 1998. Emphasis will be placed on examining how the radiometer response to SM is modified by the corn canopy scattering under different field conditions. A semi-empirical parameterization of the effective albedo will be investigated through analysis of SM and vegetation water content effects on the effective albedo.

Kurum, Mehmet↗
NASA NTRS2006

Self-Consistent Ring Current Modeling with Propagating Electromagnetic Ion Cyclotron Waves in the Presence of Heavy Ions

The self-consistent treatment of the RC ion dynamics and EMIC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt by EMIC wave scattering during a magnetic storm. That is why the modeling of EMIC waves is critical and timely issue in magnetospheric physics. To describe the RC evolution itself this study uses the ring current-atmosphere interaction model (RAM). RAM solves the gyration and bounce-averaged Boltzmann-Landau equation inside of geosynchronous orbit. Originally developed at the University of Michigan, there are now several branches of this model currently in use as describe by Liemohn namely those at NASA Goddard Space Flight Center This study will generalize the self-consistent theoretical description of RC ions and EMIC waves that has been developed by Khazanov and include the heavy ions and propagation effects of EMIC waves in the global dynamic of self-consistent RC - EMIC waves coupling. The results of our newly developed model that will be presented at GEM meeting, focusing mainly on the dynamic of EMIC waves and comparison of these results with the previous global RC modeling studies devoted to EMIC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

Khazanov, George V.↗
NASA NTRS1985

Tone excited jets. V - A theoretical model and comparison with experiment

A mathematical model of tone-excited jets is developed which consists of two major components: a mathematical description of the process by which the intrinsic instability waves of the jet are excited by the upstream tones; and the modeling of the nonlinear interaction between the mean flow of the jet, the excited large-scale instability waves or turbulence structure, and the fine scale turbulence. It is assumed that each of these jet flow components can be characterized by a few parameters, which are then related by conservation equations that are supplemented by closure models. This quasi-linear model's results are compared with experimental measurements, and good agreement is obtained over a wide range of excitation frequencies and excitation levels.

Tam, C. K. W.↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Counterions under a Surface-Adsorbed Cationic Surfactant Monolayer: Structure and Thermodynamics

The surface adsorption of ionic surfactants is fundamental for many widespread phenomena in life sciences and for a wide range of technological applications. However, direct atomic-resolution structural experimental studies of noncrystalline surface-adsorbed films are scarce. Thus, even the most central physical aspects of these films, such as their charge density, remain uncertain. Consequently, theoretical models based on contradicting assumptions as for the surface films' ionization are widely used for the description and prediction of surface thermodynamics. We employ X-ray reflectivity to obtain the angstrom Ångström-scale surface-normal structure of surface-adsorbed films of the cationic surfactant cetyltrimethylammonium bromide (CTAB) in aqueous solutions at several different temperatures and concentrations. In conjunction with published neutron reflectivity data, we determine the surface-normal charge distribution due to the dissociated surfactants' headgroups. The distribution appears to be inconsistent with the Gouy-Chapman model yet consistent with a compact Stern layer model of condensed counterions. Here, the experimental surfactant adsorption thermodynamics conforms well to classical, Langmuir and Kralchevsky, adsorption models. Furthermore, the Kralchevsky model correctly reproduces the observed condensation of counterions, allowing the values of the adsorption parameters to be resolved, based on the combination of the present data and the published surface tension measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

String Phase in an Artificial Spin Ice

Abstract One-dimensional strings of local excitations are a fascinating feature of the physical behavior of strongly correlated topological quantum matter. Here we study strings of local excitations in a classical system of interacting nanomagnets, the Santa Fe Ice geometry of artificial spin ice. We measured the moment configuration of the nanomagnets, both after annealing near the ferromagnetic Curie point and in a thermally dynamic state. While the Santa Fe Ice lattice structure is complex, we demonstrate that its disordered magnetic state is naturally described within a framework of emergent strings. We show experimentally that the string length follows a simple Boltzmann distribution with an energy scale that is associated with the system’s magnetic interactions and is consistent with theoretical predictions. The results demonstrate that string descriptions and associated topological characteristics are not unique to quantum models but can also provide a simplifying description of complex classical systems with non-trivial frustration.

36 MATERIALS SCIENCE↗