Friendly amendment: a commentary on Doyle and Ford's proposed re-definition of 'mental model'

Some amendments are proposed to a recent redefinition of the mental model concept in system dynamics. First, externalised, or articulated mental models should not be called cognitive maps; this term has a well established, alternative meaning. Second, there can be mental models of entities not yet existing beyond an individual's mind; the modelling of planned or desired systems is possible and recommended. Third, saying that mental models maintain social systems connects with some exciting research opportunities for system dynamics; however, it is probably an accidental distraction from the intended meaning of the redefinition. These minor criticisms apart, the new definition of mental model of a dynamic system is welcomed as a useful contribution to both research and practice.

Probing the 2D kinematic structure of early-type galaxies out to three effective radii

We detail an innovative new technique for measuring the two-dimensional (2D) velocity moments (rotation velocity, velocity dispersion and Gauss-Hermite coefficients h(3) and h(4)) of the stellar populations of galaxy haloes using spectra from Keck DEIMOS (Deep Imaging Multi-Object Spectrograph) multi-object spectroscopic observations. The data are used to reconstruct 2D rotation velocity maps. Here we present data for five nearby early-type galaxies to similar to three effective radii. We provide significant insights into the global kinematic structure of these galaxies, and challenge the accepted morphological classification in several cases. We show that between one and three effective radii the velocity dispersion declines very slowly, if at all, in all five galaxies. For the two galaxies with velocity dispersion profiles available from planetary nebulae data we find very good agreement with our stellar profiles. We find a variety of rotation profiles beyond one effective radius, i.e. rotation speed remaining constant, decreasing and increasing with radius. These results are of particular importance to studies which attempt to classify galaxies by their kinematic structure within one effective radius, such as the recent definition of fast- and slow-rotator classes by the Spectrographic Areal Unit for Research on Optical Nebulae project. Our data suggest that the rotator class may change when larger galactocentric radii are probed. This has important implications for dynamical modelling of early-type galaxies. The data from this study are available on-line.

Dynamics of Bose-Einstein condensates in double-well potentials

We study the dynamics of Bose-Einstein condensates in symmetric double-well potentials following a sudden change of the potential from the Mott-insulator to the superfluid regime. We introduce a continuum approximation that maps that problem onto the wave-packet dynamics of a particle in an anharmonic effective potential. For repulsive two-body interactions the visibility of interference fringes that result from the superposition of the two condensates following a stage of ballistic expansion exhibits a collapse of coherent oscillations onto a background value whose magnitude depends on the amount of squeezing of the initial state. Strong attractive interactions are found to stabilize the relative number dynamics. We visualize the dynamics of the system in phase space using a quasiprobability distribution that allows for an intuitive interpretation of the various types of dynamics.

Temperature Dynamics Investigation At Small And Shallow Lakes Using Hydrodynamic Model

A three-dimensional time-dependent hydrodynamic and heat transport model of Lake Binaba, a shallow and small dam reservoir in Ghana, emphasizing the simulation of dynamics and thermal structure has been developed. Most numerical studies of temperature dynamics in reservoirs are based on one- or two-dimensional models. These models are not applicable for reservoirs characterized with complex flow pattern and unsteady heat exchange between the atmosphere and water surface. Continuity, momentum and temperature transport equations have been solved. Proper assignment of boundary conditions, especially surface heat fluxes, has been found crucial in simulating the lake’s hydrothermal dynamics. This model is based on the Reynolds Average Navier-Stokes equations, using a Boussinesq approach, with a standard k − ε turbulence closure to solve the flow field. The thermal model includes a heat source term, which takes into account the short wave radiation and also heat convection at the free surface, which is function of air temperatures, wind velocity and stability conditions of atmospheric boundary layer over the water surface. The governing equations of the model have been solved by OpenFOAM; an open source, freely available CFD toolbox. As its core, OpenFOAM has a set of efficient C++ modules that are used to build solvers. It uses collocated, polyhedral numerics that can be applied on unstructured meshes and can be easily extended to run in parallel. A new solver has been developed to solve the hydrothermal model of lake. The simulated temperature was compared against a 15 days field data set. Simulated and measured temperature profiles in the probe locations show reasonable agreement. The model might be able to compute total heat storage of water bodies to estimate evaporation from water surface.

