In this work we study the recently developed parametrized partition function formulation and show how we can infer the thermodynamic properties of fermions based on numerical simulation of bosons and distinguishable particles at various temperatures. In particular, we show that in the three dimensional space defined by energy, temperature and the parameter characterizing parametrized partition function, we can map the energies of bosons and distinguishable particles to fermionic energies through constant-energy contours. We apply this idea to both noninteracting and interacting Fermi systems and show it is possible to infer the fermionic energies at all temperatures, thus providing a practical and efficient approach to obtain thermodynamic properties of large fermion systems with numerical simulation. As an example, we present energies for up to 50 noninteracting fermions and up to 20 interacting fermions at all temperatures and show good agreement with the analytical result for noninteracting case.

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 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.

TThis paper delves into the extension and characterization of radial positive definite functions into non-integer dimensions. We provide a thorough investigation by employing the Riemann-Liouville fractional integral and fractional Caputo derivatives, enabling a comprehensive understanding of these functions. Additionally, we introduce a secondary characterization based on the Bernstein characterization of completely monotone functions. The practical significance of our study is showcased through an examination of the positivity of the fundamental solution of the space-fractional Bessel diffusion equation, highlighting the real-world applicability of the developed concepts. Through this work, we contribute to the broader understanding of radial positive definite functions and their utility in diverse mathematical and applied contexts.

In this paper, we consider a pursuit-evasion game of inertial players, where pursuer’s control is subject to integral constraint, and evader’s control is subject to geometric constraint. In the pursuit problem, the main tool is the strategy of parallel pursuit. Sufficient conditions are obtained for the solvability of pursuit-evasion problems. Also, the main lemma describing the monotonicity of an attainability domain of the evader is proved and an explicit analytical formula for this domain is given. As one of the important results of our paper, the Isaacs Life-line game is solved for a special case.

Abstract: — management of the Continuous Stirred Tank Reactor (CSTR) operation has proven to be a difficult undertaking due to the complex and ever-changing nature of its behavior, which is both nonlinear and dynamic, which is further compounded by its instability. Attempts have been made to use Proportional-Integral-Derivative (PID) controllers for such processes, but the performance is often poor due to the inherent limitations of the controllers. In contrast, heuristic algorithms have been employed to tune the controllers, but these methods can only provide optimal performance for a restricted range of parameter variations. This paper presents a novel approach that involves using a conventional gain scheduling PID controller in combination with a fuzzy gain scheduling (FGS) controller method to calibrate set-point weighted PID controllers for the CSTR process. The paper outlines the design of the controllers and evaluates their performance by subjecting them to varying input signals within the given bandwidth.

The design of a mixed-mode proportional-integral-derivative (PID) controller circuit using current-feedback operational amplifiers (CFOAs) as active components is proposed. With the same circuit topology, the proposed configuration of three CFOAs, four resistors, and two capacitors is capable of performing the PID controller in all four operation modes: namely voltage mode, trans-admittance mode, current mode, and trans-impedance mode. Numerous mathematical analyses are conducted to determine the controller performance under ideal and non-ideal conditions. Additionally, the mixed-mode second-order lowpass filter is suggested, and also used to examine the workability of the proposed mixed-mode PID controller in a feedback control structure. The proposed PID controller is implemented with the commercially available IC-type CFOA AD844, and the simulation results are presented to illustrate the functionality of the controller and its closed-loop control system.

Ballasted railway tracks can be modelled using reduced/simplified models composed of several layers of discrete components. This paper deals with the two-layer model, which is very popular due to its computational efficiency. In order to provide some recommendations for track design, it is necessary to identify which set of parameters leads to some irregular/unexpected behavior. Such irregularity is investigated at three levels, such as: (i) critical velocity of moving constant force; (ii) instability of a single moving mass; (iii) instability of two moving masses. All results are presented in dimensionless form to cover a wide range of real parameters. Irregular cases are identified by sets of parameters leading to them and then general conclusions are drawn. Regarding the method, all results are obtained analytically or semianalytically, where “semi” refers to solving the roots of a given polynomial by predefined numerical procedures in symbolic software. No numerical integration is involved in all results presented. This means that the results are highly accurate and refer to exact values, so any kind of parametric or sensitivity analyses is readily possible.

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.

