Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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,135 records · Page 63

Significance of radiative corrections on measurements of the EMC effect

Deep inelastic scattering (DIS) from nuclear targets probes the parton distribution functions (PDFs) in nuclei. Comparisons of the PDFs from heavy nuclei and the deuteron show deviations that demonstrate a non-trivial nuclear dependence to these distributions, referred to as the EMC effect. A global analysis of the worlds data on the EMC effect reveals tensions between different extractions. Precise measurements at Jefferson Lab, studying the dependence on both the quark momentum fraction, x, and nuclear mass, show systematic discrepancies among experiments, making the extraction of the A dependence of the EMC effect sensitive to the selection of datasets. Further, by comparing various methods and assumptions used to calculate radiative corrections, we have identified differences that, while not large, significantly impact the EMC ratios and show that using a consistent radiative correction procedure resolves this discrepancy, leading to a more coherent global picture, and allowing for a more robust extraction of the EMC effect for infinite nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Large Solar Flares and Sheared Magnetic Field Configuration

This Comment gives additional information about the nature of flaring locations on the Sun described in the article "Sun unleashes Halloween storm", by R. E. Lopez, et al. What causes the large explosions from solar active regions that unleash huge magnetic storms and adverse space weather? It is now beyond doubt that the magnetic field in solar active regions harbors free energy that is released during these events. Direct measurements of the longitudinal and transverse components of active region magnetic fields with the vector magnetograph at NASA Marshall Space Flight Center (MSFC), taken on a regular basis for the last 30 years, have found key signatures of the locations of powerful flares. A vector magnetograph detects and measures the magnetic shear, which is the deviation of the observed transverse magnetic field direction from the potential field. The sheared locations possess abundant free magnetic energy for solar flares. In addition to active region NOAA 10486, the one that produced the largest flares last October, the NASA/MSFC vector magnetograph has observed several other such complex super active regions, including NOAA 6555 and 6659.

Choudhary, Debi Prasad↗

Probabilistic Path Planning of Montgolfier Balloons in Strong, Uncertain Wind Fields

Lighter-than-air vehicles such as hot-air balloons have been proposed for exploring Saturn s moon Titan, as well as other bodies with significant atmospheres. For these vehicles to navigate effectively, it is critical to incorporate the effects of surrounding wind fields, especially as these winds will likely be strong relative to the control authority of the vehicle. Predictive models of these wind fields are available, and previous research has considered problems of planning paths subject to these predicted forces. However, such previous work has considered the wind fields as known a priori, whereas in practical applications, the actual wind vector field is not known exactly and may deviate significantly from the wind velocities estimated by the model. A probabilistic 3D path-planning algorithm was developed for balloons to use uncertain wind models to generate time-efficient paths. The nominal goal of the algorithm is to determine what altitude and what horizontal actuation, if any is available on the vehicle, to use to reach a particular goal location in the least expected time, utilizing advantageous winds. The solution also enables one to quickly evaluate the expected time-to-goal from any other location and to avoid regions of large uncertainty. This method is designed for balloons in wind fields but may be generalized for any buoyant vehicle operating in a vector field. To prepare the planning problem, the uncertainty in the wind field is modeled. Then, the problem of reaching a particular goal location is formulated as a Markov decision process (MDP) using a discretized space approach. Solving the MDP provides a policy of what actuation option (how much buoyancy change and, if applicable, horizontal actuation) should be selected at any given location to minimize the expected time-to-goal. The results provide expected time-to-goal values from any given location on the globe in addition to the action policy. This stochastic approach can also provide insights not accessible by deterministic methods; for example, one can evaluate variability and risk associated with different scenarios, rather than only viewing the expected outcome.

Wolf, Michael↗

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗

Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way

Abstract We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)–like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches ≈100 cm 2 g −1 at relative velocities of ≈10 km s −1 , motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos’ gravothermal evolution, demonstrating that deviations from the collisionless R max – V max relation can be used to select deeply core-collapsed halos, where V max is a halo’s maximum circular velocity, and R max is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈10 8 M ⊙ is deeply core collapsed in our SIDM model. Core-collapsed systems form ≈10% of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above ≈10 9 M ⊙ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Plume particle collection and sizing from static firing of solid rocket motors

