Questions tagged [numerical-modelling]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
2 votes
1 answer
212 views

Calculation of the EFIE integral

I need help computing the following integral: $$ \int_{}\frac{(1+jk|\vec{r}-\vec{r}^\prime|)e^{-jk|\vec{r}-\vec{r}^\prime|}}{|\vec{r}-\vec{r}^\prime|}d\vec{r}^\prime $$ in this integral $\vec{r}$ ...
hamed's user avatar
  • 21
7 votes
1 answer
669 views

Discrete wave simulation - absorbing boundaries?

I wrote a simple 2D wave simulation using the following equations: $$\frac{\partial^2 u}{\partial t^2}=c^2\nabla^2u$$ Where $\nabla^2$ is the discrete laplace operator using a Von Neumann neighborhood ...
Bloc97's user avatar
  • 171
1 vote
2 answers
1k views

How to set up the differential equation system to speed up computation?

I've set up a system of differential equations, obtained after discretizing pde, in the following way ...
Natasha's user avatar
  • 421
6 votes
1 answer
181 views

A simple PDE solution question

I need to ask a question about partial derivatives. I want to solve this equation (steady state, one dimensional continuity equation): $$\frac{\partial (\rho u)}{\partial z}=0$$ which is equivalent to:...
ASK22's user avatar
  • 71
5 votes
1 answer
2k views

Solving for a set of coupled ODEs to get correct variable values

My question is about how I can solve a coupled system of ODE's, and print out the variables in a plot. I am solving for an q value and an e value, seen in this set of coupled ODE's below: $$ \begin{...
Nicole I.'s user avatar
5 votes
1 answer
286 views

Which optimization method can be used to do the following?

I've the following system of equations for studying information flow in the below graph, $$ \frac{d \phi}{dt} = -M^TDM\phi + \text{noise effects} \hspace{1cm} (1)$$ Here, M is the incidence ...
Natasha's user avatar
  • 421
5 votes
2 answers
741 views

Specifying ode solver options to speed up compute time

I'm specifying the 'JPattern', sparsity_pattern in the ode options to speed up the compute time of my actual system. I am sharing a sample code below to show how I ...
Natasha's user avatar
  • 421
4 votes
1 answer
3k views

Finite difference methods in cylindrical and spherical co-ordinate systems

I am quite familiar with finite difference schemes in cartesian coordinates. The key point here is that every point in the cartesian grid is treated equally as the spacing between consecutive points ...
Babaji's user avatar
  • 195
3 votes
3 answers
382 views

What numerical methods are used to model deformations in elastic physics?

What numerical methods are used to model deformations in elastic physics? For example, here's an example of a hyperelastic deformation in Ansys: Perhaps more simply than hyperelasticity, for linear ...
wyer33's user avatar
  • 757
3 votes
1 answer
872 views

Question on how MATLAB's pdepe solver works

I'm solving the following 1D transport equation in MATLAB's pdepe solver. $$\frac{\partial C}{\partial t} = D\frac{\partial ^2 C}{\partial x^2}-v\frac{\partial C}{\partial x}$$ At the inlet (left ...
Natasha's user avatar
  • 421
2 votes
2 answers
299 views

Mathematical test method for the numerical solution of PDEs?

What are some of the methods used to test for the exactitude of a numerical solution, given that the analytical solution isn't available, and the numerical solution converges ?
Sam Gomari 's user avatar
2 votes
1 answer
1k views

Numerical method for a BVP with mixed boundary conditions (MATLAB)

I've been given a second-order non-linear ODE: $$\frac{d^{2}\theta(s)}{ds^{2}} = sf_{g}\cos{\theta} + sf_{x}\cos{\phi}\sin{\theta}$$ where $f_{g}, f_{x}$ and $\phi$ are constants. The boundary ...
Kyle Massey's user avatar
2 votes
2 answers
713 views

How to implement point source or volume source in finite element implementations

I'm trying to do a simple implementation to study the advection-diffusion-reaction dynamics in a straight pipe. I have points positioned along the length of the pipe (blue dots in the image above). I ...
Natasha's user avatar
  • 421
2 votes
0 answers
106 views

Does the time-dependent 1D advection-diffusion with point sources have an analytical solution?

I am looking for the analytical solution of 1-dimensional advection-diffusion equation with several point sources, Q, along the axial length of a cylinder through which the fluid flow occurs. Neumann ...
Natasha's user avatar
  • 421
2 votes
1 answer
773 views

Analytical Solution of Transport Equation

I'm looking at the analytical solution of the convection-diffusion equation $$\frac{\partial C}{\partial t} = D\frac{\partial ^2 C}{\partial x^2}-v\frac{\partial C}{\partial x}$$ with initial ...
Natasha's user avatar
  • 421
1 vote
1 answer
456 views

RK4 integration of the three-bodies problem with C++

first of all thank you for all the answers you gave me yesterday for the integration via Symplectic Euler's method of the three-body problem. We managed to implement both Euler's and Runge Kutta 4's ...
jack23456's user avatar
  • 171
1 vote
0 answers
85 views

PDE discretization (via finite difference sheme) question

So after posting this question and reading all your comments I would like to make this new question (update). If you consider the three equations presented here: $$\frac{\partial \rho}{\partial t} +\...
ASK22's user avatar
  • 71
0 votes
1 answer
273 views

Solving differential equation by specifying jacobian pattern

This is a follow up to my previous question posted here I'm trying to construct the sparsity pattern of the jacobian matrix to speed up the computation of a large system of odes. The following is the ...
Natasha's user avatar
  • 421
0 votes
1 answer
108 views

finding boundary conditions when transforming a higher order ode to system of first order ode

given the following ODE: $$\frac{d^{4}w}{dx^{4}} + B\frac{d^{2}w}{dx^{2}} = 1$$ with boundary conditions $w(0) =0 , w(1) = 0,w'(0) = 0,w'(1) = 0$ its possible to solve analytically but I am attempting ...
Gideon Ilung's user avatar
0 votes
2 answers
415 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 ...
Haus's user avatar
  • 123