Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Composite models”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The Impact of Neutron Transfer Reactions on the Heating and Cooling of Accreted Neutron Star Crusts

Abstract Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and long-term cooling in transiently accreting systems. It was recently suggested that previously ignored neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

PDG: A Composite Method Based on the Resolution of the Identity

The Gaussian-3 (G3) composite approach for thermochemical properties is revisited in light of the enhanced computational efficiency and reduced memory costs by applying the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs) to the computationally demanding component methods in the G3 model: the energy and gradient computations via the second-order Møller–Plesset perturbation theory (MP2) and the energy computations using the coupled-cluster singles–doubles method augmented with noniterative triples corrections [CCSD(T)]. Efficient implementation of the RI-based methods is achieved by employing a hybrid distributed/shared memory model based on MPI and OpenMP. The new variant of the G3 composite approach based on the RI approximation is termed the RI-G3 scheme, or alternatively the PDG method. The accuracy of the new RI-G3/PDG scheme is compared to the “standard” G3 composite approach that employs the memory-expensive four-center ERIs in the MP2 and CCSD(T) calculations. Taking the computation of the heats of formation of the closed-shell molecules in the G3/99 test set as a test case, it is demonstrated that the RI approximation introduces negligible changes to the mean absolute errors relative to the standard G3 model (less than 0.1 kcal/mol), while the standard deviations remain unaltered. Furthermore, the efficiency and memory requirements for the RI-MP2 and RI-CCSD(T) methods are compared to the standard MP2 and CCSD(T) approaches, respectively. The hybrid MPI/OpenMP-based RI-MP2 energy plus gradient computation is found to attain a 7.5× speedup over the standard MP2 calculations. For the most demanding CCSD(T) calculations, the application of the RI approximation is found to nearly halve the memory demand, confer about a 4–5× speedup for the CCSD iterations, and reduce the computational time for the compute-intensive triples correction step by several hours.

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

Heuristic Computational Model for Predicting Lignin Solubility in Tailored Organic Solvents

Lignin is a random heteropolymer that has been extensively studied as a renewable source of aromatic precursors for high-value chemicals, biofuels, and bioplastics. A key challenge in lignin valorization is the structural and compositional heterogeneity of lignin feedstocks. Solvent-based approaches are commonly used to fractionate lignin to reduce this heterogeneity, but solvent selection can be challenging due to variability in lignin composition. In this work, we developed computational methods to predict good and poor organic solvents as a function of lignin composition. We analyzed 28 different linear pentamer structures, 18 from known libraries and 10 hypothetical polymers, and calculated their activity coefficients in 50 different organic solvents by using the conductor-like screening model for realistic solvents. We used these data to train a regression model that enabled the extensive investigation of the impact of solvent and monolignol compositions on predicted lignin solubility. The exhaustive exploration of solubility trends using model predictions revealed sets of solvents, identified using Kamlet–Taft parameters, that are predicted to promote lignin dissolution regardless of lignin composition. We further identified solvents expected to selectively isolate lignin fractions enriched in certain subunits. Furthermore, these results establish heuristic guidelines for solvent selection that can be used to tailor fractionation processes for lignin feedstocks of distinct composition or to design new processes that isolate fractions with higher proportions of selected subunits.

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

A data science approach for analysis and reconstruction of spinodal-like composition fields in irradiated FeCrAl alloys

A statistical method for the analysis of continuously distributed data representative of composition fluctuations in irradiated FeCrAl alloys acquired using Energy Dispersive X-ray Spectroscopy (EDS) method is presented. Using probability distribution functions, direct and cross-covariances between the elemental compositions, the effects of alloy composition and irradiation dose were investigated on the spatial distribution and length scale of composition fluctuations at the nanoscale. We have observed that, for neutron-irradiated FeCrAl alloys, the distribution of Fe and Cr followed a left-skewed and right-skewed distribution, respectively for all (average) alloy compositions and irradiation doses. The analysis also revealed enhanced spatial gradients in the elemental compositions at higher irradiation dose. Direct and cross-covariance estimates of the experimental data were also utilized for reconstruction of composition data through fitting it to a parametric form of the covariance functions. Linear Model of Coregionalization was used to determine the parameters of the covariance functions. Subsequently, a spectral method was utilized for simulating a realization of the alloy compositions. Close correspondence was observed between the experimental and the reconstructed data which was analyzed using probability distribution functions and covariance functions. Composition space of the experimental and reconstructed data and dislocation velocities as a function of applied stress and line directions over the entire composition maps were also examined.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase diagram of magnetic shape memory alloy Ni 50 Mn $50–x$ In $x$ , 0 < $x$ , 25 from first principles, via spin cluster expansion and phonon vibrational entropies