A unique dart system has been designed and built at the NASA Marshall Space Flight Center to collect aluminum oxide plume particles from the plumes of large scale solid rocket motors, such as the space shuttle RSRM. The capability of this system to collect clean samples from both the vertically fired MNASA (18.3% scaled version of the RSRM) motors and the horizontally fired RSRM motor has been demonstrated. The particle mass averaged diameters, d43, measured from the samples for the different motors, ranged from 8 to 11 mu m and were independent of the dart collection surface and the motor burn time. The measured results agreed well with those calculated using the industry standard Hermsen's correlation within the standard deviation of the correlation . For each of the samples analyzed from both MNASA and RSRM motors, the distribution of the cumulative mass fraction of the plume oxide particles as a function of the particle diameter was best described by a monomodal log-normal distribution with a standard deviation of 0.13 - 0.15. This distribution agreed well with the theoretical prediction by Salita using the OD3P code for the RSRM motor at the nozzle exit plane.

Sambamurthi, Jay K.↗

Extreme Weather Events and Their Relationship to Low Frequency Teleconnection Patterns

A new method for identifying the structure and other characteristics of extreme weather events is introduced and applied to both model simulations and observations. The approach is based on a linear regression model that links daily extreme precipitation amounts for a particular point on the globe to precipitation and related quantities at all other points. We present here some initial results of our analysis of extreme precipitation events over the United States, including how they are influenced by ENSO and various large-scale teleconnection patterns such as the PNA. The results are based on simulations made with the NASA/NCAR AGCM (Lin and Rood 1996). The quality of the simulated climate for the NASA/NCAR AGCM forced with observed SSTs is described in Chang et al. (2001). The runs analyzed here consist of three 20-year runs forced with idealized cold, neutral and warm ENSO SST anomalies (superimposed on the mean seasonal cycle of SST). The idealized warm or cold SST anomalies are fixed throughout each 20- year simulation and consist of the first EOF (+/- 3 standard deviations) of monthly SST data. Comparisons are made with the results obtained from a similar analysis that uses daily NOAA precipitation observations (Higgins et al. 1996) over the United States and NCEP/NCAR reanalysis data for the period 1949-1998.

Chang, Yehui↗

Real-Time Defect Correction in Large-Scale Polymer Additive Manufacturing via Thermal Imaging and Laser Profilometer

Defects can result in a failed part and are costly in terms of time and material. This cost is even greater in the context of large-scale additive manufacturing where the objects can be very large. As a result, a great deal of research has focused on defect identification and mitigation. To address defects during object construction, researchers at Oak Ridge National Laboratory’s Manufacturing Demonstration Facility investigated an in-situ control system comprised of two sensors: a thermal camera and laser profilometer. This control system adjusted material flow and build speed to mitigate three types of defects: low layer times, underfill, and overfill. Several test objects were constructed. The control system was found to adjust build parameters to handle low layer times of approximately 15 seconds and height deviations from -100% underfill (the absence of a layer) to 50% overfill. Within two layers, height deviations could be returned to within 10% of the expected layer height. Further, preliminary results suggest the system can compensate for uneven build surfaces.

Borish, Michael↗

A direct current rectification scheme for microwave space power conversion using traveling wave electron acceleration

The formation of the Vision-21 conference held three years ago allowed the present author to reflect and speculate on the problem of converting electromagnetic energy to a direct current by essentially reversing the process used in traveling wave tubes that converts energy in the form of a direct current to electromagnetic energy. The idea was to use the electric field of the electromagnetic wave to produce electrons through the field emission process and accelerate these electrons by the same field to produce an electric current across a large potential difference. The acceleration process was that of cyclotron auto-resonance. Since that time, this rather speculative ideas has been developed into a method that shows great promise and for which a patent is pending and a prototype design will be demonstrated in a potential laser power beaming application. From the point of view of the author, a forum such as Vision-21 is becoming an essential component in the rather conservative climate in which our initiatives for space exploration are presently formed. Exchanges such as Vision-21 not only allows us to deviate from the 'by-the-book' approach and rediscover the ability and power in imagination, but provides for the discussion of ideas hitherto considered 'crazy' so that they may be given the change to transcend from the level of eccentricity to applicability.

Manning, Robert M.↗

The influence of thermal inertia on Mars' seasonal pressure variation and the effect of the weather component

