Questions tagged [pde]

Use this tag for questions about partial differential equations. For questions about the Eclipse Plugin Development Environment, use eclipse-pde instead.

Partial Differential Equations are differential equations that contain unknown multivariable functions and their partial derivatives. See [http://en.wikipedia.org/wiki/Partial_differential_equation][1]

The Eclipse Plug-in Development Environment (PDE) provides tools to create, develop, test, debug, build and deploy Eclipse plug-ins, fragments, features, update sites and RCP products. PDE also provides comprehensive OSGi tooling. See [http://www.eclipse.org/pde/][2]

Questions about the Eclipse Plug-in Development Environment are better tagged as 'eclipse-pde'.

See

Partial_differential_equation

Eclipse Plug-in Development Environment

477 questions
2
votes
1 answer

Examples for using fipy.steppers.pidStepper

I'm using FiPy to simulate a moderately large system of coupled PDEs (11 equations) that calculate the time-dependent evolution of a biochemical advection-diffusion-reaction network. Currently I'm using a fixed time step, and it works fine. However,…
jmrohwer
  • 21
  • 2
2
votes
1 answer

Fixed flux boundary conditions in FiPy

I have a follow-up question for this thread Coupled non-linear equations in FyPi. I managed to set up the system and get the reasonable results when using the Neumann boundary conditions for all variables (i.e. keeping the constant concentrations of…
josefT
  • 21
  • 2
2
votes
1 answer

Non-linear HJB PDE using Fipy, PDE translation

I'm trying to solve a non-linear PDE HJB equation using FiPy, but i have some difficulties translating the PDE into the proper FiPy syntax: I tried something like : eqX = TransientTerm() == -DiffusionTerm(coeff=1) + (phi.faceGrad *…
jessi_123
  • 21
  • 1
2
votes
1 answer

Vectorising / Paralellising GEKKO equations

Using Phython's GEKKO library and following the example to solve the Wave Equation from their original paper, one needs to create a whole lot of equations. m.Equations([u[i].dt() == v[i]for i in range(npx)]) m.Equations([v[i].dt() ==…
Swift
  • 360
  • 1
  • 3
  • 12
2
votes
1 answer

GEKKO optimisation gets stuck (APM.exe)

I am trying to run a differential equation system solver with GEKKO in Python 3 (via Jupyter notebook). For bad initial parameters it of course immediately halts and says solution not found. Or immediately concludes with no changes done to the…
Swift
  • 360
  • 1
  • 3
  • 12
2
votes
0 answers

How to import non-Newtonian model in FEATOOL Multiphysics?

I am using FEATool Multiphysics. I have a MATLAB function named ex_nonnewtonian1.m, which is supported by FEATool itself and contains the definition of a power-law non-Newtonial model. This file wasn't added to FEATOOL, so I need to add it…
2
votes
0 answers

Modelling a chemical engineering diffusion problem in 1D with FiPy

I want to simulate an adsorbent bed for a chemical reactor by solving a set of partial differential equations. I want to reduce eq_1 and eq_2 to a set of ordinary differential equations via finite volume techniques for spatial terms. Ct , DL and eb…
GOKAY AVCI
  • 21
  • 1
2
votes
1 answer

How to initialize 2D-discontinuous initial data in pyhthon

I've just coded my finite difference solver in Python for the heat equation in the unit square, and in order to check an important property of the equation I need to give as initial data a discountinuous function. My domain in [0,1] x [0,1] and I'd…
slamWolfen
  • 43
  • 4
2
votes
2 answers

add resources to eclipse project / plugin

I'm new to eclipse plugin development. I've been able to create a new project type and add folders to it. What is the best way to add a starting set of files into this project? Should I have my JAVA class create these files? Or can I package them…
majestiq
  • 545
  • 8
  • 25
2
votes
1 answer

Complex PDE (Ginzburg Landau) in Julia with Pseudo-Spectral method

I want to teach myself about solving PDEs with Julia and I am trying to solve the complex Ginzburg Landau equation (CGLE) with a pseudospectral method in Julia now. However, I struggle with it and I am a bit of ideas what to try. The CGLE reads:…
2
votes
2 answers

Emacs Pretty Symbols in cperl mode: :: must be surrounded my spaces

I'm using Emacs with PDE and cperl-mode. I really want prettify-symbols-mode to work with it. When just using perl-mode this happens: Mod::thing->new( {c => 'sea'} ); becomes Mod∷thing→new( {c ⇒ 'sea'} ); however when using cperl same expression…
CharJe
  • 41
  • 6
2
votes
2 answers

Assigning to multiple sub-matrices of a matrix simultaneously. Possible optimization by vectorized indices

Is there a clever way to vectorize a for-loop which assigns elements to submatrices of a matrix? Initially, I had two for-loops: U=zeros(6*(M-2),M-2); for k=2:M-3 i=(k-1)*6+1; for j=2:M-3 …
stein
  • 23
  • 4
2
votes
0 answers

Solving a system of PDEs using implicit methods

I've got a system of partial differential equations (PDEs), specifically the diffusion-advection-reaction-equation applied to heat transfer and convection, which I solve using finite difference method. In my simulation environment I've got a…
JE_Muc
  • 5,403
  • 2
  • 26
  • 41
2
votes
2 answers

Is it possible to solve a system of PDEs using sympy in Python?

Assume that you have a system of coupled PDEs, such as 1st PDE in F(a,b) 2nd PDE in F(a,b) Following code is able to solve each one of the PDEs separately: import numpy as np import sympy as sp # definition of variables a, b = sp.symbols('a b') f…
Karlo
  • 1,630
  • 4
  • 33
  • 57
2
votes
0 answers

Numerically solving a pde with changing domain in Matlab

I'm trying to implement a changing domain into my Matlab code which solves a pde numerically. In the given problem, I am modelling the baking of bread. It is a 1D problem where there is a single heat source from above. I am using the heat diffusion…
Dedavond
  • 21
  • 2