| 1. |
- Bonito, Andrea, et al.
(författare)
-
Numerical Simulations of Surface Quasi-Geostrophic Flows on Periodic Domains
- 2021
-
Ingår i: SIAM Journal on Scientific Computing. - : Society for Industrial and Applied Mathematics. - 1064-8275 .- 1095-7197. ; 43:2, s. B405-B430
-
Tidskriftsartikel (refereegranskat)abstract
- We propose a novel algorithm for the approximation of surface quasi-geostrophic (SQG) flows modeled by a nonlinear partial differential equation coupling transport and fractional diffusion phenomena. The time discretization consists of an explicit strong-stability-preserving three-stage Runge-Kutta method while a flux-corrected-transport (FCT) method coupled with Dunford-Taylor representations of fractional operators is advocated for the space discretization. Standard continuous piecewise linear finite elements are employed, and the algorithm does not have restrictions on the mesh structure or on the computational domain. In the inviscid case, we show that the resulting scheme satisfies a discrete maximum principle property under a standard CFL condition and observe, in practice, its second order accuracy in space. The algorithm successfully approximates several benchmarks with sharp transitions and fine structures typical of SQG flows. In addition, theoretical Kolmogorov energy decay rates are observed on a freely decaying atmospheric turbulence simulation.
|
|
| 2. |
- Bouderlique, Thibault, et al.
(författare)
-
Surface flow for colonial integration in reef-building corals
- 2022
-
Ingår i: Current Biology. - : Elsevier. - 0960-9822 .- 1879-0445. ; 32:12, s. 2596-2609
-
Tidskriftsartikel (refereegranskat)abstract
- Reef-building corals are endangered animals with a complex colonial organization. Physiological mechanisms connecting multiple polyps and integrating them into a coral colony are still enigmatic. Using live imaging, particle tracking, and mathematical modeling, we reveal how corals connect individual polyps and form integrated polyp groups via species-specific, complex, and stable networks of currents at their surface. These currents involve surface mucus of different concentrations, which regulate joint feeding of the colony. Inside the coral, within the gastrovascular system, we expose the complexity of bidirectional branching streams that connect individual polyps. This system of canals extends the surface area by 4-fold and might improve communication, nutrient supply, and symbiont transfer. Thus, individual polyps integrate via complex liquid dynamics on the surface and inside the colony.
|
|
| 3. |
|
|
| 4. |
- Dao, Tuan Anh, et al.
(författare)
-
A High-Order Residual-Based Viscosity Finite Element Method for the Ideal MHD Equations
- 2022
-
Ingår i: Journal of Scientific Computing. - : Springer Nature. - 0885-7474 .- 1573-7691. ; 92:3
-
Tidskriftsartikel (refereegranskat)abstract
- We present a high order, robust, and stable shock-capturing technique for finite element approximations of ideal MHD. The method uses continuous Lagrange polynomials in space and explicit Runge-Kutta schemes in time. The shock-capturing term is based on the residual of MHD which tracks the shock and discontinuity positions, and adds sufficient amount of viscosity to stabilize them. The method is tested up to third order polynomial spaces and an expected fourth-order convergence rate is obtained for smooth problems. Several discontinuous benchmarks such as Orszag-Tang, MHD rotor, Brio-Wu problems are solved in one, two, and three spacial dimensions. Sharp shocks and discontinuity resolutions are obtained.
|
|
| 5. |
- Dao, Tuan Anh, et al.
(författare)
-
A nodal based high order nonlinear stabilization for finite element approximation of Magnetohydrodynamics
- 2024
-
Ingår i: Journal of Computational Physics. - 0021-9991 .- 1090-2716. ; 512, s. 113146-113146
-
Tidskriftsartikel (refereegranskat)abstract
- We present a novel high-order nodal artificial viscosity approach designed for solving Magnetohydrodynamics (MHD) equations. Unlike conventional methods, our approach eliminates the need for ad hoc parameters. The viscosity is mesh-dependent, yet explicit definition of the mesh size is unnecessary. Our method employs a multimesh strategy: the viscosity coefficient is constructed from a linear polynomial space constructed on the fine mesh, corresponding to the nodal values of the finite element approximation space. The residual of MHD is utilized to introduce high-order viscosity in a localized fashion near shocks and discontinuities. This approach is designed to precisely capture and resolve shocks. Then, high-order Runge-Kutta methods are employed to discretize the temporal domain. Through a comprehensive set of challenging test problems, we validate the robustness and high-order accuracy of our proposed approach for solving MHD equations.
|
|
| 6. |
- Dao, Tuan Anh, et al.
(författare)
-
A structure preserving numerical method for the ideal compressible MHD system
- 2024
-
Ingår i: Journal of Computational Physics. - 0021-9991 .- 1090-2716. ; 508, s. 113009-113009
-
Tidskriftsartikel (refereegranskat)abstract
- We introduce a novel structure-preserving method in order to approximate the compressible ideal Magnetohydrodynamics (MHD) equations. This technique addresses the MHD equations using a non-divergence formulation, where the contributions of the magnetic field to the momentum and total mechanical energy are treated as source terms. Our approach uses the Marchuk-Strang splitting technique and involves three distinct components: a compressible Euler solver, a source-system solver, and an update procedure for the total mechanical energy. The scheme allows for significant freedom on the choice of Euler's equation solver, while the magnetic field is discretized using a curl-conforming finite element space, yielding exact preservation of the involution constraints. We prove that the method preserves invariant domain properties, including positivity of density, positivity of internal energy, and the minimum principle of the specific entropy. If the scheme used to solve Euler's equation conserves total energy, then the resulting MHD scheme can be proven to preserve total energy. Similarly, if the scheme used to solve Euler's equation is entropy-stable, then the resulting MHD scheme is entropy stable as well. In our approach, the CFL condition does not depend on magnetosonic wave-speeds, but only on the usual maximum wavespeed from Euler's system. To validate the effectiveness of our method, we solve a variety of ideal MHD problems, showing that the method is capable of delivering second-order accuracy in space for smooth problems, while also offering unconditional robustness in the shock hydrodynamics regime as well.
|
|
| 7. |
- Dao, Tuan Anh, et al.
(författare)
-
Energy stable and accurate coupling of finite element methods and finite difference methods
- 2022
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 449
-
Tidskriftsartikel (refereegranskat)abstract
- We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are present. The proposed coupling technique requires minimal changes in the existing schemes while maintaining strict stability, accuracy, and energy conservation. Results are demonstrated on linear and nonlinear scalar conservation laws in two spatial dimensions.
|
|
| 8. |
- Dao, Tuan Anh, et al.
(författare)
-
Monolithic parabolic regularization of the MHD equations and entropy principles
- 2022
-
Ingår i: Computer Methods in Applied Mechanics and Engineering. - : Elsevier. - 0045-7825 .- 1879-2138. ; 398
-
Tidskriftsartikel (refereegranskat)abstract
- We show at the PDE level that the monolithic parabolic regularization of the equations of ideal magnetohydrodynamics (MHD) is compatible with all the generalized entropies, fulfills the minimum entropy principle, and preserves the positivity of density and internal energy. We then numerically investigate this regularization for the MHD equations using continuous finite elements in space and explicit strong stability preserving Runge–Kutta methods in time. The artificial viscosity coefficient of the regularization term is constructed to be proportional to the entropy residual of MHD. It is shown that the method has a high order of accuracy for smooth problems and captures strong shocks and discontinuities accurately for non-smooth problems.
|
|
| 9. |
|
|
| 10. |
- Guermond, Jean-Luc, et al.
(författare)
-
Finite element-based invariant-domain preserving approximation of hyperbolic systems : Beyond second-order accuracy in space
- 2024
-
Ingår i: Computer Methods in Applied Mechanics and Engineering. - : Elsevier. - 0045-7825 .- 1879-2138. ; 418
-
Tidskriftsartikel (refereegranskat)abstract
- This paper proposes an invariant-domain preserving approximation technique for nonlinear conservation systems that is high-order accurate in space and time. The algorithm mixes a highorder finite element method with an invariant-domain preserving low-order method that uses the closest neighbor stencil. The construction of the flux of the low-order method is based on an idea from Abgrall et al. (2017). The mass flux of the low-order and the high-order methods are identical on each finite element cell. This allows for mass preserving and invariant-domain preserving limiting.
|
|
