Comparison of actual surgical outcomes and 3-dimensional surgical simulations

PURPOSE: The advent of imaging software programs has proved to be useful for diagnosis, treatment planning, and outcome measurement, but precision of 3-dimensional (3D) surgical simulation still needs to be tested. This study was conducted to determine whether the virtual surgery performed on 3D models constructed from cone-beam computed tomography (CBCT) can correctly simulate the actual surgical outcome and to validate the ability of this emerging technology to recreate the orthognathic surgery hard tissue movements in 3 translational and 3 rotational planes of space. MATERIALS AND METHODS: Construction of pre- and postsurgery 3D models from CBCTs of 14 patients who had combined maxillary advancement and mandibular setback surgery and 6 patients who had 1-piece maxillary advancement surgery was performed. The postsurgery and virtually simulated surgery 3D models were registered at the cranial base to quantify differences between simulated and actual surgery models. Hotelling t tests were used to assess the differences between simulated and actual surgical outcomes. RESULTS: For all anatomic regions of interest, there was no statistically significant difference between the simulated and the actual surgical models. The right lateral ramus was the only region that showed a statistically significant, but small difference when comparing 2- and 1-jaw surgeries. CONCLUSIONS: Virtual surgical methods were reliably reproduced. Oral surgery residents could benefit from virtual surgical training. Computer simulation has the potential to increase predictability in the operating room.

Contributions to the theory of magnetorotational instability and waves in a rotating plasma

The one-fluid magnetohydrodynamic (MHD) theory of magnetorotational instability (MRI) in an ideal plasma is presented. The theory predicts the possibility of MRI for arbitrary 0, where 0 is the ratio of the plasma pressure to the magnetic field pressure. The kinetic theory of MRI in a collisionless plasma is developed. It is demonstrated that as in the ideal MHD, MRI can occur in such a plasma for arbitrary P. The mechanism of MRI is discussed; it is shown that the instability appears because of a perturbed parallel electric field. The electrodynamic description of MRI is formulated under the assumption that the dispersion relation is expressed in terms of the permittivity tensor; general properties of this tensor are analyzed. It is shown to be separated into the nonrotational and rotational parts. With this in mind, the first step for incorporation of MRI into the general theory of plasma instabilities is taken. The rotation effects on Alfven waves are considered.

MAGNETOHYDRODYNAMIC SIMULATIONS OF RECONNECTION AND PARTICLE ACCELERATION: THREE-DIMENSIONAL EFFECTS

Magnetic fields can change their topology through a process known as magnetic reconnection. This process in not only important for understanding the origin and evolution of the large-scale magnetic field, but is seen as a possibly efficient particle accelerator producing cosmic rays mainly through the first-order Fermi process. In this work we study the properties of particle acceleration inserted in reconnection zones and show that the velocity component parallel to the magnetic field of test particles inserted in magnetohydrodynamic (MHD) domains of reconnection without including kinetic effects, such as pressure anisotropy, the Hall term, or anomalous effects, increases exponentially. Also, the acceleration of the perpendicular component is always possible in such models. We find that within contracting magnetic islands or current sheets the particles accelerate predominantly through the first-order Fermi process, as previously described, while outside the current sheets and islands the particles experience mostly drift acceleration due to magnetic field gradients. Considering two-dimensional MHD models without a guide field, we find that the parallel acceleration stops at some level. This saturation effect is, however, removed in the presence of an out-of-plane guide field or in three-dimensional models. Therefore, we stress the importance of the guide field and fully three-dimensional studies for a complete understanding of the process of particle acceleration in astrophysical reconnection environments.

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF MAGNETIZED WINDS OF SOLAR-LIKE STARS

By means of self-consistent three-dimensional magnetohydrodynamics (MHD) numerical simulations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-beta parameter (the ratio between thermal and magnetic energy densities). This is the first study to perform such analysis solving the fully ideal three-dimensional MHD equations. We adopt in our simulations a heating parameter described by gamma, which is responsible for the thermal acceleration of the wind. We analyze winds with polar magnetic field intensities ranging from 1 to 20 G. We show that the wind structure presents characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for all cases is similar, presenting a configuration of helmet streamer-type, with zones of closed field lines and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. For the maximum magnetic intensity simulated of 20 G and solar coronal base density, the wind velocity reaches values of similar to 1000 km s(-1) at r similar to 20r(0) and a maximum temperature of similar to 6 x 10(6) K at r similar to 6r(0). The increase of the field intensity generates a larger ""dead zone"" in the wind, i.e., the closed loops that inhibit matter to escape from latitudes lower than similar to 45 degrees extend farther away from the star. The Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density and maintaining B(0) = 20 G the system recover back to slower and cooler winds. For a fixed gamma, we show that the key parameter in determining the wind velocity profile is the beta-parameter at the coronal base. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same. This degeneracy, however, can be removed if we compare other physical parameters of the wind, such as the mass-loss rate. We analyze the influence of gamma in our results and we show that it is also important in determining the wind structure.

Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Evaluation of organ-specific peripheral doses after 2-dimensional, 3-dimensional and hybrid intensity modulated radiation therapy for breast cancer based on Monte Carlo and convolution/superposition algorithms: implications for secondary cancer risk assessment.

