MATLAB Code Research Articles

Chaotic hidden attractor in a fractional order system modeling the interaction between dark matter and dark energy

In this paper the dynamics of a fractional order system modeling the interaction between dark matter and dark energy is analytically and numerically studied. It is shown for the first time that systems modeling the interaction between dark matter and dark energy, chaotic hidden attractors can be present. The chaotic attractor coexists with two asymptotically stable equilibria. Equilibria of the linearized system exhibit a center-like behavior. The numerical integration is done by means of the Adams–Bashforth–Moulton scheme and the finite Lyapunov exponents are numerically determined with a dedicated Matlab code. The 3D representation of the chaotic hidden attractor reveals the fact it is not connected with the equilibria, being “hidden” somewhere in the considered phase space.

An Efficient Code for the Multi-Material Topology Optimization of 2D/3D Continuum Structures Written in Matlab

This paper presents an efficient and compact Matlab code for 2D and 3D topology optimization of multi-materials. We extend a classical 88-line-based educational code to the multi-material problem using the mapping-based interpolation function, with which each material is modeled equally and can obtain a clear 0, 1 result for each material of the optimized structures that can be manufactured easily. A solution of topology optimization of multi-materials with minimum compliance under volume constraints is documented, including the details of the optimization model, filtering, projection, and sensitivity analysis procedures. Several numerical examples are also conducted to demonstrate the effectiveness of the code, and the influence of the parameter on the optimized results is also analyzed. Complete 2D and 3D Matlab codes are provided.

Image-based gradation and aggregate characterisation: case of cement-stabilised quarry fines

