Starting with a formula like (changing to American symbols because I'm a boorish American; N is turns, etc.),
$$V = 4.44 F N B A_e$$
(Incidentally, 4.44 is actually \$\sqrt{2} \pi\$, the peak area under a sine wave, amplitude given as RMS. The equivalent factor for a square wave is simply 4 (so using either factor for either waveform, doesn't actually make a huge difference), and note that square wave peak = RMS so we aren't measuring the square ambiguously if we leave off the suffix.)
We can substitute out N to arrive at a bulk formula for transformer design.
I'm not going to do this, actually; rather I'll do the definitional equivalent: normalize it to 1 turn, and eliminate V and I.
To do this, we need to know a few geometric factors for the core, and we need to know how much copper we're stuffing into the thing. Which means losses are a term (copper has resistance, so its amount and power dissipation matter). And then we should be mindful of core losses, but we don't have a formula handy for this. We can just guess, and assume it's alright. (The modified Steinmetz formula can in fact be used here, to give a minimum-loss design equation in closed form. For more information, see §15.4.2, Fundamentals of Power Electronics, Erickson and Maksimovic (2001).)
Geometry matters, because the core can be long and narrow, fitting many short turns, or wide and short, fitting few but longer turns. The latter are generally preferred for pulse transformers and RF baluns, where winding length limits upper bandwidth cutoff (example: P "pot" cores, binocular cores, etc.). A balanced design is preferred for power conversion, where small size and low cost are demanded, which is more or less where the cuts and stacks of conventional laminated iron come from (especially "wasteless" E-I cores, the geometry of which (and probably explanation as well) I believe I flipped past in your reference). Toroids are also popular (geometry is flexible by way of core strip width and height versus inner diameter, and assembly is relatively cheap).
So we have the winding window \$A_W\$, a fill factor \$k\$ (how much of it we can fill with copper), and some current density in that copper. We can freely assume one turn at equal current density (this of course wouldn't happen at your frequencies for a transformer of this size, were it truly solid copper, but this assumption is validated by needing many turns -- or use litz wire). Thus,
$$ I = i k A_W $$
is the total current through one turn. Note that a transformer has primary and secondary so we must halve the area and dedicate it to each winding.
We have V (for N = 1, the voltage per turn) and I, which we can multiply:
$$ VI = S = 2.22 i k A_W F B A_e $$
S is not quite power, but is the apparent power capacity (in VAs) of the transformer. When PF = 1, it is power transformed.
Cores are sometimes tabulated by area product, a sort of scaling factor. We can rearrange for that:
$$ \frac{S}{2.22 i k F B} = A_W A_e $$
Cores in a given family/series have the same aspect ratio, so it's meaningful to take the square root of this.
$$ \sqrt{ \frac{S}{2.22 i k F B} } = G A_e $$
introducing G as a geometric factor, since we don't have a family to read the Aw/Ae off of.
We might further assume k = 0.7, a typical figure for round wire and tape insulation between layers. (This can be 0.3 or less for SMPS with litz wire on bobbins, up to maybe 0.9 for carefully designed thin- or no-bobbin construction using square or flat wire.) I'll also push the coefficient to the top.
$$ \sqrt{ \frac{0.64 S}{i F B} } = G A_e $$
Further pushing around coefficients, we could move a 1/10 to the denominator, and change B from tesla to gauss, and then we have the same form as the reference. Evidently their geometric factor G is around 0.0676.
You'll have to read the context, to see what range of cores, geometry and loss density they have assumed, to arrive at such a figure. Or they might not be discussed at all, but implicit in the figures (\$i\$ and \$B_m\$), in which case you need to look up material properties to see what those figures do in turn.
Applying this for typical ferrite parameters (S = 500 VA, \$i\$ = 5 A/mm^2, F = 500Hz, B = 300mT), we need \$G A_e\$ = 653 mm^2. A typical ferrite E core (one of the few shapes available in large sizes), E65/32/27, has Ae = 540 mm^2, Aw = 394 mm^2, which seems reasonable as a start. You may need the next size up.
FYI, 500Hz is well within the range of sheet steel; if you can find a source for very fine material (under 0.1mm thickness is available), quite reasonable losses can be made, while keeping build size small (B up to 1.2T may be reasonable -- though it might also be cooking pretty good by then; look up loss curves). Grain oriented material (GOSS) is preferred.
Also, cores of this size, would normally be used as several kW at high frequency. It's entirely possible your application might benefit from a more complex solution: converting to DC, to high frequency for conversion, and back to (LF) AC. This is standard fare for inverters (regular or grid-tie), motor drivers, etc. Semiconductors are so cheap these days that production unit costs would be dominated by such a component -- assuming, of course, you need enough of them that the engineering complexity of such a design can be amortized.