A numerical computation of the non-dimensional form of a non-linear hydrodynamic model in a uniform reservoir

A mathematical model is used to simulate the water current and the elevation in a uniform reservoir. A non-linear hydrodynamic model that provides the velocity field and elevation of the water flow is considered. In the simulating process, the Lax–Wendroff technique is used to approximate the solutions. The numerical solution can be the input data for a water-quality model that is applicable for the optimal control of water treatment in the system to achieve minimum cost.     

A stochastic Keller–Segel model of chemotaxis

We introduce stochastic models of chemotaxis generalizing the deterministic Keller–Segel model. These models include fluctuations which are important in systems with small particle numbers or close to a critical point. Following Dean’s approach, we derive the exact kinetic equation satisfied by the density distribution of cells. In the mean field limit where statistical correlations between cells are neglected, we recover the Keller–Segel model governing the smooth density field. We also consider hydrodynamic and kinetic models of chemotaxis that take into account the inertia of the particles and lead to a delay in the adjustment of the velocity of cells with the chemotactic gradient. We make the connection with the Cattaneo model of chemotaxis and the telegraph equation.     

Two-level Galerkin–Lagrange multipliers method for the stationary Navier–Stokes equations

In this paper, we combine the Galerkin–Lagrange multiplier (GLM) method with the two-level method to solve the stationary Navier–Stokes equations in order to avoid the time-consuming process and the construction of zero-divergence elements. Different quadrilateral partitions are used for approximating the velocity and the pressure. Then some error estimates are obtained and some numerical results of the GLM method and the two-level GLM method are given. The results show that the two-level method based on the GLM method is more efficient than the GLM method under the convergence rate of same order.     

Numerical modelling of subcritical open channel flow using the K-ε turbulence model and the penalty function finite element technique

A numerical model has been developed that employs the penalty function finite element technique to solve the vertically averaged hydrodynamic and turbulence model equations for a water body using isoparametric elements. The full elliptic forms of the equations are solved, thereby allowing recirculating flows to be calculated. Alternative momentum dispersion and turbulence closure models are proposed and evaluated by comparing model predictions with experimental data for strongly curved subcritical open channel flow. The results of these simulations indicate that the depth-averaged two-equation k-ε turbulence model yields excellent agreement with experimental observations. In addition, it appears that neither the streamline curvature modification of the depth-averaged k-ε model, nor the momentum dispersion models based on the assumption of helicoidal flow in a curved channel, yield significant improvement in the present model predictions. Overall model predictions are found to be as good as those of a more complex and restricted three-dimensional model.     

The dynamic complexity of a host–parasitoid model with a Beddington–DeAngelis functional response

This paper investigates the complex dynamics in a discrete-time model of predator–prey interaction with a Beddington–DeAngelis functional response. Local stability analysis of this model is carried out and many forms of complexities are observed using ecology theories and numerical simulation of the global behavior. Furthermore, the existence of a strange attractor and computation of the largest Lyapunov exponent also demonstrate the chaotic dynamic behavior of the model. The results show that the system exhibits rich complexity features such as stable, periodic and chaotic dynamics.     

A 3D hydrodynamic numerical model of the Río de la Plata and Montevideo’s coastal zone

The 3D hydrodynamic numerical model MOHID was applied in the Río de la Plata and Montevideo coastal zone in order to represent the main dynamics and to study its complex circulation pattern. The hydrodynamic model was calibrated and validated considering the following main forces: fresh water flow, astronomical and meteorological tides in the oceanic boundary, and wind acting on the water surface. A series of water levels measured at six coastal stations and vertical profiles of current velocity measured at four different locations in the estuarine zone of the Río de la Plata were used for calibrating and validating the hydrodynamic model. The calibration process was carried out in two steps. First the astronomical waves propagation was calibrated comparing harmonic constants of observed and computed sea surface elevation data. Next, both the astronomical and meteorological wave propagation was calibrated. Direct comparison of scatter plot and root-mean square errors of model results and field data were used when evaluating the calibration quality. The calibrated model shows good agreement with the measured water surface level in the entire domain with mean error values being minor than 20% of the measured data and correlation factors higher than 0.74. Also, the intensity and velocity direction observed in the currents data are well represented by the model in both bottom and surface levels with errors similar to 30% of the currents data components. Using the 3D calibrated model the bottom and surface residual circulation for a four month period of time was analyzed.     

Numerical solution of the nonlinear Klein–Gordon equation using radial basis functions

The nonlinear Klein–Gordon equation is used to model many nonlinear phenomena. In this paper, we propose a numerical scheme to solve the one-dimensional nonlinear Klein–Gordon equation with quadratic and cubic nonlinearity. Our scheme uses the collocation points and approximates the solution using Thin Plate Splines (TPS) radial basis functions (RBF). The implementation of the method is simple as finite difference methods. The results of numerical experiments are presented, and are compared with analytical solutions to confirm the good accuracy of the presented scheme.     

Numerical modelling of residual flow and salinity in the Río de la Plata

