Pressure boundary condition for the time‐dependent incompressible Navier–Stokes equations

In Gresho and Sani (Int. J. Numer. Methods Fluids 1987; 7 :1111–1145; Incompressible Flow and the Finite Element Method, vol. 2. Wiley: New York, 2000) was proposed an important hypothesis regarding the pressure Poisson equation (PPE) for incompressible flow: Stated there but not proven was a so‐called equivalence theorem (assertion) that stated/asserted that if the Navier–Stokes momentum equation is solved simultaneously with the PPE whose boundary condition (BC) is the Neumann condition obtained by applying the normal component of the momentum equation on the boundary on which the normal component of velocity is specified as a Dirichlet BC, the solution (u, p) would be exactly the same as if the ‘primitive’ equations, in which the PPE plus Neumann BC is replaced by the usual divergence‐free constraint (? · u = 0), were solved instead. This issue is explored in sufficient detail in this paper so as to actually prove the theorem for at least some situations. Additionally, like the original/primitive equations that require no BC for the pressure, the new results establish the same thing when the PPE approach is employed. Copyright © 2005 John Wiley & Sons, Ltd.     

Fluid forces on a square cylinder in oscillating flows with non‐zero‐mean velocities

The unsteady forces on a square cylinder in sinusoidally oscillating flows with non‐zero‐mean velocities are investigated numerically by using a weakly compressible‐flow method with three‐dimensional large eddy simulations. The major parameters in the analysis are Keulegan–Carpenter number (KC) and the ratio between the amplitude and the mean velocities of the approaching flow (A_R). By varying the values of KC and A_R the resulting drag and lift of the cylinders are analyzed systematically at two selected approaching‐flow attack angles (0 and 22.5^°). In the case of the non‐zero attack angle, results show that both the drag and lift histories can be adequately described by Morison equations. However, Morison equations fail to correctly describing the lift history as the attack angle is zero. In addition, when the ratio of A_R/KC is near the Strouhal number of the bluff‐body flow, the resulting drag is promoted due to the occurrence of resonance. Based on the results of systematic analyses, finally, the mean and inertia force coefficients at the two selected attack angles are presented as functions of KC and A_R based on the Morison relationships. Copyright © 2008 John Wiley & Sons, Ltd.     

The effects of control‐volume distributed multi‐point flux approximation (CVD‐MPFA) on upscaling—A study using the CVD‐MPFA schemes

A family of flux‐continuous, locally conservative, finite‐volume schemes has been developed for solving the general geometry‐permeability tensor (petroleum reservoir‐simulation) pressure equation on structured and unstructured grids and are control‐volume distributed (textit Comput. Geo. 1998; 2 :259–290; Comput. Geo. 2002; 6 :433–452). The schemes are applicable to diagonal and full tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir‐simulation schemes (two‐point flux approximation) when applied to full tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization (Int. J. Numer. Meth. Fluids 2006; 51 :1177–1203). Improved convergence (for two‐ and three‐dimensional formulation) using the quadrature parameterization has been observed for the family of flux‐continuous control‐volume distributed multi‐point flux approximation (CVD‐MPFA) schemes (Ph.D. Thesis, University of Wales, Swansea, U.K., 2007). In this paper family of flux‐continuous (CVD‐MPFA) schemes are used as a part of numerical upscaling procedure for upscaling the fine‐scale grid information (permeability) onto a coarse grid scale. A series of data‐sets (SPE, 2001) are tested where the upscaled permeability tensor is computed on a sequence of grid levels using the same fixed range of quadrature points in each case. The refinement studies presented involve:

• (i) Refinement comparison study: In this study, permeability distribution for cells at each grid level is obtained by upscaling directly from the fine‐scale permeability field as in standard simulation practice.
• (ii) Refinement study with renormalized permeability: In this refinement comparison, the local permeability is upscaled to the next grid level hierarchically, so that permeability values are renormalized to each coarser level. Hence, showing only the effect of increased grid resolution on upscaled permeability, compared with that obtained directly from the fine‐scale solution.
• (iii) Refinement study with invariant permeability distribution: In this study, a classical mathematical convergence test is performed. The same coarse‐scale underlying permeability map is preserved on all grid levels including the fine‐scale reference solution.

The study is carried out for the discretization of the scheme in physical space. The benefit of using specific quadrature points is demonstrated for upscaling in this study and superconvergence is observed. Copyright © 2010 John Wiley & Sons, Ltd.     

