Numerical Simulation of Saint-Venant Equations with Thermal Energy Dependency: Applications on Global Warming

Since the Industrial Revolution, humanity has been intensifying the burning of fossil fuels and as a consequence, the average temperature on Earth has been increasing. The 20th century was the warmest and future prospects are not favorable, that is, even higher temperatures are expected. This demonstrates the importance of studies on the subject, mainly to predict possible environmental, social and economic consequences. The objective of this work was to identify the interference of the increase in ambient temperature in the dynamics of fluids, such as ocean waves advancing over the continent. For this, thermal energy was considered in the Saint-Venant equations and computational implementations were performed via Lax-Friedrichs and Adams-Moulton methods. The results indicated that, in fact, depending on the amount of thermal energy transferred to the fluid, the advance of water towards the continent can occur, even in places where such a phenomenon has never been observed.

Modified Saint-Venant equations for flow simulation in tidal rivers

Flow in tidal rivers periodically propagates upstream or downstream under tidal influence. Hydrodynamic models based on the Saint-Venant equations （the SVN model） are extensively used to model tidal rivers. A force-corrected term expressed as the combination of flow velocity and the change rate of the tidal fevel was developed to represent tidal effects in the SVN model. A momentum equation incorporating with the corrected term was derived based on Newton＇s second law. By combing the modified momentum equation with the continuity equation, an improved SVN model for tidal rivers （the ISVN model） was constructed. The simulation of a tidal reach of the Qiantang River shows that the ISVN model performs better than the SVN model. It indicates that the corrected force derived for tidal effects is reasonable; the ISVN model provides an appropriate enhancement of the SVN model for flow simulation of tidal rivers.

Saint-Venant Equations and Friction Law for Modelling Self-Channeling Granular Flows: From Analogue to Numerical Simulation

Rock avalanches are catastrophic events involving important granular rock masses (>106 m3) and traveling long distances. In exceptional cases, the runout can reach up to tens of kilometers. Even if they are highly destructive and uncontrollable events, they give important insights to understand interactions between the displaced masses and landscape conditions. However, those events are not frequent. Therefore, the analogue and numerical modelling gives fundamental inputs to better understand their behavior. The objective of the research is to understand the propagation and spreading of granular mass released at the top of a simple geometry. The flow is unconfined, spreading freely along a 45° slope and deposit on a horizontal surface. The evolution of this analogue rock avalanche was measured from the initiation to its deposition with high speed camera. To simulate the analogue granular flow, a numerical model based on the continuum mechanics approach and the solving of the shallow water equations was used. In this model, the avalanche is described from a eulerian point of view within a continuum framework as single phase of incompressible granular material. The interaction of the flowing layer with the substratum follows a Mohr-Coulomb friction law. Within same initial conditions (slope, volume, basal friction, height of fall and initial velocity), results obtained with the numerical model are similar to those observed in the analogue. In both cases, the runout of the mass is comparable and the size of both deposits matches well. Moreover, both analogue and numerical modeling gave same magnitude of velocities. In this study, we highlighted the importance of the friction on a flowing mass and the influence of the numerical resolution on the propagation. The combination of the fluid dynamic equation with the frictional law enables the self-channelization and the stop of the granular mass.

Study of the 1D Saint-Venant Equations and Application to the Simulation of a Flood Problem

In this paper we are concerned with the mathematical and numerical analysis of the one-dimensional Saint-Venant equations. Thus, we prove the existence of a weak solution for any fixed time and with low regularity on the data. For the numerical approach we use the Rusanov scheme to approximate the flux and the hydrostatic reconstruction method which consists of decentering the source term at the interface. A numerical test of the proposed resolution is performed on a non-uniform topography.

Numerical solution of the Saint-Venant equations by an efficient hybrid finite-volume/finite-difference method

A computationally efficient hybrid finite-volume/finite-difference method is proposed for the numerical solution of SaintVenant equations in one-dimensional open channel flows. The method adopts a mass-conservative finite volume discretization for the continuity equation and a semi-implicit finite difference discretization for the dynamic-wave momentum equation. The spatial discretization of the convective flux term in the momentum equation employs an upwind scheme and the water-surface gradient term is discretized using three different schemes. The performance of the numerical method is investigated in terms of efficiency and accuracy using various examples, including steady flow over a bump, dam-break flow over wet and dry downstream channels, wetting and drying in a parabolic bowl, and dam-break floods in laboratory physical models. Numerical solutions from the hybrid method are compared with solutions from a finite volume method along with analytic solutions or experimental measurements. Comparisons demonstrates that the hybrid method is efficient, accurate, and robust in modeling various flow scenarios, including subcritical, supercritical, and transcritical flows. In this method, the QUICK scheme for the surface slope discretization is more accurate and less diffusive than the center difference and the weighted average schemes.

Some Modified Equations of the Sine-Hilbert Type

Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.

Data-Driven Ai-and Bi-Soliton of the Cylindrical Korteweg-de Vries Equation via Prior-Information Physics-Informed Neural Networks

