Automated mixed dimensional modelling from 2D and 3D CAD models

The motivation for this paper is to present procedures for automatically creating idealised finite element models from the 3D CAD solid geometry of a component. The procedures produce an accurate and efficient analysis model with little effort on the part of the user. The technique is applicable to thin walled components with local complex features and automatically creates analysis models where 3D elements representing the complex regions in the component are embedded in an efficient shell mesh representing the mid-faces of the thin sheet regions. As the resulting models contain elements of more than one dimension, they are referred to as mixed dimensional models. Although these models are computationally more expensive than some of the idealisation techniques currently employed in industry, they do allow the structural behaviour of the model to be analysed more accurately, which is essential if appropriate design decisions are to be made. Also, using these procedures, analysis models can be created automatically whereas the current idealisation techniques are mostly manual, have long preparation times, and are based on engineering judgement. In the paper the idealisation approach is first applied to 2D models that are used to approximate axisymmetric components for analysis. For these models 2D elements representing the complex regions are embedded in a 1D mesh representing the midline of the cross section of the thin sheet regions. Also discussed is the coupling, which is necessary to link the elements of different dimensionality together. Analysis results from a 3D mixed dimensional model created using the techniques in this paper are compared to those from a stiffened shell model and a 3D solid model to demonstrate the improved accuracy of the new approach. At the end of the paper a quantitative analysis of the reduction in computational cost due to shell meshing thin sheet regions demonstrates that the reduction in degrees of freedom is proportional to the square of the aspect ratio of the region, and for long slender solids, the reduction can be proportional to the aspect ratio of the region if appropriate meshing algorithms are used.

Electroconvection in a homeotropic bent-rod nematic liquid crystal beyond the dielectric inversion frequency

We characterize the structural transitions in an initially homeotropic bent-rod nematic liquid crystal excited by ac fields of frequency f well above the dielectric inversion point f(i). From the measured principal dielectric constants and electrical conductivities of the compound, the Carr-Helfrich conduction regime is anticipated to extend into the sub-megahertz region. Periodic patterned states occur through secondary bifurcations from the Freedericksz distorted state. An anchoring transition between the bend Freedericksz (1317) and degenerate planar (DP) states is detected. The BF state is metastable well above the Freedericksz threshold and gives way to the DP state, which persists in the field-off condition for several hours. Numerous +1 and -1 umbilics form at the onset of BF distortion, the former being largely of the chiral type. They survive in the DP configuration as linear defects, nonsingular in the core. In the BF regime, not far from fi, periodic Williams-like domains form around the umbilics; they drift along the director easy axis right from their onset. With increasing f, the wave vector of the periodic domains switches from parallel to normal disposition with respect to the c vector. Well above fi, a broadband instability is found.

Drifting periodic structures in a degenerate-planar bent-rod nematic liquid crystal beyond the dielectric inversion frequency

We report on the electric-field-generated effects in the nematic phase of a twin mesogen formed of bent-core and calamitic units, aligned homeotropically in the initial ground state and examined beyond the dielectric inversion point. The bend-Freedericksz (BF) state occurring at the primary bifurcation and containing a network of umbilics is metastable; we focus here on the degenerate planar (DP) configuration that establishes itself at the expense of the BF state in the course of an anchoring transition. In the DP regime, normal rolls, broad domains, and chevrons (both defect-mediated and defect-free types) form at various linear defect-sites, in different regions of the frequency-voltage plane. A significant novel aspect common to all these patterned states is the sustained propagative instability, which does not seem explicable on the basis of known driving mechanisms.

QR decomposition-based matrix inversion for high performance embedded MIMO receivers

Real-time matrix inversion is a key enabling technology in multiple-input multiple-output (MIMO) communications systems, such as 802.11n. To date, however, no matrix inversion implementation has been devised which supports real-time operation for these standards. In this paper, we overcome this barrier by presenting a novel matrix inversion algorithm which is ideally suited to high performance floating-point implementation. We show how the resulting architecture offers fundamentally higher performance than currently published matrix inversion approaches and we use it to create the first reported architecture capable of supporting real-time 802.11n operation. Specifically, we present a matrix inversion approach based on modified squared Givens rotations (MSGR). This is a new QR decomposition algorithm which overcomes critical limitations in other QR algorithms that prohibits their application to MIMO systems. In addition, we present a novel modification that further reduces the complexity of MSGR by almost 20%. This enables real-time implementation with negligible reduction in the accuracy of the inversion operation, or the BER of a MIMO receiver based on this.

Anisotropic seismic inversion using a multigrid Monte Carlo approach

