SEARCH · Engineering Papers

Results for “variational principles”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Evaluation of hydrograph separation techniques with uncertain end-member composition

Hydrograph separation is one of many approaches used to analyse shifts in source water contributions to stream flow resulting from climate change in remote watersheds. Understanding these shifts is vital, as shifts in source water contributions to a stream can shape water management decisions. Because remote watersheds are often inaccessible and have poorly characterized contributing water sources, or end-members, it is critical to understand the implications of using different hydrograph separation techniques in these data-limited environments. To explore the uncertainty associated with different techniques, results from two hydrograph separation techniques, mass balance and principle component analysis, were compared using 3 years of aqueous geochemical data from the East River watershed located in the Elk Mountains of Central Colorado. Solute concentrations of the end-members were characterized by both a limited set of direct chemical measurements of different sources and detailed seasonal instream chemistry to examine the influences of uncertain end-member compositions in a data-limited environment. Annual volumetric end-member contributions to stream flow had relatively good agreement across separation techniques. Large variations in time were observed in the hydrograph separations, depending on the end-member type, and estimated flow contributions varied between the selected solutes. End-member concentrations characterized by stream chemistry showed several limitations including a reduced number of distinguishable end-members and differences in timing of flow contributions. Here the results highlight the benefits of using multiple hydrograph separation techniques by providing a ‘weight-of-evidence’ approach to environments with limited end-member concentration data.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Ferromagnetism with in-plane magnetization, Dirac spin-gapless semiconducting properties, and tunable topological states in two-dimensional rare-earth metal dinitrides

Since the successful synthesis of bulk single crystals MoN 2 and ReN 2 , which have a layered structure, transition-metal dinitrides have attracted considerable attention in recent years. Here, we focus on rare-earth metal (Rem) elements, and propose seven stable Rem dinitride monolayers with a 1T structure, namely, 1T-RemN 2 . We use first-principles calculations, and find that these monolayers have a ferromagnetic ground state with in-plane magnetization. Without spin-orbit coupling (SOC), the band structures are spin-polarized with Dirac points at the Fermi level. Remarkably, the 1T-LuN 2 monolayer exhibits an isotropic magnetocrystalline anisotropy energy in the xy plane with in-plane magnetization, indicating easy tunability of the magnetization direction. Additionally, when rotating the magnetization vector in the xy plane, we propose a model that accurately describes the variation of the SOC band gap and the two possible topological states (Weyl-like semimetal and Chern insulator states) whose properties are tunable. The Weyl-like semimetal state is a critical point between the two Chern insulator states with opposite sign of the Chern numbers (±1). The nontrivial band gap (up to 60.3 meV) and the Weyl-like semimetal state are promising for applications in spintronic devices.

2-dimensional systems↗

NASA NTRS1993

Space Shuttle Main Engine performance analysis

For a number of years, NASA has relied primarily upon periodically updated versions of Rocketdyne's power balance model (PBM) to provide space shuttle main engine (SSME) steady-state performance prediction. A recent computational study indicated that PBM predictions do not satisfy fundamental energy conservation principles. More recently, SSME test results provided by the Technology Test Bed (TTB) program have indicated significant discrepancies between PBM flow and temperature predictions and TTB observations. Results of these investigations have diminished confidence in the predictions provided by PBM, and motivated the development of new computational tools for supporting SSME performance analysis. A multivariate least squares regression algorithm was developed and implemented during this effort in order to efficiently characterize TTB data. This procedure, called the 'gains model,' was used to approximate the variation of SSME performance parameters such as flow rate, pressure, temperature, speed, and assorted hardware characteristics in terms of six assumed independent influences. These six influences were engine power level, mixture ratio, fuel inlet pressure and temperature, and oxidizer inlet pressure and temperature. A BFGS optimization algorithm provided the base procedure for determining regression coefficients for both linear and full quadratic approximations of parameter variation. Statistical information relative to data deviation from regression derived relations was also computed. A new strategy for integrating test data with theoretical performance prediction was also investigated. The current integration procedure employed by PBM treats test data as pristine and adjusts hardware characteristics in a heuristic manner to achieve engine balance. Within PBM, this integration procedure is called 'data reduction.' By contrast, the new data integration procedure, termed 'reconciliation,' uses mathematical optimization techniques, and requires both measurement and balance uncertainty estimates. The reconciler attempts to select operational parameters that minimize the difference between theoretical prediction and observation. Selected values are further constrained to fall within measurement uncertainty limits and to satisfy fundamental physical relations (mass conservation, energy conservation, pressure drop relations, etc.) within uncertainty estimates for all SSME subsystems. The parameter selection problem described above is a traditional nonlinear programming problem. The reconciler employs a mixed penalty method to determine optimum values of SSME operating parameters associated with this problem formulation.

