Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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,333 records · Page 74

Modelling time variable and total eclipses of the millisecond pulsar PSR 1744-24A

The 11 ms pulsar PSR 1744-24A in the globular clusters Terzan 5 is the first and only known example of a binary eclipsing pulsar showing time-variable eclipses with irregular and nonsymmetric eclipse patterns. Occasionally, the pulsar eclipse becomes total, and no radio emission is detected for an extended period of time. We present here numerical results of a hydrodynamical model aimed at explaining the remarkable eclipse properties of PSR 1744-24A. We show that the hydrodynamic time scale for the partial or total pulsar enshrouding is of order of the orbital period for a relatively large mass loss rate, in agreement with observations. The time-variable eclipse egress positions imply that the mass loss rate from the companion varies within one orbital period. Hydrodynamic modelling of the mass outflow of PSR 1744-24A is of crucial importance for the understanding of a class of millisecond pulsars whose radio emission might be occasionally or permanently quenched by mass outflows from their companions. The properties of PSR 1744-24A suggest a new way of studying millisecond pulsars involving radio, X-ray, gamma-ray, and optical observations.

Tavani, Marco↗

Application of solar max ACRIM data to analyze solar-driven climatic variability on Earth

Langley, in the late nineteenth century, attempted to measure solar irradiance over an extended period of time in order to detect changes. The problem with this and other early attempts was that ground based measurements are not sufficiently accurate to measure solar irradiance fluctuations, which are on the order of 0.1%. It was not until the Active Cavity Radiometer Irradiance Monitor (ACRIM) experiment on the NASA Solar Maximum Mission (SMM) was launched in 1980 that continuous data with precision 0.1% was available.

Hoffert, M. I.↗

Spectroscopic classification, variability, and SED of the Fermi -detected CSS 3C 286: the radio-loudest NLS1 galaxy?

ABSTRACT 3C 286 is a well-known calibrator source in radio astronomy. It is also one of very few compact steep-spectrum (CSS) sources detected in γ-rays. Here, we perform a detailed spectroscopic and variability analysis and present the first quasi-simultaneous optical to X-ray spectral energy distribution (SED) in order to reveal physical mechanisms that dominate its emission at different wavelengths, and arrive at a reliable optical source classification. The first main result of our study reveals several pitfalls when applying simple broad- or narrow-line Seyfert 1 (BLS1 or NLS1) classification criteria that only look at the [O iii]–Hβ complex. [O iii] and Hβ can be dominated by the same outflow components, in which case full width at half-maximum (Hβ) is no reliable classification criterion, and extinction by intrinsic or intervening material can make the highest velocity Hβ component undetectable. After careful combination of all information from ultraviolet (UV)–optical spectra along with multiwavelength data, we confirm that 3C 286 can be classified as NLS1 galaxy, with line properties and supermassive black hole mass (of order 108 M⊙ and accreting near the Eddington limit) close to the BLS1 regime, making it an important borderline object. The quasi-simultaneous SED taken with Swift shows a sharp rise in the optical–UV, implying the presence of a strong accretion-disc component with extreme UV excess, consistent with emission-line diagnostics. Finally, we report the discovery of X-ray variability of 3C 286, plausibly dominated by jet emission, and variable by at least a factor of ∼4. This result suggests to exercise caution when using 3C 286 as radio calibrator in high-resolution radio very long baseline interferometry observations.

Yao, Su↗

First-Order Frameworks for Managing Models in Engineering Optimization

Approximation/model management optimization (AMMO) is a rigorous methodology for attaining solutions of high-fidelity optimization problems with minimal expense in high- fidelity function and derivative evaluation. First-order AMMO frameworks allow for a wide variety of models and underlying optimization algorithms. Recent demonstrations with aerodynamic optimization achieved three-fold savings in terms of high- fidelity function and derivative evaluation in the case of variable-resolution models and five-fold savings in the case of variable-fidelity physics models. The savings are problem dependent but certain trends are beginning to emerge. We give an overview of the first-order frameworks, current computational results, and an idea of the scope of the first-order framework applicability.

Alexandrov, Natlia M.↗

An Assessment of Global Organic Carbon Flux Along Continental Margins

This project was designed to use real-time and historical SeaWiFS and AVHRR data, and real-time MODIS data in order to estimate the global vertical carbon flux along continental margins. This required construction of an empirical model relating surface ocean color and physical variables like temperature and wind to vertical settling flux at sites co-located with sediment trap observations (Santa Barbara Basin, Cariaco Basin, Gulf of California, Hawaii, and Bermuda, etc), and application of the model to imagery in order to obtain spatially-weighted estimates.

Thunell, Robert↗

