I don't mean to use differential equations or things like that, just
apply node laws or Kirchhoff's laws and thus be able to calculate the
values of the electrical components
The first step in analyzing this circuit is to simplify. R1 and R2 can be combined into a single resistor. C2 is small enough that at 50~60 Hz it will have negligible effect. D1..D4 can be replaced with a rectified sine wave and single blocking diode. We will assume that D5 and Q1 B-E junction have a voltage drop of ~0.7 V, and that the optocoupler LED has a voltage drop of 1.1 V. The simplified circuit then looks like this:-
simulate this circuit – Schematic created using CircuitLab
Now it is easier to see what role each component plays and what effect its value has on the zero-cross detection timing.
R1+R2 and R3 form a voltage divider that reduces the rectified mains voltage, and R1+R2 in parallel with R3 is the effective series resistance that determines how fast C1 charges (use Thévenin's theorem to derive the equivalent voltage and resistance).
When the mains voltage is close to zero the rectified input voltage drops below the voltage on C2, causing Q1 to turn on and switch R4 and LED1 across the capacitor which discharges it faster, until the input voltage rises at the end of zero crossing and turns Q1 off.
If Q1 never turned on the voltage on C1 would center around the average voltage of the rectified input waveform (2Vpk/π), less 1 diode voltage drop through D5/Q1. However with Q1 turning on during zero crossing the average capacitor voltage drops due to faster discharge through R4 and LED1. Capacitor voltage is needed to determine the value of R4 required for the desired LED current, but the calculation is complicated by interaction with the pulse ratio (and you wanted to avoid differential equations).
While you could develop a formula that includes all the required parameters, it will probably be a lot easier and less error-prone to simulate it using a program such as LTspice.