Glass Phase Research Articles

Magnetic properties of Co1-xFexCr2S4 (x = 0–0.4) solid solutions

The magnetic properties of the Co1-xFexCr2S4 solid solutions for compositions adjacent to CoCr2S4 were investigated using various techniques, including magnetization measurements and AC susceptibility. The results revealed that substituting Fe for Co in the CoCr2S4 compound led to changes in the magnetic behavior, such as the presence of ferrimagnetic ordering and the emergence of a cluster spin glass phase. The study found that the temperature of the ferrimagnetic transition (TC) decreased with increasing Fe concentration, indicating the influence of iron on the magnetic properties of the solid solutions. Additionally, the discovery of a cluster spin glass phase for specific compositions (x = 0, 0.1, 0.2, 0.4) suggests complex magnetic behavior in these materials. In summary, the study offers valuable insights into the magnetic properties of Co1-xFexCr2S4 solid solutions and highlights the potential for tuning their magnetic behavior through compositional variations.

Investigation of phase and microstructural transformations in different alumina-silica fibers at varied heat treatment temperatures: Exploring the mechanism of Silica's influence

The study investigated the phase transformation and microstructural comparison of two types of continuous alumina fibers with varying silica contents. These fibers were subjected to heat treatment at temperatures of 1100, 1200, 1300, and 1400 °C in air. By comparing the experimental results of the two different fibers, relevant influencing mechanisms were derived. The appearance of mullite phase in the fibers occurred when the heat treatment temperature reached 1200 °C. Variations in silicon content led to inconsistent pinning effects, resulting in differences in the occurrence of irregular structures on the surface of the two types of fibers. Following heat treatment at 1300 °C, the fibers exhibited scale-like and mosaic-like structures, accompanied by phase aggregation. Subsequent heat treatment at 1400 °C resulted in the fibers displaying layered structures composed of α-Al2O3 phase and undergoing significant changes in microstructure. Notably, the two types of fibers exhibited a significant difference in the presence or absence of a glass phase, and the addition of silicon oxide to a certain extent was found to inhibit the crystal transformation of Al2O3. Additionally, the experimental results were analyzed in terms of influencing mechanisms based on theories such as activation energy barriers and total free energy.

Enhanced mechanical and thermophysical properties of Mg2Ca, Al2Ca, and Al4Ca bulk metallic glasses in comparison to crystalline alloys for bone grafting applications: A molecular dynamics investigation

We investigate the glassy-state properties of Mg2Ca, Al2Ca, and Al4Ca from the grafting application viewpoint. We employed classical molecular dynamics to examine the phase transition, structural, thermodynamic, transport, and mechanical properties in the amorphous state. All properties suggest successful simulations of the glass phase at and below the glass transition temperature, ranging between 550 and 689 K for Mg2Ca, Al2Ca, and Al4Ca. Computed results are compared and discussed with the reported findings and known mechanical and thermal properties of the various parts of the human bones and biocomposites. The comparison establishes that the mechanical, thermal, and transport properties significantly improve in the glass phase compared to its crystalline alloy form. At 300 K, studied glasses have densities in close agreement with human bone density. Structural analysis and heat capacity show the second-order phase transition, verifying the formation of the glass structure. The targeted glasses exhibit excellent thermal conductivity and thermal diffusivity compared to other commonly used biocomposites for bone grafting. Furthermore, the simulated elastic properties, viz., the Poisson ratio, G/B ratio, Cauchy's pressure, and yield strength, are in close agreement with the mechanical properties of various parts of human bone. The predicted ductility nature, contrary to the brittle character of Mg2Ca, Al2Ca, and Al4Ca crystalline alloys, proves the superiority of the glassy form for the implant's functioning. The minimum enthalpy of formation and thermodynamic stability of studied compounds benefit the synthesis process; hence, we propose that the studied glasses are persuasive materials for experimental synthesis aimed at bone grating applications.

Subnanoscale Ion Beam Modification-Assisted Smoothing of Heterostructure Surfaces.

