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 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.

Growth of vortical disturbances entrained in the entrance region of a circular pipe

2021 ◽  
Vol 932 ◽  
Author(s):  
Pierre Ricco ◽  
Claudia Alvarenga
Keyword(s):  
Reynolds Number ◽  
Initial Conditions ◽  
Maximum Growth ◽  
Streamwise Velocity ◽  
Velocity Profiles ◽  
Boundary Region ◽  
Base Flow ◽  
Entrance Region ◽  
Circular Pipe ◽  
Large Reynolds Number

The development and growth of unsteady three-dimensional vortical disturbances entrained in the entry region of a circular pipe is investigated by asymptotic and numerical methods for Reynolds numbers between $1000$ and $10\,000$ , based on the pipe radius and the bulk velocity. Near the pipe mouth, composite asymptotic solutions describe the dynamics of the oncoming disturbances, revealing how these disturbances are altered by the viscous layer attached to the pipe wall. The perturbation velocity profiles near the pipe mouth are employed as rigorous initial conditions for the boundary-region equations, which describe the flow in the limit of low frequency and large Reynolds number. The disturbance flow is initially primarily present within the base-flow boundary layer in the form of streamwise-elongated vortical structures, i.e. the streamwise velocity component displays an intense algebraic growth, while the cross-flow velocity components decay. Farther downstream the disturbance flow occupies the whole pipe, although the base flow is mostly inviscid in the core. The transient growth and subsequent viscous decay are confined in the entrance region, i.e. where the base flow has not reached the fully developed Poiseuille profile. Increasing the Reynolds number and decreasing the frequency causes more intense perturbations, whereas small azimuthal wavelengths and radial characteristic length scales intensify the viscous dissipation of the disturbance. The azimuthal wavelength that causes the maximum growth is found. The velocity profiles are compared successfully with available experimental data and the theoretical results are helpful to interpret the only direct numerical dataset of a disturbed pipe-entry flow.

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 The development of biofilm architecture

2016 ◽  
Vol 472 (2188) ◽  
pp. 20150798 ◽  
Author(s):  
A. C. Fowler ◽  
T. M. Kyrke-Smith ◽  
H. F. Winstanley
Keyword(s):  
Bacterial Biofilm ◽  
Governing Equations ◽  
Biofilm Growth ◽  
One Dimensional ◽  
Direct Numerical Solution ◽  
The Common ◽  
Common Observation ◽  
The Stability ◽  
The One ◽  
A Free Boundary Problem

We extend the one-dimensional polymer solution theory of bacterial biofilm growth described by Winstanley et al . (2011 Proc. R. Soc. A 467 , 1449–1467 ( doi:10.1098/rspa.2010.0327 )) to deal with the problem of the growth of a patch of biofilm in more than one lateral dimension. The extension is non-trivial, as it requires consideration of the rheology of the polymer phase. We use a novel asymptotic technique to reduce the model to a free-boundary problem governed by the equations of Stokes flow with non-standard boundary conditions. We then consider the stability of laterally uniform biofilm growth, and show that the model predicts spatial instability; this is confirmed by a direct numerical solution of the governing equations. The instability results in cusp formation at the biofilm surface and provides an explanation for the common observation of patterned biofilm architectures.

Stable Stationary Solitons of the One-Dimensional Modified Complex Ginzburg-Landau Equation

2007 ◽  
Vol 62 (7-8) ◽  
pp. 368-372
Author(s):  
Woo-Pyo Hong
Keyword(s):  
Landau Equation ◽  
Model Parameters ◽  
Ginzburg Landau ◽  
One Dimensional ◽  
Paraxial Ray Approximation ◽  
New Family ◽  
Ginzburg Landau Equation ◽  
The Stability ◽  
The One ◽  
Paraxial Ray

We report on the existence of a new family of stable stationary solitons of the one-dimensional modified complex Ginzburg-Landau equation. By applying the paraxial ray approximation, we obtain the relation between the width and the peak amplitude of the stationary soliton in terms of the model parameters. We verify the analytical results by direct numerical simulations and show the stability of the stationary solitons.

scholarly journals The acoustic impedance of a laminar viscous jet through a thin circular aperture

2019 ◽  
Vol 864 ◽  
pp. 5-44 ◽  
Author(s):  
David Fabre ◽  
Raffaele Longobardi ◽  
Paul Bonnefis ◽  
Paolo Luchini
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Axisymmetric Flow ◽  
Classical Problem ◽  
Base Flow ◽  
Large Reynolds Number ◽  
Linear Perturbation ◽  
Circular Aperture ◽  
Vena Contracta ◽  
Numerical Resolution

The unsteady axisymmetric flow through a circular aperture in a thin plate subjected to harmonic forcing (for instance under the effect of an incident acoustic wave) is a classical problem first considered by Howe (Proc. R. Soc. Lond. A, vol. 366, 1979, pp. 205–223), using an inviscid model. The purpose of this work is to reconsider this problem through a numerical resolution of the incompressible linearized Navier–Stokes equations (LNSE) in the laminar regime, corresponding to $Re=[500,5000]$. We first compute a steady base flow which allows us to describe the vena contracta phenomenon in agreement with experiments. We then solve a linear problem allowing us to characterize both the spatial amplification of the perturbations and the impedance (or equivalently the Rayleigh conductivity), which is a key quantity to investigate the response of the jet to acoustic forcing. Since the linear perturbation is characterized by a strong spatial amplification, the numerical resolution requires the use of a complex mapping of the axial coordinate in order to enlarge the range of Reynolds number investigated. The results show that the impedances computed with $Re\gtrsim 1500$ collapse onto a single curve, indicating that a large Reynolds number asymptotic regime is effectively reached. However, expressing the results in terms of conductivity leads to substantial deviation with respect to Howe’s model. Finally, we investigate the case of finite-amplitude perturbations through direct numerical simulations (DNS). We show that the impedance predicted by the linear approach remains valid for amplitudes up to order $10^{-1}$, despite the fact that the spatial evolution of the perturbations in the jet is strongly nonlinear.

scholarly journals Chaos and complexity in a simple model of production dynamics

2000 ◽  
Vol 5 (3) ◽  
pp. 179-187 ◽  
Author(s):  
I. Katzorke ◽  
A. Pikovsky
Keyword(s):  
Simple Model ◽  
Dynamical Behavior ◽  
Production Costs ◽  
Order Flow ◽  
One Dimensional ◽  
Dynamical Behaviors ◽  
The Stability ◽  
The One ◽  
Production Dynamics ◽  
Complex Dynamical Behavior

We consider complex dynamical behavior in a simple model of production dynamics, based on the Wiendahl’s funnel approach. In the case of continuous order flow a model of three parallel funnels reduces to the one-dimensional Bernoulli-type map, and demonstrates strong chaotic properties. The optimization of production costs is possible with the OGY method of chaos control. The dynamics changes drastically in the case of discrete order flow. We discuss different dynamical behaviors, the complexity and the stability of this discrete system.

Sign in / Sign up

Export Citation Format

Share Document