Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “residual minimization”

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 55 records · Page 3

Regional Seismic Travel Time v.3.2.0

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗

Studying the effects of overlapping objects in dark energy

Observing the clustering of galaxies allows us to calculate cosmological parameters necessary for understanding dark energy. However, as the density of observed objects increases, the probability of these objects blending likewise increases, causing multiple galaxies to be observed as one. This affects the inferred values of parameters such as the galaxy bias (b) and the matter energy density (ΩM). To see whether the bias from incorrectly inferring the galaxy count is significant, we compare the correlation functions in simulated data for “true” and “observed” data sets with one-to-one and multiple-to-one correspondences, respectively. For each data set, we create two correlation functions: one “measured” function directly relying on the galaxies’ positions using the TreeCorr python library, and one “model” derived mathematically from the galaxies’ power spectrum using the Cosmological Core Library (CCL). By minimizing the residual between these two functions, we compute the ideal values for b and ΩM across the various possible redshifts that position the galaxies in three dimensional space. This minimization is done with an Markov chain Monte Carlo (MCMC) estimate that finds one value of ΩM and ten values for b corresponding to the ten redshift bins ranging from z = 0.2 to z = 1.2. We find that neither b nor ΩM is particularly affected by inclusion of blended galaxies. Though there is room for improvement, the data suggests that the fluctuations we found are a result of noise or limitations on the modeling rather than blending explicitly.

79 ASTRONOMY AND ASTROPHYSICS↗

Simultaneous optimization of hatching orientations and lattice density distribution for residual warpage reduction in laser powder bed fusion considering layerwise residual stress stacking

The reduction in the residual warpage generated via fabrication is an emerging issue in metals’ laser powder bed fusion additive manufacturing (AM). One solution from a design for additive manufacturing (DfAM) perspective is to arrange a hatching strategy that utilizes the asymmetry of the residual stress generated from a molten pool. Another approach is to arrange the geometry of the parts, such as forming an inner lattice structure to arrange the effective stiffness. In this study, the simultaneous optimization of the laser hatching orientation and lattice density distribution was conducted to minimize the residual warpage in lattice infill structures and confirm their synergetic effect. The mutual effect of each layer was highlighted when considering the mechanical background of the optimal results because an imbalance in the residual stress caused the warpage of the AM parts between the layers. The proposed methodology was validated via experimental verification using a quasi-2D plate, 3D bracket, and 3D connecting-rod design problems. The residual warpage in the quasi-2D study, measured as the edge’s vertical deformation, was reduced from 23.4% to 39.4% compared with the benchmark design. In the 3D studies, the warpage reductions ranged from 13.1% to 20.7%.

36 MATERIALS SCIENCE↗

Residual elastic strain evolution due to thermal cycling of a ceramic-metal composite (WC-Cu) via high energy X-ray diffraction and analytical modeling

Residual stress, when superimposed with in-service loading, can significantly reduce the lifetime and performance of a component. Ceramic-metal composites are susceptible to residual stresses due to the thermal expansion mismatch of the ceramic and metallic phases. The WC-Cu composite explored in the present study provides a promising combination of thermal conductivity and strength properties, while exhibiting counterintuitive improvements in strength and ductility after thermal cycling. Further, this work quantifies the evolution of the residual elastic strains as a result of processing and cyclic thermal loading in a co-continuous WC-Cu composite through experimental high energy X-ray diffraction and kinetics-based modeling. Both analyses indicate that processing-induced residual tensile stress in the copper phase is relieved upon subsequent thermal cycling, with kinetics modeling revealing the cyclic-dependent nature of the active power-law creep mechanisms. The results indicate that, through stress relaxation, this material system maintains structural stability during thermal cycling. The illustrated kinetics of relaxation can inform general material processors and designers of ceramic-metal composites to minimize detrimental residual stress and improve performance of these material systems.

36 MATERIALS SCIENCE↗

Frhodo: A program for simulating chemical kinetic measurements and optimizing kinetic mechanisms

