An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.

Atmospheric plasma treatment of porous polymer constructs for tissue engineering applications

Porous 3D polymer scaffolds prepared by TIPS from PLGA (53:47) and PS are intrinsically hydrophobic which prohibits the wetting of such porous media by water. This limits the application of these materials for the fabrication of scaffolds as supports for cell adhesion/spreading. Here we demonstrate that the interior surfaces of polymer scaffolds can be effectively modified using atmospheric air plasma (AP). Polymer films (2D) were also modified as control. The surface properties of wet 2D and 3D scaffolds were characterised using zeta-potential and wettability measurements. These techniques were used as the primary screening methods to assess surface chemistry and the wettability of wet polymer constructs prior and after the surface treatment. The surfaces of the original polymers are rather hydrophobic as highlighted but contain acidic functional groups. Increased exposure to AP improved the water wetting of the treated surfaces because of the formation of a variety of oxygen and nitrogen containing functions. The morphology and pore structure was assessed using SEM and a liquid displacement test. The PLGA and PS foam samples have central regions which are open porous interconnected networks with maximum pore diameters of 49 μm for PLGA and 73 μm for PS foams. (Figure Presented) © 2007 Wiley-VCH Verlag GmbH & Co. KGaA.

Supervised learning of graph structure

Graph-based representations have been used with considerable success in computer vision in the abstraction and recognition of object shape and scene structure. Despite this, the methodology available for learning structural representations from sets of training examples is relatively limited. In this paper we take a simple yet effective Bayesian approach to attributed graph learning. We present a naïve node-observation model, where we make the important assumption that the observation of each node and each edge is independent of the others, then we propose an EM-like approach to learn a mixture of these models and a Minimum Message Length criterion for components selection. Moreover, in order to avoid the bias that could arise with a single estimation of the node correspondences, we decide to estimate the sampling probability over all the possible matches. Finally we show the utility of the proposed approach on popular computer vision tasks such as 2D and 3D shape recognition. © 2011 Springer-Verlag.

Responsiveness of voltage-gated calcium channels in SH-SY5Y human neuroblastoma cells on quasi-three-dimensional micropatterns formed with poly (l-lactic acid)

Introduction: In this study, quasi-three-dimensional (3D) microwell patterns were fabricated with poly (l-lactic acid) for the development of cell-based assays, targeting voltage-gated calcium channels (VGCCs). Methods and materials: SH-SY5Y human neuroblastoma cells were interfaced with the microwell patterns and found to grow as two dimensional (2D), 3D, and near two dimensional (N2D), categorized on the basis of the cells’ location in the pattern. The capability of the microwell patterns to support 3D cell growth was evaluated in terms of the percentage of the cells in each growth category. Cell spreading was analyzed in terms of projection areas under light microscopy. SH-SY5Y cells’ VGCC responsiveness was evaluated with confocal microscopy and a calcium fluorescent indicator, Calcium GreenTM-1. The expression of L-type calcium channels was evaluated using immunofluorescence staining with DM-BODIPY. Results: It was found that cells within the microwells, either N2D or 3D, showed more rounded shapes and less projection areas than 2D cells on flat poly (l-lactic acid) substrates. Also, cells in microwells showed a significantly lower VGCC responsiveness than cells on flat substrates, in terms of both response magnitudes and percentages of responsive cells, upon depolarization with 50 mM K+. This lower VGCC responsiveness could not be explained by the difference in L-type calcium channel expression. For the two patterns addressed in this study, N2D cells consistently exhibited an intermediate value of either projection areas or VGCC responsiveness between those for 2D and 3D cells, suggesting a correlative relation between cell morphology and VGCC responsiveness. Conclusion: These results suggest that the pattern structure and therefore the cell growth characteristics were critical factors in determining cell VGCC responsiveness and thus provide an approach for engineering cell functionality in cell-based assay systems and tissue engineering scaffolds.

Gene delivery for the treatment of prostate cancer