Santi, L. Michael↗

NASA NTRS2002

The Influence of Tropospheric Processes in Modeling the Middle Atmosphere with Gravity Waves

Our Numerical Spectral Model (NSM) extends from the ground up into the thermosphere and has a vertical grid point resolution of about 0.5 km to resolve the interactions of gravity waves (GWs) described with Hines' Doppler Spread Parameterization (DSP). This model produces in the stratosphere and mesosphere the major features of QBO, SAO, tides, and planetary waves. The purpose of this paper is to discuss results from an initial study with our 3D model that shows how certain tropospheric processes can affect the dynamics of the middle atmosphere. Under the influence of tropospheric heating, and augmented by GW interactions, two distinct but related processes can be identified. (1) A meridional circulation develops in the stratosphere, with rising motions at low latitudes that are in magnitude comparable to the downward propagation of the QBO. As Dunkerton pointed out, a larger GW source is then required to reproduce the observed QBO, which tends to move us closer to the values recommended for the DSP. This has significant consequences for our model results that describe the upper mesosphere, considering the general importance of GWs for this region and in influencing planetary waves (e.g., 2-day wave) and tides in particular. (2) Tropospheric heating produces zonal jets near the tropopause that are related to latitudinal variations in pressure and reversing temperature variations (resembling the dynamical conditions near the mesopause), which in turn is conducive to generate baroclinic instability. Modeling results show that our ability to generate the QBO critically depends on the magnitude of the temperature reversal that is a measure of this instability. Planetary waves are generated in this process, which can apparently interfere with or augment the GW interactions. As originally demonstrated by Lindzen and Holton, the eastward propagating Kelvin waves and westward propagating Rossby gravity waves (generated by tropospheric convection) can in principle provide the acceleration to influence the QBO, and we were able to confirm this with our 3D model.

Mayr, H. G.↗

NASA NTRS2004

Generalization of Equivalent Crystal Theory to Include Angular Dependence

In the original Equivalent Crystal Theory, each atomic site in the real crystal is assigned an equivalent lattice constant, in general different from the ground state one. This parameter corresponds to a local compression or expansion of the lattice. The basic method considers these volumetric transformations and, in addition, introduces the possibility that the reference lattice is anisotropically distorted. These distortions however, were introduced ad-hoc. In this work, we generalize the original Equivalent Crystal Theory by systematically introducing site-dependent directional distortions of the lattice, whose corresponding distortions account for the dependence of the energy on anisotropic local density variations. This is done in the spirit of the original framework, but including a gradient term in the density. This approach is introduced to correct a deficiency in the original Equivalent Crystal Theory and other semiempirical methods in quantitatively obtaining the correct ratios of the surface energies of low index planes of cubic metals (100), (110), and (111). We develop here the basic framework, and apply it to the calculation of Fe (110) and Fe (111) surface energy formation. The results, compared with first principles calculations, show an improvement over previous semiempirical approaches.

Ferrante, John↗

NASA NTRS1992

Scientific issues and potential remote-sensing requirements for plant biochemical content

Application of developments in imaging spectrometry to the study of terrestrial ecosystems, which began in 1983, demonstrate the potential to estimate lignin and nitrogen concentrations of plant canopies by remote-sensing techniques. Estimation of these parameters from the first principles of radiative transfer and the interactions of light with plant materials is not presently possible, principally because of lack of knowledge about internal leaf scattering and specific absorption involving biochemical compounds. From the perspective of remote-sensing instrumentation, sensors are needed to support derivative imaging spectroscopy. Biochemical absorption features tend to occur in functional groupings throughout the 1100- to 2500-nm region. Derivative spectroscopy improves the information associated with the weaker, narrower absorption features of biochemical absorption that are superimposed on the strong absolute variations due to foliar biomass, pigments, and leaf water content of plant canopies. Preliminary sensor specifications call for 8-nm bandwidths at 2-nm centers in four spectral regions (about 400 bands total) and a signal-to-noise performance of at least 1000:1 for 20 percent albedo targets in the 2000-nm region.