The metamagnetic shape memory Heusler alloy Ni 50 Mn $50–x$ In $x$ exhibits a rich phase diagram featuring competing magnetic states, coupled magnetic–structural phase transitions, and strong compositional sensitivity. Existing first-principles approaches struggletocapturetheintertwinedchemical, magnetic, andvibrationaleffectsinthesealloys, necessitating a more integrated modeling framework. We develop a spin cluster expansion (spin-CE) framework augmented by a quasi-harmonic phonon model to capture both configurational (chemical and magnetic) and vibrational contributions to the free energy of Ni 50 Mn $50–x$ In $x$ over the full compositional range 0 ≤ x ≤25. The spin-CE includes both chemical clusters and composition-dependent Ising spin interactions, with parameters fit to a first-principles density functional theory (DFT) dataset. Using this approach, we predict the complete magnetostructural phase diagram and transformation temperatures of Ni 50 Mn $50–x$ In $x$ across the composition space. We find that vibrational entropy alone is insufficient to reproduce the martensitic transformation in the magnetic shape memory alloy regime, highlighting the essential role of magnetism. Incorporating both magnetic and vibrational contributions allows us to reproduce all experimentally known phases, including the disappearance of the stable martensite phase at a critical In concentration and the Curie temperature of the austenite phase. The method also captures the transition with increasing In in martensite from antiferromagnetic to ferromagnetic order and predicts re-entrant ferromagnetism, though the latter occurs at higher In content than reported experimentally. We discuss possible sources of this discrepancy and highlight the broader applicability of the method to other magnetostructurally complex materials, where it may offer mechanistic insight and predictive design capabilities.

Cluster expansion↗
DOE OSTI2023

In situ probing the origin of interfacial instability of Na metal anode

The unstable and fragile solid electrolyte interphase (SEI) has restricted the application of Na metal anodes. Despite numerous research efforts being put into understanding its chemical composition and physical properties, direct observation of its formation remains a challenge due to the lack of temporal and spatial resolution. Here, through combined in situ probing techniques, we exhibit two pivotal stages associated with SEI instability during the often -neglected formation process. It is revealed that Na metal that is not uniformly passivated at the initial (passivating) stage will trigger unrestricted electrolyte decomposition and homogeneous components distribution during the subsequent (growing) stage. SEI with homogenously distributed components is found to have higher solubility than that with a layered structure evolved from a compact passivation layer. Furthermore, through demonstrating an SEI dissolution model that is closely related to its formation process and compositional distribution, this work sheds light on an uncharted territory of Na metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Query relaxation for portable brick-based applications

Semantic metadata standards pave the way for interoperability by providing building operators and application developers with common schemes to describe building resources. Applications can query building metadata models to retrieve the set of entities and relationships they need to operate, instead of hard-coding references to specific points and objects from the underlying data sources. Currently, querying such models requires the developer to be very specific when formulating queries in order to obtain meaningful answers (or any answer). The developer is inevitably expected to be familiar with the systems and components of the buildings being queried, as well as the schema used to represent them. The variety of buildings - both in the composition of their subsystems and in how they happen to be modeled - means that the developer will need to use multiple queries in order to retrieve necessary results. This is complex, time-consuming and error-prone. To address this limitation, we investigate query relaxation as a technique to facilitate discovery of meaningful building resources in a collection of ontology-based buildings data. We evaluate our query relaxation approach over a set of Brick models and demonstrate its use in the context of real-world building applications.

Brick↗
DOE OSTI2020

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Real-time cure monitoring of fiber-reinforced polymer composites using infrared thermography and recursive Bayesian filtering

Despite the desirable characteristics of fiber-reinforced polymer (FRP) composites, their utilization in high volume production industries is limited by a lack of efficient manufacturing techniques. Monitoring the curing process of these composites can help with improving the quality and efficiency of the manufacturing process. This article discusses a method that uses Kalman-filter-based fusion of the information obtained from infrared thermography (surface temperature measurements) and a heat conduction model to estimate degrees of cure and internal temperatures of curing FRP composite parts in real time. The effectiveness of the methodology is demonstrated by successfully monitoring the curing of an FRP composite part in a laboratory experiment. Finally, the proposed methodology is a crucial step towards identifying anomalies in the curing process that negatively impact the quality of FRP composite parts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Evaluation of RMI Ejecta Sourcing from Reconstructed Interfaces in 2D FLAG Flyer-Plate Simulations