Experiment Design and Visualization Techniques for an X-59 Low-boom Variability Study

This presentation outlines the design of experiments approach and data visualization techniques for a simulation study of sonic booms from NASA’s X-59 supersonic aircraft. The X-59 will soon be flown over communities across the contiguous USA as it produces a low-loudness sonic boom, or low-boom. Survey data on human perception of low-booms will be collected to support development of potential future commercial supersonic aircraft noise regulatory standards. The macroscopic atmosphere plays a critical role in the loudness of sonic booms. The extensive sonic boom simulation study presented herein was completed to assess climatological, geographical, and seasonal effects on the variability of the X-59’s low-boom loudness and noise exposure region size in order to inform X-59 community test planning. The loudness and extent of the noise exposure region make up the “sonic boom carpet.” Two spatial and temporal resolutions of atmospheric input data to the simulation were investigated. A Fast Flexible Space-Filling Design was used to select the locations across the USA for the two spatial resolutions. Analysis of simulated X-59 low-boom loudness data within a regional subset of the northeast USA was completed using a bootstrap forest to determine the final spatial and temporal resolution of the countrywide simulation study. Atmospheric profiles from NOAA’s Climate Forecast System Version 2 database were used to generate over one million simulated X-59 carpets at the final selected 138 locations across the USA. Effects of aircraft heading, season, geography, and climate zone on low-boom levels and noise exposure region size were analyzed. Models were developed to estimate loudness metrics throughout the USA for X-59 supersonic cruise overflight, and results were visualized on maps to show geographical and seasonal trends. These results inform regulators and mission planners on expected variations in boom levels and carpet extent from atmospheric variations. Understanding potential carpet variability is important when planning community noise surveys using the X-59.

X-59↗

Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)

97 MATHEMATICS AND COMPUTING↗

Placement and Sizing of Inverter-Based Renewable Systems in Multi-Phase Distribution Networks

This study develops a tractable formulation for optimal placement and sizing of inverter-based renewable systems in multi-phase distribution networks. The goal of the formulation is to minimize the cost of inverter installation, average power import, and average distributed generation curtailment. Threephase and single-phase inverter models are presented that preserve the underlying mappings between renewable uncertainty to power injection. The uncertainty of distributed generators (DGs) and loads are characterized by a finite set of scenarios. Linear multi-phase power flow approximations are used in conjunction with scenario reduction techniques to arrive at a tractable twostage stochastic formulation for optimal DG placement and sizing. First-stage decisions are locations for DG deployment and capacity sizes, and second-stage decisions include DG real power curtailment, reactive power support, as well as feeder voltage profile. The resulting formulation is a mixed-integer second-order cone program and can be solved efficiently either by existing optimization solvers or by relaxing the binary variables to the [0,1] interval. Simulation studies on standard multi-phase IEEE test feeders promise that optimal stochastic planning of DGs reduces costs during validation, compared to a scheme where uncertainty is only represented by its average value.

24 POWER TRANSMISSION AND DISTRIBUTION↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

Role of Satellite Rainfall Information in Improving Understanding of the Dynamical Link Between the Tropics and Extratropics Prospects of Improved Forecasts of Weather and Short-Term Climate Variability on Sub-Seasonal Time Scales

The tropics and extratropics are two dynamically distinct regimes. The coupling between these two regimes often defies simple analytical treatment. Progress in understanding of the dynamical interaction between the tropics and extratropics relies on better observational descriptions to guide theoretical development. However, global analyses currently contain significant errors in primary hydrological variables such as precipitation, evaporation, moisture, and clouds, especially in the tropics. Tropical analyses have been shown to be sensitive to parameterized precipitation processes, which are less than perfect, leading to order-one discrepancies between estimates produced by different data assimilation systems. One strategy for improvement is to assimilate rainfall observations to constrain the analysis and reduce uncertainties in variables physically linked to precipitation. At the Data Assimilation Office at the NASA Goddard Space Flight Center, we have been exploring the use of tropical rain rates derived from the TRMM Microwave Imager (TMI) and the Special Sensor Microwave/ Imager (SSM/I) instruments in global data assimilation. Results show that assimilating these data improves not only rainfall and moisture fields but also related climate parameters such as clouds and radiation, as well as the large-scale circulation and short-range forecasts. These studies suggest that assimilation of microwave rainfall observations from space has the potential to significantly improve the quality of 4-D assimilated datasets for climate investigations (Hou et al. 2001). In the next few years, there will be a gradual increase in microwave rain products available from operational and research satellites, culminating to a target constellation of 9 satellites to provide global rain measurements every 3 hours with the proposed Global Precipitation Measurement (GPM) mission in 2007. Continued improvements in assimilation methodology, rainfall error estimates, and model parameterizations are needed to ensure that we derive maximum benefits from these observations.

