High efficacy LED lamps combined with adaptive lighting control and greenhouse integrated photovoltaics (PV) could enable the concept of year-round cultivation and become a feasible option even in the harsh climate of the Nordic countries. Meteorological satellite data of this region was analyzed in a parametric study to evaluate the potential of these technologies. The generated maps showed monthly average temperatures fluctuating from -20°C to 20°C throughout the year. The natural photoperiod and light intensity also changed drastically, resulting in monthly average daily light integral (DLI) levels ranging from 45-50 mol·m-2·d-1 in summer and contrasting with 0-5 mol·m-2·d-1 during winter. To compensate, growth room cultivation independent from outdoor conditions could be used in winter. Depending on the efficacy of the lamps, the electricity required for sole-source lighting at 300 µmol·m-2·s-1 for 16 hours would be between 1.4 and 2.4 kWh·m-2·d-1. Greenhouses with supplementary lighting could help start the cultivation earlier in spring and extend it further into autumn. The energy required for lighting highly depends on several factors such as the natural light transmittance, the light threshold settings and the lighting control protocol, resulting in electric demands between 0.6 and 2.4 kWh·m-2·d-1. Integrating PV on the roof or wall structures of the greenhouse could offset some of this electricity, with specific energy yields ranging from 400 to 1120 kWh·kWp-2·yr-1 depending on the region and system design.

In this paper we obtain some inequalities of Hermite-Hadamard type for composite convex functions. Applications for AG, AH-convex functions, GA, GG, GH-convex functions and HA, HG, HH-convex function are given. Applications for p, r-convex and LogExp functions are presented as well.

The purpose of the research is to investigate the unique solvability of a nonlocal boundary value problem for a loaded equation of parabolic-hyperbolic type in a special domain. Methods: Using representations of a general regular solution, the existence and uniqueness of the problem posed is proved. Results: A new method for proving the unique solvability of nonlocal boundary value problems for one loaded equation of mixed type in a special domain has been formed and developed. Conclusions: The obtained results make it possible to solve some problems of genetics, immunology, transonic gas dynamics and thermal physics, which are associated with displacement.

In this study, we first discover some new concept coordinated UD-convex mappings with fuzzy-number values. After that, we look into Hermite-Hadamard type inequalities via fuzzy-number-valued coordinated UD-convex fuzzy-number-valued mapping (coordinated UD-convex FNVM). In the case of coordinated UD-convex FNVM, novel conclusions are derived by making particular decisions in recently proven inequalities. Additionally, it is demonstrated that the recently discovered inequalities are expansions of comparable findings in the literature. It is important to note that the main outcomes are validated by nontrivial examples.

This study presents the solution of the second type of a two-dimensional nonlinear integral equation in Banach space. Also, the existence and uniqueness of this equation’s solution are discussed. We utilize a numerical approach involving hybrid and block-pulse functions to obtain the approximate solution of a two-dimensional nonlinear integral equation. Nonlinear integral equation in two dimensions is reduced numerically to a system of nonlinear algebraic equations that can be solved using numerical methods. This study focuses on showing the convergence analysis for the numerical approach and obtaining an error estimate. Some numerical examples have been provided to demonstrate the approach’s viability and efficacy

In this paper, the notion of θ∗-weak contraction is introduced, which is utilized to prove some fixed point results. These results are helpful to give a positive response to certain open question raised by Kannan [Amer. Math. Monthly 76:1969] and Rhoades [Contemp. Math. 72:1988] on the existence of contractive definition which does not force the mapping to be continuous at the fixed point. Some illustrative examples are also given to support our results. As applications of our result, we investigate the existence and uniqueness of a solution of non-linear matrix equations and integral equations of Volterra type as well.

Recently, the conformable derivative and its properties have been introduced. In this paper, we propose and prove some new results on conformable Boundary Value Problems. First, we introduce a conformable version of classical Sturm´s separation, and comparison theorems. For a conformable Sturm-Liouville problem, Green's function is constructed, and its properties are also studied. In addition, we propose the applicability of the Green´s Function in solving conformable inhomogeneous linear differential equations with homogeneous boundary conditions, whose associated homogeneous boundary value problem has only trivial solution. Finally, we prove the generalized Hyers-Ulam stability of the conformable inhomogeneous boundary value problem.