Conformational features of cepacian: the exopolysaccharide produced by clinical strains of Burkholderia cepacia

Conformational energy calculations and molecular dynamics investigations, both in water and in dimethyl sulfoxide, were carried out on the exopolysaccharide cepacian produced by the majority of the clinical strains of Burkholderia cepacia, an opportunistic pathogen causing serious lung infection in patients affected by cystic fibrosis, the investigation was aimed at defining the structural and conformational features, which might be relevant for clarification of the structure-function relationships of the polymer. The molecular dynamics calculations were carried out by Ramachandran-type energy plots of the disaccharides that constitute the polymer repeating unit. The dynamics of an oligomer composed of three repeating units were investigated in water and in Me2SO, a non-aggregating solvent. Analysis of the time persistence of hydrogen bonds showed the presence of a large number of favourable interactions in water, which were less evident in Me2SO. The calculations on the cepacian chain indicated that polymer conformational features in water were affected by the lateral chains, but were also largely dictated by the presence of solvent. Moreover, the large number of intra-chain hydrogen bonds in water disappeared in Me2SO solution, increasing the average dimension of the polymer chains. (c) 2005 Elsevier Ltd. All rights reserved.

Mandibular skeletal growth and modelling between 10 and 15 years of age

This study pertains to a random sample of untreated French-Canadian adolescents (79 females and 107 males) evaluated at 10 and again at 15 years of age. Superimpositions on natural reference structures were performed to describe condylar growth and modelling of 11 mandibular landmarks. Superimpositions on natural cranial/cranial base reference structures were performed to describe mandibular displacement and true rotation.The results showed significant superior and posterior growth/modelling of the condyle and ramus. Males underwent significantly (P < 0.01) greater condylar growth and ramus modelling than females. With the exception of point B, which showed significant superior drift, modelling changes for the corpus landmarks were small and variable. The mandible rotated forward 2-3.3 degrees and was displaced 9.6-12.7 mm inferiorly and 1.9-2.7 mm anteriorly. Individual differences in ramus growth and modelling, both amount and direction, can be explained by mandibular rotation and displacements. Multivariate assessments revealed that superior condylar growth and ramus modelling were most closely associated with forward rotation and inferior mandibular displacement. Posterior growth and modelling were most closely correlated with anterior mandibular displacement and forward rotation. Modelling of the lower anterior border was independent of rotation and displacement.

Modelling frequency distributions of 5 minute-averaged solar radiation indexes using Beta probability functions

Five minute-averaged values of sky clearness, direct and diffuse indices, were used to model the frequency distributions of these variables in terms of optical air mass. From more than four years of solar radiation observations it was found that variations in the frequency distributions of the three indices of optical air mass for Botucatu, Brazil, are similar to those in other places, as published in the literature. The proposed models were obtained by linear combination of normalized Beta probability functions, using the observed distributions derived from three years of data. The versatility of these functions allows modelling of all three irradiance indexes to similar levels of accuracy. A comparison with the observed distributions obtained from one year of observations indicate that the models are able to reproduce the observed frequency distributions of all three indices at the 95% confidence level.

Dynamics of bone healing after osteotomy with piezosurgery or conventional drilling - histomorphometrical, immunohistochemical, and molecular analysis

