Simultaneous perturbation methods for adaptive labor staffing in service systems

We consider the problem of optimizing the workforce of a service system. Adapting the staffing levels in such systems is non-trivial due to large variations in workload and the large number of system parameters do not allow for a brute force search. Further, because these parameters change on a weekly basis, the optimization should not take longer than a few hours. Our aim is to find the optimum staffing levels from a discrete high-dimensional parameter set, that minimizes the long run average of the single-stage cost function, while adhering to the constraints relating to queue stability and service-level agreement (SLA) compliance. The single-stage cost function balances the conflicting objectives of utilizing workers better and attaining the target SLAs. We formulate this problem as a constrained parameterized Markov cost process parameterized by the (discrete) staffing levels. We propose novel simultaneous perturbation stochastic approximation (SPSA)-based algorithms for solving the above problem. The algorithms include both first-order as well as second-order methods and incorporate SPSA-based gradient/Hessian estimates for primal descent, while performing dual ascent for the Lagrange multipliers. Both algorithms are online and update the staffing levels in an incremental fashion. Further, they involve a certain generalized smooth projection operator, which is essential to project the continuous-valued worker parameter tuned by our algorithms onto the discrete set. The smoothness is necessary to ensure that the underlying transition dynamics of the constrained Markov cost process is itself smooth (as a function of the continuous-valued parameter): a critical requirement to prove the convergence of both algorithms. We validate our algorithms via performance simulations based on data from five real-life service systems. For the sake of comparison, we also implement a scatter search based algorithm using state-of-the-art optimization tool-kit OptQuest. From the experiments, we observe that both our algorithms converge empirically and consistently outperform OptQuest in most of the settings considered. This finding coupled with the computational advantage of our algorithms make them amenable for adaptive labor staffing in real-life service systems.

Conformal perturbation theory and higher spin entanglement entropy on the torus

We study the free fermion theory in 1+1 dimensions deformed by chemical potentials for holomorphic, conserved currents at finite temperature and on a spatial circle. For a spin-three chemical potential mu, the deformation is related at high temperatures to a higher spin black hole in hs0] theory on AdS(3) spacetime. We calculate the order mu(2) corrections to the single interval Renyi and entanglement entropies on the torus using the bosonized formulation. A consistent result, satisfying all checks, emerges upon carefully accounting for both perturbative and winding mode contributions in the bosonized language. The order mu(2) corrections involve integrals that are finite but potentially sensitive to contact term singularities. We propose and apply a prescription for defining such integrals which matches the Hamiltonian picture and passes several non-trivial checks for both thermal corrections and the Renyi entropies at this order. The thermal corrections are given by a weight six quasi-modular form, whilst the Renyi entropies are controlled by quasi-elliptic functions of the interval length with modular weight six. We also point out the well known connection between the perturbative expansion of the partition function in powers of the spin-three chemical potential and the Gross-Taylor genus expansion of large-N Yang-Mills theory on the torus. We note the absence of winding mode contributions in this connection, which suggests qualitatively different entanglement entropies for the two systems.

Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man(9)GlcNAc(2)). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct beta sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

Ultra-sensitive pressure dependence of bandgap of rutile-GeO2 revealed by many body perturbation theory

The reported values of bandgap of rutile GeO2 calculated by the standard density functional theory within local-density approximation (LDA)/generalized gradient approximation (GGA) show a wide variation (similar to 2 eV), whose origin remains unresolved. Here, we investigate the reasons for this variation by studying the electronic structure of rutile-GeO2 using many-body perturbation theory within the GW framework. The bandgap as well as valence bandwidth at Gamma-point of rutile phase shows a strong dependence on volume change, which is independent of bandgap underestimation problem of LDA/GGA. This strong dependence originates from a change in hybridization among O-p and Ge-(s and p) orbitals. Furthermore, the parabolic nature of first conduction band along X-Gamma-M direction changes towards a linear dispersion with volume expansion. (C) 2015 AIP Publishing LLC.

Perturbation of nucleo-cytoplasmic transport affects size of nucleus and nucleolus in human cells