This study investigates the performance of RMI ejecta sourcing from a reconstructed interface in two-dimensional FLAG ejecta flyer-plate simulations. Enabling the RMI ejecta model to operate on reconstructed interfaces has been a longstanding challenge in FLAG, limiting its applicability to experiments in which interface reconstruction is not required. Addressing this limitation is therefore an important step toward expanding the model’s predictive capability. The first phase of this effort focused on identifying a flyer-plate/target configuration that produces large ejecta clouds from two distinct shocks using the standard RMI ejecta model. A series of simulations examined the effects of flyer composition, layer thickness, and mesh resolution on shock structure and ejecta formation. VISAR traces were used to evaluate shock timing and identify configurations that produced a clear double-shock response. A copper middle-layer flyer with a doubled tin layer yielded the most favorable shock characteristics and was selected as the baseline configuration. After confirming this configuration performed as expected when running a reference FLAG version, it became the basis for three simulations: (1) a standard two-dimensional ejecta flyer-plate calculation using a reference FLAG version, (2) the same calculation with the reconstructed interface option enabled, and (3) a calculation employing ALE strategies to force the ejecta model to operate on the reconstructed interface. While the baseline simulation produced the expected ejecta response, the reconstructed interface approach did not result in measurable ejecta, indicating that this capability remains under development. These results highlight both the progress made and the remaining challenges in extending the RMI ejecta model to reconstructed interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Direct observation on supersonic microprojectile penetration of carbon fiber composites with ultrafast synchrotron X-ray phase contrast imaging

High-speed penetration into carbon fiber composites is of fundamental importance to materials science and impact engineering, but research along this line suffers considerably from the lack of direct experimental observations. Here we investigate such penetration dynamics of a unidirectional carbon fiber reinforced epoxy (UCFRE) composite, with a combination of in situ, ultrafast, synchrotron phase contrast imaging and finite element (FE) analysis. The experiments yield the first direct observation on projectile trajectories and fiber-scale deformation and damage in the UCFRE composites during supersonic microprojectile penetration, for different fiber orientations (0 degrees - 90 degrees from the impact direction) and projectile velocities (600 - 850 m.s -1 ), at unprecedented temporal (~ 100 ps) and spatial (5 μm) resolutions. The maximum penetration depth decreases with increasing fiber orientation angles, as a result of anisotropic damage evolution in the composite sample. Strain localizations are prone to develop along a direction perpendicular to the fiber orientation, while the damage or cavity region, along the fiber direction. FE modeling with a three-dimensional Hashin criterion yields consistent projectile trajectory and cavity morphology with the experimental results. Finally, with increasing fiber orientation angles, damage analyses show a transition in the damage mode from fiber compression to matrix compression damage, in line with the increasing maximum penetration depth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantifying the impacts of marine aerosols over the southeast Atlantic Ocean using a chemical transport model: implications for aerosol–cloud interactions

The southeast Atlantic region, characterized by persistent stratocumulus clouds, has one of the highest uncertainties in aerosol radiative forcing and significant variability across climate models. In this study, we analyze the seasonally varying role of marine aerosol sources and identify key uncertainties in aerosol composition at cloud-relevant altitudes over the southeast Atlantic using the GEOS-Chem chemical transport model. We evaluate simulated aerosol optical depth (AOD) and speciated aerosol concentrations against those collected from ground observations and aircraft campaigns such as LASIC, ORACLES, and CLARIFY, conducted during 2017. The model consistently underestimates AOD relative to AERONET, particularly at remote locations like Ascension Island. However, when compared with aerosol mass concentrations from aircraft campaigns during the biomass burning period, it performs adequately at cloud-relevant altitudes, with a normalized mean bias (NMB) between -3.5 % (CLARIFY) and -7.5 % (ORACLES). At these altitudes, in the model, organic aerosols (63 %) dominate during the biomass burning period, while sulfate (41 %) prevails during austral summer, when dimethylsulfide (DMS) emissions peak in the model. Our findings indicate that marine sulfate can account for up to 69 % of total sulfate during the high-DMS period. Sensitivity analyses indicate that refining DMS emissions and oxidation chemistry may increase sulfate aerosol produced from marine sources, highlighting that there remains large uncertainty as to the role of DMS emissions in the marine boundary layer. Additionally, we find marine primary organic aerosol emissions may substantially increase total organic aerosol concentrations, particularly during austral summer. This study underscores the imperative need to refine marine emissions and their chemical transformations, as aerosols from marine sources are a major component of total aerosols at cloud-relevant altitudes and may impact uncertainties in aerosol radiative forcing over the southeast Atlantic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Self-organized patterns at the plasma–liquid anode interface in a helium glow discharge: temporal development and mechanisms

