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

learn more… | top users | synonyms

2
votes
1answer
47 views

Help understanding and implementing fast multipole method for N-body

I've been trying to understand the Fast Multipole Method but not really getting anywhere. It seems like the fastest mainstream N-body simulation algorithm, and I would like to implement it in an ...
3
votes
1answer
74 views

$(1+x^M)^{1/M}$ need to be able to calculate any order derivatives vs. $x$ and $M$ for $x\ge 0$ and $M\gt 2$

I need to be able to calculate any order derivatives of $(1+x^M)^{1/M}$ with respect to $x$ and $M$ for $x\ge 0$ and $M\gt 2$. Should I develop this into a Taylor series at a small $x$ with the order ...
0
votes
0answers
44 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 ...
6
votes
3answers
224 views

What do C, C++ and Java have that Fortran 2003 don't? [closed]

As you may know, Fortran is widely and longly (?) used in science. I've heard a lot of people saying that we natural scientists should all migrate to C++ or Java. I've seen people increasingly ...
2
votes
1answer
39 views

Is it common or to use objects for calculation of pair interactions in physics and chemistry simulations?

When simulation atoms, molecules, colloids, autc., are there programs that define each interacting unity as an object in the OOP sense? In my own case, I've been simulating magnetic nanoparticles in ...
2
votes
2answers
69 views

Is model-view-controller useful pattern useful to build scientific simulation programs?

I work with scientific simulation, but what I like the most in programming is creating graphical user interfaces and 2D or 3D outputs to visualize the evolution of the simulations. I know most input ...
1
vote
0answers
67 views

Numerical inversion involved confluent hypergeometric (1F1) (or Kummer function)

Edit: the question is solved ! Thank you for your time ! This problem arises when I tried to compute the valua of Asian call otions using Inverse Laplace transform method. Let $r=\mu = 0.15; \sigma ...
0
votes
0answers
15 views

Converting events with irregular arrival times into a continous signal (1st, 2nd and 3rd order)

I got events coming in at irregular intervals. From these I need to emit the number of events arrived and the expected first and second derivative over some specified measurement period P. So ...
0
votes
2answers
66 views

Simple java gravity simulation

I'm trying to write a class that uses Newton's law of Gravitation to work out the field of a planet. I've tested my code by inputting values for Earth, for purely vertical motion, so I should get g = ...
3
votes
1answer
71 views

How to model the mixing of two fluids in a container?

In the context of a research project of mine, I am faced with the difficult task of modeling the mixing of two fluids in a container. I would like to achieve the following: Given a container ...
2
votes
0answers
21 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 ...
3
votes
1answer
187 views

Finite difference scheme for 2D sound propagation

