Numerical study of Sivashinsky equation using a splitting scheme based on Crank‐Nicolson method

The mathematical modeling of a planar solid‐liquid interface in the solidification of a dilute binary alloy is formulating by one of nonintegrable, nonlinear evolution equation known as Sivashinsky equation. In the first part of this paper, the mathematical modeling of Sivashinsky equation is briefly discussed. Since, the exact solutions of this equation is yet unknown, obtaining its numerical solution plays an important role to simulate its behavior. Therefore, in the second part, a second‐order splitting finite difference scheme, based on Crank‐Nicolson method, is investigated to approximate the solution of the Sivashinsky equation with homogeneous boundary conditions. We prove the solvability of the present scheme and establish the error estimate of the numerical scheme.     

A second‐order backward difference time‐stepping scheme for penalized Navier–Stokes equations modeling filtration through porous media

We investigate the stability and convergence of a fully implicit, linearly extrapolated second‐order backward difference time‐stepping scheme for the penalized Navier–Stokes equations modeling filtration through porous media. In the penalization approach, an extended Navier–Stokes equation is used in the entire computational domain with suitable resistance terms to mimic the presence of porous medium. It is widely used as an alternative to the heterogeneous approach in which different types of partial differential equations (PDEs) are used in fluid and porous subregions along with suitable continuity conditions at the interface. However, the introduction of extra resistance terms makes the penalized Navier–Stokes equations more nonlinear. We prove that the linearly extrapolated scheme is unconditionally stable and derive optimal order error estimates without any stability condition. To show feasibility and applicability of the approach, it is used to numerically solve a passive control problem in which flow around a solid body is controlled by adding porous layers on the surface. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 681–705, 2016     

Global sharp interface limit of the Hele–Shaw–Cahn–Hilliard system

In this paper, we focus on a diffuse interface model named by Hele–Shaw–Cahn–Hilliard system, which describes a two‐phase Hele–Shaw flow with matched densities and arbitrary viscosity contrast in a bounded domain. The diffuse interface thickness is measured by ? , and the mobility coefficient (the diffusional Peclet number) is ? ^α . We will prove rigorously that the global weak solutions of the Hele–Shaw–Cahn–Hilliard system converge to a varifold solution of the sharp interface model as ? →0 in the case of 0≤α  < 1. Copyright © 2016 John Wiley & Sons, Ltd.     

Non‐local approach in mathematical problems of fluid–structure interaction

Three‐dimensional mathematical problems of interaction between elastic and scalar oscillation fields are investigated. An elastic field is to be defined in a bounded inhomogeneous anisotropic body occupying the domain Ω¯₁⊂ℝ³ while a physical (acoustic) scalar field is to be defined in the exterior domain Ω¯₂=ℝ³\Ω₁ which is filled up also by an anisotropic (fluid) medium. These two fields satisfy the governing equations of steady‐state oscillations in the corresponding domains together with special kinematic and dynamic transmission conditions on the interface ∂Ω₁. The problems are studied by the so‐called non‐local approach, which is the coupling of the boundary integral equation method (in the unbounded domain) and the functional‐variational method (in the bounded domain). The uniqueness and existence theorems are proved and the regularity of solutions are established with the help of the corresponding Steklov–Poincaré type operators and on the basis of the Gårding inequality and the Lax–Milgram theorem. In particular, it is shown that the physical fluid–solid acoustic interaction problem is solvable for arbitrary values of the frequency parameter. Copyright © 1999 John Wiley & Sons, Ltd.     

Criterion of hyperbolicity for non‐conservative quasilinear systems admitting a partially convex conservation law

A system of conservation laws admitting an additional convex conservation law can be written as a symmetric t‐hyperbolic in the sense of Friedrichs system. However, in mathematical modeling of complex physical phenomena, it is customary to use non‐conservative hyperbolic models. We generalize the Godunov–Friedrichs–Lax approach to this new class of models. Copyright © 2011 John Wiley & Sons, Ltd.     

First‐Order system least squares for generalized‐Newtonian coupled Stokes‐Darcy flow

