Skip to main content
Share Your Experience: Take the 2024 Developer Survey

Questions tagged [simulation]

Simulation is the process of mimicking the behaviour of a real physical process or system using a representative model.

82 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
7 votes
0 answers
203 views

Is a complete bacteria simulation with an exascale supercomputer possible?

Will it be possible to simulate a complete (at least simple) bacteria atom by atom on an exascale supercomputer? or is it possible already today with the largest systems? Here, I've read that ...
Open Food Broker's user avatar
6 votes
0 answers
341 views

FEM : energy minimization VS PDE solving

Engineering FEM When I studied engineering, I learned the traditional approach for finite elements for elasticity. The point was to solve the PDE $-div(\sigma)=f$ as: Multiply your PDE with a test ...
Txnda's user avatar
  • 61
6 votes
0 answers
1k views

How does GAMG in OpenFOAM really work?

I use OpenFOAM for CFD simulations. A very popular preconditioner is GAMG which needs a low number of iterations per a time step in SIMPLE or PISO solvers that are used to simulate the fluid flow. I ...
user avatar
5 votes
0 answers
164 views

How can Navier--Stokes equations have asymmetric solutions such as Karman vortex streets

The Navier--Stokes equations are axially symmetric, so with symmetric boundary conditions, how can features such as Karman vortex streets develop? I understand that in reality symmetry does never ...
Bananach's user avatar
  • 799
5 votes
0 answers
268 views

Sound Waves Simulation in 3D Environment

I want to do a simulation of sound waves including wave propagation, absorption, and reflection in 3D space. I did some research and I found this question in stackoverflow but it talks about ...
eyadof's user avatar
  • 151
5 votes
0 answers
121 views

What are the governing equations solved in coupled atmosphere-ocean models?

In (hydrostatic) atmospheric general circulation models, for example the so-calle Primitive Equations, consisting of the horizontal momentum equation, the hydrostatic balance, the continuity equation, ...
Dilaton's user avatar
  • 151
4 votes
0 answers
86 views

How to do a parametric study of arbitrary 2D surface?

I need to do a parametric study of the performance of a room heater on different rooms by simulating the temperature distribution in there. The problem here is that the rooms are not simple rectangles ...
Ken Grimes's user avatar
4 votes
0 answers
105 views

Any literature that extensively discusses the ability to strongly inter-/extrapolate computationally on little empirical data?

Any literature that extensively discusses the ability to strongly inter-/extrapolate computationally on little empirical data? This topic has fascinated me, but I find that it seems a bit novel. ...
mavavilj's user avatar
  • 427
4 votes
0 answers
93 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 ...
beables's user avatar
  • 41
4 votes
0 answers
105 views

Tracking for two meshes

I'm dealing with physical simulation (position based dynamic). Now I'm trying to realize tracking for two meshes. In order to explain what does it mean, let's assume that we have two similar meshes: ...
Stepan Loginov's user avatar
3 votes
0 answers
104 views

Verlet integration on grids or how to get better stability in hyperelastic simulation

I am using MLS-MPM to simulate both solid and fluids. It works, but the amount of time steps I must do for hyperelastic solids is absurd. To give you some perspective, I am able to simulate just the ...
Makogan's user avatar
  • 263
3 votes
0 answers
70 views

Spectral solver on em-pic

I'm recently studying for the spectral solver to implement EM-PIC code. I read an article and have some questions. Many PIC codes uses spectral solver to overcome numerical artifacts on FDTD. In the ...
asdgaaa1123's user avatar
3 votes
0 answers
360 views

Error in MD simulation with Lennard-Jones potential in C, resulting in a lack of conservation of energy

I'm working on a MD simulation using the Lennard-Jones potential and the Velocity Verlet algorithm. The code compiles and runs, but the energy for each time step is not conserved; it decreases ...
Luciano's user avatar
  • 151
3 votes
0 answers
119 views

Is it normal to expect the error of simulation of a damped harmonic oscillator to decrease as the damping factor decreases?

I am simulating a damped harmonic oscillator using the RK4 method of numerical integration. I am comparing the simulated results with the analytical ones (for the free evolution case) and obtaining ...
turnip's user avatar
  • 263
3 votes
0 answers
61 views

Discretizing boundary conditions for vortex methods