To make a comprehensive evaluation of organ-specific out-of-field doses using Monte Carlo (MC) simulations for different breast cancer irradiation techniques and to compare results with a commercial treatment planning system (TPS). Three breast radiotherapy techniques using 6MV tangential photon beams were compared: (a) 2DRT (open rectangular fields), (b) 3DCRT (conformal wedged fields), and (c) hybrid IMRT (open conformal+modulated fields). Over 35 organs were contoured in a whole-body CT scan and organ-specific dose distributions were determined with MC and the TPS. Large differences in out-of-field doses were observed between MC and TPS calculations, even for organs close to the target volume such as the heart, the lungs and the contralateral breast (up to 70% difference). MC simulations showed that a large fraction of the out-of-field dose comes from the out-of-field head scatter fluence (>40%) which is not adequately modeled by the TPS. Based on MC simulations, the 3DCRT technique using external wedges yielded significantly higher doses (up to a factor 4-5 in the pelvis) than the 2DRT and the hybrid IMRT techniques which yielded similar out-of-field doses. In sharp contrast to popular belief, the IMRT technique investigated here does not increase the out-of-field dose compared to conventional techniques and may offer the most optimal plan. The 3DCRT technique with external wedges yields the largest out-of-field doses. For accurate out-of-field dose assessment, a commercial TPS should not be used, even for organs near the target volume (contralateral breast, lungs, heart).

First results from a 3-dimensional middle atmosphere model

Initial results are presented from a middle atmosphere extension to a version of the European Centre For Medium Range Weather Forecasting tropospheric model. The extended version of the model has been developed as part of the UK Universities Global Atmospheric Modelling Project and extends from the ground to approximately 90 km. A comprehensive solar radiation scheme is included which uses monthly averaged climatological ozone values. A linearised infrared cooling scheme is employed. The basic climatology of the model is described; the parametrization of drag due to orographically forced gravity waves is shown to have a dramatic effect on the simulations of the winter hemisphere.

Inferring the structure of the solar corona and inner heliosphere during the Maunder minimum using global thermodynamic magnetohydrodynamic simulations

Observations of the Sun’s corona during the space era have led to a picture of relatively constant, but cyclically varying solar output and structure. Longer-term, more indirect measurements, such as from 10Be, coupled by other albeit less reliable contemporaneous reports, however, suggest periods of significant departure from this standard. The Maunder Minimum was one such epoch where: (1) sunspots effectively disappeared for long intervals during a 70 yr period; (2) eclipse observations suggested the distinct lack of a visible K-corona but possible appearance of the F-corona; (3) reports of aurora were notably reduced; and (4) cosmic ray intensities at Earth were inferred to be substantially higher. Using a global thermodynamic MHD model, we have constructed a range of possible coronal configurations for the Maunder Minimum period and compared their predictions with these limited observational constraints. We conclude that the most likely state of the corona during—at least—the later portion of the Maunder Minimum was not merely that of the 2008/2009 solar minimum, as has been suggested recently, but rather a state devoid of any large-scale structure, driven by a photospheric field composed of only ephemeral regions, and likely substantially reduced in strength. Moreover, we suggest that the Sun evolved from a 2008/2009-like configuration at the start of the Maunder Minimum toward an ephemeral-only configuration by the end of it, supporting a prediction that we may be on the cusp of a new grand solar minimum.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.

Effect Of Platform Connection And Abutment Material On Stress Distribution In Single Anterior Implant-supported Restorations: A Nonlinear 3-dimensional Finite Element Analysis.

Although various abutment connections and materials have recently been introduced, insufficient data exist regarding the effect of stress distribution on their mechanical performance. The purpose of this study was to investigate the effect of different abutment materials and platform connections on stress distribution in single anterior implant-supported restorations with the finite element method. Nine experimental groups were modeled from the combination of 3 platform connections (external hexagon, internal hexagon, and Morse tapered) and 3 abutment materials (titanium, zirconia, and hybrid) as follows: external hexagon-titanium, external hexagon-zirconia, external hexagon-hybrid, internal hexagon-titanium, internal hexagon-zirconia, internal hexagon-hybrid, Morse tapered-titanium, Morse tapered-zirconia, and Morse tapered-hybrid. Finite element models consisted of a 4×13-mm implant, anatomic abutment, and lithium disilicate central incisor crown cemented over the abutment. The 49 N occlusal loading was applied in 6 steps to simulate the incisal guidance. Equivalent von Mises stress (σvM) was used for both the qualitative and quantitative evaluation of the implant and abutment in all the groups and the maximum (σmax) and minimum (σmin) principal stresses for the numerical comparison of the zirconia parts. The highest abutment σvM occurred in the Morse-tapered groups and the lowest in the external hexagon-hybrid, internal hexagon-titanium, and internal hexagon-hybrid groups. The σmax and σmin values were lower in the hybrid groups than in the zirconia groups. The stress distribution concentrated in the abutment-implant interface in all the groups, regardless of the platform connection or abutment material. The platform connection influenced the stress on abutments more than the abutment material. The stress values for implants were similar among different platform connections, but greater stress concentrations were observed in internal connections.

