2D Slab Research Articles

TMET-20. A METABOLIC BIOMARKER IN GLIOBLASTOMA: A PROMISING APPLICATION OF CO-POLARIZED [1-13C]PYRUVATE AND [1-13C]DEHYDROASCORBATE

Abstract BACKGROUND Noninvasive biomarkers for treatment response in glioblastoma are needed to improve patient outcomes. Our understanding of redox biology in treatment resistance is rapidly evolving, and interrogation of tumor metabolism represents an ideal opportunity to meet this need. In this study, the utility of co-hyperpolarized [1-13C]pyruvate and [1-13C]dehydroascorbate (HP PA/DHA) to simultaneously evaluate brain tumor metabolism and oxidative stress in orthotopically implanted mice is investigated. METHODS To characterize HP PA/DHA, mice were treated with 0-4mM diethyl maleate and 1D slab dynamic acquisitions were performed over the mouse brain. As an initial evaluation of the utility of HP DHA/PA in brain tumors, athymic nude mice (n=8) underwent orthotopic U87 glioblastoma implantation. After treatment with 8Gy in one fraction, mice were injected with HP DHA/PA dissolved in D2O. T2-weighted and 13C Chemical Shift Imaging sequences were obtained on a 3T MRI using a quadrature double-tuned 1H/13C volume coil. RESULTS HP DHA/PA detected differential lactate and vitamin C generation over the mouse brain. Treatment with diethyl maleate, which depletes the cofactor glutathione, diminished fractional vitamin C generation (0.1 vs 0.03, p=0.0009). Elevation in tumor lactate generation was seen in untreated (125.1±17.6µM, p=0.005) and radiated mice (149.4±32.6µM, p=0.01). Treatment with radiation generated higher vitamin C concentrations in tumors relative to contralateral brain (46.8±4.8µM, p=0.002). Radiated tumors were distinguished by higher ratio of lactate to vitamin C compared to untreated tumor controls (3.7 vs 2.5, p=0.07), and the same relationship was preserved between treated and untreated non-tumor brain. CONCLUSION The ratio of vitamin C/lactate detects the metabolic effects of radiation in brain tumors, inspiring its’ potential as a biomarker for treatment response. To further characterize these findings, examination of survival and modulation of radiosensitivity are ongoing.

A Universal Target Plate Design Scheme for Stellarators: theoretical basis and its application to heat load control

Abstract This study introduces a novel divertor target design scheme for stellarators, grounded in mathematical treatments and tailored to control toroidal heat load distributions. Initially, a differential equation characterizing toroidally uniform heat load distribution has been established in a two-dimensional (2D) slab configuration, and its analytic solution has been obtained. Subsequently, a numerical scheme has been developed to adapt the analytic solution into the 3D surface shape of stellarator target. The effectiveness of this design scheme has been validated through simulations of the Chinese First Quasi-axisymmetric Stellarator (CFQS) using a suite of codes including HINT, FLARE and EMC3-EIRENE, where a toroidally uniform heat load distribution has been achieved with an island configuration. Further, the effects of input parameters on the target shape and heat load distribution have been studied. The robustness of the designed target has been investigated by simulation results with varying magnetic island configurations, confirming that the toroidal uniformity of heat load distribution is insensitive to changes in island configurations. Moreover, the designed target has been assessed with gas puffing of neon, which shows that neon injections effectively reduce the heat loads without altering the toroidal uniformity of heat load distributions. The proposed scheme highlights the importance of theoretical and mathematical foundations of target design, offering an advantageous alternative/complement to the traditional numerical schemes.

Fission source convergence diagnosis in Monte Carlo eigenvalue calculations by skewness and kurtosis estimation methods

In this study, a skewness estimation method (SEM) and kurtosis estimation method (KEM) are introduced to determine the number of inactive cycles in Monte Carlo eigenvalue calculations. The SEM and KEM can determine the number of inactive cycles on the basis that fully converged fission source distributions may follow normal distributions without asymmetry or outliers. Two convergence criteria values and a minimum cycle length for the SEM and KEM were determined from skewness and kurtosis analyses of the AGN-201K benchmark and 1D slab problems. The SEM and KEM were then applied to two OECD/NEA slow convergence benchmark problems to evaluate the performance and reliability of the developed methods. Results confirmed that the SEM and KEM provide appropriate and effective convergence cycles when compared to other methods and fission source density fraction trends. Also, the determined criterion value of 0.5 for both ε1 and ε2 was concluded to be reasonable. The SEM and KEM can be utilized as a new approach for determining the number of inactive cycles and judging whether Monte Carlo tally values are fully converged. In the near future, the methods will be applied to various practical problems to further examine their performance and reliability, and optimization will be performed for the convergence criteria and other parameters as well as for improvement of the methodology for practical usage.