In this work, the problem is considered for eliminating a non-uniform rectilinear smearing of an image by mathematical and computer-based way, for example, a picture of several cars taken with a fixed camera and moving with different speeds. The problem is described by a set of 1-dimensional Fredholm integral equations (IEs) of the first kind of convolution type with a 1-dimensional point spread function (PSF) at uniform smearing and by a set of new 1-dimensional IEs of a general type (not convolution type) with a 2-dimensional PSF at non-uniform smearing. The problem is also described by one 2-dimensional IE of convolution type with a 2-dimensional PSF at uniform smearing and by a new 2-dimensional IE of a general type with a 4-dimensional PSF at non-uniform smearing. The problem for solving Fredholm IE of the first kind is ill-posed (unstable). Therefore, IEs of convolution type are solved by the Fourier transform (FT) method and Tikhonov's regularization (TR), and IEs of general type are solved by the quadrature/cubature and TR methods. Moreover, the magnitude of the image smear Δ is determined by the original “spectral method”, which increases the accuracy of image restoration. It is shown that the use of a set of 1-dimensional IEs is preferable to one 2-dimensional IE in the case of non-uniform smearing. In the inverse problem (image restoration), the Gibbs effect (the effect of false waves) in the image may occur. It can be edge or inner. The edge effect is well suppressed by the proposed technique “diffusing the edges". In the case of an inner effect, it is eliminated with difficulty, but the image smearing itself plays the role of diffusing and suppresses the inner Gibbs effect to a large extent. It is shown (in the presence of impulse noise in an image) that the well-known Tukey median filter can distort the image itself, and the Gonzalez adaptive filter also distorts the image (but to a lesser extent). We propose a modified adaptive filter. A software package was developed in MatLab and illustrative calculations are performed.

In this paper, we investigate a new $q$-analogue of the higher order degenerate Changhee polynomials and numbers, which are called the Witt-type formula for the $q$-analogue of degenerate Changhee polynomials of order $r$. We can derive some new interesting identities related to the degenerate $(h,q)$-Changhee polynomials and numbers.

In this paper, we investigate a new $q$-analogue of the higher order degenerate Changhee polynomials and numbers, which are called the Witt-type formula for the $q$-analogue of degenerate Changhee polynomials of order $r$. We can derive some new interesting identities related to the degenerate $(h,q)$-Changhee polynomials and numbers.

In the paper, the authors establish two identities to express the generating function of the Chebyshev polynomials of the second kind and its higher order derivatives in terms of the generating function and its derivatives each other, deduce an explicit formula and an identities for the Chebyshev polynomials of the second kind, derive the inverse of an integer, unit, and lower triangular matrix, present several identities of the Catalan numbers, and give some remarks on the closely related results including connections of the Catalan numbers respectively with the Chebyshev polynomials, the central Delannoy numbers, and the Fibonacci polynomials.

In the full load operating range, wind turbines typically use pitch control to maintain rotor speed and output power at their rated values to guarantee the whole system’s safety. Blade pitch actuators are subject to physical limitations and safety requirements, which imposes strict constraints on pitch angle’s amplitude and rate of change. These constraints will bring windup problem to the integral part of the pitch controller. It leads to deterioration of the pitch control system’s performance or even the instability problem, especially for large-scale turbines operating in extreme uncertain wind conditions or suffering from cyber attacks. This paper designs an anti-windup robust pitch angle control strategy to deal with pitch rate constraint issue to enhance the security of control system. First, to facilitate controller design, a filtered tracking error technique is employed to transform the nonaffine form into an affine one. Subsequently, a feedback robust controller based on UDE is developed to handle the model’s uncertainty and external disturbances. To address the issue of integral saturation in the pitch system and guarantee its safety, an elliptical bounded constraint is integrated into the designed UDE strategy. This bounded UDE controller can improve the stability of power generation quality, reduces the mechanical loads on components and enhance the whole system’s safety. Finally, the effectiveness of the proposed scheme is verified on the Wind Turbine Blockset platform in Matlab/Simulink. It can achieve better performance than the traditional method.

The article considers an internal boundary value problem of the distribution of an electrostatic field in a lens formed by two identical semi-infinite circular cylinders coaxially located inside an infinite external cylinder. The problem is reduced to solving a system of singular Wiener-Hopf integral equations, which is further solved by the Wiener-Hopf method using factorized Bessel functions. Solutions to the problem for each region inside the infinite outer cylinder are presented as exponentially converging series in terms of eigenfunctions and eigenvalues.

