Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Cauchy-Lagrangian method for numerical analysis of Euler flow
218
0
0.0
(
0
)
Added by
Vladislav Zheligovsky A.
Publication date
2015
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A novel semi-Lagrangian method is introduced to solve numerically the Euler equation for ideal incompressible flow in arbitrary space dimension. It exploits the time-analyticity of fluid particle trajectories and requires, in principle, only limited spatial smoothness of the initial data. Efficient generation of high-order time-Taylor coefficients is made possible by a recurrence relation that follows from the Cauchy invariants formulation of the Euler equation (Zheligovsky & Frisch, J. Fluid Mech. 2014, 749, 404-430). Truncated time-Taylor series of very high order allow the use of time steps vastly exceeding the Courant-Friedrichs-Lewy limit, without compromising the accuracy of the solution. Tests performed on the two-dimensional Euler equation indicate that the Cauchy-Lagrangian method is more - and occasionally much more - efficient and less prone to instability than Eulerian Runge-Kutta methods, and less prone to rapid growth of rounding errors than the high-order Eulerian time-Taylor algorithm. We also develop tools of analysis adapted to the Cauchy-Lagrangian method, such as the monitoring of the radius of convergence of the time-Taylor series. Certain other fluid equations can be handled similarly.
rate research
Read More
This Response is concerned with the recent Comment of Ruiz-Herrera, Limitations of the Method of Lagrangian Descriptors [arXiv:1510.04838], criticising the method of Lagrangian Descriptors. In spite of the significant body of literature asserting the contrary, Ruiz-Herrera claims that the method fails to reveal the presence of stable and unstable manifolds of hyperbolic trajectories in incompressible systems and in almost all linear systems. He supports this claim by considering the method of Lagrangian descriptors applied to three specific examples. However in this response we show that Ruiz-Herrera does not understand the proper application and interpretation of the method and, when correctly applied, the method beautifully and unambiguously detects the stable and unstable manifolds of the hyperbolic trajectories in his examples.
We study the dynamics of a droplet moving on an inclined rough surface in the absence of inertial and viscous stress effects. In this case, the dynamics of the droplet is a purely geometric motion in terms of the wetting domain and the capillary surface. Using a single graph representation, we interpret this geometric motion as a gradient flow on a Hilbert manifold. We propose unconditionally stable first/second order numerical schemes to simulate this geometric motion of the droplet, which is described using motion by mean curvature coupled with moving contact lines. The schemes are based on (i) explicit moving boundaries, which decouple the dynamic updates of the contact lines and the capillary surface, (ii) a semi-Lagrangian method on moving grids and (iii) a predictor-corrector method with a nonlinear elliptic solver upto second order accuracy. For the case of quasi-static dynamics with continuous spatial variable in the numerical schemes, we prove the stability and convergence of the first/second order numerical schemes. To demonstrate the accuracy and long-time validation of the proposed schemes, several challenging computational examples - including breathing droplets, droplets on inhomogeneous rough surfaces and quasi-static Kelvin pendant droplets - are constructed and compared with exact solutions to quasi-static dynamics obtained by desingularized differential-algebraic system of equations (DAEs).
In the past decades, boreal summers have been characterized by an increasing number of extreme weather events in the Northern Hemisphere extratropics, including persistent heat waves, droughts and heavy rainfall events with significant social, economic and environmental impacts. Many of these events have been associated with the presence of anomalous large-scale atmospheric circulation patterns, in particular persistent blocking situations, i.e., nearly stationary spatial patterns of air pressure. To contribute to a better understanding of the emergence and dynamical properties of such situations, we construct complex networks representing the atmospheric circulation based on Lagrangian trajectory data of passive tracers advected within the atmospheric flow. For these Lagrangian flow networks, we study the spatial patterns of selected node properties prior to, during and after different atmospheric blocking events in Northern Hemisphere summer. We highlight the specific network characteristics associated with the sequence of strong blocking episodes over Europe during summer 2010 as an illustrative example. Our results demonstrate the ability of the node degree, entropy and harmonic closeness centrality based on outgoing links to trace important spatio-temporal characteristics of atmospheric blocking events. In particular, all three measures capture the effective separation of the stationary pressure cell forming the blocking high from the normal westerly flow and the deviation of the main atmospheric currents around it. Our results suggest the utility of further exploiting the Lagrangian flow network approach to atmospheric circulation in future targeted diagnostic and prognostic studies.
Global spectral analysis (GSA) is used as a tool to test the accuracy of numerical methods with the help of canonical problems of convection and convection-diffusion equation which admit exact solutions. Similarly, events in turbulent flows computed by direct numerical simulation (DNS) are often calibrated with theoretical results of homogeneous isotropic turbulence due to Kolmogorov, as given in Turbulence -U. Frisch, Cambridge Univ. Press, UK (1995). However, numerical methods for the simulation of this problem are not calibrated, as by using GSA of convection and/or convection-diffusion equation. This is with the exception in A critical assessment of simulations for transitional and turbulence flows-Sengupta, T.K., In Proc. of IUTAM Symp. on Advances in Computation, Modeling and Control of Transitional and Turbulent Flows, pp 491-532, World Sci. Publ. Co. Pte. Ltd., Singapore (2016), where such a calibration has been advocated with the help of convection equation. For turbulent flows, an extreme event is characterized by the presence of length scales smaller than the Kolmogorov length scale, a heuristic limit for the largest wavenumber present without being converted to heat. With growing computer power, recently many simulations have been reported using a pseudo-spectral method, with spatial discretization performed in Fourier spectral space and a two-stage, Runge-Kutta (RK2) method for time discretization. But no analyses are reported to ensure high accuracy of such simulations. Here, an analysis is reported for few multi-stage Runge-Kutta methods in the Fourier spectral framework for convection and convection-diffusion equations. We identify the major source of error for the RK2-Fourier spectral method using GSA and also show how to avoid this error and specify numerical parameters for achieving highest accuracy possible to capture extreme events in turbulent flows.
We derive a numerical method for Darcy flow, hence also for Poissons equation in mixed (first order) form, based on discrete exterior calculus (DEC). Exterior calculus is a generalization of vector calculus to smooth manifolds and DEC is one of its discretizations on simplicial complexes such as triangle and tetrahedral meshes. DEC is a coordinate invariant discretization, in that it does not depend on the embedding of the simplices or the whole mesh. We start by rewriting the governing equations of Darcy flow using the language of exterior calculus. This yields a formulation in terms of flux differential form and pressure. The numerical method is then derived by using the framework provided by DEC for discretizing differential forms and operators that act on forms. We also develop a discretization for spatially dependent Hodge star that varies with the permeability of the medium. This also allows us to address discontinuous permeability. The matrix representation for our discrete non-homogeneous Hodge star is diagonal, with positive diagonal entries. The resulting linear system of equations for flux and pressure are saddle type, with a diagonal matrix as the top left block. The performance of the proposed numerical method is illustrated on many standard test problems. These include patch tests in two and three dimensions, comparison with analytically known solution in two dimensions, layered medium with alternating permeability values, and a test with a change in permeability along the flow direction. We also show numerical evidence of convergence of the flux and the pressure. A convergence experiment is also included for Darcy flow on a surface. A short introduction to the relevant parts of smooth and discrete exterior calculus is included in this paper. We also include a discussion of the boundary condition in terms of exterior calculus.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
