Search

Abstract: In (2012), Leonetti and Siepe [10] considered solutions to boundary value problems of some anisotropic elliptic equations of the type Under some suitable conditions, they obtained an integrability result, which shows that higher integrability of the boundary datum forces solutions to have higher integrability as well. In the present paper, we consider -obstacle problems of the nonhomogeneous anisotropic elliptic equations under some controllable growth and monotonicity conditions. We obtain an integrability result, which can be regarded as a generalization of the result due to Leonetti and Siepe. [Copyright &y& Elsevier]

Abstract: The aim of this paper is twofold: First, we characterize an essentially optimal class of boundary operators Θ which give rise to self-adjoint Laplacians in with (nonlocal and local) Robin-type boundary conditions on bounded Lipschitz domains , , . Second, we extend Friedlander''s inequalities between Neumann and Dirichlet Laplacian eigenvalues to those between nonlocal Robin and Dirichlet Laplacian eigenvalues associated with bounded Lipschitz domains Ω, following an approach introduced by Filonov for this type of problems. [Copyright &y& Elsevier]

Abstract: We compare the above-mentioned article with the content of a previous publication. [Copyright &y& Elsevier]

In this paper we consider the vanishing viscosity limit of solutions to the initial boundary value problem for the compressible viscoelastic equations in the half space. When the initial deformation gradient does not degenerate and there is no vacuum initially, we establish the uniform regularity estimates of solutions to the initial-boundary value problem for the three-dimensional compressible viscoelastic equations in the Sobolev spaces. Then we justify the vanishing viscosity limit of solutions of the compressible viscoelastic equations based on the uniform regularity estimates and the compactness arguments. Both the no-slip boundary condition and the Navier-slip type boundary condition on velocity are addressed in this paper. On the one hand, for the corresponding vanishing viscosity limit of the compressible Navier-Stokes equations with the no-slip boundary condition, it is impossible to derive such uniform energy estimates of solutions due to the appearance of strong boundary layers. Consequently, our results show that the deformation gradient can prevent the formation of strong boundary layers. On the other hand, our results also provide two different kinds of boundary conditions suitable for the well-posedness of the initial-boundary value problem of the elastodynamic equations via the method of vanishing viscosity. Finally, it is worth noting that we take advantage of the Lagrangian coordinates to study the vanishing viscosity limit for the fixed boundary problem in this paper. [ABSTRACT FROM AUTHOR]

Ordinary and partial differential equations with operators containing the p -Laplace operator, the so called ϕ -Laplace operator, operators with variable exponents and the double phase operators have been studied thoroughly for at least the last 20 years and the amount of papers dealing with these subjects is huge. In this paper we deal with the problem of the existence of periodic solutions to some system of ODE's involving a fairly general operator generalizing in this form corresponding results for the areas of research just mentioned. In our results neither the function generating the differential operator nor the right hand side function in the boundary value problem considered are necessarily equal to the gradient of a function. Our main result is a here is a continuation theorem based on the Leray-Schauder degree theory which allows us to obtain the existence of solutions for the boundary value problem considered. [ABSTRACT FROM AUTHOR]

