scholarly journals Numerical Counterexamples of Lorenz System in Implicit Time Scheme

Author(s):  
Xuan Chen

In nonlinear self-consistent system, Lorenz system (Lorenz equations) is a classic case with chaos solutions which are sensitively dependent on the initial conditions. As it is difficult to get the analytical solution, the numerical methods and qualitative analytical methods are widely used in many studies. In these papers, Runge-Kutta method is the one most often used to solve these differential equations. However, this method is still a method based on explicit time scheme, which would be the main reason for the chaotic solutions to Lorenz system. In this work, numerical experiments based on implicit time scheme and explicit scheme are setup for comparison, the results show that: in implicit time scheme, the numerical solutions (counterexamples) are without chaos; for an original volume, the volume shrinks exponentially fast to 0 in common.

2021 ◽  
Vol 7 (2) ◽  
pp. 97-101
Author(s):  
H. P. Kandel ◽  
J. Kafle ◽  
L. P. Bagale

Many physical problems, such as heat transfer and wave transfer, are modeled in the real world using partial differential equations (PDEs). When the domain of such modeled problems is irregular in shape, computing analytic solution becomes difficult, if not impossible. In such a case, numerical methods can be used to compute the solution of such PDEs. The Finite difference method (FDM) is one of the numerical methods used to compute the solutions of PDEs by discretizing the domain into a finite number of regions. We used FDMs to compute the numerical solutions of the one dimensional heat equation with different position initial conditions and multiple initial conditions. Blacksmiths fashioned different metals into the desired shape by heating the objects with different temperatures and at different position. The numerical technique applied here can be used to solve heat equations observed in the field of science and engineering.


2018 ◽  
Vol 13 (5) ◽  
pp. 47 ◽  
Author(s):  
Fatima Z. Ahmed ◽  
Mayada G. Mohammed ◽  
Dmitry V. Strunin ◽  
Duc Ngo-Cong

We present numerical solutions of the semi-empirical model of self-propagating fluid pulses (auto-pulses) through the channel simulating an artificial artery. The key mechanism behind the model is the active motion of the walls in line with the earlier model of Roberts. Our model is autonomous, nonlinear and is based on the partial differential equation describing the displacement of the wall in time and along the channel. A theoretical plane configuration is adopted for the walls at rest. For solving the equation we used the One-dimensional Integrated Radial Basis Function Network (1D-IRBFN) method. We demonstrated that different initial conditions always lead to the settling of pulse trains where an individual pulse has certain speed and amplitude controlled by the governing equation. A variety of pulse solutions is obtained using homogeneous and periodic boundary conditions. The dynamics of one, two, and three pulses per period are explored. The fluid mass flux due to the pulses is calculated.


2015 ◽  
Vol 07 (01n02) ◽  
pp. 1550004 ◽  
Author(s):  
Gregori J. Clarke ◽  
Samuel S. P. Shen

This study uses the Hilbert–Huang transform (HHT), a signal analysis method for nonlinear and non-stationary processes, to separate signals of varying frequencies in a nonlinear system governed by the Lorenz equations. Similar to the Fourier series expansion, HHT decomposes a data time series into a sum of intrinsic mode functions (IMFs) using empirical mode decomposition (EMD). Unlike an infinite number of Fourier series terms, the EMD always yields a finite number of IMFs, whose sum is equal to the original time series exactly. Using the HHT approach, the properties of the Lorenz attractor are interpreted in a time–frequency frame. This frame shows that: (i) the attractor is symmetric for [Formula: see text] (i.e. invariant for [Formula: see text]), even though the signs on [Formula: see text] and [Formula: see text] are changed; (ii) the attractor is sensitive to initial conditions even by a small perturbation, measured by the divergence of the trajectories over time; (iii) the Lorenz system goes through “windows” of chaos and periodicity; and (iv) at times, a system can be both chaotic and periodic for a given [Formula: see text] value. IMFs are a finite collection of decomposed quasi-periodic signals, starting from the highest to lowest frequencies, providing detection of the lower frequency signals that may have otherwise been “hidden” by their higher frequency counterparts. EMD decomposes the original signal into a family of distinct IMF signals, the Hilbert spectra are a “family portrait” of time–frequency–amplitude interplay of all IMF members. Together with viewing the IMF energy, it is easy to discern where each IMF resides in the spectra in relation to one another. In this study, the majority of high amplitude signals appear at low frequencies, approximately 0.5–1.5. Although our numerical experiments are limited to only two specific cases, our HHT analyses of time–frequency, marginal spectra, and energy and quasi-periodicity of each IMF provide a novel approach to exploring the profound and phenomena-rich Lorenz system.


2021 ◽  
Vol 19 (1) ◽  
pp. 351-370
Author(s):  
Zulqurnain Sabir ◽  
◽  
Muhammad Asif Zahoor Raja ◽  
Abeer S. Alnahdi ◽  
Mdi Begum Jeelani ◽  
...  