Glass ceramic (GC) is the most promising material for objective lenses for extreme ultraviolet lithography that must meet the subnanometer precision, which is characterized by low values of high spatial frequency surface roughness (HSFR). However, the HSFR of GC is typically degraded during ion beam figuring (IBF). Herein, a developed method for constructing molecular dynamics (MD) models of GC was presented, and the formation mechanisms of surface morphologies were investigated. The results indicated that the generation of the dot-like microstructure was the result of the difference in the erosion rate caused by the difference in the intrinsic properties between ceramic phases (CPs) and glass phases (GPs). Further, the difference in the microstructure of the IBF surface under different beam angles was mainly caused by the difference in the two types of sputtering. Quantum mechanical calculations showed that the presence of interstitial atoms would result in electron rearrangement and that the electron localization can lead to a reduction in CP stability. To obtain a homogeneous surface, the effects of beam parameters on the heterogeneous surface were systematically investigated based on the proposed MD model. Then, a novel ion beam modification (IBM) method was proposed and demonstrated by TEM and GIXRD. The range of ion beam smoothing parameters that could effectively converge the HSFR of the modified surface was determined through numerous experiments. Using the optimized beam parameters, an ultrathin homogeneous modified surface within 3 nm was obtained. The HSFR of GC smoothed by ion beam modification-assisted smoothing (IBMS) dropped from 0.348 to 0.090 nm, a 74% reduction. These research results offer a deeper understanding of the morphology formation mechanisms of the GC surfaces involved in ion beam processing and may point to a new approach for achieving ultrasmooth heterostructure surfaces down to the subnanometer scale.

Glass and nanocrystalline phase formation in CuZrAg alloys: Insights from combinatorial thin film libraries studied by mapping synchrotron X-ray diffraction

We report detailed diffraction peak analysis and parameters indicative of the transition between fully amorphous microstructures and nanocrystal formation in sputter deposited thin film metallic glasses. Specifically, we fabricated ternary CuZrAg alloy films (2µm thickness) with a large compositional gradient on flexible, X-ray transparent polymer substrates in high vacuum conditions. The combinatorial libraries were subjected to point-by-point structural (XRD synchrotron radiation) and chemical (X-ray fluorescence spectroscopy) analysis, characterizing roughly 172 alloys. A strong correlation was found between peak symmetry and width, and early stages of crystalline phase formation. These nanocrystals are difficult to detect from the evolution of the peak position. In view of current literature, our data questions the recently claimed universality of the peak width as indicator of high glass forming ability. The addition of Ag significantly improves the poor glass formation of binary CuZr (36 alloys) under identical deposition conditions, as fully amorphous alloys are found for the majority of the investigated composition space. A thorough understanding of the microstructure of CuZrAg alloys on flexible substrates is highly relevant in view of favourable antibacterial properties and potential application as coatings for biomedical surfaces.

Exploring the corrosion mechanism of oxide glasses using advanced solid-state nuclear magnetic resonance spectroscopy

Despite over 70 years of study, the mechanism of glass corrosion, particularly the formation of the amorphous alteration layer known as the gel layer, remains a subject of controversy. A critical debate revolves around the origin of the abrupt change in chemical concentration near the interface between the gel layer and the glass phase. The prevailing concept attributes this shift to the dissolution of glass compositions, followed by their precipitation onto the glass surface, ultimately forming the gel layer. In this study, we utilized state-of-the-art solid-state nuclear magnetic resonance (SSNMR) techniques to track the atomic-scale network evolution of phosphate glass and borosilicate glass during corrosion. The SSNMR results demonstrate that water partially disrupts the glass network at the surface but hardly disrupts the glass network of the hydrated glass phase. This distinct interaction leads to different structures at the surface compared to the hydrated glass phase. A significant exchange of alkali metal ions in the glass occurs with H species (H+ or H3O+) in the surrounding water. Our findings demonstrate that the ion-selective diffusion, combined with the structural difference between the gel layer and hydrated glass phase, rather than the dissolution-reprecipitation of glass compositions, is responsible for the formation of a gel layer on the glass surface and the abrupt change in chemical concentration near the interface between the gel layer and the hydrated glass phase. These new findings provide valuable and novel insights into the corrosion mechanism of oxide glasses.

Multiscale investigation on the formation path of the apatite phase in bioactive glasses