Thermal modeling of subduction zones with prescribed and evolving 2D and 3D slab geometries

The determination of the temperature in and above the slab in subduction zones, using models where the top of the slab is precisely known, is important to test hypotheses regarding the causes of arc volcanism and intermediate-depth seismicity. While 2D and 3D models can predict the thermal structure with high precision for fixed slab geometries, a number of regions are characterized by relatively large geometrical changes over time. Examples include the flat slab segments in South America that evolved from more steeply dipping geometries to the present day flat slab geometry. We devise, implement, and test a numerical approach to model the thermal evolution of a subduction zone with prescribed changes in slab geometry over time. Our numerical model approximates the subduction zone geometry by employing time dependent deformation of a Bézier spline that is used as the slab interface in a finite element discretization of the Stokes and heat equations. We implement the numerical model using the FEniCS open source finite element suite and describe the means by which we compute approximations of the subduction zone velocity, temperature, and pressure fields. We compute and compare the 3D time evolving numerical model with its 2D analogy at cross-sections for slabs that evolve to the present-day structure of a flat segment of the subducting Nazca plate.

Convergence of alkaline earth and transition metal oxide nanotube properties toward the 2D slab limit, computational investigation for energy conversion implications

Nanoscale devices developed from low-dimensional materials (one and two-dimensional systems) have recently attracted great attention due to their critical applications in micro/nano electromechanics. However, the practical fabrication of 2D systems is often more challenging than that of their 1D counterparts. In the current investigation, the convergence of 1D nanotube properties toward the 2D slab limit has been verified and confirmed. In this regard, a variety of geometrical, electronic, and response (mechanical and piezoelectric) properties are computed for both zigzag (n,0) metal oxide MO (M = Be, Mg, Ca, Zn, Cd) nanotubes and MO 2D surfaces. The variation of such properties as a function of the tube index n (ranging from 6 to 48, corresponding to 24 to 192 atoms per cell, respectively) is highlighted. Additionally, the electronic and nuclear contributions to all aforementioned response properties have been discussed and analyzed. Interestingly, CdO (48,0) nanotube introduces considerable direct and converse piezoelectric constants: e11 = 6.04 |e|×bohr and d11 = 22.88 pm/V. These electromechanical coefficients converge well to the values of the 2D surface, while the response is entirely dominated by the ionic contribution, confirming the flexibility and softness. Such findings make these MO nanotubes good candidates for nanoscale energy conversion piezoelectric applications.

Water Self‐Dissociation is Insensitive to Nanoscale Environments

AbstractNanoconfinement effects on water dissociation and reactivity remain controversial, despite their importance to understand the aqueous chemistry at interfaces, pores, or aerosols. The pKw in confined environments has been assessed from experiments and simulations in a few specific cases, leading to dissimilar conclusions. Here, with the use of carefully designed ab initio simulations, we demonstrate that the energetics of bulk water dissociation is conserved intact to unexpectedly small length‐scales, down to aggregates of only a dozen molecules or pores of widths below 2 nm. The reason is that most of the free‐energy involved in water autoionization comes from breaking the O−H covalent bond, which has a comparable barrier in the bulk liquid, in a small droplet of nanometer size, or in a nanopore in the absence of strong interfacial interactions. Thus, dissociation free‐energy profiles in nanoscopic aggregates or in 2D slabs of 1 nm width reproduce the behavior corresponding to the bulk liquid, regardless of whether the corresponding nanophase is delimited by a solid or a gas interface. The present work provides a definite and fundamental description of the mechanism and thermodynamics of water dissociation at different scales with broader implications on reactivity and self‐ionization at the air‐liquid interface.

Water Self-Dissociation is Insensitive to Nanoscale Environments.