Positive solutions to （n-1,1） m-point boundary value problems with dependence on the first order derivative

Using the extension of Krasnoselskii＇s fixed point theorem in a cone, we prove the existence of at least one positive solution to the nonlinear nth order m-point boundary value problem with dependence on the first order derivative. The associated Green＇s function for the nth order m-point boundary value problem is given, and growth conditions are imposed on the nonlinear term f which ensures the existence of at least one positive solution. A simple example is presented to illustrate applications of the obtained results.     

A mass‐fraction‐based interface‐capturing method for multi‐component flow

An interface‐capturing method based on mass fraction is developed to solve the Riemann problem in multi‐component compressible flow. Equations of mass fraction with modified form, which is derived from conservative equations of mass, are employed here to capture the interface. By introducing mass fraction into Euler equations system, as well as other conservative coefficients, a quasi‐conservative numerical model is created. Numerical examples show that the mass fraction model performs well not only in multi‐component fluids modeled by simple stiffened gas equation of state (EOS) but also in that modeled by complex Mie–Grüneisen EOS. Moreover, the mass fraction model is applied to Riemann problem with piecewise EOS; the expression of which depends on density. It is found that the mass fraction model can well adapt to the analytic change in piecewise EOS and produce accuracy solutions with fewer unknown quantities, and the model can be easily extended to m‐component fluid mixture by using only m + 4 equations with no additional conditions. Copyright © 2013 John Wiley & Sons, Ltd.     

A comparison of the GWCE and mixed P − P1 formulations in finite‐element linearized shallow‐water models

The appearance of spurious pressure modes in early shallow‐water (SW) models has resulted in two common strategies in the finite element (FE) community: using mixed primitive variable and generalized wave continuity equation (GWCE) formulations of the SW equations. One FE scheme in particular, the P ? P₁ pair, combined with the primitive equations may be advantageously compared with the wave equation formulations and both schemes have similar data structures. Our focus here is on comparing these two approaches for a number of measures including stability, accuracy, efficiency, conservation properties, and consistency. The main part of the analysis centres on stability and accuracy results via Fourier‐based dispersion analyses in the context of the linear SW equations. The numerical solutions of test problems are found to be in good agreement with the analytical results. Copyright © 2011 John Wiley & Sons, Ltd.     

Least‐squares finite element processes in h,p, k mathematical and computational framework for a non‐linear conservation law

This paper considers numerical simulation of time‐dependent non‐linear partial differential equation resulting from a single non‐linear conservation law in h, p, k mathematical and computational framework in which k=(k₁, k₂) are the orders of the approximation spaces in space and time yielding global differentiability of orders (k₁?1) and (k₂?1) in space and time (hence k‐version of finite element method) using space–time marching process. Time‐dependent viscous Burgers equation is used as a specific model problem that has physical mechanism for viscous dissipation and its theoretical solutions are analytic. The inviscid form, on the other hand, assumes zero viscosity and as a consequence its solutions are non‐analytic as well as non‐unique (Russ. Math. Surv. 1962; 17 (3):145–146; Russ. Math. Surv. 1960; 15 (6):53–111). In references (Russ. Math. Surv. 1962; 17 (3):145–146; Russ. Math. Surv. 1960; 15 (6):53–111) authors demonstrated that the solutions of inviscid Burgers equations can only be approached within a limiting process in which viscosity approaches zero. Many approaches based on artificial viscosity have been published to accomplish this including more recent work on H(Div) least‐squares approach (Commun. Pure Appl. Math. 1965; 18 :697–715) in which artificial viscosity is a function of spatial discretization, which diminishes with progressively refined discretizations. The thrust of the present work is to point out that: (1) viscous form of the Burgers equation already has the essential mechanism of viscosity (which is physical), (2) with progressively increasing Reynolds (Re) number (thereby progressively reduced viscosity) the solutions approach that of the inviscid form, (3) it is possible to compute numerical solutions for any Re number (finite) within hpk framework and space–time least‐squares processes, (4) the space–time residual functional converges monotonically and that it is possible to achieve the desired accuracy, (5) space–time, time marching processes utilizing a single space–time strip are computationally efficient. It is shown that viscous form of the Burgers equation without linearizing provides a physical and viablemechanism for approaching the solutions of inviscid form with progressively increasing Re. Numerical studies are presented and the computed solutions are compared with published work. Copyright © 2008 John Wiley & Sons, Ltd.     