We prove Lp estimates of a class of generalized Marcinkiewicz integral operators with mixed homogeneity on product domains. By using these estimates along with an extrapolation argument we obtain the boundedness of our operators under very weak conditions on the kernel functions. Our results in this paper improve and extend several known results on both generalized Marcinkiewicz integrals and parabolic Marcinkiewicz integrals on product domains.

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 the paper, a study of the performance of unmanned ground vehicle (UGV), establishes numerical simulation model of the chassis system by three-dimensional software and finite element analysis software. Based on the numerical modelling and analysis, the work focuses on study of the strength and stiffness of the vehicle's chassis and researched the distribution of stress and deformation under the extremal condition. The paper also presents an autopilot design with new cascade control system for an autonomous motion of an unmanned ground vehicle based on Proportional-Integral-Derivative (PID) and Feed-Forward (FF) control. The PID-FF controller is a part of dedicated for presented UGV used in a hybrid control system for precise control and stabilization, necessary to increase the vehicle motion stability and maneuvers precision. The hybrid PID-FF control system proposed for the ground vehicle model gives a satisfactory quality of control, while maintaining the simplicity of the control system. Presented tests performed in area of mechanical design and control analysis gives good results and proves usefulness of designed unmanned device.

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.

In the paper, we develop some numerical schemes including finite difference scheme and finite volume scheme for the fractional Laplacian operator, and apply the resulted numerical schemes to solve some fractional diffusion equations. First, the fractional Laplacian operator can be characterized as the weak singular integral by a integral operator with zero boundary condition. Second, because of the solution of fractional diffusion equation is usually singular near the boundary, we apply finite difference scheme to discrete fractional Laplacian operator and fractional diffusion equation by fractional interpolation function. Moreover, it is find that the differential matrix of the above scheme is symmetric matrix, and strictly row-wise diagonally dominant in special fractional interpolation function. Third, we show finite volume scheme to discrete fractional diffusion equation by fractional interpolation function, and analyze the properties of differential matrix. Finally, the numerical experiments are given, and verify the correctness of theoretical results and the efficiency of the schemes.

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 contemporary construction, the prevalence of vibration serviceability issues in lightweight and slender structures has become increasingly common, owing to advancements in building materials and construction methods. While these structures often meet the criteria for ultimate limit states, they can still elicit complaints due to excessive vibrations induced by human activity. To address this challenge, the Integral Resonant Control (IRC) technique has emerged as a favored approach for actively damping vibrations in various systems. This study introduces a fresh perspective on implementing a multi-input multi-output (MIMO) IRC scheme for active vibration control (AVC) specifically tailored for pedestrian structures utilizing inertial mass actuators. Building upon a common framework and design methodology outlined in previous research, this work presents a novel application of MIMO IRC for AVC. The designed controller is rigorously tested and implemented on a laboratory floor structure to validate its efficacy.

Amplitudes of quantum transitions containing time zigzags are considered. The discussion is carried out in the framework of the Minkowski metric and the standard quantum mechanics without adding new postulates. We argue that time zigzags are not suppressed at the quantum level, but their contribution to the amplitude is zero. The result is valid for a single particle and a non-interacting scalar field.

We consider a classial case of irrational integrals containing a square root of a quadratic polynomial. It is well known that they can be expressed in terms of elementary functions by one of three Euler’s substitutions. It is less known that the Euler substittutions have a beautiful geometric interpretation. In the framework of this interpretation one can see that the number 3 is not the most suitable. We show that it is natural to introduce the fourth Euler substitution. By the way, it is not clear who was the first to attribute these three substitutions to Euler. In his original treatise Leonhard Euler uses two substitutions which are sufficient to cover all cases.

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.

We delve into nonrelativistic quantum electrodynamics within a background magnetic field, ensuring gauge invariance via a vector potential. Our exploration extends to the Lagrangian, which incorporates electron self-interactions and electromagnetic field interactions. Through the path integral formalism, we elucidate the effective action, with a specific focus on the photon propagator and screening effects. Examining density and current components reveals Laurent expansions in the hydrodynamic limit. The equations of motion in real space are derived, leading to discussions on phenomena such as Poiseuille-like flow and the Navier-Stokes equation. To ground our theoretical framework at practical contexts, we consider applications such as the fluid dynamics at binary star systems. Depending on the velocity of the expelled fluid, our analysis allows for both nonrelativistic and relativistic treatments, providing a versatile tool for understanding the intricacies of fluid behavior at astrophysical scenarios. Additionally, we recognize that at the presence of a magnetic field, magnetohydrodynamics becomes crucial. While our current focus is at nonrelativistic quantum electrodynamics, the insights gained here contribute to a broader understanding, offering a comprehensive foundation for quantum electrodynamics analysis at diverse physical systems.

