I'd like to just summarize the Lagrangian description of fluids quickly, because some texts have horrible typos and errors.
We begin with describing the position of fluid parcels using a function parametrized by position at time . Time evolution is described by the function
such that (locally) is invertible; i.e., we can write
Now we may introduce the Fluid Velocity in the Lagrangian description as
The Eulerian picture may be obtained by writing
Example 1. Consider a flow on the unit disc, with trajectories described by
Observe that and . We find we can invert this equation to determine the initial positions from the trajectory:
which permits us to write the Lagrangian parameters (initial positions) in terms of the current position of the fluid parcel. We find the time derivative of the trajectories
Plugging in the equations relating and in terms of and gives us (after some simple matrix multiplication):
This is a steady flow (since time doesn't explicit appear in the equation for or ), and just a rotation. If we want to generalize this, we could replace with a dimensionless, strictly increasing function of time , which would usually not be steady (e.g., if for positive constants , , then and ).
Exercise 1. Find the density for the flow from example 1 by solving the continuity equation, or some other way. If you do it some other way, prove your solution satisfies the equation of continuity.
Exercise 2. Is the flow from example 1 compressible or incompressible?
Example 2 (Childress, Example 2.2). Consider the one-dimensional flow with . Suppose we want to find the Lagrangian description of the flow, then we can rewrite this as
with initial condition . This has the obvious solution
This illustrates how to translate a solution in the Euler description to find the trajectories in the Lagrangian description.
Exercise 3. Solve the Euler flow equations and continuity equation for both the pressure and density. "Someone on the internet" asserts is a valid solution: is it? [Hint: use method of characteristics to determine density.]
No comments:
Post a Comment