Questions related to the application of Mathematica to problems in physics.

learn more… | top users | synonyms

1
vote
1answer
175 views

My code uses `ClebschGordan` but `Mathematica` is using `ThreeJSymbol`

I'm using a function that calculate CG coefficients with the function ClebschGordan but instead, I've got the following error because Mathematica is using ...
3
votes
1answer
114 views

Solving 2D Laplace equation for irregular boundaries [on hold]

I want to solve Laplace equation in 2D with $x$ and $y$ coordinates with the boundary conditions $V(0,y)=-180$; $V(x,0)=50x-180$; $V(x,6)=50x-180$; $V(453.595-\sqrt{450.05^2 - (y - 3)^2},y)=0$ The ...
2
votes
3answers
791 views

draw part of parabola using graphics

today I was doing my physics homework. one of the questions was draw velocity-time and position-time graphs from the acceleration-time graph. So for the first step I draw acceleration-time using ...
0
votes
0answers
92 views

Animating ballistic movement

I was tasked with creating a 3D animation of the movement of a ball through space, while taking into account influences such as the magnus effect, air drag, and rotation. In order for this to happen a ...
13
votes
2answers
266 views

Heat convection differential equations from 1952 - Mathematica “fails to converge”

I am trying to solve a fundamental problem in analytical convective heat transfer: laminar free convection flow and heat transfer from a flat plate parallel to the direction of the generating body ...
6
votes
1answer
115 views

How to numerically solve a 1D time-independent Schrödinger equation for two interacting particles

Solving the 1D single-electron time-independent Schrödinger equation has been demonstrated using NDEigensystem here. There, the single-electron Schrödinger equation ...
2
votes
2answers
259 views

Ion motion in uniform axial magnetic field and radial electric field [closed]