Pullback attractor of 2D nonautonomous <Emphasis Type="Italic">g</Emphasis>-Navier-Stokes equations with linear dampness

The pullback attractors for the 2D nonautonomous g-Navier-Stokes equations with linear dampness are investigated on some unbounded domains. The existence of the pullback attractors is proved by verifying the existence of pullback D-absorbing sets with cocycle and obtaining the pullback D-asymptotic compactness. Furthermore, the estimation of the fractal dimensions for the 2D g-Navier-Stokes equations is given.     

An optimized Schwarz method with two‐sided Robin transmission conditions for the Helmholtz equation

Optimized Schwarz methods are working like classical Schwarz methods, but they are exchanging physically more valuable information between subdomains and hence have better convergence behaviour. The new transmission conditions include also derivative information, not just function values, and optimized Schwarz methods can be used without overlap. In this paper, we present a new optimized Schwarz method without overlap in the 2d case, which uses a different Robin condition for neighbouring subdomains at their common interface, and which we call two‐sided Robin condition. We optimize the parameters in the Robin conditions and show that for a fixed frequency an asymptotic convergence factor of 1 – O(h^1/4) in the mesh parameter h can be achieved. If the frequency is related to the mesh parameter h, h = O(1/ω^γ) for γ?1, then the optimized asymptotic convergence factor is 1 – O(ω^(1–2γ)/8). We illustrate our analysis with 2d numerical experiments. Copyright © 2007 John Wiley & Sons, Ltd.     

On coupling the Reynolds‐averaged Navier–Stokes equations with two‐equation turbulence model equations

Two methods for coupling the Reynolds‐averaged Navier–Stokes equations with the q–ω turbulence model equations on structured grid systems have been studied; namely a loosely coupled method and a strongly coupled method. The loosely coupled method first solves the Navier–Stokes equations with the turbulent viscosity fixed. In a subsequent step, the turbulence model equations are solved with all flow quantities fixed. On the other hand, the strongly coupled method solves the Reynolds‐averaged Navier–Stokes equations and the turbulence model equations simultaneously. In this paper, numerical stabilities of both methods in conjunction with the approximated factorization‐alternative direction implicit method are analysed. The effect of the turbulent kinetic energy terms in the governing equations on the convergence characteristics is also studied. The performance of the two methods is compared for several two‐ and three‐dimensional problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Uniformly valid solution of limit cycle of the Duffing–van der Pol equation

The limit cycle of the Duffing–van der Pol equation is studied. By considering the product of the frequency ω of the limit cycle and the coefficient ε as an independent parameter μ=εω, an equivalent equation is obtained and then solved by Liao’s homotopy analysis method. The frequency ω is deduced as a function of μ and δ. This function provides us with an algebraic equation for ω, according to which we have an analytical approximation for the frequency. Numerical examples show that the attained approximation is very accurate. More importantly, the results are uniformly valid for all positive values of ε.     

Positive‐definite q‐families of continuous subcell Darcy‐flux CVD(MPFA) finite‐volume schemes and the mixed finite element method

A new family of locally conservative cell‐centred flux‐continuous schemes is presented for solving the porous media general‐tensor pressure equation. A general geometry‐permeability tensor approximation is introduced that is piecewise constant over the subcells of the control volumes and ensures that the local discrete general tensor is elliptic. A family of control‐volume distributed subcell flux‐continuous schemes are defined in terms of the quadrature parametrization q (Multigrid Methods. Birkhauser: Basel, 1993; Proceedings of the 4th European Conference on the Mathematics of Oil Recovery, Norway, June 1994; Comput. Geosci. 1998; 2 :259–290), where the local position of flux continuity defines the quadrature point and each particular scheme. The subcell tensor approximation ensures that a symmetric positive‐definite (SPD) discretization matrix is obtained for the base member (q=1) of the formulation. The physical‐space schemes are shown to be non‐symmetric for general quadrilateral cells. Conditions for discrete ellipticity of the non‐symmetric schemes are derived with respect to the local symmetric part of the tensor. The relationship with the mixed finite element method is given for both the physical‐space and subcell‐space q‐families of schemes. M‐matrix monotonicity conditions for these schemes are summarized. A numerical convergence study of the schemes shows that while the physical‐space schemes are the most accurate, the subcell tensor approximation reduces solution errors when compared with earlier cell‐wise constant tensor schemes and that subcell tensor approximation using the control‐volume face geometry yields the best SPD scheme results. A particular quadrature point is found to improve numerical convergence of the subcell schemes for the cases tested. Copyright © 2007 John Wiley & Sons, Ltd.     