Peterson, David L.↗

NASA NTRS2008

Deconvoluting the Friction Stir Weld Process for Optimizing Welds

In the friction stir welding process, the rotating surfaces of the pin and shoulder contact the weld metal and force a rotational flow within the weld metal. Heat, generated by the metal deformation as well as frictional slippage with the contact surface, softens the metal and makes it easier to deform. As in any thermo-mechanical processing of metal, the flow conditions are critical to the quality of the weld. For example, extrusion of metal from under the shoulder of an excessively hot weld may relax local pressure and result in wormhole defects. The trace of the weld joint in the wake of the weld may vary geometrically depending upon the flow streamlines around the tool with some geometry more vulnerable to loss of strength from joint contamination than others. The material flow path around the tool cannot be seen in real time during the weld. By using analytical "tools" based upon the principles of mathematics and physics, a weld model can be created to compute features that can be observed. By comparing the computed observations with actual data, the weld model can be validated or adjusted to get better agreement. Inputs to the model to predict weld structures and properties include: hot working properties ofthe metal, pin tool geometry, travel rate, rotation and plunge force. Since metals record their prior hot working history, the hot working conditions imparted during FSW can be quantified by interpreting the final microstructure. Variations in texture and grain size result from variations in the strain accommodated at a given strain rate and temperature. Microstructural data from a variety of FSWs has been correlated with prior marker studies to contribute to our understanding of the FSW process. Once this stage is reached, the weld modeling process can save significant development costs by reducing costly trial-and-error approaches to obtaining quality welds.

Schneider, Judy↗

DOE OSTI2020

First-principles calculation of spin and orbital contributions to magnetically ordered moments in Sr 2 IrO 4

Here, we show how an accurate first-principles treatment of the canted-anti-ferromagnetic ground state of Sr 2 IrO 4 , a prototypical 5d correlated spin-orbit coupled material, can be obtained without invoking any free parameters, such as the Hubbard $\textit{U}$ or tuning the spin-orbit coupling strength. Our theoretically predicted iridium magnetic moment of 0.250$μ_B$, canted by 12.6° off the $\textit{a}$ axis, is in accord with experimental results. By resolving the magnetic moments into their spin and orbital components, we show that our theoretically obtained variation of the magnetic scattering amplitude $\langle M_m \rangle$ as a function of the polarization angle is consistent with recent nonresonant magnetic x-ray scattering measurements. The computed value of the band gap (55 meV) is also in line with the corresponding experimental values. A comparison of the band structure to that of the cuprates suggests the presence of incommensurate charge-density wave phases in Sr 2 IrO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Vacancy-driven variations in the phonon density of states of fast neutron irradiated nuclear graphite

Research examining radiation damage and its effects in graphite began in the 1940’s during the development of moderated nuclear reactors. Interest in this topic is expanding because of emerging applications associated with fullerenes and carbon nanostructures as well as its long-standing use as a fission moderator. In this work, we report the measurements of the full phonon density of states of irradiated nuclear graphite using inelastic neutron scattering experiments at room temperature. The samples were previously exposed to different levels of neutron doses and irradiation temperatures. The phonon density of states of perfect and defected graphite supercells including different configurations of vacancies and interstitials were calculated using first-principles direct method. Despite high neutron irradiation doses, the scattering measurements show that the induced damage is localized, and the layered structure is preserved. In conclusion, a comparison of the measured phonon densities of states of irradiated samples with those calculated of defected supercells indicates that the main changes observed are mainly attributed to formation of vacancies.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Revealing the Chemical Bonding in Adatom Arrays via Machine Learning of Hyperspectral Scanning Tunneling Spectroscopy Data

