Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compositional simulation”

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,495 records · Page 83

Application of Probabilistic Analysis to Aircraft Impact Dynamics

Full-scale aircraft crash simulations performed with nonlinear, transient dynamic, finite element codes can incorporate structural complexities such as: geometrically accurate models; human occupant models; and advanced material models to include nonlinear stressstrain behaviors, laminated composites, and material failure. Validation of these crash simulations is difficult due to a lack of sufficient information to adequately determine the uncertainty in the experimental data and the appropriateness of modeling assumptions. This paper evaluates probabilistic approaches to quantify the uncertainty in the simulated responses. Several criteria are used to determine that a response surface method is the most appropriate probabilistic approach. The work is extended to compare optimization results with and without probabilistic constraints.

Lyle, Karen H.↗

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗

Low-velocity impact tests on fibrous composite sandwich structures

The effect of low-velocity projectile impact on the load-carrying ability of the composite sandwich structural components is investigated experimentally, the impact simulating the damage caused by runway debris and the accidental dropping of hand tools during servicing on secondary aircraft structures made with composites. The sandwich-type beam specimens were fabricated with graphite/epoxy face sheets, aluminum honeycomb core, and a steel (back) plate. A four-point beam-loading apparatus was used, and the ultimate strength, ultimate strain, and residual strength of the composites were determined. A faired curve is presented indicating the lower bound of the failure threshold for each of the laminate configurations tested in compression and tension as a function of the projectile impact energy. It is shown that strength degradation due to impact is dependent on the laminate configuration and the fiber/matrix combination. The laminates having more angle plies near the impact surface and unidirectional plies elsewhere seem to show extensive interply and intraply fiber delaminations at failure relative to the laminates with a cross-ply on the impact surface.

Sharma, A. V.↗

Can Fractional Crystallization of a Lunar Magma Ocean Produce the Lunar Crust?

New techniques enable the study of Apollo samples and lunar meteorites in unprecedented detail, and recent orbital spectral data reveal more about the lunar farside than ever before, raising new questions about the supposed simplicity of lunar geology. Nevertheless, crystallization of a global-scale magma ocean remains the best model to account for known lunar lithologies. Crystallization of a lunar magma ocean (LMO) is modeled to proceed by two end-member processes - fractional crystallization from (mostly) the bottom up, or initial equilibrium crystallization as the magma is vigorously convecting and crystals remain entrained, followed by crystal settling and a final period of fractional crystallization [1]. Physical models of magma viscosity and convection at this scale suggest that both processes are possible. We have been carrying out high-fidelity experimental simulations of LMO crystallization using two bulk compositions that can be regarded as end-members in the likely relevant range: Taylor Whole Moon (TWM) [2] and Lunar Primitive Upper Mantle (LPUM) [3]. TWM is enriched in refractory elements by 1.5 times relative to Earth, whereas LPUM is similar to the terrestrial primitive upper mantle, with adjustments made for the depletion of volatile alkalis observed on the Moon. Here we extend our earlier equilibrium-crystallization experiments [4] with runs simulating full fractional crystallization

Rapp, Jennifer F.↗

Effect of spatial density variation and O+ concentration on the growth and evolution of electromagnetic ion cyclotron waves