In this paper we obtain some inequalities of Hermite-Hadamard type for composite convex functions. Applications for AG, AH-h-convex functions, GA, GG, GH-h-convex functions and HA, HG, HH-h-convex function are given.

In this paper, a multi-patch weakly singular isogeometric boundary element method (WSIGABEM) for two-dimensional elastostatics is proposed. Since the method based on the weakly singular boundary integral equation, the quadrature techniques, dedicated for the weakly singular and regular integrals, are applied in the method. A new formula for the generation of collocation points is suggested to take the full advantage of the multi-patch technique. The generated collo-cation points are essentially inside the patches without any correction. If the boundary conditions are assumed to be continuous in every patch, no collocation point lies on the discontinuous boundaries, thus simplifying the implementation. The multi-patch WSIGABEM is verified by the simple examples with analytical solutions. The feature of the present multi-patch WSIGABEM is investigated by the comparison with the traditional IGABEM. Furthermore, the combination of the present multi-patch WSIGABEM and the particle swarm optimization algorithm results in a shape optimization method in two-dimensional elastostatics. By changing some specific control points and their weights, the shape optimizations of the fillet corner, the spanner and the arch bridge are verified to be effective.

Our aim in this paper is to establish some new Hermite-Hadamard- Mercer type integral inequalities by utilizing the fractional proportional-integral operators.For this purpose, Hermite-Hadamard-Mercer inequalities for di¤er- antiable mappings whose derivatives in absolute value are convex.

This contribution considers an improved control scheme for three-phase two-stage grid-tied photovoltaic (PV) power system based on integral sliding mode control (ISMC) theory. The proposed control scheme consists of maximum power point tracking (MPPT), DC-Link voltage regulation and grid currents synchronization. A modified voltage-oriented maximum power point tracking (VO-MPPT) method based on ISMC theory is proposed for design of an enhanced MPPT under irradiation changes. Moreover, a novel DC-Link voltage control based on ISMC theory is proposed in order to achieve good regulation of DC-Link voltage over its reference. To inject the generated PV power into the grid with high quality, a voltage oriented control based on space vector modulation (SVM) and ISMC (VOC-ISMC-SVM) has been developed to control the grid currents synchronization. Numerical simulations are performed in Matlab/SimulinkTM environment in order to evaluate the proposed control strategy. In comparison with conventional control scheme, the developed control strategy provides an accurate MPP tracking with less power oscillation as well as a fast and an accurate DC-Link regulation under climatic conditions variations. Moreover, the transfer of the extracted power into the grid is achieved with high quality.

Covariant integral quantization with rotational symmetry SO(3) is established for the quantum motion on this group manifold, or, alternatively, for Gabor signal analysis on this group. We revisit the action of the related (non-unitary) Weyl-Gabor operator on the Hilbert space of square integrable functions on SO(3) and disclose a set of various properties. By selecting weight functions on the corresponding discrete-continuous phase space, one derives related coherent states and Wigner transform.

In our previous published papers, considering 3 large atomic gravitational constants assumed to be associated with weak, strong and electromagnetic interactions, we have proposed the existence of a nuclear charge of magnitude, e_n=2.95e and developed a nuclear mass formula associated with strong and weak interactions having 4 simple terms and only one energy coefficient. Two important assumptions are, there exists a weak fermion of rest energy 585 GeV and strong coupling constant is the squared ratio of electromagnetic charge and nuclear charge. The aim of this paper is associated with understanding the mystery of quantum of magnetic flux, Planck’s quantum radiation constant and Reduced Planck’s constant. Proceeding further, quark charges, strong coupling constant, nuclear stability, nuclear binding energy, medium and heavy atomic X-ray levels and celestial magnetic moments can be understood in a unified approach. It may also be noted that, by considering integral nature of elementary particle masses, it seems possible to understand the discreteness of angular momentum. Considering our proposed e_n=2.95e=3e as a characteristic nuclear charge, it seems possible to understand the integral nature of quarks electromagnetic charge. With this idea, neutron, proton and pions decay can be understood very easily. In all the cases, up quark of charge (±2e) seems to play a crucial role with the internal transformation of down quark of charge (±e) and external observable elementary basic elementary particles. It needs further study at fundamental level.

