Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large-eddy 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 235 records · Page 13

Comparison of wind-farm control strategies under realistic offshore wind conditions: wake quantities of interest

Wind-farm control strategies aim to increase the efficiency, and therefore lower the levelized cost of energy, of wind farms. This is done by using turbine settings such as the yaw angle, blade pitch angles, or generator torque to manipulate the wake that negatively affects downstream turbines in the farm. Two inherently different wind-farm control methods have been identified in the literature: wake steering (WS) and active wake mixing (AWM). As one of two companion papers focused on understanding practical aspects of these two wind-farm control strategies using large-eddy simulation (LES), we below analyze the wake quantities of interest for a single wind turbine performing WS and AWM, while the companion article (Frederik et al., 2025) focuses on turbine quantities of interest including power and structural loads for the same computational setup and also includes two-turbine arrays with full and partial wake overlap. The simulations, which are based in the LES solver AMR-Wind, are tailored to have inflow conditions representative of measurements from a site off the East Coast of the US, including with strong veer and low turbulence. The turbine, which is modeled in OpenFAST and coupled to the LES, is the IEA 15 MW, an open-source offshore design. After presenting an overview of the wake recovery for the different wake-control cases, the analysis probes the fluid-dynamic causes for the different performance of the arrays reported in the companion article by examining control volumes around the wakes and the budget of the mean-flow kinetic energy (MKE) within these volumes. In the high veer environment considered, the MKE recovery is dominated by mean convection, and this is shown to especially benefit the WS strategy when a neighboring turbine is directly downstream: there is ≈65 % more available power for a downstream turbine than in the baseline case, and this power is gained primarily through mean convection on the left-tip and top-tip faces of the control volume. However, the case with imperfect knowledge of the exact wind direction favors the pulse-type AWM strategy, largely because of ≈9 % increased turbulent entrainment from aloft versus the baseline that could be related to an apparent resistance to skewing in the pulsed wake. The general reduced effectiveness of helix-type and other individual-pitch-based AWM strategies for inflow with high veer and low turbulence as reported in the companion paper is due, in part, to low magnitudes of phase-averaged turbulent entrainment. Two main findings of this study are thus that veer has a significant impact on the effectiveness of different wake-control strategies and that pulse-type AWM may be a useful strategy when the objective is power maximization in realistic, offshore flow environments with imperfect knowledge of the exact wake overlap position on the downstream turbine.

17 WIND ENERGY↗

Wake-Resolving Acoustic Tomography: Advances through Numerical Covariance Methods

Acoustic tomography offers path-integrated measurements of atmospheric velocity and temperature fluctuations with high spatial resolution. Classical implementations of time-dependent stochastic inversion rely on homogeneous, isotropic covariance models that are poorly suited to the anisotropic structure of wind turbine wakes. By directly estimating heterogeneous covariances from large-eddy simulations (LESs) into the time-dependent stochastic inversion operator, we relax implicit assumptions in the analytical models used historically. Retrievals using these LES-informed models improve agreement with true fields in variance, turbulent kinetic energy, and spectral content compared to analytical and precursor-based covariance models. The results indicate that LES-informed covariance models can enhance the accuracy of acoustic tomography retrievals in complex, anisotropic flows such as wind turbine wakes in some cases and highlight instances where analytical models still offer competitive performance, despite their simplifying assumptions.

17 WIND ENERGY↗

Power output of turbines mounted on tension-leg platforms subjected to fully developed ocean gravity waves

A concern in the deployment of large wind turbines on ocean floating platforms is the effect of floating-platform motions on their electrical power generation. Further, it is not clear how floating motions influence waking, which might affect the combined power generation of collections of turbines. We examine the average power output of a single and a collection of NREL 5 MW wind turbines mounted on a tension-leg platform (TLP) under the action of fully developed ocean wave motions, coupling floating motions with the large-eddy simulation (LES) of atmospheric and rotor dynamics. The ocean dynamics enter as fully developed waves derived from the Pierson–Moskowitz spectrum. To assess the influence of ocean motions, we performed simulations over the full range of wind speeds in the operational range of the turbine, reporting comparisons of average power output when the platforms are allowed to move to when they are held rigidly in place. In all simulations, we find that the effects of the TLP floating-platform's induced motions have a minor effect on single and multiple turbine power production and wake deficits. Even when using coherent and large amplitude harmonic-floating-induced perturbations, any significant wake modifications from floating motions are confined to the near-wake region, where downstream turbines are unlikely to be located. The relatively small amplitude of TLP motions relative to pre-existing turbulent fluctuations are the primary reason for low wake and power modifications downstream.