Mesoporous bioactive glasses (MBGs) have been reported as promising biomaterials to stimulate and promote the growth of bones. When exposed to Simulated Body Fluid (SBF), MBGs develop an apatite phase called carbonate hydroxyapatite (H(C)A) on their surfaces that closely mimics the mineral phase of bones. In order to gain deeper insights into the key stages of the surface reactions of two different types of MBGs (58S and 70S30C) soaked in SBF, we have used high energy total X-ray scattering coupled to Pair Distribution Function (PDF) analysis, along with Raman spectroscopy and scanning electron microscopy. This type of coupled analytical approach can ultimately help to unravel the details of the reaction kinetics in the complex calcium phosphate environment around the BGs surfaces. The data obtained in our in-vitro study point to the formation of other intermediate phases like amorphous calcium phosphate (ACP) and calcite (CC), that transform after a specific time into H(C)A depending on the composition of the MBGs, their surface properties and their interaction times with SBF.

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro–Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density–temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Colouring mechanism and firing technology of white porcelain from the Sui Dynasty (581–618) from the Xiangzhou kiln, Anyang, China

The Xiangzhou kiln was one of the first significant kilns in China to fire white porcelain and as well as the earliest one to adopt the technique of Make-up clay. It plays a very important role in the history of Chinese ceramics. Here we utilize ten pieces of Sui (581–618) white porcelain from the Xiangzhou kiln as the research objects. The chemical composition, micro-morphology, firing temperature, firing atmosphere, micro-area composition and other characteristics of white porcelain glaze were tested by energy dispersive X-ray fluorescence spectroscopy (EDXRF), laser Raman spectroscopy (LRS), optical microscope (OM), thermal dilatometer (TD), focused ion beam scanning electron microscope combined with X-ray energy spectrometer (SEM-EDS), X-ray photoelectron spectroscopy (XPS) and spectrophotometer (UV–Vis–NIR). Finally, the firing technology and colouring mechanism of the Xiangzhou white porcelain are revealed. Our study demonstrates that the relatively high content of Fe2O3 and TiO2 in the white porcelain glaze of the Xiangzhou kiln results in a milky white glaze with a slight greenish tinge. The high content of CaO and K2O in the glaze and the almost entirely glass phase of glaze layer lead to the fine grain, good transparency and gloss. Most of the white porcelain of the Xiangzhou kiln is fired at high temperature under reducing atmosphere. The Xiangzhou kiln white porcelain glaze belongs to the crystalline glaze, glaze disorderly distribution of some calcium feldspar precipitation crystals and a very small amount of quartz, calcite residual crystals, glaze bubbles are larger, the distribution of relatively dispersed. The overall colour effect of the Xiangzhou kiln white porcelain results from the combined action of the chemical colour of the body and glaze, along with structural colour.

High temperature resistant composite adhesive with a remarkable bonding strength in a wide temperature range from 25 °C to 1200 °C

Ensuring a consistently high-strength connection in superalloy components of space shuttles from 25 to 1200 °C is essential yet challenging. To achieve the objective, we developed an adhesive using silicon-boron modified phenolic resin (Si-BPF) as the matrix and inorganic fillers as additives. The continuous Si–O–Si–O–B–O–B skeleton endowed the Si-BPF with excellent thermal stability. In addition, the inorganic fillers inhibited the thermal decomposition of Si-BPF. The generation of low-melting-point glass phases compensated for the defects in the adhesive. Therefore, within the “300–600 °C weak strength interval”, the synergistic effect of Si-BPF and inorganic additives increased the adhesive strength to 17.26 MPa. The formation of intermetallic compounds and the ceramization of the adhesive resulted in a gradual increase in the bonding strength of the composite adhesive with increasing temperature, reaching 34.23 MPa after 1200 °C. The adhesive enabled high-strength bonding for superalloys throughout the wide temperature range, enhancing the application of superalloy components in space shuttles.

Understanding the molecular origin of aging in the three topological phases of network glasses

