Stability of a non-local kinetic model for cell migration with density-dependent speed

2020 ◽  
Author(s):  
Nadia Loy ◽  
Luigi Preziosi
Keyword(s):  
Kinetic Model ◽  
One Dimensional ◽  
Dirac Delta ◽  
Finite Distance ◽  
Non Local ◽  
Sensing Radius ◽  
The Mean ◽  
The Stability ◽  
The One ◽  
Direction Of Polarization

Abstract The aim of this article is to study the stability of a non-local kinetic model proposed by Loy & Preziosi (2020a) in which the cell speed is affected by the cell population density non-locally measured and weighted according to a sensing kernel in the direction of polarization and motion. We perform the analysis in a $d$-dimensional setting. We study the dispersion relation in the one-dimensional case and we show that the stability depends on two dimensionless parameters: the first one represents the stiffness of the system related to the cell turning rate, to the mean speed at equilibrium and to the sensing radius, while the second one relates to the derivative of the mean speed with respect to the density evaluated at the equilibrium. It is proved that for Dirac delta sensing kernels centered at a finite distance, corresponding to sensing limited to a given distance from the cell center, the homogeneous configuration is linearly unstable to short waves. On the other hand, for a uniform sensing kernel, corresponding to uniformly weighting the information collected up to a given distance, the most unstable wavelength is identified and consistently matches the numerical solution of the kinetic equation.

scholarly journals Stability of a one-dimensional morphoelastic model for post-burn contraction

2021 ◽  
Vol 83 (3) ◽  
Author(s):  
Ginger Egberts ◽  
Fred Vermolen ◽  
Paul van Zuijlen
Keyword(s):  
Wound Healing ◽  
Residual Stresses ◽  
Truncation Error ◽  
Signaling Molecules ◽  
Discrete Problem ◽  
One Dimensional ◽  
Stability Constraint ◽  
Biological Interpretation ◽  
The Stability ◽  
The One

AbstractTo deal with permanent deformations and residual stresses, we consider a morphoelastic model for the scar formation as the result of wound healing after a skin trauma. Next to the mechanical components such as strain and displacements, the model accounts for biological constituents such as the concentration of signaling molecules, the cellular densities of fibroblasts and myofibroblasts, and the density of collagen. Here we present stability constraints for the one-dimensional counterpart of this morphoelastic model, for both the continuous and (semi-) discrete problem. We show that the truncation error between these eigenvalues associated with the continuous and semi-discrete problem is of order $${{\mathcal {O}}}(h^2)$$ O ( h 2 ) . Next we perform numerical validation to these constraints and provide a biological interpretation of the (in)stability. For the mechanical part of the model, the results show the components reach equilibria in a (non) monotonic way, depending on the value of the viscosity. The results show that the parameters of the chemical part of the model need to meet the stability constraint, depending on the decay rate of the signaling molecules, to avoid unrealistic results.

scholarly journals One-Dimensional Organic–Inorganic Material (C6H9N2)2BiCl5: From Synthesis to Structural, Spectroscopic, and Electronic Characterizations

2021 ◽  
Vol 22 (4) ◽  
pp. 2030
Author(s):  
Hela Ferjani ◽  
Hammouda Chebbi ◽  
Mohammed Fettouhi
Keyword(s):  
Hydrogen Bonding ◽  
Density Functional ◽  
Crystal Packing ◽  
Diffuse Reflectance Spectrum ◽  
Enrichment Ratio ◽  
One Dimensional ◽  
Organic Part ◽  
Fukui Indices ◽  
The Stability ◽  
The One

The new organic–inorganic compound (C6H9N2)2BiCl5 (I) has been grown by the solvent evaporation method. The one-dimensional (1D) structure of the allylimidazolium chlorobismuthate (I) has been determined by single crystal X-ray diffraction. It crystallizes in the centrosymmetric space group C2/c and consists of 1-allylimidazolium cations and (1D) chains of the anion BiCl52−, built up of corner-sharing [BiCl63−] octahedra which are interconnected by means of hydrogen bonding contacts N/C–H⋯Cl. The intermolecular interactions were quantified using Hirshfeld surface analysis and the enrichment ratio established that the most important role in the stability of the crystal structure was provided by hydrogen bonding and H···H interactions. The highest value of E was calculated for the contact N⋯C (6.87) followed by C⋯C (2.85) and Bi⋯Cl (2.43). These contacts were favored and made the main contribution to the crystal packing. The vibrational modes were identified and assigned by infrared and Raman spectroscopy. The optical band gap (Eg = 3.26 eV) was calculated from the diffuse reflectance spectrum and showed that we can consider the material as a semiconductor. The density functional theory (DFT) has been used to determine the calculated gap, which was about 3.73 eV, and to explain the electronic structure of the title compound, its optical properties, and the stability of the organic part by the calculation of HOMO and LUMO energy and the Fukui indices.

A TWELFTH-ORDER FOUR-STEP FORMULA FOR THE NUMERICAL INTEGRATION OF THE ONE-DIMENSIONAL SCHRÖDINGER EQUATION

