Mixed-dimensional coupling in finite element models

In many finite element analysis models it would be desirable to combine reduced- or lower-dimensional element types with higher-dimensional element types in a single model. In order to achieve compatibility of displacements and stress equilibrium at the junction or interface between the differing element types, it is important in such cases to integrate into the analysis some scheme for coupling the element types. A novel and effective scheme for establishing compatibility and equilibrium at the dimensional interface is described and its merits and capabilities are demonstrated. Copyright (C) 2000 John Wiley & Sons, Ltd.

Atomistic Insights into Rhodopsin Activation from a Dynamic Model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Modelling the mechanical response of the Reed-mouthpiece-lip system of a clarinet. Part I. A one-dimensional distributed model

The motion of a clarinet reed that is clamped to a mouthpiece and supported by a lip is simulated in the time-domain using finite difference methods. The reed is modelled as a bar with non-uniform cross section, and is described using a one-dimensional, fourth-order partial differential equation. The interactions with the mouthpiece Jay and the player's lip are taken into account by incorporating conditional contact forces in the bar equation. The model is completed by clamped-free boundary conditions for the reed. An implicit finite difference method is used for discretising the system, and values for the physical parameters are chosen both from laboratory measurements and by accurate tuning of the numerical simulations. The accuracy of the numerical system is assessed through analysis of frequency warping effects and of resonance estimation. Finally, the mechanical properties of the system are studied by analysing its response to external driving forces. In particular, the effects of reed curling are investigated.

Modelling the mechanical response of the reed-mouthpiece-lip system of a clarinet. Part II: A lumped model approximation

A non-linear lumped model of the reed-mouthpiece-lip system of a clarinet is formulated, in which the lumped parameters are derived from numerical experiments with a finite-difference simulation based on a distributed reed model. The effective stiffness per unit area is formulated as a function of the pressure signal driving the reed, in order to simulate the effects of the reed bending against the lay, and mass and damping terms are added as a first approximation to the dynamic behaviour of the reed. A discrete-time formulation is presented, and its response is compared to that of the distributed model. In addition, the lumped model is applied in the simulation of clarinet tones, enabling the analysis of the effects of using a pressure-dependent stiffness per unit area on sustained oscillations. The analysed effects and features are in qualitative agreement with players' experiences and experimental results obtained in prior studies.

Passive intermodulation in finite lengths of printed microstrip lines

This paper addresses the theoretical aspects of passive intermodulation (PIM) generation in printed transmission lines. In order to elucidate the mechanisms of PIM generation, a new model of the transmission line length with distributed nonlinearity is proposed. The developed model has been validated by the near-field measurements of PIM product distributions along the microstrip lines. The contributions of nonlinear mixing, power dissipation, and load matching to PIM products have been analyzed in detail. The obtained results reveal the fundamental properties of PIM generation in finite lengths of printed lines with distributed non-linearity and identify possible means for PIM mitigation. It was shown for the first time that the reverse PIM products in a matched transmission line with distributed nonlinearity are generated due to nonlinear scattering. © 2008 IEEE.

An explanatory model for food-web structure and evolution

Food webs are networks describing who is eating whom in an ecological community. By now it is clear that many aspects of food-web structure are reproducible across diverse habitats, yet little is known about the driving force behind this structure. Evolutionary and population dynamical mechanisms have been considered. We propose a model for the evolutionary dynamics of food-web topology and show that it accurately reproduces observed food-web characteristics in the steady state. It is based on the observation that most consumers are larger than their resource species and the hypothesis that speciation and extinction rates decrease with increasing body mass. Results give strong support to the evolutionary hypothesis. (C) 2005 Elsevier B.V. All rights reserved.

Dynamic pricing model and algorithm for perishable products with fuzzy demand (ABS 2* Journal)

This paper studies the dynamic pricing problem of selling fixed stock of perishable items over a finite horizon, where the decision maker does not have the necessary historic data to estimate the distribution of uncertain demand, but has imprecise information about the quantity demand. We model this uncertainty using fuzzy variables. The dynamic pricing problem based on credibility theory is formulated using three fuzzy programming models, viz.: the fuzzy expected revenue maximization model, a-optimistic revenue maximization model, and credibility maximization model. Fuzzy simulations for functions with fuzzy parameters are given and embedded into a genetic algorithm to design a hybrid intelligent algorithm to solve these three models. Finally, a real-world example is presented to highlight the effectiveness of the developed model and algorithm.

Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

Exact theory for localized envelope modulated electrostatic wavepackets in space and dusty plasmas

Abundant evidence for the occurrence of modulated envelope plasma wave packets is provided by recent satellite missions. These excitations are characterized by a slowly varying localized envelope structure, embedding the fast carrier wave, which appears to be the result of strong modulation of the wave amplitude. This modulation may be due to parametric interactions between different modes or, simply, to the nonlinear (self-)interaction of the carrier wave. A generic exact theory is presented in this study, for the nonlinear self-modulation of known electrostatic plasma modes, by employing a collisionless fluid model. Both cold (zero-temperature) and warm fluid descriptions are discussed and the results are compared. The (moderately) nonlinear oscillation regime is investigated by applying a multiple scale technique. The calculation leads to a Nonlinear Schrodinger-type Equation (NLSE), which describes the evolution of the slowly varying wave amplitude in time and space. The NLSE admits localized envelope (solitary wave) solutions of bright(pulses) or dark- (holes, voids) type, whose characteristics (maximum amplitude, width) depend on intrinsic plasma parameters. Effects like amplitude perturbation obliqueness (with respect to the propagation direction), finite temperature and defect (dust) concentration are explicitly considered. Relevance with similar highly localized modulated wave structures observed during recent satellite missions is discussed.