ABSTRACT Network glasses composed of chalcogenide and modified oxides display three topological phases: the flexible, intermediate, and stressed-rigid. Modulated-DSC experiments, Raman scattering, and molar volume measurements along with topological constraint theory have shown that glass compositions in the flexible and stressed-rigid phases display glass transitions that age with time. The change of the enthalpy of relaxation (ΔHnr ) at the glass transition temperature (Tg ) steadily evolves with waiting time which can be traced to the presence of Open Degrees of Freedom (ODF). In the isostatically constrained intermediate phase (IP), glasses possess an enthalpy of relaxation which is substantially reduced in aging due to a lack of ODF. IP glasses are composed of self-organised Isostatically Rigid Local Structures (ISRLSs) which have no ODF. The melt fragility index (m) across the three topological phases displays a Gaussian-like global minimum of m = 15 near the IP centre composition, illustrating a super-strong melt behaviour leading to the realisation of the ideal Phillips glass condition where glass forming tendency is globally optimised as demonstrated in the present work. For melt compositions away from the IP centre and into the flexible and stressed-rigid phases, the fragility index steadily increases due to the increased presence of non-ISRLSs. The non-ISRLSs admix with the ISRLSs formed near the mid IP composition, thus promoting the melt networks to relax, age, and diffuse as one goes away from IP centre composition.

Hydrogen-annealing modulated magnetocaloric effect and critical behavior of Weyl semimetallic Co3Sn2S2

Modulating the electron interaction and magnetic domains is an essential strategy to evoke novel physical properties, both for fundamental research and for related magnetic applications. The emerging annealing process attracted attention due to its capability of tuning the magnetic performance, but it has been rarely studied. Herein, we conducted hydrogen annealing on the Weyl semimetallic Co3Sn2S2 and investigated the magnetic interaction, critical behavior, and magnetocaloric effect (MCE). Hydrogen-annealed Co3Sn2S2 (HA-Co3Sn2S2) exhibits a distinct anisotropy-dependent magnetic behavior, where a typical second-order transition with an easy magnetization process is observed under H//c. However, a frustrated magnetic state composed of ferromagnetic (FM), antiferromagnetic (AFM), and/or spin glass phases is observed when H//ab. Critical analyses illustrate that the ground state of HA-Co3Sn2S2 lies between the 3D-Ising model and the Mean-filed model, indicating the presence of both long-range and short-range magnetic interactions in the system. Additionally, we found that HA-Co3Sn2S2 exhibits an enhanced saturation magnetization (Ms) and magnetic refrigeration capacity, both showing anisotropic dependence, compared to the pristine Co3Sn2S2. Our work develops hydrogen annealing to modulate electron coupling and magnetic domains, paving the way for further studies of anisotropic magnetic behavior and applications in magnetic refrigeration.

In‐depth study of a speiss/matte sample from Castillo de Huarmey, North Coast of Peru, and its implications for the pre‐Columbian production of arsenic bronze in the Central Andes

AbstractThis study aims to characterize the phase composition and chemistry of the speiss/matte sample from the Metallurgist's Burial at Castillo de Huarmey and to use the information derived from these analyses to infer the temperatures, furnace conditions, and ores associated with the smelting processes, which created the speiss/matte sample. For this purpose, a number of geochemical analyses were performed on the spies/matte fragment: analysis of the general chemical composition (handheld X‐ray fluorescence spectrometry [hhXRF], X‐ray photoelectron spectroscopy [XPS]), analysis of the chemical composition in the micro area (field emission scanning electron microscope with an energy dispersive spectroscopy detector [FE‐SEM‐EDS], field emission electron probe microanalysis [FE‐EPMA]), analysis of the mineral composition (X‐ray diffraction [XRD]), and analysis of the phase composition (Raman spectroscopy). Chemical and mineralogical analyses of the speiss/matte specimen determined that the specimen is composed of distinct arsenide, arsenate, sulfide, and glass phases. During the smelting process, the charge material consisted mainly of Cu, Fe, and As sulfides. Arsenopyrite is the most likely candidate as the mineral source of arsenic. In addition, temperatures of at least 1200°C were achieved during the smelting process, with smelting occurring over a relatively short timeframe given that effective density separation of speiss and matte phases was not achieved.

The Modulation of Compositional Heterogeneity for Controlling Shear Banding in Co-P Metallic Nanoglasses.