In this paper we deal with the following boundary value problem { − Δ p u + g (u) | ∇ u | p = h (u) f in Ω, u ≥ 0 in Ω, u = 0 on ∂Ω, in a domain Ω ⊂ R N (N ≥ 2) , where 1 ≤ p < N , g is a positive and continuous function on [ 0 , ∞) , and h is a continuous function on [ 0 , ∞) (possibly blowing up at the origin). We show how the presence of regularizing terms h and g allows to prove existence of finite energy solutions for nonnegative data f only belonging to L 1 (Ω). [ABSTRACT FROM AUTHOR]

In this paper we consider time-harmonic acoustic wave propagation in a half-plane filled by inhomogeneous periodic medium. If the refractive index depends on the horizontal coordinate only, we define upward and downward radiating modes by solving a one-dimensional Sturm-Liouville eigenvalue problem with a complex-valued periodic coefficient. The upward and downward radiation conditions are introduced based on a generalized Rayleigh series. Using the variational method, we then prove uniqueness and existence for the scattering of an incoming wave mode by a grating located between an upper and lower half plane with such inhomogeneous periodic media. The solution operator of the quasiperiodic boundary value problem yields a well-defined scattering operator, the numerical approximation of which is nothing else but the S-matrix in scattering matrix algorithms like Rigorous Coupled Wave Analysis (RCWA) or Fourier Modal Method (FMM). [ABSTRACT FROM AUTHOR]

The Hunter-Saxton equation models the propagation of weakly nonlinear orientation waves in a massive director field of the nematic liquid crystal. In this paper, we study the vanishing viscosity limit for an initial boundary value problem of the Hunter-Saxton equation with the characteristic boundary condition. By the formal multiscale analysis, we first derive the characteristic boundary layer profile, which satisfies a nonlinear parabolic equation. On the base of the Galerkin method along with a compactness argument, we then establish the global well-posedness of the boundary layer equation. Finally, we prove the global stability of the boundary layer profiles together with the optimal convergence rate of the vanishing viscosity limit by the energy method. [ABSTRACT FROM AUTHOR]

Inverse scattering problem for the operator representing sum of the operator of the third derivative on semi-axis and of the operator of multiplication by a real function is studied in this paper. Properties of Jost solutions of such an operator are studied and it is shown that these Jost solutions are solutions of the Riemann boundary value problem on a system of rays. The main system of linear singular integral equations is derived. This system is equivalent to the solution of inverse scattering problem. [ABSTRACT FROM AUTHOR]

The Hamiltonian systems considered in this paper are obtained by weakly coupling two systems having completely different behaviors. The first one satisfies the twist assumptions usually considered for the application of the Poincaré–Birkhoff Theorem, while the second one presents the existence of some well-ordered lower and upper solutions. In the higher dimensional case, we also treat a coupling situation where the classical Hartman condition is assumed. [ABSTRACT FROM AUTHOR]

This paper is concerned with boundary regularity estimates in the homogenization of elliptic equations with rapidly oscillating and high-contrast coefficients. We establish uniform nontangential-maximal-function estimates for the Dirichlet, regularity, and Neumann problems with L 2 boundary data in a periodically perforated Lipschitz domain. [ABSTRACT FROM AUTHOR]

In this paper, we study the vanishing dissipation limit of solutions to the initial boundary value problem of three-dimensional incompressible viscous magnetohydrodynamic equations in the upper half space where the magnetic field is assumed to be transverse to the boundary. When the tangential component of the magnetic field is imposed the zero Neumann boundary condition, we can establish the uniform energy estimates of the solutions to the initial boundary value problem in the conormal Sobolev spaces even if the velocity field satisfies the no-slip boundary condition. Then the vanishing dissipation limit of solutions to the initial boundary value problem can be achieved in L ∞ sense based on some compactness arguments. Compared with the previous results, our results show that the strong boundary layer can still be prevented by the transverse magnetic field even with the magnetic diffusion. [ABSTRACT FROM AUTHOR]

In this paper, we deal with the free boundary value problem of the three-dimensional full primitive equations for the large-scale ocean with heat insulation in a bounded domain. The domain is periodic in the horizontal direction and finite in the vertical direction, and its upper boundary is a free surface driven by the fluid and bottom boundary is fixed. We show that the unique strong solution of the free boundary value problem for the three-dimensional full primitive equations exists globally in time and converges to the equilibrium state with exponential decay rate, when the initial data is a small perturbation of the steady state in regular Sobolev space. [ABSTRACT FROM AUTHOR]

This paper is concerned with a chemotaxis-fluid model with realistic boundary conditions matching the experiments of Kessler et al. Different from the initial and boundary value problem, the existence of time periodic solution is proved under no any smallness assumptions on chemotactic coefficient. Subsequently, we investigate the reason of pattern formation. It is generally known that there is no pattern formation for the homogeneous boundary value problem. Therefore, the pattern formation is partly caused by non-homogeneous boundary conditions, but chemotaxis is also an indispensable condition. In fact, it is proved that the density of bacteria will tend to be evenly distributed when the chemotactic intensity tends to 0. At last, the idea of finding time periodic pattern is also applied to prove the steady-state pattern formation of chemotaxis-Navier-Stokes model. [ABSTRACT FROM AUTHOR]

In this paper, we are concerned with a diffusive system modeling the interaction of mussel and algae in the water layer overlying the mussel bed, where the algae are the main food source for mussels, and the advection pushes algaes in one direction but not out of the domain. We present a threshold result on the global extinction and persistence of mussels. It shows that the condition for persistence depends on the principal eigenvalue of a scalar boundary value problem, which is related to the diffusion, the speed of the tidal flow, the conversion rate of algae to mussel production, and mussel mortality rate. We further investigate the asymptotic profile of the positive steady state when it exists. [ABSTRACT FROM AUTHOR]

Morse index theory in the classical framework provides an equality between the spectral properties of the second variation of the Lagrangian functional and the oscillation properties of the space of solutions of the associated boundary value problem. For carrying over a similar result in the non-Hamiltonian context which can be useful for investigating the dynamical properties of dissipative systems, we introduce a new topological invariant, the so-called "degree-index" defined in terms of the Brouwer degree of a suitable determinant map of a boundary matrix, which provides one possible substitute of the Maslov index in this non-Hamiltonian framework and finally we prove the equality between the Morse index and the degree-index in this non-selfadjoint setting through a new abstract trace formula. In the last part of the paper we apply our theoretical results to some 1D reaction-diffusion systems. [ABSTRACT FROM AUTHOR]

This paper considers two point boundary value problems for conservative systems defined in multiple coordinate systems, and develops a flexible a posteriori framework for computer assisted existence proofs. Our framework is applied to the study collision and near collision orbits in the circular restricted three body problem. In this case the coordinate systems are the standard rotating coordinates, and the two Levi-Civita coordinate systems regularizing collisions with each of the massive primaries. The proposed framework is used to prove the existence of a number of orbits which have long been studied numerically in the celestial mechanics literature, but for which there are no existing analytical proofs at the mass and energy values considered here. These include transverse ejection/collisions from one primary body to the other, Strömgren's asymptotic periodic orbits (transverse homoclinics for L 4 , 5), families of periodic orbits passing through collision, and orbits connecting L 4 to ejection or collision. [ABSTRACT FROM AUTHOR]

The paper is concerned with the stability property under perturbation Q → Q ˜ of different spectral characteristics of a boundary value problem associated in L 2 ([ 0 , 1 ] ; C 2) with the following 2 × 2 Dirac type equation (0.1) L U (Q) y = − i B − 1 y ′ + Q (x) y = λ y , B = ( b 1 0 0 b 2 ) , b 1 < 0 < b 2 , y = col (y 1 , y 2) , with a potential matrix Q ∈ L p ([ 0 , 1 ] ; C 2 × 2) and subject to the regular boundary conditions U y : = { U 1 , U 2 } y = 0. If b 2 = − b 1 = 1 this equation is equivalent to one dimensional Dirac equation. Our approach to the spectral stability relies on the existence of the triangular transformation operators for system (0.1) with Q ∈ L 1 , which was established in our previous works. The starting point of our investigation is the Lipshitz property of the mapping Q → K Q ± , where K Q ± are the kernels of transformation operators for system (0.1). Namely, we prove the following uniform estimate: ‖ K Q ± − K Q ˜ ± ‖ X ∞ , p 2 + ‖ K Q ± − K Q ˜ ± ‖ X 1 , p 2 ⩽ C ⋅ ‖ Q − Q ˜ ‖ p , Q , Q ˜ ∈ U p , r 2 × 2 , p ∈ [ 1 , ∞ ] , on balls U p , r 2 × 2 in L p ([ 0 , 1 ] ; C 2 × 2). It is new even for Q ˜ = 0. Here X ∞ , p 2 , X 1 , p 2 are the special Banach spaces naturally arising in such problems. We also obtained similar estimates for Fourier transforms of K Q ±. Both of these estimates are of independent interest and play a crucial role in the proofs of all spectral stability results discussed in the paper. For instance, as an immediate consequence of these estimates we get the Lipshitz property of the mapping Q → Φ Q (⋅ , λ) , where Φ Q (x , λ) is the fundamental matrix of the system (0.1). Assuming the spectrum Λ Q = { λ Q , n } n ∈ Z of L U (Q) to be asymptotically simple, denote by F Q = { f Q , n } | n | > N a sequence of corresponding normalized eigenvectors, L U (Q) f Q , n = λ Q , n f Q , n. Assuming boundary conditions (BC) to be strictly regular , we show that the mapping Q → Λ Q − Λ 0 sends L p ([ 0 , 1 ] ; C 2 × 2) either into ℓ p ′ or into the weighted ℓ p -space ℓ p ({ (1 + | n |) p − 2 }) ; we also establish its Lipshitz property on compact sets in L p ([ 0 , 1 ] ; C 2 × 2) , p ∈ [ 1 , 2 ]. The proof of the second estimate involves as an important ingredient inequality that generalizes classical Hardy-Littlewood inequality for Fourier coefficients. It is also shown that the mapping Q → F Q − F 0 sends L p ([ 0 , 1 ] ; C 2 × 2) into the space ℓ p ′ (Z ; C ([ 0 , 1 ] ; C 2) of sequences of continuous vector-functions, and has the Lipshitz property on compacts sets in L p ([ 0 , 1 ] ; C 2 × 2) , p ∈ [ 1 , 2 ]. Certain modifications of these spectral stability results are also proved for balls U p , r 2 × 2 in L p ([ 0 , 1 ] ; C 2 × 2) , p ∈ [ 1 , 2 ]. Note also that the proof of the Lipshitz property of the mapping Q → F Q − F 0 involves the deep Carleson-Hunt theorem for maximal Fourier transform, while the proof of this property for the mapping Q → Λ Q − Λ 0 relies on the estimates of the classical Fourier transform and is elementary in character. We apply our previous results to the damped string equation to establish the Riesz basis property and the asymptotic behavior of the eigenvalues of the corresponding dynamic generator under the assumptions d ∈ L 1 [ 0 , ℓ ] , ρ ∈ W 1 , 1 [ 0 , ℓ ] on the damping coefficient and the density of the string, that are weaker than previously treated in the literature. We also establish Lipshitz dependence on d and ρ in ℓ p -spaces of the remainders in the asymptotic formula for the eigenvalues. [ABSTRACT FROM AUTHOR]

Direct and inverse scattering problem for an operator with non-local potential is solved in the paper. The method is based on the Riemann boundary value problem on a bundle of three straight lines. Description of scattering problem data is given. [ABSTRACT FROM AUTHOR]

In the paper we established the L 2 estimates for commutators of the Dirichlet-to-Neumann map generated by a divergence form parabolic operator, defined on R + n + 2 , with real, symmetric and t , λ − independent coefficient matrix, or more generally, a small complex L ∞ perturbation of such. The major new challenge, compared to our previous work in elliptic setting, is to handle the first order derivatives with respect to the time variable t which not only fall on solutions to parabolic equations but also on auxiliary functions. [ABSTRACT FROM AUTHOR]

The initial boundary value problems for compressible Navier-Stokes-Poisson and Euler-Poisson equations of plasma are considered on exterior domains in this paper. With the radial symmetry assumption, the global existence of solutions to compressible Navier-Stokes-Poisson equations with the large initial data on a domain exterior to a ball in R n (n ≥ 1) is proved. Moreover, without any symmetry assumption, the global existence of smooth solutions near a given constant steady state for both compressible Navier-Stokes-Poisson and Euler-Poisson equations on an exterior domain in R 3 with physical boundary conditions is also established with the exponential stability. A key issue addressed in this paper is on the global-in-time regularity of solutions near physical boundaries. This is in particular so for the 3-D compressible Navier-Stokes-Poisson equations to which global smooth solutions of initial boundary value problems are seldom found in literature to the best of knowledge. [ABSTRACT FROM AUTHOR]

Dynamics of the macroscopic motion of a compressible, viscous and heat-conducting, one-dimensional monatomic ionized gas can be described by a compressible Navier-Stokes type system. Compared with the corresponding compressible Navier-Stokes type systems modeling a viscous and heat conducting ideal polytropic gas [7,10,13–16,18,19,23,24,26,28,31,32] , the real gases [17] , nonlinear thermoviscoelasticity [5,6] , and a viscous radiative and reactive gas [11,12,20,21,30] , the dependence of the state of an ionized gas on the degree of ionization leads to the loss of concavity of the physical entropy in some small bounded domain. Such a property of the ionized gas makes it impossible to deduce the dissipative estimates on the first-order spatial derivative of the bulk velocity and the absolute temperature from the basic energy estimates obtained through the normalized entropy and, as a consequence, makes the construction of its global large solutions and the study of their large time behavior more difficult. By making full use of the intrinsic structure of equations, the existence of global large smooth solutions to two types of initial boundary value problems of one-dimensional, compressible Navier-Stokes system for a viscous and heat-conducting ionized gas is established in [22] for a class of density and/or temperature dependent transport coefficients, but since the estimates obtained there depend on the time variable, the problem on their large time behavior remains unsolved and the main purpose of this paper is to give the precise description of the large time behavior of the global solutions constructed in [22]. The key ingredient in our analysis is to deduce the uniform-in-time positive lower and upper bounds on the specific volume and the absolute temperature. [ABSTRACT FROM AUTHOR]

This paper concerns the Ginzburg-Landau equation and the Meissner equations on a bounded multiply-connected domain in R 3. The novelty of these equations is that both of them contain an unknown potential and an unknown Neumann field as the Lagrange multipliers. The minimal and non-minimal solutions for each equation are obtained by variational methods, and the effect of the applied magnetic field on the Meissner states is examined. The Meissner A-system with potential is also examined and a solution is obtained by using the Mountain Pass Lemma of Ambrosetti and Rabinowitz. In our analysis the estimates of the natural boundary value problem of the Maxwell-Stokes system play an important role. [ABSTRACT FROM AUTHOR]

Developing surrogate models for costly mathematical models representing physical systems is challenging since it is typically not possible to generate large training data sets, i.e. to create a large experimental design. In such cases, it can be beneficial to constrain the surrogate approximation to adhere to the known physics of the model. This paper presents a novel methodology for the construction of physics-informed polynomial chaos expansions (PCE) that combines the conventional experimental design with additional constraints from the physics of the model represented by a set of differential equations and specified boundary conditions. A computationally efficient means of constructing physically constrained PCEs, termed PC2, are proposed and compared to the standard sparse PCE. Algorithms are presented for both full-order and sparse PC2 expansions and an iterative approach is proposed for addressing nonlinear differential equations. It is shown that the proposed algorithms lead to superior approximation accuracy and do not add significant computational burden over conventional PCE. Although the main purpose of the proposed method lies in combining training data and physical constraints, we show that the PC2 can also be constructed from differential equations and boundary conditions alone without requiring model evaluations. We further show that the constrained PCEs can be easily applied for uncertainty quantification through analytical post-processing of a reduced PCE by conditioning on the deterministic space-time variables. Several deterministic examples of increasing complexity are provided and the proposed method is demonstrated for uncertainty quantification. • PC2 – a novel framework for physically constrained polynomial chaos expansions is proposed. • An efficient algorithm based on constrained least squares and sparse solver is developed. • Analytical uncertainty quantification of approximated PDEs by PC2 is performed. [ABSTRACT FROM AUTHOR]

The Immersed Boundary (IB) method of Peskin (1977) [1] is useful for problems that involve fluid-structure interactions or complex geometries. By making use of a regular Cartesian grid that is independent of the geometry, the IB framework yields a robust numerical scheme that can efficiently handle immersed deformable structures. Additionally, the IB method has been adapted to problems with prescribed motion and other PDEs with given boundary data. IB methods for these problems traditionally involve penalty forces which only approximately satisfy boundary conditions, or they are formulated as constraint problems. In the latter approach, one must find the unknown forces by solving an equation that corresponds to a poorly conditioned first-kind integral equation. This operation can therefore require a large number of iterations of a Krylov method, and since a time-dependent problem requires this solve at each step in time, this method can be prohibitively inefficient without preconditioning. This work introduces a new, well-conditioned IB formulation for boundary value problems, called the Immersed Boundary Double Layer (IBDL) method. In order to lay the groundwork for similar formulations of Stokes and Navier-Stokes equations, this paper focuses on the Poisson and Helmholtz equations to introduce the methodology and to demonstrate its efficiency over the original constraint method. In this double layer formulation, the equation for the unknown boundary distribution corresponds to a well-conditioned second-kind integral equation that can be solved efficiently with a small number of iterations of a Krylov method without preconditioning. Furthermore, the iteration count is independent of both the mesh size and spacing of the immersed boundary points. The method converges away from the boundary, and when combined with a local interpolation, it converges in the entire PDE domain. Additionally, while the original constraint method applies only to Dirichlet problems, the IBDL formulation can also be used for Neumann boundary conditions. • Reformulated Immersed Boundary method for double layer integral equations. • Efficient IB constraint method for PDEs with prescribed boundary values. • Well-conditioned IB method resulting from second-kind integral equation formulation. • IB method for Neumann boundary conditions. [ABSTRACT FROM AUTHOR]

In this paper, we study the asymptotic behavior of the solutions to an initial boundary value problem of the Hunter-Saxton equation in the half space when the viscosity tends to zero. By means of the asymptotic analysis with multiple scales, we first formally derive the equations for boundary layer profiles. Next, we study the well-posedness of the equations for the boundary layer profiles by using the compactness argument. Moreover, we construct an accurate approximate solution and use the energy method to obtain the convergence results of the vanishing viscosity limit. [ABSTRACT FROM AUTHOR]

This paper is concerned with the local well-posedness for the free boundary value problem of smooth solutions to the cylindrical symmetric Euler equations with damping and related models, including the compressible Euler equations and the Euler-Poisson equations. The free boundary is moving in the radial direction with the radial velocity, which will affect the angular velocity but does not affect the axial velocity. However, the compressible Euler equations or Euler-Poisson equations with damping become a degenerate system at the moving boundary. By setting a suitable weighted Sobolev space and using Hardy's inequality, we successfully overcome the singularity at the center point and the vacuum occurring on the moving boundary, and obtain the well-posedness of local smooth solutions. We also summarize the recent related results on the free boundary value problem for the Euler equations with damping, compressible Euler equations and Euler-Poisson equations. [ABSTRACT FROM AUTHOR]

Highlights • Pipeline fixing method may significantly influence the course of water hammer. • For boundary condition (BC) at pipeline support the stiffness matrix is required. • Flexibility matrices of flat springs used for actual pipe supports were determined. • Analytical, experimental and numerical results were received and compared. • BC at support was modeled for singular as well as nonsingular flexibility matrix. Abstract Elastic mounting elements are often used in mechanical systems for reduction of undesired loads. In many engineering applications one-dimensional (1D) approximation of absorber reaction is an acceptable approach however in general, a six order stiffness matrix of an elastic element should be known. This can also be the case of water hammer with fluid-structure interaction in a pipeline fixed with elastic supports. For the standard model of this phenomenon the general boundary condition (BC) at the support should take into account all degrees of freedom of the pipe element. In this paper various methods of flexibility matrix determination of actual, self manufactured flat springs (FS) used as supports at a laboratory pipeline are discussed. Analytical calculations are based on modeling of elastic strain energy of a bar under external loads and formulas being developed for the FS flexibility coefficients are presented. Experimental verification of selected coefficients was performed and the results were consistent with theoretical findings. A general experimental method and a data processing scheme based on linear least squares (LS) method is also proposed and shortly discussed. Numerical investigations are based on FEM modeling of the FS structure and computations of their load-displacement characteristics with the ANSYS software. Flexibility matrix was determined with that same LS data processing scheme. Analytical, experimental and numerical results are compared and discussed. For non-singular flexibility matrix the stiffness matrix used in the BC can be easily determined as its reverse. Otherwise, for a singular case, a special treatment is required. A solution to this problem is developed and presented in the paper. It has been evidenced that singularity results in dimension reduction of the BC. [ABSTRACT FROM AUTHOR]

Highlights • This paper presents an original updating method that enables to update the model parameters one-at-a-time. • For each model parameter, the model outputs are transformed into a new output which is sensitive to this parameter only. • The sensitivity is quantified using Sobol indices. • The proposed approached is robust with respect to modelling errors. Abstract This paper is interested in model updating problems which consists in identifying optimal values of model parameters by comparing the model outputs with the experimental outputs. Such a problem generally yields a challenging multivariate inverse problem to be solved in high dimension. The high-dimensionality requires the use of a global optimization algorithm in order the explore efficiently the parameters space. In this paper we propose an alternative algorithm which allows each model parameters to be identified separately and sequentially by solving separated univariate inverse problems. For each parameter, a devoted inverse problem is introduced by identifying an output which is sensitive to this parameter only, the sensitivity being quantified using Sobol indices. The proposed method is illustrated through a three-storey structure for which experimental measurements are collected. [ABSTRACT FROM AUTHOR]

Abstract This paper proposes a numerical method to model wave scattering by local inhomogeneities in elastic open waveguides. The damaged zone of the waveguide and its vicinity are described by a finite-element model. This model is coupled on the cross-section boundaries to a modal representation of the field propagating along the waveguide axis in the undamaged parts, yielding transparent boundary conditions. However, open waveguides are unbounded in the transverse direction, which complicates both the numerical resolution and the physical analysis. Theoretically, the wave fields are described by a discrete sum of trapped modes and two continuous sums on radiation modes. The latter is sometimes approximated by a discrete sum of leaky modes, which grow at infinity in the transverse direction. In this paper, a Perfectly Matched Layer (PML) of finite thickness is introduced to absorb outgoing waves in the transverse direction, yielding three types of discrete modes: trapped modes, leaky modes and PML modes (PML modes oscillate mainly inside the layer and are non-intrinsic to the physics). Two numerical test cases are considered: the scattering at the junction between a closed and an open cylindrical waveguides and the scattering by an axisymmetrical notch. Good agreement with literature results is found. In particular, the influence of PML modes on the scattered field is highlighted through numerical tests. Due to the lack of power orthogonality of leaky modes, it is also shown that the modal cross power can become significant, which complicates the scattering analysis of open waveguides. Finally, the generality of the proposed method is discussed through a three-dimensional test case considering an embedded bar with an oblique break. Graphical abstract Image 1 [ABSTRACT FROM AUTHOR]

Abstract It is well-known that the convective wave equation with a uniform mean flow can be transformed into a standard wave equation without flow by a Prandtl-Glauert-Lorentz type transformation. This paper examines the boundary condition to be used in the transformed coordinates for acoustic scattering by rigid bodies. Recently, based on the energy conservation equation for acoustic waves propagating in a uniform flow, a new solid wall boundary condition, dubbed the Zero Energy Flux (ZEF) solid surface boundary condition, has been proposed. Instead of the commonly used solid surface boundary condition that the normal acoustic velocity be zero, the ZEF condition requires that the acoustic energy flux be zero. In this paper, we point out that when formulated in the acoustic velocity potential and under the ZEF boundary condition, the normal derivative of the velocity potential remains to be zero when computed in the transformed coordinates. As such, numerical methods developed for solving the standard wave equation can be directly applied to the convective wave equation through a use of the Prandtl-Glauert-Lorentz transformation. Utilization of the transformation is particularly beneficial for numerical methods based on the boundary integral formulation because of the simplified kernel function in the transformed coordinates, for both the time domain and the frequency domain approaches. Useful relations for applying the transformation are derived and numerical examples that demonstrate the effectiveness of the transformation and ZEF boundary condition are presented. [ABSTRACT FROM AUTHOR]

Abstract In this paper, we investigate the heat transport across Germanium Silicon interface using an enhanced ballistic-diffusive equation (EBDE), where we introduce the temperature jump boundary condition coupled with a thermal boundary resistance (TBR) and an effective thermal conductivity model. This paper focuses on the thermal transport of sub-23 nm Ge/Si thin films. The present model is aimed to describe the ballistic-diffusive phonon transport across Ge/Si interface. We have found that the temperature jump occurs in the interface due to phonon-boundary interactions. In addition, the interfacial heat transport is influenced by the surface roughness effect. The prediction of the suggested EBDE model are in good accordance with analytical method reported in the literature. The proposed model shows excellent agreement with the phonon Boltzmann transport equation (BTE) approach and Monte Carlo simulation (MC). Further, the analytic model for the effective thermal conductivity (ETC) is in strong agreement with experimentally based approach and also the theoretical model. Graphical abstract Image Highlights • An enhanced ballistic-diffusive model is developed for heat transport across the interface scattering between nanomaterials. • An excellent agreement with experimental results and Monte Carlo simulation. • Investigation of phonon-surface scattering in nanostructures. • This work help research optimization of interface quality. [ABSTRACT FROM AUTHOR]

Abstract In this paper we study bifurcation solutions from the unique radial solution of a free boundary problem modeling stationary state of tumors with angiogenesis. This model comprises two elliptic equations describing the distribution of the nutrient concentration σ = σ (x) and the inner pressure p = p (x). Unlike similar tumor models that have been intensively studied in the literature where Dirichlet boundary condition for σ is imposed, in this model the boundary condition for σ is a Robin boundary condition. Existence and uniqueness of a radial solution of this model have been successfully proved in a recently published paper [20]. In this paper we study existence of nonradial solutions by using the bifurcation method. Let { γ k } k = 2 ∞ be the sequence of eigenvalues of the linearized problem. We prove that there exists a positive integer k ⁎ ⩾ 2 such that in the two dimension case for any k ⩾ k ⁎ , γ k is a bifurcation point, and in the three dimension case for any even k ⩾ k ⁎ , γ k is also a bifurcation point. [ABSTRACT FROM AUTHOR]

The aim of this paper is to study the problem { u t t − Δ u + P ( x , u t ) = f ( x , u ) in ( 0 , ∞ ) × Ω , u = 0 on ( 0 , ∞ ) × Γ 0 , u t t + ∂ ν u − Δ Γ u + Q ( x , u t ) = g ( x , u ) on ( 0 , ∞ ) × Γ 1 , u ( 0 , x ) = u 0 ( x ) , u t ( 0 , x ) = u 1 ( x ) in Ω ‾ , where Ω is a bounded open subset of R N with C 1 boundary ( N ≥ 2 ), Γ = ∂ Ω , Γ 1 is relatively open on Γ, Δ Γ denotes the Laplace–Beltrami operator on Γ, ν is the outward normal to Ω, and the terms P and Q represent nonlinear damping terms, while f and g are nonlinear perturbations. In the paper we establish local and global existence, uniqueness and Hadamard well-posedness results when source terms can be supercritical or super-supercritical. [ABSTRACT FROM AUTHOR]

This paper investigates the use of non-linear structural dynamics computational models with multiple levels of fidelity for the calibration of input-dependent model parameters. Non-linear behavior complicates the calibration of model parameters that are input-dependent (i.e., functions of temperature, loading, etc.). Different types of models may also be available for the estimation of unmeasured system properties, with different levels of physics fidelity, mesh resolution and boundary condition assumptions. To infer these system properties, Bayesian calibration can be used to combine information from multiple sources (such as experimental data and prior knowledge) and comprehensively quantify the uncertainty in the model parameters. Estimating the posterior distributions is done using Markov Chain Monte Carlo sampling, which requires a large number of computations, thus making the use of a high-fidelity model for calibration prohibitively expensive. On the other hand, use of a low-fidelity model could lead to significant error in calibration and prediction. Therefore, this paper pursues a multi-fidelity approach for input-dependent model parameter calibration with a surrogate of the low-fidelity model corrected using higher fidelity simulation data. The effects of non-linear behavior and input-dependence of model parameters and measurements are systematically organized and distinguished using a Bayesian network. The methodology is illustrated with a curved panel, subjected to acoustic and thermal loading, where the damping properties of the panel could be functions of the acoustic loading and temperature magnitudes. Two models (a frequency response analysis and a full time history analysis) are combined to estimate the damping characteristics of the panel. The effect of the temperature on the performance of the strain gages and therefore on the calibration results is also studied. [ABSTRACT FROM AUTHOR]

The paper is concerned with the interconnection of the boundary behaviour of the solutions of the exterior Dirichlet, Neumann and Robin problems of harmonic analysis for the unit ball in R 3 with the corresponding behaviour of the associated ergodic inverse problems for the entire space. Rates of approximation play a basic role. The solutions themselves are evaluated by means of Fourier expansions with respect to spherical harmonics. In case of the first two problems, the basis for the investigation of the approximation and ergodic behaviour is the theory of semigroups of linear operators mapping a Banach space X into itself. The connection between the semigroup property and the major premise of Huygens' principle is emphasized. Another tool is a Drazin-like inverse operator A ad for the infinitesimal generator A of a semigroup that arises naturally in ergodic theory. This operator is a closed, not necessarily bounded, operator. It was introduced in a paper with U. Westphal (Butzer and Westphal, 1970/71) and extended to a generalized setting with J.J. Koliha (Butzer and Koliha, 2009). Unlike the latter two problems, the solution of Robin's problem does not have the semigroup property and therefore the semigroup methods applied to Dirichlet's and Neumann's problem do not work. The authors give several hints how to overcome these difficulties. [ABSTRACT FROM AUTHOR]

The main goal of this paper is to propose an approach to inverse spectral problems for functional-differential operators (FDO) with involution. For definiteness, we focus on the second-order FDO with involution-reflection. Our approach is based on the reduction of the problem to the matrix form and on the solution of the inverse problem for the matrix Sturm-Liouville operator by developing the method of spectral mappings. The obtained matrix Sturm-Liouville operator contains the weight, which causes qualitative difficulties in the study of the inverse problem. As a result, we show that the considered FDO with involution is uniquely specified by five spectra of certain regular boundary value problems. [ABSTRACT FROM AUTHOR]

In this paper we study the stability of the energy associated to an initial boundary value problem involving a semilinear wave equation. The domain is an unbounded subset of R 2 with unbounded boundary. On the function responsible by the semilinearity we consider a very general decay rate control. The main tools are the Trudinger–Moser inequality combined with microlocal analysis arguments and the techniques of Lasiecka and Tataru. [ABSTRACT FROM AUTHOR]

This paper deals with a two-species and a three-species predator–prey reaction diffusion systems where the diffusion coefficients are density dependent and time-delays are involved in the reaction functions. The diffusion terms are of porous medium type which are degenerate and the boundary conditions are of Neumann type. The aim of the paper is to investigate the asymptotic behavior of the time-dependent solution in relation to various steady-state solutions of the system. This includes the existence of a unique positive classical solution, global attraction of a steady-state solution, and stability or instability of various semitrivial solutions and the positive steady-state solution. [ABSTRACT FROM AUTHOR]

This paper deals with the chemotaxis–haptotaxis system with signal-dependent sensitivity { u t = Δ u − ∇ ⋅ ( χ ( v ) u ∇ v ) − ξ ∇ ⋅ ( u ∇ w ) + μ u ( 1 − u − w ) , x ∈ Ω , t > 0 , v t = Δ v − v + u , x ∈ Ω , t > 0 , w t = − v w , x ∈ Ω , t > 0 under homogeneous Neumann boundary conditions and initial conditions, where Ω ⊂ R n ( n ≥ 3 ) is a bounded domain with smooth boundary, ξ , μ > 0 are constants and χ is a function satisfying some conditions. In the case that χ is a constant it is known that the above system possesses a global classical solution under some conditions (Cao [4] , Tao [10] , Tao and Winkler [11] ); however, in the case that χ is a function, the above system has not been studied. The purpose of this paper is to establish global existence and boundedness in the above system. [ABSTRACT FROM AUTHOR]

In this paper we study a free boundary problem modeling the growth of vascularized tumors. The model is a modification to the Byrne–Chaplain tumor model that has been intensively studied during the past two decades. The modification is made by replacing the Dirichlet boundary value condition with the Robin condition, which causes some new difficulties in making rigorous analysis of the model, particularly on existence and uniqueness of a radial stationary solution. In this paper we successfully solve this problem. We prove that this free boundary problem has a unique radial stationary solution which is asymptotically stable for large surface tension coefficient, whereas unstable for small surface tension coefficient. Tools used in this analysis are the geometric theory of abstract parabolic differential equations in Banach spaces and spectral analysis of the linearized operator. [ABSTRACT FROM AUTHOR]

Abstract This paper addresses the question of the impact of the boundary on the dynamical behaviour of finite Boolean automata networks on Z2. The evolution over discrete time of such networks is governed by a specific stochastic threshold non-linear transition rule derived from the classical rule of formal neural networks. More precisely, the networks considered in this paper are finite but the study is done for arbitrarily large sizes. Moreover, the boundary impact is viewed as a classical definition of a phase transition in probability theory, characterising in our context the fact that a network admits distinct asymptotic behaviours when different boundary instances are assumed. The main contribution of this paper is the highlight of a formula for a necessary condition for boundary sensitivity, whose sufficiency and necessity are entirely proven with natural constraints on interaction potentials. [ABSTRACT FROM AUTHOR]

We consider a class of tumor growth models under the combined effects of density-dependent pressure and cell multiplication, with a free boundary model as its singular limit when the pressure-density relationship becomes highly nonlinear. In particular, the constitutive law connecting pressure p and density ρ is p ( ρ ) = m m − 1 ρ m − 1 , and when m ≫ 1 , the cell density ρ may evolve its support according to a pressure-driven geometric motion with sharp interface along its boundary. The nonlinearity and degeneracy in the diffusion bring great challenges in numerical simulations. Prior to the present paper, there is lack of standard mechanism to numerically capture the front propagation speed as m ≫ 1 . In this paper, we develop a numerical scheme based on a novel prediction-correction reformulation that can accurately approximate the front propagation even when the nonlinearity is extremely strong. We show that the semi-discrete scheme naturally connects to the free boundary limit equation as m → ∞ . With proper spatial discretization, the fully discrete scheme has improved stability, preserves positivity, and can be implemented without nonlinear solvers. Finally, extensive numerical examples in both one and two dimensions are provided to verify the claimed properties in various applications. [ABSTRACT FROM AUTHOR]

The finite-difference time-domain (FDTD) method has been extensively used to simulate varieties of electromagnetic interaction problems. However, because of its Courant–Friedrich–Levy (CFL) condition, the maximum time step size of this method is limited by the minimum size of cell used in the computational domain. So the FDTD method is inefficient to simulate the electromagnetic problems which have very fine structures. To deal with this problem, the Hybrid Implicit Explicit (HIE)-FDTD method is developed. The HIE-FDTD method uses the hybrid implicit explicit difference in the direction with fine structures to avoid the confinement of the fine spatial mesh on the time step size. So this method has much higher computational efficiency than the FDTD method, and is extremely useful for the problems which have fine structures in one direction. In this paper, the basic formulations, time stability condition and dispersion error of the HIE-FDTD method are presented. The implementations of several boundary conditions, including the connect boundary, absorbing boundary and periodic boundary are described, then some applications and important developments of this method are provided. The goal of this paper is to provide an historical overview and future prospects of the HIE-FDTD method. [ABSTRACT FROM AUTHOR]

It has been 40 years since Lawson and Osserman introduced the three minimal cones associated with Dirichlet problems in their 1977 Acta paper [13] . The first cone was shown area-minimizing by Harvey and Lawson in the celebrated paper [10] . In this paper, we confirm that the other two are also area-minimizing. In fact, we show that every Lawson–Osserman cone of type ( n , p , 2 ) constructed in [26] is area-minimizing. [ABSTRACT FROM AUTHOR]

In the paper a numerical approach to the four equation model of water hammer (WH) with fluid–structure interaction (FSI) is presented. An algorithm for numerical solution of that model in time domain based on the method of characteristics (MOC) is proposed. Special attention is paid to modeling of the influence of viscoelastic pipe supports. The boundary condition at the support is formulated as a differential equation of junction motion which is solved numerically concurrently with the MOC compatibility equations. Numerical simulations were done with the use of an own computer program. The tests were conducted for a model of real pipeline built in the lab and fastened with specific, complex, supporting system. The numerical results were compared to the experimental records and after some modeling and calibration quite a good agreement was achieved. Basic behaviors in the pressure records were identified and existed discrepancies were discussed and explained. The other physical model tested and preliminary analyzed in the paper is a straight pipeline fixed to the floor only with viscoelastic supports. Such a design created the possibility of significant displacements of the pipeline on the supports and thus allowing for effective testing of the influence of supports stiffness and damping properties on the WH behaviors. The specific initial conditions for that case were formulated, solved and implemented in a computer code. It was found that proper selection of support parameters may produce significant reduction of WH pressure amplitudes, mainly due to effective energy absorption and dissipation at the supports. [ABSTRACT FROM AUTHOR]

In this paper we investigate boundary blow-up solutions of the problem { − Δ p ( x ) u + f ( x , u ) = ± K ( x ) | ∇ u | m ( x ) in Ω , u ( x ) → + ∞ as d ( x , ∂ Ω ) → 0 , where Δ p ( x ) u = div ( | ∇ u | p ( x ) − 2 ∇ u ) is called the p ( x ) -Laplacian. Our results extend the previous work [25] of Y. Liang, Q.H. Zhang and C.S. Zhao from the radial case to the non-radial setting, and [46] due to Q.H. Zhang and D. Motreanu from the assumption that K ( x ) | ∇ u ( x ) | m ( x ) is a small perturbation, to the case in which ± K ( x ) | ∇ u | m ( x ) is a large perturbation. We provide an exact estimate of the pointwise different behavior of the solutions near the boundary in terms of d ( x , ∂ Ω ) and in terms of the growth of the exponents. Furthermore, the comparison principle is no longer applicable in our context, since f ( x , ⋅ ) is not assumed to be monotone in this paper. [ABSTRACT FROM AUTHOR]

This paper presents a robust numerical solution to the electromagnetic scattering problem involving multiple multi-layered cavities in both transverse magnetic and electric polarizations. A transparent boundary condition is introduced at the open aperture of the cavity to transform the problem from an unbounded domain into that of bounded cavities. By employing Fourier series expansion of the solution, we reduce the original boundary value problem to a two-point boundary value problem, represented as an ordinary differential equation for the Fourier coefficients. The analytical derivation of the connection formula for the solution enables us to construct a small-scale system that includes solely the Fourier coefficients on the aperture, streamlining the solving process. Furthermore, we propose accurate numerical quadrature formulas designed to efficiently handle the weakly singular integrals that arise in the transparent boundary conditions. To demonstrate the effectiveness and versatility of our proposed method, a series of numerical experiments are conducted. • Innovative Transparent Boundary Conditions: A novel transparent boundary condition is designed to enhance the accuracy and efficiency of simulations. • Enhanced Numerical Quadratures: Highly efficient and accurate numerical quadrature techniques are developed to ensure reliable computations. • Comprehensive Numerical Examples: An array of carefully chosen numerical examples are presented to demonstrate the method's effectiveness and versatility. [ABSTRACT FROM AUTHOR]

Suppose G is a compact Lie group, H is a closed subgroup of G , and the homogeneous space G / H is connected. The paper investigates the Ricci flow on a manifold M diffeomorphic to [ 0 , 1 ] × G / H . First, we prove a short-time existence and uniqueness theorem for a G -invariant solution g ( t ) satisfying the boundary condition II ( g ( t ) ) = F ( t , g ∂ M ( t ) ) and the initial condition g ( 0 ) = g ˆ . Here, II ( g ( t ) ) is the second fundamental form of ∂ M , g ∂ M is the metric induced on ∂ M by g ( t ) , F is a smooth map and g ˆ is a metric on M . Second, we study Perelman's F -functional on M . Our results show, roughly speaking, that F is non-decreasing on a G -invariant solution to the modified Ricci flow, provided that this solution satisfies boundary conditions inspired by a paper of Gianniotis. [ABSTRACT FROM AUTHOR]

In this paper, we study a free boundary model describing growth of tumors under action of drugs. To our knowledge, in theoretical discussion for free boundary problems, the proliferation rate in tumor models discussed in previous bifurcation results is a linear function of nutrients and inhibitors. Whereas in this paper we consider the net proliferation rate as a nonlinear function depending on both nutrients and drugs. First, the existence and the uniqueness of radially symmetric stationary solutions are obtained. Second, we prove that symmetry-breaking solutions bifurcate from the radially symmetric stationary solutions when the concentration of drug on the boundary of tumor is less than one in the rescaled model. [ABSTRACT FROM AUTHOR]