Questions tagged [pde]

Partial differential equations (PDEs) are equations that relate the partial derivatives of a function of more than one variable. This tag is intended for questions on modeling phenomena with PDEs, solving PDEs, and other related aspects.

Filter by
Sorted by
Tagged with
0
votes
0answers
22 views

Overflow in Solving numerically the 1D Kuramoto-Sivashinsky Equation using spectral methods

I'm having trouble with the following implementation of the KS model (see below) found on Solving numerically the 1D Kuramoto-Sivashinsky Equation using spectral methods When tN > 300 an overflow ...
0
votes
1answer
114 views

Best method to solve this system of PDEs?

I have a system of PDEs constituting an initial value problem (IVP) consisting of three coupled PDEs: \begin{align} \partial_t \rho + \partial_x(\rho v) &= \left(k_A (1-\phi) + k_B \phi \right)...
3
votes
0answers
41 views

Appropiate Artificial Boundary Conditions for the radial part of the Klein Gordon equation?

I am trying to simulate the following equation using FDTD $ \left(- \partial^2_t + \partial^2_x + V(x) \right) \psi(x,t) =0 $ subjected to the initial conditions $\psi(x,0) = f(x),~ \partial_t \psi(x,...
0
votes
1answer
71 views

1D wave equation using Finite difference method MATLAB

I have the wave equation $$u_{tt} = 4 u_{xx}$$ with the boundary conditions $$u(0,t) = u(L,t) = 0\,,\quad x \leq 0 \leq 2\pi \,,\quad t\geq 0$$ and initial conditions $$\begin{align} &u(x,0)=\...
6
votes
0answers
107 views

FEM : energy minimization VS PDE solving

Engineering FEM When I studied engineering, I learned the traditional approach for finite elements for elasticity. The point was to solve the PDE $-div(\sigma)=f$ as: Multiply your PDE with a test ...
3
votes
1answer
109 views

Wave equation, wave is bouncing off Neuman boundary

Wave equation. Mixed BC. Applied Neuman boundary condition ($\frac{\partial u}{\partial x}\big|_N=0$) to the RHS of the domain You may observe the sharp edge in the middle of the string in the image ...
2
votes
1answer
63 views

Finite element modelling of thermal expansion of 3D solid bodies

I want to solve the thermal expansion of solid by using FEM approach. When I developed the model based on the the principle the minimum potential energy, the solutions for thermal expansion are not ...
1
vote
1answer
141 views

Non-reflective boundary condition

I'm currently solving incompressible Navier-Stokes system of equations with periodic flow and high viscosity. Is there any outlet boundary types that avoids the reflection of flow from the outlet back ...
2
votes
2answers
298 views

Semi-infinite domain transformation

Question is mostly related to literature or suggestions. Given a semi infinite domain: $x=[0; +\infty);y=[0; +\infty)$. Willing to transform it to computational domain of: $[0,1]\times[0,1]$. I did ...
1
vote
2answers
89 views

What is known about C0 triangular finite elements with nonstandard mesh point placement?

I'm curious about the general case, but for ease of explaining lets just take the case of a $P^2(\Omega)$ approximation. For simplicity, let's also just consider the reference element $(0,0), (0,1), (...
3
votes
1answer
79 views

"Optimal" domain partitioning in domain decomposition algorithms

When solving a PDE numerically by domain decomposition methods, what is the "optimal way" to split the domain? Are there any results stating that a particular partition of the domain yields &...
2
votes
1answer
185 views

Solving stiff ODEs: Dealing with Jacobian terms which take too long to compute with finite differences

I have a system of PDEs describing atmospheric chemistry and transport. I use finite-differences to make my system of PDEs into a system of ~10,000 ODEs. I then integrate the ODEs forward in time with ...
3
votes
1answer
90 views

Stability analysis simplification for PDE

I have the nonlinear PDE $$\frac{\partial U(z,t)}{\partial t} + A(U)\frac{\partial U(z,t)}{\partial z} + B(U)U(z,t) + C(z,t) = 0,$$ where $A(U)$ and $B(U)$ are guaranteed to be real and positive. I ...
0
votes
0answers
49 views

What is the most common loss function used with collocation methods for differential equations

I was looking at the Cheney and Kincaid book (6th edition) on numerical methods, with respect to collocation method for differential equations. Now for linear systems of ODES, collocation is just a ...
0
votes
1answer
27 views

How to decrease error in (FTCS) forward time centered space method?

I am using the FTCS method for solving differential equations. I know that the condition for stable output is $$ \frac{\alpha \Delta t}{\Delta x ^2} < \frac{1}{2} $$ But when I use the distance ...
1
vote
0answers
50 views

How to numerically solve PDE that governs the free vortex wake model?

Crossposted at Math SE I am reading a paper on the free vortex wake model for a helicopter rotor blade, which is described by the following PDE: $$\frac{\partial \vec{r}}{\partial \psi} (\psi, \zeta) ...
0
votes
0answers
44 views

introduction to non-overlapping domain decomposition in 1D

I am new to domain decomposition and have been searching in google for an introduction in 1D which goes over the complete procedure from the continuous to the discrete problem as well as the algorithm ...
0
votes
0answers
53 views

Discretizing Multi-species Ion Exchange Equations by Finite Volume Method

I'm solving a system of multispecies ion exchange equations (diffusion+drift fluxes) in 1-d spherical domain using finite volume method to obtain the ion concentrations at the next time step. After ...
1
vote
2answers
181 views

Continuation of solution to $-\nabla\cdot (k(x,y)\nabla u)=f$

I'm trying to solve the following problem, I had previously opened another discussion for the implementation and well, it seems that it has turned out well, it can be found here. I need to calculate ...
3
votes
0answers
82 views

Strange Picard iteration

I am interested in solving the equation $$ \begin{aligned} \nabla \cdot\left(\nabla \phi-\frac{\nabla \phi}{|\nabla \phi|}\right) &=0 & & \text { in } \Omega \\ \phi &=0 & & \...
3
votes
0answers
75 views

Solving PDEs: What is the best way to deal with non-banded/dense jacobians?

I have a system of PDEs describing atmospheric chemistry and transport. I use finite-differences to make my system of PDEs into a system of ~10,000 ODEs. I then integrate the ODEs forward in time with ...
0
votes
2answers
356 views

Implementing routine for $-\nabla\cdot (k(x,y) \nabla u)=f$ in Matlab

I am solving the Poisson Equation for 2D given by the following expression: $$-\nabla\cdot (k(x,y) \nabla u)=f$$ in a rectangle with Dirichlet conditions on the boundary using Matlab. In principle I ...
0
votes
1answer
138 views

Help me choose a book on the numerical integration of PDEs

For ODEs I have these books: Griffiths, David, Higham, Desmond J., Numerical Methods for Ordinary Differential Equations, 2010 Alfio Quarteroni, Riccardo Sacco, Fausto Saleri, Numerical Mathematics, ...
5
votes
1answer
256 views

What is the weak form of a vector type Laplace equation?

For a scalar variable $u$, the weak form of the following Laplacian: $$\nabla^2 u =0 $$ with the assumption that $v$ vanishes at the boundary is $$\int \nabla v . \nabla u \, d\Omega = 0 $$ which is ...
0
votes
2answers
85 views

Stability of the Forward-Time Central Space method, section 9.3 in LeVeque

I am reading section 9.3 in Leveque about stability. The explanations are very short and brief so I am asking for some help to understand and elaborations. The result of the section is that for FTCS ...
2
votes
1answer
143 views

How to discretize a non-linear PDE with boundary conditions and intial value

Consider this non linear PDE: $$u_t + c(u^2)_x =\alpha u_{xx} \text{ , } -1<x<1 , t>0 $$ with $$u(-1,t) = g_L(t) \text{ , } u(1,t) = g_R(t) \text{ and } u(x,0)=f(x) $$ where the 3 functions(...
2
votes
1answer
56 views

How to demonstrate the order of convergence of FTBS method for solving a hyperbolic PDE

consider the Purely hyperbolic model problem $$u_t+au_x=0$$ $$u(-1,t)=u(1,t) \text{ (periodic boundary)}$$ $$u(x,0)=f(x)$$ with $f(y)=\sin(2\pi y)$. Furthermore the exact solution is given by $u(x,t)=...
2
votes
3answers
183 views

Solving a PDE with 2 variables, with one variable whose derivative with respect to space is only known

I am trying to solve a PDE of the form $$\frac{\partial u}{\partial t} = D \frac{\partial^2u}{\partial x^2} + k\ \ \ (1)$$ where only $k$'s first derivative with respect to $x$ is known $$\frac{\...
1
vote
0answers
51 views

Stability plot of upward difference implicit time

I am analyzing the stability of 1D convection (advection) equation as shown in the picture. When I derive the equations as shown I want to get rid of the complex number. I`m asking if those stability ...
0
votes
0answers
122 views

Error $L_{2}$ convergence in Finite Element for Poisson Equation

I have written a Matlab code to solve the equation $-u'' = f$ with conditions $u(0) = u'(1) = 0$ on the domain $x \in [0,1]$. I tested the code with $f(x) = -2, \forall x\in [0,1]$. I check the plot ...
3
votes
1answer
156 views

Finite element method for high-frequency electromagnetics

I am writing a project about the Finite element method for use in high-frequency solutions of Maxwell's equations. This could be for use in antenna design and similar. I have some trouble ...
1
vote
0answers
42 views

discretization of advection diffusion with variable coefficients

I am looking for help to find a somewhat stable FD numerical scheme for the advection diffusion equation posed on a curve $(x(r),y(r))$. The equation becomes $$u_t=\alpha(r) u_r +\beta(r) u_{rr}+f(r,t)...
0
votes
1answer
81 views

High Running Time and Suboptimal Accuracy of 2D Wave Equation Solver with Finite Differences

Im trying to solve the following 2D wave equation: $$u_{tt} = u_{xx} + u_{yy}, \hspace{3mm} u(x,y,0) = \cos(4 \pi x) \sin(4 \pi y), \hspace{3mm} u_t(x,y,0) = 0$$ with the periodic boundary condition ...
0
votes
1answer
172 views

Trouble Implementing 1d Wave Equation Finite Difference Solver

Im trying to solve the 1d Wave Equation on $x \in \mathbb{R}, t > 0$: $$u_{tt} = c^2u_{xx}, \hspace{5mm} u(x,0) = \cos(4 \pi x), \hspace{5mm} u_t(x,0) = 0$$ with $c = 1$ and a periodic boundary ...
2
votes
0answers
67 views

Divergence on wave equation simulation

I'm currenly working on my own PDE solver for non-linear simulations in python. I've done succesfully simulations for KdV and Fisher's equation, but now I'm playing with second order derivatives in ...
1
vote
0answers
60 views

Solving Laplace equation with constraint on boundary

I have found the following PDE problem in a paper: Essentially, we have a rectangular domain where there is a unknown interface $z=\xi(x)$ (liquid-air interface) separating the domain into two medium ...
0
votes
0answers
48 views

Cauchy problem ill-posed?

Find the solution to the Cauchy problem consisting of the wave equation : $$u_{xx}-u_{yy}=0$$ together with initial conditions: $$ u(x,0)=0,$$ $$u_{y}(x,0)=g(x)$$ for some known initial datum $g$. Is ...
0
votes
0answers
42 views

Output solution vector in Petsc

I am using petsc to solve a linear elasticity problem discretized by finite elements.The initial mesh is read by a mesh file and the distribution in each processor is done using METIS.I am using only ...
0
votes
0answers
54 views

When to stop iterations in SOR solver for 3D Poisson equation

I'm writing a solver (in C) for 3D incompressible fluids, using the finite-differences method, and I'm finding a somewhat surprising behaviour: the solver provides "good-looking" solutions, ...
1
vote
0answers
59 views

Integrating a wavelike equation with absorbing boundary conditions

I am trying to numerically solve the following equation: $\frac{\partial^{2} \phi}{\partial t^{2}}-\frac{\partial^{2} \phi}{\partial x^{2}}+V(x) \phi(x, t)=0$ On some domain, with: $\phi(x, 0) = I(x)$ ...
0
votes
1answer
81 views

Acoustic Simulation, how are boundaries handled?

I don't have a background in numerical modeling so this question is rather broad. What I am interested in is modeling the propagation of an ultrasonic acoustic wave in 3d space. The basic 3d wave ...
0
votes
0answers
38 views

Devising Convergent Numerical Scheme for PDE

I'm currently looking at the PDE \begin{align*} u_t + \left[x(1-y) - (1-x)\right]u_x - (1-y) u_y + (z-xy) u_z = (z-xy) u_{xy} - (1-x)u& \\ \end{align*} with \begin{align*} u(x,y,z,0) = 1& \\ ...
4
votes
1answer
749 views

FEM Python book

Is there any book or site available with Finite element Method for partial differential equations with python code apart from Fenics?
0
votes
1answer
65 views

why do i receive this error?

I want to plot the 100,200 and 400 iterations of this function of non homogeneous parabolic pde ...
0
votes
1answer
139 views

Why this error occurs in my code for Lax Wendroff?

I want to implement the Lax Wendroff method for a non linear advection equation which is $$\frac{u_{i}^{n+1}-u_{i}^{n}}{t} + \frac{f(u_{i+1}^{n})-f(u_{i-1}^{n}) }{2h} -\frac{t}{2h} \left( F_{i+1/2}^{n}...
2
votes
1answer
115 views

Developing a meshfree contouring algorithm

I would like to find the contours of a scalar function $u(x,y)$ available as a discrete set of function values $u_i = u(x_i,y_i)$ over a scattered set of points $(x_i,y_i), i=1,...,N$. In my case, the ...
1
vote
1answer
102 views

Which 2D PDE with an exact solution can I use to test/verify my FEM-PDE code?

I have created a program to solve 2D, time-dependent PDEs with the finite element method and get reasonable looking results for the 2D acoustic wave equation. Now I would like to go further and solve ...
0
votes
1answer
77 views

Are spatial boundary conditions required for PDEs discretized with Method of Lines?

As far as I understand, you need to define boundary conditions in time and space to select a unique solution to a PDE and make it solvable. However, in ODEs I only need to specific the initial value ...
1
vote
1answer
110 views

What is wrong in the code for this upwind method?

I want to implement the upwind method in following advection equation problem : $$ u_{t}+2 u_{x} =0 ,$$ for $0\leq t \leq 1,$ $0\leq x \leq 1 $ $$ u(0,x) = u_{0}(x) = \begin{cases} 10^4 (0.1-...
1
vote
1answer
90 views

Finite Difference for Advection Equation With Source

I'm trying to find a convergent finite difference scheme for the PDE \begin{equation} \begin{split} u_t + (x-1) u_x &= (x-1)u, \hspace{.5cm} x \in [0,1] \\ u(x,0) &= 1 \\ u(1,t) &= 1. \\ \...

1
2 3 4 5
16