The tag has no wiki summary.

learn more… | top users | synonyms

0
votes
0answers
12 views

How to solve a nonlinear diffusion equation-comsol?

Consider a thin film with a perpendicular applied magnetic field ${H_{a}}$ (A/m) in z-axis. In fact, By increasing ${H_{a}}$ the magnetic field penetrates gradually the film. The equation for the time ...
0
votes
0answers
46 views

How to solve a nonlinear diffusion equation?

Consider a thin film with a perpendicular applied magnetic field Ha in z-axis. The nonlocal relation between Ha, the self-field Hself (generated by the eddy current J) and the local magnetic flux ...
4
votes
1answer
44 views

Eikonal Equation solver with different grid densities

The Fast Marching Method, Fast Iterative Method, and Fast Sweeping Method are three ways of solving the Eikonal Equation on a discrete grid, essentially just a wavefront spreading out from initial ...
3
votes
0answers
89 views

numerical analysis of a partial integro-differential equation

I have to numerically solve a nonlinear partial integro-differential equation. This is my equation, $$\frac{\partial y(x,t)}{\partial t}=\int_{-\infty}^\infty K_0(|x-u|) \frac{\partial^2 ...
1
vote
1answer
61 views

OpenMP threaded nonlinear solver for complex numbers

Problem: I have translated Jacobian-Free Newton-Krylov solver written by C. T. Kelley to Fortran and now want to parallelize it on shared-memory system with OpenMP. In addition, I want to precondition ...
6
votes
2answers
133 views

Strong coupling of a non-linear multiphysic problem: failure with Newton Raphson method

I am trying to solve a multiphysic problem using finite elements and a Newton Raphson solution scheme. I have two non-linear subsystems that are coupled bi-directionally. The first subsystem includes ...
1
vote
1answer
41 views

solving for unknown inside an expectation

I have a an equation that need to find its root. The function is the following $f(\theta) \equiv E[R(\theta;\eta)]=0$ for some unknown $\theta$ which is deterministic, while the expectation is taken ...
3
votes
0answers
52 views

Finding self-kissing points on a plane curve?

I have a curve in the complex plane given by $$ f(t) = \sum_k r_k\exp(2\pi\mathrm{i}(t+\varphi_k)p_k). $$ Some of the parameters are specially chosen: $r_k>0$, $\sum_k r_k=1$, $p_k\in\mathbb{Z}$, ...
3
votes
2answers
81 views

Bessel EVP and fzero

I am trying to solve the Eigenvalue problem $$ x^2 y''+ x y' + x^2 y = \lambda^2 y\,,\quad x\in(0,1)\,,\quad y(0)=0\,,\quad y'(1)=y(1) $$ The differential equation is the Bessel equation. The solution ...
1
vote
1answer
61 views

methods for a peculiar BVP system

Consider the following system defined on the open interval (-1, 1): $y_1' = c y_3 \\ y_2' = c y_4 \\ y_3' = -f(y_1, y_2)y_2 \\ y_4' = f(y_1, y_2)y_1 $ given $ y_3(-1) = 0 = y_3(1) \\ ...
0
votes
1answer
101 views

Implement Robin Boundary Condition

This is a follow up to this question. I have a nonlinear BVP on $x=0$ to $L$: $$ (T^2)\frac{\partial^2 T}{\partial x^2} + T \left( \frac{\partial T}{\partial x}\right)^2 + Q = 0 $$ to which I apply a ...
4
votes
2answers
134 views

Large-scale nonlinear optimization problem

I want to solve a nonlinear optimization problem of the following form \begin{equation} min( \sum_i d^{x_i}c_{i})\\ 0 \leq x_{i} \leq a\\ \sum_{i} x_{i} \leq b \end{equation} $a$, $b$, $0.95 < d ...
2
votes
1answer
188 views

Solving a Nonlinear BVP using Finite Difference Method

I am trying to write a code to solve a nonlinear BVP using the Finite Difference Method. The BVP is: $(T^2)\frac{\partial^2 T}{\partial x^2} + T \left( \frac{\partial T}{\partial x}\right)^2 + Q = 0$ ...
1
vote
1answer
341 views

Nonlinear ODE to solve Duffing's equation

I am trying to solve the Duffing's equation in MATLAB. $ m\ddot{y}+c\dot{y}+ky+k_{3}y^{3} = f(t) $ where $ f(t) = A \sin{\omega t}$ To do that I wrote a function to be given to the ode45. ...
0
votes
0answers
31 views

Solving a nonlinear equation with a Markov process and RVs

Assume that we have the following equation and the following assumption. The scope is to solve for some particular variables expressed later. Update $$E_{t}\left[ b(A_{t+1})^{1-\gamma} ...
0
votes
0answers
47 views

Solvers for nonlinear parabolic PDEs

Could you please advise some programs or libraries for solving parabolic PDEs (or its systems) in 1D, 2D and 3D, for example, with the method of lines? The system of parabolic PDEs can be nonlinear in ...
2
votes
1answer
136 views

Solving a nonlinear equation with random variable

I would like to solve an equation that looks like this UPDATE $E[(R^{1-\gamma})(r_k+\theta-r_z)]=0$ , where $R=\phi r_z+(1-\phi)(r_k+\theta)$ and $\phi\in[0,1]$, $\theta$, is a random variable ...
2
votes
1answer
135 views

Solving large, non-linear systems of ODEs numerically: what do I need to consider in order to figure out which solver to use?

I would prefer recommendations that don't require the use of proprietary tools (such as Matlab). I know of two ODE solving options for the Python ecosystem: PyDSTool (Dopri, Radau, other Runge-Kutta ...
0
votes
1answer
29 views

Why might the time taken to compute the solution of an ODE system over some interval increase non-linearly with increasing size of interval?

Currently, my problem requires me to solve a system a large system of non-linear ODEs (up to ~5000). So far, I have been using scipy.integrate.odeint as my ...
1
vote
1answer
33 views

Possible to reduce effort needed to solve non-linear ODEs by taking some coefficients/parameters as constant over small time intervals?

So far, I have been using scipy.integrate.odeint as my "workhorse" ODE solver. My current problem requires that I solve a large system (up to ~5000) ODEs. Here's ...
1
vote
1answer
103 views

PETSc or Trilinos for GPU?

I have to choose between the PETSc and Trilinos libraries for a project that involves the solution of non-linear systems of equations. From their web sites they both mention support for GPUs, ...
0
votes
0answers
33 views

suggestion needed for solution algorithim non-linear differntial algebraic system

I have a nonlinear differential algebraic system in time arising from the weak formulation of coupled transient non-linear 1D problem. The system roughly looks like. ...
6
votes
1answer
87 views

Numerical methods for boundary-value ODEs with a jump condition

I want to solve a non linear system of equations of a particular kind. I find it hard to formulate clearly so I directly give a simple example. $ f''=A(f,g)\\ g''=B(f,g) $ with the boundary ...
6
votes
2answers
718 views

Is the shooting method the only general numerical method for solving nonlinear boundary value ODEs?

During my wandering in Mathematica.se, I gradually noticed that a certain kind of differential equation solving problem is "troubling" us all the time, that is, the boundary value problem (BVP) of ...
5
votes
2answers
185 views

Implementation of nonlinear term in FEM

Although there are similar questions, I am also struggling with the implementation of the following term in "my own code" by Finite Element Method, namely, $\nabla \phi \cdot \nabla \phi$. $\phi$ is ...
1
vote
1answer
156 views

Solving a system of nonlinear equations with an ODE solver is faster than with the Newton method?

This is somehow unexpected, but my recent experience with solving a system of nonlinear equations is that treating them as the right hand side of a system of ordinary equations and then evolve the ...
-1
votes
1answer
92 views

solving a non-linear equation with integrals involved

I would like to solve the following equation, wrt $n(e)$ $$f(n(e))=g(n(e)) + \int_{\alpha}^{e} w(n(x))dx $$ The integral there it confuses me. Any suggestion on how I can implement this on a the ...
5
votes
2answers
118 views

What are the numerical methods for huge polynomial systems?

Let a system of $n$ polynomial equations of degree $d$ with $m$ variables. I'm interested in a sparse system with $d = 3$, $n \sim 2000000$, $m \sim 50000$ and integer coefficients. What techniques ...
0
votes
1answer
26 views

Non linear system of equations with discretization on k-space

I want to numerically solve the following system of differential equations at the steady state: \begin{equation} \begin{aligned} \frac{\partial \rho_{11\mathbf{k}}}{\partial t} =& ...
0
votes
1answer
66 views

Finding roots without knowing much about the function

Consider solving numerically for roots: $( x_0, y_0): f(x_0, y_0) = 0, g(x_0, y_0) = 0$ where you only know that f, g continuously differentiable but the theoretical differentiation is not a ...
3
votes
1answer
127 views

Time Integration of a nonlinear reaction-diffusion system

I want to solve the following system of nonlinear reaction-diffusion equations (Schnakenberg Turing) using FEM methods (such as deal.ii): $$ \partial_{t} u = \Delta u + \gamma\left(a-u+u²v\right)$$ ...
0
votes
0answers
20 views

Model actuators in non-linear beam element structure

I have written a non-linear finite element code using beam elements, and I'd like to represent hydraulic actuators in the model. In real terms, this would mean that the distance between two nodes at a ...
3
votes
2answers
95 views

Non-linear root finding with positive definite Jacobian

I am dealing with a system of non-linear equations: $$ f(\boldsymbol{x}) = \boldsymbol{y}, \;\;\; \boldsymbol{x}, \boldsymbol{y} \in \mathbb{R}^d. $$ And I know that the Jacobian $J(\boldsymbol{x})$ ...
2
votes
0answers
121 views

Integration of nonlinear PIDE via spectral methods

At the mean-field level, the dynamics of a polariton condensate can be described by a type of nonlinear Schrodinger equation (Gross-Pitaevskii-type), for a classical (complex-number) wavefunction ...
4
votes
1answer
280 views

Coupled nonlinear PDEs with time dependence on the RHS

I would like to numerically solve the following system of 2 coupled partial differential equations for the unknown functions $\psi_X(x,y,t)$ and $\psi_C(x,y,t)$: $\partial_t \psi_X = -i\psi_C - ...
0
votes
1answer
114 views

Solver for large non-linear system of equations

I am curently using R package nleqslv for solving a non-linear system of equations with 300 variables. I need to scale this to the system with ~50k variables and naturally this does not scale very ...
4
votes
1answer
117 views

Fast way to repeatedly solve a small nonlinear equation system

A small nonlinear equation system (sizes around 12 ✕ 12) needs to be solved repeatedly (millions of times); each time with some variation in parameters/coefficients (although the equation set is ...
0
votes
1answer
83 views

MATLAB Newton non-linear equation

I have the following non-linear equation: where $w0=0.25,w0=0.5,w0=0.75$. I have to prove that if $k$ is a root, then also $−k$ is a root and that there exists only one $k∈(0,1)$ root, but my ...
0
votes
2answers
174 views

Systems of nonlinear equations

Consider the nonlinear system of equations $$ (1) \quad qk^2a_1^2E^2+wna_0a_1AE+pnka_0^2a_1E+rn^2a_0^2A^2-rn^2a_0^3A^2+qk^2a_0a_1ABE-qk^2na_0^2E^2=0, $$ $$ (2) \quad ...
1
vote
1answer
358 views

System of nonlinear equations in MATLAB

I've got some problems solving (numerically) this system of equations. \begin{array}{l} 40 \cdot \cos (2t) + 105 \cdot \cos ({\theta _3}) - 75 \cdot \cos ({\theta _4}) - 91.924 \cdot \cos ({337.62}) ...
0
votes
0answers
106 views

Solving a non linear equation and iterating for various values in python

problem:Unable to solve a unknown for multiple known values in a non-linear equation ...
5
votes
2answers
130 views

Convergence issues for a non-linear system

I have a nasty system of coupled integral equations, which I managed to discretize and recast a non-linear system, i.e. something like: $$ \vec{w} = F \left( \vec{w} \right) \hspace{32pt} w \in ...
6
votes
1answer
169 views

Non-linear optimization using approximate gradient

I'm working with non-linear optimization for imaging, such as MRI and CT. Our problem is of the form $\|Ax-b \|_2^2+\lambda \|Wx\|_1$. $A$ is never formed explicitly, so we're limited to approaches ...
7
votes
2answers
130 views

For a non-linear PDEs should the source term be discretised at $u_j$ or averaged over $(u_{j+1} + u_{j-1})/2$?

The non-linear Poisson equation in one-dimension, $$ 0 = \frac{\partial^2u}{\partial x^2} - f(u) $$ can be discretised as to give, $$ u_{j-1} -2u_{j} + u_{j+1} = h^2 f(u_j) $$ where $h$ is the ...
3
votes
1answer
134 views

System of quadratic algebraic equations

I have this problem $H_i(x_1,x_2,\dots, x_N) = a_{ijk} x_j x_k + b_{ij} x_j + c_i = 0 \quad 1\leq i \leq N$ And I need to show that applying Newton-Raphson can fail to find even one real solution ...
2
votes
1answer
68 views

Methods to solve this equation on finite fields?

Is there any analytical (exact, closed-form solution) or numerical method to solve an equation such as $p(x) = r^x$ where $p(x)$ is a polynomial whose coefficients are drawn from a finite field, ...
1
vote
2answers
295 views

finding wave function for anharmonic oscillator

I'd like to find the normalized ground state wavefunction for the anharmonic oscillator (Duffing) whose potential for which there is no analytic solution; an oscillator with a quartic potential, in ...
3
votes
2answers
268 views

Slow convergence of Newton's method for finite elements

The application is a simple non-linear advection diffusion problem (steady state) using DGFEM. My error at each iteration is given by $$ e_{n+1} = ...
10
votes
1answer
772 views

When is Newton-Krylov not an appropriate solver?

Recently I have been comparing different non-linear solvers from scipy and was particularly impressed with the Newton-Krylov example in the Scipy Cookbook in which they solve a second order ...
3
votes
0answers
247 views

Finding roots of systems of equations with a Jacobian that is singular everywhere

Given $a\in\mathbb{R}^{mn\times n}$, find a $C\in\mathbb{R}^{n}$, $x\in\mathbb{R}^{m\times n}$ such that $$ 0 = f_{k}(\boldsymbol{C}, \boldsymbol{x}):=\sum_{i=1}^{m} C_{i} \left(\prod_{j=1}^{n} ...