Nanoconfinement effects on water dissociation and reactivity remain controversial, despite their importance to understand the aqueous chemistry at interfaces, pores, or aerosols. The pKw in confined environments has been assessed from experiments and simulations in a few specific cases, leading to dissimilar conclusions. Here, with the use of carefully designed ab initio simulations, we demonstrate that the energetics of bulk water dissociation is conserved intact to unexpectedly small length-scales, down to aggregates of only a dozen molecules or pores of widths below 2 nm. The reason is that most of the free-energy involved in water autoionization comes from breaking the O-H covalent bond, which has a comparable barrier in the bulk liquid, in a small droplet of nanometer size, or in a nanopore in the absence of strong interfacial interactions. Thus, dissociation free-energy profiles in nanoscopic aggregates or in 2D slabs of 1 nm width reproduce the behavior corresponding to the bulk liquid, regardless of whether the corresponding nanophase is delimited by a solid or a gas interface. The present work provides a definite and fundamental description of the mechanism and thermodynamics of water dissociation at different scales with broader implications on reactivity and self-ionization at the air-liquid interface.

A Monte Carlo Thermal Radiative Transfer Solver with Nonlinear Elimination

In this paper, we present a new Monte Carlo method for solving the thermal radiative transfer (TRT) equations via the method of nonlinear elimination (NLEM). This method is inspired by the previous application of NLEM to thermal radiation diffusion. Our approach, called diffusion accelerated Implicit Monte Carlo (DAIMC), is a hybrid technique which combines a Monte Carlo method for solving a purely-absorbing transport equation and a diffusion solution that accounts for effective scattering, or absorption–reemission. The method aims to improve the implicitness of the traditional implicit Monte Carlo (IMC) method. We derive DAIMC generally for 3D Cartesian geometries, but in this paper, we present results and analysis in 1D slab geometry. These preliminary results indicate that DAIMC implementations may provide more accurate and robust TRT solutions than IMC in certain test problems.

A numerical study on the role of instabilities on multi-wavelength emission signatures of blazar jets

Context. Blazars, a class of active galaxies whose jets are relativistic and collimated flows of plasma directed along the line of sight, are prone to a slew of magnetohydrodynamic (MHD) instabilities. These jets show characteristic multi-wavelength and multi-timescale variabilities. Aims. We aim to study the interplay of radiation and particle acceleration processes in regulating the multi-band emission and variability signatures from blazars. In particular, the goal is to decipher the impact of shocks arising due to MHD instabilities in driving the long-term variable emission signatures from blazars. Methods. To this end, we performed relativistic MHD (RMHD) simulations of a representative section of a blazar jet. The jet was evolved using a hybrid Eulerian-Lagrangian framework to account for radiative losses due to synchrotron process as well as particle acceleration due to shocks. Additionally, we incorporated and validated radiative losses taking into consideration the external Compton (EC) process that is relevant for blazars. We further compared the effects of different radiation mechanisms through numerical simulation of 2D slab jet as a validation test. Finally, we carried out a parametric study to quantify the effect of magnetic fields and external radiation field characteristics by performing 3D simulations of a plasma column. The synthetic light curves and spectral energy distribution (SEDs) were analyzed to qualitatively understand the impact of instability driven shocks. Results. We observed that shocks produced with the evolution of instabilities give rise to flaring signatures in the high-energy band. The impact of such shocks is also evident from the instantaneous flattening of the synchrotron component of the SEDs. At later stages, we observed the transition in X-ray emission from the synchrotron process to that dominated by EC. The inclusion of the EC process also gives rise to γ-ray emission and shows signatures of mild Compton dominance that is typically seen in low-synchrotron peaked blazars.

Droplet entrainment and its role in the combustion of HTPB/paraffin fuels

Droplet entrainment offers high regression rate of HTPB/paraffin fuels, but also leads to more complicated mass and heat transfer mechanisms during combustion. The droplet entrainment of the fuels in a windowed two-dimensional (2D) slab burner were investigated by using a high-speed camera, and the convective heat transfer coefficient was modified. The results indicated that the entrained droplets could be easily generated below 1.29 MPa of combustion pressure (pc), which was dominant in the regression of fuels, whereas the regression rate decreased by approximately 0.9 mm when pc increased from 0.32 MPa to 1.29 MPa due to the negative correlation of droplet entrainment rate with pc. However, the regression at higher pc mainly depended on weak evaporation and thermal decomposition rather than droplet entrainment, which was insensitive to pc and decreased by only about 0.05 mm s−1 with the increase in pc from 1.85 MPa to 2.52 MPa. The carbonized layer could be easily formed on the burning surface at higher pc, which hindered the droplet entrainment. Additionally, the improvement of convective heat transfer coefficient with the presence of droplet entrainment was evaluated, and the blocking effect was weakened, especially under high oxidizer mass flux.