Prostate cancer is one of the most common cancers diagnosed in men. Whilst treatments for early-stage disease are largely effective, current therapies for metastatic prostate cancer, particularly for bone metastasis, offer only a few months increased lifespan at best. Hence new treatments are urgently required. Small interfering RNA (siRNA) has been investigated for the treatment of prostate cancer where it can ‘silence’ specific cancer-related genes. However the clinical application of siRNA-based gene therapy is limited due to the absence of an optimised gene delivery vector. The optimisation of such gene delivery vectors is routinely undertaken in vitro using 2D cell culture on plastic dishes which does not accurately simulate the in vivo bone cancer metastasis microenvironment. The goal of this thesis was to assess the potential of two different targeted delivery vectors (gold or modified β-cyclodextrin derivatives) to facilitate siRNA receptor-mediated uptake into prostate cancer cells. Furthermore, this project aimed to develop a more physiologically relevant 3D in vitro cell culture model, to mimic prostate cancer bone metastasis, which is suitable for evaluating the delivery of nanoparticulate gene therapeutics. In the first instance, cationic derivatives of gold and β-cyclodextrin were synthesized to complex anionic siRNA. The delivery vectors were targeted to prostate cancer cells using the anisamide ligand which has high affinity for the sigma receptor that is overexpressed by prostate cancer cells. The gold nanoparticle demonstrated high levels of uptake into prostate cancer PC3 cells and efficient gene silencing when transfection was performed in serum-free media. However, due to the absence of a poly(ethylene glycol) (PEG) stabilising group, the formulation was unsuitable for use in serum-containing conditions. Conversely, the modified β-cyclodextrin formulation demonstrated enhanced stability in the presence of serum due to the inclusion of a PEG chain onto which the anisamide ligand was conjugated. However, the maximum level of gene silencing efficacy from three different prostate cancer cell lines (DU145, VCaP and PC3 cells) was 30 %, suggesting that further optimisation of the formulation would be required prior to application in vivo. In order to develop a more physiologically-relevant in vitro model of prostate cancer bone metastasis, prostate cancer cells (PC3 and LNCaP cells) were cultured in 3D on collagenbased scaffolds engineered to mimic the bone microenvironment. While the model was suitable for assessing nanoparticle-mediated gene knockdown, prostate cancer cells demonstrated a phenotype with lower invasive potential when grown on the scaffolds relative to standard 2D cell culture. Hence, prostate cancer cells (PC3 and LNCaP cells) were subsequently co-cultured with bone osteoblast cells (hFOB 1.19 cells) to enhance the physiological relevance of the model. Co-cultures secreted elevated levels of the MMP9 enzyme, a marker of prostate cancer metastasis, relative to prostate cancer cell monocultures (2D and 3D) indicating enhanced physiological relevance of the model. Furthermore, the coculture model proved suitable for investigating nanoparticle-mediated gene silencing. In conclusion, the work outlined in this thesis identified two different sigma receptor-targeted gene delivery vectors with potential for the treatment of prostate cancer. In addition, a more physiologically relevant model of prostate cancer bone metastasis was developed with the capacity to help optimise gene delivery vectors for the treatment of prostate cancer.

On the context-dependent scaling of consumer feeding rates

The stability of consumer-resource systems can depend on the form of feeding interactions (i.e. functional responses). Size-based models predict interactions - and thus stability - based on consumer-resource size ratios. However, little is known about how interaction contexts (e.g. simple or complex habitats) might alter scaling relationships. Addressing this, we experimentally measured interactions between a large size range of aquatic predators (4-6400 mg over 1347 feeding trials) and an invasive prey that transitions among habitats: from the water column (3D interactions) to simple and complex benthic substrates (2D interactions). Simple and complex substrates mediated successive reductions in capture rates - particularly around the unimodal optimum - and promoted prey population stability in model simulations. Many real consumer-resource systems transition between 2D and 3D interactions, and along complexity gradients. Thus, Context-Dependent Scaling (CDS) of feeding interactions could represent an unrecognised aspect of food webs, and quantifying the extent of CDS might enhance predictive ecology.

The Status of Multi-Dimensional Core-Collapse Supernova Models

Models of neutrino-driven core-collapse supernova explosions have matured considerably in recent years. Explosions of low-mass progenitors can routinely be simulated in 1D, 2D, and 3D. Nucleosynthesis calculations indicate that these supernovae could be contributors of some lighter neutron-rich elements beyond iron. The explosion mechanism of more massive stars remains under investigation, although first 3D models of neutrino-driven explosions employing multi-group neutrino transport have become available. Together with earlier 2D models and more simplified 3D simulations, these have elucidated the interplay between neutrino heating and hydrodynamic instabilities in the post-shock region that is essential for shock revival. However, some physical ingredients may still need to be added/improved before simulations can robustly explain supernova explosions over a wide range of progenitors. Solutions recently suggested in the literature include uncertainties in the neutrino rates, rotation, and seed perturbations from convective shell burning. We review the implications of 3D simulations of shell burning in supernova progenitors for the ‘perturbations-aided neutrino-driven mechanism,’ whose efficacy is illustrated by the first successful multi-group neutrino hydrodynamics simulation of an 18 solar mass progenitor with 3D initial conditions. We conclude with speculations about the impact of 3D effects on the structure of massive stars through convective boundary mixing.

