Questions tagged [differential-equations]
The differential-equations tag has no usage guidance.
100
questions
2
votes
2answers
52 views
Methods for solving discrete PDEs using algorithmic differentiation results
I'm looking for a method to solve a 20000 variable, 20000 residual non-linear PDE with a Galerkin method.
I have Fortran subroutines for:
The residuals: $\vec{r}(\vec{x})$;
Their Jacobian multiplied ...
0
votes
0answers
39 views
Solver for large dense BVP system in python
I have a large system of boundary value problems of the form
$$ \frac{d^2 y }{dt^2} = C(t) y + b(t), $$
where the variable $y$ is a vector that has anywhere from 50 to around 500 components, $C$ is a ...
2
votes
2answers
87 views
Finite difference method having a discontinuity
I am trying to understand the FDM which is a widely used method solving differential equations by using approximation below.
$$\dfrac{\partial u}{\partial x}=\dfrac{u(i+1)-u(i-1)}{2\Delta x}$$
How can ...
0
votes
1answer
86 views
Two RK4 method in one program
I want to solve this integral using RK4 by coding in Fortran:
$$R=∫1/a(t) dt → dR/dt=1/a(t) =f(t)$$
Initial point: t=0 (or a=0.001) and R=0
And I have to get a(t) by solving another ...
-1
votes
0answers
47 views
Analytical vs Numerical solution to differential equation with final value
I have the following differential equaiotn: $x' = -x-2$ with a final value x(2)=0. I solved it analytically and got: $x(t)=2(e^{2-t} -1)$. It satisfies the final condition x(2)=0. Also, at t=0, x(0)=...
0
votes
0answers
26 views
Simulating the response of nonlinear system with stiff differential equations
I want to simulate the response of a nonlinear system given in the following form:
$$ \dot{x_1} = f_1(\bar{x_1})+g_1(\bar{x_1})x_2, \ x_1(0) = 0.2 $$
$$ \dot{x_2} = f_2(\bar{x_2})+g_2(\bar{x_2})x_3, \...
1
vote
1answer
93 views
Partial differential equation FEM application
I have a PDE which looks like Helmholtz wave equation on one dimensional domain.
$$\dfrac{d^2u(x)}{dx^2}+\pi^2u(x)=f(x)$$
where $-\infty <x<\infty $
Also, $f(x)= 1$ for $-0.25<x<0.25$, I ...
1
vote
1answer
56 views
How to write a code of 2D ADI method in matlab?
I tried to write a code for the alternating direction implicit (ADI) method in 2D, but I got stuck.
My equation is:
$$\frac{\partial U(t,x,y)}{\partial t} = 2\Delta U(t,x,y) -10(\frac{\partial U(t,x,y)...
0
votes
1answer
33 views
MATLAB ode45 doesn't start at initial conditions
I wrote a code in MATLAB to solve a system of differential equations, but my solution doesn't seem to take into consideration the initial conditions I specified. I am not sure how to interpret this ...
1
vote
3answers
123 views
Runge-Kutta method for an ODE with initial value which is root of denominator
I wrote a code in Fortran to solve this differential equation using RK4 method:
$$
\frac{dy}{dx}=A\sqrt{\frac{B}{y}+\frac{C}{y^2}}
$$
$A$, $B$, and $C$ are some known constants. The problem is that ...
0
votes
0answers
61 views
Applying boundary Conditions on FEM
I have a partial differential equatons as shown below.
$$\dfrac{d}{dx}((1+x)\dfrac{du(x)}{dx})=0$$
With the following boundary conditions.
$$u(0)=0, u(3)=10$$
To solve it using FEM, I multiplied the ...
0
votes
0answers
11 views
Differential parameterized inequalities
Let $H$ be an Hamiltonian and denote $\vec{H}$ the associated Hamiltonian vector field.
I am interested in solving numerically the following problem
$$
\dot z(t) = \vec{H}(z(s),t_1,\ldots, t_p) \...
2
votes
0answers
70 views
Numerical method for harmonic oscillator with jumping constant
Let $k_1 \neq k_2$ be positive reals, $t_0 > 0$ and consider the following Cauchy problem in $[0,+\infty)$:
\begin{cases}
y(t) + k(t)y''(t) = 0 \newline
y(0) = 1/\sqrt{k_1} \newline
y'(0) = 0,
\end{...
2
votes
1answer
128 views
Solving numerically a linear ODE
I start by saying that I do not have a strong background in numerical analysis, so I may miss some basic things or make trivial mistakes.
Motivated by some problems in digital signal processing, I ...
2
votes
1answer
114 views
Lambdifying a symbolic matrix in Julia
If I have a symbolic matrix defined as T below, is there any way to lambdify this as function of variables, say σ..., and return ...
2
votes
0answers
47 views
Modelling of Stefan Maxwell equation
I am trying to solve Maxwell Stefan's equation over a membrane to get the transient mole fraction distribution over the membrane thickness 'z'. But somehow I am not able to code it using ODE45, more ...
0
votes
0answers
50 views
How to solve odd-order differential equations in FEM? Petrov-Galerkin?
I've recently learned about using weighted residuals with the Galerkin method to numerically approximate even-order differential equations (for linear elements, I'm still a beginner). It seems for odd-...
1
vote
0answers
20 views
Optimality conditions for optimal control: BVP - DAE
I am solving an optimal control problem of the form
$$
\min_u \qquad\int_0^T \langle u(t), u(t) \rangle \, \mathrm{d}t \\
s.t. \quad \dot{x} = \tilde{f}(x) + u, \quad x(0)=x_0 \\
\qquad \tilde{\Phi}(...
4
votes
2answers
65 views
MInimizing cost function using iterative search for a minimum method
I want to estimated the parameters $\ \hat{\theta} $ of a model using an iterative search for the minimum of a cost function. The cost function is defined as follows:
$$ V_N(\hat{\theta}) = \frac{1}{...
1
vote
2answers
102 views
Solving a parameter estimation problem using trajectory optimization
This is a follow-up to my previous question here
I've the following system of equations for studying information flow in the below graph,
$$ \frac{d \phi}{dt} = -M^TDM\phi + \text{noise ...
2
votes
1answer
80 views
Using Implicit Euler with second order differential equations
We can numerically integrate first order differential equations using Euler method like this:
$$y_{n+1} = y_n + hf(t_n, y_n)$$
And with Implicit Euler like this:
$$y_{n+1} = y_n + hf(t_{n+1},y _{n+...
1
vote
1answer
42 views
How to properly compute weights for Weighted Least Squares (WLS)?
I want to apply the weighted least squares method in order to identify parameters of a dynamic process. The process is described by a second order differential equation of the form:
$$ \ddot{y}+a_1\...
2
votes
0answers
52 views
Numerical errors due to terms of the form $\frac{1}{r}$ (r goes to 0 at the boundary) while using finite difference method
I am trying to solve a system of differential equations using finite difference method.
There are few terms of the form $\frac{A(r)}{r}$, both $A(r)$ and r go to zero at the boundary. Analytically ...
1
vote
0answers
63 views
Are linear, CTCS codes always stable?
I would like to solve some equations which basically look like this
$$\frac{\partial u}{\partial x}=F\left(v,\frac{\partial v}{\partial y},\frac{\partial^2 v}{\partial y^2}\right),$$
$$\frac{\partial ...
5
votes
2answers
823 views
Solving coupled differential equations in Python, 2nd order
I have a system of coupled differential equations, one of which is second-order. I am looking for a way to solve them in Python. I would be extremely grateful for any advice on how can I do that!
$k$...
1
vote
1answer
52 views
Dealing with boundary conditions using Fourier spectral methods
I am currently working on a project where I need to use Fourier spectral methods to solve the KS equation. I found this code which is using the Fourier spectral methods to solve the classic 1D heat ...
0
votes
0answers
272 views
Numerically solving a partial differential equation in python with Runge Kutta 4
I'm supposed to solve the following partial differential equation in python using Runge-Kutta 4 method in time.
$$
\frac{\partial}{\partial t}v(y,t)=Lv(t,y)
$$
where $L$ is the following linear ...
2
votes
1answer
294 views
Numerical Solution to Rayleigh Plesset Equation in Python
I have been trying to numerically solve the Rayleigh-Plesset equation for a sonoluminescence bubble in Python. You can read about this phenomenon here: https://iopscience.iop.org/article/10.1088/0143-...
1
vote
0answers
54 views
Stably solve transport equation with source term
I am trying to solve a transport equation of the form for the variable $\psi(t,r)$
\begin{equation}
\partial_t\psi-\alpha(r)\partial_r\psi-\beta(r)^2\psi-f(t,r)=0
,
\end{equation}
where I am solving ...
1
vote
0answers
20 views
Using nondimensionalization to solve an ode in MATLAB [duplicate]
I am trying to solve an ode that uses some extremely large numbers and some extremely small numbers, namely
$$
e = 1.6\times 10^{-19}\\
E = 10^6\\
\tau = 6\times 10^{-24}\\
m = 9.1\times 10^{-31}\\
c ...
2
votes
1answer
92 views
MATLAB's ode45 not dealing with initial conditions well [RESOLVED]
*Concern highlighted in yellow
*Solution at bottom
I have a differential equation to solve for the motion of an electron:
$$
\frac{d^2v}{dt^2} = \frac{1}{\gamma^6}\left( \frac{eE}{\tau m} - \left( \...
1
vote
1answer
97 views
Crank-Nicholson for diffusion-advection vs diffusion equation
Let's consider the following 1D diffusion equation:
$\frac{\partial u}{\partial t} = xk \frac{\partial}{\partial x}(\frac{1}{x}\frac{\partial u}{\partial x})$
where we assume that the diffusion ...
1
vote
1answer
117 views
Numerical integration in 2D
I would like to solve the following problem
$$
\vec{v}(x,y)= k\, \nabla \theta(x,y)
$$
with respect to the unknown function $\theta$. Parameter $k$ is just a real constant quantity.
I have two ...
3
votes
1answer
163 views
Solving an SDE with time-dependent parameter in R
I am trying to solve a system of SDEs in R using the Diffeqr package.
Let's reduce the system to a simple ODE:
...
1
vote
1answer
84 views
How to Break Coupled ODEs down to first order for Runge-Kutta
My question might seem a bit simple. I am trying to solve a system of ODEs using Runge-Kutta method. I am having difficulty breaking down the equations into a system of first order ones required ...
1
vote
1answer
77 views
Numerical methods for non-linear diffusion
I have the following non-linear diffusion equation, for $\ z(x,t)$:
$\ z_t = -C(\sin(\omega t))^m x^{hm}(hm x^{-1}(z_x)^n + n z_{xx} (z_x)^{n-1}) $
Any advice for numerical (or analytical) solutions?...
2
votes
0answers
32 views
Book Recommendation: Analysis and design of mechanistic models - such as pharmacokinetics or hydrology models
I have been looking at an interesting book "Pharmacokinetic-Pharmacodynamic Modeling and Simulation" by Peter Bonate on pharmacokinetic models: the models of how medical drugs work their way through ...
2
votes
0answers
64 views
Numerical integration of SDE: choice of $dt$ and algorithm
I am working on the following Stochastic Differential Equation (SDE) in the Quantum Mechanics context:
$$dX_{t} = a X_{t} dt + b X_{t} dW$$
where $X_{t}$ is my stochastic varible, $dt$ is my ...
0
votes
0answers
70 views
How one can simulate a system given by differential equation?
I want to simulate a diffusion environment given by the differential equation
$$\frac{\partial u(x,y,t)}{\partial t}=D\left(\frac{\partial^2 u(x,y,t)}{\partial x^2}+\frac{\partial^2 u(x,y,t)}{\...
1
vote
0answers
55 views
Governing equations vs Transport equations
This is a basic question. But I did not find any explanations for this. How are governing equations, like mass, momentum, energy conservations equations, different from 'Transport equation'?. Is a ...
0
votes
0answers
57 views
What are the differences between these different forms of equation?
What are the differences between Conservative differential form, Non-conservative differential form, Conservative Integral form and Non-conservative integral form of differential equations? I know ...
1
vote
0answers
81 views
Solved : Damped spring-mass system, wrong position, correct speed and acceleration
I am modulating a spring-mass system with gravitation and aero drag, with python programming. The spring is hanging vertically and attached a weight. The user then selects a length to drag it down ...
4
votes
1answer
139 views
Fast and free server for computing
I have to calculate a huge differential equation. With my laptop, it's going to be computed for several days.
Is there a free (I need just for 3 days) fast server for scientific calculations?
My ...
2
votes
1answer
71 views
Symplectic linear multistep method?
I'm doing a gravitational n-body simulator and I'm thinking of implementing linear multistep methods like Adam-Bashforth. But is there any symplectic multistep methods?
1
vote
1answer
83 views
Online Parameter Estimation using steepest descent
I have a first order system which is described by the following differential equation:
dx/dt = -a*x + b*u where u is the input <...
3
votes
2answers
112 views
Finite difference for a highly nonlinear equation - The wind within the forest
Based on the Navier-Stokes equations and a few parameterizations, the horizontal steady-state wind $u(z)$ within a forest of height H satisfies:
$$
a\left(\frac{du}{dz}\right)^2 + b\frac{du}{dz} \...
3
votes
2answers
124 views
Damped Harmonic Oscillation. Efficient algorithm to find the parameters resulting in threshold oscillation amplitude
Let's assume, that we have damped harmonic oscillation of a body in the form of a cone, immersed in a liquid.
Equilibrium condition of the body is:
$$m\overrightarrow{a} = \overrightarrow{F_\text{...
1
vote
0answers
44 views
Why does the correlation function of this stochastic differential equation starts at different points?
I am working with the following differential equation:
The equation is $$x=\beta +\sqrt{2D} \xi(t)$$ where $\xi(t)$ is a white noise term, with a reflecting wall boundary conditions.
After solving ...
2
votes
0answers
87 views
How to implement adaptive step size Runge-Kutta Cash-Karp?
Trying to implement an adaptive step size Runge-Kutta Cash-Karp but failing with this error:
...
1
vote
1answer
387 views
Solving differential equation in Python with variable coefficients (I just know the coefficients numerically)
I am trying to implement a routine to solve a differential equation in Python. Basically the kind of equation that I am interested in solving is of the form:
$\displaystyle \frac{d}{dx^2} \left(x y(x)...
