skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Physical Pendulum without Small Angle Approximation
The simple pendulum usually studied by analytic methods invokes the small angle approximation (SAA) so that one can easily reduce the equation of motion to Hooke’s law and thereby obtain the period of the pendulum and other associated quantities. If the approximation is relaxed, the problem becomes analytically intractable and one must resort to computational methods. In this exercise, the pendulum with and without the SAA are compared to allow students to discover what happens to the temporal behavior for larger angles of displacement. The students will also be able to obtain quantitative estimates of what a “small angle” means, and the limits of validity of the SAA. In the computations, students will learn that the second order differential equation that describes the motion of the pendulum can be reduced to two coupled first order differential equations, which can then be solved by the Euler-Cromer algorithm. In this specialized exercise set, students will also learn to use the ODE45 package in MATLAB to solve differential equations and its advantages over Euler-Cromer method, as well as ‘findpeaks’ command in MATLAB.  more » « less
Award ID(s):
2021209
PAR ID:
10554044
Author(s) / Creator(s):
;
Publisher / Repository:
American Association of Physics Teachers
Date Published:
Format(s):
Medium: X
Location:
Indiana University
Institution:
Partnership for the Integration of Computation into Undergraduate Physics
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    A reliable, accurate, and yet simple dynamic model is important to analyzing, designing, and controlling hybrid rigid–continuum robots. Such models should be fast, as simple as possible, and user-friendly to be widely accepted by the ever-growing robotics research community. In this study, we introduce two new modeling methods for continuum manipulators: a general reduced-order model (ROM) and a discretized model with absolute states and Euler–Bernoulli beam segments (EBA). In addition, a new formulation is presented for a recently introduced discretized model based on Euler–Bernoulli beam segments and relative states (EBR). We implement these models in a Matlab software package, named TMTDyn, to develop a modeling tool for hybrid rigid–continuum systems. The package features a new high-level language (HLL) text-based interface, a CAD-file import module, automatic formation of the system equation of motion (EOM) for different modeling and control tasks, implementing Matlab C-mex functionality for improved performance, and modules for static and linear modal analysis of a hybrid system. The underlying theory and software package are validated for modeling experimental results for (i) dynamics of a continuum appendage, and (ii) general deformation of a fabric sleeve worn by a rigid link pendulum. A comparison shows higher simulation accuracy (8–14% normalized error) and numerical robustness of the ROM model for a system with a small number of states, and computational efficiency of the EBA model with near real-time performances that makes it suitable for large systems. The challenges and necessary modules to further automate the design and analysis of hybrid systems with a large number of states are briefly discussed. 
    more » « less
  2. We consider the problem of stabilizing what we call a pendulum skate, a simple model of a figure skater developed by Gzenda and Putkaradze. By exploiting the symmetry of the system as well as taking care of the part of the symmetry broken by the gravity, the equations of motion are given as nonholonomic Euler–Poincaré equation with advected parameters. Our main interest is the stability of the sliding and spinning equilibria of the system. We show that the former is unstable and the latter is stable only under certain conditions. We use the method of Controlled Lagrangians to find a control to stabilize the sliding equilibrium. 
    more » « less
  3. This paper proves that the motion of small-slope vorticity fronts in the two-dimensional incompressible Euler equations is approximated on cubically nonlinear timescales by a Burgers–Hilbert equation derived by Biello and Hunter (2010) using formal asymptotic expansions. The proof uses a modified energy method to show that the contour dynamics equations for vorticity fronts in the Euler equations and the Burgers–Hilbert equation are both approximated by the same cubically nonlinear asymptotic equation. The contour dynamics equations for Euler vorticity fronts are also derived. 
    more » « less
  4. Quantum computers can produce a quantum encoding of the solution of a system of differential equations exponentially faster than a classical algorithm can produce an explicit description. However, while high-precision quantum algorithms for linear ordinary differential equations are well established, the best previous quantum algorithms for linear partial differential equations (PDEs) have complexity p o l y ( 1 / ϵ ) , where ϵ is the error tolerance. By developing quantum algorithms based on adaptive-order finite difference methods and spectral methods, we improve the complexity of quantum algorithms for linear PDEs to be p o l y ( d , log ⁡ ( 1 / ϵ ) ) , where d is the spatial dimension. Our algorithms apply high-precision quantum linear system algorithms to systems whose condition numbers and approximation errors we bound. We develop a finite difference algorithm for the Poisson equation and a spectral algorithm for more general second-order elliptic equations. 
    more » « less
  5. In this paper, we study ultra-weak discontinuous Galerkin methods with generalized numerical fluxes for multi-dimensional high order partial differential equations on both unstructured simplex and Cartesian meshes. The equations we consider as examples are the nonlinear convection-diffusion equation and the biharmonic equation. Optimal error estimates are obtained for both equations under certain conditions, and the key step is to carefully design global projections to eliminate numerical errors on the cell interface terms of ultra-weak schemes on general dimensions. The well-posedness and approximation capability of these global projections are obtained for arbitrary order polynomial space based on a wide class of generalized numerical fluxes on regular meshes. These projections can serve as general analytical tools to be naturally applied to a wide class of high order equations. Numerical experiments are conducted to demonstrate these theoretical results. 
    more » « less