Water diffusion in tissues is generally restricted and often anisotropic. Neural tissue is of particular interest, since it is well known that injury alters diffusion in a characteristic manner. Both Monte Carlo simulations and approximate analytical models have previously been reported in attempts to predict water diffusion behavior in the central nervous system. These methods have relied on axonal models, which assume simple geometries (e.g., ellipsoids, cylinders, and square prisms) and ignore the thickness of the myelin sheath. The current work describes a method for generating models using synthetic images. The computations are based on a 3D finite difference (FD) approximation of the diffusion equation. The method was validated with known analytic solutions for diffusion in a cylindrical pore and in a hexagonal array of cylinders. Therefore, it is envisioned that, by exploiting histologic images of neuronal tissues as input model, current method allows investigating the water diffusion behavior inside biological tissues and potentially as- Diffusion-sensitized MRI provides a means of probing the architectural features of structures that are much smaller than a voxel. In the CNS the degree of diffusion anisotropy is sensitive to local differences in nerve fiber orientation (1-3). Various mechanisms have been proposed to explain changes in the apparent diffusion coefficient (ADC) during pathologic states. For example, in the study of cat brain, cellular swelling in the initial stages of ischemia may account for the decreased ADC since water protons originally in the faster-diffusing extracellular space migrate into the intracellular space (4) In addition, the cellular swelling reduces the extracellular space and, thus, increases the tortuosity of the diffusion paths (5). In the rat spinal cord diffusion was found to be anisotropic (6 -8). Spinal cord injury resulted in reduced diffusion anisotropy (9) which has also been attributed to axonal swelling or loss and inflammation (10,11). Changes in the myelin sheath may also affect the ADC. For example, the demyelination resulting from mechanical trauma (12) and the varying levels of remyelination during recovery (13) may cause changes in the permeability of the barrier between the intracellular and extracellular compartments. Furthermore, regenerating and sprouting axon fibers are typically tortuous, small, and unmyelinated (14). The effect of these tissue parameters on the ADC may potentially be exploited to assess the extent of neural damage and repair in vivo. Reduced diffusion anisotropy in multiple sclerosis lesions have also been attributed to loss of axons and myelin (15,16). Diffusion-weighted imaging may thus potentially provide insight into the pathophysiologic changes in MS lesions if the relationship between morphology and the diffusion-sensitized MR signal are better understood.Previous attempts to predict the effects of morphologic changes have relied on Monte Carlo (MC) simulations (11,17) or approximate analytic expressions (18). Exact solutions...