京都大学 / 新制・課程博士 / 博士(理学) / 甲第23717号 / 理博第4807号 / 新制||理||1688(附属図書館) / 京都大学大学院理学研究科化学専攻 / (主査)教授 谷村 吉隆, 教授 林 重彦, 教授 渡邊 一也 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DGAM

<p>L^p is supposed to be L with a superscript lower case 'p.'</p> / <p>This thesis investigates a 1-dimensional Gross-Pitaevskii (GP) equation from the viewpoint of a system of Hamiltonian partial differential equations (PDEs). A theorem on Birkhoff normal forms is a particularly important goal of this study. The resulting system is a perturbed system of a completely resonant system, which we analyze, using several forms of perturbation theory.</p> <p>In chapter two, we study estimates 011 integrals of products of four Hermite functions, which represent coefficients of mode coupling, and play an important role in the proof of the Birkhoff normal form theorem. This is a basic problem, which has a close relationship with a problem of Besicovitch, namely the behavior of the L^p norms of L² -normalized Hermite functions.</p> <p>In chapter three we carefully reconsider the linear Schrodinger equation with a harmonic potential, and we introduce a family of Hilbert spaces for studying the GP equation, which generalize the traditional energy spaces in which one works. One unexpected fact is that these function spaces have a close relationship with the former works for the tempered distributions, in particular the N-representation theory due to B. Simon, and V. Bargmann's theory, which uncovers relationship between the tempered distributions and his function spaces through the so-called Segal-Bargmann transformation. In addition, our function spaces have a nice relationship with the Sobolev spaces. In this chapter, a few other questions regarding these function spaces are discussed.</p> <p>In chapter four the proof of the Birkhoff normal form theorem on spaces we have introduced are provided. The analysis is divided into two cases according to the regularity of the related function space. After proving the Birkhoff normal form theorem, we made an analysis of the impact of the perturbation on the main part of the GP system, which we remark is completel:y resonant.</p> / Doctor of Philosophy (PhD)

yes / Traditional rotor dynamics mainly focuses on the steady- state behavior of the rotor and shafting. However, for systems such as hydro turbine generating sets (HTGS) where the control and regulation is frequently applied, the shafting safety and stabilization in transient state is then a key factor. The shafting transient state inevitably involves multiparameter domain, multifield coupling, and coupling dynamics. In this paper, the relative value form of the Lagrange function and its equations have been established by defining the base value system of the shafting. Takingthe rotation angle and the angular speed of the shafting as a link, the shafting lateral vibration and generator equations are integrated into the framework of generalized Hamiltonian system. The generalized Hamiltonian control model is thus established. To make the model more general, additional forces of the shafting are taken as the input excitation in proposed model. The control system of the HTGS can be easily connected with the shafting model to form the whole simulation system of the HTGS. It is expected that this study will build a foundation for the coupling dynamics theory using the generalized Hamiltonian theory to investigate coupling dynamic mechanism among the shafting vibration, transient of hydro turbine generating sets, and additional forces of the shafting. / National Natural Science Foundation of China under Grant Nos. 51179079 and 50839003

No / Under local multi-machine condition, voltage and current changes of generator d and q axis are important parameters that mirror the coupling characteristics between local generator and others. Fifth order generator Hamiltonian model included electromagnetic transient of the d and q axis is established in this paper, which provides a foundation that reveals inner dynamics mechanism of multi-machine coupling and interaction. Based on electromechanical analysis dynamics theory, each of subsystem Lagrange function, dissipative function and generalized forces can be derived from fundamental energy relationship of the generator, then the Lagrange-Maxwell equations of the generator can be derived. The equations are transformed into the generalized Hamiltonian model by defining the generalized momentum, in which its structure is clear, its energy flow is consistent with actual physical system ones. The Hamiltonian model of generator is improved by the transition to that with practical parameters resulting from local simplified method. The conception of the local multi-machine system of power station is proposed, implicit multi-machine system Hamiltonian model described by measurable parameters on the generator port is derived. Analysis show that proposed model is well in structure.

Yes / A diesel engine is a typical dynamic system. In this paper, a dynamics method is proposed to establish the Hamiltonian model of the diesel engine, which solves the main difficulty of con-structing a Hamiltonian function under the multi-field coupling condition. Furthermore, the control method of Hamiltonian model structure modification is introduced to study the control of a diesel engine. By means of the principle of energy-shaping and Hamiltonian model structure modification theories, the modified energy function is constructed, which is proved to be a quasi-Lyapunov function of the closed-loop system. Finally, the control laws are derived, and the simulations are carried out. The study reveals the dynamic mechanism of diesel engine operation and control and provides a new way to research the modeling and control of a diesel engine system. / This research was funded by the National Natural Science Foundation of China, grant number 51869007.