Background: Previous papers have developed a statistical mechanics of neocortical interactions (SMNI) fit to short-term memory and EEG data (Ingber, 2018). Adaptive Simulated Annealing (ASA) was used for all fits to data. A numerical path-integral for quantum systems, qPATHINT, was used. Objective: The quantum path-integral for Calcium ions was used to derive a closed-form analytic solution at arbitrary time. The quantum effects is parameterized here, whereas the previous 2018 paper applied a nominal ratio of 1/2 to these effects. Method: Methods of mathematical-physics for optimization and for path integrals in classical and quantum spaces are used. The quantum path-integral is used to derive a closed-form analytic solution at arbitrary time, and is used to calculate interactions with classical-physics SMNI interactions among scales. Results: The mathematical-physics and computer parts of the study are successful, in that three cases with Subjects (blind to this author) after 1,000,000 visits to the cost function gave: Subject-07 = 0.04, Subject-08 = 0.55, and Subject-09 = 1.00. All other 9 Subjects gave 0.

Ship wave is of great interest for wave drag, coastal erosion and ship detection. In this paper, a highly-paralleled numerical scheme is proposed for simulating three-dimensional (3-D) nonlinear Kelvin ship waves effectively. First, a numerical model for nonlinear ship waves is established based on potential flow theory, boundary integral method and Jacobian-free Newton-Krylov (JFNK) method. In order to improve computational efficiency and reduce data storage of JFNK method, a banded preconditioner method is then developed by formulating the optimal bandwidth selection rule. After that, a Graphics Process Unit (GPU) based parallel computing framework is designed, and a GPU solver is developed by using Compute Unified Device Architecture (CUDA) language. Finally, numerical simulations of 3-D nonlinear ship waves under multiple scales are performed by using the GPU and CPU solvers. Simulation results show that the proposed GPU solver is more efficient than the CPU solver with the same accuracy. More than 66% GPU memory can be saved and the computational speed can be accelerated up to 20 times. Hence, the computation time for Kelvin ship waves simulation can be significantly reduced by applying the GPU parallel numerical scheme, which lays a solid foundation for practical ocean engineering.

We study the tight focusing properties of generic bright and dark HyOPS beams in the presence of primary coma. The role of the polarization singularity index and handedness of the polarization of the HyOPS beams on the focused structure has been discussed. Results have been presented for the total intensity, component intensities, and component phase distributions for left-, and right-handed bright and dark star and lemon types singularities. The presence of primary coma distorted the focal plane intensity distributions for both positive and negative index generic C-points. Coma is known to disturb the circular symmetry of the focal plane intensity distribution. Similarly in tight focusing polarization is known to disturb the symmetry. Therefore, a beam with structured and inhomogeneous polarization distribution tightly focused under the influence of coma is a fit case to study. It is found that the presence of primary coma aberration in the focusing system produces a positional shift of the high-intensity peaks and a reduction of the intensity on one side of the center. As the strength of the primary coma increases, the focal plane intensity distributions shift more and more toward the right from the initial position. Unlike the scalar vortex case, in the case of HyOPS beams, the focal plane intensity distribution undergoes rotation, as the helicity of the HyOPS beams is inverted, in addition to shifting. All the component phase distributions are found to be embedded with phase vortices of charge ±1.

Ensuring precise electric locomotive stopping is pivotal for achieving intelligent and unmanned auxiliary transportation systems. Currently, human drivers' expertise plays a central role in the stopping process, posing challenges for automation and cargo location accuracy, particularly in coal mining. This paper focuses on achieving accurate electric locomotive stopping under diverse conditions by utilizing permanent magnet synchronous motor-driven locomotives. This technique introduces a novel stopping control method that combines a fuzzy Proportional-Integral-Derivative (F-PID) controller and a vector control model for permanent magnet synchronous motors (PMSM). After this, we design the F-PID using the PMSM technique based on the new fuzzy rules for each sub-system. In the end, extensive simulations and real-world experiments are conducted on an electric locomotive stopping test bed, which demonstrate the effectiveness of the proposed control method. From the simulation results, it confirming that the ability to achieve precise stopping under various working conditions.