Lp-CMFD acceleration schemes in multi-energy group 2D Monte Carlo transport

The linear prolongation flux update scheme is extended to both regular CMFD acceleration, as well as partial CMFD acceleration in 2D multi energy group Monte Carlo k-eigenvalue neutron transport problems. The acceleration performance of these CMFD variants were investigated in simple 2D slab geometries, first with a monoenergetic case and then with a three group problem on the same geometry based on the monoenergetic cross sections. Flux convergence was determined via an on-the-fly convergence diagnostics developed by Ueki and Brown. It is found that on top of providing better acceleration in general, the linear prolongation scheme is also able to correct for instabilities in the CMFD scheme. Overall, the lp-pCMFD scheme employing a maximum history length is found to have the best performance across the cases presented.

Immiscible Viscous Fingering: The Simulation of Tertiary Polymer Displacements of Viscous Oils in 2D Slab Floods.

Immiscible viscous fingering in porous media occurs when a high viscosity fluid is displaced by an immiscible low viscosity fluid. This paper extends a recent development in the modelling of immiscible viscous fingering to directly simulate experimental floods where the viscosity of the aqueous displacing fluid was increased (by the addition of aqueous polymer) after a period of low viscosity water injection. This is referred to as tertiary polymer flooding, and the objective of this process is to increase the displacement of oil from the system. Experimental results from the literature showed the very surprising observation that the tertiary injection of a modest polymer viscosity could give astonishingly high incremental oil recoveries (IR) of ≥100% even for viscous oils of 7000 mPa.s. This work seeks to both explain and predict these results using recent modelling developments. For the 4 cases (µo/µw of 474 to 7000) simulated in this paper, finger patterns are in line with those observed using X-ray imaging of the sandstone slab floods. In particular, the formation of an oil bank on tertiary polymer injection is very well reproduced and the incremental oil response and water cut drops induced by the polymer are very well predicted. The simulations strongly support our earlier claim that this increase in incremental oil displacement cannot be explained solely by a viscous “extended Buckley-Leverett” (BL) linear displacement effect; referred to in the literature simply as “mobility control”. This large response is the combination of this effect (BL) along with a viscous crossflow (VX) mechanism, with the latter VX effect being the major contributor to the recovery mechanism.

Boundary layer combustion of HTPB/paraffin fuels for hybrid propulsion applications

HTPB/paraffin fuels exhibit high regression rates and excellent mechanical properties due to the droplet entrainment and enhancement of the HTPB matrix, however, the special structure of the fuels also leads to more complex combustion mechanisms that may be quite different from that of paraffin and HTPB alone. The boundary layer combustion characteristics of HTPB/paraffin fuels with oxygen flow in a windowed 2D slab burner were investigated by using a high speed camera, Schlieren and thermocouples. The results demonstrated that the regression rate of HTPB/paraffin fuels mainly composed of paraffin was higher than that of HTPB due to the dominant role of droplet entrainment in mass transfer. The flame height and combustion boundary layer thickness of the fuel were inversely correlated with the oxidizer mass flux, and decreased by about 2.8 times and 2.5 times when the oxidizer mass flux increased from 12.40 kg⋅m⋅−2s−1 to 64.86 kg⋅m⋅−2s−1, respectively, however, the chamber pressure exhibited little effect. The addition of HTPB limited the development of combustion boundary layer thickness and hindered the generation of entrainment droplets due to an increase in liquid layer viscosity, but both the proportion of flame area in the combustion boundary layer and the gas temperature gradient in the fuel-rich zone were positively correlated with HTPB content attributed to the higher combustion efficiency.

Magnetic resonance spectroscopy shows associations between neurometabolite levels and perivascular space volume in Parkinson's disease: a pilot and feasibility study.