Self-organization at the plasma–liquid anode interface is a commonly observed phenomenon for atmospheric pressure glow plasmas, resulting in patterns with distinctive shapes such as circular ring, star-shaped, and gear-like structures, depending primarily on the discharge current and solution conductivity. Recent studies have shown that the electrode gap distance, solute used for liquid anode solution, and gas composition can also significantly impact pattern formation. Nonetheless, an overarching model or explanation of the key underlying mechanisms consistent with all experimentally observed trends is not yet reported. We propose a key underlying mechanism enabling pattern formation motivated by a detailed parametric study of pattern formation complemented by the temporal development of patterns and consistent with all observed trends. Pattern formation was observed to be on a time scale of 100 μs, similar to the time scales of gas heating and evaporation. It was found that a minimum water evaporation rate of (3.5 ± 0.5) × 10 –6 kg s –1 and reduced electric field in the positive column of 16.6 ± 0.4 Td is required for pattern formation in the investigated cases irrespective of solution conductivity and gas composition for NaCl solutions. Nonetheless, the presence of cations for which the corresponding metal atom has a low ionization energy was identified as a necessary condition for pattern formation. Furthermore, the reported results suggest that the presence of a small amount of metal atoms in the gas phase with low ionization energy enhances the overall ionization rate in the near anode region which triggers pattern formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Model Package Report: Central Plateau Vadose Zone Models

This model package report describes concisely the modeling objectives, conceptualization, implementation, uncertainty and sensitivity, configuration control, limitations, data needs, and recommendations for improvement of the Central Plateau Vadose Zone Models. This collection of models is developed to meet the vadose zone simulation needs of the Hanford Site Composite Analysis (CA), the Hanford Site Cumulative Impact Evaluation (CIE), and are expected to find other applications in Hanford Site remedial cleanup decision-making processes. This model package report describes and documents the development of the models themselves to fulfill technical approaches defined for the CA and the CIE. This report does not document any specific calculation using these models: instead, applications of individual vadose models to perform specific calculations will be documented in environmental calculation files, including inputs and results, as appropriate. Configuration management for the Central Plateau Vadose Zone Models is managed through assignment of unique sequential numerical version numbers and archival of the models in the Environmental Modeling Management Archive.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Dynamics and activation of membrane-bound B cell receptor assembly

B-cell receptor (BCR) complexes are expressed on the surface of a B-cell and are critical in antigen recognition and modulating the adaptive immune response. Even though the relevance of antibodies has been known for almost a hundred years, the antigen-dependent activation mechanism of B-cells has remained elusive. Several models have been proposed for BCR activation, including cross-linking, conformation-induced oligomerization, and dissociation activation models. Recently, the first cryo-EM structures of the human B-cell antigen receptor of the IgM and IgG isotypes have been published that validates the asymmetric organization of the BCR complex. Here, we carry out extensive molecular dynamics simulations to probe the conformational changes upon antigen binding and the influence of the membrane lipids. We identify two critical dynamical events that could be associated with antigen-dependent activation of BCR. First, antigen binding causes increased flexibility in regions distal to the antigen binding site. Second, antigen binding alters the rearrangement of IgM transmembrane helices, including the relative interaction of Igα/Igβ that mediates intracellular signaling. Furthermore, these transmembrane rearrangements lead to changes in localized lipid composition. Our work indirectly supports the conformational-change induced models of BCR activation and contributes to the understanding of the antigen-dependent activation mechanism of BCRs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

CFD modeling of near-wall combustion and unburned methane prediction in natural gas spark ignition engines