The coupled problem for a generalized Newtonian Stokes flow in one domain and a generalized Newtonian Darcy flow in a porous medium is studied in this work. Both flows are treated as a first‐order system in a stress‐velocity formulation for the Stokes problem and a volumetric flux‐hydraulic potential formulation for the Darcy problem. The coupling along an interface is done using the well‐known Beavers–Joseph–Saffman interface condition. A least squares finite element method is used for the numerical approximation of the solution. It is shown that under some assumptions on the viscosity the error is bounded from above and below by the least squares functional. An adaptive refinement strategy is examined in several numerical examples where boundary singularities are present. Due to the nonlinearity of the problem a Gauss–Newton method is used to iteratively solve the problem. It is shown that the linear variational problems arising in the Gauss–Newton method are well posed. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1150–1173, 2015     

A model for resistance welding including phase transitions and Joule heating

In this paper, we introduce a new model for solid–liquid phase transitions triggered by Joule heating as they arise in the case of resistance welding of metal parts. The main novelties of the paper are the coupling of the thermistor problem with a phase‐field model and the consideration of phase‐dependent physical parameters through a mixture ansatz. The PDE system resulting from our modeling approach couples a strongly nonlinear heat equation, a non‐smooth equation for the the phase parameter (standing for the local proportion of one of the two phases) with a quasistatic electric charge conservation law. We prove the existence of weak solutions in the three‐dimensional (3D) case, whereas the regularity result and the uniqueness of solution is stated only in the two‐dimensional case. Indeed, uniqueness for the 3D system is still an open problem. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling of cutting force prediction in machining high-volume SiCp/Al composites

Variations in cutting forces significantly influence the tool wear and part quality in machining high-volume SiC-particle-reinforced aluminum matrix (SiCp/Al) composites. Properties of the reinforcement SiC particles, such as size and volume fraction, contribute to the change in the cutting forces. This paper presents a cutting force model based on the geometrical and mechanical nature of the tool and workpiece, considering the effect of the SiC reinforcement particles. The cutting force is predicted as three components (F_z, F_x, and F_y) and the resultant cutting force F_τs. The cutting force was considered to generate three deformed zones: (a) shear deformed zone, (b) friction deformed zone on the chip–tool interface, (c) plow deformed zone. The effect of SiC reinforcement particles on friction deformed zone is analyzed emphatically. The friction force from friction deformed zone was obtained by calculating the sliding friction force and rolling friction force. To verify the feasibility and validity of the predicted model of cutting force, cutting experiments were performed with different combinations of cutting speed, feed rate, depth of cut, and tool nose radius. The predicted cutting force values demonstrate good agreement with the measured experimental cutting force values in most cutting conditions. The average percentages of the prediction error were 1.93%, 6.20%, and 10.48% for the F_z, F_x, and F_y components, respectively, thus proving the validity and accuracy of the predicted model of cutting forces.     

Parametrices of mixed elliptic problems

Mixed elliptic problems for differential operators A in a domain Ω with smooth boundary Y are studied in theform $$ Au = f \;\; {\rm in} \;\; {\rm \Omega} \, , \quad T_{\pm}u = g_{\pm} \;\; {\rm on} \;\; Y_{\pm} \, , $$ where Y_± ? Y are manifolds with a common boundary Z, such that Y_– ∪ Y₊ = Y and Y_– ∩ Y₊ = Z, with boundary conditions T_± on Y_± (with smooth coefficients up to Z from the respective side) satisfying the Shapiro–Lopatinskij condition. We consider such problems in standard Sobolev spaces and characterise natural extra conditions on the interface Z with an analogue of Shapiro–Lopatinskij ellipticity for an associated transmission problem on the boundary; then the extended operator is Fredholm. The transmission operators on the boundary with respect to Z belong to a complete pseudo‐differential calculus, a modification of the algebra of boundary value problems without the transmission property. We construct parametrices of elliptic elements in that calculus, and we obtain parametrices of the original mixed problems under additional conditions on the interface. We consider the Zaremba problem and other mixed problems for the Laplace operator, determine the number of extra conditions and calculate the index of associated Fredholm operators. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Complexity of clique‐coloring odd‐hole‐free graphs