Higher volume fraction of perivascular space (PVS) has recently been reported in Parkinson's disease (PD) and related disorders. Both elevated PVS and altered levels of neurometabolites, assayed by proton magnetic resonance spectroscopy (MRS), are suspected indicators of neuroinflammation, but no published reports have concurrently examined PVS and MRS neurometabolites. In an exploratory pilot study, we acquired multivoxel 3-T MRS using a semi-Localization by Adiabatic SElective Refocusing (sLASER) pulse-sequence (repetition time/echo time = 2810/60 ms, voxels 10 × 10 × 10 mm3) from a 2D slab sampling bilateral frontal white matter (FWM) and anterior middle cingulate cortex (aMCC). PVS maps obtained from high-resolution (0.8 × 0.8 × 0.8 mm3) T1-weighted MRI were co-registered with MRS. In each MRS voxel, PVS volume and neurometabolite levels were measured. Linear regression accounting for age, sex, and BMI found greater PVS volume for higher levels of choline-containing compounds (Cho; P = 0.047) in FWM and lower PVS volume for higher levels of N-acetyl compounds (NAA; P = 0.012) in aMCC. Since (putatively) higher Cho is associated with inflammation while NAA has anti-inflammatory properties, these observations add to evidence that higher PVS load is a sign of inflammation. Additionally, lower Montreal Cognitive Assessment scores were associated with lower NAA in aMCC (P = 0.002), suggesting that local neuronal dysfunction and inflammation contribute to cognitive impairment in PD. These exploratory findings indicate that co-analysis of PVS and MRS is feasible and may help elucidate the cellular and metabolic substrates of glymphatic and inflammatory processes in PD.

Metasurface-Dressed Two-Dimensional on-Chip Waveguide for Free-Space Light Field Manipulation.

We show that a metasurface-coated two-dimensional (2D) slab waveguide enables the generation of arbitrary complex light fields by combining the extreme versatility and freedom on the wavefront control of optical metasurfaces with the compactness of photonic integrated circuits. We demonstrated off-chip 2D focusing and holographic projection with our metasurface-dressed photonic integrated devices. This technology holds the potential for many other optical applications requiring 2D light field manipulation with full on-chip integration, such as solid-state LiDAR and near-eye AR/VR displays.

Boundary layer combustion of paraffin fuels for hybrid propulsion applications

Droplet entrainment offers paraffin fuels a regression rate of 3–4 times higher than classical polymer fuels, which is favorable for the hybrid rocket motor, however, the heterogeneous combustion reactions between droplets and oxidizer lead to a more complicated combustion mechanism. Combustion diagnosis was carried out for paraffin fuels combusted in a windowed 2D slab burner, and the boundary layer combustion was investigated by using high speed camera and schlieren in this study. The results showed that the presence of remarkable periodic flame and chamber pressure oscillation was due to the periodic accumulation and breakdown of liquid film layer; Kelvin−Helmholtz instability was a necessary but not the sufficient condition for droplet entrainment, and the critical oxidizer mass flux of Kelvin−Helmholtz instability and droplet entrainment was about 12.31 kg m−2 s−1 and 30.15 kg m−2 s−1; the thickness of turbulent boundary layer and the height of combustion flame were inversely correlated with the oxidizer mass flux, indicating the more efficient combustion under higher oxidizer mass flux; the entrainment droplets disappeared at elevated chamber pressures due to the inhibition of Kelvin−Helmholtz instability, meanwhile, the thickness of turbulent boundary layer was found to increase; entrainment droplets increased and the height of combustion flame decreased due to the acceleration of combustion reactions when the copper chromite was used as the combustion catalyst, which favored the enhancement of the regression rate of fuels and the combustion efficiency of hybrid rocket motors.

Automatic quantification of perivascular spaces in T2-weighted images at 7 T MRI.

Perivascular spaces (PVS) are believed to be involved in brain waste disposal. PVS are associated with cerebral small vessel disease. At higher field strengths more PVS can be observed, challenging manual assessment. We developed a method to automatically detect and quantify PVS. A machine learning approach identified PVS in an automatically positioned ROI in the centrum semiovale (CSO), based on -resolution T2-weighted TSE scans. Next, 3D PVS tracking was performed in 50 subjects (mean age 62.9 years (range 27-78), 19 male), and quantitative measures were extracted. Maps of PVS density, length, and tortuosity were created. Manual PVS annotations were available to train and validate the automatic method. Good correlation was found between the automatic and manual PVS count: ICC (absolute/consistency) is 0.64/0.75, and Dice similarity coefficient (DSC) is 0.61. The automatic method counts fewer PVS than the manual count, because it ignores the smallest PVS (length <2mm). For 20 subjects manual PVS annotations of a second observer were available. Compared with the correlation between the automatic and manual PVS, higher inter-observer ICC was observed (0.85/0.88), but DSC was lower (0.49 in 4 persons). Longer PVS are observed posterior in the CSO compared with anterior in the CSO. Higher PVS tortuosity are observed in the center of the CSO compared with the periphery of the CSO. Our fully automatic method can detect PVS in a 2D slab in the CSO, and extract quantitative PVS parameters by performing 3D tracking. This method enables automated quantitative analysis of PVS.