Shear banding is much dependent on the glass-glass interfaces (GGIs) in metallic nanoglasses (NGs). Nevertheless, the current understanding of the glass phase of GGIs is not well established for controlling the shear banding in NGs. In this study, Co-P NGs are investigated by molecular dynamics simulations to reveal the phenomenon of elemental segregation in the GGI regions where the content of Co is dominant. Specifically, Co segregation results in the formation of GGIs, whose atomic structures are comparatively less dense than those present in the interiors of glassy grains. It is suggested that the Co segregation significantly reduces the shear resistance of GGIs. Thus, such compositional heterogeneity influences the mechanical properties of Co-P NGs. Particularly, shear banding is much altered through enhancing the Co segregation in the GGI regions, which leads to improvements in the ductility of Co-P NGs. This study advances knowledge of the formation of the GGI phase in NGs, which could enable GGI engineering in enhancing the mechanical properties of NGs.

The changing technology trends in processing titaniferous magnetite

ABSTRACT Titaniferous magnetite deposits are a significant source of vanadium for use in high-strength steel alloying and other chemical applications as well as producing titanium slag for the chloride production of titanium pigment. Ore beneficiation is practised to generate an iron-based mineral concentrate containing much of the vanadium with low combined silica and alumina levels of less than 4%. If higher than 4% increased sodium salt is used in roasting lower vanadium recovery is achieved. Typical processing problems encountered are fine liberation size, potential over-grinding of the magnetite, and variable response to magnetic separation as a result of the weathering profile. Magnetic separation is the standard route to produce iron vanadium concentrates. Salt roasting to extract vanadium followed by water leaching has been the standard processing route for many years and this locks up the remaining titanium in a glass phase. The cost of suitable sodium salts, availability of capital and increasing energy costs and energy availability are critical factors in determining the viability of aspiring vanadium producers. In this paper, factors influencing the recovery of vanadium are reviewed and emerging processing trends are explored. The growing demand for vanadium electrolytes in vanadium redox flow batteries (VRFB) now adds another value-added product, which can assist vanadium project economics.

Effect of MgO addition on the structural evolution, crystallization and properties of self-glazed ceramics sintered from granite sludge

Recycling of granite sludge has both economic and environmental benefits. In this paper, high glossiness self-glazed ceramics were prepared from granite sludge with the addition of TiO2, Ca(OH)2 and MgO by combining two-step sintering and low-temperature crystallization. The effect of MgO addition on the phase evolution, crystallization and properties of the glazed ceramics was studied to reveal its role in the diopside precipitation. With increasing MgO content, the reduced glass viscosity accelerates the dissolution of quartz and anorthite into the liquid glass phase, resulting in a decrease in flexural strength above 1 % MgO. Meanwhile, the addition of MgO increases the non-bridging oxygen to lower the glass transition temperature. The increase in the Al/Si ratio of the glass phase and the mixed alkaline earth effect contribute to an increase in the crystallization temperature. After thermal treatment at 950 °C, the precipitated diopside increases the crystallinity of the glazed glass-ceramic from 6.4 % to 19.3 % at 3 % MgO, resulting in a surface glossiness of 99.7 GU and a flexural strength of 64.1 MPa. High glossiness self-glazed ceramics with an aesthetic appearance could be used in decorative building tiles for scalable utilization of granite sludge.

Compositional Characterization of Glassy Volcanic Material From VSWIR and MIR Spectra Using Partial Least Squares Regression Models

AbstractThe glass phase in volcanic rocks presents a challenge to obtaining compositional data from visible and short‐wave‐infrared (VSWIR) and mid‐infrared (MIR) spectral data of remote surfaces due to its amorphous structure and variable composition. Nonetheless, glass is a common phase in volcanic materials because it forms via the rapid quench of magma and can constitute up to the entirety of a volcanic deposit. Use of partial least squares regression (PLS) to predict glass contents creates models that are insensitive to viewing geometry and sample conditions such as grain size and spectrally inactive compositional variables, enhancing the ability to detect glasses with remote sensing. PLS models are used here to predict crystallinity and oxide composition of samples from VSWIR and MIR spectral data using training spectra from natural volcanic rocks and geologically relevant synthetic samples. Three spectral resolutions of VSWIR and MIR spectra (1, 10, and 100 nm/band, and 1.9, 19, and 190 cm−1/band, respectively) were tested to assess the effects of collection configuration on different spectrometers. PLS models trained on 1 nm and 1.9 cm−1 data sets have the lowest uncertainties of glass modal abundance for VSWIR and MIR, respectively. MIR models predicting sample wt. % SiO2 and FeO, and VSWIR models of wt. % FeO provide accurate estimates (e.g., RMSE‐P of 3.4 wt. % FeO) at all spectral resolutions. Results are based on training data sets skewed to mafic compositions, which affects model accuracies.

