اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدA hyperbolic reformulation of the Serre-Green-Naghdi model for general bottom topographies
96
0
0.0
(
0
)
تأليف
Caterina Bassi
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We present a novel hyperbolic reformulation of the Serre-Green-Naghdi (SGN) model for the description of dispersive water waves. Contrarily to the classical Boussinesq-type models, it contains only first order derivatives, thus allowing to overcome the numerical difficulties and the severe time step restrictions arising from higher order terms. The proposed model reduces to the original SGN model when an artificial sound speed tends to infinity. Moreover, it is endowed with an energy conservation law from which the energy conservation law associated with the original SGN model is retrieved when the artificial sound speed goes to infinity. The governing partial differential equations are then solved at the aid of high order ADER discontinuous Galerkin finite element schemes. The new model has been successfully validated against numerical and experimental results, for both flat and non-flat bottom. For bottom topographies with large variations, the new model proposed in this paper provides more accurate results with respect to the hyperbolic reformulation of the SGN model with the mild bottom approximation recently proposed in C. Escalante, M. Dumbser and M.J. Castro. An efficient hyperbolic relaxation system for dispersive non-hydrostatic water waves and its solution with high order discontinuous Galerkin schemes, Journal of Computational Physics 2018.
قيم البحث
اقرأ أيضاً
A new modified Galerkin / Finite Element Method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low-order Lagrange finite element spaces, despite the fact that the
system contains third order spatial partial derivatives for the depth averaged velocity of the fluid. After studying the efficacy and the conservation properties of the new numerical method, we proceed with the validation of the new numerical model and boundary conditions by comparing the numerical solutions with laboratory experiments and with available theoretical asymptotic results.
In this work we present an adaptive boundary element method for computing the electromagnetic response of wave interactions in hyperbolic metamaterials. One unique feature of hyperbolic metamaterial is the strongly directional wave in its propagating
cone, which induces sharp transition for the solution of the integral equation across the cone boundary when wave starts to decay or grow exponentially. In order to avoid a global refined mesh over the whole boundary, we employ a two-level a posteriori error estimator and an adaptive mesh refinement procedure to resolve the singularity locally for the solution of the integral equation. Such an adaptive procedure allows for the reduction of the degree of freedom significantly for the integral equation solver while achieving desired accuracy for the solution. In addition, to resolve the fast transition of the fundamental solution and its derivatives accurately across the propagation cone boundary, adaptive numerical quadrature rules are applied to evaluate the integrals for the stiff matrices. Finally, in order to formulate the integral equations over the boundary, we also derive the limits of layer potentials and their derivatives in the hyperbolic media when the target points approach the boundary.
Kaniadakis and Quarati (1994) proposed a Fokker--Planck equation with quadratic drift as a PDE model for the dynamics of bosons in the spatially homogeneous setting. It is an open question whether this equation has solutions exhibiting condensates in
finite time. The main analytical challenge lies in the continuation of exploding solutions beyond their first blow-up time while having a linear diffusion term. We present a thoroughly validated time-implicit numerical scheme capable of simulating solutions for arbitrarily large time, and thus enabling a numerical study of the condensation process in the Kaniadakis--Quarati model. We show strong numerical evidence that above the critical mass rotationally symmetric solutions of the Kaniadakis--Quarati model in 3D form a condensate in finite time and converge in entropy to the unique minimiser of the natural entropy functional at an exponential rate. Our simulations further indicate that the spatial blow-up profile near the origin follows a universal power law and that transient condensates can occur for sufficiently concentrated initial data.
The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov--Poisson--Lenard--Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is
decomposed as $f=mathcal{E}[boldsymbol{rho}_{f}]+g$, where $mathcal{E}$ is a local equilibrium distribution, depending on the macroscopic moments $boldsymbol{rho}_{f}=int_{mathbb{R}}boldsymbol{e} fdv=langleboldsymbol{e} frangle_{mathbb{R}}$, where $boldsymbol{e}=(1,v,frac{1}{2}v^{2})^{rm{T}}$, and $g$, the microscopic distribution, is defined such that $langleboldsymbol{e} grangle_{mathbb{R}}=0$. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for $boldsymbol{rho}_{f}$ and $g$. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the $langleboldsymbol{e} grangle_{mathbb{R}}=0$ constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. We present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for $f$.
We present both modeling and computational advancements to a unique three-dimensional discontinuous Petrov-Galerkin finite element model for the simulation of laser amplification in a fiber amplifier. Our model is based on the time-harmonic Maxwell e
quations, and it incorporates both amplification via an active dopant and thermal effects via coupling with the heat equation. As a full vectorial finite element simulation, this model distinguishes itself from other fiber amplifier models that are typically posed as an initial value problem and make significantly more approximations. Our model supports co-, counter-, and bi-directional pumping configurations, as well as inhomogeneous and anisotropic material properties. The longer-term goal of this modeling effort is to study nonlinear phenomena that prohibit achieving unprecedented power levels in fiber amplifiers, along with validating typical approximations used in lower-fidelity models. The high-fidelity simulation comes at the cost of a high-order finite element discretization with many degrees of freedom per wavelength. This is necessary to counter the effect of numerical pollution due to the high-frequency nature of the wave simulation. To make the computation more feasible, we have developed a novel longitudinal model rescaling, using artificial material parameters with the goal of preserving certain quantities of interest. Our numerical tests demonstrate the applicability and utility of this scaled model in the simulation of an ytterbium-doped, step-index fiber amplifier that experiences laser amplification and heating. We present numerical results for the nonlinear coupled Maxwell/heat model with up to 240 wavelengths.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