We propose a new approach for the inversion of anisotropic P-wave data based on Monte Carlo methods combined with a multigrid approach. Simulated annealing facilitates objective minimization of the functional characterizing the misfit between observed and predicted traveltimes, as controlled by the Thomsen anisotropy parameters (epsilon, delta). Cycling between finer and coarser grids enhances the computational efficiency of the inversion process, thus accelerating the convergence of the solution while acting as a regularization technique of the inverse problem. Multigrid perturbation samples the probability density function without the requirements for the user to adjust tuning parameters. This increases the probability that the preferred global, rather than a poor local, minimum is attained. Undertaking multigrid refinement and Monte Carlo search in parallel produces more robust convergence than does the initially more intuitive approach of completing them sequentially. We demonstrate the usefulness of the new multigrid Monte Carlo (MGMC) scheme by applying it to (a) synthetic, noise-contaminated data reflecting an isotropic subsurface of constant slowness, horizontally layered geologic media and discrete subsurface anomalies; and (b) a crosshole seismic data set acquired by previous authors at the Reskajeage test site in Cornwall, UK. Inverted distributions of slowness (s) and the Thomson anisotropy parameters (epsilon, delta) compare favourably with those obtained previously using a popular matrix-based method. Reconstruction of the Thomsen epsilon parameter is particularly robust compared to that of slowness and the Thomsen delta parameter, even in the face of complex subsurface anomalies. The Thomsen epsilon and delta parameters have enhanced sensitivities to bulk-fabric and fracture-based anisotropies in the TI medium at Reskajeage. Because reconstruction of slowness (s) is intimately linked to that epsilon and delta in the MGMC scheme, inverted images of phase velocity reflect the integrated effects of these two modes of anisotropy. The new MGMC technique thus promises to facilitate rapid inversion of crosshole P-wave data for seismic slownesses and the Thomsen anisotropy parameters, with minimal user input in the inversion process.

Temporal characterization and in vitro comparison of cell survival following the delivery of 3D-conformal, intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT)

A phantom was designed and implemented for the delivery of treatment plans to cells in vitro. Single beam, 3D-conformal radiotherapy (3D-CRT) plans, inverse planned five-field intensity-modulated radiation therapy (IMRT), nine-field IMRT, single-arc volumetric modulated arc therapy (VMAT) and dual-arc VMAT plans were created on a CT scan of the phantom to deliver 3 Gy to the cell layer and verified using a Farmer chamber, 2D ionization chamber array and gafchromic film. Each plan was delivered to a 2D ionization chamber array to assess the temporal characteristics of the plan including delivery time and 'cell's eye view' for the central ionization chamber. The effective fraction time, defined as the percentage of the fraction time where any dose is delivered to each point examined, was also assessed across 120 ionization chambers. Each plan was delivered to human prostate cancer DU-145 cells and normal primary AGO-1522b fibroblast cells. Uniform beams were delivered to each cell line with the delivery time varying from 0.5 to 20.54 min. Effective fraction time was found to increase with a decreasing number of beams or arcs. For a uniform beam delivery, AGO-1552b cells exhibited a statistically significant trend towards increased survival with increased delivery time. This trend was not repeated when the different modulated clinical delivery methods were used. Less sensitive DU-145 cells did not exhibit a significant trend towards increased survival with increased delivery time for either the uniform or clinical deliveries. These results confirm that dose rate effects are most prevalent in more radiosensitive cells. Cell survival data generated from uniform beam deliveries over a range of dose rates and delivery times may not always be accurate in predicting response to more complex delivery techniques, such as IMRT and VMAT.

Conventional 3D staging PET/CT in CT simulation for lung cancer: impact of rigid and deformable target volume alignments for radiotherapy treatment planning

Objective: Positron emission tomography (PET)/CT scans can improve target definition in radiotherapy for non-small cell lung cancer (NSCLC). As staging PET/CT scans are increasingly available, we evaluated different methods for co-registration of staging PET/CT data to radiotherapy simulation (RTP) scans.

Methods: 10 patients underwent staging PET/CT followed by RTP PET/CT. On both scans, gross tumour volumes (GTVs) were delineated using CT (GTVCT) and PET display settings. Four PET-based contours (manual delineation, two threshold methods and a source-to-background ratio method) were delineated. The CT component of the staging scan was co-registered using both rigid and deformable techniques to the CT component of RTP PET/CT. Subsequently rigid registration and deformation warps were used to transfer PET and CT contours from the staging scan to the RTP scan. Dice’s similarity coefficient (DSC) was used to assess the registration accuracy of staging-based GTVs following both registration methods with the GTVs delineated on the RTP PET/CT scan.

Results: When the GTVCT delineated on the staging scan after both rigid registration and deformation was compared with the GTVCT on the RTP scan, a significant improvement in overlap (registration) using deformation was observed (mean DSC 0.66 for rigid registration and 0.82 for deformable registration, p50.008). A similar comparison for PET contours revealed no significant improvement in overlap with the use of deformable registration.

Conclusions: No consistent improvements in similarity measures were observed when deformable registration was used for transferring PET-based contours from a staging PET/CT. This suggests that currently the use of rigid registration remains the most appropriate method for RTP in NSCLC.

A comparison of nonstaggered compact FDTD schemes for the 3D wave equation

This paper aims at providing a better insight into the 3D approximations of the wave equation using compact finite-difference time-domain (FDTD) schemes in the context of room acoustic simulations. A general family of 3D compact explicit and implicit schemes based on a nonstaggered rectilinear grid is analyzed in terms of stability, numerical error, and accuracy. Various special cases are compared and the most accurate explicit and implicit schemes are identified. Further considerations presented in the paper include the direct relationship with other numerical approaches found in the literature on room acoustic modeling such as the 3D digital waveguide mesh and Yee's staggered grid technique.