Core Reactor Calculation Using the Adaptive Remeshing with a Current Error Estimator

With the objective to improve the reactor physics calculation on a 2D and 3D nuclear reactor via the Diffusion Equation, an adaptive automatic finite element remeshing method, based on the elementary area (2D) or volume (3D) constraints, has been developed. The adaptive remeshing technique, guided by a posteriori error estimator, makes use of two external mesh generator programs: Triangle and TetGen. The use of these free external finite element mesh generators and an adaptive remeshing technique based on the current field continuity show that they are powerful tools to improve the neutron flux distribution calculation and by consequence the power solution of the reactor core even though they have a minor influence on the critical coefficient of the calculated reactor core examples. Two numerical examples are presented: the 2D IAEA reactor core numerical benchmark and the 3D model of the Argonauta research reactor, built in Brasil.

Coupled CFD Shape Optimization for aerodynamic profiles

The present document deals with the optimization of shape of aerodynamic profiles -- The objective is to reduce the drag coefficient on a given profile without penalising the lift coefficient -- A set of control points defining the geometry are passed and parameterized as a B-Spline curve -- These points are modified automatically by means of CFD analysis -- A given shape is defined by an user and a valid volumetric CFD domain is constructed from this planar data and a set of user-defined parameters -- The construction process involves the usage of 2D and 3D meshing algorithms that were coupled into own- code -- The volume of air surrounding the airfoil and mesh quality are also parametrically defined -- Some standard NACA profiles were used by obtaining first its control points in order to test the algorithm -- Navier-Stokes equations were solved for turbulent, steady-state ow of compressible uids using the k-epsilon model and SIMPLE algorithm -- In order to obtain data for the optimization process an utility to extract drag and lift data from the CFD simulation was added -- After a simulation is run drag and lift data are passed to the optimization process -- A gradient-based method using the steepest descent was implemented in order to define the magnitude and direction of the displacement of each control point -- The control points and other parameters defined as the design variables are iteratively modified in order to achieve an optimum -- Preliminary results on conceptual examples show a decrease in drag and a change in geometry that obeys to aerodynamic behavior principles

Conversor de energia das ondas em energia elétrica com dispositivo de coluna de água oscilante: simulação numérica e estudo geométrico

Os oceanos representam um dos maiores recursos naturais, possuindo expressivo potencial energético, podendo suprir parte da demanda energética mundial. Nas últimas décadas, alguns dispositivos destinados à conversão da energia das ondas dos oceanos em energia elétrica têm sido estudados. No presente trabalho, o princípio de funcionamento do conversor do tipo Coluna de Água Oscilante, do inglês Oscillating Water Colum, (OWC) foi analisado numericamente. As ondas incidentes na câmara hidro-pneumática da OWC, causam um movimento alternado da coluna de água no interior da câmara, o qual produz um fluxo alternado de ar que passa pela chaminé. O ar passa e aciona uma turbina a qual transmite energia para um gerador elétrico. O objetivo do presente estudo foi investigar a influência de diferentes formas geométricas da câmara sobre o fluxo resultante de ar que passa pela turbina, que influencia no desempenho do dispositivo. Para isso, geometrias diferentes para o conversor foram analisadas empregando modelos computacionais 2D e 3D. Um modelo computacional desenvolvido nos softwares GAMBIT e FLUENT foi utilizado, em que o conversor OWC foi acoplado a um tanque de ondas. O método Volume of Fluid (VOF) e a teoria de 2ª ordem Stokes foram utilizados para gerar ondas regulares, permitindo uma interação mais realista entre o conversor, água, ar e OWC. O Método dos Volumes Finitos (MVF) foi utilizado para a discretização das equações governantes. Neste trabalho o Contructal Design (baseado na Teoria Constructal) foi aplicado pela primeira vez em estudos numéricos tridimensionais de OWC para fim de encontrar uma geometria que mais favorece o desempenho do dispositivo. A função objetivo foi a maximização da vazão mássica de ar que passa através da chaminé do dispositivo OWC, analisado através do método mínimos quadrados, do inglês Root Mean Square (RMS). Os resultados indicaram que a forma geométrica da câmara influencia na transformação da energia das ondas em energia elétrica. As geometrias das câmaras analisadas que apresentaram maior área da face de incidência das ondas (sendo altura constante), apresentaram também maior desempenho do conversor OWC. A melhor geometria, entre os casos desse estudo, ofereceu um ganho no desempenho do dispositivo em torno de 30% maior.

Random path sampling applied to vector models of magnetism