Quarry fines are secondary materials that can be reused for road construction, commonly after improving their properties through stabilising. To better understand and characterise these recycled materials, we created an image-based, two-dimensional representation of the particles that compose the material. This computational representation enables (1) the quantification of the (area-based) gradation from images, (2) a detailed particle-by-particle characterisation and (3) potentially, random rearrangement of the particles into reconstructed specimens for further computational analyses. The images of cement-stabilised quarry fines were obtained using a scanning electron microscope. The particles were isolated and characterised using image processes and a computational approach (MATLAB code). The techniques and code presented herein can be applied as well to other granular-based materials. The MATLAB code and associated files are available as an open resource in GitHub (https://github.com/danielcastillo-beta/particle-analysis).

Signal Processing for the Condition-Based Maintenance of Rotating Machines via Vibration Analysis: A Tutorial.

One of the common methods for implementing the condition-based maintenance of rotating machinery is vibration analysis. This tutorial describes some of the important signal processing methods existing in the field, which are based on a profound understanding of the component's physical behavior. Furthermore, this tutorial provides Python and MATLAB code examples to demonstrate these methods alongside explanatory videos. The goal of this article is to serve as a practical tutorial, enabling interested individuals with a background in signal processing to quickly learn the important principles of condition-based maintenance of rotating machinery using vibration analysis.

Simulink Methods to Simulate COVID-19 Outbreak Using Fractional Order Models

The objective of this work is to solve a coronavirus transmission model using simulink, a platform for model based design that facilitates simulation and design at the system level. The simulation is divided into two sections: the first, deals with setting up the parameters, and the second deals with computing the fractional derivatives. The fundamental SIR, SEIR, SEIQR, and SEIARM models were used in this study. An effective, quick, easy, and visually appealing method is used to simulate pandemic outbreaks. To accurately follow the progress of the infection, we contrast the simulation findings with those acquired by MATLAB code. The applications can be used for research projects and also as a teaching tool.

One-dimensional pump geometry prediction modeling for energy loss analysis of pumps working as turbines

Pumps as turbine (PAT) for small-scale hydropower, pressure relief valves, and hydro-pumping applications have recently gained popularity. In contrast to conventional turbines, pumps do not have movable guiding vanes to accommodate changes in flow rate and head. This limits their use to the single best operating point. As a result, the precise estimation of reverse mode operating parameters is critical. One of the most successful techniques in this regard is the energy loss model. However, it necessitates extremely detailed pump geometry and experimental data. Although it is claimed to be more accurate, the lack of data makes it impractical. This calls for the devising of a pump geometrical parameter prediction model from limited manufacturer information. Using 137 sets of commercial pump operating data, geometry redesign was carried out by MATLAB code. Following that, multivariate linear regression analysis was carried out, and a one-dimensional geometrical prediction model was developed. Then the result was validated with computational analysis and measurement. The model results were 99 %, 100 %, and 95 % correlated for pump inlet, volute, and throat inlet diameters, respectively, with the experimental results. Similarly, the model results were also 99 %, 71 %, and 89 % correlated for pump inlet, outlet, and volute blade width, respectively, with experimental results. The model gave a more accurate geometrical result than the analytical and numerical methods.

MFCDFT and impedance characteristic-based adaptive technique for fault and power swing discrimination

ABSTRACT In today’s fast-changing power systems, it is vital to have a protection system that can distinguish between power swing and fault conditions. With advanced protective devices (Digital Relays), an appropriate protection strategy can be established to identify power swing and fault conditions individually. The proposed algorithm uses the sliding window concept, which requires very little computational data unless any abnormal conditions exist in the power system. Modified Full Cycle Discrete Fourier Transform (MFCDFT) is more effective in the extraction of the fundamental phasors from the original complex signal by effectively removing the harmonics and DC components. An adaptive impedance characteristic-based MFCDFT is utilized for the discrimination of power swing and fault in this work. Several fault scenarios and power swings are generated on the test system concurrently and independently from the relay station at varied distances. PSCAD/EMTDC is utilized on a single-machine infinite bus system to analyze various test conditions. The simulation is utilized in this case to collect data from various test circumstances. After obtaining the data, the MFCDFT algorithm along with adaptive impedance estimation is performed using MATLAB code. After examining various test results, it is discovered that the suggested technique correctly recognizes power swings and various fault scenarios. It can also identify the difference between stable and unstable power swings. According to the results, the proposed technique is proficient at spotting fault and power swing conditions. In addition, the algorithm ensures dependable operation with quick detection and avoids maloperation under a typical swing condition.

Heat and mass transfer under MHD mixed convection in a four‐sided lid‐driven square cavity

AbstractThis paper investigates the heat and mass transfer under magnetohydrodynamic mixed convection flow of a binary gas mixture in a four‐sided lid‐driven square cavity. The enclosure's left wall is sinusoidally heated and acts as a source term, while the right wall functions as a sink. The cavity's horizontal walls are adiabatic and impermeable to mass transfer. The governing equations under Boussinesq approximation and stream function‐vorticity formulation are solved using the alternating‐direction‐implicit scheme, a finite‐difference method. The numerical scheme's consistency and stability are demonstrated using the matrix method. The MATLAB code is written, validated against some existing studies, and used to perform numerical simulations. The numerical solutions are graphically examined by visualizing the streamline, isotherm, and concentration contours for nondimensional parameters, such as Hartmann number , heat absorption or generation coefficient , Richardson number , and buoyancy ratio . The magnetic field modifies the temperature and concentration distribution in the cavity, depending on the convection mode. The magnetic field forces the fluid to stagnate in different regions of the cavity, depending on the mode of convection. It was found that the difference between the maximum and minimum temperature and concentration at the cavity's midpoint increases up to 13 and 10 times, respectively, in the natural convection compared with the forced convection. The average Nusselt number on the vertical walls of the cavity is maximum in natural convection in the absence of a magnetic field but reaches a minimum value at in forced and mixed convection. The average Sherwood number on the cavity's vertical walls decreases with the magnetic field in mixed and natural convection.

Design and Position Control of a Robot with 5 Degrees of Freedom

In recent decades, robotics and artificial intelligence have gained significant importance for their involvement in various industrial processes currently at their peak. Nevertheless, advanced robots are not designed explicitly for lemon supply. This research aims to develop a manipulator robot with 5 degrees of freedom and control its trajectories for lemon supply purposes. To achieve this goal, kinematic and dynamic calculations of the manipulator robot were performed, along with the development of programming code in Matlab to determine its trajectories, positions, speeds, and accelerations. In addition, the Proportional Integral Derivative (PID) tuner was used to obtain the optimal controller parameters and ensure accurate joint trajectory generation.

Cyclostationarity of Communication Signals in Underwater Acoustic Channels

The effect of underwater acoustic propagation on the cyclostationary features of communication signals is modeled and analyzed. Two kinds of channels are considered: the multiscale-multilag channel, over which mobile and wideband acoustic systems usually communicate, and the dispersive channel resulting from low-frequency modal propagation in shallow water. It is shown that multiscale-multilag channels transform cyclostationary signals into a sum of velocity and acceleration-dependent time-warped cyclostationary processes. This time-warping is carefully taken into account to efficiently recover the cyclostationary features. On the other hand, it is found that low-frequency dispersive channels preserve the original periodicity but attenuate the shorter cycles and spread the correlations. To illustrate the theoretical results, applications with simulated and real data are also presented. Specifically, the problem of estimating time-varying Doppler scales is addressed for multiscale-multilag channels as well as the detection of signals with unique cyclostationary signatures. The example of blind symbol-rate estimation applied to covert communications in dispersive channels is also discussed. Special attention is paid to phase-shift keying (PSK), quadrature amplitude modulation (QAM), orthogonal frequency-division multiplexing (OFDM), and direct-sequence spread-spectrum (DSSS) signals. Accompanying supplementary material provides the MATLAB code used for the estimation and detection examples.

Experimental study and modeling of water film thickness in aero-engine under air–water mist flow

The water film not only plays an important role in the mass, momentum, and energy transfer between the air–water-surface but also determines the on-wing washing effect of the aero-engine. In view of this, air–water mist flow visualization experiments have been conducted at different gas velocities and water-to-air ratios in a compressor cascade, and the microscopic water film images have been analyzed to extract the transient water film thickness data by the Matlab code. It was found that the transient water film thickness fluctuation has no obvious association with the gas velocity, and the water film fluctuation is more affected by the water-to-air ratio. As the water-to-air ratio increases, the fluctuation magnitude of the water film thickness increases. The average water film thickness has been studied in relation to gas velocity and water flow rate, i.e., the average thickness of water film decreases with increasing gas velocity and increase with the increasing water flow rate. On the basis of the water film flow equation and taking the droplet collection efficiency into account, i.e., from the perspective of the physical mechanism of water film formed, a new model for predicting the water film thickness of a compressor blade surface under the air–water mist flow condition was proposed and validated. This model predicted, with a root mean square error and the mean absolute percentage error of 11.6% and 9.15%, respectively, under the present experimental flow conditions.

Establishing a Mouse Genotyping Core Facility Based on Automation, High-resolution Melting Analysis, and Purpose-developed MATLAB Applications for Sample Tracking, Data Analysis, and Report Generation

An in-house genotyping facility should aim to be more cost-effective than outsourced service and more reliable than genotyping performed by short-term employees or students of individual research groups. Reliable genotyping allows efficient and economical management of mice colonies and promotes accurate and reproducible research results. Here we provide a detailed description of our approach to establishing a genotyping core facility, relying on automated PCR assembly and high-resolution melting (HRM) analysis (first derivative). The workflow we devised was tightly managed by purpose-designed applications developed using MATLAB App Designer that allowed straightforward work planning, ensured sample tracking throughout the process, and provided a platform for reliable data analysis and generation of genotyping reports. We successfully transitioned PCR product analysis of more than 250 different target genes from standard gel electrophoresis to the more advanced HRM analysis. About 23% of the target genes required a redesign of primers to adapt to our protocol. The process was highly universal, and only 2% of the target genes required deviation from the standard PCR method to a more restricted protocol that reduces the amplification of nonspecific products. We currently run more than 1,000 PCR reactions weekly, of samples taken at weaning or experimental endpoint, and assemble a large variety of target genes in every PCR plate. We also showed that genotyping of blastocytes instead of embryos can serve as quality control of cryopreservation. Thus, our genotyping protocol promotes the 3Rs (Replacement, Reduction, and Refinement) principles. Our refined genotyping process facilitates cost-effective colony management, replaces tissue types as well as traditional methods with advanced ones, and provides reliable results in a timely manner. MATLAB codes and related data are available in supplementary materials and online.

An exploratory study of “large‐n” detrital zircon geochronology of the Book Cliffs, UT via rapid (3 s/analysis) U–Pb dating

AbstractDetrital zircon (DZ) U–Pb geochronology has improved the way geologists approach questions of sediment provenance and stratigraphic age. However, there is debate about what constitutes an appropriate sample size (i.e., the number of dates in a DZ sample, n), which depends on project objectives, sample complexity, and, critically, analytical budget. Additionally, there is ongoing concern about bias introduced by zircon grain size. We tested a recently developed rapid (3 s/analysis) data acquisition method by multicollector laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) that incorporates an automated selection routine and calculates two‐dimensional grain geometry from polished sample surfaces. Eleven samples were analysed from below and above the Late Cretaceous (Campanian) basal Castlegate unconformity of the Book Cliffs, Utah, in a down‐depositional‐dip transect including Price, Horse, Tusher, and Thompson canyons. 12,448 new concordant dates were generated during two measurement sessions. Results are consistent with recent studies suggesting there is no major provenance change and little time (1–2 Myr) represented across the unconformity. Grain size and sample size both exert a strong control on sample dissimilarity. Age distributions constructed from subsamples of large grains are systematically less similar to whole samples; age distributions composed of small grains are overall more similar to whole samples. As such, North American sediment sources that produce large grains such as the Grenville and Yavapi‐Mazatzal belts can bias age distributions if only large grains are analysed. A sample size of n = 100 is inadequate for characterizing age distributions as complex as those of the Book Cliffs, whereas a sample size of n = 300 provides good characterization. Sample size of n ≈ 1000 or more is unnecessary unless project objectives include scanning for subordinate age groups, such as when identifying the youngest grains for calculating a maximum depositional age (MDA). Dates used in MDA calculations acquired with rapid acquisition are best re‐analysed with longer LA‐ICP‐MS acquisition methods or isotope dilution thermal ionization mass spectrometry for increased accuracy and precision. We include new MATLAB code and open‐source software programs, DZpick and DZmda, for automated spot picking and calculating MDAs.