By the modifying loss function MSE and training area of physics-informed neural networks(PINNs),we propose a neural networks model,namely prior-information PINNs(PIPINNs).We demonstrate the advantages of PIPINNs by simulating Ai-and Bi-soliton solutions of the cylindrical Korteweg-de Vries(cKdV)equation.

Matrix Riccati Equations in Optimal Control

In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.

THE STABILITY OF BOUSSINESQ EQUATIONS WITH PARTIAL DISSIPATION AROUND THE HYDROSTATIC BALANCE

This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).

The Maxwell-Heaviside Equations Explained by the Theory of Informatons

In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.

THE SMOOTHING EFFECT IN SHARP GEVREY SPACE FOR THE SPATIALLY HOMOGENEOUS NON-CUTOFF BOLTZMANN EQUATIONS WITH A HARDPOTENTIAL

In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.

On entire solutions of some Fermat type differential-difference equations

On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].

Equivalence between the internal observability and equation

This paper is concerned with a third order in time linear Moore-Gibson-Thompson equation which describes the acoustic velocity potential in ultrasound wave program.Influenced by the work of Kaltenbacher,Lasiecka and Marchand(Control Cybernet.2011,40:971-988),we establish an observability inequality of the conservative problem,and then discuss the equivalence between the exponential stabilization of a dissipative system and the internal observational inequality of the corresponding conservative system.

Multi-soliton solutions, breather-like and bound-state solitons for complex modified Korteweg–de Vries equation in optical fibers

Under investigation in this paper is a complex modified Korteweg–de Vries(KdV) equation, which describes the propagation of short pulses in optical fibers. Bilinear forms and multi-soliton solutions are obtained through the Hirota method and symbolic computation. Breather-like and bound-state solitons are constructed in which the signs of the imaginary parts of the complex wave numbers and the initial separations of the two parallel solitons are important factors for the interaction patterns. The periodic structures and position-induced phase shift of some solutions are introduced.

ON THE STABILITY OF PERIODIC SOLUTIONS OF PIECEWISE SMOOTH PERIODIC DIFFERENTIAL EQUATIONS

In this paper,we address the stability of periodic solutions of piecewise smooth periodic differential equations.By studying the Poincarémap,we give a sufficient condition to judge the stability of a periodic solution.We also present examples of some applications.

Exact and heuristic formulae to compute the geodetic height from the ellipse equation

The conversion of the cartesian coordinates of a point to its geodetic equivalent coordinates in reference to the geodetic ellipsoid is one of the main challenges in geodesy.The ellipse equation in the meridian plane significantly influences the value of the geodetic coordinates.This research analyzes this influence and how it can contribute to their solutions.The study investigates the mathematical relation between them and presents an exact formula relating to the geodetic height and the ellipse equation.In addition,a heuristic formula for the relation between the geodetic height and the ellipse equation is proposed,which is independent of the geodetic latitude and has a relative accuracy better than 99.9 %.The calculation is stable,and the cost is low.

A deep learning method based on prior knowledge with dual training for solving FPK equation

The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macroscopic variables in the stochastic dynamic system. Traditional methods for solving these equations often struggle with computational efficiency and scalability, particularly in high-dimensional contexts. To address these challenges, this paper proposes a novel deep learning method based on prior knowledge with dual training to solve the stationary FPK equations. Initially, the neural network is pre-trained through the prior knowledge obtained by Monte Carlo simulation(MCS). Subsequently, the second training phase incorporates the FPK differential operator into the loss function, while a supervisory term consisting of local maximum points is specifically included to mitigate the generation of zero solutions. This dual-training strategy not only expedites convergence but also enhances computational efficiency, making the method well-suited for high-dimensional systems. Numerical examples, including two different two-dimensional(2D), six-dimensional(6D), and eight-dimensional(8D) systems, are conducted to assess the efficacy of the proposed method. The results demonstrate robust performance in terms of both computational speed and accuracy for solving FPK equations in the first three systems. While the method is also applicable to high-dimensional systems, such as 8D, it should be noted that computational efficiency may be marginally compromised due to data volume constraints.

Sparse-Grid Implementation of Fixed-Point Fast Sweeping WENO Schemes for Eikonal Equations

Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.

On Two Types of Stability of Solutions to a Pair of Damped Coupled Nonlinear Evolution Equations

The stability of a set of spatially constant plane wave solutions to a pair of damped coupled nonlinear Schrödinger evolution equations is considered. The equations could model physical phenomena arising in fluid dynamics, fibre optics or electron plasmas. The main result is that any small perturbation to the solution remains small for all time. Here small is interpreted as being both in the supremum sense and the square integrable sense.

GLOBAL SOLUTIONS IN THE CRITICAL SOBOLEV SPACE FOR THE LANDAU EQUATION

The Landau equation is studied for hard potential with-2≤γ≤1.Under a perturbation setting,a unique global solution of the Cauchy problem to the Landau equation is established in a critical Sobolev space H_(x)^(d)L_(v)^(2)(d>3/2),which extends the results of[11]in the torus domain to the whole space R_(x)^(3).Here we utilize the pseudo-differential calculus to derive our desired result.