Hyperbolic partial differential equations describe wave behavior.

learn more… | top users | synonyms

2
votes
0answers
42 views

Euler Equation Eigensystem with Gravity in the Energy Flux

I am modifying a conservative form of the Euler equations with gravity in the energy flux (see previous question: Energy Conservation in Conservation Laws with Source Terms) for use in a Riemann ...
1
vote
1answer
81 views

How can I solve this nonlinear, variable-coefficient hyperbolic PDE in one space dimension?

I need to solve the following hyperbolic equation in x and phi co-ordinates $$\frac{\partial \left ( -s/f \right )}{\partial \varphi }+\frac{\partial \left ( 1/f \right )}{\partial x}=0$$ $$\varphi ...
3
votes
2answers
99 views

How to set up a shock tube problem such that the solution includes a shock with a specified Mach number

One of the famous and convenient test cases for shock wave modeling is the 1D Sod's shock tube. This is a Riemann problem for the compressible Euler equations of gas dynamics. The initial set up has ...
2
votes
1answer
77 views

What is the exact formulation of compressible Euler equation of gas dynamics in polar coordinates with artificial diffusion in 2D?

The interested equation is advection-diffusion equation. One of the canonical example is Navier-Stokes equations. However, I would like to let the coefficient of diffusion constant goes to zero, ...
0
votes
0answers
39 views

Numerical solution of hyperbolic system