We simulate electromagnetic ion cyclotron (EMIC) wave growth and evolution within three regions, the plasmasphere (or plasmaspheric plume), the plasmapause, and the low-density plasmatrough outside the plasmapause. First, we use a ring current simulation with a plasmasphere model to model the particle populations that give rise to the instability for conditions observed on 9 June 2001. Then, using two different models for the cold ion composition, we do a full scale hybrid code simulation in dipole coordinates of the EMIC waves on a meridional plane at MLT = 18 and at 1900 UT within a range of L shell from L = 4.9 to 6.7. EMIC waves were observed during June 9, 2001 by Geostationary Operational Environmental Satellite (GOES) spacecraft. While an exact comparison between observed and simulated spectra is not possible here, we do find significant similarities between the two, at least at one location within the region of largest wave growth. We find that the plasmapause is not a preferred region for EMIC wave growth, though waves can grow in that region. The density gradient within the plasmapause does, however, affect the orientation of wave fronts and wave vector both within the plasmapause and in adjacent regions. There is a preference for EMIC waves to be driven in the He+ band (frequencies between the O+ and He+ gyrofrequencies) within the plasmasphere, although they can also grow in the plasmatrough. If present, H+ band waves are more likely to grow in the plasmatrough. This fact, plus L dependence of the frequency and possible time evolution toward lower frequency waves, can be explained by a simple model. Large O+ concentration limits the frequency range of or even totally quenches EMIC waves. This is more likely to occur in the plasmatrough at solar maximum. Such large O+ concentration significantly affects the H+ cutoff frequency and hence the width in frequency of the stop band above the He+ gyrofrequency. EMIC wave surfaces predicted by cold plasma theory are altered by the finite temperature of the ring current H+.

58 GEOSCIENCES↗

Analysis of Monolith Cores from an Engineering Scale Demonstration of a Prospective Cast Stone Process

The primary disposition path of Low Activity Waste (LAW) at the DOE Hanford Site is vitrification. A cementitious waste form is one of the alternatives being considered for the supplemental immobilization of the LAW that will not be treated by the primary vitrification facility. Washington River Protection Solutions (WRPS) has been directed to generate and collect data on cementitious or pozzolanic waste forms such as Cast Stone. This report documents the coring and leach testing of monolithic samples cored from an engineering-scale demonstration (ES Demo) with non-radioactive simulants. The ES Demo was performed at SRNL in October of 2013 using the Scaled Continuous Processing Facility (SCPF) to fill an 8.5 ft. diameter x 3.25 ft. high container with simulated Cast Stone grout. The Cast Stone formulation was chosen from the previous screening tests. Legacy salt solution from previous Hanford salt waste testing was adjusted to correspond to the average LAW composition generated from the Hanford Tank Waste Operation Simulator (HTWOS). The dry blend materials, ordinary portland cement (OPC), Class F fly ash, and ground granulated blast furnace slag (GGBFS or BFS), were obtained from Lafarge North America in Pasco, WA. In 2014 core samples originally obtained approximately six months after filling the ES Demo were tested along with bench scale molded samples that were collected during the original pour. A latter set of core samples were obtained in late March of 2015, eighteen months after completion of the original ES Demo. Core samples were obtained using a 2” diameter x 11” long coring bit. The ES Demo was sampled in three different regions consisting of an outer ring, a middle ring and an inner core zone. Cores from these three lateral zones were further segregated into upper, middle and lower vertical segments. Monolithic core samples were tested using the Environmental Protection Agency (EPA) Method 1315, which is designed to provide mass transfer rates (release rates) of inorganic analytes contained in monolithic material under diffusion controlled release conditions as a function of leaching time. Compressive strength measurements and drying tests were also performed on the 2015 samples. Leachability indices reported are based on analyte concentrations determined from dissolution of the dried samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation

The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Response of the North Atlantic Bloom to NAO Forcing