Size regulation of human cell nucleus and nucleolus are poorly understood subjects. 3D reconstruction of live image shows that the karyoplasmic ratio (KR) increases by 30-80% in transformed cell lines compared to their immortalized counterpart. The attenuation of nucleo-cytoplasmic transport causes the KR value to increase by 30-50% in immortalized cell lines. Nucleolus volumes are significantly increased in transformed cell lines and the attenuation of nucleo-cytoplasmic transport causes a significant increase in the nucleolus volume of immortalized cell lines. A cytosol and nuclear fraction swapping experiment emphasizes the potential role of unknown cytosolic factors in nuclear and nucleolar size regulation.

Perturbation of the stability of amyloid fibrils through alteration of electrostatic interactions.

The self-assembly of proteins and peptides into polymeric amyloid fibrils is a process that has important implications ranging from the understanding of protein misfolding disorders to the discovery of novel nanobiomaterials. In this study, we probe the stability of fibrils prepared at pH 2.0 and composed of the protein insulin by manipulating electrostatic interactions within the fibril architecture. We demonstrate that strong electrostatic repulsion is sufficient to disrupt the hydrogen-bonded, cross-β network that links insulin molecules and ultimately results in fibril dissociation. The extent of this dissociation correlates well with predictions for colloidal models considering the net global charge of the polypeptide chain, although the kinetics of the process is regulated by the charge state of a single amino acid. We found the fibrils to be maximally stable under their formation conditions. Partial disruption of the cross-β network under conditions where the fibrils remain intact leads to a reduction in their stability. Together, these results support the contention that a major determinant of amyloid stability stems from the interactions in the structured core, and show how the control of electrostatic interactions can be used to characterize the factors that modulate fibril stability.

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

Crack propagation in the power-law nonlinear viscoelastic material

An analysis on crack creep propagation problem of power-law nonlinear viscoelastic materials is presented. The creep incompressilility assumption is used. To simulate fracture behavior of craze region, it is assumed that in the fracture process zone near the crack tip, the cohesive stress sigma(f) acts upon the crack surfaces and resists crack opening. Through a perturbation method, i. e., by superposing the Mode-I applied force onto a referential uniform stress state, which has a trivial solution and gives no effect on the solution of the original problem, the nonlinear viscoelastic problem is reduced to linear problem. For weak nonlinear materials, for which the power-law index n similar or equal to 1, the expressions of stress and crack surface displacement are derived. Then, the fracture process zone local energy criterion is proposed and based on which the formulas of cracking incubation time t

Rheological effect on thermocapillary flow of a liquid film jet painted on a moving boundary

In the present paper, a liquid (or melt) film of relatively high temperature ejected from a vessel and painted on the-moving solid film is analyzed by using the second-order fluid model of the non-Newtonian fluid. The thermocapillary flow driven by the temperature gradient on the free surface of a Newtonian liquid film was discussed before. The effect of rheological fluid on thermocapillary flow is considered in the present paper. The analysis is based on the approximations of lubrication theory and perturbation theory. The equation of liquid height and the process of thermal hydrodynamics of the non-Newtonian liquid film are obtained, and the case of weak effect of the rheological fluid is solved in detail.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220perturbation breaking the equatorial symmetry. The mechanism of development of a saddle point in the meridional plane at higher Re number and its role in the formation of two-vortex flow are analyzed.

Free surface oscillation of thermocapillary convection in liquid bridge of half-floating-zone

Free surface deformations of thermocapillary convection in a small liquid bridge of half floating-zone are studied in the present paper. The relative displacement and phase difference of free surface oscillation are experimentally studied, and the features of free surface oscillation for various applied temperature differences are obtained. It is discovered that there is a sort of surface waves having the character of small perturbation, and having a wave mode of unusually large amplitude in one corner region of the liquid bridge.

Singularity criteria for perturbation series

This paper is concerned with some mathematical aspects of the Van Dyke method inperturbation theory, i.e. the singularity criteria of perturbation series. The author suggestsa sign criterion and a Domb-syke plot for the cases with complex conjugate singularities, thussucceeding in extending the conclusions of Van Dyke's. Subsequently. effects of singularitiesof the lower order upon the criteria are taken into account. In addition, a method of locat-ing singular points is developed by analysing the new perturbation series derived by the Eulertransformation.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.