PRECESSING JETS AND X-RAY BUBBLES FROM NGC 1275 (3C 84) IN THE PERSEUS GALAXY CLUSTER: A VIEW FROM THREE-DIMENSIONAL NUMERICAL SIMULATIONS

The Perseus galaxy cluster is known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths; both morphologies suggest that the active galactic nucleus (AGN) jet is subject to precession. In this work, we performed three-dimensional hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, whose dynamics are coupled to a Navarro-Frenk-White dark matter gravitational potential. The AGN jet inflates cavities that become buoyantly unstable and rise up out of the cluster core. We found that under certain circumstances precession can originate multiple pairs of bubbles. For the physical conditions in the Perseus cluster, multiple pairs of bubbles are obtained for a jet precession opening angle >40 degrees acting for at least three precession periods, reproducing both radio and X-ray maps well. Based on such conditions, assuming that the Bardeen-Peterson effect is dominant, we studied the evolution of the precession opening angle of this system. We were able to constrain the ratio between the accretion disk and the black hole angular momenta as 0.7-1.4. We were also able to constrain the present precession angle to 30 degrees-40 degrees, as well as the approximate age of the inflated bubbles to 100-150 Myr.

Placental Volumes Measured by 3-Dimensional Ultrasonography in Normal Pregnancies From 12 to 40 Weeks` Gestation

Objective. The purpose of this study was to construct nomograms of placental volumes according to gestational age and estimated fetal weight. Methods. From March to November 2007, placental volumes were prospectively measured by ultrasonography in 295 normal pregnancies from 12 to 40 weeks` gestation and correlated with gestational age and estimated fetal weight. Inclusion criteria were healthy women, singleton pregnancies with normal fetal morphologic characteristics on ultrasonography, and confirmed gestational age by first-trimester ultrasonography. Results. The mean placental volume ranged from 83 cm(3) at 12 weeks to 427.7 cm(3) at 40 weeks. Linear regression yielded the following formula for the expected placental volumes (ePV) according to gestational age (GA): ePV` (cm(3)) = -64.68 + 12.31 x GA (r = 0.572; P < .001). Placental volumes also varied according to estimated fetal weight (EFW), and the following mathematical equation was also obtained by linear regression: ePV = 94.19 + 0.09 x EFW (r = 0.505; P < 0.001). Conclusions. Nomograms of placental volumes according to gestational age and estimated fetal weight were constructed, generating reference values.

Predictable drill-free screw positioning with a graduated 3-dimensional radiographic-surgical guide: A preliminary report

Introduction: Mini-implants are placed in restricted sites, requiring an accurate surgical technique. However, no systematic study has quantified technique accuracy to reliably predict the surgical risks. Therefore, a graduated 3-dimensional radiographic-surgical guide (G-RSG) was proposed, and its inaccuracy and risk index (RI) were estimated. Methods: The sample consisted of 6 subjects (4 male, 2 female), who used mini-implant anchorage. Ten drill-free screws (DFS) were placed by using the G-RSG. The central point of the mesiodistal septum width (SW) was the selected implant site on the presurgical radiograph. The distances between DFS and the adjacent teeth (5-DFS and 6-DFS) were measured to evaluate screw centralization and inaccuracy degree (ID). These distances were statistically compared by independent t tests, and inaccuracy was determined by the expression ID = (5-DFS-6-DFS)/2, which represents deviation of the mini-implant`s final position regarding the central point initially selected. Then SW, ID, and screw diameter (SO) were combined to estimate the surgical risk with RI expressed by RI = SO/SW-ID. Results: The 5-DFS and 6-DFS distances were not significantly different. The ID of the G-RSG was 0.17 mm. The low ID ensured a safe RI (<1) in spite of the restricted SW. Conclusions: The G-RSG accuracy allowed fine prediction of the final DFS position in the inter-radicular septum, with a low RI, which is a helpful tool to estimate surgical risks. (Am J Orthod Dentofacial Orthop 2009; 136: 722-35)