<abstract> <p>These investigations are to find the numerical solutions of the nonlinear smoke model to exploit a stochastic framework called gudermannian neural works (GNNs) along with the optimization procedures of global/local search terminologies based genetic algorithm (GA) and interior-point algorithm (IPA), i.e., GNNs-GA-IPA. The nonlinear smoke system depends upon four groups, temporary smokers, potential smokers, permanent smokers and smokers. In order to solve the model, the design of fitness function is presented based on the differential system and the initial conditions of the nonlinear smoke system. To check the correctness of the GNNs-GA-IPA, the obtained results are compared with the Runge-Kutta method. The plots of the weight vectors, absolute error and comparison of the results are provided for each group of the nonlinear smoke model. Furthermore, statistical performances are provided using the single and multiple trial to authenticate the stability and reliability of the GNNs-GA-IPA for solving the nonlinear smoke system.</p> </abstract>


2014 ◽  
Vol 746 ◽  
pp. 332-367 ◽  
Author(s):  
Lawrence K. Forbes

AbstractA previous study by Forbes (ANZIAM J., vol. 53, 2011, pp. 87–121) has argued that, when a light fluid is injected from a point source into a heavier ambient fluid, the interface between them is most unstable to perturbations at the lowest spherical mode. This means that, regardless of initial conditions, the outflow from a point source eventually becomes a one-sided jet. However, two-sided (bi-polar) outflows are nevertheless often observed in astrophysics, in apparent contradiction to this prediction. While there are many possible explanations for this fact, the present paper considers the effect of a straining flow in the ambient fluid. In addition, solid-body rotation in the inner fluid is also accounted for, in a Boussinesq viscous model. It is shown analytically that there are circumstances under which straining flow alone is sufficient to convert the one-sided jet into a genuine bi-polar outflow, in linearized theory. This is confirmed in a numerical solution of a viscous model of the flow, based on a spectral solution technique that accounts for nonlinear effects. Rotation can also generate flows that are two-sided, and this is likewise revealed through an asymptotic analysis and numerical solutions of the nonlinear equations.


2017 ◽  
Vol 13 (2) ◽  
pp. 4657-4670
Author(s):  
W. S. Amer

This work touches two important cases for the motion of a pendulum called Sub and Ultra-harmonic cases. The small parameter method is used to obtain the approximate analytic periodic solutions of the equation of motion when the pivot point of the pendulum moves in an elliptic path. Moreover, the fourth order Runge-Kutta method is used to investigate the numerical solutions of the considered model. The comparison between both the analytical solution and the numerical ones shows high consistency between them.


2012 ◽  
Vol 26 (29) ◽  
pp. 1250146 ◽  
Author(s):  
BHASKAR SEN GUPTA ◽  
SHANKAR P. DAS

The renormalized dynamics described by the equations of nonlinear fluctuating hydrodynamics (NFH) treated at one loop order gives rise to the basic model of the mode coupling theory (MCT). We investigate here by analyzing the density correlation function, a crucial prediction of ideal MCT, namely the validity of the multi step relaxation scenario. The equilibrium density correlation function is calculated here from the direct solutions of NFH equations for a hard sphere system. We make first detailed investigation for the robustness of the correlation functions obtained from the numerical solutions by varying the size of the grid. For an optimum choice of grid size we analyze the decay of the density correlation function to identify the multi-step relaxation process. Weak signatures of two step power law relaxation is seen with exponents which do not match predictions from the one loop MCT. For the final relaxation stretched exponential (KWW) behavior is seen and the relaxation time grows with increase of density. But apparent power law divergences indicate a critical packing fraction much higher than the corresponding MCT predictions for a hard sphere fluid.


1999 ◽  
Vol 31 (4) ◽  
pp. 217-230 ◽  
Author(s):  
R. A. Lebensohn ◽  
T. Leffers

The rules for the lattice rotation during rolling (plane strain) deformation of fcc polycrystals are studied with a viscoplastic self-consistent model. Very high values of the ratesensitivity exponent are used in order to establish Sachs-type conditions with large local deviations from the macroscopic strain. The lattice rotation depends on the grain shape. For equiaxed grains the lattice rotation follows the MA rule, which is the one normally used in solid mechanics. For elongated and flat grains the lattice rotation follows a different rule, the PSA rule. In the standard version the model performs a transition from MA to PSA with increasing strain. There is avery clear difference between the textures resulting from the two different rules. MA leads to a copper-type texture, and PSA leads to a brass-type texture.


1996 ◽  
Vol 324 ◽  
pp. 163-179 ◽  
Author(s):  
A. Levy ◽  
G. Ben-Dor ◽  
S. Sorek

The governing equations of the flow field which is obtained when a thermoelastic rigid porous medium is struck head-one by a shock wave are developed using the multiphase approach. The one-dimensional version of these equations is solved numerically using a TVD-based numerical code. The numerical predictions are compared to experimental results and good to excellent agreements are obtained for different porous materials and a wide range of initial conditions.


Sign in / Sign up

Export Citation Format

Share Document