The Río de la Plata discharges into the Atlantic Ocean. The particular characteristics of the study area, the variable width and shallowness of the river, the high fluvial discharges and the dynamic processes involving interactions between river discharges, tidal currents and wind, generate complex velocity and salinity fields. We applied the hydrodynamic model RMA-10 to examine the effects of various forcing (tides, flow discharge and winds) on residual currents and salinity fields in the Río de la Plata, focusing on the outer zone of the river. The RMA-10 code, developed by Ian King, is a multiparameter finite element model representing estuarine flow in three dimensions. In this study the model has been applied in a depth-averaged-baroclinic mode and a series of observed data is used for model calibration and verification. The model result shows that it is able to simulate velocity and the salinity fields with a reasonable accuracy. The analysis of residual currents in the river, when forced by freshwater discharge and astronomical tide, shows that the flow discharge takes place mainly over the shallower areas of the river and that the saline water is advected up-river through the deeper channels. The numerical simulations show that the winds from the South-West and North-East quadrants have a great influence over the salinity and velocity fields.     

A sigma coordinate 3D k– model for turbulent free surface flow over a submerged structure

A fully three-dimensional unsteady flow model is developed to simulate free surface flow over a submerged structure. A new sigma coordinate is used to map the physical domain containing the wavy free surface and uneven bottom to a rectangular prism, and to keep the size of the submerged block unchanged in the sigma coordinate system. The numerical difficulty encountered in the conventional sigma coordinate system in which the block changes dynamically due to the time varying free surface is thus eliminated. A split operator scheme is used in the numerical solution so that different numerical schemes can be purposely chosen to deal with the distinctive mathematical and physical characteristics of the phenomena at different steps. k– model is used in the parameterization of turbulence due to its efficiency and reasonable performance. The model is applied to simulate the propagation of a solitary wave with good results. It is subsequently used to simulate a free surface flow against a submerged cube with one face perpendicular (or 45° inclined) to the flow. The numerical results compare favorably with the experimental measurements. In particular, no excessive turbulent kinetic energy is accumulated at the impingement regions.     

Algebraic multigrid preconditioners for the bidomain reaction–diffusion system

The so-called bidomain system is possibly the most complete model for the cardiac bioelectric activity. It consists of a reaction–diffusion system, modeling the intra, extracellular and transmembrane potentials, coupled through a nonlinear reaction term with a stiff system of ordinary differential equations describing the ionic currents through the cellular membrane. In this paper we address the problem of efficiently solving the large linear system arising in the finite element discretization of the bidomain model, when a semiimplicit method in time is employed. We analyze the use of structured algebraic multigrid preconditioners on two major formulations of the model, and report on our numerical experience under different discretization parameters and various discontinuity properties of the conductivity tensors. Our numerical results show that the less exercised formulation provides the best overall performance on a typical simulation of the myocardium excitation process.     

Identification of a free boundary in Norton–Hoff flows with thermal effects

This work deals with a free boundary identification problem in a steady viscoplastic flow. We provide a novel identification model based on a non-linear optimization. The fluid motion is governed by the incompressible Norton–Hoff model coupled with the heat equation. The viscosity of the fluid is modeled by the non-linear Arrhenius law. Our point of view is to treat the problem as a shape sensitivity of a cost functional formulated on the free boundary and governed by the normal component of the velocity of the flow. We analyze the mathematical statement of the forward problem. The equations related to the free boundary are simplified. Various properties of this optimization are proved. Since the state of Norton–Hoff model is not regular enough we introduce a parameter penalization. The shape gradient of the considered cost functional is given in the strong sense up to the parameter of penalization. We supply the expression of the shape gradient in a weak sense.     

On the thermodynamic limit for Hartree–Fock type models

We continue here our study [10–13] of the thermodynamic limit for various models of Quantum Chemistry, this time focusing on the Hartree–Fock type models. For the reduced Hartree–Fock models, we prove the existence of the thermodynamic limit for the energy per unit volume. We also define a periodic problem associated to the Hartree–Fock model, and prove that it is well-posed.     

Decoupled energy‐law preserving numerical schemes for the Cahn–Hilliard–Darcy system

We study two novel decoupled energy‐law preserving and mass‐conservative numerical schemes for solving the Cahn‐Hilliard‐Darcy system which models two‐phase flow in porous medium or in a Hele–Shaw cell. In the first scheme, the velocity in the Cahn–Hilliard equation is treated explicitly so that the Darcy equation is completely decoupled from the Cahn–Hilliard equation. In the second scheme, an intermediate velocity is used in the Cahn–Hilliard equation which allows for the decoupling. We show that the first scheme preserves a discrete energy law with a time‐step constraint, while the second scheme satisfies an energy law without any constraint and is unconditionally stable. Ample numerical experiments are performed to gauge the efficiency and robustness of our scheme. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 936–954, 2016     

A large-eddy-based lattice Boltzmann model for turbulent flow simulation