Results from the climatologically forced coupled ice/ocean/biogeochemical model that covers the Arctic and North Atlantic Oceans are presented and compared to the chlorophyll fields of satellite-derived ocean color measurements. Biogeochemical processes in the model are determined from the interactions among four phytoplankton functional groups (diatoms, chlorophytes, cyanobacteria and coccolithophores) and four nutrients (nitrate, ammonium, silicate and dissolved iron). The model simulates the general large-scale pattern in April, May and June, when compared to both satellite-derived and in situ observations. The subpolar North Atlantic was cool in the 1980s and warm in the latter 1990s, corresponding to the CZCS and SeaWiFS satellite observing periods, respectively. The oceanographic conditions during these periods resemble the typical subpolar upper ocean response to the NAO+ and NAO-phases, respectively. Thus, we use the atmospheric forcing composites from the two NAO phases to simulate the variability of the mid-ocean bloom during the satellite observing periods. The model results show that when the subpolar North Atlantic is cool, the NAO+ case, more nutrients are available in early spring than when the North Atlantic is warm, the NAO-case. However, the NAO+ simulation produces a later bloom than the NAO-simulation. This difference in the bloom times is also identified in SeaWiFS and CZCS satellite measurements. In the model results, we can trace the difference to the early diatom bloom due to a warmer upper ocean. The higher nutrient abundance in the NAO+ case did not provide larger total production than in the NAO- case, instead the two cases had a comparable area averaged amplitude. This leads us to conclude that in the subpolar North Atlantic, the timing of the spring phytoplankton bloom depends on surface temperature and the magnitude of the bloom is not significantly impacted by the nutrient abundance.

Mizoguchi, Ken-Ichi↗

Instrumental and observational requirements for space-based imaging of magnetospheric emissions

Simulated images of extreme ultraviolet (EUV) emissions from energetic outflowing ions have been constructed to study techniques for remotely sensing the dynamic behavior of hot plasmas in the near-Earth environment. These calculations include realistic assumptions about the energetic ion outflow from high latitudes and take into account the effects of cold plasmaspheric and ionospheric ions. The energetic ion outflow is determined from a statistical study based on five years of measurements from the Energetic Ion Composition Spectrometer on Dynamics Explorer 1. The simulated images change significantly with viewing geometry and certain spacecraft locations are clearly favorable for observing emissions from energetic ions. For example, for a near equatorial orbit, viewing locations greater than 9 Earth radii are required to observe outflowing ions above the cold plasmaspheric background. We will discuss other important considerations for magnetospheric imaging including the sensitivity requirements of the detector. In particular, we consider the performance of multi-layer optics for EUV wavelengths.

Robinson, R. M.↗

Molecular Modeling of Surfactant Interaction on Phospholipid Bilayers Mimicking Corneal Epithelium

Surfactants found in consumer products can compromise eye corneal membrane integrity upon accidental exposure. Traditional in vitro and in vivo approaches to evaluate membrane–surfactant interaction pose experimental limitations such as species variability, reproducibility, and most often do not provide the overall picture. These limitations motivate the use of in silico models to study phenomena like cellular disruption assays caused by surfactants at the molecular scale. In this work, coarse-grained molecular dynamics simulations have been employed to investigate how nonionic alcohol ethoxylate (AE) and anionic surfactant alcohol ethoxy sulfate (AES) interact with lipid bilayer liposomes that mimic corneal epithelial cell membranes. The spherical liposome is composed of 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phospho-l-serine (DOPS), and cholesterol, resembling the composition of the corneal epithelial cells’ membrane bilayer. The simulation consisted of varying degrees of representative surfactant compositions and two initial types of surfactant configurations within or outside the liposome. Our results reveal that both surfactants induce outer leaflet bulging, agreeing with membrane solubilization models. The more highly ethoxylated surfactant, AE, caused more consistent inner leaflet disruption than AES, resulting in significantly more water permeation and membrane thinning. In addition, both surfactants increase the lateral diffusion of lipids within the membrane layers, with higher ethoxylated AE showing a stronger effect than AES. This study demonstrates how surfactant structure and localization influence bilayer membrane integrity, offering mechanistic insights into the irritation potential, thus guiding the rational design of effective surfactant-based formulations.

Lipids↗

SU-E-T-380: Particle Microdosimetry Study Based On 3D-Cylindrical Silicon Radiaton Detectors