Hou, Arthur Y.↗

Soil carbon cycling proxies: Understanding their critical role in predicting climate change feedbacks

The complexity of processes and interactions that drive soil C dynamics necessitate the use of proxy variables to represent soil characteristics that cannot be directly measured (correlative proxies), or that aggregate information about multiple soil characteristics into one variable (integrative proxies). These proxies have proven useful for understanding the soil C cycle, which is highly variable in both space and time, and are now being used to make predictions of the C fate and persistence under future climate scenarios. As these proxies are used at increasingly larger scales, the C pools and processes that proxies represent must be thoughtfully considered in order to minimize uncertainties in empirical understanding, as well as in model parameters and in model outcomes. The importance of these uncertainties is further amplified by the current need to make predictions of the C cycle for the non steady state environmental conditions resulting from global climate change. To clarify the appropriate uses of proxy variables, we provide specific examples of proxy variables that could improve decision making, adaptation choices, and modeling skill, while not foreclosing on – and also encouraging – continued work on their mechanistic underpinnings. We explore the use of three common soil proxies used to study soil organic matter: metabolic quotient, clay content, and physical fractionation. We also consider emerging data types, specifically genome-sequence data, and how these serve as proxies for microbial community activities. We opine that the demand for increasing mechanistic detail, and the flood of data from new imaging and genetic techniques, does not replace the value of correlative and integrative proxies--variables that are simpler, easier, or cheaper to measure. By closely examining the current knowledge gaps and broad assumptions in soil C cycling with the proxies already in use, we can develop new hypotheses and specify criteria for new and needed proxies.

Bailey, Vanessa L.↗

Aerogel-derived nickel-iron oxide catalysts for oxygen evolution reaction in alkaline media

Anion exchange membrane water electrolyzers (AEMWEs) can generate hydrogen with a pure water feed using noble metal-free catalysts. Here, the development of highly active and stable catalysts for oxygen evolution reaction (OER) is required for improving performance of AEMWEs systems. Ni-Fe (oxy)hydroxides show high OER catalytic activity in alkaline media, but typically have low surface area. In this work, we investigate a series of Ni-Fe oxides with high surface area and disordered morphology, obtained using an aerogel synthesis method. We evaluate the impact of different synthesis variables on the OER activity and demonstrate that heat treatment at high temperatures generates more ordered structure, resulting in a decrease in OER activity. Advanced characterization reveals that maintaining highly disordered and porous structure of the aerogel is essential to achieving high OER activity, as it enables the formation of highly OER-active lamellar structures of the catalyst.

36 MATERIALS SCIENCE↗

A velocity/stress mixed stabilized nodal finite element for elastodynamics: Analysis and computations with strongly and weakly enforced boundary conditions

A new nodal mixed finite element is proposed for the simulation of linear elastodynamics and wave propagation problems in time domain. Our method is based on equal-order interpolation discrete spaces for both the velocity (or displacement) and stress (or strain) tensor variables. The mixed form is derived using either the velocity/stress or velocity/strain pair of unknowns, the latter being instrumental in extensions of the method to nonlinear mechanics. The proposed approach works equally well on hexahedral or tetrahedral grids and, for this reason, it is suitable for time-domain engineering applications in complex geometry. The peculiarity of the proposed approach is the use of the rate form of the stress update equation, which yields a set of governing equations with the structure of a non-dissipative space/time Friedrichs’ system. We complement standard traction boundary conditions for the stress with strongly and weakly enforced boundary conditions for the velocity (or displacement). Weakly enforced boundary conditions are particularly suitable when considering complex geometrical shapes, because they do not require dedicated data structures for the imposition of the boundary degrees of freedom, but, rather, they utilize the structure of the variational formulation. Here we also show how the framework of weakly enforced boundary conditions can be used to develop variational forms for multi-domain simulations of heterogeneous media. A complete analysis including stability and convergence proofs is included, in the case of a space–time variational approach. A series of computational tests are used to demonstrate and verify the performance of the proposed approach.

42 ENGINEERING↗

A review of smart building sensing system for better indoor environment control