2003 ◽  
Vol 14 (08) ◽  
pp. 1087-1105 ◽  
Author(s):  
ZHONGCHENG WANG ◽  
YONGMING DAI
Keyword(s):  
Schrödinger Equation ◽  
Numerical Integration ◽  
Schrodinger Equation ◽  
Numerical Test ◽  
New Method ◽  
Step Method ◽  
One Dimensional ◽  
The Stability ◽  
The One ◽  
Order Four

A new twelfth-order four-step formula containing fourth derivatives for the numerical integration of the one-dimensional Schrödinger equation has been developed. It was found that by adding multi-derivative terms, the stability of a linear multi-step method can be improved and the interval of periodicity of this new method is larger than that of the Numerov's method. The numerical test shows that the new method is superior to the previous lower orders in both accuracy and efficiency and it is specially applied to the problem when an increasing accuracy is requested.

Stability of Viscoelastic Channel Flow

2007 ◽  
Author(s):  
Nariman Ashrafi
Keyword(s):  
Reynolds Number ◽  
Viscosity Ratio ◽  
Weissenberg Number ◽  
Base Flow ◽  
Galerkin Projection ◽  
Large Reynolds Number ◽  
One Dimensional ◽  
Stress Effects ◽  
The Stability ◽  
The One

The nonlinear stability and bifurcation of the one-dimensional channel (Poiseuille) flow is examined for a Johnson-Segalman fluid. The velocity and stress are represented by orthonormal functions in the transverse direction to the flow. The flow field is obtained from the conservation and constitutive equations using the Galerkin projection method. Both inertia and normal stress effects are included. The stability picture is dramatically influenced by the viscosity ratio. The range of shear rate or Weissenberg number for which the base flow is unstable increases from zero as the fluid deviates from the Newtonian limit as decreases. Typically, two turning points are observed near the critical Weissenberg numbers. The transient response is heavily influenced by the level of inertia. It is found that the flow responds oscillatorily. When the Reynolds number is small, and monotonically at large Reynolds number when elastic effects are dominated by inertia.

scholarly journals Assessment of CFD turbulence models for free surface flow simulation and 1-D modelling for water level calculations over a broad-crested weir floodway

Water SA ◽  
2019 ◽  
Vol 45 (3 July) ◽  
Author(s):  
Ahmed M Helmi
Keyword(s):  
Experimental Data ◽  
Free Surface ◽  
Water Level ◽  
Dimensional Analysis ◽  
Turbulence Model ◽  
Turbulence Models ◽  
Cfd Simulations ◽  
One Dimensional ◽  
The Mean ◽  
The One

Floodways, where a road embankment is permitted to be overtopped by flood water, are usually designed as broad-crested weirs. Determination of the water level above the floodway is crucial and related to road safety. Hydraulic performance of floodways can be assessed numerically using 1-D modelling or 3-D simulation using computational fluid dynamics (CFD) packages. Turbulence modelling is one of the key elements in CFD simulations. A wide variety of turbulence models are utilized in CFD packages; in order to identify the most relevant turbulence model for the case in question, 96 3-D CFD simulations were conducted using Flow-3D package, for 24 broad-crested weir configurations selected based on experimental data from a previous study. Four turbulence models (one-equation, k-ε, RNG k-ε, and k-ω) ere examined for each configuration. The volume of fluid (VOF) algorithm was adopted for free water surface determination. In addition, 24 1-D simulations using HEC-RAS-1-D were conducted for comparison with CFD results and experimental data. Validation of the simulated water free surface profiles versus the experimental measurements was carried out by the evaluation of the mean absolute error, the mean relative error percentage, and the root mean square error. It was concluded that the minimum error in simulating the full upstream to downstream free surface profile is achieved by using one-equation turbulence model with mixing length equal to 7% of the smallest domain dimension. Nevertheless, for the broad-crested weir upstream section, no significant difference in accuracy was found between all turbulence models and the one-dimensional analysis results, due to the low turbulence intensity at this part. For engineering design purposes, in which the water level is the main concern at the location of the flood way, the one-dimensional analysis has sufficient accuracy to determine the water level.

Microscale temperature and velocity spectra in the atmospheric boundary layer

1977 ◽  
Vol 83 (3) ◽  
pp. 547-567 ◽  
Author(s):  
R. M. Williams ◽  
C. A. Paulson
Keyword(s):  
Reynolds Number ◽  
Prandtl Number ◽  
Temperature Fluctuations ◽  
Inertial Subrange ◽  
Velocity Fluctuations ◽  
One Dimensional ◽  
Temperature Spectra ◽  
The Mean ◽  
Velocity Spectra ◽  
The One

