Questions tagged [computational-physics]

Computational physics is the study and implementation of numerical algorithms to solve problems in physics for which a quantitative theory already exists.

Filter by
Sorted by
Tagged with
0
votes
1answer
81 views

Mechanics - Rotational Inertia

I am facing a problem solving a mechanical system which is shown in the scheme below. We assume that coupling of the two shafts is done with a clutch Equations: $$J1*\frac{d\omega1}{dt} + Td(\...
0
votes
1answer
119 views

Solving differential equations in MAPLE (ballistic curve problem) [closed]

TASK: Make a simple animation of an object shot from a catapult. Object is moving through a medium. The problem must be solved in MAPLE. NARROWING OF THE PROBLEM: 2D problem consider atmospheric ...
0
votes
2answers
3k views

Most efficient library to diagonalize exactly large hermitian or unitary matrices

I am working on a physics problem which requires obtaining the exact eigenvalues and eigenvectors of Hermitian and Unitary matrices numerically. Naturally I would like to ask the experts what are the ...
0
votes
1answer
65 views

How to model waterflow when only a couple of sample points available

Figure below depicts a cross section of a creek for which I am trying to measure the water flow for that section. What we have as inputs are a bunch of sample points on the river. For each sample ...
0
votes
0answers
4 views

Propagation of a Gaussian beam using FFT

I am trying to simulate the propagation of a gaussian beam through a lens using an FFT approach. I tried to implement the approach described by Couairon in this paper at page 43: https://link....
0
votes
0answers
49 views

Is it possble to do this complex symbolic calculation with Matlab?

Sorry it's bit abrupt, but recently I am caught up in some symbolic calcualtion which is tedious and almost impossible with mere human hands, so just wondering is it possible to solve the double ...
0
votes
0answers
42 views

Troullier-Martins pseudopotential with nonzero angular momentum quantum number

According to [1] the Troullier-Martins pseudopotential for quantum number $l$ is computed by $$ V^{\textrm{KB}}_{\textrm{nonlocal},l}(r) = \frac{\vert{}V_{\textrm{nonlocal},l}(r)\Phi^{\textrm{PP},0}_l(...
0
votes
0answers
137 views

Ising model simulation offset critical temperature and interal ernergy

I'm writing a code for the Ising model using WHAM (the weighted histogram analysis method),But it seems to produce critical temperature and internal energy wrong. (newest rewritten code is below) <...
0
votes
0answers
29 views

solving electron motion in undulator by Boris method

I am trying to use Boris scheme to solve the electron trajectory in undulator. The undulator field I used is: $$B_x = b_0\sin(2\pi \tfrac{z}{\lambda_u})$$ where $b_0 = \dfrac{2\pi c_{0}K}{q m_{e} \...
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
30 views

A preconditioner for self-consistent iteration

I tried to derive a preconditioner for self-consistent iteration similar to section IX in arXiv:0804.2583. For simplicity, consider here only one orbital (one or two electrons) systems. Suppose that ...
0
votes
0answers
65 views

Computing a double Integral using two Riemann sums & graphing multiple isosurfaces

I'm trying to compute the potential and electric field of a uniformly charged spherical shell and plot the results in the space outside the shell using MATLAB. I've tried everything from for loops to ...
0
votes
0answers
27 views

How can I implement the Invaded Cluster Algorithm for a network of Ising spins?

My primary concern is about finding the percolating cluster for any given network. For a lattice it is straight forward : When the size of a cluster reaches the length of the lattice, then it is said ...
0
votes
0answers
13 views

Compute cell values and isosurface constant from density values of particles

I am trying to reconstruct the surface for a fluid simulation based on a list of particles using the Marching Cubes algorithm. From different resources, such as http://paulbourke.net/geometry/...
0
votes
0answers
93 views

Double mach reflection at a inclined wedge

I am running into a strange problem when solving the 2D compressible Euler equations on a inclined wedge. To elaborate, my top boundary condition seems to emitting some type of instability. I have ...
0
votes
0answers
40 views

Time sampling changes solution

I'm currently trying to solve a problem using numerical methods. The set-up is rather long, so I apologize in advance... TL;DR: My solutions change depending on how big my steps are and I don't know ...
0
votes
0answers
32 views

3D Tollmien-Schlichting Waves Imposed in a Channel Flow (Are Physics correct?, etc)

So I am trying to do some further tests on a 2nd-order code Incompressible Navier Stokes equations, by studying transition to turbulence in a Poiseuille flow. Specifically, I'm interested to see ...
0
votes
0answers
32 views

Compute a Boltzman partition function

I'm trying to calculate total energy of a system $$ E(v, h) = -\sum a_iv_i - \sum b_jh_j - \sum_{i,j} v_ih_jw_{ij} $$ Python equivalent looks like this ...
0
votes
0answers
49 views

Determining the pseudo-time period of a system of $n$-pendulums via Kane's method in Python

We can use Kane's method to integrate the equations of motion for a system of $n$ pendulums with arbitrary masses and lengths (see derivation). In particular, if $(x_i,y_i)$ denotes the Cartesian ...
0
votes
0answers
53 views

Does radioisotope decay affect computer systems?

It is known that there exist a number of radioisotopes of elements commonly used in computer systems. Is the decay of these materials known to affect device performance over time, or is its impact so ...
0
votes
2answers
63 views

Physical meaning behind the choice of the proposal distribution in Markov Chain Monte Carlo (MCMC) methods

Let us consider the conventions on names used in the theoretical derivation of Metropolis-Hastings Monte Carlo as outlined here, for the sake of common nomenclature. What we are building is a step-...
0
votes
0answers
107 views

How to keep velocities in check in molecular dynamics simulation?

I am trying to make a very simple molecular dynamics simulator with Reflective Boundary Conditions. I am assigning the initial positions in a cube randomly while making sure they are not too close to ...
-1
votes
0answers
56 views

Simulating an object in orbit

This question is more oriented around suggestions for simulation tools and how to approach simulating an object in orbit. So high level I am trying to simulate the concept of a Sky Hook. What are the ...
-1
votes
0answers
23 views

non premixed Combustion by using oxygen droplets

I am trying to simulate a case in ANSYS FLUENT for combustion of methane and oxygen where oxygen enters the domain in the form of droplets. There is a change in temperature but there is no formation ...
-1
votes
1answer
45 views

Smoothed Particle Hydrodynamics: Weird clustering of particles. Is that normal?

I implemented a rather simple SPH simulation using a cubic-spline-kernel and a simple non-iterative pressure solver as described in this PDF in equation 9. I followed algorithm 1 of that paper (...
-1
votes
1answer
57 views

Simulating magnetic particles in a field free point generated by two opposing magnets

This is probably a long shot with such a short time, but I've been trying to get theoretical data for a project I'm working on. The project involves using a very simplified version of magnetic ...
-2
votes
1answer
246 views

Solving Schrodinger equation numerically

Here is the Schrodinger equation that is to be solved: A 1D hard wall potential in $[0, 1]$. The potential within the potential well is given by a linear combination of Gaussian dips $$v(x) = - \...