The adatom arrays on surfaces offer an ideal playground to explore the mechanisms of chemical bonding via changes in the local electronic tunneling spectra. While this information is readily available in hyperspectral scanning tunneling spectroscopy data, its analysis has been considerably impeded by a lack of suitable analytical tools. Here we develop a machine learning based workflow combining supervised feature identification in the spatial domain and unsupervised clustering in the energy domain to reveal the details of structure-dependent changes of the electronic structure in adatom arrays on the Co 3 Sn 2 S 2 cleaved surface. This approach, in combination with first-principles calculations, provides insight for using artificial neural networks to detect adatoms and classifies each based on their local neighborhood comprised of other adatoms. These structurally classified adatoms are further spectrally deconvolved. The unexpected inhomogeneity of electronic structures among adatoms in similar configurations is unveiled using this method, suggesting there is not a single atomic species of adatoms, but rather multiple types of adatoms on the Co 3 Sn 2 S 2 surface. This is further supported by a slight contrast difference in the images (or slight size variation) of the topography of the adatoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1992

Free-edge effects in laminates under extension, bending and twisting. II - Sublaminate/layer modeling and analysis

The stress-function-based variational method of Yin (1991) is extended and modified into a combined layer/sublaminate approach applicable to a laminated strip composed of a large number of differently orientated, anisotropic elastic plies. Lekhnitskii's (1963) stress functions are introduced into two interior layers adjacent to a particular interface. The remaining layers are grouped into an upper sublaminate and a lower sublaminate. The stress functions are expanded in truncated power series of the thickness coordinate, and the differential equations governing the coefficient functions are derived by using the complementary virtual work principle. The layer/sublaminate approach limits the dimension of the eigenvalue problem to a fixed number irrespective of the number of layers in the sublaminate, so that reasonably accurate solutions of the interlaminar stresses can be computed with extreme ease. For symmetric, four-layer, angle-ply and cross-ply laminates, a comparison of the previous analysis results based on the pure layer model and new results based on two different layer/sublaminate models indicates reasonable over-all agreement in the interlaminar stresses and superior agreement in the total peeling and shearing force.

Yin, Wan-Lee↗

NASA NTRS1986

The effects of mantle and anelasticity on nutations, earth tides, and tidal variations in rotation rate

The paper models the effects of mantle anelasticity on luni-solar nutations, on tidal deformation, on tidal variations in rotation rate, and on the eigenfrequency of the free core nutation. The results can be used to invert observations to solve for the anelastic contributions to the shear and bulk moduli of the upper and lower mantle. Specific anelastic models are used to numerically estimate the effects of anelasticity on these geodetic observables. The nutation estimates are compared with observational results. Among the conclusions: (1) mantle anelasticity is likely to be the most important source of damping for the free core nutation; (2) present VLBI nutation results are, in principle, accurate enough to usefully bound anelasticity at diurnal periods. But the discrepancy between the VLBI observed nutations and the 1984 IAU nutation model cannot be explained by anelasticity and is not yet well enough understood to allow anelasticity to be determined from the data.

Wahr, John↗

NASA NTRS1989

Minimal complexity control law synthesis

A paradigm for control law design for modern engineering systems is proposed: Minimize control law complexity subject to the achievement of a specified accuracy in the face of a specified level of uncertainty. Correspondingly, the overall goal is to make progress towards the development of a control law design methodology which supports this paradigm. Researchers achieve this goal by developing a general theory of optimal constrained-structure dynamic output feedback compensation, where here constrained-structure means that the dynamic-structure (e.g., dynamic order, pole locations, zero locations, etc.) of the output feedback compensation is constrained in some way. By applying this theory in an innovative fashion, where here the indicated iteration occurs over the choice of the compensator dynamic-structure, the paradigm stated above can, in principle, be realized. The optimal constrained-structure dynamic output feedback problem is formulated in general terms. An elegant method for reducing optimal constrained-structure dynamic output feedback problems to optimal static output feedback problems is then developed. This reduction procedure makes use of star products, linear fractional transformations, and linear fractional decompositions, and yields as a byproduct a complete characterization of the class of optimal constrained-structure dynamic output feedback problems which can be reduced to optimal static output feedback problems. Issues such as operational/physical constraints, operating-point variations, and processor throughput/memory limitations are considered, and it is shown how anti-windup/bumpless transfer, gain-scheduling, and digital processor implementation can be facilitated by constraining the controller dynamic-structure in an appropriate fashion.