Recently, Shehata et al. [37] introduced the r+1Rs(B;C; z) matrix function and established some properties. The aim of this study established to devote and derive certain basic properties including analytic properties, recurrence matrix relations, dierential properties, new integral representations, k- Beta transform, Laplace transform, fractional k-Fourier transform, fractional integral properties, the k-Riemann{Liouville and k-Weyl fractional integral and derivative operators an extended version of r+1Rs;k matrix function. We establish its relationships with other well known special matrix functions which have some particular cases in the context of three parametric Mittag-Leer matrix function, k- Konhauser and k-Laguerre matrix polynomials. Finally, some special cases of the established formulas are also discussed.

We show succinctly a natural origin of the beautiful Gamma function in light of the WJ distribution, linking to quite distinct fields.

This article is a study on the (k,s)-Riemann-Liouville fractional integral, a generalization of the Riemann-Liouville fractional integral. Firstly, we introduce several properties of the extended integral of continuous functions. Furthermore, we make the estimation of the Box dimension of the graph of continuous functions after the extended integral. It presents that the upper Box dimension of the (k,s)-Riemann-Liouville fractional integral for any continuous functions is no more than the upper Box dimension of the functions on the unit interval I=[0,1], which indicates that the upper Box dimension of the integrand f(x) will not be increased by the σ-order (k,s)-Riemann-Liouville fractional integral ksD−σf(x) where σ>0 on I. Additionally, we prove that the fractal dimension of ksD−σf(x) of one-dimensional continuous functions f(x) is still one.

This paper is devoted to the study of multi-dimensional integral transform with Fox H-function in the kernel in weighted spaces integrable functions in the domain R+n with positive coordinates. Mapping properties such as the boundedness, the range, the representations of the considered transformation are established.

In this paper for every number field K generated by a root α of a trinomial x⁷ + ax + b ∈ Z[x] and for every prime integer p, we calculate ν_p(i(K)), the highest power of p dividing the index i(K) of the field K. In particular, we calculate the index i(K). As application, when the index of K is not trivial, then K is not monogenic.

Body-centered cubic bismuth (Bi) is considered an enticing pressure marker, and therefore it is highly desirable to command its accurate equation of state (EOS). However, significant discrepancies are noted among the previous experimental EOSs. In the present work, the EOS up to 300 GPa is theoretically obtained by direct solving the partition function via a direct integral approach (DIA). The calculated results nearly reproduce the hydrostatic experimental measurements below 75 GPa, and the deviations from the measurements gradually get larger with increasing pressure. Based on the ensemble theory of equilibrium state, DIA works with high precision particularly in high-pressure conditions, so the hydrostatic EOS presented in this work is expected to be a reliable pressure standard.

This paper presents two control design strategies for voltage regulation in a single inductor dual output (SIDO) DC-DC Buck converter system. Based on a nominal multiple input multiple output (MIMO) plant model and performance requirements, both an LQR and a decoupled PI control laws are designed to control the power converter system, under parametric uncertainties such as voltage source variation, CPL power variation, and load resistance variation. Therefore, the SIDO converter can cascade operating, which may cause the undesired effects of voltage oscillation and reduce the system's stability margin. The control performance was assessed under both uncertainties resistive loading and power variation of a constant power load (CPL). A SIDO DC-DC Buck converter board has been developed for experimental tests. The experimental results show that the LQR strategy, compared to the PI strategy, presented a robust performance in the presence of parametric uncertainties in the resistive loads being, however, sensitive to uncertainties due to the presence of CPL loads.

The Polyurethane (PU) mixture is a high−performance engineering material with different dynamic properties compared to a hot−mixed asphalt mixture (HMA). This paper analyzed whether the Kramers−Kronig (K−K) relation and thermo−rheological simple properties were applicable to the PU mixture. Based on the result discussion, the PU mixture exhibited thermo−rheological simple properties within the test conditions. The time−temperature superposition principle (TTSP) was applicable for the PU mixture to construct the dynamic modulus master curve using the standard logistic Sigmoidal (SLS) model, the generalized logistic sigmoidal (GLS) model, and the Havriliak−Negami (HN) model. The Hilbert integral transformed SLS and GLS models for the phase angle can accurately fit the measured phase angle data with newly fitted shift factors and predict the phase angle within the viscoelastic range. The Core−Core and black space diagrams both displayed single continuous smooth curves, which can be utilized to characterize the viscoelastic property of the PU mixture. The K−K relation is applicable for the PU mixture to obtain the phase angle master curve model, storage modulus, and loss modulus from the complex modulus test results. The phase angle of the PU mixture cannot be derived from the dynamic modulus data and currently needs to be derived from the laboratory measurement.

Optical communications are described and analyzed by shaped beams; we challenge the ‎effects of parameters that impact the profile of a Flat Top Gaussian (FTG) beam. When ‎the laser beam propagates throughout the atmosphere, it can be influenced by different ‎optical phenomena including scattering, absorption, and turbulence due to changes in the ‎scintillation index and the forms of intensity that are displayed in the source and receiver ‎planes. In this project, the FTG laser beam that propagates through a weak turbulent ‎region is numerically investigated using open-source software. This simulation will be ‎performed according to a mathematical model based on the split-step beam propagation ‎method. Intensity distributions at the source plane and the received average intensity in ‎atmospheric turbulence are calculated, and additional contour is in the transducer plane. ‎The scintillation index, structure constant, source size, and other parameters, are applied ‎in the Rytov method to quantify the weak turbulent model. Moreover, these parameters ‎are analyzed in near-field propagation. Also, the effects of the beam’s scintillation and ‎beam wander are determined. All results simulated are discussed and compared with the ‎TEM00 Gaussian beam. Finally, these results are compared to measurements in the ‎experimental part of the work.

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.

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.