The recently reported Monte Carlo Random Path Sampling method (RPS) is here improved and its application is expanded to the study of the 2D and 3D Ising and discrete Heisenberg models. The methodology was implemented to allow use in both CPU-based high-performance computing infrastructures (C/MPI) and GPU-based (CUDA) parallel computation, with significant computational performance gains. Convergence is discussed, both in terms of free energy and magnetization dependence on field/temperature. From the calculated magnetization-energy joint density of states, fast calculations of field and temperature dependent thermodynamic properties are performed, including the effects of anisotropy on coercivity, and the magnetocaloric effect. The emergence of first-order magneto-volume transitions in the compressible Ising model is interpreted using the Landau theory of phase transitions. Using metallic Gadolinium as a real-world example, the possibility of using RPS as a tool for computational magnetic materials design is discussed. Experimental magnetic and structural properties of a Gadolinium single crystal are compared to RPS-based calculations using microscopic parameters obtained from Density Functional Theory.

Spectral self-action of THz emission from ionizing two-color laser pulses in gases

The spectrum of terahertz (THz) emission in gases via ionizing two-color femtosecond pulses is analyzed by means of a semi-analytic model and numerical simulations in 1D, 2D and 3D geometries taking into account propagation effects of both pump and THz fields. We show that produced THz signals interact with free electron trajectories and thus significantly influence further THz generation upon propagation, i.e., make the process inherently nonlocal. This self-action contributes to the observed strong spectral broadening of the generated THz field. Weshow that diffraction of the generated THz radiation is the limiting factor for the co-propagating low frequency amplitudes and thus for the self-action mechanism in 2D and 3D geometries.

A comparative study of navigation meshes

International audience

Contributions to gait recognition using multiple views

Esta tesis se centra en la identificación de personas a través de la forma de caminar. El problema del reconocimiento del paso ha sido tratado mediante diferentes enfoques, en los dominios 2D y 3D, y usando una o varias vistas. Sin embargo, la dependencia con respecto al punto de vista, y por tanto de la trayectoria del sujeto al caminar sigue siendo aún un problema abierto. Se propone hacer frente al problema de la dependencia con respecto a la trayectoria por medio de reconstrucciones 3D de sujetos caminando. El uso de reconstrucciones varias ventajas que cabe destacar. En primer lugar, permite explotar una mayor cantidad de información en contraste con los métodos que extraen los descriptores de la marcha a partir de imágenes, en el dominio 2D. En segundo lugar, las reconstrucciones 3D pueden ser alineadas a lo largo de la marcha como si el sujeto hubiera caminado en una cinta andadora, proporcionando así una forma de analizar el paso independientemente de la trayectoria seguida. Este trabajo propone tres enfoques para resolver el problema de la dependencia a la vista: 1. Mediante la utilización de reconstrucciones volumétricas alineadas. 2. Mediante el uso de reconstrucciones volumétricas no alineadas. 3. Sin usar reconstrucciones. Se proponen además tres tipos de descriptores. El primero se centra en describir el paso mediante análisis morfológico de los volúmenes 3D alineados. El segundo hace uso del concepto de entropa de la información para describir la dinámica del paso humano. El tercero persigue capturar la dinámica de una forma invariante a rotación, lo cual lo hace especialmente interesante para ser aplicado tanto en trayectorias curvas como rectas, incluyendo cambios de dirección. Estos enfoques han sido probados sobre dos bases de datos públicas. Ambas están especialmente diseñadas para tratar el problema de la dependencia con respecto al punto de vista, y por tanto de la dependencia con respecto a la trayectoria. Los resultados experimentales muestran que para el enfoque basado en reconstrucciones volumétricas alineadas, el descriptor del paso basado en entropa consigue los mejores resultados, en comparación con métodos estrechamente relacionados del Estado del Arte actual. No obstante, el descriptor invariante a rotación consigue una tasa de reconocimiento que supera a los métodos actuales sin requerir la etapa previa de alineamiento de las reconstrucciones 3D.

Spectral self-action of THz emission from ionizing two-color laser pulses in gases

The spectrum of terahertz (THz) emission in gases via ionizing two-color femtosecond pulses is analyzed by means of a semi-analytic model and numerical simulations in 1D, 2D and 3D geometries taking into account propagation effects of both pump and THz fields. We show that produced THz signals interact with free electron trajectories and thus significantly influence further THz generation upon propagation, i.e., make the process inherently nonlocal. This self-action contributes to the observed strong spectral broadening of the generated THz field. We show that diffraction of the generated THz radiation is the limiting factor for the co-propagating low frequency amplitudes and thus for the self-action mechanism in 2D and 3D geometries.