This thesis deals with non-linear non-quadratic optimal control problems in an autonomous system and a related iterative numerical method, the Kleinman-Newton method, for solving the problem. The thesis proves the local convergence of Kleinman-Newton method using the contraction mapping theorem and then describes how this Kleinman-Newton method may be used to numerically solve for the optimal control and the corresponding solution. In order to show the proof and the related numerical work, it is necessary to review some of earlier work in the beginning of Chapter 1 [Zhang], and to introduce the Kleinman-Newton method at the end of the chapter. In Chapter 2 we will demonstrate the proof. In Chapter 3 we will show the related numerical work and results. / Ph. D.

Ensemble and variational techniques have gained wide popularity as the two main approaches for solving data assimilation and inverse problems. The majority of the methods in these two approaches are derived (at least implicitly) under the assumption that the underlying probability distributions are Gaussian. It is well accepted, however, that the Gaussianity assumption is too restrictive when applied to large nonlinear models, nonlinear observation operators, and large levels of uncertainty. This work develops a family of fully non-Gaussian data assimilation algorithms that work by directly sampling the posterior distribution. The sampling strategy is based on a Hybrid/Hamiltonian Monte Carlo (HMC) approach that can handle non-normal probability distributions. The first algorithm proposed in this work is the "HMC sampling filter", an ensemble-based data assimilation algorithm for solving the sequential filtering problem. Unlike traditional ensemble-based filters, such as the ensemble Kalman filter and the maximum likelihood ensemble filter, the proposed sampling filter naturally accommodates non-Gaussian errors and nonlinear model dynamics, as well as nonlinear observations. To test the capabilities of the HMC sampling filter numerical experiments are carried out using the Lorenz-96 model and observation operators with different levels of nonlinearity and differentiability. The filter is also tested with shallow water model on the sphere with linear observation operator. Numerical results show that the sampling filter performs well even in highly nonlinear situations where the traditional filters diverge. Next, the HMC sampling approach is extended to the four-dimensional case, where several observations are assimilated simultaneously, resulting in the second member of the proposed family of algorithms. The new algorithm, named "HMC sampling smoother", is an ensemble-based smoother for four-dimensional data assimilation that works by sampling from the posterior probability density of the solution at the initial time. The sampling smoother naturally accommodates non-Gaussian errors and nonlinear model dynamics and observation operators, and provides a full description of the posterior distribution. Numerical experiments for this algorithm are carried out using a shallow water model on the sphere with observation operators of different levels of nonlinearity. The numerical results demonstrate the advantages of the proposed method compared to the traditional variational and ensemble-based smoothing methods. The HMC sampling smoother, in its original formulation, is computationally expensive due to the innate requirement of running the forward and adjoint models repeatedly. The proposed family of algorithms proceeds by developing computationally efficient versions of the HMC sampling smoother based on reduced-order approximations of the underlying model dynamics. The reduced-order HMC sampling smoothers, developed as extensions to the original HMC smoother, are tested numerically using the shallow-water equations model in Cartesian coordinates. The results reveal that the reduced-order versions of the smoother are capable of accurately capturing the posterior probability density, while being significantly faster than the original full order formulation. In the presence of nonlinear model dynamics, nonlinear observation operator, or non-Gaussian errors, the prior distribution in the sequential data assimilation framework is not analytically tractable. In the original formulation of the HMC sampling filter, the prior distribution is approximated by a Gaussian distribution whose parameters are inferred from the ensemble of forecasts. The Gaussian prior assumption in the original HMC filter is relaxed. Specifically, a clustering step is introduced after the forecast phase of the filter, and the prior density function is estimated by fitting a Gaussian Mixture Model (GMM) to the prior ensemble. The base filter developed following this strategy is named cluster HMC sampling filter (ClHMC ). A multi-chain version of the ClHMC filter, namely MC-ClHMC , is also proposed to guarantee that samples are taken from the vicinities of all probability modes of the formulated posterior. These methodologies are tested using a quasi-geostrophic (QG) model with double-gyre wind forcing and bi-harmonic friction. Numerical results demonstrate the usefulness of using GMMs to relax the Gaussian prior assumption in the HMC filtering paradigm. To provide a unified platform for data assimilation research, a flexible and a highly-extensible testing suite, named DATeS , is developed and described in this work. The core of DATeS is implemented in Python to enable for Object-Oriented capabilities. The main components, such as the models, the data assimilation algorithms, the linear algebra solvers, and the time discretization routines are independent of each other, such as to offer maximum flexibility to configure data assimilation studies. / Ph. D.

