Finite-volume modelling of system with compressible flow propelled and actuated body motion

A finite-volume model for the gas flow in a compressible flow driven system with two moving components is developed and applied. In this system, high pressure and temperature compressible flow propels one body along a main flow channel. When this body passes an opening in the main flow channel, some of the high pressure gas bleeds off, entering a secondary, telescoping, flow channel that empties into a plenum. The bleed-off flow causes the pressure in the plenum to increase, actuating the motion of a secondary (actuated) body that forms a bounding surface of the plenum, causing the plenum to expand and the secondary channel to lengthen.     

Numerical analysis of the weighted particle method applied to the semiconductor Boltzmann equation

Summary The semiconductor Boltzmann equation involves an integral operator, the kernel of which is a measure supported by a surface. This feature introduces some singularities of the exact solution, which makes the numerical approximation of this equation difficult. This paper is devoted to the error analysis of the weighted particle method (introduced by Mas-Gallic and Raviart [14]) applied to the space homogeneous semiconductor Boltzmann equation. The results are commented in view of the practical use of the method. This paper is closely related to [12], where results of numerical simulations on both test and real problems are given.     

Boundary layers for a fluid–particle interaction system with density-dependent viscosity and cylindrical symmetry

In this paper, we study the initial boundary value problem for a cylindrical symmetry fluid–particle interaction system in three dimensions. The boundary layer phenomena is investigated when the shear viscosity $\mu =\kappa {\rho }^{\beta }$ goes to zero. Furthermore, we establish the boundary layer thickness of the order $O\left({\kappa }^{\alpha }\right)$ for more general initial data when $0<\alpha <\frac{1}{2}$ and give the optimal boundary-layer thickness for the system with more general initial data. As a byproduct, this work improves the corresponding results in Yao et al. (2011) for isentropic compressible Navier–Stokes equations where $0<\alpha <\frac{1}{4}$.     

A competitive and cooperative co-evolutionary approach to multi-objective particle swarm optimization algorithm design

Multi-objective particle swarm optimization (MOPSO) is an optimization technique inspired by bird flocking, which has been steadily gaining attention from the research community because of its high convergence speed. On the other hand, in the face of increasing complexity and dimensionality of today’s application coupled with its tendency of premature convergence due to the high convergence speeds, there is a need to improve the efficiency and effectiveness of MOPSO. In this paper a competitive and cooperative co-evolutionary approach is adapted for multi-objective particle swarm optimization algorithm design, which appears to have considerable potential for solving complex optimization problems by explicitly modeling the co-evolution of competing and cooperating species. The competitive and cooperative co-evolution model helps to produce the reasonable problem decompositions by exploiting any correlation, interdependency between components of the problem. The proposed competitive and cooperative co-evolutionary multi-objective particle swarm optimization algorithm (CCPSO) is validated through comparisons with existing state-of-the-art multi-objective algorithms using established benchmarks and metrics. Simulation results demonstrated that CCPSO shows competitive, if not better, performance as compared to the other algorithms.     

Application of quorum response and information entropy to animal collective motion modeling

“Quorum response” is a type of social interaction in which an individual's chance of choosing an option is a nonlinear function of the number of other individuals already committing to it. This interaction has been widely used to characterize collective decision‐making in animal groups. Here, we first implement it in 1D and 2D models of collective animal movement, and find that the resulting group motion shows the characteristic behaviors which were observed in previous experimental and modeling studies. Further, the analytic form of quorum response renders us an opportunity to propose a mean field theory in 1D with globally interacting particles, so we can estimate the average time period between changes in the group direction (mean switching time). We find that the theoretical results provide an upper bound to the simulation results when the interaction radius grows from local to global. Information entropy, a concept widely used to quantify the uncertainty of a random variable, is introduced here as a new order parameter to study the evolution of systems of two cases in 2D models. The explicitly formulated probability of a particle's dynamic state in the framework of quorum response makes information entropy directly computable. We find that, besides the global order, information entropy can also capture the structural features of local order of the system which previous order parameters such as alignment cannot. © 2016 Wiley Periodicals, Inc. Complexity 21: 584–592, 2016     

Existence of local strong solutions to fluid–beam and fluid–rod interaction systems

We study an unsteady nonlinear fluid–structure interaction problem. We consider a Newtonian incompressible two-dimensional flow described by the Navier–Stokes equations set in an unknown domain depending on the displacement of a structure, which itself satisfies a linear wave equation or a linear beam equation. The fluid and the structure systems are coupled via interface conditions prescribing the continuity of the velocities at the fluid–structure interface and the action-reaction principle. Considering three different structure models, we prove existence of a unique local-in-time strong solution, for which there is no gap between the regularity of the initial data and the regularity of the solution enabling to obtain a blow up alternative. In the case of a damped beam this is an alternative proof (and a generalization to non zero initial displacement) of the result that can be found in [20]. In the case of the wave equation or a beam equation with inertia of rotation, this is, to our knowledge the first result of existence of strong solutions for which no viscosity is added. The key points consist in studying the coupled system without decoupling the fluid from the structure and to use the fluid dissipation to control, in appropriate function spaces, the structure velocity.     