In this paper we investigate the problem of clique‐coloring, which consists in coloring the vertices of a graph in such a way that no monochromatic maximal clique appears, and we focus on odd‐hole‐free graphs. On the one hand we do not know any odd‐hole‐free graph that is not 3‐clique‐colorable, but on the other hand it is NP‐hard to decide if they are 2‐clique‐colorable, and we do not know if there exists any bound k₀ such that they are all k₀ ‐clique‐colorable. First we will prove that (odd hole, codiamond)‐free graphs are 2‐clique‐colorable. Then we will demonstrate that the complexity of 2‐clique‐coloring odd‐hole‐free graphs is actually Σ₂ P‐complete. Finally we will study the complexity of deciding whether or not a graph and all its subgraphs are 2‐clique‐colorable. © 2009 Wiley Periodicals, Inc. J Graph Theory 62: 139–156, 2009     

Automated multi‐level sub‐structuring for fluid–solid interaction problems

The automated multi‐level sub‐structuring (AMLS) method is a powerful technique to determine a large number of eigenpairs with moderate accuracy of huge symmetric and definite eigenvalue problems in structural analysis. This paper is concerned with an adapted version of AMLS for eigenfrequency analysis of fluid–solid interaction systems. Although fluid–solid vibrations are governed by an unsymmetric eigenproblem, the modified AMLS method needs approximately the same computational effort. An error bound related to the eigenvalue approximations is proved. Copyright © 2010 John Wiley & Sons, Ltd.     

A decoupling two‐grid algorithm for the mixed Stokes‐Darcy model with the Beavers‐Joseph interface condition

In this article, a decoupling scheme based on two‐grid finite element for the mixed Stokes‐Darcy problem with the Beavers‐Joseph interface condition is proposed and investigated. With a restriction of a physical parameter α, we derive the numerical stability and error estimates for the scheme. Numerical experiments indicate that such two‐grid based decoupling finite element schemes are feasible and efficient. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1066–1082, 2014     

Derivation and analysis of an effective model for biofilm growth in evolving porous media

The present article deals with the growth of biofilms produced by bacteria within a saturated porous medium. Starting from the pore‐scale, the process is essentially described by attachment/detachment of mobile microorganisms to a solid surface and their ability to build biomass. The increase in biomass on the surface of the solid matrix changes the porosity and impedes flow through the pores. Using formal periodic homogenization, we derive an averaged model describing the process via Darcy's law and upscaled transport equations with effective coefficients provided by the evolving microstructure at the pore‐scale. Assuming, that the underlying pore geometry may be described by a single parameter, for example, porosity, the level set equation locating the biofilm‐liquid interface transforms into an ordinary differential equation (ODE) for the parameter. For such a setting, we state significant analytical and algebraic properties of these effective parameters. A further objective of this article is the analytical investigation of the resulting coupled PDE–ODE model. In a weak sense, unique solvability either global in time or at least up to a possible clogging phenomenon is shown. Copyright © 2016 John Wiley & Sons, Ltd.     

s‐Regular cubic graphs as coverings of the complete bipartite graph K3,3

A graph is s‐regular if its automorphism group acts freely and transitively on the set of s‐arcs. An infinite family of cubic 1‐regular graphs was constructed in [10], as cyclic coverings of the three‐dimensional Hypercube. In this paper, we classify the s‐regular cyclic coverings of the complete bipartite graph K_3,3 for each ≥ 1 whose fibre‐preserving automorphism subgroups act arc‐transitively. As a result, a new infinite family of cubic 1‐regular graphs is constructed. © 2003 Wiley Periodicals, Inc. J Graph Theory 45: 101–112, 2004     

Coupling of mixed finite element and stabilized boundary element methods for a fluid–solid interaction problem in 3D

We introduce and analyze the coupling of a mixed finite element and a boundary element for a three‐dimensional time‐harmonic fluid–solid interaction problem. We consider a formulation in which the Cauchy stress tensor and the rotation are the main variables in the elastic structure and use the usual pressure formulation in the acoustic fluid. The mixed variational formulation in the solid is completed with boundary integral equations relating the Cauchy data of the acoustic problem on the coupling interface. A crucial point in our formulation is the stabilization technique introduced by Hiptmair and coworkers to avoid the well‐known instability issue appearing in the boundary element method treatment of the exterior Helmholtz problem. The main novelty of this formulation, with respect to a previous approach, consists in reducing the computational domain to the solid media and providing a more accurate treatment of the far field effect. We show that the continuous problem is well‐posed and propose a conforming Galerkin method based on the lowest‐order Arnold–Falk–Winther mixed finite element. Finally, we prove that the numerical scheme is convergent with optimal order.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1211–1233, 2014     