This dissertation presents innovative unified approaches to understand and predict the motion between potential wells. The theoretical-computational framework, based on the tube dynamics, will reveal how the dissipative and gyroscopic forces change the phase space structure that governs the escape (or transition) from potential wells. In higher degree of freedom systems, the motion between potential wells is complicated due to the existence of multiple escape routes usually through an index-1 saddle. Thus, this dissertation firstly studies the local behavior around the index-1 saddle to establish the criteria of escape taking into account the dissipative and gyroscopic forces. In the analysis, an idealized ball rolling on a surface is selected as an example to show the linearized dynamics due to its special interests that the gyroscopic force can be easily introduced by rotating the surface. Based on the linearized dynamics, we find that the boundary of the initial conditions of a given energy for the trajectories that transit from one side of a saddle to the other is a cylinder and ellipsoid in the conservative and dissipative systems, respectively. Compared to the linear systems, it is much more challenging or sometimes impossible to get analytical solutions in the nonlinear systems. Based on the analysis of linearized dynamics, the second goal of this study is developing a bisection method to compute the transition boundary in the nonlinear system using the dynamic snap-through buckling of a buckled beam as an example. Based on the Euler-Bernoulli beam theory, a two degree of freedom Hamiltonian system can be generated via a two mode-shape truncation. The transition boundary on the Poincar'e section at the well can be obtained by the bisection method. The numerical results prove the efficiency of the bisection method and show that the amount of trajectories that escape from the potential well will be smaller if the damping of the system is increasing. Finally, we present an alternative idea to compute the transition boundary of the nonlinear system from the perspective of the invariant manifold. For the conservative systems, the transition boundary of a given energy is the invariant manifold of a periodic orbit. The process of obtaining such invariant manifold compromises two parts, including the computation of the periodic orbit by solving a proper boundary-value problem (BVP) and the globalization of the manifold. For the dissipative systems, however, the transition boundary of a given energy becomes the invariant manifold of an index-1 saddle. We present a BVP approach using the small initial sphere in the stable subspace of the linearized system at one end and the energy at the other end as the boundary conditions. By using these algorithms, we obtain the nonlinear transition tube and transition ellipsoid for the conservative and dissipative systems, respectively, which are topologically the same as the linearized dynamics. / Doctor of Philosophy / Transition or escape events are very common in daily life, such as the snap-through of plant leaves and the flipping over of umbrellas on a windy day, the capsize of ships and boats on a rough sea. Some other engineering problems related to escape, such as the collapse of arch bridges subjected to seismic load and moving trucks, and the escape and recapture of the spacecraft, are also widely known. At first glance, these problems seem to be irrelated. However, from the perspective of mechanics, they have the same physical principle which essentially can be considered as the escape from the potential wells. A more specific exemplary representative is a rolling ball on a multi-well surface where the potential energy is from gravity. The purpose of this dissertation is to develop a theoretical-computational framework to understand how a transition event can occur if a certain energy is applied to the system. For a multi-well system, the potential wells are usually connected by saddle points so that the motion between the wells generally occurs around the saddle. Thus, knowing the local behavior around the saddle plays a vital role in understanding the global motion of the nonlinear system. The first topic aims to study the linearized dynamics around the saddle. In this study, an idealized ball rolling on both stationary and rotating surfaces will be used to reveal the dynamics. The effect of the gyroscopic force induced by the rotation of the surface and the energy dissipation will be considered. In the second work, the escape dynamics will be extended to the nonlinear system applied to the snap-through of a buckled beam. Due to the nonlinear behavior existing in the system, it is hard to get the analytical solutions so that numerical algorithms are needed. In this study, a bisection method is developed to search the transition boundary. By using such method, the transition boundary on a specific Poincar'e section is obtained for both the conservative and dissipative systems. Finally, we revisit the escape dynamics in the snap-through buckling from the perspective of the invariant manifold. The treatment for the conservative and dissipative systems is different. In the conservative system, we compute the invariant manifold of a periodic orbit, while in the dissipative system we compute the invariant manifold of a saddle point. The computational process for the conservative system consists of the computation of the periodic orbit and the globalization of the corresponding manifold. In the dissipative system, the invariant manifold can be found by solving a proper boundary-value problem. Based on these algorithms, the nonlinear transition tube and transition ellipsoid in the phase space can be obtained for the conservative and dissipative systems, respectively, which are qualitatively the same as the linearized dynamics.

