To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
At first [Frodo] could see little. He seemed to be in a world of mist in which there were only shadows: the Ring was upon him. Then here and there the mist gave way and he saw many visions: small and clear as if they were under his eyes upon a table, and yet remote. There was no sound, only bright living images. The world seemed to have shrunk and fallen silent. He was sitting upon the Seat of Seeing, on Amon Hen, the Hill of the Eye of the Men of Númenor.
(J. R. R. Tolkien, The Lord of the Rings)
It is hard to imagine a time when our sight was not bombarded with images of far off people and places, when one had to actually be in the presence of someone to see their face without it being filtered through the hands of an artist (or while sitting on a Seat of Seeing). Images allow us to render the very small, the very large, and the very distant; they provide a visual record of time and place; they move us and entertain us; they replace perception with reality. Images are acquired by Earth-orbiting satellites, interplanetary probes, subterranean robots, surveillance cameras in unmanned aircraft and on nearly every street corner, and our mobile phones. They are created using optical principles and magnetic fields, and images display what is “seen” in every corner of the electromagnetic spectrum.
At a purchase cost of the order of, say, $300,000 or $400,000, the university can operate a computer for the education of its students quite economically… In a few years it is not unlikely that the computer may have settled immutably into our thought as an absolutely essential part of any university research program in the physical, psychological, and economic sciences.
(Alan J. Perlis, Computers and the World of the Future, 1962)
Before proceeding to the applications of variational methods that will occupy our attention throughout the remainder of the text, it is worthwhile to pause and consider how the solution of variational problems is accomplished when a closed form solution is not possible. We have already encountered several such scenarios, particularly when dealing with functionals with multiple independent variables that produce partial differential equations, and we will encounter numerous more situations in the balance of the text.
If an Euler equation cannot be solved in closed form, we may obtain approximate solutions using the following techniques:
Solve the differential Euler equation using numerical methods, such as finite difference methods, spectral methods, or finite-element methods, where the latter two are based on the Galerkin (or other method of weighted residual) approach.
Solve the integral variational form approximately using the Rayleigh-Ritz method or finite-element methods based on the Rayleigh-Ritz method.
The mathematician plays a game in which he himself invents the rules, while the physicist plays a game in which the rules are provided by nature, but as time goes on, it becomes increasingly evident that the rules which the mathematician finds interesting are the same as those which nature has chosen.
(Paul A. M. Dirac)
As described in Chapter 5, classical (Newtonian) mechanics is centered around Hamilton's principle and the associated Newton's second law, particularly as applied to statics and dynamics of nondeformable (discrete) and deformable (continuous) bodies. Such systems are typically of O(1) size operating at O(1) speeds. Modern physics generally refers to branches of physics that have arisen since Einstein's relativity theory in the early part of the twentieth century, including special relativity, general relativity, and quantum mechanics, for example. Modern physics addresses physical phenomena that occur at very small sizes, such as atoms, or very large sizes, such as galaxies, and/or very high speeds close to the speed of light. We focus on the theory of relativity and quantum mechanics, both of which reduce to the classical mechanics limit that has occupied our considerations thus far. Relativity reduces to classical mechanics when υ ≪ c, where c is the speed of light in a vacuum, and υ is a characteristic velocity of the system; quantum mechanics reduces to classical mechanics when length scales dominant within the system are much larger than Planck's universal constant ħ.
The miracle of the appropriateness of the language of mathematics to the formulation of the laws of physics is a wonderful gift which we neither understand nor deserve. We should be grateful for it and hope that it will remain valid in future research and that it will extend, for better or worse, to our pleasure, even though perhaps also to our bafflement, to wide branches of learning.
(Eugene Wigner)
Electromagnetic wave radiation is remarkable for two reasons: 1) it spans a vast spectrum of wavelengths providing for a wide range of physics and applications, and 2) it does not require a medium within which to travel. While all electromagnetic waves travel at the speed of light if traveling through a vacuum, they span ten orders of magnitude in wavelength, essentially corresponding to the full range of sizes that exist in the universe. Amazingly, the governing equations of electromagnetic waves, which do not require a medium, and physical waves, which do, are essentially the same. Rather than a guitar string or gas molecules, however, electromagnetic waves travel via massless photons. While the photons do not have mass, they do carry energy; it is this energy that determines their wavelength. A photon's energy is directly proportional to its frequency (and inversely proportional to the wavelength).
There bought a space of ground, which (Byrsa call'd,
From the bull's hide) they first inclos'd, and wall'd.
(Virgil, The Aeneid)
The above excerpt cryptically recounts the legend of Dido, the Queen of Carthage, from the ninth century bc. After being exiled from Tyre in Lebanon, Dido purportedly sailed to the shores of North Africa (now Tunisia) and requested that the local inhabitants give her and her party the land that could be enclosed by the hide of an ox. Not thinking that an oxhide could encompass a large portion of land, they granted her wish. She then proceeded to have the hide cut into narrow strips and extended end-to-end to form a semicircle bounding the shoreline and encompassing a nearby hill. This area became known as Carthage.
Certain branches of mathematics have arisen out of consideration of the theoretical consequences of known mathematical theorems. More often than not, however, new branches of mathematics have been developed to provide the means to address certain types of practical problems that initially are often very specific. For example, how did Dido know to arrange the oxhide in a circular shape in order to enclose the largest possible area with a given perimeter length?
We have to make every effort to understand for ourselves what the dangers are, and this points up a fundamental thing about computers: They involve more thought and not less thought. They may save certain parts of our efforts, but they do not eliminate the need for intelligence.
(Norbert Wiener, Computers and the World of the Future, 1962)
Our final application of optimization theory is in grid generation, which is an important topic in many computational fields in which we desire to obtain numerical solutions of partial differential equations with complex physics and/or on domains with complex shapes. When choosing a numerical grid on which to solve differential equations, it is necessary that the grid both faithfully represent the geometry of the domain and be sufficiently refined in order to maintain the numerical truncation errors at an acceptable level. Traditionally, this has been done in a rather ad hoc manner by choosing a grid that “looks good” in some qualitative sense. Calculus of variations allows us to optimize the choice of a grid in a formal manner by explicitly enforcing certain criteria on the generation of the grid, thereby providing a more intuitive and mathematically formal basis for grid generation. Variational grid generation is generally formulated to produce the “best” grid in a least squares sense, eliminating the trial and error necessary in other grid generation approaches, such as algebraic and elliptic grid generation.
I am an old man now, and when I die and go to heaven there are two matters on which I hope for enlightenment. One is quantum electrodynamics, and the other is the turbulent motion of fluids. And about the former I am rather optimistic.
(Sir Horace Lamb)
Fluid mechanics is one of the fields that is contributing significantly to the current resurgence of interest in variational methods. Historically, calculus of variations has not been taught or utilized as a core mathematical tool in the arsenal of the fluid mechanics practitioner; it is a rare fluid mechanics textbook that even mentions variational calculus. However, recent research developments are revolutionizing the field by reframing certain fluid mechanics phenomena within a variational framework. This is particularly the case in the areas of flow control and hydrodynamic stability.
Traditionally, flow control has been implemented in an ad hoc manner involving a significant amount of trial and error within an empirical (experimental or numerical) framework. Beginning in the 1990s, flow control is increasingly being framed within the context of optimal control theory. As will be seen in Chapter 10, the solution to optimal control formulations, particularly for large problems involving many degrees of freedom, is very computationally intensive. It is only recently that the computational resources have become available that are capable of solving realistic fluid mechanics scenarios within such an optimal control framework.
Similarly, significant advances are being made in our understanding of stability of fluid flows through the application of transient growth analysis that seeks the “optimal,” or most unstable, initial perturbation (disturbance) that results in the greatest growth of the instability.
Catbert: Wally, you can't float through life with no goals and no ambition. Wally: You misjudge me. I have my entire career planned out. My five-year plan is to avoid any sort of work in which my individual accomplishments can be measured. I'll hoard knowledge about one of our legacy systems so I seem indispensable. When I get to within four years of retirement, I'll only work on projects that have a five-year payback. I'll protect my cardiovascular system by getting plenty of naps and not caring about the quality of my work. Then I'll stick a straw in our pension fund and suck on it for the next forty years.
(Dilbert Comic by Scott Adams)
The focus of Chapters 4 through 9 has been to develop physical variational principles and their associated differential equations of mathematical physics. The number and variety of physical phenomena captured in Hamilton's principle are truly remarkable and are a testament to the fundamental nature of the law of conservation of energy on which it is based. Whereas our emphasis has been to elucidate the laws of physics using calculus of variations, we now turn our attention to optimization and control of these systems. Interestingly, the same variational methods used to formulate the governing Euler equations provide the framework for performing their optimization and control. The governing Euler-Lagrange equations of motion obtained in the previous chapters become the differential constraints for control of such dynamical, fluid mechanical, and electromagnetic systems.
Now, here is this principle, so wise, so worthy of the Supreme Being: when some change occurs in Nature, the amount of Action used for this change is always the smallest possible.
The laws of movement thus deduced [from the principle of least action], being found to be precisely the same as those observed in nature, we can admire the application of it to all phenomena, in the movement of animals, in the vegetation of plants, in the revolution of the heavenly bodies: and the spectacle of the universe becomes so much the grander, so much the more beautiful, so much more worthy of its Author, when one knows that a small number of laws, most wisely established, suffice for all movements.
(Pierre Louis Moreau de Maupertuis)
Classical mechanics encompasses those areas of physics that originated prior to the development of relativistic mechanics at the beginning of the twentieth century. Its primary focus is on the application of Newtonian mechanics to macroscopic systems. Classical mechanics provides the basis for many of the important fields in engineering, including solid mechanics, fluid mechanics, transport phenomena, and dynamics. These fields are employed broadly in the design of almost all devices that make modern life possible, including being sure that your mobile phone can withstand a fall on a hard surface, placing and maintaining communication satellites in their orbits, and both terrestrial and extraterrestrial transportation systems.
Steady water infiltration in homogeneous soils is governed by the Richards equation. This equation can be studied more conveniently by transforming to a type of Helmholtz equation. In this study, a dual-reciprocity boundary element method (DRBEM) is employed to solve the Helmholtz equation numerically. Using the solutions obtained, numerical values of the suction potential are then computed. The proposed method is tested on problems involving infiltration from different types of periodic channels in a homogeneous soil. Moreover, the method is also examined using infiltration from periodic trapezoidal channels in three different types of homogeneous soil.
This paper presents an integrated guidance and control (IGC) design method for an unmanned aerial vehicle with static stability which is described by a nonlinear six-degree-of-freedom (6-DOF) model. The model is linearized by using small disturbance linearization. The dynamic characteristics of pitching mode, rolling mode and Dutch rolling mode are obtained by analysing the linearized model. Furthermore, an IGC design procedure is also proposed in conjunction with a proportional–integral–derivative (PID) control method and fuzzy control method. A PID controller is applied in the control loop of the elevator and aileron, and the attitude angle and attitude angular velocity are used as compensation feedback, giving a simple and low-order control law. A fuzzy control method is applied to perform the cross-coupling control of rolling and yawing. Finally, the 6-DOF simulation shows the effectiveness of the developed method.
In this article we consider elliptic partial differential equations with random coefficients and/or random forcing terms. In the current treatment of such problems by stochastic Galerkin methods it is standard to assume that the random diffusion coefficient is bounded by positive deterministic constants or modeled as lognormal random field. In contrast, we make the significantly weaker assumption that the non-negative random coefficients can be bounded strictly away from zero and infinity by random variables only and may have distributions different from a lognormal one. We show that in this case the standard stochastic Galerkin approach does not necessarily produce a sequence of approximate solutions that converges in the natural norm to the exact solution even in the case of a lognormal coefficient. By using weighted test function spaces we develop an alternative stochastic Galerkin approach and prove that the associated sequence of approximate solutions converges to the exact solution in the natural norm. Hereby, ideas for the case of lognormal coefficient fields from earlier work of Galvis, Sarkis and Gittelson are used and generalized to the case of positive random coefficient fields with basically arbitrary distributions.
Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of their time about the minima of the underlying potential, occasionally transitioning into different basins of attraction. The essential idea of parallel replica dynamics is that the exit distribution from a given well for a single process can be approximated by the distribution of the first exit of N independent identical processes, each run for only 1 /N-th the amount of time. While promising, this leads to a series of numerical analysis questions about the accuracy of the exit distributions. Building upon the recent work in [C. Le Bris, T. Lelièvre, M. Luskin and D. Perez, Monte Carlo Methods Appl. 18 (2012) 119–146], we prove a unified error estimate on the exit distributions of the algorithm against an unaccelerated process. Furthermore, we study a dephasing mechanism, and prove that it will successfully complete.
In this work, the error behaviour of high-order exponential operator splitting methods for the time integration of nonlinear evolutionary Schrödinger equations is investigated. The theoretical analysis utilises the framework of abstract evolution equations on Banach spaces and the formal calculus of Lie derivatives. The general approach is substantiated on the basis of a convergence result for exponential operator splitting methods of (nonstiff) order p applied to the multi-configuration time-dependent Hartree–Fock (MCTDHF) equations, which are associated with a model reduction for high-dimensional linear Schrödinger equations describing free electrons that interact by Coulomb force. Provided that the analytical solution of the MCTDHF equations constituting a system of coupled linear ordinary differential equations and low-dimensional nonlinear partial differential equations satisfies suitable regularity requirements, convergence of order p − 1 in the H1 Sobolev norm and convergence of order p in the L2 norm is proven. An analogous result follows for the cubic nonlinear Schrödinger equation, which is also illustrated by a numerical experiment.