I am working on a fluid simulation using vortex methods. For this I must compute the vortex sheet on my boundaries given as: $$ \gamma(\mathbf{x}) - \frac{1}{\pi}\int_S \frac{\partial}{\partial\mathbf{...
Alex's user avatar
  • 31
3 votes
0 answers
260 views

N-body simulation for particles

I would like use N-body simulation for measure a torque on an object in a liquid (helium). It's a 2D study. On the site : http://www.browndeertechnology.com/docs/BDT_OpenCL_Tutorial_NBody-rev3.html I ...
Ludovic's user avatar
  • 31
3 votes
0 answers
220 views

How to estimate if a velocity field is statistically homogeneous?

For a 3D velocity field $\mathbf{u}$ obtained by direct numerical simulation. Assuming that the field is defined on a periodic domain $\mathcal{P}$ of periodicity $L_x$ in the $x$ direction, $L_y$ ...
ucsky's user avatar
  • 419
2 votes
0 answers
49 views

In electromagnetic simulation, how much does the feed model impact an antenna's directivity performance and E-field phase readings?

Question: How much would a simplified feed model in an EM simulation alter an antenna's directivity and the E-field phase reading when compared to using a more complicated/realistic feed model? I am ...
Xingda Chen's user avatar
2 votes
0 answers
107 views

How to calculate the change in water level height after each step in a water displacement simulation?

I'm trying to write a very simple simulation of volume/water displacement. Every timestep a sphere gets lowered a bit further into a cylindrical container filled with water. While the physics behind ...
Johann Gruber's user avatar
2 votes
0 answers
170 views

Poisson equation solution in a semiconductor structure

I am trying to solve the $\textbf{1-D}$ Poisson equation for a semiconductor structure at equilibrium (There is no external bias applied). $\textbf{Background}$ \begin{equation} \frac{d^2V}{dx^2} = -\...
0-0's user avatar
  • 33
2 votes
0 answers
385 views

Numerical simulation of magnetic dipole in a inhomogeneous magnetic field

Goal: I want to use Python to illustrate how a magnetic dipole with magnetic moment m2 moves in a non-homogeneous magnetic field in a 2D-Plane. This field is generated by another magnetic dipole with ...
tester931's user avatar
2 votes
0 answers
44 views

Decoupling of the fluid viscosity and incompressibility constraint in a partitioned scheme

I am new here and this is my very first question: I hope I respect all the criteria and rules. I am just getting started with partitioned approaches for solving FSI problems. I am interested in a ...
user111283's user avatar
2 votes
0 answers
55 views

What is the best methodology for physics simulators of large floating base rigid body systems?

I want to implement a physics simulator for large floating base rigid body systems from scratch. The Rigid Body Dynamics Systems (RBD) should typically have the following characteristics: About ~50 ...
Jenny Reininger's user avatar
2 votes
0 answers
227 views

Initial conditions of nbody problem

Note: I have edited the original question to try to be clearer. I'm trying to simulate a gravitational nbody problem with a Barnes-Hut algorithm, in which the bodies start in a homogeneous sphere ...
RedPointyJackson's user avatar
2 votes
0 answers
52 views

Limitations of particle in cell method for high density plasma

Are there any limitations of particle in cell (PIC) method for high density plasma? To be more specific, is modeling of a narrow channel of high density plasma possible or are there any limitations ...
ziky's user avatar
  • 121
2 votes
0 answers
94 views

Tracking the speed of 2D oscillations on a lattice

I wrote a Monte Carlo simulation of the 2D Lotka-Volterra model on a discrete lattice (with periodic boundary conditions). A video that I produced (which images the system after some number of monte ...
Loonuh's user avatar
  • 253
2 votes
0 answers
259 views

Algorithms for adding hydrogens on a molecule

I want to add hydrogens to some linear polymer molecules (polyethylenes). I know some working methods like using PyMOL internal function h_add. This method works, but hydrogens are added at distances ...
Samanosuke Akechi's user avatar
2 votes
0 answers
79 views

Can similarity criterion help reducing calculations in Finite Element Analysis (FEA) problems?

Can similarity criterion help reducing calculations in Finite Element Analysis (FEA) problems? I'm using NX (a thermal and CFD simulation software) to simulate the temperature and fluid field of a ...
McBear Holden's user avatar
1 vote
0 answers
101 views

Can I reduce my simulation error with a staggered grid, postprocessing and compatibility equation feedback?

What I did Using the finite difference method, I solved with a certain amount of error the following system of hyperbolic partial differential equations in cylindrical coordinates (the problem is ...
FriendlyNeighborhoodEngineer's user avatar
1 vote
0 answers
47 views

Monte Carlo simulation of classical Heisenberg model doesn't represent Curie curve

I've created a JavaScript file to execute and log average energy and magnetization values of 2D lattice classical Heisenberg model. I run the simulation with parameters, ...
M. Çağlar TUFAN's user avatar
1 vote
0 answers
33 views

Normalisation in tensor networks

I am trying to implement the iTEBD algorithm for the $PXP$ model, i.e, the hamiltonian is $$H = \sum_iP_{i-1}X_iP_{i+1}.$$ Here $P$ is the projector onto the ground state and $X$ is the usual pauli x ...
Souroy's user avatar
  • 11
1 vote
0 answers
38 views

Can I build an electrostatic simulation of conductive objects using just the method of images?

I wish to build a simulation of self-assembling wires. Think lots of ball bearings in a 20KV electric field and they form emergent structures. If I use the method of images to figure out what the ...
Mandrake_The_Magician's user avatar
1 vote
0 answers
55 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 ...
user_1_1_1's user avatar
1 vote
0 answers
55 views

Reproducing a paper's result for Topological Insulators

For the past weeks I have been trying to reproduce Agarwala's results but I've been unsuccessful. From this paper I am trying to reproduce the first and last columns of Fig.2, by implementing eq.2; ...
Victor RM's user avatar
1 vote
0 answers
44 views

How to make a time-parametrization slower around a point but not too slow?

For animation purposes (see below) I need to use a parametrization which 'slows down' near a specific value; More precisely, I am looking for a 'nice' monotonically increasing function $r:[0,1] \to [...
Asaf Shachar's user avatar
1 vote
0 answers
95 views

computing time scale and steady state concentration in microfluidic channels

I have been performing convection-diffusion transport studies on microfluidic channels like the following The inlet concentration is specific and I obtain the time-dependent concentration profiles of ...
Natasha's user avatar
  • 421
1 vote
0 answers
76 views

Problems with simulation of a spatial filter 4f setup (Python)

I have a question about my code which computes numerically the output field of a 4f setup with a pinhole in the middle which works as a spatial filter. My setup consists of two lenses with 50mm focal ...
Nabla94's user avatar
  • 11
1 vote
0 answers
42 views

Discretization formula for a system of two differential equations. "Solution to one of these is the initial condition of the other". In which sense?

Consider the following stochastic differential equation \begin{equation} dy=\left(A-\left(A+B\right)y\right)dt+C\sqrt{y\left(1-y\right)}dW\tag{1} \end{equation} where $A$, $B$ and $C$ are parameters ...
Strictly_increasing's user avatar
1 vote
0 answers
39 views

Tackling multiscale problem in numerical simulation

In a dusty plasma system there are more than one component with different masses, i.e, electrons, ions,neutrals and dust grains. Accordingly, there are more than one temporal and spatial scales ...
user35569's user avatar
1 vote
0 answers
83 views

How to simulate a PDF using data samples like this?

I know two methods to simulate a PDF from random data samples using MATLAB : 1) Using a histogram where I use this command histogram(data,'Normalization','pdf'), it gives PDF like bins. 2) Another ...
Userhanu's user avatar
  • 151
1 vote
0 answers
42 views

How to avoid density getting "deleted" in two way rigid body coupling with LBM CFD?

I've been reading this paper recently, which talks about using Lattice Boltzmann methods and two way coupling. Specifically, it outlines fluid solid coupling, and solid fluid coupling, and how simply ...
Krupip's user avatar
  • 161
1 vote
0 answers
195 views

Using MATLAB to simulate the Ising Model

I am using MATLAB to simulate a 1D Ising Chain. I am running into an issue where when trying to find heat capacity, my system has a tremendous amount of noise. I'll post my code and an image of the ...
MurderOfCrows's user avatar
1 vote
0 answers
77 views

Magnetostatic modelling Radia: Increasing distance between magnets a produces positive force until a certain point, beyond which force goes haywire

I have two sets of magnets. One set consists of two electromagnets (Shown below: Blue) and the other set consists of two NdFeB N40 permanent magnets. Both sets of magnets lie on the same plane. I want ...
James O'Farrell's user avatar
1 vote
0 answers
164 views

Dealing with simultaneous collisions in N-body sim

I have written an 2d N-body simulation in Python which allows collisions between the bodies. A body is modeled as a circle whose area is proportional to its mass. Each time the sim advances by one ...
math_lover's user avatar
1 vote
0 answers
32 views

Does excitation type matter in a time-domain simulation/computation of a transfer function of an LTI system?

Let's say I am running a FDTD simulation of a wave-equation to determine a transfer function of an LTI system: \begin{equation} H(f) = \dfrac{Y(f)}{X(f)}\ \end{equation} where $Y(f)$ and $X(f)$ are ...
tomi.lee.jones's user avatar
1 vote
0 answers
43 views

Simulating periodic boundary condition with long_range interaction

Consider 2 particles in a 2D plane. There is a long-range interaction between two particles: $$f (r_i-r_j)= \frac{\vec{r_i} - \vec{r_j}}{|r_i-r_j|^3}$$ $$\vec{r_i} =(x_i,y_i)$$ $r_i$ determines the ...
sara nj's user avatar
  • 111
1 vote
0 answers
463 views

Simulating a Hill-Type muscle model

I'm trying to do biologically based locomotion animation. I have a Hill-Type Muscle Model (see reference 1) that I've implemented in code. I can find the initial lengths of the tendons, then do a step ...
Phylliida's user avatar
  • 323
1 vote
0 answers
75 views

Dielectric in an external electric field

I want to numerically solve for electric field of an object in an external electric field. The object has linear displacement field such that its permitivity ($\epsilon$) is constant. The object's ...
heaven-of-intensity's user avatar
1 vote
0 answers
196 views

RCWA for Non-Periodic Structure?

I have a quick question regarding Rigorous Coupled Wave Analysis (RCWA) for numerically solving problems (in optics). I am trying to simulate a device which is non-periodic in x, y directions (but is ...
jfherrm's user avatar
  • 11
1 vote
0 answers
66 views

Solving a nonlinear problem with CDF

I'm trying to solve this problem: $\begin{cases} \partial_t E=-k\left([f(\rho)-i.\left[\delta+\frac{1}{2}a\left[\dfrac{\nabla^2_{\bot}}{4}+1-\rho^2\right]\right]]E - 2CP\right)\\ \partial_t P=-\...
alxg's user avatar
  • 166