The purpose of this thesis was to analyze President Trump's Afghanistan foreign policy and to determine if it fits the mold of one of the four historical foreign policy paradigms as described by Walter Russell Mead in his book, Special Providence: American Foreign Policy and How It Changed the World. Mead describes four U.S. foreign policy schools of thought, in which he titles after influential statesmen who embody the specific school's core principles. These paradigms include the Hamiltonians, who believe in a strong relationship between big business and government for foreign policy success. The Wilsonians, who encourage the spread of democratic principles abroad. The Jeffersonians, who favor the protection of domestic liberal democracy over other foreign policy endeavors. Lastly, the Jacksonians, who prioritize the physical and economic security of American citizens above all else. The primary research question in this thesis states, which of the four traditions of U.S. foreign policy identified by Walter Russell Mead helps explain President Trump's Afghanistan foreign policy? President Trump's rhetoric and specific foreign policy actions were analyzed. His rhetoric was examined through his August 2017 Afghanistan Strategy speech and his specific foreign policy actions were measured through various air operation metrics, U.S. aid to Afghanistan, and U.S. troop deployment trends. Overall, this thesis gave support to my hypothesis that President Trump's Afghanistan foreign policy contains various Hamiltonian and Wilsonian principles, but it has proven to be principally Jacksonian. / Master of Arts / The objective of this thesis was to further understand President Trump's Afghanistan foreign policy through the theoretical framework of Walter Russel Mead's four historical foreign policy paradigms. Mead's four historical paradigms are useful tools to examine and understand U.S. foreign policy. Mead provides in-depth historical context, goes into great detail on core principles, and also provides a surfeit of advantages and disadvantages for each school of thought. His breakdown of U.S. foreign policy into complementary yet combative paradigms is one of the most complete explanations of U.S. foreign policy to date. The primary research question in this thesis states, which of the four traditions of U.S. foreign policy identified by Walter Russell Mead helps explain President Trump's Afghanistan foreign policy? President Trump's rhetoric and specific foreign policy actions were analyzed. His rhetoric was examined through his August 2017 Afghanistan Strategy speech and his specific foreign policy actions were measured through various air operation metrics, U.S. aid to Afghanistan, and U.S. troop deployment trends. Overall, this thesis gave support to my hypothesis that President Trump's Afghanistan foreign policy contains various Hamiltonian and Wilsonian principles, but it has proven to be principally Jacksonian.

In this dissertation, the extended Melnikov's method has been applied to several nonlinear ship dynamics models, which are related to the new generation of stability criteria in the International Maritime Organization (IMO). The advantage of this extended Melnikov's method is it overcomes the limitation of small damping that is intrinsic to the implementation of the standard Melnikov's method. The extended Melnikv's method is first applied to two published roll motion models. One is a simple roll model with nonlinear damping and cubic restoring moment. The other is a model with a biased restoring moment. Numerical simulations are investigated for both models. The effectiveness and accuracy of the extended Melnikov's method is demonstrated. Then this method is used to predict more accurately the threshold of global surf-riding for a ship operating in steep following seas. A reference ITTC ship is used here by way of example and the result is compared to that obtained from previously published standard analysis as well as numerical simulations. Because the primary drawback of the extended Melnikov's method is the inability to arrive at a closed form equation, a 'best fit'approximation is given for the extended Melnikov numerically predicted result. The extended Melnikov's method for slowly varying system is applied to a roll-heave-sway coupled ship model. The Melnikov's functions are calculated based on a fishing boat model. And the results are compared with those from standard Melnikov's method. This work is a preliminary research on the application of Melnikov's method to multi-degree-of-freedom ship dynamics. In the last part of the dissertation, the method of manufactured solution is applied to systems with chaotic behavior. The purpose is to identify points with potential numerical discrepancies, and to improve computational efficiency. The numerical discrepancies may be due to the selection of error tolerances, precisions, etc. Two classical chaotic models and two ship capsize models are examined. The current approach overlaps entrainment in chaotic control theory. Here entrainment means two dynamical systems have the same period, phase and amplitude. The convergent region from control theory is used to give a rough guideline on identifying numerical discrepancies for the classical chaotic models. The effectiveness of this method in improving computational efficiency is demonstrated for the ship capsize models. / Ph. D.