Large‐time behavior of a two‐scale semilinear reaction–diffusion system for concrete sulfatation

We study the large‐time behavior of (weak) solutions to a two‐scale reaction–diffusion system coupled with a nonlinear ordinary differential equations modeling the partly dissipative corrosion of concrete (or cement)‐based materials with sulfates. We prove that as t → ∞ , the solution to the original two‐scale system converges to the corresponding two‐scale stationary system. To obtain the main result, we make use essentially of the theory of evolution equations governed by subdifferential operators of time‐dependent convex functions developed combined with a series of two‐scale energy‐like time‐independent estimates. Copyright © 2014 John Wiley & Sons, Ltd.     

Macro‐hybrid mixed variational two‐phase flow through fractured porous media

Macro‐hybrid mixed variational models of two‐phase flow, through fractured porous media, are analyzed at the mesoscopic and macroscopic levels. The mesoscopic models are treated in terms of nonoverlapping domain decompositions, in such a manner that the porous rock matrix system and the fracture network interact across rock–rock, rock–fracture, and fracture–fracture interfaces, with flux transmission conditions dualized. Alternatively, the models are scaled to a macroscopic level via an asymptotic process, where the width of the fractures tends to zero, and the fracture network turns out to be an interface system of one less spatial dimension, with variable high permeability. The two‐phase flow is characterized by a fractional flow dual mixed variational model. Augmented two‐field and three‐field variational reformulations are presented for regularization, internal approximations, and macro‐hybrid mixed finite element implementation. Also abstract proximal‐point penalty‐duality algorithms are derived and analyzed for parallel computing. Copyright © 2016 John Wiley & Sons, Ltd.     

The A4‐structure of a graph

We define the A₄‐structure of a graph G to be the 4‐uniform hypergraph on the vertex set of G whose edges are the vertex subsets inducing 2K₂, C₄, or P₄. We show that perfection of a graph is determined by its A₄‐structure. We relate the A₄‐structure to the canonical decomposition of a graph as defined by Tyshkevich [Discrete Math 220 (2000) 201–238]; for example, a graph is indecomposable if and only if its A₄‐structure is connected. We also characterize the graphs having the same A₄‐structure as a split graph.     

Clique‐coloring some classes of odd‐hole‐free graphs

We consider the problem of clique‐coloring, that is coloring the vertices of a given graph such that no maximal clique of size at least 2 is monocolored. Whereas we do not know any odd‐hole‐free graph that is not 3‐clique‐colorable, the existence of a constant C such that any perfect graph is C‐clique‐colorable is an open problem. In this paper we solve this problem for some subclasses of odd‐hole‐free graphs: those that are diamond‐free and those that are bull‐free. We also prove the NP‐completeness of 2‐clique‐coloring K₄‐free perfect graphs. © 2006 Wiley Periodicals, Inc. J Graph Theory 53: 233–249, 2006     

Integrative interaction analysis using threshold gradient directed regularization

For many complex business and industry problems, high‐dimensional data collection and modeling have been conducted. It has been shown that interactions may have important implications beyond the main effects. The number of unknown parameters in an interaction analysis can be larger or much larger than the sample size. As such, results generated from analyzing a single data set are often unsatisfactory. Integrative analysis, which jointly analyzes the raw data from multiple independent studies, has been conducted in a series of recent studies and shown to outperform single–data set analysis, meta‐analysis, and other multi–data set analyses. In this study, our goal is to conduct integrative analysis in interaction analysis. For regularized estimation and selection of important interactions (and main effects), we apply a threshold gradient directed regularization approach. Advancing from the existing studies, the threshold gradient directed regularization approach is modified to respect the “main effects, interactions” hierarchy. The proposed approach has an intuitive formulation and is computationally simple and broadly applicable. Simulations and the analyses of financial early warning system data and news‐APP (application) recommendation behavior data demonstrate its satisfactory practical performance.