信息能力国际比较的主成分分析法

本文利用主成分分析方法对世界主要国家的信息能力进行了排序比较 ,并对所得结果进行了简要的分析。     

A new formulation of regularized meshless method applied to interior and exterior anisotropic potential problems

This paper proposes a new formulation of regularized meshless method (RMM), which differs from the traditional RMM in that the traditional formulation generates the diagonal elements of influence matrix via null-field integral equations, while our new one directly employs the boundary integral equations at the domain point to evaluate the diagonal elements. We test the present RMM formulation to two-dimensional anisotropic potential problems in finite and infinite domains in comparison with the traditional RMM. Numerical results show that the present RMM sharply outperforms the traditional RMM in the solution of interior problems, while the latter is clearly superior for exterior problems. A rigorous theoretical analysis of circular domain case also corroborates such numerical experiment observations and is provided in the appendix of this paper.     

A post-optimization method for the routing and wavelength assignment problem applied to scheduled lightpath demands

We consider here a NP-hard problem related to the Routing and Wavelength Assignment (RWA) problem in optical networks, dealing with Scheduled Lightpath Demands (SLDs). An SLD is a connection demand between two nodes of the network, during a certain time. Given a set of SLDs, we want to assign a lightpath, i.e. a routing path and a wavelength, to each SLD, so that the total number of required wavelengths is minimized. The constraints are the following: a same wavelength must be assigned all along the edges of the routing path of any SLD; at any time, a given wavelength on a given edge of the network cannot be used to satisfy more than one SLD. To solve this problem, we design a post-optimization method improving the solutions provided by a heuristic. The experimental results show that this post-optimization method is quite efficient to reduce the number of necessary wavelengths.     

An augmented galerkin procedure for the boundary integral method applied to two-dimensional screen and crack problems

Here we present a new solution procedure for Helm-holtz and Laplacian Dirichlet screen and crack problems in IR² via boundary integral equations of the first kind having as an unknown the jump of the normal derivative across the screen or a crack curve T. Under the assumption of local finite energy we show the equivalence of the integral equations and the original boundary value problem. Via the method of local Mellin transform in [5]-[lo] and the calculus of pseudodifferential operators we derive existence, uniqueness and regularity results for the solution of our boundary integral equations together with its explicit behaviour near the screen or crack tips.With our integral equations we set up a Galerkin scheme on T and obtain high quasi-optimal convergence rates by using special singular elements besides regular splines as test and trial functions.     

Local and pointwise error estimates of the local discontinuous Galerkin method applied to the Stokes problem

We prove local and pointwise error estimates for the local discontinuous Galerkin method applied to the Stokes problem in two and three dimensions. By using techniques originally developed by A. Schatz [Math. Comp., 67 (1998), 877-899] to prove pointwise estimates for the Laplace equation, we prove optimal weighted pointwise estimates for both the velocity and the pressure for domains with smooth boundaries.

     

The decomposition method applied to systems of partial differential equations and to the reaction–diffusion Brusselator model

In this work, systems of linear and nonlinear partial differential equations and the reaction–diffusion Brusselator model are handled by applying the decomposition method. The advantage of this work is twofold. Firstly, the decomposition method reduces the computational work. Secondly, in comparison with existing techniques, the decomposition method is an improvement with regard to its accuracy and rapid convergence. The decomposition method has the advantage of being more concise for analytical and numerical purposes.     

Monotone iterative method of upper and lower solutions applied to a multilayer combustion model in porous media

This paper presents a new nonlinear reaction–diffusion–convection system coupled with a system of ordinary differential equations that models a combustion front in a multilayer porous medium. The model includes heat transfer between the layers and heat loss to the external environment. A few assumptions are made to simplify the model, such as incompressibility; then, the unknowns are determined to be the temperature and fuel concentration in each layer. When the fuel concentration in each layer is a known function, we prove the existence and uniqueness of a classical solution for the initial and boundary value problem for the corresponding system. The proof uses a new approach for combustion problems in porous media. We construct monotone iterations of upper and lower solutions and prove that these iterations converge to a unique solution for the problem, first locally and then, in time, globally.     

Using abell's narrative method to build a theory: The case of the theory of distributive justice and its generalization to the theory of comparison processes