Good day. I have the following system $\begin{cases} \frac{d^2A}{d z^2}+{\Delta}_\perp A - \frac{1}{c^2}\frac{d^2A}{d t^2} =\frac{4\pi}{c}j_\perp, \\ \frac{d n}{d t}+\frac{d}{d ...
0
votes
0answers
8 views

Wave Equation with Constant Boundary Conditions [migrated]

I need to find a formal solution to \begin{eqnarray} &u_{tt} &= c^2 u+{xx}, \;\;\;0<x<1, \mathrm{and \;}t>0\\ &u(x,0)&=x+1,\\ &u_t(x,0)&=x(1-x), \;\;\;\;0 \leq x \leq ...
2
votes
1answer
123 views

Energy Conservation in Conservation Laws with Source Terms

I'm wondering if anyone can help me understand energy conservation when using conservation law methods (i.e. Riemann solver, High-Resolution Wave-Propagation Methods) with the addition of source ...
5
votes
2answers
121 views

Advice on numerical solution for 2D hyperbolic PDE with zero flux boundary conditions

I would like to numerically solve a hyperbolic PDE of the form $\frac{\partial\theta_t}{\partial t}(x,y)+\frac{\partial\left[\theta_t \gamma_t^x\right]}{\partial x}(x,y)+\frac{\partial\left[\theta_t ...
0
votes
2answers
93 views

WENO reconstruction of flux involving derivative terms

I have a set of modified compressible Euler equations that I would like to solve using a WENO method. The issue is that the modified flux function involves derivative and filtering terms and I'm not ...
3
votes
1answer
125 views

Euler's equations for a tube with varying cross-section

$\def\pd{\partial}$ $\def\l{\left}\def\r{\right}$ $\def\mdot{{\dot{m}}}$ $\def\eps{\varepsilon}$ Consider a tube with longitudinal coordinate $x$ from $0$ to $l$ and varying cross-section $A(x)$. ...
6
votes
1answer
189 views

Conservative finite-difference expression for the advection equation

Following on from the earlier question I am trying to derive a finite-difference scheme for the advection equation which is conservative. It was suggested that for advection equation with variable ...
7
votes
0answers
140 views

How does positivity preservation fit into the implication chain from monotone to monotonicity preserving?

I know from "Numerical Methods for Conservation Laws" by Randall J. LeVeque that there is an implication chain of properties of methods for conservation laws: monotone $\Rightarrow$ $L^1$-contractive ...
8
votes
1answer
314 views

Higher order Lax-Wendroff type scheme?

Suppose we want to solve a hyperbolic conservation law $u_t+f(u)_x=0$. I really like to use Lax-Wendroff, which reads $u_j^{n+1} = u_j^n -\frac{\Delta t}{\Delta ...
1
vote
2answers
340 views

complexity of flux limiter techniques

My question is not related to any particular problem, rather, I am looking at the equations of the form $$u_t+c(t,x)u_x=0$$ and attempt to solve it numerically. According to ...
4
votes
1answer
190 views

mathematical statement of “open” boundary condition

For your information, the original equation comes from here. Note: You DON'T have to read the paper. I will make the question as self-contained as possible. The central equation to solve is equation ...
3
votes
1answer
456 views

amplification factor for the Crank Nicolson scheme for the advection equation

I will try one more time being more detailed and careful. Consider the transport equation of the form $$u_t+au_x=0, t\in[0,T],x\in \mathbb{R}, a>0$$ and initial condition $u(0,x)=u_0(x)$. I would ...
1
vote
1answer
170 views

amplification factor of some schemes for the transport equation [closed]

Any source of FDM schemes for the transport equation starts by explaining that explicit central differences for the equation of the type $u_t+au_x=0$ can cause oscillations and unconditionally ...
0
votes
1answer
239 views

ENO/WENO vs monotone Hermite interpolation

I have see the method PCHIP in matlab that implements the monotone Hermite interpolation method which was originally proposed by Carlson in 1980s. It seem to accomplish the goal of preventing the ...
1
vote
1answer
127 views

energy norm for transport equation

I asked this question before but did not have any luck with an answer. It might be a student level question but I need to understand that with possibly some help. I am considering the hyperbolic ...
8
votes
1answer
306 views

Nonlinear wave equation - Finite element or finite difference

I would like to know the which is more advantageous when it comes to solving nonlinear hyperbolic equations, Finite Element or Finite difference methods? Which method will be better in capturing ...
1
vote
0answers
110 views

Enforcing continuity conditions in pseudospectral domain decomposition methods for time dependent PDEs

I have a partial differential equation of the form $$ \frac{d}{dt}f(x,t) = \Theta(x) f(x,t) \qquad \Theta(x) \sim \left[\frac{d^2}{dx^2} + k^2(x)\right] $$ subject to $f(x,t=0) = f_0(x)$, and ...
6
votes
1answer
117 views

interpolation combined with methods of characteristics can cause oscillations for the transport equation?

I would like to know about the effect of using a higher order interpolator for the methods of characteristics. I am solving $$u_t+a(x,t)u_x=0$$ with some nonsmooth initial data $u_0(x)$ by the method ...
6
votes
2answers
316 views

Numerical Green functions for a nonlinear wave equation

I am trying to put down some code to get numerically the solution of the following PDE: $$ \partial^2_t\phi-\partial^2_x\phi+\lambda\phi^3=\delta(x)\delta(t). $$ Of course, there are several ...
10
votes
2answers
488 views

Finite difference scheme for “wave equation”, method of characteristics

Consider the following problem $$ W_{uv} = F $$ where the forcing term can depend on $u,v$ (see Edit 1 below for the formulation), and $W$ and its first derivatives. This is a 1+1 dimensional wave ...
7
votes
1answer
674 views

Finite difference coordinate transformation for spherical polar coordinates

I have part of a problem that is described by the momentum conservation equation: $\frac{\partial \rho}{\partial t} + \frac{1}{\sin\theta} \frac{\partial}{\partial \theta}(\rho u \sin \theta) =0$ ...
7
votes
1answer
176 views

Numerical solution of hyperbolic PDEs with nonconvex flux

In some hyperbolic PDEs the flux is nonconvex. One example is equations in magnetohydrodynamics. What are the complications in the wave structures of such problems? What general precautions one should ...
3
votes
2answers
286 views

Entropy fix for godunov scheme

For non linear system of hyperbolic PDE, The finite volume methods work well (because of inherent conservation). Godunov scheme is a very elegant solution philosophy. For linear system, it is nothing ...
15
votes
3answers
3k views

Recommendation for Finite Difference Method in Scientific Python

For a project I am working on (in hyperbolic PDEs) I would like to get some rough handle on the behavior by looking at some numerics. I am, however, not a very good programmer. Can you recommend ...
3
votes
2answers
2k views

Can't understand a simple wave equation matlab code

I'm trying to figure out how to draw a wave equation progress in a 2D graph with Matlab. I found this piece of code which effectively draw a 2D wave placing a droplet in the middle of the graph (I ...
8
votes
1answer
188 views

Adaptive mesh refinement with perfectly matched layers?

We have an adaptive mesh refinement (AMR) code for solving the elastic wave equation with frictional fault interfaces (based on Chombo for those that are interested). One of the things that we have ...
3
votes
1answer
526 views

How to obtain an implicit finite difference scheme for the wave equation?

Suppose I had the following problem: $U_{tt}=U_{xx}+U_{yy}$ in $\Omega=[0,1]\times[0,1]$ $U(x,y,0)=f(x,y)$ $U_{t}(x,y,0)=g(x,y)$ $U=0$ on $\partial \Omega$ I know that there is an explicit ...
12
votes
2answers
266 views

Which time-integration methods should we use for hyperbolic PDEs?

If we employ the Method of Lines for discretization (separate time and space discretization) of hyperbolic PDEs we obtain after spatial discretization by our favorite numerical method (fx. Finite ...
10
votes
1answer
496 views

What are possible methods to solve compressible Euler equations

I would like to write my own solver for compressible Euler equations, and most importantly I want it to work robustly in all situations. I would like it to be FE based (DG is ok). What are the ...
8
votes
2answers
301 views

Coupling FEM DG methods to Riemann solvers

Are there any good papers and or codes that couple discontinuous galerkin finite element solvers with Riemann solvers? I need to explore coupling elliptic and hyperbolic problems but most splitting ...
10
votes
1answer
169 views

How to construct well-balanced finite volume and discontinuous Galerkin methods for hyperbolic PDEs with source terms?

Source terms, such as those due to bathymetry in the shallow water equations, need to be integrated in a special way in order to preserve physical steady states. Is there a general way to construct ...
15
votes
4answers
375 views

Which methods can ensure that physical quantities remain positive throughout a PDE simulation?

Physical quantities like pressure, density, energy, temperature, and concentration should always be positive, but numerical methods sometimes compute negative values during the solution process. This ...
20
votes
2answers
677 views

Why is local conservation important when solving PDEs?

Engineers often insist on using locally conservative methods such as finite volume, conservative finite difference, or discontinuous Galerkin methods for solving PDEs. What can go wrong when using a ...
11
votes
1answer
626 views

When should implicit methods be used in the integration of hyperbolic PDEs?

Numerical methods for solving PDEs (or ODEs) fall into two broad categories: explicit and implicit methods. Implicit methods allow larger stable timesteps but require more work per step. For ...
10
votes
2answers
508 views

How can I choose a good Riemann solver when numerically solving a system of hyperbolic PDEs?

Many numerical methods for hyperbolic PDEs are based on the use of Riemann solvers. Such solvers are essential for accurately capturing shock waves. There are a range of such solvers available for ...