Chemical kinetic simulations are frequently used to extract kinetic and mechanistic information from experimental data by fitting simulations to the data. For simple systems of consecutive single channel reactions this is a relatively straightforward task. However, as the complexity of the system increases and particularly for systems involving competitive multi-channel reactions manually optimizing a mechanism to obtain the best fits to a range of experimental data quickly becomes a time-intensive and challenging process, especially when rate coefficients have pressure and temperature dependencies. There is considerable potential for exploiting automated optimization methods to rapidly screen reaction mechanisms and optimize the fit by adjusting rate coefficients within well-defined constraints. A new chemical kinetics simulation program, Frhodo, has been developed for this purpose. The program allows either manual simulation of individual experiments or fully automated optimization against a range of experimental data by adjusting user selected reactions within user defined constraints. Frhodo incorporates machine learning-based optimization routines that work by either minimizing standardized residuals or through Bayesian parameter estimation. Frhodo's optimization capabilities are demonstrated by reexamining two previously studied systems of relevance to combustion. Dissociation of diacetyl followed by recombination of methyl radicals is an example of a sequence of single channel, consecutive reactions. Pyrolysis of 2-methyl furan exemplifies multi-channel, unimolecular reactions that are often encountered and have temperature and pressure dependencies in the rate coefficients of the channels and branching ratios between them. For both systems the optimization routines resulted in solutions similar to the original studies. Finally, comments are made on the strengths and limitations of the approach.

0D simulation↗

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. Furthermore, we demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion–reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.

Data-driven model reduction↗

Record-Low Thermal Boundary Resistance between Diamond and GaN-on-SiC for Enabling Radiofrequency Device Cooling

The implementation of 5G-and-beyond networks requires faster, high-performance, and power-efficient semiconductor devices, which are only possible with materials that can support higher frequencies. Gallium nitride (GaN) power amplifiers are essential for 5G-and-beyond technologies since they provide the desired combination of high frequency and high power. These applications along with terrestrial hub and backhaul communications at high power output can present severe heat removal challenges. The cooling of GaN devices with diamond as the heat spreader has gained significant momentum since device self-heating limits GaN’s performance. However, one of the significant challenges in integrating polycrystalline diamond on GaN devices is maintaining the device performance while achieving a low diamond/GaN channel thermal boundary resistance. In this study, we achieved a record-low thermal boundary resistance of around 3.1 ± 0.7 m 2 K/GW at the diamond/Si 3 N 4 /GaN interface, which is the closest to theoretical prediction to date. The diamond was integrated within ~1 nm of the GaN channel layer without degrading the channel’s electrical behavior. Furthermore, we successfully minimized the residual stress in the diamond layer, enabling more isotropic polycrystalline diamond growth on GaN with thicknesses >2 μm and a ~1.9 μm lateral grain size. More isotropic grains can spread the heat in both vertical and lateral directions efficiently. Using transient thermoreflectance, the thermal conductivity of the grains was measured to be 638 ± 48 W/m K, which when combined with the record-low thermal boundary resistance makes it a leading-edge achievement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics-Informed Neural Network Method for Forward and Backward Advection-Dispersion Equations

Advection-dispersion equations (ADEs) are commonly used to describe transport phenomena in porous media. Even though mature discretization-based numerical methods for ADEs exist, some challenges still remain, especially when it comes to solving advection-dominated forward ADEs and diffusion-dominated backward ADEs. The latter problem usually arises in the source identification context and leads to numerically unstable grid-based solutions that require a form of regularization or should be treated as an inverse problem that is computationally more expensive because it requires solving the forward problem multiple times. In this study, we propose a discretization-free approach based on the physics-informed neural network (PINN) method for solving coupled ADE and Darcy flow equations with space-dependent hydraulic conductivity. In this approach, the hydraulic conductivity, hydraulic head, and concentration fields are approximated with deep neural networks (DNNs). We assume that the conductivity field is given by its values on a grid, and we use these values to train the conductivity DNN. The head and concentration DNNs are trained by minimizing the residuals of the flow equation and ADE and using the initial and boundary conditions as additional constraints. The PINN method is applied to one- and two-dimensional forward ADE problems, where its performance for various P\'{e}clet numbers ($Pe$) is compared with the analytical and numerical solutions. We find that the PINN method is accurate with errors of less than 1\% and outperforms some conventional discretization-based methods for $Pe$ larger that 100. Next, we demonstrate that the PINN method remains accurate for the backward ADEs, with the relative errors in most cases staying under 5\% compared to the reference concentration field. Finally, we show that when available, the concentration measurements can be easily incorporated in the PINN method and significantly improve (by more than 50\% in the considered cases) the accuracy of the PINN solution of the backward ADE.

