Questions tagged [numerical-modelling]

The tag has no usage guidance.

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

Mineral dissolution and solute transport around a solid

I am trying to simulate solute transport of acid (HCl) and consequent mineral dissolution around a grain (calcite). The governing equation for transport is the advection-diffusion equation, given as: ...
3
votes
1answer
105 views

Solving and Plotting Mutualism Model in Python

I am a beginner in programming. I need to program a mutualism model of two species in python that would solve and graph using the following equations: $$ \frac{dN_1}{dt} = N_1(r_1 - e_1N_1 + \alpha _{...
2
votes
2answers
128 views

What is the most suitable numerical approach for modelling multiphase flow with particle interactions?

If I want to build a solver for this following problem: 1. There is stagnant water governed by the Navier-Stokes equation in the domain. ...
-1
votes
0answers
33 views

Modeling of gestation period using Root library

The "gestation period" for humans averages 40 weeks, with a standard deviation of 10 days. Using CERN's root package, calculate and plot the probability that a woman will give birth tomorrow,...
2
votes
0answers
63 views

Solute transport around a solid obstacle

I am a newbie in CFD and single/multiphase flow and transport in general. As part of my quest to learn, I am trying to model solute transport around a solid object in the center of a 2D domain. The ...
2
votes
0answers
95 views

Floquet theory for periodic delay differential equations: current numerical routines

I would like to determine the stability of a system of periodic delay differential equations (a seasonal host-parasite model). I've tried to implement the method described in Lemma 2.5 in this paper: ...
1
vote
1answer
81 views

Discrete model of cell - cell communication

I am trying to understand how cell to cell communication is studied using a discrete modelling framework. Could someone please suggest suitable references or libraries which already have ...
2
votes
1answer
63 views

Finite element modelling of thermal expansion of 3D solid bodies

I want to solve the thermal expansion of solid by using FEM approach. When I developed the model based on the the principle the minimum potential energy, the solutions for thermal expansion are not ...
1
vote
0answers
18 views

How can this attempt to implement a Milstein numerical approximation for an Ito process of multiple components be fixed?

I have been reading Kloeden and Platen's Numerical Solution of Stochastic Differential Equations, and have more or less been trying to systematically complete the various exercises therein as I go ...
4
votes
0answers
74 views

Identifying an unknown P.D.E. from solution data

I have a black-box simulation that produces the time evolution of a probability density function p(x, t) in 1 dimension from arbitrary initial conditions p(x, 0). The underlying simulation occurs on a ...
0
votes
1answer
61 views

Does the leap-frog algorithm conserve energy for n-body problems?

The leap-frog algorithm is able to conserve to a certain extent the energy of a system, which flucutates as a cosine around a stable value. Is this true if we apply the algorithm to a n-body ...
1
vote
0answers
37 views

Procedure to convert continuous equations of motion to discrete version

Let's take a mobile robot and suppose we know its continuous equations of motion, for example this car-like simple model. Now if I am simulating this robot in a continuous 2D plane, then coding the ...
0
votes
0answers
53 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 ...
2
votes
1answer
113 views

How to couple the vibro-acoustic equations by Mortar method for non-matching meshes?

Assume we have two domains $\Omega_a$ a acoustic domain with boundary $\Gamma_a$ and $\Omega_s$ a domain of a solid body with boundary $\Gamma_s$. $\Omega_a$ and $\Omega_s$ have the common interface $\...
2
votes
0answers
50 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 ...
2
votes
2answers
304 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 ...
0
votes
2answers
356 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 ...
0
votes
0answers
74 views

How to distinguish primary hosts (stars) and orbiting satellites (planets) and tertiary bodies (moons) by their mass and trajectory?

I posted this question in the astronomy stackexchange. There are no responses, and it was suggested that I pose the question here. The "too long, didn't read" was taken from a comment, and ...
5
votes
2answers
563 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 ...
0
votes
1answer
67 views

How error scales with numerical precision in molecular dynamics?

In terms of time-step, numerical error in molecular dynamics scales with square, i.e. $error \approx dt^2$. But how it look for numerical precision ? E.g. how much bigger will be numerical error when ...
0
votes
0answers
55 views

How to efficiently perform this 2D integral in Quadpy?

I need to integrate a function defined in 2Dims (z and radius r), for which I don't have an expression. I can just query the ...
0
votes
1answer
55 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 ...
2
votes
2answers
144 views

Generating particles from a distribution function using Monte Carlo

I have been given a 4D ($x, y, v_x, v_y$) distribution function, $f(x,y,v_x, v_y)$, generated by an external code. I want to generate a set of particles from this distribution function, say 10k ...
-2
votes
1answer
48 views

Solving differential equation by setting vectorization `on` in MATLAB

This is a follow up to my previous question posted here. I've set up an ode system in MATLAB and I'm trying to vectorize the code to increase the speed of computation. The follow is the code for my ...
0
votes
0answers
37 views

What is an integrated modelling tool?

In this document they describe their code as an integrated modelling tool. I am trying to understand how this is different from a regular modelling tool? Most codes include different modules that do ...
0
votes
1answer
85 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 ...
1
vote
1answer
56 views

solving differential equations with jacobian pattern

I'm trying to compare the simulation time for solving a system of differential equations with and without jacobian pattern for a toy model using ode15s in MATLAB. ...
0
votes
0answers
46 views

Estimating the Jacobian in Harmonic Balance Method

I am trying to solve a set of ODEs using the Harmonic Balance method. In order to do this, I need to compute the Jacobian of the set of equations. However I am very confused regarding the dimensions ...
0
votes
1answer
81 views

Acoustic Simulation, how are boundaries handled?

I don't have a background in numerical modeling so this question is rather broad. What I am interested in is modeling the propagation of an ultrasonic acoustic wave in 3d space. The basic 3d wave ...
4
votes
2answers
132 views

Obstacle too thin in fluid-structure interaction, so I consider it as membrane

I need to simulate a 3D fluid-structure problem where the obstacle is an elastic and very thin structure (then I want to consider this structure as a surface). I need to solve this problem using ...
6
votes
0answers
105 views

A Question About a Claim from 1991 Computational EM paper about the Cancellation of certain Boundary Terms

Please let me know if this is not the appropriate site for this question. I found questions regarding EFIE/MFIE/CFIE on this site, so I thought my question might fit. I am studying the paper by Putnam ...
0
votes
0answers
46 views

Unsteady Stokes equations in ALE framework

I'm trying to solve Unsteady Stokes equations on a moving domain, using an ALE formulation, that is $$\frac{\partial \mathbf{u}}{\partial t} - \mathbf{w}\cdot \nabla\mathbf{u} = \nu\Delta\mathbf{u} - \...
1
vote
2answers
444 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 ...
3
votes
1answer
139 views

Choosing an appropriate time step for a discrete & continuous dynamics simulation

I have inherited of a flight dynamics simulation in C++ which represents a small drone with it's autopilot, actuator dynamics and a solid state IMU. Hence, it is composed of a few models, some ...
3
votes
1answer
223 views

Dirichlet Boundary Condition finite difference method using sparse-matrix $Ax = b$ system

I am trying to solve the boundary value problem for heat equation: $$ u_{xx} + u_{yy} = f(x,y) $$ where the solution $u(x,y) \in [0,1] \times [0,1]$ and the Dirichlet boundary condition $u(x,y) = ...
6
votes
2answers
131 views

What makes a good computational grid?

Most computational methods for solving PDEs are grid-based. What makes a computational grid "good", other than being sufficiently fine to resolve features of numerical solutions? Are grids ...
1
vote
0answers
86 views

How can I practice multivariable root-finding?

Recently, I've been reading up on various root-finding / optimization algorithms such as the Levenberg-Marquardt method, Gauss-Newton, Conjugate Gradient, trust-region and trust-region-dogleg. I've ...
2
votes
2answers
61 views

In a dynamical system, what might be a good reason why periodicity in an object's velocities is important?

I'm studying periodic motions in a dynamical system and, as a newbie, I narrowly think of an object's periodicity in its spatial x-y coordinates, but what might be a good reason why the existence of ...
1
vote
1answer
84 views

Would modeling a fluid flow that decreases in magnitude violate conservation of mass?

If I make a quasi-steady assumption in a model such as keeping the fluid density constant and assuming the flow is incompressible, does modeling a flow with decreasing velocity / magnitude, say, as ...
0
votes
0answers
65 views

The error in SOR algorithm suddenly falls to zero when it reaches 1e-7 range

I am solving the Poisson equation for heterojunction using Fortran90. I use the SOR algorithm to arrive at the potential profile. I see the weird behavior where the error (the difference between the $...
1
vote
1answer
487 views

Reading VTK file into C++ for analysis

I apologize in advance if this post is at all ignorant or elementary, I am a pure mathematician who is newly getting into the world of scientific computing. For my research, my advisor would like me ...
2
votes
0answers
54 views

Solving Stokes Equations in 3D - Do I need to treat pressure-velocity coupling iteratively?

I need to develop a code to solve Stokes Equations in 3D in cubic geometries (structured grid, uniform mesh spacing). My code needs to take a pressure gradient in one direction as a BC (pinlet=p1, ...
0
votes
0answers
31 views

At what l/d ratio will a frame element start to behave as a shell element?

I'm working in ETABS. There are few columns of dimension 300mmx1400mm. The height of building is 36.6 meter above ground level and the dimension of building is 26mx68m. I'm getting the time period of ...
3
votes
0answers
89 views

How to construct a Fortin Operator for Crouzeix-Raviart Element?

I want to prove the LBB condition for the Stokes Equations discretised by the Crouzeix-Raviart element. The continuous Stokes Equation in the weak formulation is Find $u \in H_0^1(\Omega, \mathbb{R}^...
4
votes
1answer
1k 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 ...
3
votes
1answer
169 views

Suggestions for scientific computing projects to undertake to sharpen core skills

Currently I am completing all exercises in books like "Introduction to Python for Science and Engineering, David Pine" and "Guid to Scientific Computing in C++, Pitt-Francis, Whitley". I am looking ...
1
vote
0answers
66 views

Finding isocontour lines, populating points on said lines (in MATLAB)

I'm trying to extract data from a potential field. I have an example code of a potential field, using meshgrid in MATLAB. I've pasted the code below. What I'd like to do next is 'convert' this field ...
1
vote
0answers
41 views

Advice for a topic in a seminar

I am a master student in Mathematics, and I have to prepare a seminar for a course in mathematical methods for applied sciences. I have a good background in numerical analysis for ODEs, PDEs and hence ...
5
votes
1answer
281 views

What kind of a researcher am I?

So far, I've worked a bit in modeling, simulations and simple lab experiments, and I've really enjoyed all three research methods to approach a single research question. I can write tricky (in terms ...
0
votes
0answers
73 views

Numerical flow visualization in 2D for a moving boundary,

I have a rigid body that moves according to a set of governing ODEs, and I'd like to numerically visualize the vortices that are shed by this object. How could I proceed? I've been reading up on the ...