Lagrangian Mechanics and The Heat Equation
I have been noticing a lot of discussion on Lucia’s blog about an energy Balance Climate Model.
Two Box Model: Algebra for ver. 1 test.
Two Box Model: Now assuming surface temperature are ‘mixed’ values.
Two Box Model: Rough idea how to obtain parameters.
Two Box Models & The 2nd Law of Thermodynamics
Which started based on a post by Tamino, called “Two Box”
The biggest criticism I’ve scene about these models is that there is very little physics which is used to form these models and consequently they won’t capture much of the nonlinear dynamics.
The basic premises of the model is that the time lags are associated with the heat capacity of the ocean. Consequently, the ocean and some of the atmosphere is partitioned into boxes and the model tries to determine how these boxes are coupled.
The end result though is the time legs end up being associated with eignvectors which, project onto both boxes. Consequently the modes of the system do not coincide with the boxes. Moreover, it is desirable to minimize the number of model parameters which need to be fit and therefore principles of physics, should be used to derive a more realistic “Two box model” or should I say more generally “Two mode model” as the modes do not coincide with the boxes.
The most simple equation with regards to the transfer of energy is the heat equation:
1)
The heat equation is essentially a diffusion equation based on Brownian Motion. But we can adjust the constant to try and account for other heat exchange processes.
In one dimension the heat equation can be written as:
2)
In this equation the time derivative only depends on the spatial derivatives. Keeping a mind on the dynamics, the higher the frequency that the climate forcing is, the more shallow it will penetrate into the ocean because the more ocean you include the greater the thermal inertia (heat capacity) and consequently the greater the damping. Thus a good fit for the temperature response at a given frequency as a function of depth, may be an exponential function of the depth.
Thus, from this observation and because the mathematics are fairly simple consider a set of basis functions which the special component is exponentially decreasing from the surface. That is let:
3)
Plugging this into the heat equation one gets:
4)
Where:
The constant represents how quickly heat diffuses in the ocean at a given depth. Equating terms on the heat equation gives:
5)
Now let there be a function , and let the quanity being diffused in the heat equation be the temperature. Then the energy constraint is given by:
6)
This constraint makes one of the generalized coordinates redundant.
Differentiating (5) with respect to time, gives:
7)
is determined by the forcing. Now some notes, with regards to the coordinate system:
1)I suggest taking the as the redundant basis function, the one that has the smallest value .
2)Since this model is a linearization (or approximation), I suggest using as the temperature variable at a given depth, the amount by which it exceeds the average temperature at that depth. That is model the temperature anomaly instead of the actual temperature.
3)Given two then we are looking at the change in the heat transfer from the mean instead of the actual heat transfer.
4)These are only suggestions.
Now the Lagrangian, is essentially a cost function. In our cost function, we want to minimize the error in the derivatives (the square of equation (5) integrated over y) and as an additional constraint, I want to minimize the change in entropy with respect to y.
When a parcel of water rises because it is hotter then the surrounding and therefore less dense it will expand with constant entropy (adiabatic expansion) if it exchanges no heat with the surroundings. Therefore, the heat induced ocean currents will act in order to try and minimize the entropy change between adjacent regions. (Note the actual effect of this process is to increase entropy).
From thermodynamics:
8 )
Therefore:
9)
from which we can get:
10)
Substituting into (10) the expression for temperature (equation (3)) equation (10) gives:
11)
The entropy cost function is:
12)
The Lagrangian is given by:
13)
Where is a stiffness parameter for the entropy.
Equation 13) is subject to the constraints for some (These are the energy constraints mentioned above.
14)
15)
Now the Euler Lagrange Equation is used to obtain the dynamics (differential equations) from the Lagrangian
16)
This will give a second order differential equation and it might be necessary to use some linear algebra to rearrange the equation. You can convert this into a first order differential equation by using the Hamiltonian form. If only three basis functions are used then one gets “two mode model”. There are no restrictions on the number of bais functions used.
7 Comments »
Leave a Reply

Recent
 Laplace Transform Via Limits
 log(CO2) and Scary Graphs
 Numeric Solutions to The Heat Equation
 Coriolis Forces in Hopkins and Simmons Vorticity Equation
 The Cross Product in Non Orthogonal Coordinate Systems
 Lagrangian Mechanics and The Heat Equation
 Laplace Transform of f(t) Related to smoothed f(t)?
 Coriolis Forces
 Vector Operations in Hoskins and Simmons Coordinates
 API/Object Viewers/Memory Mapping/
 Defining a Microsoft access Datasource
 Fractal Modeling of Turbulence

Links

Archives
 July 2012 (1)
 September 2009 (5)
 August 2009 (19)
 March 2009 (2)

Categories

RSS
Entries RSS
Comments RSS
I saw a link for vertical ocean mixing posted on Lucia’s blog. It looked interesting:
http://oceanworld.tamu.edu/resources/ocng_textbook/chapter08/chapter08_05.htm
Also, I haven’t taken the above equations beyond this post but it looks like if I include my entropy condition in the Lagrangian then the equations become nonlinear but if I leave the entropy condition out the equations are linear.
Comment by s243a  September 1, 2009 
Just as a side note, I realized after posting this thread that the two box model is also an approximation to the heat equation. I discuss the case of overlapping boxes here:
Overlapping Boxes
Also of interest is a post I saw on the CrankNicolson method which is the standard technique for solving the heat equation:
Approximations used in CrankNicolson method for solving PDEs numerically
Comment by s243a  September 8, 2009 
I just wanted to share a good post with regards to modeling that I read over at the blackboard:
http://rankexploits.com/musings/2009/arthurscase2ithink/
Comment by s243a  September 10, 2009 
I’m just going to post this link here for future reference. It is an old post I made about approximating nonlinear systems using the matrix exponential
http://www.physicsforums.com/showthread.php?t=325690
Comment by s243a  September 10, 2009 
Here is some good links on ocean dynamics I’ll want to look at later:
Ocean Dynamics
Comment by s243a  September 12, 2009 
[…] suggests my choice of a negative exponential basis in my post (Lagrangian Mechanics and the Heat Equation) was not two bad a choice although, Guasian functions will decay slightly faster then negative […]
Pingback by Numeric Solutions to The Heat Equation « Earth Cubed  September 15, 2009 
I don’t know If I said it already but …I’m so glad I found this site…Keep up the good work I read a lot of blogs on a daily basis and for the most part, people lack substance but, I just wanted to make a quick comment to say GREAT blog. Thanks, 🙂
A definite great read..Tony Brown
Comment by Tony Brown  September 23, 2009 