Performance evaluation of gamma-type Stirling engine using combined Schmidt and mechanical loss model

This paper focuses on the study of gamma-type Stirling engine prototype using a combined Schmidt closed-form and mechanical loss analysis. Not restricted to optimizing the indicated power as classic Schmidt theory is set to, this analysis allows to maximize the shaft power due to the mechanical loss in power transmission. For this purpose, MATLAB code was developed to calculate the indicated and the shaft powers of the engine at different operating parameters. The results showed that shaft power peaks at swept volume ratios smaller than those of indicated power at different values of mechanism effectiveness. Within the range of engine mechanism effectiveness typically between 0.7 and 0.9, it was found that maximum shaft power for this particular engine can be achieved at different optimum values of swept volume ratio between 0.75 and 0.95 and phase angle between 80° and 90°. However, an optimum swept volume ratio was found to be <i>k</i>=0.55 of the same engine size for different scenarios of operation. Also, the developed model can be used as a design tool in the preliminary stage to find the optimum geometry of the engine. The new engine design parameters including the stroke, the crank radius and power piston bore, and engine alteration were presented.

Tool path planning for hole-flanging process using single point incremental forming

Hole-flanging is a process used to create flanges around holes in sheet metal parts. It traditionally relies on dies and punches, posing cost and time challenges for small quantities or prototypes. While Incremental forming does not require dedicated sets of tools and dies, and it provides more flexibility and adaptability with enhanced formability, and emerged as a better alternative to conventional press forming. However, precise and efficient tool path planning is essential to achieve high-quality hole-flanged components. Precise toolpath planning is one of the crucial and important part in designing of hole flanging process for both through numerical simulations and experiments. The objective of the present work is to generate tool paths as well as associated G-codes using MATLAB code for different geometries of incremental hole flanging process using SPIF. Six different geometries of hole flanges were considered to generate tool path using MATLAB code. The proposed methodology was successfully applied for different shapes of hole flanges. Since the proposed methodology provides unified programming technique in the current MATLAB environment, which performs dual duties, creating toolpath coordinates against time data and outputting a G-code file and it saves computational time and programming efforts. Hence, proposed methodology is found to be better for incremental hole flanging in contrast to previously existed methodology for development of toolpath for incremental forming cases. Thus, it was concluded that this proposed methodology can be effective for capturing any intricate shapes for incremental hole flanging. You should leave 8 mm of space above the abstract and 10 mm after the abstract. The heading Abstract should be typed in bold 9-point Arial. The body of the abstract should be typed in normal 9-point Times in a single paragraph, immediately following the heading. The text should be set to 1 line spacing. The abstract should be centred across the page, indented 17 mm from the left and right page margins and justified. It should not normally exceed 200 words.