I am simulating the sound wave propagation in non-rectangular and asymetric spaces using finite-difference method. I presume linear acoustics equations to be enough (i.e. $\Box p = 0$, $\Box \vec{v} = ...
1
vote
2answers
172 views

Good desktop PC for molecular dynamics simulations

I want to build a new desktop PC for MD simulations of relatively simple protein/carbohydrate systems using mainly GROMACS and Amber force field. No ab initio or QM-MM simulations are intended for ...
0
votes
0answers
90 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 ...
0
votes
0answers
38 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 ...
0
votes
0answers
20 views

validating results of a hypersonic flow

I am trying to simulate the air flow behind a strong shockwave. The air mixture in nonequilibrium contains only 5 species(N2, O2, NO, O and N) where the three molecules are in vibrational ...
0
votes
0answers
71 views

Color Coding Gradient in Matlab

I am trying to visualize bond strength via a color gradient in matlab (such that a high energy would be in red, low in blue). The underlying data structure, bond locations, and bond energies are ...
3
votes
1answer
61 views

Solve implicit ODE numerically in orbit simulation

I'm trying to plot the orbit of a compact binary star system where general relativistic effects become important. I'm using post-Newtonian approximation and I want to solve the orbit numerically based ...
0
votes
1answer
112 views

How to calculate divergence and vorticity from a velocity field using finite elements

I am in the process of writing a finite element solver for the Navier-Stokes equations and am having trouble computing things like the divergence and vorticity correctly. Currently to compute the ...
5
votes
1answer
128 views

Implementing the pressure correction method using finite elements

Ok so I am nearing the completion of my finite element Navier-Stokes solver that uses the $\theta$-method for time stepping and the pressure correction method for the pressure. I am following the ...
1
vote
0answers
58 views

ANSYS Fluent: Defining External Force on a Cell Zone [closed]

I am using ANSYS Fluent 15 to simulate a case of fluid flow. In my case, I have three cell zones in my simulation, let's call them top, middle and bottom. My requirement is to put a vertical constant ...
2
votes
1answer
110 views

Pressure projection method boundary conditions

When using the pressure projection method to solve the incompressible Navier-Stokes equations do we apply Neumann boundary conditions for pressure only where there are associated no-slip velocity ...
2
votes
1answer
75 views

CFD turbulence modelling mean pressures vs peak pressures

I have been experimenting with using Autodesk CFD to investigate facade/ cladding pressures on a (rectangular) building, comparing results with cladding/ facades pressures pressures from design codes ...
-1
votes
1answer
224 views

How to initiate spirals in this model?

I am trying to reproduce Tang & Othmer paper which is related to excitations and oscillations in G-protein model in Dictyostelium discoideum, an amoeba species. The mathematical model in the paper ...
4
votes
2answers
108 views

Can other running processes harm the performance of my simulation?

I don't have an high-end machine on which run my simulation (Python + C extensions, based on numpy for calculations) Thus I intend to run the program on my PC, for several hours. I won't use the ...
0
votes
0answers
15 views

How to do numerical simulation of Impulsive Differential system?

I am working on a research paper on Plant-Pest-Natural enemy model. I want to plot graphs for my system as given below. I know how to simulate differential system using ODE45 command and then using ...
0
votes
2answers
116 views

Relationship between FEM solutions of PDE with different spatial resolutions

I use FEM to simulate deformations of elastic objects for animation applications in computer graphics. The governing equation is generally with the form: $$ \mathbf{M}\ddot{\mathbf{u}} ...
0
votes
0answers
51 views

modify COMSOL codes

I am a physics student and I am planning to simulate microfluidics and capillary waves on a fluid film and I decided to use COMSOL. I am not so much aware of COMSOL abilities. I wonder if I can ...
1
vote
0answers
56 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 ...
2
votes
2answers
209 views

GPU-enabled Lattice Boltzmann solvers?

Is anybody aware of any GPU-enabled Lattice Boltzmann solvers (preferably on C++/OpenCL and open-source) that would be recommended? I have found Advanced Simulation Library, but it seems to be very ...
0
votes
1answer
73 views

How to generate movie of the simulation done in LAMMPS? [closed]

I simulated nanowire deformation using LAMMPS and visualized it using VMD in Centos 6.5. But I'm not able to generate .jpg image file or the simulation movie. How can I proceed??
0
votes
0answers
58 views

How come in distributed computing 10X 32GB systems doesn't equal 1 320GB system

is there a mathematical formula to relate how much memory per node is lost when you use parallel computing? for example, when running COMSOL (FEA, MUMPS) in parallel, a problem that would require ...
3
votes
2answers
135 views

Should a colorbar have constant values throughout a simulation?

The default option for matplotlib's pcolormesh (also true for VisIt and I think also true for Matlab) is to allow the maximum and minimum of a colorbar vary. Which ...
0
votes
0answers
90 views

Modelling a Quantum Cascade Laser (QCL - GaAs atop Copper Heat Sink) with COMSOL

I am trying to model this: GaAs on-top of a heat-sink made of Copper. It's trying to simulate a type of Quantum Cascade Laser, I'm looking at the effect the temperature has by running a 10 V ...
1
vote
1answer
176 views

Hypersonic flow in Fluent

I am working in my graduating project about the hypersonic flow, particularly the properties of the air behind a strong shock wave, I try to do a simulation for that process in Fluent but I don't know ...
0
votes
0answers
31 views

Book recommendation on simulation of linear and nonlinear electronic circuits

I am looking for a book on electronic circuit simulation explaining the numerical modeling using transient and frequency based methods from a computational point of view. It should also have an ...
1
vote
2answers
130 views

Simulating a Simple Pendulum - Increasing amplitude on each swing?

I am trying to animate a dampened pendulum using RK4 for a highschool project. The equations that describe the dampened system are as follows: (from ...
3
votes
1answer
105 views

Surface charge boundary conditon for Poisson-Boltzmann solver

I want to set up a surface charge boundary condition for the simulation of semiconductor and electrolyte interface. The 2D-Poisson example and surface charge boundary condition are shown in following ...
1
vote
2answers
136 views

How to interpret/convert density in 2D simulations?

Density has units of $kg/m^3$. In simulations often scaling is done to deal with non-dimensional values, so the physical density is converted to the non-dimensional density so that it can be used in ...
1
vote
1answer
50 views

Value of density when there are no or very few neighbours in SPH simulation

Sorry for the noob question. I am trying to implement SPH using the directions shown in this paper. The density needs to be updated using the formula The smoothing kernel is If there are no ...
3
votes
0answers
72 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 ...
3
votes
0answers
53 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 ...
2
votes
1answer
113 views

Poisson solver diverged

My 2D-Poisson solver is build for simulation of semiconductor. The algorithm is Gauss-Seidel iterative method. If I use simple PN junction for simulation testing, it diverged while I applied high ...
4
votes
3answers
498 views

What are the differences between CFD simulations and realistic ocean/atmosphere model simulations?

The field of computational fluid dynamics (CFD) is dedicated to solving the Navier-Stokes equations (or some simplification of them). A subset of CFD, ocean and atmospheric models numerically solve ...
2
votes
2answers
465 views

Using Trilinos and/or Petsc in Windows with Visual Studio?

I am currently working on a project of building a simulator. Up till now I have mostly played around with very simple examples, but with regards to the future I have been thinking that it might be a ...
1
vote
1answer
98 views

Confusion in modeling and computing [closed]

I come from a computer graphics background without any professional in computational science. But I find my job(game engine development) being related to this field. May I ask some general questions ...
3
votes
1answer
406 views

Numerically solve a PDE in Python with a term calculated by coarse-graining

I'm trying to solve a PDE in Python of the form, $\dfrac{\partial c(\mathbf{x}, t)}{\partial t} = \mathrm{D} \nabla^2 c(\mathbf{x}, t) -\gamma \rho(\mathbf{x}, t) c(\mathbf{x}, t)$ where $c$ ...
2
votes
2answers
97 views

Simulating interaction of Smart Fluids with Solids

Is it possible to simulate interaction of Smart Fluids with Solids? Is there a software capable of doing so? What I Know: The two software that I know of are Autodesk Simulation Mechanical and ...
5
votes
1answer
195 views

Lagrange Multipliers in Multi-body Finite Element Code

I'm working on a Multibody dynamics code using the finite element method to simulate the behaviour of flexible beams (using this paper if anyone is interested/ it is relevant). I'd like to model ...
2
votes
2answers
183 views

Which processor should I consider buying for a desktop built for scientific computing?

I run Python/C simulations for biological problems, so the range of computational tasks can be broad. Here are some examples: 1) continuum model simulations: FEM type 2) agent-based simulations: ...