Using a Leighton-Murray type diurnal and seasonal Mars thermal model, we found that it is possible to reproduce the seasonal variation in daily-averaged pressures (approximately 680-890 Pa) measured by Viking Lander 1 (VL1), during years without global dust storms, with a standard deviation of less than 5 Pa. In this simple model, surface CO2, frost condensation, and sublimation rates at each latitude are determined by the net effects of radiation, latent heat, and heat conduction in subsurface soil layers. An inherent assumption of our model is that the seasonal pressure variation is due entirely to the exchange of mass between the atmosphere and polar caps. However, the results of recent Mars GCM modeling have made it clear that there is a significant dynamical contribution to the seasonal pressure variation. This 'weather' component is primarily due to large-scale changes in atmospheric circulation, and its magnitude depends somewhat on the dust content of the atmosphere. The overall form of the theoretical weather component at the location of VL1, as calculated by the AMES GCM, remains the same over the typical range of Mars dust opacities.

Wood, S. E.↗

Evaluation and Windspeed Dependence of MODIS Aerosol Retrievals Over Open Ocean

The Maritime Aerosol Network (MAN) data set provides high quality ground-truth to validate the MODIS aerosol product over open ocean. Prior validation of the ocean aerosol product has been limited to coastal and island sites. Comparing MODIS Collection 5 ocean aerosol retrieval products with collocated MAN measurements from ships shows that MODIS is meeting the pre-launch uncertainty estimates for aerosol optical depth (AOD) with 64% and 67% of retrievals at 550 nm, and 74% and 78% of retrievals at 870 nm, falling within expected uncertainty for Terra and Aqua, respectively. Angstrom Exponent comparisons show a high correlation between MODIS retrievals and shipboard measurements (R= 0.85 Terra, 0.83 Aqua), although the MODIS aerosol algorithm tends to underestimate particle size for large particles and overestimate size for small particles, as seen in earlier Collections. Prior analysis noted an offset between Terra and Aqua ocean AOD, without concluding which sensor was more accurate. The simple linear regression reported here, is consistent with other anecdotal evidence that Aqua agreement with AERONET is marginally better. However we cannot claim based on the current study that the better Aqua comparison is statistically significant. Systematic increase of error as a function of wind speed is noted in both Terra and Aqua retrievals. This wind speed dependency enters the retrieval when winds deviate from the 6 m/s value assumed in the rough ocean surface and white cap parameterizations. Wind speed dependency in the results can be mitigated by using auxiliary NCEP wind speed information in the retrieval process.

Kleidman, Richard G.↗

Incipient failure detection (IFD) of SSME ball bearings

Because of the immense noise background during the operation of a large engine such as the SSME, the relatively low level unique ball bearing signatures were often buried by the overall machine signal. As a result, the most commonly used bearing failure detection technique, pattern recognition using power spectral density (PSD) constructed from the extracted bearing signals, is rendered useless. Data enhancement techniques were carried out by using a HP5451C Fourier Analyzer. The signal was preprocessed by a Digital Audio Crop. DAC-1024I noise cancelling filter in order to estimate the desired signal corrupted by the backgound noise. Reference levels of good bearings were established. Any deviation of bearing signals from these reference levels indicate the incipient bearing failures.

Source record↗

Search for Higgs bosons produced in association with a high-energy photon via vector-boson fusion and decaying to a pair of b -quarks in the ATLAS detector

A search for Standard Model Higgs bosons produced in association with a high-energy photon and decaying to $b\bar{b}$ is performed using 133 fb −1 of $\sqrt{s}$ = 13 TeV pp collision data collected with the ATLAS detector at the Large Hadron Collider at CERN. The photon requirement reduces the multijet background, and the H → $b\bar{b}$ decay is the dominant decay mode. Event selection requirements target events produced by vector-boson fusion, the dominant production mode in this channel. Several improvements enhance the search sensitivity compared to previous measurements. These improvements include better background modelling and characterization, the use of a neural-network classifier, and an updated signal extraction strategy adopting a direct binned-likelihood fit to the classifier output. With these improvements, the Higgs boson signal strength is measured to be 0.2 ± 0.7 relative to the Standard Model prediction. This corresponds to an observed significance of 0.3 standard deviations, compared to an expectation of 1.5 standard deviations assuming the Standard Model.

ATLAS experiment↗

Predictive Coordinated and Cooperative Voltage Control for Systems With High Penetration of PV

