渗流问题灰色数值模型的解法研究

灰色数值模型的求解是研究灰色数值模型的一个重要问题 .本文根据灰集合、灰数及其灰色运算规则 ,在渗流系统的基本灰色数值模型的基础上 ,分析了求解这类模型的一整套灰色数值算法 ,并对灰色数值算法、普通算法和经典数值方法的计算结果进行了全面比较 ,论证了灰色数值算法对灰信息传递的正确性和对渗流系统描述的合理性 .     

The Modeling Dynamics of HIV and CD4$$^{}$$ T-cells During Primary Infection in Fractional Order: Numerical Simulation

This article is devoted to introduce a numerical treatment using Adams–Bashforth–Moulton method of the fractional model of HIV-1 infection of CD4\(^{+}\) T-cells. We study the effect of the changing the average number of viral particles N with different sets of initial conditions on the dynamics of the presented model. The fractional derivative is described in Caputo sense. Special attention is given to present the local stability of the proposed model using fractional Routh–Hurwitz stability criterion. Qualitative results show that the model has two equilibria: the disease-free equilibrium and the endemic equilibrium points. We compare our numerical solutions with those numerical solutions using fourth-order Runge–Kutta method (RK4). The obtained numerical results of the proposed model show the simplicity and the efficiency of the proposed method.     

Local discontinuous Galerkin method for solving an N‐carrier system

In this article, we present local discontinuous Galerkin (LDG) method for solving a model of energy exchanges in an N ‐carrier system with Neumann boundary conditions. This model extends the concept of the well‐known parabolic two‐step model for microheat transfer to the energy exchanges in a generalized N ‐carrier system with heat sources. The energy norm stability and error estimate of the LDG method is proved for solving N ‐carrier system. Some numerical examples are given. The numerical results when compared with the exact solution and other numerical results indicate that the present method is seen to be a very good alternative to some existing techniques for realistic problems. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

Application of Tikhonov’s regularization method to solve inverse problems for two population models

Tikhonov’s regularization method is applied to numerical solution of inverse problems for two population models. For the first model we solve the inverse problem that involves simultaneous determination of the mortality rate and the initial distribution of individuals given supplementary information on population density. For the second model we determine the growth rate of the individuals given additional information about their density. Examples of numerical solution are presented for both inverse problems. __________ Translated from Prikladnaya Matematika i Informatika, No. 23, pp. 5–14, 2006.     

Finite element method for viscoelastic fluids flows: approximation of the Maxwell model

We discuss about the finite element approximation of viscoelastic fluid flow. We consider a fluid obeying the Oldroyd model and particularly we study the purely viscoelastic case, the so-called Maxwell model, important in practice for the applications. We examine two kinds of methods used for the approximation of the Maxwell model: method using a splitting technique and finite element method satisfying inf-sup conditions relating tensor and velocity. We present numerical results for these methods and we discuss about their stability. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Oscillation analysis of θ‐methods for the Nicholson's blowflies model

This paper is concerned with the oscillation of numerical solution for the Nicholson's blowflies model. Using two kinds of θ‐methods, namely, the linear θ‐method and the one‐leg θ‐method, several conditions under which the numerical solution oscillates are derived. Moreover, it is shown that every non‐oscillatory numerical solution tends to equilibrium point of the original continuous‐time model. Finally, numerical experiments are provided to illustrate the analytical results. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical predictions of three-dimensional flow in a ventilated room using turbulence models

A three-dimensional recirculation flow in a ventilated room was predicted by the numerical methods in which the turbulence models are applied. The predicted results are compared with the experimental results obtained in a model room in order to estimate the practical utilities of such methods from the viewpoint of engineering. Taking account of the practicability of prediction method which the engineers regard as important, two turbulence models were selected and they were incorporated into the numerical prediction methods respectively. One is the two-equation model, in which transport equations of turbulence energy and its rate of dissipation are adopted. The other is the Deardoff's model, in which the subgrid scale eddy coefficient is utilized. The prediction was made by each numerical method. Consequently, no noticeable difference is recognized between both predicted results. Each result is compared with the experimental results. Generally speaking, each agreement is good with regard to the mean velocity. Thus we can conclude that the numerical method using the two-equation model has more practical utility than that using Deardoff's model, because it can give the solutions in a shorter computer time.     

Characteristic-based schemes for dispersive waves I. The method of characteristics for smooth solutions