Establishing a Mouse Genotyping Core Facility Based on Automation, High-resolution Melting Analysis, and Purpose-developed MATLAB Applications for Sample Tracking, Data Analysis, and Report Generation.

An in-house genotyping facility should aim to be more cost-effective than outsourced service and more reliable than genotyping performed by short-term employees or students of individual research groups. Reliable genotyping allows efficient and economical management of mice colonies and promotes accurate and reproducible research results. Here we provide a detailed description of our approach to establishing a genotyping core facility, relying on automated PCR assembly and high-resolution melting (HRM) analysis (first derivative). The workflow we devised was tightly managed by purpose-designed applications developed using MATLAB App Designer that allowed straightforward work planning, ensured sample tracking throughout the process, and provided a platform for reliable data analysis and generation of genotyping reports. We successfully transitioned PCR product analysis of more than 250 different target genes from standard gel electrophoresis to the more advanced HRM analysis. About 23% of the target genes required a redesign of primers to adapt to our protocol. The process was highly universal, and only 2% of the target genes required deviation from the standard PCR method to a more restricted protocol that reduces the amplification of nonspecific products. We currently run more than 1,000 PCR reactions weekly, of samples taken at weaning or experimental endpoint, and assemble a large variety of target genes in every PCR plate. We also showed that genotyping of blastocytes instead of embryos can serve as quality control of cryopreservation. Thus, our genotyping protocol promotes the 3Rs (Replacement, Reduction, and Refinement) principles. Our refined genotyping process facilitates cost-effective colony management, replaces tissue types as well as traditional methods with advanced ones, and provides reliable results in a timely manner. MATLAB codes and related data are available in supplementary materials and online.