In this paper, we propose a predictive coordinated and cooperative voltage control method in a power distribution system with high penetration of photovoltaic (PV) units. First, an integrated coordinated voltage control of voltage regulators (VRs) tap positions and cooperative distributed control of the reactive power output from PV inverters are used to maintain system voltages within an appropriate bandwidth. Next, solar power forecasting is applied to predict voltage changes, which are used to set the VR tap positions and capacitor switch status to prevent large voltage fluctuations. The fine tuning of voltage adjustment is then achieved by cooperative control of PV inverters to maintain a uniform voltage profile across the system. The proposed method is tested on a modified IEEE 123-node test feeder with high penetration of PVs using real measurement data and compared with the base case. Simulation results demonstrate the effectiveness of the integrated voltage control, as well as the enhancement from the predictive control through solar power forecasting-enabled voltage change estimates. Comparison to previous work in the literature shows significant improvement in terms of voltage deviation and reduction in excessive tap changes.

14 SOLAR ENERGY↗

Evaluation of Observed and Modelled Aerosol Lifetimes Using Radioactive Tracers of Opportunity and an Ensemble of 19 Global Models

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (Cs-137) and xenon-133 (Xe-133) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. Cs-137 size distribution measurements taken close to the power plant suggested that accumulation mode (AM) sulfate aerosols were the main carriers of cesium. Hence, Cs-137 can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas Xe-133 behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of Cs-137that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and Xe-133 emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled Cs-137and Xe-133 concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime e, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95% confidence interval 13.1-15.7 days). The equivalent modelled e lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.42.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (Xe-133) concentrations in the Arctic as well but to a smaller extent than for the aerosol (Cs-137) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

Radioactive isotopes↗

Measurement of the $t\bar{t} t\bar{t}$ production cross section in pp collisions at $ \sqrt{s} $ = 13 TeV with the ATLAS detector

A measurement of four-top-quark production using proton-proton collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS detector at the Large Hadron Collider corresponding to an integrated luminosity of 139 fb -1 is presented. Events are selected if they contain a single lepton (electron or muon) or an opposite-sign lepton pair, in association with multiple jets. The events are categorised according to the number of jets and how likely these are to contain b-hadrons. A multivariate technique is then used to discriminate between signal and background events. The measured four-top-quark production cross section is found to be 26$^{+17}_{-15}$ fb, with a corresponding observed (expected) significance of 1.9 (1.0) standard deviations over the background-only hypothesis. The result is combined with the previous measurement performed by the ATLAS Collaboration in the multilepton final state. The combined four-top-quark production cross section is measured to be 24$^{+7}_{-6}$ fb, with a corresponding observed (expected) signal significance of 4.7 (2.6) standard deviations over the background-only predictions. It is consistent within 2.0 standard deviations with the Standard Model expectation of 12.0 ± 2.4 fb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Role of Bottlebrush Additives on the Structure of Block Copolymers in the Bulk and Thin Films

Blending block copolymers (BCP) with additives is a useful approach for controlling BCP morphology and properties. In athermal systems, blends of BCPs with polymer additives having very high molecular (M n ) mass generally result in macrophase separation. Bottlebrush polymers, which consist of a linear backbone and grafted side chains, present an interesting alternative where the overall M n of the system can be very large but the low M n side chains may drive miscibility with the BCP. Here, in this study, a bottlebrush with a polynorbornene backbone and polystyrene (PS) side chains is blended with PS-b-poly(methyl methacrylate) (PS-b-PMMA) of varying M n , and the resulting morphologies are examined in both the bulk and thin films. Two different M n of PS-b-PMMA were used in the bulk study, and the analysis of small-angle X-ray scattering data shows that the blends were miscible and lamellar at all concentrations. This deviates from reference series of both low and high M n linear polymer additives, which either showed morphological transitions from lamellae to cylinders (low M n ) or were immiscible at all mass fractions studied (high M n ). The relative molecular mass of the side chain (N SC ) and the corresponding component in the BCP (N A ) dictate the distribution of the bottlebrush throughout the BCP, analogous to BCP/linear blends or grafted nanoparticles in a homopolymer matrix. The studies on thin films show a thickness dependence for bottlebrush mass fractions at or above 0.17, a behavior which may be driven by conformational changes of the bottlebrush upon confinement.

36 MATERIALS SCIENCE↗