He, Qizhi↗

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-tool optimization for computer controlled optical surfacing

With the rapid development of precision technologies, the demand of high-precision optical surfaces has drastically increased. These optical surfaces are mainly fabricated with computer controlled optical surfacing (CCOS). In a CCOS process, a target surface removal profile is achieved by scheduling the dwell time for a set of machine tools. The optimized dwell time should be positive and smooth to ensure convergence to the target while considering CNC dynamics. The total run time of each machine tool is also expected to be balanced to improve the overall processing efficiency. In the past few decades, dwell time optimization for a single machine tool has been extensively developed. While the methods are applicable to multi-tool scenarios, they fail to consider the overall contributions of multiple tools simultaneously. In this paper, we conduct a systematic study on the strategies for multi-tool dwell time optimization and propose an innovative method for simultaneously scheduling dwell time for multiple tools for the first time. First, the influential factors to the positiveness and smoothness of dwell time solutions for a single machine tool are analyzed. The compensation strategies that minimize the residual while considering the CNC dynamics limit are then proposed. Afterwards, these strategies are extended to the proposed multi-tool optimization that further balances the run time of machine tools. Finally, the superiority of each strategy is carefully studied via simulation and experiment. The experiment is performed by bonnet polishing a 60 mm × 60 mm mirror with three tools of different diameters ( i.e., 12 mm, 8 mm, and 5 mm). The figure error of the mirror is reduced from 45.42 nm to 11.18 nm root mean square in 13.28 min. Moreover, the measured polishing result well coincides with the estimation, which proves the effectiveness of the proposed method.

36 MATERIALS SCIENCE↗

Age at release affects developmental physiology and sex-specific phenotypic diversity of hatchery steelhead trout (Oncorhynchus mykiss)

Most steelhead trout hatcheries increase growth rate during rearing to produce and release yearling smolts for harvest augmentation, but natural steelhead exhibit variable age of smoltification, so this common rearing practice may not be ideal for programs focused on recovering imperiled wild stocks; therefore, it is important to investigate and compare alternative hatchery rearing methods that promote life history diversity. Over six consecutive years, the Winthrop National Fish Hatchery on the Methow River, WA reared and released paired groups of age-1 (S1) and age-2 (S2) steelhead smolts. To understand how the two rearing methods affected developmental ontogeny and life-history, fish were sampled prior to hatchery release for factors associated with smoltification (size, gill Na+/K+ ATPase activity, and a qualitative smolt phenotype) and sexual maturation (sex, pituitary and testis mRNA transcripts, gonadosomatic index, and plasma 11-ketotestosterone). Our objectives were to quantify levels of smoltification and male maturation during hatchery rearing, combine metrics to estimate residualism (failure to migrate upon release), and compare the treatments by sex. Overall, S2 rearing produced 7.8% more smolts and 44-fold (4.4 vs. 0.1%) more precociously mature males than S1 rearing. Conversely, S1 rearing produced 31.6% more residuals than S2 rearing. While the proportion of total male residuals was comparable between treatments, the S1 treatment produced approximately five-fold more female residuals (20.6 vs. 4.2%). Because residuals contribute minimally to adult returns and the number of returning adult females is critical to the success of salmonid supplementation efforts, developing rearing techniques that maximize migration in females is a management priority. Physiological assessments are useful for characterizing and quantifying the effects and risks of different hatchery rearing regimes on steelhead life-history, in addition to providing sex-specific guidance to inform and optimize conservation management goals in supplementation programs.

