Questions tagged [differential-geometry]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
4
votes
3answers
507 views

Trying to define the Lie bracket of two vector fields

I am trying to define in the simplest possible way (only one coordinate system, no checking that variables are vectors, etc.) the Lie bracket of two vector fields in 3-space. What is wrong with the ...
3
votes
1answer
100 views

Solving a System of Differential Equations for Pendant Drop Application [duplicate]

So I'm trying to solve the system of differential equations describing a pendant drop. The system is as follows: ...
0
votes
1answer
49 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
2answers
124 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 ...
8
votes
2answers
227 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 ...
0
votes
0answers
73 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 ...
9
votes
2answers
306 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}$$ ...
1
vote
1answer
35 views

Partial differentiation second order

I have some rules for differentiation: ...
3
votes
1answer
135 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
1answer
91 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 ...
5
votes
1answer
771 views

Expand wedge product

How can I force mathematica to expand for example this expression $$(\cos (\theta ) dr-r d\theta \sin (\theta ))\wedge (\sin (\theta ) dr+r d\theta \cos (\theta ))$$ into what is should be, that ...
45
votes
6answers
34k views

How to calculate scalar curvature, Ricci tensor and Christoffel symbols in Mathematica?

I am seeking a convenient and effective way to calculate such geometric quantities. I've used packages like TensoriaCalc, but they don't work at all time. Sometimes,...
3
votes
2answers
156 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 ...
12
votes
3answers
5k views

Computing Gaussian curvature

Can Gaussian curvature $K$ be computed from WolframAlpha or any other available Mathematica program? Please indicate the program or its reference. If input ...
-1
votes
1answer
163 views

Can I numerically solve these equations in Mathematica? [closed]

I have this couple of equations : $ \partial_\mu \partial^\mu z^i + G^{i\bar{p}} (\partial_j G_{k\bar{p}} ) \partial_\mu z^j \partial^\mu z^k + G^{i\bar{j}} (\partial_{\bar{j}} G_{k\bar{l}} ) \...
4
votes
1answer
93 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
1answer
151 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 ...
6
votes
1answer
186 views

Comparing unit normal definition in calculus with FrenetSerretSystem

I'm trying to compare the unit normal definition in calculus texts (i.e., $\vec N=\vec T'/\|\vec T'\|$) where $\vec T$ is the unit tangent vector, with the unit normal vector returned by ...
0
votes
0answers
47 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: ...
2
votes
1answer
266 views

Principal Curvature of a monkey saddle using ParametricPlot3D but graph is showing blank

Edit: I'm working through a textbook by Alfred Gray, Modern Differential Geometry of Curves and Surfaces with Mathematica. I'm trying to plot the principal curvature of a monkey saddle, which I've ...
3
votes
1answer
352 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 ...
2
votes
1answer
166 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
4answers
331 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
3answers
335 views

How can I use procedure, circle3D, to make an animation of the osculating circle of a parametric curve?

I try to use circle3D procedure, shown below, to make an animation of the osculating circle of a parametric curve. ...
1
vote
0answers
60 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
1answer
233 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? ...
4
votes
1answer
354 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
0answers
45 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\...
1
vote
3answers
136 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: ...
18
votes
3answers
415 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 ...
8
votes
3answers
358 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 ...
1
vote
1answer
101 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 ...
1
vote
0answers
73 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 ...
8
votes
1answer
277 views

How to calculate specific area on surface of sphere?

I am trying to calculate the solid angle subtended by arbitrary-shaped loops on a sphere's surface. First, I parametrize circular loops by: $$\theta(t,k_{x0},r) = k_{x_0} + r \cos(t);$$ $$\phi(t,k_{...
1
vote
0answers
98 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
2answers
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 ...
2
votes
1answer
118 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
1answer
368 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
0answers
107 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
1answer
487 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 <...
5
votes
2answers
392 views

Can Mathematica solve nonlinear, coupled differential equations?

I've got two equations that describe a Geodesic on a sphere. $$ \frac{\mathrm d^2 u}{\mathrm d\lambda^2} - \cos u \sin u \frac{\mathrm dv}{\mathrm d\lambda} \frac{\mathrm dv}{\mathrm d\lambda} = 0 \\ \...
2
votes
1answer
168 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\...
4
votes
3answers
433 views

How to calculate arbitrary area on surface of sphere?

I am trying to calculate the solid angle subtended by arbitrary-shaped loops on a sphere's surface. First, I parametrize circular loops by: $$\theta(t,k_{x0},r) = k_{x_0} + r \cos(t);$$ $$\phi(t,k_{...
5
votes
1answer
229 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+...
5
votes
2answers
259 views

Finding inflection points of 2D BSplineFunction

This code creates a jagged line in two dimensions and fits a BSpline function to it. ...
5
votes
3answers
495 views

Calculating curvature of a contour

I have an equation of a scalar field in the form $$f(x, y) = x^2 + y^2 + xy + c$$ I want to find the curvature of the contour of the curve at $f_c = f(0.5, 0.5)$. So I need to calculate the ...
1
vote
1answer
122 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
1answer
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
3answers
573 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 ...
0
votes
0answers
80 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 ...