Computational discovery of two-dimensional copper chalcogenidesCuX(X= S, Se, Te)

We investigate the energy landscape of 2D CuS, CuSe, and CuTe compounds using a global evolutionary algorithm in combination with density functional theory. Four low-lying energy ${\mathrm{Cu}}_{2}{X}_{2}$ (X=S, Se, Te) slabs with $P\overline{3}m1$, $P4/mmm$, $P4/nmm$, and $Pmmn$ symmetries have been identified on their respective potential energy surfaces. Three structures present tetrahedral ${\mathrm{Cu}X}_{4}$ motifs which pave the 2D slabs, while the latter is based on shared ${\mathrm{Cu}}_{4}{X}_{2}$ octahedra. The viability of each phase was examined by looking at dynamical, thermodynamic, and thermal stability criteria. We find that 2D copper monosulfide crystallizes in the $P\overline{3}m1$ phase, reminiscent of the covalent slab found in the $\mathrm{Ca}{\mathrm{Al}}_{2}{\mathrm{Si}}_{2}$ prototype. Two polymorphs of 2D copper selenide with symmetries $P\overline{3}m1$ and $P4/mmm$ are close in energy. The latter ${\mathrm{Cu}}_{2}{\mathrm{Se}}_{2}$ slab is made of fused ${\mathrm{Cu}}_{4}{\mathrm{Se}}_{2}$ square bipyramids with electronegative Se atoms in apical positions. Finally, the ground-state 2D CuTe has a $Pmmn$ space group containing distorted ${\mathrm{CuTe}}_{4}$ tetrahedra with some Te--Te bonding. The electronic and bond analyses show that all of these predicted 2D CuX phases are metallic with ionocovalent bonds.

Evolution of anisotropic turbulence in the fast and slow solar wind: Theory and Solar Orbiter measurements

Aims. Solar Orbiter (SolO) was launched on February 9, 2020, allowing us to study the nature of turbulence in the inner heliopshere. We investigate the evolution of anisotropic turbulence in the fast and slow solar wind in the inner heliosphere using the nearly incompressible magnetohydrodynamic (NI MHD) turbulence model and SolO measurements. Methods. We calculated the two dimensional (2D) and the slab variances of the energy in forward and backward propagating modes, the fluctuating magnetic energy, the fluctuating kinetic energy, the normalized residual energy, and the normalized cross-helicity as a function of the angle between the mean solar wind speed and the mean magnetic field (θUB), and as a function of the heliocentric distance using SolO measurements. We compared the observed results and the theoretical results of the NI MHD turbulence model as a function of the heliocentric distance. Results. The results show that the ratio of 2D energy and slab energy of forward and backward propagating modes, magnetic field fluctuations, and kinetic energy fluctuations increases as the angle between the mean solar wind flow and the mean magnetic field increases from θUB = 0° to approximately θUB = 90° and then decreases as θUB → 180°. We find that solar wind turbulence is a superposition of the dominant 2D component and a minority slab component as a function of the heliocentric distance. We find excellent agreement between the theoretical results and observed results as a function of the heliocentric distance.

Ion temperature anisotropy in the tokamak scrape-off layer

Understanding tokamak exhaust-power heat loads on divertor plates depends critically on having a realistic model of the scrape-off layer (SOL) plasma. The Braginskii fluid model is often solved to understand the SOL plasma behavior. This model is based on the collisional limit for transport along the magnetic field . The ions and electron gyrofrequencies are assumed to be much larger than the Coulomb collision frequencies, which are nonetheless, sufficiently large to yield common parallel and perpendicular temperatures for each species, i.e. the temperatures are assumed to be isotropic. In certain circumstances such as encountered for the tokamak H-mode, the ion temperature can be quite anisotropic. In this work, the anisotropy effects are implemented in the two-dimensional (2D) transport code UEDGE. Various geometries (1D slab, 2D slab and a toroidal tokamak geometry) are used to study the 2D structure of ion temperature anisotropy and its effects on plasma transport in detail. Results show that the effects of ion temperature anisotropy on the plasma parallel transport are substantial near the magnetic X-point, which leads to different steady state density profiles in the divertor regions. The extra mirror force introduced by ion temperature anisotropy can be one of the main forces contributing to the plasma flow in the SOL.