The dropper stress characteristics under different moving load speeds for a high-speed railway

Dropper failure seriously threatens the operation safety of a high-speed railway. In this work, for a simple chain suspension catenary, one span with five droppers is performed to establish a model and thus the effects of the moving load speed on dropper stress are investigated. First, the partial differential vibration equation of dropper is obtained through the mechanical analysis and converted into the finite difference equation. Then, we consider contact line as a beam element to obtain its motion equation. Furthermore, the boundary and initial conditions of five droppers are determined. Finally, the stresses of five droppers are numerically calculated and the effects of the moving load speed on dropper stress are investigated by writing a MATLAB code. The results suggest that the dropper location significantly affects its stress. Compared with other droppers, droppers II and IV have much more severe vibration amplitudes. Different moving load speeds could cause different stress change of each dropper. With the increasing speed, dropper experiences longer bending compression stage and the bending amplitude increases. The impact of the moving load speed on dropper stress is significant.

An iterative algorithm to estimate the energy spectrum of an electron beam from PDD curves

Electron beam central-axis percentage depth dose (PDD) curves in water phantom are routinely employed to evaluate the electron beam energy at the phantom surface, in particular the mean and most probable energies from the values of R 50 (half-value range) and Rp (practical range). However, these two quantities are not enough to evaluate important details of the energy distribution, such as the FWHM (Full Width Half Maximum) and the possible presence of a low-energy tail. This paper presents a numerical method that allows estimating the shape of the energy spectrum from a PDD curve. The algorithm uses a database consisting of a set of depth dose curves for monochromatic beams computed by FLUKA in the range 0.1-7.0 MeV by steps of 0.1 MeV and, using an adaptive iterative Monte Carlo process, reconstructs the incident energy spectrum by minimizing the distance between the measured PDD and the computed one. Applications of a MATLAB code based on this algorithm to simulated and real measurements of electron beams done at APAM lab (ENEA Frascati) are presented. This approach represents a strong simplification with respect to energy analysis based on the use of a magnetic spectrometer.