Bernstein, Dennis S.↗

NASA NTRS2000

Retinex Preprocessing for Improved Multi-Spectral Image Classification

The goal of multi-image classification is to identify and label "similar regions" within a scene. The ability to correctly classify a remotely sensed multi-image of a scene is affected by the ability of the classification process to adequately compensate for the effects of atmospheric variations and sensor anomalies. Better classification may be obtained if the multi-image is preprocessed before classification, so as to reduce the adverse effects of image formation. In this paper, we discuss the overall impact on multi-spectral image classification when the retinex image enhancement algorithm is used to preprocess multi-spectral images. The retinex is a multi-purpose image enhancement algorithm that performs dynamic range compression, reduces the dependence on lighting conditions, and generally enhances apparent spatial resolution. The retinex has been successfully applied to the enhancement of many different types of grayscale and color images. We show in this paper that retinex preprocessing improves the spatial structure of multi-spectral images and thus provides better within-class variations than would otherwise be obtained without the preprocessing. For a series of multi-spectral images obtained with diffuse and direct lighting, we show that without retinex preprocessing the class spectral signatures vary substantially with the lighting conditions. Whereas multi-dimensional clustering without preprocessing produced one-class homogeneous regions, the classification on the preprocessed images produced multi-class non-homogeneous regions. This lack of homogeneity is explained by the interaction between different agronomic treatments applied to the regions: the preprocessed images are closer to ground truth. The principle advantage that the retinex offers is that for different lighting conditions classifications derived from the retinex preprocessed images look remarkably "similar", and thus more consistent, whereas classifications derived from the original images, without preprocessing, are much less similar.

Thompson, B.↗

NASA NTRS1990

Derivatives of buckling loads and vibration frequencies with respect to stiffness and initial strain parameters

A uniform variational approach to sensitivity analysis of vibration frequencies and bifurcation loads of nonlinear structures is developed. Two methods of calculating the sensitivities of bifurcation buckling loads and vibration frequencies of nonlinear structures, with respect to stiffness and initial strain parameters, are presented. A direct method requires calculation of derivatives of the prebuckling state with respect to these parameters. An adjoint method bypasses the need for these derivatives by using instead the strain field associated with the second-order postbuckling state. An operator notation is used and the derivation is based on the principle of virtual work. The derivative computations are easily implemented in structural analysis programs. This is demonstrated by examples using a general purpose, finite element program and a shell-of-revolution program.

Haftka, Raphael T.↗

NASA NTRS1994

Feedback Control of Plasma Instabilities with Charged Particle Beams and Study of Plasma Turbulence

A new non-perturbing technique for feedback control of plasma instabilities has been developed in the Columbia Linear Machine (CLM). The feedback control scheme involves the injection of a feedback modulated ion beam as a remote suppressor. The ion beam was obtained from a compact ion beam source which was developed for this purpose. A Langmuir probe was used as the feedback sensor. The feedback controller consisted of a phase-shifter and amplifiers. This technique was demonstrated by stabilizing various plasma instabilities to the background noise level, like the trapped particle instability, the ExB instability and the ion-temperature-gradient (ITG) driven instability. An important feature of this scheme is that the injected ion beam is non-perturbing to the plasma equilibrium parameters. The robustness of this feedback stabilization scheme was also investigated. The principal result is that the scheme is fairly robust, tolerating about 100% variation about the nominal parameter values. Next, this scheme is extended to the unsolved general problem of controlling multimode plasma instabilities simultaneously with a single sensor-suppressor pair. A single sensor-suppressor pair of feedback probes is desirable to reduce the perturbation caused by the probes. Two plasma instabilities the ExB and the ITG modes, were simultaneously stabilized. A simple 'state' feedback type method was used where more state information was generated from the single sensor Langmuir probe by appropriate signal processing, in this case, by differentiation. This proof-of-principle experiment demonstrated for the first time that by designing a more sophisticated electronic feedback controller, many plasma instabilities may be simultaneously controlled. Simple theoretical models showed generally good agreement with the feedback experimental results. On a parallel research front, a better understanding of the saturated state of a plasma instability was sought partly with the help of feedback. A plasma instability is usually observed in its saturated state and appears as a single feature in the frequency spectrum with a single azimuthal and parallel wavenumbers. The physics of the non-zero spectral width was investigated in detail because the finite spectral width can cause "turbulent" transport. One aspect of the "turbulence" was investigated by obtaining the scaling of the linear growth rate of the instabilities with the fluctuation levels. The linear growth rates were measured with the established gated feedback technique. The research showed that the ExB instability evolves into a quasi-coherent state when the fluctuation level is high. The coherent aspects were studied with a bispectral analysis. Moreover, the single spectral feature was discovered to be actually composed of a few radial harmonics. The radial harmonics play a role in the nonlinear saturation of the instability via three-wave coupling.

