skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: A remark on the invariant energy quadratization (IEQ) method for preserving the original energy dissipation laws
In this letter, we revisit the invariant energy quadratization (IEQ) method and provide a new perspective on its ability to preserve the original energy dissipation laws. The IEQ method has been widely used to design energy stable numerical schemes for phase-field or gradient flow models. Although there are many merits of the IEQ method, one major disadvantage is that the IEQ method usually respects a modified energy law, where the modified energy is expressed in the auxiliary variables. Still, the dissipation laws in terms of the original energy are not guaranteed by the IEQ method. Using the widely-used Cahn-Hilliard equation as an example, we demonstrate that the Runge-Kutta IEQ method indeed can preserve the original energy dissipation laws for certain situations up to arbitrary high-order accuracy. Interested readers are encouraged to extend this idea to more general cases and apply it to other thermodynamically consistent models.  more » « less
Award ID(s):
2111479
PAR ID:
10337445
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Electronic Research Archive
Volume:
30
Issue:
2
ISSN:
2688-1594
Page Range / eLocation ID:
701 to 714
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Computers & Mathematics with Applications)
    In this paper, we propose a linear second-order numerical method for solving the Allen-Cahn equation with general mobility. The fully-discrete scheme is achieved by using the Crank-Nicolson formula for temporal integration and the central difference method for spatial approximation, together with two additional stabilization terms. Under mild constraints on the two stabilizing parameters, the proposed numerical scheme is shown to unconditionally preserve the discrete maximum bound principle and the discrete original energy dissipation law. Error estimate in the 𝐿∞ norm is successfully derived for the proposed scheme. Finally, some numerical experiments are conducted to verify the theoretical results and demonstrate the performance of the proposed scheme in combination with an adaptive time-stepping strategy. 
    more » « less
  2. ; (, ESAIM: Mathematical Modelling and Numerical Analysis)
    We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization Method (IEQ). The model consists of the Navier–Stokes equations for the fluid flow, coupled to the Q-tensor gradient flow describing the liquid crystal molecule alignment. By using the Invariant Energy Quadratization Method, we obtain a linearly implicit scheme, accelerating the computational speed. However, this introduces an auxiliary variable to replace the bulk potential energy and it isa prioriunclear whether the reformulated system is equivalent to the Beris-Edward system. In this work, we prove stability properties of the scheme and show its convergence to a weak solution of the coupled liquid crystal system. We also demonstrate the equivalence of the reformulated and original systems in the weak sense. 
    more » « less
  3. ; (, Proceedings of the European Wave and Tidal Energy Conference)
    Wave energy converters typically use various control methods to extract energy from ocean waves. The objective of the control system is to optimize the energy extraction process, taking into account the dynamics of the system and the wave conditions. The task of deriving the optimal control laws of wave energy converter arrays for regular and irregular waves using the Pontryagin minimum principle was previously investigated in the literature. The result is a combination between the singular arc and bang-bang control laws. For irregular waves, some complexity arises due to the radiation state-space representation, which requires ignoring the hydrodynamic coupling terms related to the added mass and radiation-damping coefficients; this reduces the computational complexity of the control force but adversely affects the solution's accuracy. Also, the derived control laws are specific to a particular wave condition. Recently, the optimal control of a flexible buoy wave energy converter was derived using the convolution representation for the radiation force. In this work, the optimal control laws of flexible buoy wave energy converters are modified to simulate wave energy converter arrays; then, the results are compared to those obtained by dropping the hydrodynamic radiation coupling terms. Although using a convolution representation adds computational complexity to the optimal control problem, it generates an equation that is generic to any wave condition, can be used with any wave spectrum, and provides an expression for the switching condition. 
    more » « less
  4. ; ; ; ; ; ; (, Journal of Physical Oceanography)
    Abstract Oceanic mixing, mostly driven by the breaking of internal waves at small scales in the ocean interior, is of major importance for ocean circulation and the ocean response to future climate scenarios. Understanding how internal waves transfer their energy to smaller scales from their generation to their dissipation is therefore an important step for improving the representation of ocean mixing in climate models. In this study, the processes leading to cross-scale energy fluxes in the internal wave field are quantified using an original decomposition approach in a realistic numerical simulation of the California Current. We quantify the relative contribution of eddy–internal wave interactions and wave–wave interactions to these fluxes and show that eddy–internal wave interactions are more efficient than wave–wave interactions in the formation of the internal wave continuum spectrum. Carrying out twin numerical simulations, where we successively activate or deactivate one of the main internal wave forcing, we also show that eddy–near-inertial internal wave interactions are more efficient in the cross-scale energy transfer than eddy–tidal internal wave interactions. This results in the dissipation being dominated by the near-inertial internal waves over tidal internal waves. A companion study focuses on the role of stimulated cascade on the energy and enstrophy fluxes. 
    more » « less
  5. ; (, Proceedings of the Royal Society of Edinburgh: Section A Mathematics)
    null (Ed.)
    Abstract Smooth solutions of the incompressible Euler equations are characterized by the property that circulation around material loops is conserved. This is the Kelvin theorem. Likewise, smooth solutions of Navier–Stokes are characterized by a generalized Kelvin's theorem, introduced by Constantin–Iyer (2008). In this note, we introduce a class of stochastic fluid equations, whose smooth solutions are characterized by natural extensions of the Kelvin theorems of their deterministic counterparts, which hold along certain noisy flows. These equations are called the stochastic Euler–PoincarĂ© and stochastic Navier–Stokes–PoincarĂ© equations respectively. The stochastic Euler–PoincarĂ© equations were previously derived from a stochastic variational principle by Holm (2015), which we briefly review. Solutions of these equations do not obey pathwise energy conservation/dissipation in general. In contrast, we also discuss a class of stochastic fluid models, solutions of which possess energy theorems but do not, in general, preserve circulation theorems. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email