A new type of ordinary differential equation is introduced and discussed, namely, the time-dependent order ordinary differential equations. These equations can be solved via fractional calculus and are mapped into Volterra integral equations of second kind with singular integrable kernel. The solutions of the time-dependent order differential equations smoothly deforms solutions of the classical integer order ordinary differential equations into one-another, and can generate or remove singularities. An interesting symmetry of the solution in relation to the Riemann zeta function and Harmonic numbers was also proved.

The main goal of the present paper is to construct some families of the Carlitz’s q-Bernoulli polynomials and numbers. We firstly introduce the modified Carlitz’s q-Bernoulli polynomials and numbers with weight (α, β) and investigate their some explicit properties and identities arising from the bosonic q-Volkenborn integral on ℤ_p. We then define the modified degenerate Carlitz’s q-Bernoulli polynomials and numbers with weight (α, β) and obtain some recurrence relations and other identities. Moreover, we derive some correlations with the modified Carlitz’s q-Bernoulli polynomials with weight (α, β), the modified degenerate Carlitz’s q-Bernoulli polynomials with weight (α, β), the Stirling numbers of the first kind and second kind.

This paper presents a symmetrical topology for the design of solid-state transformer, made up of power switching converters, to replace conventional bulky transformers. The proposed circuitry not only reduces the overall size but also provides power flow control with the ability to be interfaced with renewable energy resources (RESs) to fulfill the future grid requirements at consumer end. Solid state transformer provides bidirectional power flow with variable voltage and frequency operation and has the ability to maintain unity power factor, and current total harmonic distortion (THD) for any type of load within defined limits of IEEE standard. Solid State Transformer offers much smaller size as compared to that of the conventional iron core transformer. MATLAB/Simulink platform is adopted to test the validity of the proposed circuit for different scenarios by providing the simulation results evaluated at 25 kHz switching frequency.

We obtain some generalised Minkowski-type integral formulas for compact Riemannian (resp. spacelike) hypersurfaces in Riemannian (resp. Lorentzian) manifolds admitting an arbitrary vector field that we assume to be timelike in the case where the ambient space is Lorentzian. Some of these formulas generalize existing formulas in the case of conformal and Killing vector fields. We apply these integral formulas to obtain interesting results concerning the characterization of such hypersurfaces in some particular cases such as when the ambient space is Einstein admitting an arbitrary (in particular, conformal or Killing) vector field, and when the hypersurface has constant mean curvature.

In this article, the multitouch option, also called the n  touch option (or the “baseball” option when n  3 ) is analyzed and valued in closed form. This is a kind of barrier option that sets a gradual knockout / knock-in mechanism based on the number of times the underlying asset has crossed a predefined barrier in various time intervals before expiry. The higher the number of predefined time intervals during which the barrier has been touched, the lower the value of a knock-out contract at expiry, and conversely for a knock-in one. Multitouch options can be viewed as an extension of step barrier options, preserving the ability of the latter to adjust the exposure to risk over time, while eliminating the notorious danger of “sudden death” that holders of step barrier options are faced with. Unlike occupation time derivatives, the payoff at expiry does not depend on the amount of time spent outside the authorized range, but on the number of passages beyond the authorized range.

In the paper we show that the general relativity in recent Einstein-Palatini formulation is equivalent to a gauge field. We begin with a bit of information of the Einstein-Palatini formulation and derive Einstein field equations from it. In the next section, we consider general relativity with a positive cosmological constant in terms of the corrected curvature. We show that in terms of the corrected curvature general relativity takes the form typical for a gauge field. Finally, we give a geometrical interpretation of the corrected curvature.

In the paper, by virtue of a Lévy–Khintchine representation and an alternative integral representation for the weighted geometric mean, the authors establish a Lévy–Khintchine representation and an alternative integral representation for the Toader–Qi mean. Moreover, the authors also collect an probabilistic interpretation and applications in engineering of the Toader–Qi mean.

Despite all efforts, the treatment of breast cancer (BC) cannot be considered to be a success story. The advances in surgery, chemotherapy and radiotherapy have not been sufficient. Indeed, the accumulated experience clearly indicates that new perspectives and non-main stream approaches are needed to better characterize the etiopathogenesis and treatment of this disease. This contibution deals with how the new pH-centric anticancer paradigm plays a fundamental role in reaching a more integral understanding of the etiology, etiopathogenesis and treatment of this multifactorial disease. For the first time the armamentarium available for the treatment of the different types and phases of BC is approached here from a Unitarian perspective based upon the hydrogen ion dynamics of cancer. The wide-ranged pH-related molecular, biochemical and metabolic model is able to embrace most the fields and subfields of breast cancer pathology. This single and integrated approach allows to advance a unidirectional program to treatment. Further efforts in this line are likely to first improve the therapeutics of each subtype of this tumor, then every phase of the disease and finally every individual patient.

