In one such instance I've been studying for years, the WKB approximation, I've realized two things: 1) the approximations are not well understood and 2) the mathematics are quite complicated, but these points notwithstanding the equations are used in experiments. You can read the few-page introduction in "Lectures on the Geometry of Quantization" by Bates & Weinstein <https://math.berkeley.edu/~alanw/GofQ.pdf> to see some of this, in particular the subsection "Quantization and the classical limit". I'll just quote the relevant paragraph:
> Although there remain some unsettled issues connected with the question, “How can ℏ become small?” the answer is essentially the following. For any particular mechanical system, there are usually characteristic distances, masses, velocities, . . . from which a unit of action appropriate to the system can be derived, and the classical limit is applicable when ℏ divided by this unit is much less than 1.
But remember, this is just one approach to the subject. Another heuristic is this: h has dimensions energy x time, which means it converts frequency into energy, e.g. E = hf. In the Fourier transform, the character is exp(2πihx·ξ), where ξ is the frequency. The effect of h -> 0 would be to dampen high-energy waves. Irregularity comes from high frequencies (think of it like this: a sum of sines of large periods would not have many kinks.) When you "iron out" the irregularity of the quantum solution, you end up with a classical one.
Again disclaimer: not a physicist, nor an expert.
(Many books claim that Planck's constant is a quantum of "action". This is extremely wrong, because action is not a quantized physical quantity, so it cannot have quanta, and moreover the ratio between energy and frequency is not an action, but an angular momentum. This ridiculous mistake is caused by the failure to understand that the unit of plane angle is a base unit that cannot be derived from any of the units of the other physical quantities and forgetting to write the plane angle in the dimensional relationships between physical quantities leads to very serious errors. Angular momentum is the ratio between action and plane angle, while energy is the ratio between action and time and the linear momentum is the ratio between action and length. Frequency is not the inverse of time, but the ratio between plane angle and time, as the old name of "Hertz", i.e. "cycles per second" was making this obvious.)
Dividing or multiplying by combinations of Planck's constant with the speed of light and the elementary charge only changes the system of units, between the traditional units and natural units. There are several variants of "natural" systems of units and the difference between using them and using SI is that when using natural systems of units there are much less "universal" constants in the relationships between physical quantities.
Like electric charge, angular momentum is one of the quantities that are discrete, not continuous (though angular momentum is not necessarily discrete like electric charge; depending on the physical system it may be either discrete or continuous).
Like the elementary charge is the quantum of electric charge (when all quarks are bound in hadrons), a half of Planck's constant is the quantum of angular momentum (there are multiple possible definitions of the angular momentum, depending on the unit chosen for plane angle, which lead to multiple possible numeric values for Planck's constant).
The importance of the c, e and h constants is due to them being respectively the limit value for velocity and the quanta of two fundamental discrete quantities (electric charge and angular momentum), which makes them appear in many relationships between physical quantities, unless the system of units is changed to a "natural" system of units, when the "universal" constants become "1", so they disappear from the formulae.
When a natural system of units is not used in quantum mechanics, then all the equations that contain both energy or momentum and time or length, like the equation of Schroedinger, will contain Planck's constant, possibly combined with other universal constants. This will have as a consequence the appearance of those constants also in the expressions that give the pairs of quantities related by Fourier transforms.
The fact that hbar is the quantum of angular momentum can be derived from the boundary conditions of a wavefunction with a closed dimension (such as an angle), combined with the fact that one h of action phases a wavefunction by 360º.
In particular, for a wavefunction with only angular dependence:
boson: \psi(2pi) = \psi(0) fermion: \psi(4pi) = \psi(0)
Because the value has to be the same when you get back to where you started.
Therefore, (for a boson) psi must decompose into eigenfunctions of the form \psi(\theta) = exp(i * n * \theta) where n is an integer.
These eigenfunctions satisfy the boundary condition stated above and are a complete basis
We can take the fourier transform of this (closed) function and we get a dirac delta at n. I.e. the fourier transforms of this basis are just dirac deltas at integers in the fourier domain.
But to convert this to physical units, we need to use the embedding from integers to reals, which is done by multiplication by planck's constant.
So, some related questions that I think are further down the causal hierarchy:
* Why does one h of action phase a wavefunction by 360º
* Why does h necessarily show up in the fourier transform