Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice accretion”

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 19 records

Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach

This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A rigorous framework for an improved Messinger/Myers model of ice accretion under conditions of variable property and unsteady aircraft icing

We analyse the Messinger/Myers model by critically evaluating simplifying assumptions through a rigorous formulation of the rime ice accretion process. We explore the effects of both constant and variable ice density and thermal conductivity, along with the effects of sublimation from the ice surface. The effects of key factors such as droplet impact rate, ambient temperature relative to the freezing temperature and the temperature difference between the ambient air and the airfoil surface are examined. Under these varying conditions, the present rigorous formulation is used to assess the significance of unsteady effects, variable ice properties and sublimation. We observe that the Myers model performs remarkably well in certain icing situations and analyse the reasons for this strong performance. We also show that partially relaxing the model’s assumptions can lead to poorer performance. The Myers model can lead to overprediction of ice surface temperature and correspondingly underprediction of transition time under conditions of relatively weak sublimation and surface cooling. A modified Myers model is presented, which can be used to recover near-perfect results under widely varying icing conditions of relevance. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗

Stokes-dependent droplet collection efficiency on a NACA 0012 airfoil from droplet-informed simulations with statistical overloading

Accurate modelling of ice accretion on aircraft wings requires analysing droplet impingement on the surface to optimize the design of ice-protection systems. We perform Euler–Lagrange simulations of a droplet-laden flow impinging on a NACA 0012 airfoil. Our study includes water droplets with eight discrete sizes ranging from 1 to 160 microns. We vary the free-stream velocity of the incoming airflow in the range 60 ≤ U ≤ 240 m s −1 and the chord length of the airfoil in the range 0.5 ≤ c ≤ 2 m. Due to the dilute nature of supercooled clouds, one-way coupling is used in the simulations. The effects of droplet breakup and collision are also neglected. To reduce the computational cost, we employ statistical overloading of droplets, allowing us to simulate millions of impinging droplets in a time span on the order of milliseconds. Our results show that the droplet collection efficiency, which measures the likelihood of droplet impingement on the airfoil surface, increases with droplet size and free-stream velocity but decreases with airfoil size. We demonstrate that collection efficiency, impingement velocity and impingement angle are primarily dictated by a single non-dimensional parameter, the droplet Stokes number. We also identify a critical stagnation-streamline Stokes number below which impingements do not occur and use it to estimate the minimum droplet size for impingement. In addition, we observe droplet behaviour to become Stokes number independent at large values of the Stokes number. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

Aircraft Icing Severity Evaluation

Aircraft icing refers to the ice buildup on the surface of an aircraft flying in icing conditions. The ice accretion on the aircraft alters the original aerodynamic configuration and degrades the aerodynamic performances and may lead to unsafe flight conditions. Evaluating the flow structure, icing mechanism and consequences is of great importance to the development of an anti/deicing technique. Studies have shown computational fluid dynamics (CFD) and machine learning (ML) to be effective in predicting the ice shape and icing severity under different flight conditions. CFD solves a set of partial differential equations to obtain the air flow fields, water droplets trajectories and ice shape. ML is a branch of artificial intelligence and, based on the data, the self-improved computer algorithms can be effective in finding the nonlinear mapping relationship between the input flight conditions and the output aircraft icing severity features.

Li, Sibo (ORCID:000000021705844X)↗

High-efficiency frost and ice control via sensing-assisted nanovibrational slippery surfaces

Frost and ice accretions on surfaces pose persistent challenges across numerous industrial, residential and transportation systems. While various removal strategies exist, they often suffer from limited effectiveness or high energy consumption, such as frosting delay, ice crack generation, and Joule heating. Here, in this work, we report a novel integrated approach combining vibrational quasi-liquid surface (QLS) and capacitive sensing for efficient condensate, frost, and ice management. Compared to Joule heating, our approach does not rely on complete melting and evaporation for removal, resulting in 68% and 95% energy savings for frost and ice removal, respectively. Our QLS coating significantly reduces surface retention forces, achieving 91% and 87% less residual mass compared to hydrophilic surfaces for frost and ice removal through surface nanovibration, respectively. The integrated capacitive sensor provides real-time detection of different phase states, enabling on-demand removal in precise timeframes. This sensor-assisted approach showed 3.8 times lower energy consumption compared to conventional Joule heating for defrosting. This synergistic integration of surface engineering, nanovibration, and intelligent sensing represents a significant advancement in phase change processes, offering an energy-efficient solution for frost and ice mitigation in energy-intensive systems.