In order to embark on the development of numerical schemes for stiff problems, we have studied a model of relaxing heat flow. To isolate those errors unavoidably associated with discretization, a method of characteristics is developed, containing three free parameters depending on the stiffness ratio. It is shown that such “decoupled” schemes do not take into account the interaction between the wave families and hence result in incorrect wave speeds. We also demonstrate that schemes can differ by up to two orders of magnitude in their rms errors even while maintaining second-order accuracy. We show that no method of characteristics solution can be better than second-order accurate. Next, we develop “coupled” schemes which account for the interactions, and here we obtain two additional free parameters. We demonstrate how coupling of the two wave families can be introduced in simple ways and how the results are greatly enhanced by this coupling. Finally, numerical results for several decoupled and coupled schemes are presented, and we observe that dispersion relationships can be a very useful qualitative tool for analysis of numerical algorithms for dispersive waves. © 1993 John Wiley & Sons, Inc.     

Exact controllability of a spherical cap: Numerical implementation of HUM

This paper deals with the numerical implementation of the exact boundary controllability of the Reissner model for shallow spherical shells (Ref. 1). The problem is attacked by the Hilbert uniqueness method (HUM, Refs. 2–4), and we propose a semidiscrete method for the numerical approximation of the minimization problem associated to the exact controllability problem. The numerical results compare well with the results obtained by a finite difference and conjugate gradient method in Ref. 5.This work was done when the first two authors were at CNR-IAC, Rome, Italy as Graduate Students.     

Numerical simulation two phase flows of casting filling process using SOLA particle level set method

Mould filling process is a typical gas–liquid metal two phase flow phenomenon. Numerical simulation of the two phase flows of mould filling process can be used to properly predicate the back pressure effect, the gas entrapment defects, and better understand the complex motions of the gas phase and the liquid phase. In this paper, a novel sharp interface incompressible two phase numerical model for mould filling process is presented. A simple ghost fluid method like discretization method and a density evaluation method at face centers of finite difference staggered grid are proposed to overcome the difficulties when solving two phase Navier–Stokes equations with large-density ratio and large-viscosity ratio. A new mass conservation particle level set method is developed to capture the gas–liquid metal phase interface. The classical pressure-correction based SOLA algorithm is modified to solve the two phase Navier–Stokes equations. Two numerical tests including the Zalesak disk problem and the broken dam problem are used to demonstrate the accuracy of the present method. The numerical method is then adopted to simulate three mould filling examples including two high speed CCD camera imaging water filling experiments and an in situ X-ray imaging experiment of pure aluminum filling. The simulation results are in good agreement with the experiments.     

Hopf Bifurcation Analysis for a Delayed Business Cycle Model-The Equivalence of Multiple Time Scales Versus Center Manifold Reduction Methods

In this paper, we study the Hopf bifurcation of a model with a second order term, which is the business cycle model with delay. Multiple time scales method, which is mainly used by the engineering researchers, and center manifold reduction method, which is mainly used by researchers from mathematical society, are used to derive the two types of normal forms near the Hopf critical point. A comparison between the two methods shows that the two normal forms are equivalent. Scholars can derive the normal form by choosing appropriate methods according to their actual demands. Moreover, bifurcation analysis and numerical simulations are given to verify the analytical predictions.     

A robust multigrid approach for variational image registration models

Variational registration models are non-rigid and deformable imaging techniques for accurate registration of two images. As with other models for inverse problems using the Tikhonov regularization, they must have a suitably chosen regularization term as well as a data fitting term. One distinct feature of registration models is that their fitting term is always highly nonlinear and this nonlinearity restricts the class of numerical methods that are applicable. This paper first reviews the current state-of-the-art numerical methods for such models and observes that the nonlinear fitting term is mostly ‘avoided’ in developing fast multigrid methods. It then proposes a unified approach for designing fixed point type smoothers for multigrid methods. The diffusion registration model (second-order equations) and a curvature model (fourth-order equations) are used to illustrate our robust methodology. Analysis of the proposed smoothers and comparisons to other methods are given. As expected of a multigrid method, being many orders of magnitude faster than the unilevel gradient descent approach, the proposed numerical approach delivers fast and accurate results for a range of synthetic and real test images.     

Transient analysis of two queues in parallel with jockeying

Recently, Tarabia (Appl. Math. Model., 2008, 802) studied the steady-state probabilities of two parallel queues with jockeying and restricted capacities, using the matrix-analytical technique. In this paper, the differential–difference equations which describe the transient state case are derived. Using the fourth order Runge–Kutta method and randomization methods, transient-state probabilities of the Tarabia (2008) model are computed. It is shown that these two methods are closely related, but that the randomization method is superior to the Runge–Kutta method. In the transient case, a numerical comparison between Tarabia's model and Conolly's (J. Appl. Prob., 1984, 394) model is presented to highlight the effect of jockeying on the average of the queue length and the waiting time. Finally, some illustrative numerical results are provided, and conclusions are presented.     

A lattice Boltzmann model for the eddy–stream equations in two-dimensional incompressible flows

