Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 99
PROCEEDINGS OF THE ELEVENTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY Edited by: B.H.V. Topping
Paper 253
Optimal Prediction of Lyapunov Exponents for a Duffing Oscillator A. LaBryer, P.J. Attar and P. Vedula
University of Oklahoma, Norman OK, United States of America A. LaBryer, P.J. Attar, P. Vedula, "Optimal Prediction of Lyapunov Exponents for a Duffing Oscillator", in B.H.V. Topping, (Editor), "Proceedings of the Eleventh International Conference on Computational Structures Technology", CivilComp Press, Stirlingshire, UK, Paper 253, 2012. doi:10.4203/ccp.99.253
Keywords: multiscale modeling, reduced order modeling, optimal prediction, Duffing oscillator.
Summary
A novel multiscale reduced order modeling (ROM) method called optimal spatiotemporal reduced order modeling (OPSTROM), is presented in this paper, which was recently proposed by the authors [1,2,3,4], and can be applied to any nonlinear dynamical system. The governing equations are modified for an underresolved numerical simulation in space and time with an arbitrary discretization scheme. Basic filtering concepts are used to demonstrate the manner in which subgridscale dynamics arise with a coarse computational grid. Models are then developed to account for the underlying spatiotemporal structure using inclusion of statistical information into the governing equations on a multipoint stencil. These subgridscale models are designed to provide closure by accounting for the interactions between the spatiotemporal microscales and macroscales as the system evolves. Predictions for the modified system are based upon principles of meansquare error minimization, conditional expectations and stochastic estimation, thus rendering the optimal solution with respect to the chosen resolution.
The OPSTROM approach has practical applications for a wide range of multiscale problems in physics and engineering. When fully resolved simulations are not feasible, the technique can be used to expedite simulations by coarsening the computational grid in space and time while maintaining reliable predictions for the quantities of interest. Three variants of the technique are conceivable: OPTROM (which accounts for the subgrid temporal structure), OPSROM (which accounts for the subgrid spatial structure), and OPSTROM (which accounts for the complete subgrid spatiotemporal structure). The OPTROM variant was originally proposed by the authors [1,2] whereas OPSROM is analogous to the numerical method known as largeeddy simulation for turbulent flows. The primary objective of this paper is to test the ability of the OPTROM approach to improve the accuracy of Lyapunov exponents for a system in which the response may be regular or chaotic. The Duffing oscillator, has been selected for study, which contains no spatial dependence, as it enables exclusive focus on the effects of a coarse temporal grid. Maps are drawn for the Lyapunov exponents in a large forcing parameter space in which the response may undergo symmetric timeperiodicity, asymmetry, period Ntupling and chaos. The exponents are used to discriminate between regular and chaotic response types. In order to expedite simulations, coarse timesteps are used with an implicit timemarching scheme. Substantial discretization errors occur in the form of strong numerical damping, which produces artificial regularity in some of the numerical solutions. In order to improve the accuracy of the underresolved simulations, optimal subgridscale models are developed efficiently by means of direct calculation and estimation [3]. By resolving more time scales with OPTROM, underresolved predictions for the Lyapunov exponents significantly improve. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
