In the present work the simulation of water impacts is discussed. The investigation is mainly focused on the energy dissipation involved in liquid impacts in both the frameworks of the weakly compressible and incompressible models. A detailed analysis is performed using a weakly compressible Smoothed Particle Hydrodynamics (SPH) solver and the results are compared with the solutions computed by an incompressible mesh-based Level-Set Finite Volume Method (LS-FVM). Impacts are numerically studied using single-phase models through prototypical problems in 1D and 2D frameworks. These problems were selected for the conclusions to be of interest for, e.g., the numerical computation of the flow around plunging breaking waves. The conclusions drawn are useful not only to SPH or LS-FVM users but also for other numerical models, for which accurate results on benchmark test-cases are provided.
Energy loss
Finite Volume Method
Free-surface flow
Level-Set
Smoothed Particle Hydrodynamics
Water impact
MIPAS measurements on ENVISAT represent a unique database for the study of atmospheric composition and of the time variation of atmospheric constituents.
For trend studies it is important that instrumental drifts are reduced. Some of the MIPAS spectral bands are affected by time-dependent non-linearity which have been recently corrected.
In addition to this non-linearity correction, the forthcoming new version of MIPAS products (Version 7) contains several other improvements: a new approach for handling continuum capable of making the retrieval more stable, a new selection of spectral intervals selected for the analysis of the full resolution measurements aimed to reduce the bias between full resolution and optimized resolution measurements, the regularization of the H2O profiles, the products of five new species (HCFC-22, CFC-14, HCN, COF2, CCl4), which makes equal to 15 the total number of retrieved species.
The latest improvements implemented in the ESA processor and some of the results on trend will be presented and discussed.
Eugenio Elia Levi è stato uno dei più grandi matematici italiani del '900 (come del resto il fratello Beppo). Morì nell'ottobre del 1917, ucciso da un cecchino, nelle fasi iniziali della disfatta di Caporetto.
Attivo interventista, allo scoppio della prima guerra mondiale si era arruolato volontariamente nel Genio Zappatori e fu poi promosso capitano per meriti di guerra. Eugenio Elia Levi si era laureato in Matematica nel 1904, studiando alla "Normale" di Pisa dove ebbe come maestri Luigi Bianchi e Ulisse Dini. Nel 1909 ottenne la cattedra di Analisi infinitesimale presso l'Università di Genova. La sua produzione scientifica fu tanto profonda quanto differenziata - ha riguardato i gruppi di Lie, le equazioni alle derivate parziali e la teoria delle funzioni di più variabili complesse - e venne immediatamente apprezzata negli ambienti matematici internazionali.
Benché fosse di soli due anni più anziano di Mauro Picone (fondatore dell'IAC), quest'ultimo lo considerò sempre come un suo grande maestro, alla pari di Luigi Bianchi e Ulisse Dini.
MIPAS measurements on ENVISAT represent a unique database for the study of atmospheric composition and of the time variation of atmospheric constituents.For trend studies it is important that instrumental drifts are reduced. Some of the MIPAS spectral bands are affected by time-dependent non-linearity that have been recently corrected.In addition to this non-linearity correction, the forthcoming new version of MIPAS products (Version 7) contains several other improvements: a new approach for handling continuum leading to a more stable retrieval, a new selection of spectral intervals for the analysis of the full resolution measurements aiming to reduce the bias between full resolution and optimized resolution measurements, the regularization of the H2O profiles. Furthermore, the implementation of the retrieval of five new species (HCFC-22, CFC-14, HCN, COF2, CCl4) leads to a total of 15 species in ESA products.The latest improvements implemented in the ESA processor, the results of the validation of the products and some preliminary results on the trend of some ozone depleting substances will be presented and discussed.
A sharp integrability condition on the right-hand side of the p-Laplace system for all its solutions
to be continuous is exhibited. Their uniform continuity is also analyzed and estimates for their modulus
of continuity are provided. The relevant estimates are shown to be optimal as the right-hand side ranges in
classes of rearrangement-invariant spaces, such as Lebesgue, Lorentz, Lorentz-Zygmund, and Marcinkiewicz
spaces, as well as some customary Orlicz spaces.
Nonlinear elliptic systems
continuity of solutions
modulus of continuity
classical Lorentz spaces
Orlicz spaces
Sobolev embeddings.
We obtain a comparison result for solutions to nonlinear fully anisotropic elliptic problems by means of
anisotropic symmetrization. As consequence we deduce a priori estimates for norms of the relevant solutions.
Anisotropic symmetrization rearrangements
A priori estimates
Dirichlet problems.
The hybrid use of exact and heuristic derivative-free methods for global unconstrained optimization problems is presented. Many real-world problems are modeled by computationally expensive functions, such as problems in simulationbased design of complex engineering systems. Objective-function values are often provided by systems of partial differential equations, solved by computationally expensive black-box tools. The objective-function is likely noisy and its derivatives are often not available. On the one hand, the use of exact optimization methods might be computationally too expensive, especially if asymptotic convergence properties are sought. On the other hand, heuristic methods do not guarantee the stationarity of their final solutions. Nevertheless, heuristic methods are usually able to provide an approximate solution at a reasonable computational cost, and have been widely applied to real-world simulation-based design optimization problems. Herein, an overall hybrid algorithm combining the appealing properties of both exact and heuristic methods is discussed, with focus on Particle Swarm Optimization (PSO) and line search-based derivative-free algorithms. The theoretical properties of the hybrid algorithm are detailed, in terms of limit points stationarity. Numerical results are presented for a specific test function and for two real-world optimization problems in ship hydrodynamics.
Derivative-free optimization
Global optimization
Hybrid optimization algorithm
Line search algorithm
Particle swarm optimization
Ship design
Simulationbased design
The book describes a computational model of the immune system reaction, C-ImmSim, built along the lines of the computer model known as the Celada-Seiden model (CS-model). The computational counterpart of the CS-model is called IMMSIM which stands for IMMune system SIMulator. IMMSIM was written in 1992 by the physicist Phil E. Seiden and the immunologist Franco Celada. This model was built around the idea of developing a computerized system to perform experiments similar in vivo experiments; a tool developed to help biologists testing theories and hypothesis about how the immune system works.
C-ImmSim is best viewed as a collection of models in a single program. It incorporates the principal core facts of today's immunological knowledge, such as the diversity of specific elements, MHC restriction, clonal selection, thymic education of T cells, antigen processing and presentation (both the cytosolic and endocytic pathways are implemented), cell-cell cooperation, homeostasis of cells created by the bone marrow, hyper mutation of antibodies, maturation of the cellular and humoral response, and memory. Besides, an antigen can represent a bacterium, a virus, or an allergen or a tumor cell. C-ImmSim has been recently customized to simulate the HIV-1 infection. Moreover, it can simulate the immunotherapy for cancer. These features are all present in the code and people can choose to turn them on and off at compiling time.
The book presents the basic model as well as the various customizations to implement the description of different diseases and the way they have been used in practice to produce new knowledge either from hypothesis or from lab-experiment data. In this respect, the book can be used as a practical guide to implement a computational model with which to study a specific disease and to try to address realistic clinical questions.
TRANSLATION OF THE ORIGINAL DOCUMENT "NETWORK LITERACY" AVAILABLE AT https://sites.google.com/a/binghamton.edu/netscied/teaching-learning/network-concepts
Mentre il nostro mondo diventa sempre più connesso attraverso l'uso di reti, o network, che rendono le comunicazioni e la diffusione di informazioni pressoché istantanee, il livello di comprensione di come queste reti funzionino avrà un ruolo importante nel determinare quanto la società trarrà beneficio da questa connettività accresciuta.
In breve, una società connessa richiede un'alfabetizzazione sul concetto di rete, ossia una conoscenza di base di cosa siano le reti, come possano essere utilizzate come strumento per la scoperta e per i processi decisionali, nonché sulle problematiche e i potenziali vantaggi resi accessibili a tutti coloro che vivono nel mondo interconnesso di oggi.
In this paper we present a general model of drug release from a drug delivery device and the subsequent transport in biological tissue. The model incorporates drug diffusion, dissolution and solubility in the polymer coating, coupled with diffusion, convection and reaction in the biological tissue. Each layer contains bound and free drug phases so that the resulting model is a coupled two-phase two-layer system of partial differential equations. One of the novelties is the generality of the model in each layer. Within the drug coating, our model includes diffusion as well as three different models of dissolution. We show that the model may also be used in cases where dissolution is rapid or not relevant, and additionally when drug release is not limited by its solubility. Within the biological tissue, the model can account for nonlinear saturable reversible binding, with linear reversible binding and linear irreversible binding being recovered as special cases. The generality of our model will allow the simulation of the release from a wide range of drug delivery devices encompassing many different applications. To demonstrate the efficacy of our model we simulate results for the particular application of drug release from arterial stents.
Program summary
We present the open-source computer program JETSPIN, specifically designed to simulate the electro-spinning process of nanofibers. Its capabilities are shown with proper reference to the underlying model, as well as a description of the relevant input variables and associated test-case simulations. The various interactions included in the electrospinning model implemented in JETSPIN are discussed in detail. The code is designed to exploit different computational architectures, from single to parallel processor workstations. This paper provides an overview of JETSPIN, focusing primarily on its structure, parallel implementations, functionality, performance, and availability.
Electrospinning
Jet dynamics
Viscoelasticity
Nanofibers
Coarse grained model
Lagrangian model
We investigate the effects of dissipative air drag on the dynamics of electrified jets in the initial stage of the electrospinning process. The main idea is to use a Brownian noise to model air drag effects on the uniaxial elongation of the jets. The developed numerical model is used to probe the dynamics of electrified polymer jets at different conditions of air drag force, showing that the dynamics of the charged jet is strongly biased by the presence of air drag forces. This study provides prospective beneficial implications for improving forthcoming electrospinning experiments. (C) 2015 Elsevier Ltd. All rights reserved.
Electrospinning
Air drag
Viscoelasticity
Nanofibers
We present a nonlinear Langevin model to investigate the early-stage dynamics of electrified polymer jets in electrospinning experiments. In particular, we study the effects of air drag force on the uniaxial elongation of the charged jet, right after ejection from the nozzle. Numerical simulations show that the elongation of the jet filament close to the injection point is significantly affected by the nonlinear drag exerted by the surrounding air. These results provide useful insights for the optimal design of current and future electrospinning experiments.
In this paper we present a model of drug release from a drug eluting-stent and the subsequent drug transport in the arterial wall. In order to study the complete process, a two-phase mathematical model describing the transport of a drug between two coupled media of different properties and dimensions is presented. A system of partial differential equations describes both the solid-liquid transfer (dissolution) and diffusion processes in the polymeric substrate as well as diffusion, convection and reaction in the tissue layer. We adopt a semi-analytical approach in solving the four resulting equations. The model serves as a useful tool for providing insight into the influence of solid-liquid mass transfer in drug-eluting stent systems.
Drug-Eluting Stents
Drug Delivery
Solid-Liquid Mass Transfer
Reaction-Diffusion-Convection Equations
Two-Phase Model
Coupled Partial Differential Equations
Modelling the glycocalyx-endothelium-erythrocyte interaction in the microcirculation: a computational study
Pontrelli Giuseppe
;
Halliday Ian
;
Spencer Tim J
;
Koenig Carola S
;
Collins Michael W
A novel, coarse-grained, single-framework 'Eulerian' model for blood flow in the microvascular circulation is presented and used to estimate the variations in flow properties that accrue from all of the following: (i) wall position variation, associated with the endothelial cells' (ECs) shape, (ii) glycocalyx layer (GL) effects and (iii) the particulate nature of blood. We stress that our new model is fully coupled and uses only a single Eulerian computational framework to recover complex effects, dispensing altogether with the need for, e.g. re-meshing and advected sets of Lagrangian points. Physically, blood is modelled as a two-component, incompressible fluid - the plasma and corpuscular elements dispersed in it. The latter are modelled as deformable liquid droplets of increased viscosity. Interfacial membrane effects are present to mimic key blood properties and to avoid droplets' coalescence. The model is encapsulated within a multi-component lattice Boltzmann method that uses a sub-lattice 'wavy wall' closure to represent the ECs. Between this boundary and the flow domain, the model incorporates a coarse-grained representation of the endothelial GL, which is known to cover microvessel walls. The endothelial glycocalyx is modelled as a medium of variable and adaptive porosity, with approaching droplets being subject to a repulsive elastic force. Numerical simulations are presented to show the combined and simultaneous influence on fundamental flow properties of the EC wall undulation, the glycocalyx compression and repulsion and the particulate nature of blood. Several characteristic hemodynamical features of microvessel flow are successfully reproduced, including the deformability of particulates and the Fahraeus-Lindqvist effect. Moreover, the importance of modelling the GL is manifest in the magnitude of and the temporal variations in the flow rate and wall shear stresses.
In this paper we propose a discrete in continuous mathematical model for the morphogenesis of the posterior lateral line system in zebrafish. Our model follows closely the results obtained in recent biological experiments. We rely on a hybrid description: discrete for the cellular level and continuous for the molecular level. We prove the existence of steady solutions consistent with the formation of particular biological structure, the neuromasts. Dynamical numerical simulations are performed to show the behavior of the model and its qualitative and quantitative accuracy to describe the evolution of the cell aggregate.
We consider two models which were both designed to describe the movement of eukaryotic cells responding to chemical signals. Besides a common standard parabolic equation for the diffusion of a chemoattractant, like chemokines or growth factors, the two models differ for the equations describing the movement of cells. The first model is based on a quasilinear hyperbolic system with damping, the other one on a degenerate parabolic equation. The two models have the same stationary solutions, which may contain some regions with vacuum. We first explain in details how to discretize the quasilinear hyperbolic system through an upwinding technique, which uses an adapted reconstruction, which is able to deal with the transitions to vacuum. Then we concentrate on the analysis of asymptotic preserving properties of the scheme towards a discretization of the parabolic equation, obtained in the large time and large damping limit, in order to present a numerical comparison between the asymptotic behavior of these two models. Finally we perform an accurate numerical comparison of the two models in the time asymptotic regime, which shows that the respective solutions have a quite different behavior for large times.
In this paper we study a semilinear hyperbolic-parabolic system modeling biological phenomena evolving on a network composed of oriented arcs. We prove the existence of global (in time) smooth solutions to this problem. The result is obtained by using energy estimates with suitable transmission conditions at nodes.
nonlinear hyperbolic systems
networks
transmission conditions
global existence of solutions
chemotaxis