A lattice Boltzmann model for two-dimensional incompressible flows with eddy–stream equations is proposed. By using two kinds of distribution functions and employing several higher-order moments of equilibrium distribution functions, the eddy equation and stream function equation with the second-order truncation error are obtained. In the numerical examples, we compared the numerical results of this scheme with those obtained by other classical method. The numerical results agree well with the classical ones.     

Numerical simulation of the endolymph flow in a semicircular canal

We present a concept for the simulation of the fluid flow in the semicircular canals of the inner ear. Based on the temporal dynamics of an idealized model configuration we formulate the governing equations and devise a strategy for their numerical solution. We analyze the proposed method and find a numerical instability. This instability is characterized by an error bound which can serve as a guideline for tuning the numerical method to the specific boundary conditions. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Local preconditioners for two‐level non‐overlapping domain decomposition methods

We consider additive two‐level preconditioners, with a local and a global component, for the Schur complement system arising in non‐overlapping domain decomposition methods. We propose two new parallelizable local preconditioners. The first one is a computationally cheap but numerically relevant alternative to the classical block Jacobi preconditioner. The second one exploits all the information from the local Schur complement matrices and demonstrates an attractive numerical behaviour on heterogeneous and anisotropic problems. We also propose two implementations based on approximate Schur complement matrices that are cheaper alternatives to construct the given preconditioners but that preserve their good numerical behaviour. Through extensive computational experiments we study the numerical scalability and the robustness of the proposed preconditioners and compare their numerical performance with well‐known robust preconditioners such as BPS and the balancing Neumann–Neumann method. Finally, we describe a parallel implementation on distributed memory computers of some of the proposed techniques and report parallel performances. Copyright © 2001 John Wiley & Sons, Ltd.     

Fitted Galerkin spectral method to solve delay partial differential equations

In this paper, we consider a class of parabolic partial differential equations with a time delay. The first model equation is the mixed problems for scalar generalized diffusion equation with a delay, whereas the second model equation is a delayed reaction‐diffusion equation. Both of these models have inherent complex nature because of which their analytical solutions are hardly obtainable, and therefore, one has to seek numerical treatments for their approximate solutions. To this end, we develop a fitted Galerkin spectral method for solving this problem. We derive optimal error estimates based on weak formulations for the fully discrete problems. Some numerical experiments are also provided at the end. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical Modeling of Laser-Induced Cavitation Bubbles with a Finite Volume Method

We report on the numerical modeling of laser-induced cavitation bubbles with the finite volume method using the open source software package OpenFOAM. The numerical model is validated by comparison to experimental data for the two cases of a bubble collapsing under normal ambient conditions in an unbounded liquid as well as close to a solid wall. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Fourth-Order Symmetric DIRK Methods for Periodic Stiff Problems

New symmetric DIRK methods specially adapted to the numerical integration of first-order stiff ODE systems with periodic solutions are obtained. Our interest is focused on the dispersion (phase errors) of the dominant components in the numerical oscillations when these methods are applied to the homogeneous linear test model. Based on this homogeneous test model we derive the dispersion conditions for symmetric DIRK methods as well as symmetric stability functions with real poles and maximal dispersion order. Two new fourth-order symmetric methods with four and five stages are obtained. One of the methods is fourth-order dispersive whereas the other method is symplectic and sixth-order dispersive. These methods have been applied to a number of test problems (linear as well as nonlinear) and some numerical results are presented to show their efficiency when they are compared with the symplectic DIRK method derived by Sanz-Serna and Abia (SIAM J. Numer. Anal. 28 (1991) 1081–1096).     

Hybrid simulation of space plasmas: Models with massless fluid representation of electrons. II. Slow and intermediate shocks

This is a second article in a series of reviews on hybrid simulation of low-frequency processes in space plasmas. A hybrid model is described with ions represented by particles and electrons by a massless fluid. The main numerical schemes for the implementation of this model are described: the generalized Ohm law scheme and the predictor-corrector scheme. The first part of the article provides basic back-ground information: MHD models (ideal, resistive, and Hall model); the Rankine-Hugoniot relationship for MHD discontinuities; the Hoffman-Teller coordinate system; and a classification of discontinuities. The review part of the article surveys the literature on simulation of slow shocks (including switch-off shocks) and intermediate shocks. The survey of literature on hybrid simulation of intermediate shocks is concluded with a review of studies that use two different numerical codes (the hybrid model and the resistive Hall MHD model). The computation results produced by the two codes are compared. The concluding part presents some remarks concerning the existence of intermediate shocks and their relationship with rotational discontinuities in various numerical models (ideal MHD, resistive MHD, the hybrid model). Translated from Prikladnaya Matematika i Informatika, No. 2, pp. 5–33, 1999.