Strong-field approximation for Coulomb explosion of H-2(+) and D-2(+) by short intense laser pulses

We present a simple quantum mechanical model to describe Coulomb explosion of H-2(+) and D-2(+) by short, intense infrared laser pulses. The model is based on the length gauge version of the molecular strong-field approximation and is valid when the process of dissociation prior to ionization is negligible. The results are compared with recent experimental data for the proton kinetic energy spectrum [Th. Ergler , Phys. Rev. Lett. 95, 093001 (2005); D. S. Murphy , J. Phys. B 40, S359 (2007)]. Using a Franck-Condon distribution over initial vibrational states, the theory reproduces the overall shape of the spectrum with only a small overestimation of slow protons. The agreement between theory and experiment can be made perfect by using a non-Frank-Condon initial distribution characteristic for H-2(+) (D-2(+)) targets produced by strong-field ionization of H-2 (D-2). For comparison, we also present results obtained by two different tunneling models for this process.

Ambiguity, social opinion and the use of common property resources

In this paper we argue that ambiguity, combined with social opinion formation, can be used as the foundation of a game-theoretic equilibrium concept that transcends the standard Nash equilibrium concept, applied to a model of the tragedy of the commons. Our approach sheds light on the international environmental crisis and the relevant ongoing international negotiations. We conclude that social opinion formation in most cases has a significant impact on equilibrium common property resource usage.

Redistribution of vibrational population in a molecular ion with nonresonant strong-field laser pulses

We present an experimental demonstration of nonresonant manipulation of vibrational states in a molecule by an intense ultrashort laser pulse. A vibrational wave packet is generated in D-2(+) through tunnel ionization of D-2 by a few-cycle pump pulse. A similar control pulse is applied as the wave packet begins to dephase so that the dynamic Stark effect distorts the electronic environment of the nuclei, transferring vibrational population. The time evolution of the modified wave packet is probed via the D-2(+) photodissociation yield that results from the application of an intense probe pulse. Comparing the measured yield with a quasiclassical trajectory model allows us to determine the redistribution of vibrational population caused by the control pulse. ©

Finite difference time domain modeling of phase grating diffusion

In this paper, a method for modeling diffusion caused by non-smooth boundary surfaces in simulations of room acoustics using finite difference time domain (FDTD) technique is investigated. The proposed approach adopts the well-known theory of phase grating diffusers to efficiently model sound scattering from rough surfaces. The variation of diffuser well-depths is attained by nesting allpass filters within the reflection filters from which the digital impedance filters used in the boundary implementation are obtained. The presented technique is appropriate for modeling diffusion at high frequencies caused by small surface roughness and generally diffusers that have narrow wells and infinitely thin separators. The diffusion coefficient was measured with numerical experiments for a range of fractional Brownian diffusers.

Modal Representation of the Resonant Body within a Finite Difference Framework for Simulation of String Instruments

This paper investigates numerical simulation of a string coupled
transversely to a resonant body. Starting from a complete nite
difference formulation, a second model is derived in which the
body is represented in modal form. The main advantage of this hybrid form is that the body model is scalable, i.e. the computational
complexity can be adjusted to the available processing power. Numerical results are calculated and discussed for simplied models
in the form of string-string coupling and string-plate coupling.

How can cells sense the elasticity of a substrate? An analysis using a cell tensegrity model

A eukaryotic cell attaches and spreads on substrates, whether it is the extracellular matrix naturally produced by the cell itself, or artificial materials, such as tissue-engineered scaffolds. Attachment and spreading require the cell to apply forces in the nN range to the substrate via adhesion sites, and these forces are balanced by the elastic response of the substrate. This mechanical interaction is one determinant of cell morphology and, ultimately, cell phenotype. In this paper we use a finite element model of a cell, with a tensegrity structure to model the cytoskeleton of actin filaments and microtubules, to explore the way cells sense the stiffness of the substrate and thereby adapt to it. To support the computational results, an analytical 1D model is developed for comparison. We find that (i) the tensegrity hypothesis of the cytoskeleton is sufficient to explain the matrix-elasticity sensing, (ii) cell sensitivity is not constant but has a bell-shaped distribution over the physiological matrix-elasticity range, and (iii) the position of the sensitivity peak over the matrix-elasticity range depends on the cytoskeletal structure and in particular on the F-actin organisation. Our model suggests that F-actin reorganisation observed in mesenchymal stem cells (MSCs) in response to change of matrix elasticity is a structural-remodelling process that shifts the sensitivity peak towards the new value of matrix elasticity. This finding discloses a potential regulatory role of scaffold stiffness for cell differentiation.