Skip to main content
Log in

Three-dimensional ray tracing through curvilinear interfaces with application to laser Doppler anemometry in a blood analogue fluid

  • Original Article
  • Published:
Medical & Biological Engineering & Computing Aims and scope Submit manuscript

Abstract

Prediction of the effects of refractive index (RI) mismatch on laser Doppler anemometer (LDA) measurements within a curvilinear cavity (an artificial ventricle) was achieved by developing a general technique for modelling the paths of the convergent beams of the LDA system using 3D vector geometry. Validated by ray tracing through CAD drawings, the predicted maximum tolerance in RI between the solid model and the working fluid was ± 0.0005, equivalent to focusing errors commensurate with the geometric and alignment uncertainties associated with the flow model and the LDA arrangement. This technique supports predictions of the effects of refraction within a complex geometry. Where the RI mismatch is unavoidable but known, it is possible not only to calculate the true position of the measuring volume (using the probe location and model geometry), but also to estimate degradation in signal quality arising from differential displacement and refraction of the laser beams.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8

Similar content being viewed by others

Notes

  1. These latter locations were (0, ± r z , z) when the probe was parallel to the x-axis, and (± r z , 0, z) when the probe was parallel to the y-axis.

  2. A capital Z represents an approximating function (in the zr plane), as opposed to the Cartesian coordinate.

  3. Nevertheless, the model can be readily extended should more interfaces be required.

  4. The parameters pertaining to the other three beams are readily obtained by invoking symmetry about the optic axis; this gives rise to three sets of equations analogous to Eqs. 36, featuring the subscripts P2, P3, P4.

Abbreviations

α :

Half angle of convergence of the laser beams (in air)

B :

Basepoint vector

γ :

Scalar variable

d :

y displacement of the probe objective lens from the near wall of the VAD

d min :

The minimum separation between two rays

f :

Focal length of the LDA probe (in air)

h :

x displacement of the optic axis from the origin of the VAD coordinates

k :

z displacement of the optic axis from the origin of the VAD coordinates

μ :

Fluid viscosity

n :

Refractive index

N :

Normal vector

r :

Radial coordinate (normal to the outflow axis)

ρ :

Fluid density

R :

Equatorial radius of the VAD

s :

Displacement of a beam from the optic axis at the probe objective lens

s sub :

Displacement from the optic axis of an intersection between a laser beam and a model interface (‘sub’ = F or P)

t :

Minimum thickness of the PMMA wall

T :

Transmission vector

x :

Cartesian coordinate in the plane of the diaphragm–housing junction, parallel to the projection of the inflow axis

y :

Cartesian coordinate in the plane of the diaphragm–housing junction, normal to the projection of the inflow axis

z :

Cartesian coordinate parallel to the outflow axis

Z :

Function of r describing the profile of the axially symmetric part of the VAD cavity (subscript C indicates circular arc; subscript P indicates polynomial segment)

F:

The PMMA:fluid interface

O:

outflow tract

P:

The air:PMMA interface

S:

intersection of VAD profile functions

1:

The y-most beam in the xy plane

2:

the y-least beam in the xy plane

3:

the z-most beam in the yz plane

4:

the z-least beam in the yz plane

References

  1. Baldwin JT, Tarbell JM, Deutsch S, Geselowitz DB (1989) Mean flow velocity patterns within a ventricular-assist device. ASAIO Trans 35:429–433

    Article  Google Scholar 

  2. Dybbs A, Edwards RV (1987) Refractive index matching for difficult situations. In: 2nd international conference on laser anemometry—advances and applications, paper I1. Strathclyde, UK, 21–23 September, pp 1–22

  3. Gardavsky J, Hrbek J, Chara Z, Severa M (1989) Refraction corrections for LDA measurements in circular tubes within rectangular optical boxes. Dantec Information No. 08. Dantec Dynamics A/S, Skovlunde

  4. Jin W, Clark C (1993) Experimental investigation of unsteady flow within a sac-type ventricular-assist device (VAD). J Biomech 26(6):697–707

    Article  Google Scholar 

  5. Narrow TL, Yoda M, Abdel-Khalik SI (2000) A simple model for the refractive index of sodium iodide aqueous solutions. Exp Fluids 28:282–283

    Article  Google Scholar 

  6. Nugent AH (2005) Fluid dynamical investigation of a ventricular-assist device. Ph.D. thesis, University of New South Wales

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Allen H. Nugent.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Nugent, A.H., Bertram, C.D. Three-dimensional ray tracing through curvilinear interfaces with application to laser Doppler anemometry in a blood analogue fluid. Med Biol Eng Comput 48, 147–156 (2010). https://doi.org/10.1007/s11517-009-0511-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11517-009-0511-7

Keywords

Navigation