Questions tagged [differential-geometry]
The differential-geometry tag has no usage guidance.
132
questions
0
votes
1
answer
61
views
Covariant Derivative
I'm trying to create a function which receives as inputs a coordinate system $x^\mu$, metric tensor $g_{\mu\nu}$, a general tensor $T$ of rank $n$, and a list $I$ of length $n$ that looks like this:
$$...
0
votes
2
answers
117
views
Is it how you do a projection? Is it a Projection? I don’t know what this is
I don't really know, if this is a projection. Maybe you can tell me. . . From the parametric equations
...
1
vote
1
answer
79
views
How to compute the divergence of a four-vector?
I have a quadri-vector which is given by
u = {(E^(-φ0[r]))*(1 - ε δφ[t,
r]), (E^(-φ0[r])) D[ε ξ[t, r], t], 0,
0}
and a quantity n which is given by
<...
0
votes
1
answer
102
views
How to solve the einstein field equation symbolically? [closed]
How can I use Mathematica to solve the disturbed Einstein field equation? Is there any notebook that introduces this subject or package like diffgeo?
0
votes
0
answers
50
views
How to compute the covariant derivative of a vector? [duplicate]
I want to compute a derivative of a vector, but I'm new to mathematica (I programmed in other languages). I was advised to work with the diffgeo.m package.
What I want to do is calculate
using ...
1
vote
2
answers
269
views
How do I get the Schwarzschild solution from the Einstein Equations?
Consider the Schwartzchild-like metric $$ds^2=-A(r)dt^2+B(r)dr^2+r^2(d\theta^2+\sin^2\theta d\phi^2).$$ The Einstein field equations for this metric reduce to $$R_{\mu\nu}=0,$$ which is also known as ...
0
votes
1
answer
305
views
How to do Einstein Summation in Mathematica?
I am working with the de Sitter metric which takes the form
$$g=\tau^{-2}\left(-\left(\Lambda / 3-\tau^{2}\right)^{-1} d \tau^{2}+\left(\Lambda / 3-\tau^{2}\right) d t^{2}+g_{\mathbb{S}^{2}}\right)$$
...
1
vote
3
answers
104
views
What is the most efficient way to turn a metric formula into a metric tensor?
I have a metric formula:
ds=(-dt^2)*(c3+a3*t)^2+(dxC^2*(t0^2+t1^2)^2)/(4*t0^4)+
(dxM^2*(t0^2+t1^2)^2)/(4*t0^4)+(dxY^2*(t0^2+t1^2)^2)/(4*t0^4)
How do I turn this ...
0
votes
0
answers
97
views
Christoffel symbol of the first kind [duplicate]
Suppose that we are given a metric $$ds^2=-\left(1-\frac{r_s}{r}\right)c^2dt^2+\left(1-\frac{r_s}{r}\right)^{-1}dr^2+r^2d \theta^2+r^2\sin^2(\theta)d \phi^2.$$ Given the Catalogs of Spacetimes pdf we ...
0
votes
0
answers
52
views
Finding transformation matrix for a coordinate transform
I have a matrix in two different bases. Suppose they are called $g_{\mu \nu}$ in one basis and $\eta_{ab}$ in the other basis. If the transformation between the bases is represented by $T_{\mu}^{\text{...
2
votes
1
answer
88
views
Problem with ArcCurvature
When comparing the outputs of
...
8
votes
0
answers
212
views
Shortest Distance between two points on a 2D surface
I am working on an interesting problem, that consists of trying to compute the shortest distance between two points on a 2D surface. This 2D surface can be represented in a 3D space by the function $f(...
0
votes
1
answer
54
views
Define Matrix Function in a For loop
Let $g\in\mathcal R^{3\times 3}$ be a given known rank 2 tensor function, $gInv$ its inverse, and $dg\in\mathcal R^{3\times 3\times 3}$ its gradiant:
...
1
vote
2
answers
152
views
Plot or draw a differential geometry graph in Mathematica?
I want to draw a differential geometry diagram, but I don’t know how to draw it. I can only draw graphics like the first picture. There are many differential geometry diagrams in textbooks. What tools ...
0
votes
0
answers
84
views
Computations on differential geometry
I have to do some really long computations of basic riemannian geometry. I am trying to use Mathematica to check them. This is the kind of computations I want to do. Let $M$ be a riemannian manifold ...
1
vote
1
answer
37
views
Partial differentiation second order
I have some rules for differentiation:
...
8
votes
2
answers
248
views
What is the best source to learn how to use tensor operations (exterior algebra) in Mathematica?
I'm specifically interested in the TensorProduct,TensorWedge, HodgeDual and certain build in functions to do tensor arithmetic like TensorReduce, TensorExpand.
I would like to do exterior algebra ...
1
vote
1
answer
120
views
Einstein field equations for Bondi-Sachs formalism
I'm trying to re-derive the results of Bondi-Sachs formalism. The metric is given in the form
\begin{array}{c}g_{a b} d x^{a} d x^{b}=-\frac{V}{r} e^{2 \beta} d u^{2}-2 e^{2 \beta} d u d r+r^{2} h_{A ...
3
votes
2
answers
180
views
Calculate surface normals at the boundary of a Graphics3D object
How do I go about calculating and plotting the surface normals at the boundary of a Graphics3D object?
For example, consider this custom-defined ...
4
votes
1
answer
104
views
Differentiation by indexed variable in equation of Christoffel Symbols
I am very new to Mathematica, and am trying to use it to compute Christoffel symbols for a certain manifold. All this requires is taking some indexed sums of derivatives, but it builds up on several ...
5
votes
1
answer
193
views
Finding unit tangent, normal, and binormal vectors for interpolated function
As an extension of this question, is it possible to find the unit tangent, normal, and binormal vectors for an interpolated function? eg
...
9
votes
2
answers
319
views
Prove $R^\top R = I_3$ and find skew-symmetric matrix $R^\top \dot R$
Consider the following operator defined over unit vectors of $\mathbb{R}^3$:
$$R(u,v) = (u\cdot v)I_3 + hat(u\times v) + \dfrac{(u\times v)\otimes (u\times v)}{1+u\cdot v}$$
...
0
votes
0
answers
54
views
How do I interpret this table of Christoffel symbols?
So I found a code that allows me to compute the covariant derivative of some vector, here it is:
...
3
votes
1
answer
180
views
Covariant derivative of a vector [duplicate]
I have a code that gives me the Christoffel symbols of a metric. How do I take the covariant derivative of a vector?
It does not necessarily have to build upon my code, but this is what I have used so ...
1
vote
0
answers
76
views
Product of manifolds with non−zero non−diagonal boxes in the metric
I'm trying to construct in xAct a metric like this
where $g_{\mu\nu}$ is 4-dimentional, and $g_{MN}$ - 5-dimentional, $A_\mu$ - 4-vector and $\phi$ is a scalar field.
I already tried to do it like ...
2
votes
1
answer
275
views
How can I draw a sphere in the Minkowski space?
If the equation of the circle in the Minkowski 3 space is given as
$S_1^2 = \{x \in E_1^3:- x_1^2 + x_2^2 + x_3^2 \}$, how can I replace it in the following code?
...
0
votes
0
answers
50
views
Parametrization of twisted pseudospheres types 2 and 3
Among the three types of rotationally symmetric Pseudospherical surfaces $K=-1$ (Beltrami central, type 2 and type 3)
(http://xahlee.info/surface/gallery.html)
we have Dini twist addition term $ b\...
2
votes
1
answer
209
views
A doubt on ParametricPlot3D, RevolutionPlot3D, ListPlots and NIntegrate: can I build an "RevolutionListPlot3D"?
First of all: this is question lies in the context of Surfaces and Embbedings on differential geometry. More precisely in the context of Kruskal coordinates and how to plot a 3D dynamical ...
7
votes
4
answers
349
views
Why do isolated large values of WorkingPrecision fail in NDSolve?
Executive Summary:
Getting an accurate answer often depends on setting the WorkingPrecision high enough. Once it is high enough, though, I would expect that its ...
1
vote
0
answers
84
views
Problem with the analytical study of transient processes in nonlinear systems with linear dynamic links
I ask for advice and help.
I am having difficulties of this nature. There is a nonlinear system of the following type:
I need to analyze analytically the transient process in such a system. The ...
1
vote
1
answer
115
views
DSolve does not work
I am trying to solve this coupled nonlinear pdes for $\kappa(x,t)$ and $\tau(x,t)$:
where $\zeta_1 = \kappa(x,t)$, $\zeta_2 = 0$. I used this code
...
18
votes
3
answers
441
views
Coordinate-free derivative
Given the function
\begin{align*}
f \colon \mathbb{R}^n &\to \mathbb{R}^n\\
v&\mapsto \dfrac{v}{\|v\|},
\end{align*}
I would like to compute the derivative of $f$, that is $df(v)$. It is ...
1
vote
3
answers
152
views
Cannot derive `Norm` or `Normalize` when recreating Frenet Serret equations
I'm trying to calculate the torsion of a curve at a point using the following code:
...
1
vote
0
answers
105
views
How Can I draw a curve on a sphere in Minkowski 3 space?
There is a curve, which evolves with time, in the the Euclidean space. And the solution of its evolution equations has been used to draw it using this code:
...
8
votes
3
answers
388
views
How is Grad defined for array particularly in non-Cartesian coordinates?
This question can be viewed as a follow-up of
What is the definition of Curl in Mathematica?
First argument of Grad can be an array, but what definition does ...
2
votes
1
answer
142
views
How to identify specific Christoffel symbols and Riemann Tensor components from a general solution [closed]
I want to learn how to identify specific Christoffel Symbols and Riemman Tensor components from the general solution provided by Mathematica. Let us work out an example to see what I mean clearly.
...
4
votes
1
answer
488
views
Computing Christoffel symbols of the second kind [duplicate]
I want to compute the Christoffel-symbol for a given metric. I am using the code here, but I am missing something.
The Chrisfoffel-symbol formula is
$\Gamma^{\mu}_{\phantom{\mu}\nu\sigma}=\frac{1}{2}g^...
1
vote
0
answers
113
views
Computing wedge product of vector fields
This is a question of actually linear algebra. Say I have a vector spaces spanned on x3,x4,..., x9 I am trying to find wedge product of two elements of this vector space. My numbers a3, ..., a9 are ...
1
vote
1
answer
614
views
Parametric Plot 3D: A curve and a surface plot together
I am new to Mathematica and
I tried to make a superimposed plot of a curve on a half sphere and need help with some graphics.
This is the code line that i used to generate the following graphs
<...
2
votes
1
answer
182
views
How can I create this graph in Mathematica?
I'm interested in illustrating a hyperbolic plane for a report I'm writing. Here's the metric that I'm using:
$$ds^2=dr^2+\sqrt{\lvert k\rvert}^{-1}\sinh\left(r\sqrt{\lvert k\rvert}\right)\left( d\...
1
vote
1
answer
125
views
Frenet frame in Pseudo Galilean Space [duplicate]
If I have the Frenet frames for an admissible curve say $\alpha$ in a pseudo Galilean space, which is given by: $ t'(x) = \kappa(x) n(x)$, $n'(x) = \tau(x) b(x)$ and $b'(x) = tau(x) n(x)$, where tau ...
2
votes
1
answer
66
views
I need help solving this hyperbolic equation [duplicate]
I have some data and I'd like to calculate the radius of curvature. The formula is:
$$R_{oc}\space Sinh\left[\frac{D_{LSS}}{R_{OC}}\right]=\frac{s_*}{\theta_*}$$
Noting that $s_*$ is sh, $\theta_*$ ...
4
votes
1
answer
408
views
Find geodesics given two points
Suppose I want to find the geodesic on a paraboloid going through two points located on the paraboloid, using Mathematica 12.
I chose the paraboloid parametrization as follows:
...
0
votes
0
answers
82
views
How to compute trajectories normal to field lines?
I have function [Psi][r,z] found from solution of a Grad-Shafranov equation. Magnetic field is expressed through [Psi][r,z] as follows
...
1
vote
0
answers
102
views
Dirac Operator in de Sitter background
I am attempting to write the Dirac operator in a curved background, and eventually solve the equation as a second order PDE, since I am attempting to bring it into Klein-Gordon form. Essentially what ...
3
votes
1
answer
521
views
Difficulties on Mathematica code to solve Christoffel Symbols of a particular metric
I) The Problem
There's a particular metric $[1],[2]$ in general relativity which is written as:
$$ds^{2} = -[c^2-v_{s}^2f(r_{s})^2]dt^2+v_{s}f(r_{s})dtdx+v_{s}f(r_{s})dxdt+ dy^2+dz^2 \tag{1}$$
So ...
5
votes
1
answer
555
views
Fourth Order Tensor rotation [duplicate]
What is the easiest way to perform Rotation for Higher Order Tensors in Mathematica ? For Instance 4th order tensor
$C_{ijkl} = \lambda_{im}\lambda_{jn}\lambda_{ko}\lambda_{lp} C_{mnop}$
4
votes
3
answers
634
views
It seems Eigensystem[m] returns vectors that are not eigenvectors
I am new to here so please forgive me if I do something wrong carelessly. I have faced a serious problem in eigensystem method, or more particular, eigenvalue. It seems that the following codes that ...
5
votes
1
answer
236
views
Archimedean spiral from curvature
I am trying to reconstruct an Archimedean spiral from its curvature
$$\kappa (\text{s$\_$})\text{:=}\frac{s^2+2}{\left(s^2+1\right)^{3/2}};$$
eqns:
$$\left(
\begin{array}{c}
t'(s)=\frac{\left(s^2+...
8
votes
2
answers
1k
views
The time-like geodesics (orbits) in the Schwarzschild spacetime
I am trying to plot Schwarzschild's orbit without invoking the geodesic equation.
As a reference I am using Chandrasekhar's Book (The Mathematical theory of Black Holes, Oxford University Press). In ...