Purpose: A new design of a solid-state-microdetector based on silicon 3D microfabrication and its performance to characterise Lineal energy, Specific Energy, dose, LET and other microdosimetric variables required for modelling particle relative biological effectiveness (RBE) is presented. Methods: A microdosimeter formed by a matrix of independent sensors with well-defined micrometric cylindrical shape and with a volume similar to those of cellular dimensions is used to measure microdosimetric variables. Each sensor measures the radiation deposited energy which, divided by the mean cord length of the sensors, provides us with the Linear Energy (y) of the radiation as well as its energy distribution, and frequencymean. Starting from the these distributions in different points of a proton beam, we generate biophysical data (e.g. Linear Energy Transfer (LET), Specific Energy (z), etc…) needed for relative biological effectiveness (RBE) calculations radiation effect models used in particle radiotherapy treatment planning. In addition, a Tissue Equivalent Proportional Counter (TEPC) will be used as baseline to calibrate the “y” magnitude of the microdosimeter unit-cells. Results: The experimental measurements will soon be carried out at the Perelman Center for Advanced Medicine (University of Pennsylvania), which provides proton beam for clinical research proposals. The results of distributions measured of the microdosimetric variables from the first tests developed in the proton facility will be presented and compared with Monte Carlo simulations using the Geant4 code. Conclusion: The use of 3D microdosimeters such as the one presented here will enhance the accuracy of RBE calculations normally affected by the inherent uncertainty of monte carlo simulations due to the approximation of material composition and energy dependent physical laws involved in such calculations. The effect of such approximations will be quatified by comparison with absolute measurement of radiation quality parameters.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗

The effects of radiation on the interlaminar fracture toughness of a graphite/epoxy composite

An experimental study is made of the effect of electron irradiation (10 to the 10th rad), simulating a 30-year geosynchronous orbit exposure, on the fracture toughness of a graphite/epoxy composite, T300/934. The double cantilever beam (DBC) test is used to determine Mode I (peel) critical strain energy release rate and the edge delamination tension (EDT) test is used to determine mixed Mode I and II (peel and shear) critical strain energy release rate. It is found that the electron interaction of the epoxy matrix material enhances the fracture toughness properties of the composite and that the test temperature has a significant effect on the fracture toughness of both baseline and irradiated material.

Funk, J. G.↗

SENSITIVITY ANALYSIS OF FUEL PHYSICAL PROPERTY EFFECTS ON SPARK IGNITION ENGINE PERFORMANCE

Alternative fuels are of interest to automakers and regulators due to their potential to reduce net greenhouse gas emissions from transportation sources. Alternative fuels also have fuel properties which may enable advanced combustion modes with higher engine thermal efficiencies. There has been previous work to identify the relationship between various fuel properties and engine performance, but most of this work has been experiments or simulations where the change in properties was obtained through changing the fuel composition, making isolating the effects of individual fuel properties difficult. In this study, numerical simulations have been used to investigate the effects of individual fuel physical properties such as viscosity or heat of vaporization (HoV) on engine performance. Simulations have been performed of two different engine platforms, the first an optical, single-cylinder research engine and the second a multi-cylinder production engine. Both engines are direct-injection spark-ignition engines with pent-roof heads and are designed for automotive applications. Each engine was run at a different operating condition, one stable and one knock-limited. Different base fuels provided a variety of simulated conditions. Up to six different fuel properties were varied as part of Global Sensitivity Analyses performed for each of the engines with multiple performance targets including thermal efficiency, combustion efficiency and combustion phasing. Results show trends that are largely consistent with previous experimental findings using multiple fuels. The engine thermal efficiency was primarily sensitive to the fuel's HoV, with other fuel physical properties having smaller effects. For optical engine results, the magnitude of the effect was greater in this study than expected based on previous experimental results were many fuel physical and chemical properties were varied simultaneously. However, for the multi-cylinder production engine, the relationship between thermal efficiency and HoV was slightly smaller.

engine simulation↗

Quantitative interpretation of molecular dynamics simulations for X-ray photoelectron spectroscopy of aqueous solutions