A compact numerical algorithm for solving the time‐dependent mild slope equation

The mild slope equation has been widely used to describe combined wave refraction and diffraction. In this study, a new numerical algorithm is developed to solve the time‐dependent mild slope equation in a second‐order hyperbolic form. The numerical algorithm is based on a compact and explicit finite difference method that is second‐order accurate in both time and space. The algorithm has the similar structure to the leap‐frog method but is constructed on three time levels for the second‐order time derivative term. The numerical model has the capability of simulating transient wave motion by correctly predicting the speed of wave energy propagation, which is important for the real‐time forecast of the arrival time of storm waves generated in the far field. The model is validated against analytical solution for wave shoaling and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope (Coastal Eng. 1982; 6 :255). Lastly, the realistic scale Homma's island (Geophys. Mag. 1950; 21 :199) is studied with the use of various wave periods of T = 720s, T = 120 s, and T = 24 s. These wave periods correspond to long, intermediate, and short waves for the given topography, respectively. Comparisons are made between numerical results and existing analytical solutions in terms of the wave amplification around the island, which serves as the indicator for the potential wave runup. Excellent agreements are obtained. The model runs on a PC (Pentium IV 1.8GHz) and the computer capacity allows the computation of a mesh system up to 3000 × 3000, which is equivalent to about 150 × 150 waves or a large area of 540km × 540km for a wave train with the period of T = 60 s. Copyright 2004 John Wiley & Sons, Ltd.     

On parallel pre‐conditioners for pressure Poisson equation in LES of complex geometry flows

This paper presents an assessment of fast parallel pre‐conditioners for numerical solution of the pressure Poisson equation arising in large eddy simulation of turbulent incompressible flows. Focus is primarily on the pre‐conditioners suitable for domain decomposition based parallel implementation of finite volume solver on non‐uniform structured Cartesian grids. Bi‐conjugate gradient stabilized method has been adopted as the Krylov solver for the linear algebraic system resulting from the discretization of the pressure Poisson equation. We explore the performance of multigrid pre‐conditioner for the non‐uniform grid and compare its performance with additive Schwarz pre‐conditioner, Jacobi and SOR(k) pre‐conditioners. Numerical experiments have been performed to assess the suitability of these pre‐conditioners for a wide range of non‐uniformity (stretching) of the grid in the context of large eddy simulation of a typical flow problem. It is seen that the multigrid preconditioner shows the best performance. Further, the SOR(k) preconditioner emerges as the next best alternative. Copyright © 2016 John Wiley & Sons, Ltd.     

On the modeling and simulation of a laser‐induced cavitation bubble

Single cavitation bubbles exhibit severe modeling and simulation difficulties. This is due to the small scales of time and space as well as due to the involvement of different phenomena in the dynamics of the bubble. For example, the compressibility, phase transition, and the existence of a noncondensable gas inside the bubble have strong effects on the dynamics of the bubble. Moreover, the collapse of the bubble involves the occurrence of critical conditions for the pressure and temperature. This adds extra difficulties to the choice of equations of state. Even though several models and simulations have been used to study the dynamics of the cavitation bubbles, many details are still not clearly accounted for. Here, we present a numerical investigation for the collapse and rebound of a laser‐induced cavitation bubble in liquid water. The compressibility of the liquid and vapor are involved. In addition, great focus is devoted to study the effects of phase transition and the existence of a noncondensable gas on the dynamics of the collapsing bubble. If the bubble contains vapor only, we use the six‐equation model for two‐phase flows that was modified in our previous work [A. Zein, M. Hantke, and G. Warnecke, J. Comput. Phys., 229(8):2964‐2998, 2010]. This model is an extension to the six‐equation model with a single velocity of Kapila et al. (Phys. Fluid, 13:3002‐3024, 2001) taking into account the heat and mass transfer. To study the effect of a noncondensable gas inside the bubble, we add a third phase to the original model. In this case, the phase transition is considered only at interfaces that separate the liquid and its vapor. The stiffened gas equations of state are used as closure relations. We use our own method to determine the parameters to obtain reasonable equations of state for a wide range of temperatures and make them suitable for the phase transition effects. We compare our results with experimental ones. Also our results confirm some expected physical phenomena. Copyright © 2013 John Wiley & Sons, Ltd.     