Restrepo, Juan [ORNL] (ORCID:0000000326092882)↗

Introduction to and comparison of deep learning and optimization approaches to analytical wake modeling of a tilted wind turbine

This paper introduces innovative optimization and deep learning techniques to enhance the prediction of complex wake dynamics in the downstream wind velocity of tilted wind turbines. Traditional methods for calibrating the Bastankhah wake model often lead to increased errors in wind velocity distribution due to overfitting of the local wake characteristics. To address this issue, we propose an additional global optimization step to reduce errors in wind velocity predictions with respect to various wake parameters. Despite this improvement, the Bastankhah model's axisymmetric Gaussian wake shape limits its accuracy for complex wake structures. Therefore, we also propose a deep learning approach, which demonstrates promising results by accurately modeling complex wake shapes across a broader range of tilt angles with minimal computational cost. The deep learning approach achieves near-identical predictions to high-fidelity large-eddy simulations, representing a promising advancement in wake modeling.

17 WIND ENERGY↗

Increased power gains from wake steering control using preview wind direction information

Abstract. Yaw controllers typically rely on measurements taken at the wind turbine, resulting in a slow reaction to wind direction changes and subsequent power losses due to misalignments. Delayed yaw action is especially problematic in wake steering operation because it can result in power losses when the yaw misalignment angle deviates from the intended one due to a changing wind direction. This study explores the use of preview wind direction information for wake steering control in a two-turbine setup with a wind speed in the partial load range. For these conditions and a simple yaw controller, results from an engineering model identify an optimum preview time of 90 s. These results are validated by forcing wind direction changes in a large-eddy simulation model. For a set of six simulations with large wind direction changes, the average power gain from wake steering increases from only 0.44 % to 1.32 %. For a second set of six simulations with smaller wind direction changes, the average power gain from wake steering increases from 1.24 % to 1.85 %. Low-frequency fluctuations are shown to have a larger impact on the performance of wake steering and the effectiveness of preview control, in particular, than high-frequency fluctuations. From these results, it is concluded that the benefit of preview wind direction control for wake steering is substantial, making it a topic worth pursuing in future work.

17 WIND ENERGY↗

3D Continuous Forcing Dataset from 3D Constrained Variational Analysis at SGP

The continuous 3D large-scale forcing (VARANAL3D) data set derived from 3D constrained variational analysis (3DCVA) extends the conventional constrained variational analysis method by incorporating multiple sub-columns within the analysis domain. This advancement introduces spatial variability into the large-scale forcing fields, thereby enriching the data set’s applicability. The VARANAL3D data set spans from 2004 to 2018 and covers a region of 5˚×4.5˚ domain around the ARM SGP site. The analysis domain is divided into 10×9 sub-columns with 0.5˚ resolution. The 3D large-scale forcing data provides necessary variables to drive and evaluate single-column models (SCM), cloud-resolving models (CRM) ,and large-eddy simulations (LES), as well as information for testing model sensitivity to spatial variability of the large-scale forcing data, facilitating more rigorous testing and refinement of physical processes in SCM/CRM/LES.

54 ENVIRONMENTAL SCIENCES↗

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗

Evaluation of probability density function descriptions for three-component Rayleigh–Taylor mixing