Effect of crystalline phase formed by compound flame retardant on the flame retardancy and ceramization of polyethylene composites

AbstractCeramic polyolefin composites have the capability to transform into hard ceramics when exposed to fire conditions. During the ceramization process, the formation of new crystalline phase plays a crucial role in enhancing flame‐retardant and ceramifiable properties. Consequently, ceramic polyolefin composites show great potential for the applications in fire‐resistant wires and cables. In this article, the incorporation of the compound flame retardant consisting of ammonium polyphosphate/melamine cyanurate/zinc borate (APP/MCA/ZB) was found to enhance the flame retardancy and ceramization of polyethylene/wollastonite fiber/phosphate glass frits (PE/WF/PGF) composites. The results indicated that ceramifiable flame‐retarding PE composites with compound flame retardant exhibited superior flame retardancy compared to pure PE and PE composites with a single flame retardant. Specifically, the limiting oxygen index (LOI) was significantly increased to 26.8%, and the vertical combustion test rating in UL‐94 (test for flammability of plastic materials for parts in devices and appliances) reached V‐0. During the heating process, ZB thermally decomposed to produce 2ZnO ⋅ 3B2O3, which reacted with CaSiO3 to form a silicate glass intermediate phase (CaO ⋅ SiO2 ⋅ 2ZnO ⋅ 3B2O3). APP thermally decomposed to produce (HPO3)n, which reacted with 2ZnO ⋅ 3B2O3 to form a phosphate glass intermediate phase (nP2O5 ⋅ 2ZnO ⋅ 3B2O3). These two glass phases experienced a eutectic reaction with WF, ultimately producing the formation of a new crystalline phase of calcium zinc phosphate (CZP, Ca19Zn2(PO4)14). This newly formed CZP phase made sintered ceramics more compact and had higher flexural strength. The flexural strength of ceramic residues after sintering was 11.68 MPa, meeting the requirements for practical applications.

Modeling of the temperature‐dependent tensile strength of polymer considering the effect of the two‐phase transition of glass and rubber phases

AbstractThis study explores the contribution of glass phase and rubber phase to polymer failure, and establishes a theoretical model of temperature‐dependent tensile strength (TDTS) considering the effect of two‐phase transition based on the Force‐heat Equivalence Energy Density Principle (FHEEDP). The model achieves good verification by comparison with the TDTS experimental data of various polymer materials, including epoxy‐based shape memory polymer (SMP), PEEK‐based SMP and CdS/PMMA. The results demonstrate that the model can perform an accurate prediction for the tensile strength of polymer at different temperatures by easily obtainable material parameters. In conclusion, the proposed theoretical model deepens the understanding on the two‐phase transition and their effect on the TDTS, as well as supplies a practical predicted method on the TDTS of polymer, which is crucial for polymer applied in extreme condition.

Evaluating the residual stresses and fracture behavior of glass-ceramics by numerical modeling

Residual stresses often arise after glass crystallization due to the residual glass and crystal phases' thermal and elastic properties mismatch. Using the Positional Finite Element Method, we developed a computational tool for nonlinear 2D geometric analysis to simulate residual stresses from crystal growth in a glass matrix. We implemented a damage model in which interface elements are introduced into the mesh. These elements' elastic modules were penalized after reaching a tensile stress limit to simulate crack propagation. The crystal mesh was generated separately and included in the glass mesh via kinematic compatibility without introducing degrees of freedom nor requiring mesh conformity. Results were consistent with analytical calculations using Selsing's model. The fracture mode in crystal-glass composite microstructures with different signs of thermal expansion mismatch is closely in accordance with the microstructure pictures in the literature. The adapted damage model showed excellent potential in evaluating and predicting crack patterns due to residual stresses in glass-ceramics.