I'm trying to model an ion's 3D motion in a constant/uniform axial magnetic field and a radial electric field (similar to a cylindrical magnetron, but two positive parallel center electrodes offset by ...
3
votes
1answer
153 views

Controlling scaling of multiple vector fields in VectorPlot3D

I'm trying to plot the electric and magnetic fields of a relativistically moving electric charge using VectorPlot3D, and apply "Manipulate" to change the velocity. The code I'm using is as follows: <...
8
votes
1answer
103 views

Calculating density of states from lists of momentum/energy data

I am trying to calculate the density of states for a 2D system given certain band structure data organized in a specific format. Below is an example of my procedure: Beginning with a simple 2 x 2 ...
9
votes
3answers
2k views

How to create an inverse pole figure color map

To interpret electron backscatter diffraction (EBSD) results, inverse pole figures are used. A three dimensional space of directions ([100], [110] and [111] for those who know about crystallography) ...
3
votes
1answer
119 views

How? 2-dimension visualisation of all 3d-orbitals as either a density plot or as function of maximal radius

I've been playing around with some d-orbitals and have been trying to view their maximum in two dimensions. At the moment I've progressed to the point where the d__ orbitals are superimposed onto one ...
0
votes
1answer
73 views

Lagrangian for Spring-Pendulum [closed]

I solved a Spring-Pendulum System with NDSolve. I tried to plot a graph with the interpolated data but always returns a message "An improperly formatted directive with head Symbol was encountered." <...
21
votes
4answers
3k views

Convert spectral distribution to RGB color

Presume that I have a spectrum as a function of wavelength (an example being the blackbody spectrum): I want to convert that to a single RGB color to display on-screen, i.e. the "color" of that ...
8
votes
2answers
325 views

Lagrangian to Hamiltonian

I want to go from Lagrangian description to Hamiltonian one. Using the example given by Mathematica I do something like: ...
31
votes
3answers
975 views

A better “VisibleSpectrum” function?

Since ColorData["VisibleSpectrum"] is wrong, I would like to have a more accurate function to use. Can this information be extracted from Mathematica itself?
1
vote
1answer
240 views

Why Integral does not converge?

I am trying to solve for $\Omega$ this nonlinear integral equation: $$1+\dfrac{z}{k^2}-\dfrac{z^2}{K_2(z)} \dfrac{\Omega}{k^3} \displaystyle\int_{1}^{\infty} \gamma^2\, \text{ArcTanh} \left(\sqrt{\...
2
votes
0answers
104 views

What is a good way to model fluid stream contours around airfoils?

I am an undergraduate student and am doing a project on conformal mapping for my complex analysis class. I have had one heck of a time trying to plot stream function contours to model fluid flow in ...
0
votes
0answers
51 views

Eigenelements of conductivity equation

I am trying to calculate the eigenvalues and eigenfunctions of the conductivity equation in an annulus. In particular I am looking for $(\lambda, u)$ s.t. $$ \begin{cases} \Delta u = \lambda u & \...
0
votes
0answers
34 views

Solving Time-Dependent PDE's Analytically with no non-trivial Boundary Conditions

I'm having some trouble solving a partial differential equation that is supposed to model the propagation of a scalar field in Schwarzschild spacetime using spherical coordinates. I'm using the DSolve ...
0
votes
0answers
38 views

Recycle boundary condition

I'm having the following issue that can not find a solution , by using this boundary condition u[0, t] == 2 u[1, t] that Mathematica does not understand, someone ...
0
votes
2answers
84 views

Jump boundary conditions

I'm having the following issue that can not find a solution , by using this boundary condition ...
0
votes
2answers
88 views

Free-fall air resistance problem

I´m working on a free fall problem and I got an equation with variable $k$ when the situation has air resistance. Then I compared with no air resistance equation ($-4.9$, considering t variable). So, ...
-1
votes
1answer
42 views

Solve returns empty list on system of 62 linear equations with undeclared known variables

Context: I am trying to solve a set of equations generated by hand from a complicated, but not particularly strange system via Newton's Second and Third Laws, and a few constraints. All the resulting ...
0
votes
1answer
59 views

Verify angle pi/4 maximizes launch range [closed]

I'm solving the equation of motion for a projectile in the x, y plane and got this: R = v0^2 sin(2 theta) / g where R is the range of how far the object will ...
0
votes
1answer
99 views

Plotting graphical analysis of chaotic behavior

I am trying to plot some results of chaotic behavior in a system based on a paper for one of my classes. I am a physics undergrad. I found some sample code online springy.nb which I have been trying ...
0
votes
3answers
173 views

Can I convert a 3-D plot into a 2-D plot?

I am studying the Landau Free Energy Equation. First of all the plot displayed only shown negative range. can i show the positive range too? (y-axes). Then, can I turn this plot into 2-D Plot? ...
-1
votes
1answer
57 views

NDSolve and vectors [closed]

I think I can't find solution of equation below because I need use scalars, not vectors as I have tried. But I don't know how I can do express the problem with scalars. I have read Wolfram Language ...
22
votes
1answer
288 views

How to model wooden joints with mathematica's FEM?

This is a dovetail joint: and I'd like to see the stresses and deformation on the joint.I haven't seen any modeling of disconnected regions with FEM, only connected regions, so I'm curious if you ...
12
votes
2answers
302 views

Sampling a phase space uniformly for a given energy?

The energy for a 1D harmonic oscillator can be written: $$ E = \frac{1}{2 m} p^2 + \frac{m \omega^2}{2} x^2 $$ where x is position and p is the momentum. I would like to sample phase space $(x,p)$ ...
1
vote
0answers
102 views

How to find the eigenvalues of a custom operator using NDEigensystem? [closed]

I would like to solve a Hamiltonian presented by a 2x2 matrix. Hamiltonian is given by: ...
0
votes
1answer
125 views

Animate spherical pendulum

I want to plot the position of Spherical pendulum. First I tried to plot the simple pendulum: ...
0
votes
1answer
124 views

How can I solve a Hamiltonian with numerical methods?

Hamiltonian in terms of two level atom's operator (0 is ground state and 1 is excited state) which are 2*2 matrices and cavity modes are given by $a$ and $a^\dagger$. Rest of the parameters are ...
0
votes
2answers
80 views

Determine the initial speed and launch angle of a projectile from its trajectory

I have the coordinates {x,y} needed to plot the trajectory of a projectile, which are given below: ...
6
votes
0answers
238 views

Problems with NDSolve for solving coupled Schrödinger-like PDEs with perfectly matched layer boundary condition

I'm trying to solve a full-vectorial wave equation for an arbitrarily shaped wave guide, by using NDSolve and perfectly matched layer (PML) conditions. The PML ...
2
votes
1answer
122 views

Improving the speed of evaluating numerical solution to a sytem of PDE

I am running a simulation for a kinetic reaction inside a pellet modeled by a system fo reaction-diffusion equations. ...
9
votes
1answer
148 views

Beam deformation due to localized force distribution

I would like to simulate the deformation of a beam which is fixed on one end and has some localized forces applied along the top. I have been inspired by this Mathematica example: http://www....
9
votes
3answers
675 views

Calculate inertia tensors

Moments of inertia and inertia tensors are well-known characteristics of rigid bodies in physics and applied mathematics. How to calculate them in Mathematica? Is possible to use regions which was ...
30
votes
2answers
3k views

The Orbit and Perigee of the Flamsteed comet

Historical context This year we have the 330-th anniversary of the Battle of Vienna - one of the great formative events of European history, it took place on September 12, 1683. Kara Mustafa, Grand ...
5
votes
2answers
272 views

How to see phase difference of two pendulum

i want to see the simulation of two uncoupled Simple Harmonic Motion (SHM) pendulum with different Lengths and released at different initial displacement angles. in this case adding one more pendulum ...
2
votes
2answers
89 views

How to make A simulation of 2-Dimensional Projectiles

Lets say two cannon is fired at different initial position, initial angle and initial speed. let say initial speed,angle and position for projectile 1= ...
0
votes
1answer
77 views

Making a Manipulate showing a the motion of a particle [closed]

Solving a differential equation I got to the following result -(((R[n*Pi - t/w] Cos[t/w])/(n*Pi)) Then, using this solution, I create . ...
1
vote
1answer
161 views

Boundary Problems

I am trying to study what happens to a wave as it passes different boundaries. This is the working wave equation defined on a plane ...
1
vote
1answer
188 views

Raman spectra collected at over 200 spatial points and in separate text files, how can I collate and visualize the data?

I need to do density plot and contour plot of Raman spectra with a specific peak, Raman mode, as a way to study the intensity distribution of such peak throughout the sample used. As such ...
9
votes
1answer
1k views

Schrödinger eigenvalue problem in two dimensions (Harmonic Oscillator)

I read here the discussion about how to solve a one-dimensional eigenvalue problem. I am wondering, how can one generalize these methods to two dimensions? For example, how to solve this equation: $$...
0
votes
2answers
241 views

A function generating particle-in-a-box eigenfunctions

This function sequence is defined over $\left[-\frac{L}{2},\frac{L}{2}\right]$ (all elements are valued 0 outside this interval): $$\psi_n(x)=\begin{cases}\sqrt{\frac2{L}}\sin{\left(\frac{n\pi x}{L}\...