In this paper a novel and simple large-eddy-based lattice Boltzmann model is proposed to simulate two-dimensional turbulence. Unlike existing lattice Boltzmann models for turbulent flow simulation, which were based on primitive-variables Navier–Stokes equations, the target macroscopic equations of the present model are vorticity-streamfunction equations. Thanks to the intrinsic features of vorticity-streamfunction equations, the present model is efficient, stable and simple for two-dimensional turbulence simulation. The advantages of the present model are validated by numerical experiments.     

A joint SSPDF–CMC method for simulating turbulent methane–air jet combustion

A joint single scalar probability density function and conditional moment closure (SSPDF–CMC) method is proposed for modeling a turbulent methane–air jet flame. In general, the probability density function (PDF) of passive scalar (such as mixture fraction) is non-Gaussian and not fully determined by the advecting velocity field, therefore the presumed shape of PDF of mixture fraction assumed as clipped Gaussian distribution or beta function in normal conditional moment closure (CMC) method is incorrect. In SSPDF–CMC method, the PDF of mixture fraction is obtained using a Monte-Carlo method to solve a PDF transport equation. An assumption that the averaged scalar advection is approximately equal to the averaged scalar dissipation in the wake of a grid-generated turbulence flow is adopted to model the averaged scalar dissipation. The predictions using the proposed method are compared with those using the conventional CMC method and the experimental data. It is seen that the predicted Favre conditional averaged statistics and Favre unconditional averaged statistics using the proposed method are in better agreement with the measurement data than those using the conventional CMC method. The predicted conditional or unconditional mean NO even using the SSPDF model is only in fair agreement with the experiments. It shows that the first-order closure for the conditional reaction rate of NO should be improved.     

Hybrid Simulation of Space Plasmas: Models with Massless Fluid Representation of Electrons. IV. Kelvin–Helmholtz Instability

This is a survey of the literature on hybrid simulation of the Kelvin–Helmholtz instability. We start with a brief review of the theory: the simplest model of the instability—a transition layer in the form of a tangential discontinuity; compressibility of the medium; finite size of the velocity shear region; pressure anisotropy. We then describe the electromagnetic hybrid model (ions as particles and electrons as a massless fluid) and the main numerical schemes. We review the studies on two-dimensional and three-dimensional hybrid simulation of the process of particle mixing across the magnetopause shear layer driven by the onset of a Kelvin–Helmholtz instability. The article concludes with a survey of literature on hybrid simulation of the Kelvin–Helmholtz instability in finite-size objects: jets moving across the magnetic field in the middle of the field reversal layer; interaction between a magnetized plasma flow and a cylindrical plasma source with zero own magnetic field.     

Bivariate Birnbaum–Saunders distribution and associated inference

Univariate Birnbaum–Saunders distribution has been used quite effectively to model positively skewed data, especially lifetime data and crack growth data. In this paper, we introduce bivariate Birnbaum–Saunders distribution which is an absolutely continuous distribution whose marginals are univariate Birnbaum–Saunders distributions. Different properties of this bivariate Birnbaum–Saunders distribution are then discussed. This new family has five unknown parameters and it is shown that the maximum likelihood estimators can be obtained by solving two non-linear equations. We also propose simple modified moment estimators for the unknown parameters which are explicit and can therefore be used effectively as an initial guess for the computation of the maximum likelihood estimators. We then present the asymptotic distributions of the maximum likelihood estimators and use them to construct confidence intervals for the parameters. We also discuss likelihood ratio tests for some hypotheses of interest. Monte Carlo simulations are then carried out to examine the performance of the proposed estimators. Finally, a numerical data analysis is performed in order to illustrate all the methods of inference discussed here.     

Arch model with Box–Cox transformed dependent variable

Box–Cox power transformation has been used traditionally to linearise otherwise nonlinear models. In this paper, Engle's linear ARCH specification is considered for a regression model in which the dependent variable is Box–Cox transformed. The consequent issues arising in both testing and estimation of the model are investigated. A Lagrange multiplier test is also developed to test Engle's linear ARCH model against this wider class of models. The usefulness of this generalisation is examined by applying it to the daily closing prices on the Bombay Stock Exchange Sensitive Index, and the findings strongly favour the proposed model.     

A nonconforming rectangular finite element pair for the Darcy–Stokes–Brinkman model

In this article, we consider the Darcy–Stokes–Brinkman model that can be sorted into three problems: the Darcy problems, the Stokes–Brinkman interface problems and the coupled Darcy–Stokes problems. We study finite element approximation of the model with Dirichlet boundary conditions and make a unified analysis of the three problems based on nonconforming element. Optimal error estimates for the fluid velocity and pressure are derived. Finally, we present some numerical examples verifying the theoretical predictions. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Vibration of a Reissner–Mindlin–Timoshenko plate–beam system

In this paper, we consider a plate–beam system in which the Reissner–Mindlin plate model is combined with the Timoshenko beam model. Natural frequencies and vibration modes for the system are calculated using the finite element method. The interface conditions at the contact between the plate and beams are discussed in some detail. The impact of regularity on the enforcement of certain interface conditions is an important feature of the paper.