Questions tagged [nonlinear-equations]

Solution of nonlinear systems of equations. The equations might be algebraic or differential equations.

Filter by
Sorted by
Tagged with
1
vote
0answers
50 views

Fitting data with a Voigt function

I have some data, (xrd data), that I would like peak fit with a pseudo-Voigt function, a combination of a Gaussian and a Lorentzian function. These are the functions $G(x) = I \exp\left( -\frac{4\ln(2)...
4
votes
1answer
96 views

Method to linearize highly nonlinear partial differential equation

I have a set of coupled pdes which I want to solve using finite-difference, of which one is nonlinear. The three linear pdes for quantities $T_f$, $T_s$ and $c$ are convection-diffusion-reaction-like ...
0
votes
0answers
50 views

Best way to solve system of quadratic forms

I have a system of equations that have the following structure. Let $x\in\mathbb{R}^m$ and let $x_k$ be the $k$-th element of $x$. Let $H_k\in\mathbb{R}^{m\times m}$ for $k=1,\ldots, m$. I need to ...
1
vote
0answers
90 views

Is it possible to use a fixed point iteration for solving this nonlinear system?

Consider the following differential equation \begin{align} \frac{\partial f(u)}{\partial x} &= g(x), \ \ x\in [x_{L},x_{R}] \label{Eq2.2} \\ u(x_{L}) &= g_{1} \end{align} where $f(u)$ is a ...
5
votes
1answer
125 views

Analysis of nonlinear finite element methods

I have been doing a lot of reading on the development of finite element methods and their analysis using, e.g., functional analysis. I am clear on the formulation of the weak form of a PDE and ...
1
vote
1answer
46 views

Verification of coupled system of equations for light propagation

I am trying to simulate the propagation of light in material using the non-linear schrödinger equation (NLSE): $$\partial_zE=\frac{i}{2k_0}\nabla^2_\perp E+\frac{ik_0n_2}{n_0}\vert E\vert^2E-0.5\beta^{...
3
votes
0answers
124 views

How to treat nonlinear radiation term in heat equation using Finite-element method?

I am trying to solve the time-dependant non-linear heat equation with radiation. This equation is coupled to the radiative transfer equation but for the purpose of my question, this does not matter ...
0
votes
0answers
55 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 ...
3
votes
0answers
84 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 & & \...
2
votes
1answer
149 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
0answers
101 views

How does the error work for the Strang Splitting?

We know in Strang splitting that the splitting error in the steady state solution is proportional to $h^2$. I want ask 2 things: If this error in the steady state solution is the global error? If we ...
0
votes
1answer
131 views

Discretization of a non-linear ODE using FDM isn't grid indepenent

I am trying to solve the ODE : $\frac{d^2T}{dx^2} = \omega_1 T+\omega_2 T^2$ + using different numerical methods. I have tried the following discretizations so far and none of them seem to be grid ...
1
vote
0answers
32 views

Approximation in the derivation of the Arc Length method

I am studying the proof of the Arc Length method in section 2.2 of this thesis. In equation (2.2) the author introduces the supplementary conditions $$ (\Delta {\bf u} + \delta {\bf u})^T \cdot (\...
0
votes
0answers
56 views

Proof of kinematic relationship between updated and total Lagrangian 2nd Piola-Kirchhoff stresses

On page 587 of Finite Element Procedures by Bathe the author gives the following kinematic transformations $$ {}^t\tau_{ij} = \frac{{}^t\rho}{{}^o\rho} \; {}^t_ox_{i,r} \; {}^t_oS_{rs} \; {}^t_ox_{j,...
1
vote
0answers
127 views

Solving coupled PDEs with self-consistency condition

I am figuring out how to attack a problem (the Usadel equations of superconductivity) in which I need to solve a set of nonlinear PDEs for the fields $\{G_i (r)\}$ $$ U(G_i(r), \nabla G_i(r), \Delta(r)...
1
vote
0answers
94 views

Stochastic differential equation system (SDE) : overflow encountered in double-scalars

I'm trying to integrate the following SDE system from Dekker et al. [1] $$\begin{cases} \frac{dx}{dt}=a_1x^3+a_2x+\phi+\zeta_x\\ \frac{dy}{dt}=b_1z+b_2(\kappa(x)-(y^2+z^2))y+\zeta_y\\ \frac{dz}{dt}=...
1
vote
0answers
63 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
2answers
123 views

diffusivity matrix assembly in nonlinear finite element analysis

I want to solve a diffusion analysis using finite elements. According to fick's law, governing equation is $$\frac{\partial h}{\partial t} = D \frac{\partial^{2} h}{\partial x^{2}}$$ . h is relative ...
2
votes
0answers
73 views

Linearising Nonlinear Coupled Partial Differential Equations - Alfvénic Diffusion

I am trying to solve the following coupled partial differential equations with a finite difference scheme: $$\partial_tf+v\partial_zf+\partial_z\frac{1}{W}\partial_zf=0$$ $$\partial_tW+v\partial_zW-\...
1
vote
0answers
75 views

Improve code of logarithmic quantizer

I am implementing a logarithmic quantizer which is defined as follows: $$ q(u) = \begin{cases} u_i , \frac{u_i}{1+\delta} < u < \frac{u_i}{1-\delta} \\ 0 , 0 \leq u \leq \frac{u_o}{1+\delta} \\ -...
3
votes
1answer
119 views

Nonlinear root solving libraries which accept a Jacobian in band-storage

I'm in search for a library for solving large systems of non-linear equations, similar to MINPACK, but unlike MINPACK, can accept a Jacobian in band-storage. My Jacobian is sometimes not invertible, ...
4
votes
1answer
308 views

Backwards Difference Implicit Method for Nonlinear Parabolic PDE in Python

Original Stack Overflow Question: https://stackoverflow.com/questions/65683788/indexerror-index-31-is-out-of-bounds-for-axis-1-with-size-31?noredirect=1#comment116218335_65683788 PDE: u_t = u_xx + u(...
4
votes
0answers
56 views

Continuous limit and nonlinear functional analysis

I have a kind of general question about approximation schemes in nonlinear functional analysis. Given a nonlinear map $\Phi$ from an open set (in an infinite dimensional Banach space) of functions to ...
0
votes
1answer
391 views

Solve non-linear equation in R

I need to solve the following equation for $x$ in [0, 1]. Assume $0<\alpha<1$ and $0<\lambda$. $$(1 - x)^{\alpha+1} - \lambda (x+1)^{\alpha+1} = -2\lambda (\alpha + 1) x^\alpha$$ Would very ...
3
votes
1answer
183 views

How to derive the simplified Newton iteration in the TR-BDF2 ODE integration scheme

The Problem The TR-BDF2 explained in this paper [1], is quite a popular numerical scheme used to integrate $\dot{y} = f(t,y)$, consistent of the following two stages: \begin{align} y_{n+\gamma} &...
1
vote
0answers
217 views

Optimize speed for calculating the approximate inverse of a large matrix

I am searching for a faster method to calculate an approximate inverse of a large matrix (up to 32000x32000) resulting from a discrete non-linear system of partial differential equations. I'm using C++...
2
votes
2answers
180 views

Numerical Methods of solving a non-linear ODE?

I want to solve the nonlinear equation $\frac{d^2x}{dt^2} + k\sin x = 0$, numerically. I found that solving this elliptic integral would be cumbersome, so is there a numerical method i could use to ...
5
votes
0answers
76 views

2nd-order TVD criteria for flux-limiter

Consider a nonlinear hyperbolic conservation equation: $$ \partial_{t}u = -\partial_{x}f(u) $$ The spatial derivative of $f(u)$ may approximated after a spatial discretization by $x_{j}=j\Delta x$ $$ \...
1
vote
0answers
68 views

Numerically estimate the Jacobian from a scalar time-series

I'm trying to numerically estimate the Jacobian from a time-series. Following the paper, https://journals.aps.org/pra/abstract/10.1103/PhysRevA.43.2787. Consider that I have a scalar time series $x = (...
0
votes
1answer
131 views

Interpreting multivariable root-finding results from Matlab's fsolve algorithm

Edit: So I was able to get the same value of r that's given, when coding up the sum of squares of function values directly in the script file, rather than on the Command Window. So, maybe there's a ...
2
votes
1answer
155 views

Solving system of nonlinear vector functions

I am trying to figure out how to implement a solver for a system of nonlinear equations of the form \begin{align*} u_1 &= y_n + h\left(a_{1,1}f(t_n + c_1 h, u_1) + a_{1,2}f(t_n + c_2 h, u_2)\...
1
vote
1answer
76 views

How to effectively find a starting point solving a non-linear equation?

I have the following equation (the Kurz-Giovanola-Trivedi model [1]) $$ v^2 \frac{\pi^2 \Gamma}{P^2 D^2} + v \frac{mC_0(1-k)\xi}{D[1-(1-k)Iv(P)]} + G = 0, $$ where $Iv(P)=P \cdot \exp(P) \cdot E(P)$, $...
0
votes
0answers
78 views

Derivative-free ill-conditioned non-linear least squares

I am looking for a package which can solve (non-linear) least squares problems without the use of derivatives (because of an expensive model), but which also deals with ill-conditioning well (such as ...
1
vote
2answers
711 views

Jacobians with automatic differentiation

I have an objective function F: Nx1 -> Nx1, where N>30000. There are many sparse matrix/tensor multiplications in this function, so taking an analytic Jacobian by paper and pen is cumbersome. ...
3
votes
1answer
134 views

How avoid square shape with Laplacian operator in reaction diffusion calculations?

I have used different variants of the Laplacian operator (div grad) using 4, 8, 12, 20 and 24 of the closest points. I get problems due to the chosen coordinate system and the discretization of the ...
1
vote
0answers
78 views

RK4-method starts oscillating above certain input parameters

I am trying to solve an equation of the following type $$\partial_zE(z)=-c_0J$$ with $$J=c_1\beta E^3(z)$$ using the boost::odeint-framework and a fixed time stepper, with $c_0$, $c_1$ and $\beta$ ...
1
vote
1answer
77 views

Issue solving nonlinear equation containing a quotient

I have a coupled set of PDEs that need to be solved as part of a larger problem. I am currently approaching this by computing spatial derivatives with finite differences and using PETSc's nonlinear ...
2
votes
2answers
97 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
1answer
132 views

Method of Lines Runge-Kutta nonlinear stability and behavior

I have a system of 4 nonlinear 1st-order PDEs. I want to solve them numerically by method of lines, first discretizing space. This leads to the system of $N\times 4$ coupled ODEs. $$ \mathbf{u}_{i} =...
3
votes
0answers
92 views

Is the matrix exponential and the Jordan canonical form actually useful for solving differential equations?

All of my yearlong graduate-level Linear Algebra course notes from my professor—an algebraist/representation theorist—shows his love for the exponential map $e^A$ and the Jordan canonical form—and one ...
2
votes
1answer
209 views

Numerical methods that can be written in flux conservative form

I have a system of non-linear PDEs that I expect to have shocks as well as the appearance of Gibbs phenomena (spurious oscillations that form near the shock) for 2nd-order methods or higher. I have ...
3
votes
0answers
70 views

Conjugate Gradient for nonlinear equation system

Is it possible to apply adaptions of the conjugate gradient algorithm i.e. Fletcher-Reeves, Polak-Ribere or others to systems of nonlinear equations? How should the equation system be adjusted so one ...
0
votes
1answer
125 views

FDM on nonlinear PDEs

I'm working with a 2D Navier Stokes PDE in the unstabilized version - the equation is a linear equation of the type $\frac{∂u}{∂t} = F(u,t)$. In order to perform time discretization with FDM (finite ...
0
votes
1answer
22 views

Python: Getting second output variable from minimizing a computationally intensive function on first outputs

I have a function in python that is quite computationally expensive to evaluate, of the form: ...
0
votes
1answer
39 views

Coupled pdes of the first order

May question is about possible approaches to solve the following system $$ \begin{array}{rcl} \nabla{n}&=&n\,\mathbf{E},\\ \nabla\cdot\mathbf{E}&=&1-n, \end{array} $$ in general with ...
1
vote
0answers
81 views

Bifurcation points on homotopy path by numerical continuation?

I am trying to implement an algorithm that finds (possibly) all solutions of a system of nonlinear polynomial equations $$F(X) = 0$$ I thought about using the (convex) homotopy $$H(X,T) = TF(X) + (1-T)...
4
votes
2answers
80 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
0answers
97 views

Linearize non-linear PDE with BCs to hyperbolic problem: How does linearization affect BCs?

I am working with the Shallow Water equations that is a system of non-linear PDEs that simulate water waves propagation on some domain, in my cases the $x$-axis. I have here a GIF showing the results (...
1
vote
1answer
103 views

Full approximation scheme - smoothers - literature recomendation

I would like to use full approximation scheme(FAS) for solving nonlinear PDE. I am looking into practical implementation of FAS in parallel(MPI) settings. I noticed the most common/effective smoother ...
0
votes
0answers
76 views

Solving a nonlinear problem with a very small components with finite element method

In solving nonlinear hyperelastic solid mechanics problems, to converge to the correct solution we need to do step-by-step loading which makes the deformation at each step very small (for my ...

1
2 3 4 5
7