Natural gas-powered engines play a critical role in gas drilling, compression, and transmission sectors, but methane (CH 4 ) from engine combustion slip can be significant over their lifespan, contributing to atmospheric pollution and signaling reduced engine efficiency. Here, to address this challenge, computational fluid dynamics (CFD) simulations offer valuable insights into the in-cylinder combustion process, enabling the optimization of combustion strategies and engine designs to minimize unburned CH 4 slip. This study aims to evaluate and improve combustion models for simulating the combustion process and predicting unburned CH 4 concentrations in natural gas spark-ignition (SI) engines, including engines that are part of combined reformer-engine systems. Specifically, the performance of two flamelet-based combustion models—the Extended Coherent Flame Model (ECFM) and the G-equation model—was assessed using experimental engine data collected under varying excess-air ratio (λ) conditions and fuel compositions, including natural gas and syngas blends. In addition, to enhance the predictive capabilities of the G-equation model, a flame-wall interaction (FWI) sub-model was integrated into its framework. The effects of its model parameters, such as quenching and influence distance, on combustion behavior and unburned methane predictions were analyzed in detail. The ECFM tended to predict delayed combustion phasing under diluted mixture conditions, resulting in overprediction of unburned CH 4 concentrations. In contrast, the G-equation model provided reasonable predictions of combustion pressure, while representing higher the CH 4 reduction rate across the operating condition compared to experimental data. Incorporating the FWI sub-model—with the quenching distance calculated based on a pressure-dependent relation (P -0.48 ) and a fixed influence distance of 1.5 mm—further improved the G-equation model’s accuracy in predicting CH 4 reduction rates without compromising its ability to simulate the combustion process.

Combustion model↗
DOE OSTI2023

Aligning theoretical and empirical representations of soil carbon-to-nitrogen stoichiometry with process-based terrestrial biogeochemistry models

Soil carbon-nitrogen (C:N) stoichiometry acts as a control over decomposition and soil organic matter formation and loss, making it a key soil property for understanding ecosystem dynamics and projected ecosystems responses to global environmental change. However, the controls of soil C:N and how they respond to increasing pressures from global change agents are not fully understood. The “foundational” controls on soil C:N, namely plant and microbial C:N, have been used to predict soil C:N, but fail to accurately simulate all ecosystems and may be insufficient for predictions under global environmental change. Here we present an “emerging” representation of controls of soil C:N that includes plant-microbe-mineral feedbacks that have been shown to regulate soil C:N. We argue that including representation of these emerging drivers in process-based terrestrial biogeochemistry models, which include biological N fixation, mycorrhizae, priming, root exudation of organic acids, and mineralogy (including soil texture, mineral composition, and aggregation), will improve mechanistic representation of soil C:N and associated processes. Such improvements will produce models that will better simulate a variety of ecological states and predict soil C:N when global changes modify plant-microbe-mineral interactions. Here, we align our empirical understanding of controls of soil C:N with those controls represented in models, identifying contexts where emerging drivers might be particularly important to represent (e.g., priming and root exudation in nutrient-limited conditions) and areas of future work. Additionally, we show that implementing emerging drivers of soil C:N results in different simulated outcomes at steady state and in response to elevated atmospheric CO2. Our review and preliminary simulations support the need to incorporate emerging drivers of soil C:N into process-based terrestrial biogeochemistry models, allowing for both theoretical exploration of mechanisms and potentially more accurate predictions of land biogeochemical responses to global change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Probing interspecies metabolic interactions within a synthetic binary microbiome using genome-scale modeling

Metabolic interactions within a microbial community play a key role in determining the structure, function, and composition of the community. However, due to the complexity and intractability of natural microbiomes, limited knowledge is available on interspecies interactions within a community. In this work, using a binary synthetic microbiome, a methanotroph-photoautotroph (M-P) coculture, as the model system, we examined different genome-scale metabolic modeling (GEM) approaches to gain a better understanding of the metabolic interactions within the coculture, how they contribute to the enhanced growth observed in the coculture, and how they evolve over time. Using batch growth data of the model M-P coculture, we compared three GEM approaches for microbial communities. Two of the methods are existing approaches: SteadyCom, a steady state GEM, and dynamic flux balance analysis (DFBA) Lab, a dynamic GEM. We also proposed an improved dynamic GEM approach, DynamiCom, for the M-P coculture. SteadyCom can predict the metabolic interactions within the coculture but not their dynamic evolutions; DFBA Lab can predict the dynamics of the coculture but cannot identify interspecies interactions. DynamiCom was able to identify the cross-fed metabolite within the coculture, as well as predict the evolution of the interspecies interactions over time. A new dynamic GEM approach, DynamiCom, was developed for a model M-P coculture. Constrained by the predictions from a validated kinetic model, DynamiCom consistently predicted the top metabolites being exchanged in the M-P coculture, as well as the establishment of the mutualistic N-exchange between the methanotroph and cyanobacteria. The interspecies interactions and their dynamic evolution predicted by DynamiCom are supported by ample evidence in the literature on methanotroph, cyanobacteria, and other cyanobacteria-heterotroph cocultures.

59 BASIC BIOLOGICAL SCIENCES↗