Shen, Yuchen [Univ. of Texas at Dallas, Richardson↗

Fast Evaluation of Aircraft Icing Severity Using Machine Learning Based on XGBoost

Aircraft icing represents a serious hazard in aviation which has caused a number of fatal accidents over the years. In addition, it can lead to substantial increase in drag and weight, thus reducing the aerodynamics performance of the airplane. The process of ice accretion on a solid surface is a complex interaction of aerodynamic and environmental variables. The complex relationship makes machine learning-based methods an attractive alternative to traditional numerical simulation-based approaches. In this study, we introduce a purely data-driven approach to find the complex pattern between different flight conditions and aircraft icing severity prediction. The supervised learning algorithm Extreme Gradient Boosting (XGBoost) is applied to establish the prediction framework which makes prediction based on any set of observations. The input flight conditions for the proposed prediction framework are liquid water content, droplet diameter and exposure time. The proposed approach is demonstrated in three cases: maximum ice thickness prediction, icing area prediction and icing severity level evaluation. Performance comparison studies and error analysis are also conducted to verify the effectiveness and performance of the proposed method. Results show that the proposed method has reasonable capability in evaluating aircraft icing severity.

Li, Sibo (ORCID:000000021705844X)↗

Abundant ammonia and nitrogen-rich soluble organic matter in samples from asteroid (101955) Bennu

Organic matter in meteorites reveals clues about early Solar System chemistry and the origin of molecules important to life, but terrestrial exposure complicates interpretation. Samples returned from the B-type asteroid Bennu by the Origins, Spectral Interpretation, Resource Identification, and Security–Regolith Explorer mission enabled us to study pristine carbonaceous astromaterial without uncontrolled exposure to Earth’s biosphere. Here we show that Bennu samples are volatile rich, with more carbon, nitrogen and ammonia than samples from asteroid Ryugu and most meteorites. Nitrogen-15 isotopic enrichments indicate that ammonia and other N-containing soluble molecules formed in a cold molecular cloud or the outer protoplanetary disk. We detected amino acids (including 14 of the 20 used in terrestrial biology), amines, formaldehyde, carboxylic acids, polycyclic aromatic hydrocarbons and N-heterocycles (including all five nucleobases found in DNA and RNA), along with ~10,000 N-bearing chemical species. All chiral non-protein amino acids were racemic or nearly so, implying that terrestrial life’s left-handed chirality may not be due to bias in prebiotic molecules delivered by impacts. The relative abundances of amino acids and other soluble organics suggest formation and alteration by low-temperature reactions, possibly in NH 3 -rich fluids. Bennu’s parent asteroid developed in or accreted ices from a reservoir in the outer Solar System where ammonia ice was stable.

79 ASTRONOMY AND ASTROPHYSICS↗

Comets and carbonaceous chondrites delivered noble gases to the Moon

Trapped xenon isotopes in two Apollo 17 mare basalt fragments are similar to those in primitive meteorites. Xenon would have been effectively excluded from the Moon in the aftermath of its formation in a giant impact. As such, lunar mantle xenon trapped in the mare basalts is best explained as being derived from late accretion occurring before or concurrent with the formation of the lunar crust. Here, the isotopic composition indicates that a combination of comets and carbonaceous chondrites delivered this xenon to the lunar mantle. The inferred mass of accreted cometary ice would have delivered significantly less than a part per million of water to the lunar mantle. The inferred mass of accreted carbonaceous chondrites would have supplied at least a half of a part per million of water. The data further indicate that enstatite chondrites are unlikely to have supplied the majority of late accreted mass.

Chemistry↗

A numerical simulation of CCN impacts on weather modification efficiency

Aerosols affect development of clouds and precipitation by serving as cloud condensation nuclei (CCN) and ice nuclei (IN). Considering the dramatically changing ambient aerosol concentration, it is important to examine the potential “side effect” of aerosol pollution on precipitation enhancement by weather modification. In this study, the cloud seeding was performed on a precipitation event in Beijing in the summer of 2008, which is simulated by the NSSL two-moment cloud scheme of the Weather Research and Forecasting (WRF) model. Sensitivity tests were conducted by modifying the ambient aerosol concentration and the ice crystal seeding amount to investigate the cloud seeding efficacy in different CCN concentration scenarios. There was a slight difference in the precipitation distribution between the simulations with two ambient CCN concentrations: the northern precipitation center in polluted scenario was weaker and the southern center was stronger. Compared with normal CCN scenario, the cloud liquid water mass and ice crystal mass in the severe pollution scenario is larger, and the total contents of snow and graupel were not sensitive to the CCN concentration. With the same amount of man-made ice crystals seeding, the precipitation enhancement was quite different under different CCN conditions. The higher the CCN concentration usually leads to stronger precipitation suppression. As CCN concentration increase, the deposition growth of snow, auto-conversion and accretion of ice crystals to snow were weakened, as well as the conversion of melting snow and graupel into rainwater.

54 ENVIRONMENTAL SCIENCES↗

Triple oxygen and hydrogen stable isotope composition of hydroxyl water in Orgueil (CI-type) and Tagish Lake (C2-type) carbonaceous chondrites

The primitive carbonaceous chondrites are of interest to cosmochemical science because they contain relatively large amounts of ‘water’ (H 2 O and/or OH – ) within phyllosilicate minerals. This water is evidence for the accretion of ices by their parent planetesimals, and thus represents an archive of the isotopic compositions of H 2 O in the protoplanetary environments. Here, in this study, we used thermogravimetry-enabled laser spectroscopy (TGA-IRIS) analyses of the Orgueil and Tagish Lake meteorites to make δ 2 H, δ 18 O, and Δ′ 17 O measurements of the H 2 O and OH – contained in the different hydrous minerals that comprise each meteorite. In Orgueil, we measured mass-weighted averages of OH – in the saponite and serpentine phyllosilicate matrix to be δ 2 H = 192 ‰, δ 18 O = 1.5 ‰, which are unquestionably of extraterrestrial origin with Δ′ 17 O = 1.0 ‰. For Tagish Lake, analogous values of OH– in the saponite and serpentine phyllosilicate matrix are δ 2 H = 704 ‰, δ 18 O = 11.3 ‰, and are similarly unambiguously extraterrestrial with Δ′ 17 O = 0.82 ‰. We estimate that the parent H 2 O involved in aqueous alteration of Orgueil had δ 18 O value ≥ +23 ‰. In Orgueil, we interpret the phyllosilicate petrographic relationships, and the δ 18 O values of OH – in saponite and serpentine to indicate that saponite formed first, at a lower temperature by 35 to 53 °C than serpentine. This suggests that the Orgueil parent body experienced increasing temperature during the phase of active aqueous alteration (prograde) which set the δ 18 O OH values of the serpentine and saponite. In the case of Tagish Lake, serpentine formed at a lower temperature by 32 to 60 °C than saponite, for the simplest case with constant δ 18 O H2O values. If serpentine formed first, followed by saponite formation at 32 to 60 °C °C higher temperature, this suggests that Tagish Lake sample TL1 underwent prograde aqueous alteration as the parent body heated up. We find evidence that the parent H 2 O for Orgueil, Tagish Lake sample TL1, and Murchison (based on data from a previous study) may have had Δ′ 17 O values of >0.64 ‰.

Aqueous alteration↗

Formation and Evolution Simulations of Saturn, Including Composition Gradients and Helium Immiscibility

In this paper, the formation of Saturn is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. Previous models are enhanced to include the dissolution of accreting planetesimals, composed of water ice, rock, and iron, in the gaseous envelope of the planet, leading to a nonuniform composition with depth. The immiscibility of helium in metallic hydrogen layers is also considered. The calculations start at a mass of 0.5 Earth masses and are extended to the present day. At 4.57 Gyr, the model, proceeding outward, has the following structure: (i) a central core composed of 100% heavy elements and molecules, (ii) a region with a decreasing heavy-element mass fraction, down to a value of 0.1, (iii) a layer of uniform composition with the helium mass fraction Y enhanced over the primordial value, (iv) a helium rain region with a gradient in Y, (v) an outer convective, adiabatic region with uniform composition in which Y is reduced from the primordial value, and (vi) the very outer layers where cloud condensation of the heavy elements occurs. Models of the distribution of heavy elements as a function of radius are compared with those derived to fit the observations of the Cassini mission, with rough qualitative agreement. The helium mass fraction in Saturn’s outer layers is estimated to be around 20%. Models are found that provide good agreement with Saturn’s intrinsic luminosity and radius.

79 ASTRONOMY AND ASTROPHYSICS↗

Facilitating Large‐Scale Snow Shedding from In‐Field Solar Arrays using Icephobic Surfaces with Low‐Interfacial Toughness

Abstract Large‐scale accrual of snow and ice on solar arrays in northern latitudes can cause significant power generation losses during winter. Depending on environmental conditions, snow can encompass a wide range in physical characteristics from dry snow (modulus ≈100 kPa and density ≈0.1 g cm −3 ) to bulk ice (modulus ≈8 GPa and density ≈0.9 g cm −3 ). This variation in snow morphology has made the development of a passive, broad‐spectrum, snow and ice‐shedding surface challenging. Here, the authors develop one of the first surfaces that simultaneously possesses both low‐interfacial strength ( τ˄ ice < 50 kPa) and toughness (Γ ice < 0.5 J m −2 ) with ice. These surfaces, fabricated via the addition of mobile polymer chains/oils to a thin polymeric coating, require extremely low detachment forces for ice, enabling its passive shedding at virtually any accretion length scale. Preliminary evidence that the new surfaces can shed different forms of snow and ice from field‐deployed solar arrays, over a range of subzero temperatures for several weeks, leading to significant increases in power generation is provided. The optically transparent surfaces are easily scalable and can be widely deployed by the solar industry in areas that see persistent snow. Other applications include automotive windshields, LIDAR covers for autonomous vehicles, and cold climate optical sensors.

Dhyani, Abhishek↗

Evidence that the Hot Jupiter WASP-77 A b Formed Beyond Its Parent Protoplanetary Disk’s H2O Ice Line

Idealized protoplanetary disk and giant planet formation models have been interpreted to suggest that a giant planet's atmospheric abundances can be used to infer its formation location in its parent protoplanetary disk. It has recently been reported that the hot Jupiter WASP-77 A b has subsolar atmospheric carbon and oxygen abundances with a solar C/O abundance ratio. Assuming solar carbon and oxygen abundances for its host star WASP-77 A, WASP-77 A b's atmospheric carbon and oxygen abundances possibly indicate that it accreted its envelope interior to its parent protoplanetary disk's H 2 O ice line from carbon-depleted gas with little subsequent planetesimal accretion or core erosion. We show that the photospheric abundances of carbon and oxygen in WASP-77 A are supersolar with a subsolar C/O abundance ratio, implying that WASP-77 A b's atmosphere has significantly substellar carbon and oxygen abundances with a superstellar C/O ratio. Our result possibly indicates that WASP-77 A b's envelope was accreted by the planet beyond its parent protoplanetary disk's H 2 O ice line. While numerous theoretical complications to these idealized models have now been identified, the possibility of nonsolar protoplanetary disk abundance ratios confound even the most sophisticated protoplanetary disk and giant planet formation models. We therefore argue that giant planet atmospheric abundance ratios can only be meaningfully interpreted relative to the possibly nonsolar mean compositions of their parent protoplanetary disks as recorded in the photospheric abundances of their dwarf host stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhanced Convective Microphysics Scheme and Its Impacts on Mean Climate in E3SM

Abstract To improve the representation of microphysical processes in convective clouds and their interaction with aerosol and stratiform clouds, a two‐moment convective microphysics parameterization (CMP) scheme developed by Song and Zhang (2011, https://doi.org/10.1029/2010jd014833 ) is upgraded and implemented in E3SM. The new developments include: (a) implementing a parameterization for graupel to enhance the representation of ice‐phase microphysical processes; (b) representing the impact of spatial inhomogeneity of cloud droplets in cumulus ensembles on autoconversion and accretion processes to improve the representation of warm‐rain microphysical processes; (c) implementing a comprehensive Bergeron process parameterization to better represent mixed‐phase microphysical processes; and (d) representing the interactions between ice‐phase microphysics and cloud thermodynamics. Simulations show that the cloud microphysical properties simulated by the CMP are generally in good agreement with observations. It reasonably simulates the changes in droplets effective radius related to precipitation formation in convective clouds, as identified from satellite observations. It also successfully simulates the contrast in these processes between maritime and continental clouds, demonstrating its capability to simulate the impact of aerosols on convection. Analyses of the impact of CMP on climate mean state simulation demonstrate that the CMP slightly improves the simulations of precipitation, cloud macrophysical properties, longwave cloud radiative forcing, zonal wind, and temperature. However, a degradation in shortwave cloud radiative forcing occurs.

GCM↗

Tracing the Rain Formation Pathways in Numerical Simulations of Deep Convection

Quantifying the microphysical process contributions to surface precipitation in numerical simulations can be challenging. This is due to the fact that many microphysical processes contribute to the formation and depletion of rain drops and there is almost always a spatial/temporal mismatch between where/when rain is formed and where/when it strikes the surface. In this work, we develop a tracing method that tracks the sources and sinks of raindrop mass and number as they are advected by the Weather Research and Forecasting model. Applying the method to an idealized squall line confirms that convective precipitation is dominated by warm rain processes (autoconversion and accretion) while stratiform precipitation is dominated by the melting of rimed and unrimed ice crystals. Sensitivity experiments in which the prescribed cloud drop number concentration is increased confirm the conventional wisdom that weakened autoconversion increases the fraction of raindrops originating from cold rain processes. The method also reveals that when applied to deep convection the Khairoutdinov and Kogan autoconversion scheme produces an excessive number of raindrops which are subsequently clipped in P3 microphysics to keep the rain size distribution within prescribed limits. This problem can mostly be mitigated by increasing the assumed radius for raindrops created by autoconversion.

58 GEOSCIENCES↗

Mechanism for the Efficient Homogeneous Nucleation of Ice in a Weakly Ionized, Ultracold Plasma

Abstract It is proposed that the rapid observed homogeneous nucleation of ice dust in a cold, weakly ionized plasma depends on the formation of hydroxide (OH − ) by fast electrons impacting water molecules. These OH − ions attract neutral water molecules because of the high dipole moment of the water molecules and so hydrates of the form (OH) − (H 2 O) n are formed. The hydrates continuously grow in the cold environment to become macroscopic ice grains. These ice grains are negatively charged as a result of electron impact and so continue to attract water molecules. Because hydroxide is a negative ion, unlike positive ions, it does not suffer recombination loss from collision with plasma electrons. Recombination with positive ions is minimal because positive ions are few in number (weak ionization) and slow-moving as result of being in thermal equilibrium with the cold background gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Detectability of standing accretion shock instabilities activity in supernova neutrino signals

We introduce a novel methodology for establishing the presence of standing accretion shock instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice- based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test- statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda et al. [Astrophys. J. 851, 62 (2017).] as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to ~6 kpc for IceCube and up to ~3 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