Over the past decade, energy-dependent ambient pressure X-ray photoelectron spectroscopy (XPS) has emerged as a powerful analytical probe of the ion spatial distributions at the vapor (vacuum)-aqueous electrolyte interface. These experiments are often paired with complementary molecular dynamics (MD) simulations in an attempt to provide a complete description of the liquid interface. There is, however, no systematic protocol that permits a straightforward comparison of the two sets of results. XPS is an integrated technique that averages signals from multiple layers in a solution even at the lowest photoelectron kinetic energies routinely employed, whereas MD simulations provide a microscopic layer-by-layer description of the solution composition near the interface. Here, we use the National Institute of Standards and Technology database for the Simulation of Electron Spectra for Surface Analysis (SESSA) to quantitatively interpret atom-density profiles from MD simulations for XPS signal intensities using sodium and potassium iodide solutions as examples. We show that electron inelastic mean free paths calculated from a semi-empirical formula depend strongly on solution composition, varying by up to 30% between pure water and concentrated NaI. The XPS signal thus arises from different information depths in different solutions for a fixed photoelectron kinetic energy. XPS signal intensities are calculated using SESSA as a function of photoelectron kinetic energy (probe depth) and compared with a widely employed ad hoc method. SESSA simulations illustrate the importance of accounting for elastic-scattering events at low photoelectron kinetic energies (<300 eV) where the ad hoc method systematically underestimates the preferential enhancement of anions over cations. Finally, some technical aspects of applying SESSA to liquid interfaces are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High‐Performance and Stable Semi‐Transparent Perovskite Solar Cells through Composition Engineering

Abstract Semi‐transparent perovskite solar cells (ST‐PeSCs) have tremendous potential as solar windows owing to their higher efficiency and visible transmittance. However, studies toward this application are still nascent, particularly in unraveling the interplay between how the perovskite composition impacts the achievable device performance and stability. Here, the role of A‐ and X‐site modification in APbX 3 perovskites is studied to understand their influence on these factors. Through detailed experimental and simulation work, it is found that a perovskite composition consisting of cesium (Cs) and formamidinium (FA) at the A‐site delivers the best device performance over a range of band gaps, which are tuned by changes to the X‐site anion. Using this optimized perovskite composition, power conversion efficiencies of 15.5% and 4.1% are achieved for ST‐PeSCs with average visible transmittance values between 20.7% and 52.4%, respectively. Furthermore, the CsFA‐based ST‐PeSCs show excellent long‐term stability under continuous illumination and heating. The stability of the precursor solutions across each of the studied compositions has also been considered, showing dramatic differences in the structural properties of the perovskites and their device performance for all mixed A‐site compositions possessing the archetypal methyl ammonium species, while also confirming the superior stability of the CsFA precursor solutions.

14 SOLAR ENERGY↗

Epithermal Neutron Resonance Analysis Using a Compact D-T Generator

Nuclear material identification and isotopic composition analysis are important capabilities for safeguards and the verification of nuclear arms control and nonproliferation treaties. Epithermal neutron resonance imaging has been proposed for spent fuel assay measurements, and the development of a mobile neutron imaging system could provide such a tool for on-site applications. This technique is highly sensitive to isotopes of interest for safeguards applications (e.g. 235,238 U, 238,239,240 Pu) and yields unique transmission spectra for a target containing actinides and some mid- and high-Z elements. Paired with a position-sensitive neutron detector, such a system could provide spatial information of the target isotopic composition. Recent progress in the development of high intensity neutron sources could make mobile epithermal neutron imaging designs increasingly feasible to serve as a tool for nuclear inspectors. In this work, Monte Carlo (MC) simulations were performed to further optimize a previous design for an epithermal neutron resonance-based analysis apparatus using a D-T portable neutron generator. Additionally, experimental work is being performed to test the feasibility of this method. The MC simulations of moderation time distributions for different moderator compositions and thicknesses were measured to provide a multi-group correction for neutron time-of-flight energy reconstruction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