Results from simulations of a three-component Rayleigh–Taylor (RT) mixing problem are presented. These simulations are conducted in heavy–light–heavy and heavy–intermediate–light configurations, and each of these configurations are further considered in high- and low-Reynolds-number regimes. This results in RT-unstable flow with one or both interfaces initially unstable, permitting the influence of problem configuration on the statistical description of three-component RT-driven mixing to be considered. Mass fraction covariances are observed to undergo a sign change through the mixing layer in all four configurations considered. This appears to be unique to the multi-component case and represents another way in which multi-component RT mixing differs from the two-component case. Qualitative and quantitative comparisons of joint and marginal probability density function (PDF) descriptions of species concentration are made. Three-, five-, and six-parameter model PDFs are compared against simulation data to assess how accurately they describe the mixing, and it is found that three-component mixing requires at least a five-parameter model PDF to accurately describe the mixing. Notably, the marginal distributions of three-component mixing do not appear to conform to a beta distribution, representing a departure from the classical two-component RT case. In conclusion, statistical neutrality also appears to influence the optimal choice of model PDF, which is found to be a function of problem configuration.

Large-eddy simulation↗

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗

Bridging New Observational Capabilities and Process-Level Simulation: Insights into Aerosol Roles in the Earth System

The spatial distribution of ambient aerosol particles significantly impacts aerosol–radiation–cloud interactions, which contribute to the largest uncertainty in global anthropogenic radiative forcing estimations. However, the atmospheric boundary layer and lower free troposphere have not been adequately sampled in terms of spatiotemporal resolution, hindering a comprehensive characterization of various atmospheric processes and impeding our understanding of the Earth system. To address this research data gap, we have leveraged the development of uncrewed aerial systems (UAS) and advanced measurement techniques to obtain mesoscale spatial data on aerosol microphysical and optical properties around the U.S. Southern Great Plains (SGP) atmospheric observatory. Our study also benefits from state-of-the-art laboratory facilities that include three-dimensional molecular imaging techniques enabled by secondary ion mass spectrometry and nanogram-level chemical composition analysis via micronebulization aerosol mass spectrometry. Through our study, we have developed a framework for observation–modeling integration, enabling an examination of how various assumptions about the organic–inorganic components mixing state, inferred from chemical analysis, affect clouds and radiation in observation-constrained model simulations. By integrating observational constraints (derived from offline chemical analysis of the aerosol surface using collected samples) with in situ UAS observations, we have identified a prominent role of organic-enriched nanometer layers located at the surface of aerosol particles in determining profiles of aerosol optical and hygroscopic properties over the SGP observatory. Furthermore, we have improved the agreement between predicted clouds and ground-based cloud lidar measurements. This UAS–model–laboratory integration exemplifies how these new advanced capabilities can significantly enhance our understanding of aerosol–radiation–cloud interactions.

54 ENVIRONMENTAL SCIENCES↗

Assessing turbulence and mixing parameterizations in the gray-zone of multiscale simulations over mountainous terrain during the METEX21 field experiment

Multiscale numerical weather prediction models transition from mesoscale, where turbulence is fully parameterized, to microscale, where the majority of highly energetic scales of turbulence are resolved. The turbulence gray-zone is situated between these two regimes and multiscale models must downscale through these resolutions. Here, we compare three multiscale simulations which vary by the parameterization used for turbulence and mixing within the gray-zone. The three parameterizations analyzed are the Mellor-Yamada Nakanishi and Niino (MYNN) Level 2.5 planetary boundary layer scheme, the TKE-1.5 large eddy simulation (LES) closure scheme, and a recently developed three-dimensional planetary boundary layer scheme based on the Mellor-Yamada model. The simulation domain includes complex (i.e., mountainous) terrain in Nevada that was instrumented with meteorological towers, profiling and scanning lidars, a tethered balloon, and a surface flux tower. Simulations are compared to each other and to observations, with assessment of model skill at predicting wind speed, wind direction and TKE, and qualitative evaluations of transport and dispersion of smoke from controlled releases. This analysis demonstrates that microscale predictions of transport and dispersion can be significantly influenced by the choice of turbulence and mixing parameterization in the terra incognita , particularly over regions of complex terrain and with strong local forcing. This influence may not be apparent in the analysis of model skill, and motivates future field campaigns involving controlled tracer releases and corresponding modeling studies of the turbulence gray-zone.