Background: Piezosurgery is an osteotomy system used in medical and dental surgery. Many studies have proven clinical advantages of piezosurgery in terms of quality of cut, maneuverability, ease of use, and safety. However, few investigations have tested its superiority over the traditional osteotomy systems in terms of dynamics of bone healing. Therefore, the aim of this study was to evaluate the dynamics of bone healing after osteotomies with piezosurgery and to compare them with those associated to traditional bone drilling.Methods: One hundred and ten rats were divided into two groups with 55 animals each. The animals were anesthetized and the tibiae were surgically exposed to create defects 2 mm in diameter by using piezosurgery (Piezo group) and conventional drilling (Drill group). Animals were sacrificed at 3, 7, 14, 30 and 60 days post-surgery. Bone samples were collected and processed for histological, histomorphometrical, immunohistochemical, and molecular analysis. The histological analysis was performed at all time points (n = 8) whereas the histomorphometrical analysis was performed at 7, 14, 30 and 60 days post-surgery (n = 8). The immunolabeling was performed to detect Vascular Endothelial Growth Factor (VEGF), Caspase-3 (CAS-3), Osteoprotegerin (OPG), Receptor Activator of Nuclear Factor kappa-B Ligand (RANKL), and Osteocalcin (OC) at 3, 7, and 14 days (n = 3). For the molecular analysis, animals were sacrificed at 3, 7 and 14 days, total RNA was collected, and quantification of the expression of 21 genes related to BMP signaling, Wnt signaling, inflammation, osteogenenic and apoptotic pathways was performed by qRT-PCR (n = 5).Results: Histologically and histomorphometrically, bone healing was similar in both groups with the exception of a slightly higher amount of newly formed bone observed at 30 days after piezosurgery (p < 0.05). Immunohistochemical and qRT-PCR analyses didn't detect significant differences in expression of all the proteins and most of the genes tested.Conclusions: Based on the results of our study we conclude that in a rat tibial bone defect model the bone healing dynamics after piezosurgery are comparable to those observed with conventional drilling. © 2013 Esteves et al.; licensee BioMed Central Ltd.

Interpolating scattering amplitudes between the instant form and the front form of relativistic dynamics

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Scattering Amplitudes Interpolating Between Instant Form and Front Form of Relativistic Dynamics

Among the three forms of relativistic Hamiltonian dynamics proposed by Dirac in 1949, the front form has the largest number of kinematic generators. This distinction provides useful consequences in the analysis of physical observables in hadron physics. Using the method of interpolation between the instant form and the front form, we introduce the interpolating scattering amplitude that links the corresponding time-ordered amplitudes between the two forms of dynamics and provide the physical meaning of the kinematic transformations as they allow the invariance of each individual time-ordered amplitude for an arbitrary interpolation angle. We discuss the rationale for using front form dynamics, nowadays known as light-front dynamics (LFD), and present a few explicit examples of hadron phenomenology that LFD uniquely can offer from first-principles quantum chromodynamics. In particular, model-independent constraints are provided for the analyses of deuteron form factors and the N Delta transition form factors at large momentum transfer squared Q(2). The swap of helicity amplitudes between the collinear and non-collinear kinematics is also discussed in deeply virtual Compton scattering.

Modelling the maximum potential of nitrogen deposition and requirements of lysine for broilers

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Electromagnetic gauge field interpolation between the instant form and the front form of the Hamiltonian dynamics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Reduced-order modelling, circuit-level design and SOI fabrication of microelectromechanical resonators