This paper examines the usefulness for theoretical work of the narrative method proposed by Peter Abell. Our assessment proceeds by using the narrative method to perform the two main tasks of theoretical analysis—constructing postulates and deriving predictions. Tb illustrate, we focus on the theory of distributive justice and the more general theory of comparison processes to which it led. The results of our assessment of the usefulness of Abell's narrative method for theoretical work indicate that the narrative method has far wider applicability than Abell has claimed for it. For example, (i) it is useful for all theoretical work in the sociobehavioral sciences, not only for theoretical work based on game theory, (ii) it is useful for analyzing thought‐experiments as well as narrative accounts of actual actions and events, and (iii) the events in the narrative need not be restricted to human actions but can include events not traceable to human agency. We conclude also that Abell's narrative method complements the use of mathematical analysis in theoretical work and that it may be especially valuable for theoretical derivation involving two or more theories jointly.     

Network simulation method applied to radiation and viscous dissipation effects on MHD unsteady free convection over vertical porous plate

The effects of thermal radiation and viscous dissipation on magneto-hydrodynamic (MHD) unsteady free-convection flow over a semi-infinite vertical porous plate are analysed. The fluid considered is non-gray (absorption coefficient dependent on wave length). The Network Simulation Method is used to solve the boundary-layer equations based on the finite-difference formulation; only discretization of the spatial co-ordinates is necessary, while time remains as a real continuous variable. This method provides a solution for both transient and steady-state problems at the same time, and programming does not require manipulation of the sophisticated mathematical software that is inherent in other numerical methods. The velocity, temperature, local skin-friction and local Nusselt number are studied for different parameters, including the radiation parameter, Eckert number, magnetic number and suction (or injection).     

Mathematical modeling of the dynamics of the bladder cancer and the immune response applied to a patient: Evolution and short‐term prediction

Bladder cancer is one of the most common malignant diseases in the urinary system and a highly aggressive neoplasm. The prognosis is not favorable usually, and its evolution for particular patients is very difficult to find out. In this paper, we propose a dynamic mathematical model that describes the bladder tumor growth and the immune response evolution. This model is customized for a single patient, determining appropriate model parameter values via model calibration. Due to the uncertainty of the tumor evolution, using the calibrated model parameters, we predict the tumor size and the immune response evolution over the next few months assuming three different scenarios: favorable, neutral, and unfavorable. In the former, it is not expected any trace of the cancer in the middle of September 2018 (after 16 mo). In the neutral scenario, at the same date, a 7‐ to 8‐mm tumor is expected. In the worst case, a 40‐mm tumor is expected. The patient was cited on 10 September 2018 to check the tumor size, and according to the doctors, there was no sign of recurrence. It seems that we are in the favorable scenario. The patient will be called again for follow‐up in mid‐2019.     

The discrete (G′/G)‐expansion method applied to the differential‐difference Burgers equation and the relativistic Toda lattice system

We introduce the discrete (G′/G)‐expansion method for solving nonlinear differential–difference equations (NDDEs). As illustrative examples, we consider the differential–difference Burgers equation and the relativistic Toda lattice system. Discrete solitary, periodic, and rational solutions are obtained in a concise manner. The method is also applicable to other types of NDDEs. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 28: 127‐137, 2012     

Instability of the finite‐difference split‐step method applied to the generalized nonlinear Schrödinger equation. III. external potential and oscillating pulse solutions

This is the final part of a series of articles where we have studied numerical instability (NI) of localized solutions of the generalized nonlinear Schrödinger equation (gNLS). It extends our earlier studies of this topic in two ways. First, it examines differences in the development of the NI between the case of the purely cubic NLS and the case where the gNLS has an external bounded potential. Second, it investigates how the NI is affected by the oscillatory dynamics of the simulated pulse. The latter situation is common when the initial condition is not an exact stationary soliton. We have found that in this case, the NI may remain weak when the time step exceeds the threshold quite significantly. This means that the corresponding numerical solution, while formally numerically unstable, can remain sufficiently accurate over long times, because the numerical noise will stay small. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 633–650, 2017     

The P1 and P2 approximations of the spherical harmonic method as applied to radiative heat transfer computations with allowance for strong blowing from the surface of a vehicle moving in the Jovian atmosphere

The spherical harmonics method in the P₁ and P₂ approximations is used to analyze radiative heat transfer for a space vehicle entering a planet’s atmosphere. Strong blowing of ablated materials from the vehicle surface is taken into account by using a two-layer flow model without allowance for the viscosity and thermal conduction of the gas. The heating and ablation of a multilayered thermal protection system are determined simultaneously with the vehicle’s flight trajectory, which is calculated taking into account the mass loss due to the ablation. The approach is illustrated by computing the flight of a space vehicle shaped as a spherical segment or a spherically blunted cone entering the Jovian atmosphere at a speed of 60 km/s and an entry angle of ?5°.