The Periodicity of the Accuracy of Numerical Integration Methods for the Solution of Different Engineering Problems

Newton-Cotes integration formulae have been researched for a long time, but the topic is still of interest since the correctness of the techniques has not yet been explicitly defined in a sequence for diverse engineering situations. The purpose of this paper is to give the readers an overview of the four numerical integration methods derived from Newton-Cotes formula, namely the Trapezoidal rule, Simpson's 1/3rd rule, Simpson's 3/8th rule, and Weddle's rule, as well as to demonstrate the periodicity of the most accurate methods for solving each engineering integral equation by varying the number of sub-divisions. The exact expressions by solving the numerical integral equations have been determined by Maple program and comparisons have been done using Python version 3.8.

CMMSE: Study of a new symmetric anomaly in the elliptic, hyperbolic and parabolic Keplerian motion

In the present work, we define a new anomaly, $\Psi$, termed semifocal anomaly. It is determined by the mean between the true anomaly, $f$, and the antifocal anomaly, $f^{\prime}$; Fukushima defined $f^{\prime}$ as the angle between the periapsis and the secondary around the empty focus. In this first part of the paper, we take an approach to the study of the semifocal anomaly in the hyperbolic motion and in the limit case correspoding to the parabolic movement. From here we find a relation beetween the semifocal anomaly and the true anomaly that holds independently of the movement type. We focus on the study of the two-body problem when this new anomaly is used as the temporal variable.\\ In the second part, we show the use of this anomaly —combined with numerical integration methods— to improve integration errors in one revolution. Finally, we analyze the errors committed in the integration process —depending on several values of the eccentricity— for the elliptic, parabolic and hyperbolic cases in the apsidal region.

Efficient computational modelling of smooth muscle orientation and function in the aorta

Understanding the mechanical effects of smooth muscle cell (SMC) contraction on the initiation and the propagation of cardiovascular diseases such as aortic dissection is critical. Framed by elastic lamellar sheets in the lamellar unit, there are SMCs in the media with a distinct radial tilt, which indicates their contribution to the radial strength. However, the mechanical effects of this type of anisotropy have not been fully discussed. Therefore, in this study, we propose a constitutive framework that models the passive and active mechanics of the aorta, taking into account the dispersed nature of the aortic constituents by applying the discrete fibre dispersion method. We suggest an isoparametric approach by evaluating various numerical integration methods and introducing a non-uniform discretization of the unit hemisphere to increase its computational efficiency. Finally, the constitutive parameters are fitted to layer-specific experimental data and initial computational results are briefly presented. The radial tilt of SMCs is also analysed, which has a noticeable influence on the mechanical behaviour of the aorta. In the absence of sufficient experimental data, the results indicate that the active contribution of SMCs has a remarkable impact on the mechanics of the healthy aorta.

Comparison of numerical integration methods for calculating diffraction of a plane electromagnetic wave by a rectangular aperture

We discuss features of the calculation of a Fraunhofer integral by traditional quadrature numerical integration methods and a special collocation Levin method when calculating the diffraction of a plane electromagnetic wave by a rectangular aperture. For the quadrature numerical integration methods, a criterion for the assessment of the integration step is derived depending on the screen size and required calculation accuracy. Advantages of the use of the special collocation Levin method in comparison with the traditional quadrature numerical integration methods are shown.

Numerical Integration Implementation Using Trapezoidal Rule Method To Calculate Aproximation Area Of West Java Province

An area can be shaped into a regular shape or an irregular shape. There is an area of irregular shape which is restricted by an unknown function, to determine that area must use a numerical integration. One of numerical integration methods is Trapezoidal Rule by replacing (????) with an integral approach function which can be evaluated, then let the (????) approximated by a linear polynomial in the certain interval, denoted as closed interval . This study is going to calculate the area of West Java Province by using this method with several different number of partitions in each quadrant such as, 9 partitions, 11 partitions, and 36 partitions in for different quadrants. This study provides the final result of the approximate area which will be compared with the actual area based in the error of result. The main finding is the approximate total area will be closer to the actual area followed by the increasing number of partitions.

Lumped, Constrained Cable Modeling With Explicit State-Space Formulation Using An Elastic Version of Baumgarte Stabilization

This paper presents a modeling approach for efficient simulation of slender structures, such as wires, cables and ropes. Lumped structural elements are connected using constraints. These are solved explicitly, using an elastic version of Baumgarte stabilization. This avoids singularities in the matrix inversions. The resulting explicit state-space formulation filters the higher order dynamics and can be solved using simple numerical integration methods. Constraints are demonstrated for modeling different aspects: Internal cable forces, one cable sliding along another cable and contact between cable and seabed. Also, a cable initialization routine is presented for rapid building of different interconnected cable geometries, ranging from cases in offshore crane operations to in-sea equipment such as seismic cables. Two case studies are presented to illustrate the effectiveness and the robustness of the proposed modeling approach; the first one being a test of two connected, sinking cables, and the last one being a larger case demonstrating the use of the cable library in an offshore seismic survey case.

Pair correlations of Halton and Niederreiter Sequences are not Poissonian

Niederreiter and Halton sequences are two prominent classes of higher-dimensional sequences which are widely used in practice for numerical integration methods because of their excellent distribution qualities. In this paper we show that these sequences—even though they are uniformly distributed—fail to satisfy the stronger property of Poissonian pair correlations. This extends already established results for one-dimensional sequences and confirms a conjecture of Larcher and Stockinger who hypothesized that the Halton sequences are not Poissonian. The proofs rely on a general tool which identifies a specific regularity of a sequence to be sufficient for not having Poissonian pair correlations.

Benchmarking of numerical integration methods for ODE models of biological systems

Ordinary differential equation (ODE) models are a key tool to understand complex mechanisms in systems biology. These models are studied using various approaches, including stability and bifurcation analysis, but most frequently by numerical simulations. The number of required simulations is often large, e.g., when unknown parameters need to be inferred. This renders efficient and reliable numerical integration methods essential. However, these methods depend on various hyperparameters, which strongly impact the ODE solution. Despite this, and although hundreds of published ODE models are freely available in public databases, a thorough study that quantifies the impact of hyperparameters on the ODE solver in terms of accuracy and computation time is still missing. In this manuscript, we investigate which choices of algorithms and hyperparameters are generally favorable when dealing with ODE models arising from biological processes. To ensure a representative evaluation, we considered 142 published models. Our study provides evidence that most ODEs in computational biology are stiff, and we give guidelines for the choice of algorithms and hyperparameters. We anticipate that our results will help researchers in systems biology to choose appropriate numerical methods when dealing with ODE models.

Adaptive Generalized Synchronization between Circuit and Computer Implementations of the Rössler System

The synchronization between chaotic systems implemented in similar ways—e.g., computer models or circuits—is a well-investigated topic. Nevertheless, in many practical applications, such as communication, identification, machine sensing, etc., synchronization between chaotic systems of different implementation types—e.g., between an analog circuit and computer model—might produce fruitful results. In this research, we study the synchronization between a circuit modeling the Rössler chaotic system and a computer model using the same system. The theoretical possibility of this kind of synchronization is proved, and experimental evidence of this phenomenon is given with special attention paid to the numerical methods for computer model simulation. We show that synchronization between a circuit with uncertain parameters and a computer model is possible, and the parameters obtained from the synchronized computer model are in high correspondence with the circuit element specification. The obtained results establish the possibility of using adaptive generalized synchronization for the parameter identification of real systems. It was also found that Heun’s method yielded the most accurate results in synchronization among second-order numerical integration methods. The best among the first-order methods appears to be the Euler–Cromer method, which can be of interest in embedded applications.