High-frequency fluctuations in temperature and velocity were measured at a height of 2 m above a harvested, nearly level field of rye grass. Conditions were both stably and unstably stratified. Reynolds numbers ranged from 370000 to 740000. Measurements of velocity were made with a hot-wire anemometer and measurements of temperature with a platinum resistance element which had a diameter of 0[sdot ]5 μm and a length of 1 mm. Thirteen runs ranging in length from 78 to 238 s were analysed.Spectra of velocity fluctuations are consistent with previously reported universal forms. Spectra of temperature, however, exhibit an increase in slope with increasing wavenumber as the maximum in the one-dimensional dissipation spectrum is approached. The peak of the one-dimensional dissipation spectrum for temperature fluctuations occurs at a higher wavenumber than that of simultaneous spectra of the dissipation of velocity fluctuations. It is suggested that the change in slope of the temperature spectra and the dissimilarity between temperature and velocity spectra may be due to spatial dissimilarity in the dissipation of temperature and velocity fluctuations. The temperature spectra are compared with a theoretical prediction for fluids with large Prandtl number, due to Batchelor (1959). Even though air has a Prandtl number of 0[sdot ]7, the observations are in qualitative agreement with predictions of the theory. The non-dimensional wavenumber at which the increase in slope occurs is about 0[sdot ]02, in good agreement with observations in the ocean reported by Grantet al. (1968).For the two runs for which the stratification was stable, the normalized spectra of the temperature derivative fall on average slightly below the mean of the spectra of the remaining runs in the range in which the slope is approximately one-third. Hence the Reynolds number may not have always been sufficiently high to satisfy completely the conditions for an inertial subrange.Universal inertial-subrange constants were directly evaluated from one-dimensional dissipation spectra and found to be 0[sdot ]54 and 1[sdot ]00 for velocity and temperature, respectively. The constant for velocity is consistent with previously reported values, while the value for temperature differs from some of the previous direct estimates but is only 20% greater than the mean of the indirect estimates. This discrepancy may be explained by the neglect in the indirect estimates of the divergence terms in the conservation equation for the variance of temperature fluctuations. There is weak evidence that the one-dimensional constant, and hence the temperature spectra, may depend upon the turbulence Reynolds number, which varied from 1200 to 4300 in the observations reported.

The problem of determining the one-dimensional kernel of viscoelasticity equation with a source of explosive type

2020 ◽  
Vol 28 (1) ◽  
pp. 43-52
Author(s):  
Durdimurod Kalandarovich Durdiev ◽  
Zhanna Dmitrievna Totieva
Keyword(s):  
Inverse Problem ◽  
Hyperbolic Equations ◽  
Stability Estimate ◽  
Continuous Functions ◽  
System Of Integral Equations ◽  
One Dimensional ◽  
Weighted Norms ◽  
The Stability ◽  
The One ◽  
Global Unique Solvability

AbstractThe integro-differential system of viscoelasticity equations with a source of explosive type is considered. It is assumed that the coefficients of the equations depend only on one spatial variable. The problem of determining the kernel included in the integral terms of the equations is studied. The solution of the problem is reduced to one inverse problem for scalar hyperbolic equations. This inverse problem is replaced by an equivalent system of integral equations for unknown functions. The principle of constricted mapping in the space of continuous functions with weighted norms to the latter is applied. The theorem of global unique solvability is proved and the stability estimate of solution to the inverse problem is obtained.

scholarly journals Tomographic retrievals of ozone with the OMPS Limb Profiler: algorithm description and preliminary results

2018 ◽  
Vol 11 (4) ◽  
pp. 2375-2393 ◽  
Author(s):  
Daniel J. Zawada ◽  
Landon A. Rieger ◽  
Adam E. Bourassa ◽  
Douglas A. Degenstein
Keyword(s):  
Polar Vortex ◽  
Retrieval Algorithm ◽  
Vertical Profiles ◽  
Two Dimensional ◽  
One Dimensional ◽  
Preliminary Results ◽  
Horizontal Gradients ◽  
The Mean ◽  
The One

Abstract. Measurements of limb-scattered sunlight from the Ozone Mapping and Profiler Suite Limb Profiler (OMPS-LP) can be used to obtain vertical profiles of ozone in the stratosphere. In this paper we describe a two-dimensional, or tomographic, retrieval algorithm for OMPS-LP where variations are retrieved simultaneously in altitude and the along-orbital-track dimension. The algorithm has been applied to measurements from the center slit for the full OMPS-LP mission to create the publicly available University of Saskatchewan (USask) OMPS-LP 2D v1.0.2 dataset. Tropical ozone anomalies are compared with measurements from the Microwave Limb Sounder (MLS), where differences are less than 5 % of the mean ozone value for the majority of the stratosphere. Examples of near-coincident measurements with MLS are also shown, and agreement at the 5 % level is observed for the majority of the stratosphere. Both simulated retrievals and coincident comparisons with MLS are shown at the edge of the polar vortex, comparing the results to a traditional one-dimensional retrieval. The one-dimensional retrieval is shown to consistently overestimate the amount of ozone in areas of large horizontal gradients relative to both MLS and the two-dimensional retrieval.

Sign in / Sign up

Export Citation Format

Share Document