Questions tagged [simulation]

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

53 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
7
votes
0answers
139 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 ...
5
votes
0answers
126 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 ...
5
votes
0answers
238 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 ...
5
votes
0answers
115 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, ...
5
votes
0answers
996 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 ...
4
votes
0answers
100 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: ...
3
votes
0answers
49 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 ...
3
votes
0answers
283 views

Simulation of surface tension-dominated interfaces

I want to simulate the shape of the free surface in a small fuel tank in microgravity, which is very slowly being emptied. The tank is not symmetric, the geometry is given by CAD (e.g. step file). ...
3
votes
0answers
102 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 ...
3
votes
0answers
55 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{...
3
votes
0answers
244 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 ...
3
votes
0answers
215 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$ ...
3
votes
1answer
428 views

General case Kutta condition

I'm working on a 2D inviscid fluid simulation using a "panel method", with Potential being used to enforce the no-through boundary condition. I'm trying to incorporate the Kutta condition, which says ...
2
votes
0answers
45 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 ...
2
votes
0answers
34 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 ...
2
votes
0answers
45 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 ...
2
votes
0answers
288 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 ...
2
votes
0answers
195 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 ...
2
votes
0answers
42 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 ...
2
votes
0answers
92 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 ...
2
votes
0answers
169 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 ...
2
votes
0answers
78 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 ...
1
vote
0answers
38 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 [...
1
vote
0answers
51 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 ...
1
vote
0answers
57 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 ...
1
vote
0answers
37 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 ...
1
vote
0answers
31 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 ...
1
vote
0answers
26 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 ...
1
vote
0answers
37 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 ...
1
vote
0answers
69 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 ...
1
vote
0answers
56 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 ...
1
vote
0answers
65 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 ...
1
vote
0answers
23 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 ...
1
vote
0answers
39 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 ...
1
vote
0answers
340 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 ...
1
vote
0answers
71 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 ...
1
vote
0answers
85 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 ...
1
vote
0answers
63 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=-\...
1
vote
0answers
918 views

BiCopter simulation in Matlab

I derived a bicopter dynamical model with two servos and two BLDC motors. And now am trying to simulate it using Matlab. As base for simulation I used this paper and this code Unfortunately, the ...
1
vote
0answers
55 views

Software to simulate behavioral economic models?

I am looking for a way to simulate behavioral economic models to illustrate the effects of changing parameters. I am familiar with statistical analysis (especially R), but I only have few programming ...
1
vote
0answers
113 views

Conservation at grid interface in adaptive mesh refinement

I am using adaptive mesh refinement to solve one dimensional inviscid Burgers equation. However I am facing some difficulty to handle grid interfaces which are not uniform (coarse-fine grid interface)....
1
vote
0answers
41 views

Spectra of the energy flux VS transfer spectra?

In direct numerical simulation of turbulence what is the advantage to look at the spectra of the energy flux rather than the transfer spectra? The transfer spectra corresponds to the triple-velocity ...
0
votes
0answers
52 views

Why is perfect sampling not used in large-scale lattice model simulations?

The statistical physics literature is replete with papers describing simulations of lattice models, such as the Ising model. Typically, these are done through Monte Carlo methods, such as the ...
0
votes
0answers
34 views

Effect of Peclet number on concentration profiles

I have been performing convection-diffusion transport studies on microfluidic channels like the following I came across an article that illustrates the effect fo Peclet number on concentration ...
0
votes
0answers
33 views

Numerical simulation for a bounded process. Is slight deviation a “normal” fact?

Suppose I have to numerically simulate a process $\{y_t\}$ such that $y_t\geq0$ $\forall t\in\mathbb{N}$, with $t$ denoting time-step. Let's suppose I use MonteCarlo with $\mathscr{N}$ simulation ...
0
votes
0answers
58 views

Stress during unloading in FEM simulation of an elasto-plastic truss element

I am simulating the behaviour of a single bar composed of two node points by fixing one point and applying a force on the oder node with $F_x = 1$ and fixing $u_y=0$. I am assuming the material to ...
0
votes
0answers
72 views

How one can simulate a system given by differential equation?

I want to simulate a diffusion environment given by the differential equation $$\frac{\partial u(x,y,t)}{\partial t}=D\left(\frac{\partial^2 u(x,y,t)}{\partial x^2}+\frac{\partial^2 u(x,y,t)}{\...
0
votes
0answers
48 views

Simulate friction arround steady-state point

When simulating a friction dynamic system, I have come to a problem when approaching steady-state. In the example below, a box moves along a surface with a certain friction coefficient: ...
0
votes
0answers
85 views

References on the topic of DEM and XDEM

DEM: discrete element method. XDEM: extended discrete element method. For my current project of furnace simulation with granular materials, I am interested in the methods mentioned above. I have not ...
0
votes
0answers
59 views

Energy anomolies in many body simulation

I am trying to simulate the gravitational interaction between many bodies. I am using a direct PP force calculation and a 4th order symplectic integrator with a variable step size. The energy of ...