Search

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 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]

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]

Various techniques are developed for addressing the existence, uniqueness and numerical calculation of Wiener path integral (WPI) most probable path solutions. Specifically, the WPI technique for determining the stochastic response of diverse nonlinear dynamical systems treats the system response joint transition probability density function as a functional integral over the space of all possible paths connecting the initial and the final states of the response vector. This functional integral is evaluated, ordinarily, by resorting to an approximate approach that considers the contribution only of the most probable path. The most probable path corresponds to an extremum of the functional integrand and is determined by solving a functional minimization problem that takes the form of a deterministic boundary value problem (BVP). In this paper, first, it is shown that for the commonly considered case of the system nonlinearity being of polynomial form, there exist globally optimal solutions corresponding to the most probable path BVP. Further, relying on algebraic geometry concepts and tools, a condition is derived for determining if the BVP for the most probable path exhibits a unique solution over a specific region. Furthermore, a novel solution approach is developed for the BVP by relying on Sylvester's dialytic method of elimination. Notably, the method reduces the complexity of the BVP system of coupled multivariate polynomial equations by eliminating one or more variables. Various numerical examples pertaining to diverse nonlinear oscillators are included for demonstrating the capabilities of the developed techniques. [ABSTRACT FROM AUTHOR]