INVESTIGATING ENERGY ABSORPTION CAPABILITY OF ADDITIVELY MANUFACTURED CUBOCTAHEDRAL LATTICE STRUCTURES VIA TAGUCHI’S METHOD: EFFECTS OF PROCESS PARAMETERS

This research aims to investigate the influence of process parameters of fused filament deposition additive manufacturing technique on the energy absorption capacity of cuboctahedral lattice structures using Taguchi’s method. Four process parameters (i.e., print temperature, print speed, layer thickness, and line width) were considered with three levels of parameter values for each, which resulted in nine combinations in Taguchi’s L9 orthogonal array. The lattice structure was fabricated with each of nine combinations of process parameters and tested under a compression load to obtain the energy absorption experimentally. Signal-to-noise ratio analysis was conducted for the results obtained from the experiments of the nine Taguchi sets. Two more lattice specimens were fabricated with the parameter values which resulted in the best and worst energy absorption results and tested. The optimum values of the print temperature, print speed, layer thickness, and line width were determined to be 225 ºC, 30 mm/s, 0.12 mm, and 0.35 mm, respectively. The specific energy absorption (SEA) of the lattice specimen fabricated by using the optimum process parameter values presented 9.4% improvement compared to the highest SEA obtained from the lattices in the Taguchi’s L9 array.

Possibility of Quenching of Limit Cycles in Multi Variable Nonlinear Systems with Special Attention to 3X3 Systems

The present work proposes novel methods of Quenching self-sustained oscillations in the event of the existence of limit cycles (LC) in 3x3 non-linear systems. It explores the possibility of Stabilising/Quenching the LC by way of signal stabilization using high frequency dither signals both deterministic and random when 3X3 systems exhibit such self-sustained nonlinear oscillations under autonomous state. The present work also explores the suppression limit cycles of 3X3 systems using state feedback by either arbitrary pole placement or optimal selection of pole placement. The complexity involved, in implicit non-memory type nonlinearity for memory type nonlinearities, it is extremely difficult to formulate the problem. Under this circumstance, the harmonic linearization/harmonic balance reduces the complexity considerably. Furthermore, the method is made simpler assuming the whole 3X3 system exhibits the LC predominantly at a single frequency. It is equally a formidable task to make an attempt to suppress the limit cycles for 3X3 systems with memory type nonlinearity in particular. Backlash is one of the nonlinearities commonly occurring in physical systems that limit the performance of speed and position control in robotics, the automation industry, and other occasions of modern applications. The proposed methods are well illustrated through examples and substantiated by digital simulation (a program developed using MATLAB CODES) and the use of the SIMULINK Toolbox of MATLAB software.