This thesis deals with two important research aspects concerning radio frequency (RF) microresonators and switches. First, a new approach for compact modeling and simulation of these devices is presented. Then, a combined process flow for their simultaneous fabrication on a SOI substrate is proposed. Compact models for microresonators and switches are extracted by applying mathematical model order reduction (MOR) to the devices finite element (FE) description in ANSYS c° . The behaviour of these devices includes forms of nonlinearities. However, an approximation in the creation of the FE model is introduced, which enables the use of linear model order reduction. Microresonators are modeled with the introduction of transducer elements, which allow for direct coupling of the electrical and mechanical domain. The coupled system element matrices are linearized around an operating point and reduced. The resulting macromodel is valid for small signal analysis around the bias point, such as harmonic pre-stressed analysis. This is extremely useful for characterizing the frequency response of resonators. Compact modelling of switches preserves the nonlinearity of the device behaviour. Nonlinear reduced order models are obtained by reducing the number of nonlinearities in the system and handling them as input to the system. In this way, the system can be reduced using linear MOR techniques and nonlinearities are introduced directly in the reduced order model. The reduction of the number of system nonlinearities implies the approximation of all distributed forces in the model with lumped forces. Both for microresonators and switches, a procedure for matrices extraction has been developed so that reduced order models include the effects of electrical and mechanical pre-stress. The extraction process is fast and can be done automatically from ANSYS binary files. The method has been applied for the simulation of several devices both at devices and circuit level. Simulation results have been compared with full model simulations, and, when available, experimental data. Reduced order models have proven to conserve the accuracy of finite element method and to give a good description of the overall device behaviour, despite the introduced approximations. In addition, simulation is very fast, both at device and circuit level. A combined process-flow for the integrated fabrication of microresonators and switches has been defined. For this purpose, two processes that are optimized for the independent fabrication of these devices are merged. The major advantage of this process is the possibility to create on-chip circuit blocks that include both microresonators and switches. An application is, for example, aswitched filter bank for wireless transceiver. The process for microresonators fabrication is characterized by the use of silicon on insulator (SOI) wafers and on a deep reactive ion etching (DRIE) step for the creation of the vibrating structures in single-crystal silicon and the use of a sacrificial oxide layer for the definition of resonator to electrode distance. The fabrication of switches is characterized by the use of two different conductive layers for the definition of the actuation electrodes and by the use of a photoresist as a sacrificial layer for the creation of the suspended structure. Both processes have a gold electroplating step, for the creation of the resonators electrodes, transmission lines and suspended structures. The combined process flow is designed such that it conserves the basic properties of the original processes. Neither the performance of the resonators nor the performance of the switches results affected by the simultaneous fabrication. Moreover, common fabrication steps are shared, which allows for cheaper and faster fabrication.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Modelling coupled physical-biogeochemical processes in ice-covered oceans

The last decades have seen a large effort of the scientific community to study and understand the physics of sea ice. We currently have a wide - even though still not exhaustive - knowledge of the sea ice dynamics and thermodynamics and of their temporal and spatial variability. Sea ice biogeochemistry is instead largely unknown. Sea ice algae production may account for up to 25% of overall primary production in ice-covered waters of the Southern Ocean. However, the influence of physical factors, such as the location of ice formation, the role of snow cover and light availability on sea ice primary production is poorly understood. There are only sparse localized observations and little knowledge of the functioning of sea ice biogeochemistry at larger scales. Modelling becomes then an auxiliary tool to help qualifying and quantifying the role of sea ice biogeochemistry in the ocean dynamics. In this thesis, a novel approach is used for the modelling and coupling of sea ice biogeochemistry - and in particular its primary production - to sea ice physics. Previous attempts were based on the coupling of rather complex sea ice physical models to empirical or relatively simple biological or biogeochemical models. The focus is moved here to a more biologically-oriented point of view. A simple, however comprehensive, physical model of the sea ice thermodynamics (ESIM) was developed and coupled to a novel sea ice implementation (BFM-SI) of the Biogeochemical Flux Model (BFM). The BFM is a comprehensive model, largely used and validated in the open ocean environment and in regional seas. The physical model has been developed having in mind the biogeochemical properties of sea ice and the physical inputs required to model sea ice biogeochemistry. The central concept of the coupling is the modelling of the Biologically-Active-Layer (BAL), which is the time-varying fraction of sea ice that is continuously connected to the ocean via brines pockets and channels and it acts as rich habitat for many microorganisms. The physical model provides the key physical properties of the BAL (e.g., brines volume, temperature and salinity), and the BFM-SI simulates the physiological and ecological response of the biological community to the physical enviroment. The new biogeochemical model is also coupled to the pelagic BFM through the exchange of organic and inorganic matter at the boundaries between the two systems . This is done by computing the entrapment of matter and gases when sea ice grows and release to the ocean when sea ice melts to ensure mass conservation. The model was tested in different ice-covered regions of the world ocean to test the generality of the parameterizations. The focus was particularly on the regions of landfast ice, where primary production is generally large. The implementation of the BFM in sea ice and the coupling structure in General Circulation Models will add a new component to the latters (and in general to Earth System Models), which will be able to provide adequate estimate of the role and importance of sea ice biogeochemistry in the global carbon cycle.