Middleton, Mollie A. (ORCID:0009000905577865)↗

MULTISTEP AND CONTINUOUS PHYSICS-INFORMED NEURAL NETWORK METHODS FOR LEARNING GOVERNING EQUATIONS AND CONSTITUTIVE RELATIONS

In this work, we investigate the applicability and relative merit of discrete and continuous versions of physics-informed neural network (PINN) methods for learning unknown governing equations or constitutive relations in a nonlinear dynamical system. In the case of unknown dynamics, entire right-hand-side (RHS) equations of the ordinary differential equations are unknown. In the case of unknown constitutive relations, however, the RHS equations are known up to the specification of constitutive relations (that may depend on the state of the system). We use a deep neural network to model unknown governing equations or constitutive relations. The discrete PINN approach combines classical multistep discretization methods for dynamical systems with neural-network-based machine learning methods. On the other hand, the continuous versions utilize deep neural networks to minimize the residual function for the continuous governing equations. We use the case of a fedbatch bioreactor system to study the effectiveness of these approaches and discuss conditions for their applicability. Our results indicate that the accuracy of the trained neural network models is much higher for the cases where we only have to learn a constitutive relation instead of all dynamics. This finding corroborates the well-known fact from scientific computing that building as much structural information as is available into an algorithm can enhance its efficiency and/or accuracy.

97 MATHEMATICS AND COMPUTING↗

Molecular mechanism of secreted amyloid-β precursor protein in binding and modulating GABABR1a

A recent phenomenal study discovered that the extension domain of secreted amyloid-β precursor protein (sAPP) can bind to the intrinsically disordered sushi 1 domain of the γ-aminobutyric acid type B receptor subunit 1a (GABA B R1a) and modulate its synaptic transmission. The work provided an important structural foundation for the modulation of GABA B R1a; however, the detailed molecular interaction mechanism, crucial for future drug design, remains elusive. Here, we further investigated the dynamical interactions between sAPP peptides and the natively unstructured sushi 1 domain using all-atom molecular dynamics simulations, for both the 17-residue sAPP peptide (APP 17-mer) and its minimally active 9 residue segment (APP 9-mer). We then explored mutations of the APP 9-mer with rigorous free energy perturbation (FEP) calculations. Our in silico mutagenesis studies revealed key residues (D4, W6, and W7) responsible for the binding with the sushi 1 domain. More importantly, one double mutation based on different vertebrate APP sequences from evolution exhibited a stronger binding (ΔΔG = -1.91 ± 0.66 kcal mol -1 ), indicating a potentially enhanced GABA B R1a modulator. These large-scale simulations may provide new insights into the binding mechanism between sAPP and the sushi 1 domain, which could open new avenues in the development of future GABA B R1a-specific therapeutics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geometrical effects on the measured neutron residual stress in additively manufactured cold spray AA6061

Aluminum cold spray additive builds were made with varying wall widths and geometries to explore the effects of build geometry on residual stress development in the additively sprayed material. Deposition rates of over 750 g/hr were achieved without any cracking or delamination of the additively manufactured components. Neutron residual stress measurements revealed a maximum tensile residual stress of 41 MPa at the interface of the substrate and a maximum compressive residual stress of −35 MPa in the cold sprayed material. While there was not found to be significant residual stress changes between the cold sprayed walls, the cold sprayed cylinder had varying residual stress throughout the diameter with the outer diameter having a higher tensile residual stress than the interior. The low residual stress and minimal geometrical dependency demonstrate further part size is possible without concerns for component failure in AA6061. The application of common coating residual stress models to full scale additive builds is evaluated and discussed. The lack of residual stress evolution was studied by in-situ coating property measurement that revealed that the deposition stress of the AA6061 sprays reduced with increased layer thickness. In conclusion, this study provides insight into the evolution of residual stress in cold spray additive components and how both the mechanics and material responses differ from coatings.