Tham, Philip Kin-Wah↗

DOE OSTI2023

Optimizing the Synthetic Potential of O 2 : Implications of Overpotential in Homogeneous Aerobic Oxidation Catalysis

Molecular oxygen is the quintessential oxidant for organic chemical synthesis, but many challenges continue to limit its utility and breadth of applications. Extensive historical research has focused on overcoming kinetic challenges presented by the ground-state-triplet electronic structure of O 2 and the various reactivity and selectivity challenges associated with reactive oxygen species derived from O 2 reduction. This Perspective will analyze thermodynamic principles underlying catalytic aerobic oxidation reactions, borrowing concepts from the study of the oxygen reduction reaction (ORR) in fuel cells. This analysis is especially important for "oxidase"-type liquid-phase catalytic aerobic oxidation reactions, which proceed by a mechanism that couples two sequential redox half-reactions: (1) substrate oxidation, and (2) oxygen reduction, typically affording H 2 O 2 or H 2 O. The catalysts for these reactions feature redox potentials that lie between the potentials associated with the substrate oxidation and oxygen reduction reactions, and changes in the catalyst potential lead to variations in effective overpotentials for the two half reactions. Catalysts that operate at low ORR overpotential retain more thermodynamic driving force for the substrate oxidation step, enabling O 2 to be used in more challenging oxidations. While catalysts that operate at high ORR overpotential have less driving force available for substrate oxidation, they often exhibit different or improved chemoselectivity relative to the high-potential catalysts. The concepts are elaborated in a series of case studies to highlight their implications for chemical synthesis. Examples include comparisons of (a) NO x /oxoammonium and Cu/nitroxyl catalysts, (b) high-potential quinones and amine oxidase biomimetic quinones, and (c) Pd aerobic oxidation catalysts, with or without NO x cocatalysts. In addition, we show how reductive activation of O 2 provides a means to access potentials not accessible with conventional oxidase-type mechanisms. Altogether, this analysis highlights the central role of catalyst overpotential in guiding the development of aerobic oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Vibrational properties and thermal transport in quaternary chalcogenides: The case of Te-based compositions

Vibrational thermal properties of CuZn 2 InTe 4 , AgZn 2 InTe 4 , and Cu 2 CdSnTe 4 , derived from binary II-VI zinc-blendes, are reported based on first-principles calculations. While the chalcogenide atoms in these materials have the same lattice positions, the cation atom arrangements vary, resulting in different crystal symmetries and subsequent properties. The compositional differences have important effects on the vibrational thermal characteristics of the studied materials, which demonstrate that low-frequency optical phonons hybridize with acoustic phonons and lead to enhanced phonon-phonon scattering and low lattice thermal conductivities. The phonon density of states, mode Grüneisen parameters, and phonon scattering rates are also calculated, enabling deeper insight into the microscopic thermal conduction processes in these materials. Compositional variations drive differences among the three materials considered here; nonetheless, their structural similarities and generally low thermal conductivities (0.5–4 W/m K at room temperature) suggest that other similar II-VI zinc-blende derived materials will also exhibit similarly low values, as also corroborated by experimental data. Finally, this, combined with the versatility in designing a variety of motifs on the overall structure, makes quaternary chalcogenides interesting for thermal management and energy conversion applications that require low thermal conductivity.

36 MATERIALS SCIENCE↗