Let be a power function $f_{r,M}(s)$ defined for every $s$ within the finite set $M$ as follows $$f_{r,M}(s)= \begin{cases} s^r, \ &s\in M,\\ 0, \ &\mathrm{otherwise}. \end{cases} $$ Let a discrete convolution of $f_{r,M}(s)$ be denoted as follows $\mathrm{Conv}_{r,M}[n]=(f_{r,M}*f_{r,M})[n]$. Let a real coefficients $A_{m,j}$ be given by the following recurrence $$ A_{m,j} = \begin{cases} 0, & \mathrm{if } \ j<0 \ \mathrm{or } \ j>m, \\ (2j+1)\binom{2j}{j} \sum_{d=2j+1}^{m} A_{m,d} \binom{d}{2j+1} \frac{(-1)^{d-1}}{d-j} B_{2d-2j}, & \mathrm{if } \ 0 \leq j < m, \\ (2j+1)\binom{2j}{j}, & \mathrm{if } \ j=m. \end{cases} $$ In this paper we show that for every $n>0$ the following odd-power identities involving coefficients $A_{m,j}$ and convolution transform $\mathrm{Conv}_{r,M}[n]$ hold $$ \begin{split} n^{2m+1}+1&=\sum_{r=0}^{m}A_{m,r}\mathrm{Conv}_{r,\mathbb{N}}[n],\\ n^{2m+1}-1&=\sum_{r=0}^{m}A_{m,r}\mathrm{Conv}_{r,\mathbb{Z}_{>0}}[n],\\ n^{2m+1}&=\sum_{r=0}^{m}A_{m,r}\sum_{k=1}^{n} k^r(n-k)^r\\ &=\sum_{r=0}^{m}A_{m,r}\sum_{k=0}^{n-1} k^r(n-k)^r. \end{split} $$

An advanced driver-assistance system (ADAS) is critical to driver-vehicle-interaction systems. Driving behaviour modelling and control significantly improves the global performance of ADAS. A haptic assistive system assists the driver by providing a specific torque on the steering wheel according to the driving-vehicle-road profile to improve the steering control. However, the main problem is designing a compensator dealing with the high-level uncertainties in different driving scenarios with haptic driver assistance, where different personalities and diverse perceptions of drivers are considered. If not properly accounted for, these differences can lead to poor driving performance. This paper focuses on designing a data-driven model-free compensator considering various driving behaviours with a haptic feedback system. A back propagation neural network (BPNN) models driving behaviour based on real driving data (speed, acceleration, vehicle orientation, and current steering angle). Then using the genetic algorithm (GA) and the integral time absolute error (ITEA) criterion to optimize multiple PID compensation parameters for various driving behaviour (such as speeding/braking, lane-keeping and turning), which are then be combined by the fuzzy logic to provide different driving commands. An experiment was conducted with 5 participants in a driving simulator. During the second experiment, seven participants drove in the simulator to evaluate the robustness of the proposed combined GA-PID offline and the Fuzzy-PID controller applied online. The experiment results evaluated the ITAE of lateral displacement and yaw angle, during various driving behaviour. The results validated the proposed method by significantly enhancing the driving performance.

The aim of this article is to establish a new class of unified integrals associated with the generalized Lommel-Wright functions, which are expressed in terms of Wright Hypergeometric function.Some integrals involving trigonometric,generalized Bessel function and Struve functions are also indicated as special cases of our main results.further, with the help of our main results and their special cases, we obtain two reduction formulas for the Wright hypergeometric function.

In the paper, by induction and recursively, the author proves that the generating function of multivariate logarithmic polynomials and its reciprocal are a Bernstein function and a completely monotonic function respectively, establishes a Lévy-Khintchine representation for the generating function of multivariate logarithmic polynomials, deduces an integral representation for multivariate logarithmic polynomials, presents an integral representation for the reciprocal of the generating function of multivariate logarithmic polynomials, computes real and imaginary parts for the generating function of multivariate logarithmic polynomials, derives two integral formulas, and denies the uniform convergence of a known integral representation for Bernstein functions.