Here this paper aims to provide a systemic review of how indoor sensors influence in managing optimal energy saving, thermal comfort, visual comfort, and indoor air quality in the built environment. The optimal management of energy saving and occupant comfort plays a vital role in the built environment because the occupant's productivity and health are highly influenced by indoor environmental quality. In order to do this, there must be a functional sensing system that connects the environment variables (e.g., temperature) with building environmental control systems such as the heating, ventilation, and air-conditioning system. This paper starts with an overview of the importance of energy saving and occupant comfort in the built environment. It then discusses sensors and their importance in the built environment and reviews the different types of sensors, which explains them in terms of how they influence the indoor built environment and occupant productivity. The paper further explores the application of sensors in the built environment and analyzes this in terms of energy saving, thermal comfort, visual comfort, and indoor air quality. Following this, the data analysis is discussed in terms of data, information, and knowledge accrued from the sensors. Lastly, the paper discusses the future challenges for the improvement of building indoor environmental quality and energy saving by the application of sensors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Scalar conservation and boundedness in simulations of compressible flow

With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.

97 MATHEMATICS AND COMPUTING↗

Estimated capital costs of fish exclusion technologies for hydropower facilities

Hydropower is a reliable source of renewable energy, and its future expansion is likely to be in the form of either smaller new stream development (NSD) projects or powering existing non-powered dams. Thresholds for entrainment risk to fish and the requirements for fish exclusion at hydropower facilities often differ depending on the species involved, the characteristics of the facility, and the goals of stakeholders, but little quantitative information is present within the literature regarding the specific costs of fish exclusion measures. Cost data associated with protection, mitigation, and enhancement (PM&E) measures related to positive barrier screening were identified using keyword searches of an existing environmental mitigation cost data set and manual extraction from regulatory licensing documents available in the Federal Energy Regulatory Commission (FERC) eLibrary. This approach yielded a total of 50 p.m.&E mitigation measures with estimated capital construction costs pertaining to positive barrier screens and represented <10% of the 171 total FERC project dockets available in the data set. These data were highly skewed toward conventional relicensing projects, as <7% were associated with NSD projects. Results indicate highly variable costs are associated with fish screening, with flow-normalized costs one to two orders of magnitude higher for screening with the highest exclusion capability (≤0.09 in. spacing) compared with coarser screening (1–2 in.). These data provide an initial baseline for estimating exclusion costs for hydropower development and may help developers consider options for more fish-friendly generation technologies, though gaps remain relating to a lack of data, particularly for NSD projects.

13 HYDRO ENERGY↗

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. Here, these observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of uniaxial single-ion anisotropy on the magnetic and thermal properties of Heisenberg antiferromagnets within unified molecular field theory

Here, the influence of uniaxial single-ion anisotropy –DS 2 z on the magnetic and thermal properties of Heisenberg antiferromagnets (AFMs) is investigated. The uniaxial anisotropy is treated exactly and the Heisenberg interactions are treated within unified molecular field theory (MFT), where thermodynamic variables are expressed in terms of directly measurable parameters. The properties of collinear AFMs with ordering along the z axis (D>0) in applied field H z = 0 are calculated versus D and temperature T, including the ordered moment μ, the Néel temperature T N , the magnetic entropy, internal energy, heat capacity, and the anisotropic magnetic susceptibilities χ ∥ and χ ⊥ in the paramagnetic (PM) and AFM states. The high-field average magnetization per spin μ z (H z ,D,T) is found, and the critical field H c (D,T) is derived at which the second-order AFM to PM phase transition occurs. The magnetic properties of the spin-flop (SF) phase are calculated, including the zero-field properties T N (D) and μ(D,T). The high-field μ z (H z ,D,T) is determined, together with the associated spin-flop field H SF (D,T) at which a second-order SF to PM phase transition occurs. The free energies of the AFM, SF, and PM phases are derived from which H z –T phase diagrams are constructed. For f J =–1 and –0.75, where f J = θ pJ /T NJ and θ pJ and T NJ are the Weiss temperature in the Curie-Weiss law and the Néel temperature due to exchange interactions alone, respectively, phase diagrams in the H z –T plane similar to previous results are obtained. However, for f J = 0 we find a topologically different phase diagram where a spin-flop bubble with PM and AFM boundaries occurs at finite H z and T. Also calculated are properties arising from a perpendicular magnetic field, including the perpendicular susceptibility χ ⊥ (D,T), the associated effective torque at low fields arising from the –DS 2 z term in the Hamiltonian, the high-field perpendicular magnetization μ ⊥ , and the perpendicular critical field H c⊥ at which the second-order AFM to PM phase transition occurs. In addition to the above results for D > 0, the T N (D) and ordered moment μ(T,D) for collinear AFM ordering along the x axis with D < 0 are determined. In order to compare the properties of the above spin systems with those of noninteracting systems with –DS 2 z uniaxial anisotropy with either sign of D, Supplemental Material is provided in which results for the thermal and magnetic properties of such noninteracting spin systems are given.

36 MATERIALS SCIENCE↗