I'm a bit of a student of color space theory (my father worked at tektronix during my childhood on color printers and display devices) also I worked in media streaming for a while and got familiar with the limitations and eccentricities of computer color spaces such as CCIR-601, the most popular YCrCb type color space and Rec 709 and their many many relatives:
Note that human perceivable color is mostly represented by what's referred to as CIE 1931* which is at present a superset of colors and shades representable by display devices, which is why much of this is inexact, up to interpretation or depends on uncontrollable environmental factors such as interior lighting or the weather, which is crudely modeled in the form of a white point coordinate.
I wasn't able to retrieve t_convert.html but I gleaned the following from the other two:
The specific matrix constants use appear to come from here (at least):
So it can be explained that while the CSS TR converts first to CIEXYZ at D65 white point* (Note that white point affects human color perception) then to CIEXYZ at D50 white point, then to linear rgb, the d3 code skips a step by cribbing a matrix approximation of CIEXYZ directly to linear rgb at D50. Therefore, I believe the D3 code is using (implicitly) illuminant d65 as specified for the standard ciexyz coordinate system.
My read of the lab2xyz and xyz2lab functions in the d3 code is that they precisely mirror the function described on the CIE description on wikipedia* Note the definitions of t2 and t3 equal 3δ² and δ³.
rgb2lrgb is converting sRGB (the gamma corrected, "standard" sRGB* color space to linear RGB space with no gamma correction). Their exact representation appears to match the wikipedia entry. This conversion is necessary because CIEXYZ is a linear color space and standard RGB is not.
Regarding CIELUV, I might be missing the context, but cylindrical CIELUV appears to be a simple polar conversion in the XZ CIELUV plane, and not hue corrected as in relatives of the conic HLS space. My read is that CIELUV does not bias color coordinates by supposing perceptual overlap between primaries as the CIELAB space attempts to, making it simpler to reason about (since each primary axis is separately offset in the direction of the white point's color temperature).
Hope some of this helps and it isn't tl;dr or missing the point (or stuff you already researched).