Additive manufacturing↗

Performance of pumped counterflow virtual impactors to study aerosol interactions with laboratory generated warm clouds

Pumped Counterflow Virtual Impactors (PCVI) are designed to separate aerosols based on aerodynamic diameter, which is particularly useful for isolating cloud droplets and ice crystals from smaller particles. However, the PCVI transmission efficiency (TE) values reported in the literature show considerable variability, and little information is available on the TEs for cloud liquid droplets. Here, we determined the optimal flow conditions for PCVI sampling for different activation ratios from the MTU Pi-cloud chamber, highlighting the conditions that maximize droplet residual sampling while minimizing interstitial transmission. Even for a lower limit cloud activation ratio = 1:10, an add flow of 1.5 LPM achieves a residual fraction greater than 0.80, with a droplet TE of ∼15% for typical Pi-Chamber droplets with diameters between 3.5 to 10 µm. This framework can be adapted for field cloud measurements based on the cloud conditions of interest. TE for cloud droplets under flow conditions at which ∼99% of the unactivated submicron particles were removed was substantially lower between 8 and 16% than for supermicron dry particles (∼30 to 40% for 3 µm polystyrene latex spheres). In addition, we conducted a comparison of the performances of three similar PCVI units to assess repeatability for dry aerosols and cloud droplets. All three PCVIs performed similarly in rejecting submicron particles smaller than the desired cutoff diameter. One of the units showed lower TEs due to a misalignment of the internal orifice. We discuss a procedure that improved alignment and performance.

54 ENVIRONMENTAL SCIENCES↗

Consumable development to tailor residual stress in parts fabricated using directed energy deposition processes

Distortion and residual stresses are major challenges that limit the ability to fabricate large scale structures using Additive Manufacturing (AM). Researchers worldwide are evaluating techniques to induce compressive residual stress in the parts via intermittent rolling. While reasonable success has been documented, the idea of lowering the martensite start temperature to induce compressive stresses has not been evaluated in the context of AM, despite demonstrated success by the welding community. This study validates the hypothesis that, by a proper selection of materials and process parameters, one may effectively reduce distortion and induce a compressive residual stress in AM parts. Using neutron diffraction to measure residual stresses in parts, we demonstrate that, in addition to selection of the correct materials, the inter pass temperature plays a major role in controlling the residual stress evolution. The observations relating to the residual stresses are rationalized based on a microstructural evolution in these samples. Based on this preliminary study, a strategy to fabricate large structures with minimal distortion and residual stress is outlined.

36 MATERIALS SCIENCE↗

Mechanism-guided engineering of a minimal biological particle for genome editing

The widespread application of genome editing to treat and cure disease requires the delivery of genome editors into the nucleus of target cells. Enveloped delivery vehicles (EDVs) are engineered virally derived particles capable of packaging and delivering CRISPR-Cas9 ribonucleoproteins (RNPs). However, the presence of lentiviral genome encapsulation and replication proteins in EDVs has obscured the underlying delivery mechanism and precluded particle optimization. Here, we show that Cas9 RNP nuclear delivery is independent of the native lentiviral capsid structure. Instead, EDV-mediated genome editing activity corresponds directly to the number of nuclear localization sequences on the Cas9 enzyme. EDV structural analysis using cryo-electron tomography and small molecule inhibitors guided the removal of ~80% of viral residues, creating a minimal EDV (miniEDV) that retains full RNP delivery capability. MiniEDVs are 25% smaller yet package equivalent amounts of Cas9 RNPs relative to the original EDVs and demonstrated increased editing in cell lines and therapeutically relevant primary human T cells. These results show that virally derived particles can be streamlined to create efficacious genome editing delivery vehicles with simpler production and manufacturing.

59 BASIC BIOLOGICAL SCIENCES↗