We introduce the Symplectic Structure of Information Geometry based on Souriau’s Lie Group Thermodynamics model, with a covariant definition of Gibbs equilibrium via invariances through co-adjoint action of a group on its moment space, defining physical observables like energy, heat, and moment as pure geometrical objects. Using Geometric (Planck) Temperature of Souriau model and Symplectic cocycle notion, the Fisher metric is identified as a Souriau Geometric Heat Capacity. Souriau model is based on affine representation of Lie Group and Lie algebra that we compare with Koszul works on G/K homogeneous space and bijective correspondence between the set of G-invariant flat connections on G/K and the set of affine representations of the Lie algebra of G. In the framework of Lie Group Thermodynamics, an Euler-Poincaré equation is elaborated with respect to thermodynamic variables, and a new variational principal for thermodynamics is built through an invariant Poincaré-Cartan-Souriau integral. The Souriau-Fisher metric is linked to KKS (Kostant-Kirillov-Souriau) 2-form that associates a canonical homogeneous symplectic manifold to the co-adjoint orbits. We apply this model in the framework of Information Geometry for the action of an affine Group for exponentiel families, and provide some illustrations of use cases for multivariate Gaussian densities. Information Geometry is presented in the context of seminal work of Fréchet and his Clairaut-Legendre equation. Souriau model of Statistical Physics is validated as compatible with Balian gauge model of thermodynamics. We recall the precursor work of Casalis on affine group invariance for natural exponential families.

In this paper, the 1st level general fractional derivatives of arbitrary order are defined and investigated for the first time. We start with a generalization of the Sonin condition for the kernels of the general fractional integrals and derivatives and then specify a set of the kernels that satisfy this condition and posses an integrable singularity of power law type at the origin. The 1st level general fractional derivatives of arbitrary order are integro-differential operators of convolution type with the kernels from this set. They contain both the general fractional derivatives of arbitrary order of the Riemann-Liouville type and the regularized general fractional derivatives of arbitrary order considered in the literature so far. For the 1st level general fractional derivatives of arbitrary order, some important properties including the 1st and the 2nd fundamental theorems of Fractional Calculus are formulated and proved.

A new framework in quantum chemistry has been proposed recently ("An approach to first principles electronic structure calculation by symbolic-numeric computation'' by A. Kikuchi). It is based on the modern technique of computational algebraic geometry, viz. the symbolic computation of polynomial systems. Although this framework belongs to molecular orbital theory, it fully adopts the symbolic method. The analytic integrals in the secular equations are approximated by the polynomials. The indeterminate variables of polynomials represent the wave-functions and other parameters for the optimization, such as atomic positions and contraction coefficients of atomic orbitals. Then the symbolic computation digests and decomposes the polynomials into a tame form of the set of equations, to which numerical computations are easily applied. The key technique is Grröbner basis theory, by which one can investigate the electronic structure by unraveling the entangled relations of the involved variables. In this article, at first, we demonstrate the featured result of this new theory. Next, we expound the mathematical basics concerning computational algebraic geometry, which are necessitated in our study. We will see how highly abstract ideas of polynomial algebra would be applied to the solution of the definite problems in quantum mechanics. We solve simple problems in "quantum chemistry in algebraic variety'' by means of algebraic approach. Finally, we review several topics related to polynomial computation, whereby we shall have an outlook for the future direction of the research.

Recent technological advancements in produced water for the oil and gas industry have a long history that spreads over centuries and has roots with many mathematicians, physicists, engineers, and inventors. Working to understand the many ideas and questions has helped to bring the produced water technologies to where we are today. Fluid mechanics is an essential factor in the world of engineering because it comes with significant benefits when designing the necessary equipment to overcome the challenging endeavors of daily operations. Many of us believe certain ideas or products were developed within a specific year not fully understanding the idea or concept might be hundreds of years old. For example, when we think of fluid mechanics, we do not think of Archimedes' immersing King Hiero II crown in water. Archimedes did this to determine the displacement of the water as he knew gold was denser than silver. Archimedes was one of the pioneers who was part of the beginning of fluid mechanics. Another interesting fact is the first valve invented was in 1 AD by Heron of Alexandria. This valve can be described as the very first hydraulic valve. The evolution of technology over hundreds of years has made considerable innovations benefiting engineers all over the world. These innovations come in the form of PLCs, RTUs, PIDs, mathematical formulas, processing methods, and SCADA systems.