54 ENVIRONMENTAL SCIENCES↗

Surface‐Atmosphere Decoupling Prolongs Cloud Lifetime Under Warm Advection Due To Reduced Entrainment Drying

Abstract An initially well‐mixed stratocumulus deck can remain overcast for several tens of hours under warm‐advection conditions, although moisture supply is cut off from the ocean due to surface‐atmosphere decoupling (stabilization of the surface‐atmosphere interface). In this study, a set of idealized large‐eddy simulations were performed to investigate the physical mechanism of how warm‐air advection impacts the evolution of a pre‐existing stratocumulus deck. To mimic warm‐air advection, we decrease the sea surface temperature linearly over time in a doubly periodic domain. Given the same initial conditions, the stratocumulus deck is more persistent when experiencing warm‐air advection than cold‐air advection. This persistence is caused by reduced cloud‐top entrainment drying due to decoupling, a process more influential than the decoupling‐induced cutoff of moisture supply. This mechanism is more notable when the free troposphere becomes more humid. The relevance of the mechanism to previous observations of less low‐level cloudiness under warm‐advection conditions is discussed.

54 ENVIRONMENTAL SCIENCES↗

Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence

ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().

17 WIND ENERGY↗

Designing a Convection‐Cloud Chamber for Collision‐Coalescence Using Large‐Eddy Simulation With Bin Microphysics

Abstract Collisional growth of cloud droplets is an essential yet uncertain process for drizzle and precipitation formation. To improve the quantitative understanding of this key component of cloud‐aerosol‐turbulence interactions, observational studies of collision‐coalescence in a controlled laboratory environment are needed. In an existing convection‐cloud chamber (the Pi Chamber), collisional growth is limited by low liquid water content and short droplet residence times. In this work, we use numerical simulations to explore various configurations of a convection‐cloud chamber that may intensify collision‐coalescence. We employ a large‐eddy simulation (LES) model with a size‐resolved (bin) cloud microphysics scheme to explore how cloud properties and the intensity of collision‐coalescence are affected by the chamber size and aspect ratio, surface roughness, side‐wall wetness, side‐wall temperature arrangement, and aerosol injection rate. Simulations without condensation and evaporation within the domain are first performed to explore the turbulence dynamics and wall fluxes. The LES wall fluxes are used to modify the Scalar Flux‐budget Model, which is then applied to demonstrate the need for non‐uniform side‐wall temperature (two side walls as warm as the bottom and the two others as cold as the top) to maintain high supersaturation in a tall chamber. The results of LES with full cloud microphysics reveal that collision‐coalescence is greatly enhanced by employing a taller chamber with saturated side walls, non‐uniform side‐wall temperature, and rough surfaces. For the conditions explored, although lowering the aerosol injection rate broadens the droplet size distribution, favoring collision‐coalescence, the reduced droplet number concentration decreases the frequency of collisions.

54 ENVIRONMENTAL SCIENCES↗

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗

Generating a stratocumulus-like cloud top in a convection-cloud chamber

Stratocumulus-topped boundary layers play a crucial role in influencing daily weather and earth energy balance. Entrainment at the stratocumulus cloud top affects the cloud’s lifetime, precipitation, and radiative properties, but our understanding remains limited due to the lack of resolution in both field observations and numerical simulations. A recently proposed convection-cloud chamber with detailed control of sidewall temperatures can provide a unique opportunity to explore this mechanism in a laboratory setting. In this work, we use numerical simulations to demonstrate that this design can produce a cloud top that mimics the entrainment interfacial layer in a stratocumulus cloud. Our results show that a steady-state cloud can be formed by cooling the lower portions of the sidewalls and warming the bottom surface, while a temperature inversion at the cloud top can be generated by keeping the upper sidewalls and top surface warmer than the bottom. The turbulent kinetic energy profile and budget are similar to those found in a convective boundary layer, and inhomogeneous mixing near the cloud top can be observed. These findings significantly enhance the scientific value of constructing the tall convection-cloud chamber.

54 ENVIRONMENTAL SCIENCES↗