Mathematical analysis of the well known model of a piezoelectric energy harvester is presented. The harvester is designed as a cantilever Timoshenko beam with piezoelectric layers attached to its top and bottom faces. Thin, perfectly conductive electrodes are covering the top and bottom faces of the piezoelectric layers. These electrodes are connected to a resistive load. The model is governed by a system of three partial differential equations. The first two of them are the equations of the Timoshenko beam model and the third one represents Kirchhoff’s law for the electric circuit. All equations are coupled due to the piezoelectric effect. We represent the system as a single operator evolution equation in the Hilbert state space of the system. The dynamics generator of this evolution equation is a non-selfadjoint matrix differential operator with compact resolvent. The paper has two main results. Both results are explicit asymptotic formulas for eigenvalues of this operator, i.e., the modal analysis for the electrically loaded system is performed. The first set of the asymptotic formulas has remainder terms of the order O ( 1 n ) , where n is the number of an eigenvalue. These formulas are derived for the model with variable physical parameters. The second set of the asymptotic formulas has remainder terms of the order O ( 1 n 2 ) , and is derived for a less general model with constant parameters. This second set of formulas contains extra term depending on the electromechanical parameters of the model. It is shown that the spectrum asymptotically splits into two disjoint subsets, which we call the α -branch eigenvalues and the θ -branch eigenvalues. These eigenvalues being multiplied by “i” produce the set of the vibrational modes of the system. The α -branch vibrational modes are asymptotically located on certain vertical line in the left half of the complex plane and the θ -branch is asymptotically close to the imaginary axis. By having such spectral and asymptotic results, one can derive the asymptotic representation for the mode shapes and for voltage output. Asymptotics of vibrational modes and mode shapes is instrumental in the analysis of control problems for the harvester.
more »
« less
Asymptotic distribution of the eigenvalues of the bending‐torsion vibration model with fully nondissipative boundary feedback
Asymptotic and spectral results on the initial boundary‐value problem for the coupled bending‐torsion vibration model (which is important in such areas of engineering sciences as bridge and tall building designs, aerospace and oil pipes modeling, etc.) are presented. The model is given by a system of two hyperbolic partial differential equations equipped with a three‐parameter family of non‐self‐adjoint (linear feedback type) boundary conditions modeling the actions of self‐straining actuators. The system is rewritten in the form of the first‐order evolution equation in a Hilbert space of a four‐component Cauchy data. It is shown that the dynamics generator is a matrix differential operator with compact resolvent, whose discrete spectrum splits asymptotically into two disjoint subsets called the α‐branch and the β‐branch, respectively. Precise spectral asymptotics for the eigenvalues from each branch as the number of an eigenvalue tends to ∞ have been derived. It is also shown that the leading asymptotical term of the α‐branch eigenvalue depends only on the
- NSF-PAR ID:
- 10420587
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Studies in Applied Mathematics
- Volume:
- 150
- Issue:
- 4
- ISSN:
- 0022-2526
- Page Range / eLocation ID:
- p. 996-1025
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)The distribution of natural frequencies of the Euler–Bernoulli beam resting on elastic foundation and subject to an axial force in the presence of several damping mechanisms is investigated. The damping mechanisms are: ( i ) an external or viscous damping with damping coefficient ( − a 0 ( x )), ( ii ) a damping proportional to the bending rate with the damping coefficient a 1 ( x ). The beam is clamped at the left end and equipped with a four-parameter (α, β, κ 1 , κ 2 ) linear boundary feedback law at the right end. The 2 × 2 boundary feedback matrix relates the control input (a vector of velocity and its spacial derivative at the right end) to the output (a vector of shear and moment at the right end). The initial boundary value problem describing the dynamics of the beam has been reduced to the first order in time evolution equation in the state Hilbert space of the system. The dynamics generator has a purely discrete spectrum (the vibrational modes). Explicit asymptotic formula for the eigenvalues as the number of an eigenvalue tends to infinity have been obtained. It is shown that the boundary control parameters and the distributed damping play different roles in the asymptotical formulas for the eigenvalues of the dynamics generator. Namely, the damping coefficient a 1 and the boundary controls κ 1 and κ 2 enter the leading asymptotical term explicitly, while damping coefficient a 0 appears in the lower order terms.more » « less
-
Classical statistical mechanics has long relied on assumptions such as the equipartition theorem to understand the behavior of the complicated systems of many particles. The successes of this approach are well known, but there are also many well-known issues with classical theories. For some of these, the introduction of quantum mechanics is necessary, e.g., the ultraviolet catastrophe. However, more recently, the validity of assumptions such as the equipartition of energy in classical systems was called into question. For instance, a detailed analysis of a simplified model for blackbody radiation was apparently able to deduce the Stefan–Boltzmann law using purely classical statistical mechanics. This novel approach involved a careful analysis of a “metastable” state which greatly delays the approach to equilibrium. In this paper, we perform a broad analysis of such a metastable state in the classical Fermi–Pasta–Ulam–Tsingou (FPUT) models. We treat both the α-FPUT and β-FPUT models, exploring both quantitative and qualitative behavior. After introducing the models, we validate our methodology by reproducing the well-known FPUT recurrences in both models and confirming earlier results on how the strength of the recurrences depends on a single system parameter. We establish that the metastable state in the FPUT models can be defined by using a single degree-of-freedom measure—the spectral entropy (η)—and show that this measure has the power to quantify the distance from equipartition. For the α-FPUT model, a comparison to the integrable Toda lattice allows us to define rather clearly the lifetime of the metastable state for the standard initial conditions. We next devise a method to measure the lifetime of the metastable state tm in the α-FPUT model that reduces the sensitivity to the exact initial conditions. Our procedure involves averaging over random initial phases in the plane of initial conditions, the P1-Q1 plane. Applying this procedure gives us a power-law scaling for tm, with the important result that the power laws for different system sizes collapse down to the same exponent as Eα2→0. We examine the energy spectrum E(k) over time in the α-FPUT model and again compare the results to those of the Toda model. This analysis tentatively supports a method for an irreversible energy dissipation process suggested by Onorato et al.: four-wave and six-wave resonances as described by the “wave turbulence” theory. We next apply a similar approach to the β-FPUT model. Here, we explore in particular the different behavior for the two different signs of β. Finally, we describe a procedure for calculating tm in the β-FPUT model, a very different task than for the α-FPUT model, because the β-FPUT model is not a truncation of an integrable nonlinear model.more » « less
-
more » « less
Abstract Some physical processes that occur during a star's main-sequence evolution also affect its post-main-sequence evolution. It is well known that stars with masses above approximately 1.1
M ⊙have well-mixed convective cores on the main sequence; however, the structure of the star in the neighborhood of the convective core regions is currently underconstrained. We use asteroseismology to study the properties of the stellar core, in particular convective boundary mixing through convective overshoot, in such intermediate-mass stars. These core regions are poorly constrained by the acoustic (p ) mode oscillations observed for cool main-sequence stars. Consequently, we seek fossil signatures of main-sequence core properties during the subgiant and early first-ascent red giant phases of evolution. During these stages of stellar evolution, modes of mixed character that sample the deep interior can be observed. These modes sample the parts of the stars that are affected by the main-sequence structure of these regions. We model the global and near-core properties of 62 subgiant and early first-ascent red giant branch stars observed by theKepler , K2, and TESS space missions. We find that the effective overshoot parameter,α ov,eff, increases fromM = 1.0M ⊙toM = 1.2M ⊙before flattening out, although we note that the relationship betweenα ov,effand mass will depend on the incorporated modeling choices of internal physics and nuclear reaction network. We also situate these results within existing studies of main-sequence convective core boundaries. -
more » « less
Abstract Quartz crystals with zircon inclusions were synthesized using a piston-cylinder apparatus to experimentally evaluate the use of inclusions in “soft” host minerals for elastic thermobarometry. Synthesized zircon inclusion strains and, therefore, pressures (
P inc) were measured using Raman spectroscopy and then compared with the expected inclusion strains and pressures calculated from elastic models. Measured inclusion strains and inclusion pressures are systematically more tensile than the expected values and, thus, re-calculated entrapment pressures are overestimated. These discrepancies are not caused by analytical biases or assumptions in the elastic models and strain calculations. Analysis shows that inclusion strain discrepancies progressively decrease with decreasing experimental temperature in the α-quartz field. This behavior is consistent with inelastic deformation of the host–inclusion pairs induced by the development of large differential stresses during experimental cooling. Therefore, inclusion strains are more reliable for inclusions trapped at lower temperature conditions in the α-quartz field where there is less inelastic deformation of the host–inclusion systems. On the other hand, entrapment isomekes of zircon inclusions entrapped in the β-quartz stability field plot along the α–β quartz phase boundary, suggesting that the inclusion strains were mechanically reset at the phase boundary during experimental cooling and decompression. Therefore, inclusions contained in soft host minerals can be used for elastic thermobarometry and inclusions contained in β-quartz may provide constraints on theP –T at which the host–inclusion system crossed the phase boundary during exhumation.
