On CML model for study of spatiotemporal chaos

A new coupled map lattic(CML)model is given by using some stability analysis for the related difference equations.Numerical results show that the new model is.an effective one of studying spatiotemporal chaos,especially for strongly coupled systems.     

Periodic solutions in one-dimensional coupled map lattices

IntroductionDynamicsofcoupledmaplattices (CML )haverecentlyattractedconsiderableattention[1- 4].IthasbeenfoundthatCMLsexhibitavarietyofspace_timepatterns:space_timeintermittency ,spatiotemporalchaos,frontsandinterfaces,nonlinearperiodicsolutions.Inthispaper,anone_dimensionalinfinitechainofnonlinearcoupledmapsisconsideredandtheexistenceofperiodicsolutionsisprovedwhenthecouplingstrengthαissmallenough .Ourapproachistheanti_integrablelimitintroducedinRefs.[5-7] .Intheanti_integrablelimit(coup…     

二维Logistic映射的分岔与分形

理论分析了二维Logistic映射的分岔,并采用相图、分岔图、功率谱、Lyapunov指数和分维数计算的方法,揭示出：二维Logistic映射可按倍周期分岔和Hopf分岔走向混沌;在倍周期分岔过程中,系统在参数空间和相空间中都表现出自相似性和尺度变换下的不变性．对二维Logistic映射的吸引盆及其Mandelbrot-Julia集(简称M-J集)的研究表明：吸引盆中周期和非周期区域之间的边界是分形的,这意味着无法预测相平面上点运动的归宿;M-J集的结构由控制参数决定,且它们的边界是分形的．     

Space temporal maps for vortical flow field construction

A class of models to generate a two-dimensional complex flow field, deriving from the coupled map lattice dynamics, is presented here. It automatically satisfies continuity equation for an incompressible fluid. The method is numerically implemented on a square lattice and some results relatively to a fully deterministic and a semi stochastic evolution are presented here. The qualitative similarity with two dimensional hydrodynamic turbulence is encouraging. In view of these first results, directions for advances are proposed.

---

Sommario Si presenta una classe di modelli per generare un campo di moto complesso, derivante dai modelli dinamici a mappe accoppiate, in modo tale che l'equazione di continuita' per un fluido incomprimibile sia automaticamente soddisfatta. Il metodo e'implementato numericamente su una griglia quadrata. Si presentano alcuni risultati relativi ad un'evoluzione deterministica e semi stocastica. La somiglianza qualitativa con campi di moto idrodinamici turbolenti bidimensionali e'incoraggiante. Alla luce dei risultati ottenuti si propongono inoltre possibili direzioni per ulteriori sviluppi.

     

Large-Scale Flows in a Turbulent Convective Layer with an Immersed Moving Thermal Insulator

Large-scale motions developing against the background of developed convective turbulence in a rectangular vessel heated from below, in which a horizontal disk floats freely at a fixed depth, are studied experimentally. It is found that regular or chaotic oscillations may develop in the vessel. Transitional regimes characterized by interruptions of the motion or frequency variations are also possible. A regime map in the Rayleigh number — relative disk depth plane is constructed. It is shown that in the regular-oscillation regime the oscillation frequency is a linear function of the Rayleigh number.     

Bifurcations and chaos in parametrically excited single-degree-of-freedom systems

The behavior of single-degree-of-freedom systems possessing quadratic and cubic nonlinearities subject to parametric excitation is investigated. Both fundamental and principal parametric resonances are considered. A global bifurcation diagram in the excitation amplitude and excitation frequency domain is presented showing different possible stable steady-state solutions (attractors). Fractal basin maps for fundamental and principal parametric resonances when three attractors coexist are presented in color. An enlargement of one region of the map for principal parametric resonance reveals a Cantor-like set of fractal boundaries. For some cases, both periodic and chaotic attractors coexist.     

A coupled multigrid-domain-splitting technique for simulating incompressible flows in geometrically complex domains

This paper describes a domain decomposition numerical procedure for solving the Navier-Stokes equations in regions with complex geometries. The numerical method includes a modified version of QUICK (quadratic upstream interpolation convective kinematics) for the formulation of convective terms and a central difference scheme for the diffusion terms. A second-order-accurate predictor-corrector scheme is employed for the explicit time stepping. Although the momentum equations are solved independently on each subdomain, the pressure field is computed simultaneously on the entire flow field. A multigrid technique coupled with a Schwarz-like iteration method is devised to solve the pressure equation over the composite domains. The success of this strategy depends crucially on appropriate methods for specifying intergrid pressure boundary conditions on subdomains. A proper method for exchanging information among subdomains during the Schwarz sweep is equally important to the success of the multigrid solution for the overall pressure field. These methods are described and subsequently applied to two forced convection flow problems involving complex geometries to demonstrate the power and versatility of the technique. The resulting pressure and velocity fields exhibit excellent global consistency. The ability to simulate complex flow fields with this method provides a powerful tool for analysis and prediction of mixing and transport phenomenon.     

Experimental investigation of convective flows in a plane horizontal fluid layer heated from below and rotating about the vertical axis

The convective coherent structures in a plane horizontal fluid layer, heated from below and capable of rotation about the vertical axis, are experimentally investigated. It is shown that with increase in the supercriticality the time it takes for the convective structures to be formed decreases sharply. Rotation and an increase in the layer thickness-to-diameter ratio lead to an increase in the steady-state attainment time.     

A control volume finite‐element method for numerical simulating incompressible fluid flows without pressure correction

This paper presents a numerical model to study the laminar flows induced in confined spaces by natural convection. A control volume finite‐element method (CVFEM) with equal‐order meshing is employed to discretize the governing equations in the pressure–velocity formulation. In the proposed model, unknown variables are calculated in the same grid system using different specific interpolation functions without pressure correction. To manage memory storage requirements, a data storage format is developed for generated sparse banded matrices. The performance of various Krylov techniques, including Bi‐CGSTAB (Bi‐Conjugate Gradient STABilized) with an incomplete LU (ILU) factorization preconditioner is verified by applying it to three well‐known test problems. The results are compared to those of independent numerical or theoretical solutions in literature. The iterative computer procedure is improved by using a coupled strategy, which consists of solving simultaneously the momentum and the continuity equation transformed in a pressure equation. Results show that the strategy provides useful benefits with respect to both reduction of storage requirements and central processing unit runtime. Copyright © 2006 John Wiley & Sons, Ltd.     

基于力热耦合材料模型的Steven试验数值模拟方法

在力热耦合材料模型中,增加炸药自热放能模型,建立了Steven试验的力-热-化耦合的数值模拟方法。数值计算模型中,应力应变关系采用双线性硬化弹塑性模型,炸药受力后的热作用采用各向同性热材料模型,炸药的化学反应采用Arrhenius反应率函数,同时还考虑了升温和熔化对材料力学、热学性能的影响。针对标准Steven试验,通过数值分析得到了靶板的变形情况和炸药点火的速度阈值,将计算结果与实验数据进行了比较,两者符合较好。表明该方法可以较好地模拟Steven试验,而且与以往的分析模型和方法相比,本文的方法不需要增加经验性的点火准则和判据,具有更广泛的适用性,可以为研究低速撞击条件下炸药的力热响应和局域化点火问题提供参考。     

Coexistence of the chaos and the periodic solutions in planar fluid flows

This paper discusses the dynamic behavior of the Kelvin-Stuart cat’s eye flow underperiodic perturbations.By means of the Melnikov method the conditions to havebifurcations to subharmonics of even order for the oscillating orbits and to have bifurcationsto subharmonics of any order for the rotating orbits are given,and further,the coexistencephenomena of the chaotic motions and periodic solutions are presented.     

Chaos in the Convective Flow of a Fluid Mixture in a Porous Medium

The results of the study of the global behaviour of the convective flow of a binary mixture in a porous medium are presented. Bifurcation diagram, fixed points, periodic, chaotic solutions, stable and unstable manifolds, and basins of attraction have been calculated. Different behaviours (chaos, undecidable behaviour, etc.) have been found.     

气液两相流动与固壁相互作用耦合求解的研究

气液两相流动与固壁相互作用的研究是液滴撞击壁面运动研究的重要基础.以结合了VOF和Level Set两种方法优点的用于气液相界面追踪的复合Level Set-VOF方法和利用唯象分析方法建立的能够反映接触角滞后性及壁面性质对润湿过程影响的壁面润湿模型为基础,提出了气液两相流动与固壁相互作用耦合求解流程,给出了气液两相流动与固壁相互作用耦合求解过程中接触线速度的计算方法及边界条件的确定方法.通过与已有实验结果的对比,对提出的气液两相流动与固壁相互作用耦合求解方法的有效性进行了验证.     

SMALE HORSESHOES AND CHAOS IN DISCRETIZEDPERTURBED NLS SYSTEMS(Ⅰ)—POINCARÉ MAP

The existence of Smale horseshoes for a certain discretized perturbed nonlinear Schroedinger （NLS） equations was established by using n-dimensional versions of the Conley-Moser conditions. As a result, the discretized perturbed NLS system is shown to possess an invadant set A on which the dynamics is topologically conjugate to a shift on four symbols.     

A fully implicit scheme for simulating ionized gas flows using the gas dynamics electrodynamics coupled system

The flow of ionized gases under the influence of electromagnetic fields is governed by the coupled system of the compressible flow equations and the Maxwell equations. In this system, coupling of the flow with the electromagnetic field is obtained through nonlinear and stiff source terms, which may cause difficulties with the numerical solution of the coupled system. The discontinuous Galerkin finite element method is used for the numerical solution of this system. For the magnetic field vector, discontinuous Galerkin discretization is performed using a divergence‐free vector base for the magnetic field to preserve zero divergence in the element and retain the implicit constraint of a divergence‐free magnetic field vector down to very low level both globally and locally. To circumvent difficulties resulting from the presence of the stiff source terms, implicit time marching is used for the fully coupled system to avoid wrong wave shapes and propagation speeds that are obtained when the coupling source terms are lagged in time or by using splitting iterative schemes. Numerical solutions for benchmark problems computed on collocated meshes for the flow and electromagnetic field variables with this fully coupled monolithic approach showed good agreement with other numerical solutions and exact results. Copyright © 2014 John Wiley & Sons, Ltd.     

A simple method for simulating viscoelastic fluid flows in slit channel

A low-cost semi-analysis finite element technique, named the finite piece method (FPM) is presented in this article. It aims to solve three-dimensional (3D) viscoelastic slit flows. The viscoelastic stress of the fluid is modelled using an K-BKZ integral constitutive equation of the Wagner type. Picard iteration is used to solve non-linear equations. The FPM is tested on flow problems in both planar and contraction channels. The accuracy of the method is assessed by comparing flow distributions and pressure with results obtained by 3D finite element method (FEM). It shows that the solution accuracy is excellent and a substantial amount of computing time and memory requirement can be saved.     

Performance of RNG Turbulence Modelling for Forced Convective Heat Transfer in Ducts

This investigation concerns numerical calculation of turbulent forced convective heat transfer and fluid flow in straight ducts using the RNG (Re-Normalized Group) turbulence method.

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts with different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with the RNG κ?ε model and the RNG non-linear κ-ε model of Speziale. The turbulent heat fluxes are modeled by the simple eddy diffusivity (SED) concept, GGDH and WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models arc implemented for an arbitrary three dimensional duct.

Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non-staggered grid arrangement. The pressure-velocity coupling is handled by using the SIMPLEC-algorithm. The convective terms are treated by the QUICK, scheme while the diffusive terms are handled by the central-difference scheme. The hybrid scheme is used for solving the κ and ε equations.

The overall comparison between the models is presented in terms of friction factor and Nusselt number. The secondary flow generation is also of major concern.     

The analytical predictive criteria for chaos and escape in nonlinear oscillators: A survey

The purpose of this paper is to provide a brief summary of the various analytical predictive criteria in order for strange phenomena to occur in a class of softening nonlinear oscillators, oscillators which may exhibit escape from a potential well. Implications of Melnikov's criteria are discussed first and transient chaos in the twin-well potential oscillator is illustrated. Three different heuristic criteria for steady state chaos or escape solution, proposes by F. Moon, G. Schmidt and W. Szempliskia-Stupnicka, are then presented and compared to computer simulation results.     

Route to Chaos for Moderate Prandtl Number Convection in a Porous Layer Heated from Below

The route to chaos for moderate Prandtl number gravity driven convection in porous media is analysed by using Adomian's decomposition method which provides an accurate analytical solution in terms of infinite power series. The practical need to evaluate numerical values from the infinite power series, the consequent series truncation, and the practical procedure to accomplish this task, transform the otherwise analytical results into a computational solution achieved up to a desired but finite accuracy. The solution shows a transition to chaos via a period doubling sequence of bifurcations at a Rayleigh number value far beyond the critical value associated with the loss of stability of the convection steady solution. This result is extremely distinct from the sequence of events leading to chaos in low Prandtl number convection in porous media, where a sudden transition from steady convection to chaos associated with an homoclinic explosion occurs in the neighbourhood of the critical Rayleigh number (unless mentioned otherwise by 'the critical Rayleigh number' we mean the value associated with the loss of stability of the convection steady solution). In the present case of moderate Prandtl number convection the homoclinic explosion leads to a transition from steady convection to a period-2 periodic solution in the neighbourhood of the critical Rayleigh number. This occurs at a slightly sub-critical value of Rayleigh number via a transition associated with a period-1 limit cycle which seem to belong to the sub-critical Hopf bifurcation around the point where the convection steady solution looses its stability. The different regimes are analysed and periodic windows within the chaotic regime are identified. The significance of including a time derivative term in Darcy's equation when wave phenomena are being investigated becomes evident from the results.