Pseudo‐characteristic formulation and dynamic boundary conditions for computational aeroacoustics

The present paper deals with the use of the pseudo‐characteristic formulation of the Navier–Stokes and Euler equations recently introduced by Sesterhenn (Comput. Fluid. 2001; 30 :37–67) for the simulation of acoustic wave propagation. The emphasis is put on the formulation of an efficient method on structured curvilinear grids, along with the definition and implementation of efficient boundary conditions. The cases of inflow, outflow, rigid/compliant walls and walls with prescribed impedance are addressed. The proposed boundary conditions are assessed on generic cases. The pseudo‐characteristic formulation enables a straightforward and optimal use of high‐order upwind dispersion‐relation‐preserving schemes, yielding an efficient method. Copyright © 2006 John Wiley & Sons, Ltd.     

The Extinction Behavior of the Solutions for a Class of Reaction-Diffusion Equations

1　ＩｎｔｒｏｄｕｃｔｉｏｎａｎｄＬｅｍｍａｓＴｈｅｒｅａｒｅｍａｎｙｒｅｓｕｌｔｓａｂｏｕｔｅｘｉｓｔｅｎｃｅ (ｇｌｏｂａｌｏｒｌｏｃａｌ)ａｎｄａｓｙｍｐｔｏｔｉｃｂｅｈａｖｉｏｒｏｆｓｏｌｕｔｉｏｎｓｆｏｒｒｅａｃｔｉｏｎ＿ｄｉｆｆｕｓｉｏｎｅｑｕａｔｉｏｎｓ[1- 9].Ｂｙｔｈｅａｉｄｓｏｆｒｅｓｕｌｔｓ[2 ,3]ｏｆｅｑｕａｔｉｏｎ ｕ/ ｔ=Δｕ-λ｜ｕ｜γ- 1ｕｗｉｔｈｉｎｉｔｉａｌ＿ｂｏｕｎｄａｒｙｖａｌｕｅｓ,ｐａｐｅｒ [4 ]ｓｔｕｄｉｅｄｔｈｅｐｒｏｂｌｅｍｏｆ ｕ/ ｔ=Δｕ-λ｜ｅβｔｕ｜γ- …     

An analysis of the three-dimensional elastic solid with internal rectangular crack

In this paper, the three-dimensional elastic solid with internal rectangular crack is considered. Let the crack surfaces be subjected to the equal and opposite normal tractionsp ₀. This problem is reduced by means of Fourier transforms to the standard set of dual integral equations with two variables. Then the formulas of analytic solution of the displacements on the crack surfaces and of the stress-intensity factors of crack border are obtained.     

An upwind finite‐volume element scheme and its maximum‐principle‐preserving property for nonlinear convection–diffusion problem

For a class of nonlinear convection–diffusion equation in multiple space dimensions, a kind of upwind finite‐volume element (UFVE) scheme is put forward. Some techniques, such as calculus of variations, commutating operators and prior estimates, are adopted. It is proved that the UFVE scheme is unconditionally stable and satisfies maximum principle. Optimal‐order estimates in H¹‐norm are derived to determine the error in the approximate solution. Numerical results are presented to observe the performance of the scheme. Copyright © 2007 John Wiley & Sons, Ltd.     

Inviscid evolution of incompressible swirling flows in pipes: The dependence of the flow structure upon the inlet velocity field

This paper analyses the influence of the inlet swirl on the structure of incompressible inviscid flows in pipes. To that end, the inviscid evolution along a pipe of varying radius with a central body situated inside the pipe is studied for three different inlet swirling flows by solving the Bragg–Hawthorne equation both asymptotically and numerically. The downstream structure of the flow changes abruptly above certain threshold values of the swirl parameter (L). In particular, there exist a value L_r above which a near-wall region of flow reversal is formed downstream, and a critical value L_f above which the axial vortex flow breaks down. It is shown that the dependence upon the pipe geometry of these critical values of the swirl parameter varies strongly with the inlet azimuthal velocity profile considered. An excellent agreement between asymptotic and numerical results is found.