Quasi-separation of the biharmonic partial differential equation

In this paper, we consider analytical and numerical solutions to the Dirichlet boundary-value problem for the biharmonic partial differential equation on a disc of finite radius in the plane. The physical interpretation of these solutions is that of the harmonic oscillations of a thin, clamped plate. For the linear, fourth-order, biharmonic partial differential equation in the plane, it is well known that the solution method of separation in polar coordinates is not possible, in general. However, in this paper, for circular domains in the plane, it is shown that a method, here called quasi-separation of variables, does lead to solutions of the partial differential equation. These solutions are products of solutions of two ordinary linear differential equations: a fourth-order radial equation and a second-order angular differential equation. To be expected, without complete separation of the polar variables, there is some restriction on the range of these solutions in comparison with the corresponding separated solutions of the second-order harmonic differential equation in the plane. Notwithstanding these restrictions, the quasi-separation method leads to solutions of the Dirichlet boundary-value problem on a disc with centre at the origin, with boundary conditions determined by the solution and its inward drawn normal taking the value 0 on the edge of the disc. One significant feature for these biharmonic boundary-value problems, in general, follows from the form of the biharmonic differential expression when represented in polar coordinates. In this form, the differential expression has a singularity at the origin, in the radial variable. This singularity translates to a singularity at the origin of the fourth-order radial separated equation; this singularity necessitates the application of a third boundary condition in order to determine a self-adjoint solution to the Dirichlet boundary-value problem. The penultimate section of the paper reports on numerical solutions to the Dirichlet boundary-value problem; these results are also presented graphically. Two specific cases are studied in detail and numerical values of the eigenvalues are compared with the results obtained in earlier studies.

The General Differential Operators Generated by a Quasi-Differential Expressions with their Interior Singular Points

The general ordinary quasi-differential expression M of n-th order with complex coefficients and its formal adjoint M + are considered over a regoin (a, b) on the real line, −∞ ≤ a < b ≤ ∞, on which the operator may have a finite number of singular points. By considering M over various subintervals on which singularities occur only at the ends, restrictions of the maximal operator generated by M in L2|w (a, b) which are regularly solvable with respect to the minimal operators T0 (M ) and T0 (M + ). In addition to direct sums of regularly solvable operators defined on the separate subintervals, there are other regularly solvable restrications of the maximal operator which involve linking the various intervals together in interface like style.

On the Stabilization of the Wave Equation by the Boundary

* Partially supported by CNPq (Brazil)

Asymptotic Property of Eigenvalues and Eigenfunctions of the Laplace Operator in Domain with a Perturbed Boundary

2000 Mathematics Subject Classification: 35J05, 35C15, 44P05

Anti-Periodic Boundary Value Problem for Impulsive Fractional Integro Differential Equations

MSC 2010: 34A37, 34B15, 26A33, 34C25, 34K37

Crack Theory with Singularities at the Boundary

2002 Mathematics Subject Classification: 35S15, 35J70, 35J40, 38J40

Leak identification in saturated unsteady flow via a Cauchy problem

This work is an initial study of a numerical method for identifying multiple leak zones in saturated unsteady flow. Using the conventional saturated groundwater flow equation, the leak identification problem is modelled as a Cauchy problem for the heat equation and the aim is to find the regions on the boundary of the solution domain where the solution vanishes, since leak zones correspond to null pressure values. This problem is ill-posed and to reconstruct the solution in a stable way, we therefore modify and employ an iterative regularizing method proposed in [1] and [2]. In this method, mixed well-posed problems obtained by changing the boundary conditions are solved for the heat operator as well as for its adjoint, to get a sequence of approximations to the original Cauchy problem. The mixed problems are solved using a Finite element method (FEM), and the numerical results indicate that the leak zones can be identified with the proposed method.

Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.

An Optimal Order Interior Penalty Discontinuous Galerkin Discretization of the Compressible Navier-Stokes Equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior-penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L_2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi & Rebay, cf. [11], the standard SIPG method outlined in [25], and an NIPG variant of the new scheme will be undertaken.

Unfolding Operator Method for Thin Domains with a Locally Periodic Highly Oscillatory Boundary

We analyze the behavior of solutions of the Poisson equation with homogeneous Neumann boundary conditions in a two-dimensional thin domain which presents locally periodic oscillations at the boundary. The oscillations are such that both the amplitude and period of the oscillations may vary in space. We obtain the homogenized limit problem and a corrector result by extending the unfolding operator method to the case of locally periodic media.

Homogenization of the p-Laplacian with nonlinear boundary condition on critical size particles: Identifying the strange term for the some non smooth and multivalued operators

We extend previous papers in the literature concerning the homogenization of Robin type boundary conditions for quasilinear equations, in the case of microscopic obstacles of critical size: here we consider nonlinear boundary conditions involving some maximal monotone graphs which may correspond to discontinuous or non-Lipschitz functions arising in some catalysis problems.

Second-order elliptic PDE with discontinuous boundary data

We shall consider the weak formulation of a linear elliptic model problem with discontinuous Dirichlet boundary conditions. Since such problems are typically not well-defined in the standard H^1-H^1 setting, we will introduce a suitable saddle point formulation in terms of weighted Sobolev spaces. Furthermore, we will discuss the numerical solution of such problems. Specifically, we employ an hp-discontinuous Galerkin method and derive an L^2-norm a posteriori error estimate. Numerical experiments demonstrate the effectiveness of the proposed error indicator in both the h- and hp-version setting. Indeed, in the latter case exponential convergence of the error is attained as the mesh is adaptively refined.

Migração lateral da desembocadura do Rio Itapocú, SC, Brasil: evolução morfológica e condicionantes físicas

Desembocaduras são ambientes bastante dinâmicos e sujeitos à complexa interação entre fatores estabilizadores e desestabilizadores. Dependendo dessa interação, desembocaduras podem apresentar a tendência de migração ao longo de barreiras arenosas. Um dos mecanismos mais eficientes de transporte de sedimento paralelo à costa, e consequentemente migração de canais, são as correntes longitudinais geradas pelas ondas se aproximando obliquamente à costa. A motivação do presente trabalho é entender o comportamento morfodinâmico do sistema de desembocadura do rio Itapocú, localizado no centro-norte de Santa Catarina (SC), frente aos processos forçantes que atuam na sua migração ao longo da linha de costa. A morfologia dos pontais arenosos foi obtida a partir de levantamentos morfológicos com o uso de DGPS. Para analisar a refração de ondas foi utilizado o modelo numérico MIKE 21 SW, sendo considerados como condições de contorno os dados de ondas referentes ao ano de 2002 e os dados de ondas previstos referentes ao período de coleta. Os dados de saída do modelo foram utilizados para estimar a deriva litorânea potencial na região. Os resultados morfológicos obtidos demonstraram uma migração da desembocadura para o norte durante o período analisado, sendo mais intenso durante o inverno e o verão. Ondas incidentes do quadrante sul sofreram mais o fenômeno da refração e as ondas de leste apresentaram menor variação angular ao se aproximarem à costa. A deriva litorânea potencial anual para os dados de ondas de 2002 apresentou sentido norte-sul, com inversão de sentido durante o outono. Utilizando os dados de ondas previstas para o período dos levantamentos, a deriva litorânea potencial estimada apresentou sentido sul-norte, concordando com a migração observada. Na região próxima a desembocadura, nos pontais arenosos, a deriva potencial apresentou direção para o norte durante todas as estações. Os dados de descarga fluvial não apresentaram influência na migração do canal, porém apresentaram uma relação com a largura do mesmo sazonalmente.Os dados de morfologia juntamente com os dados de deriva litorânea referentes às ondas de 2004/2005 mostraram claramente a migração do canal para o norte sendo a deriva a principal contribuinte para a migração da desembocadura.

Brazilian offshore wave climate based on NWW3 reanalysis

This paper provides a description of the wave climate off the Brazilian coast based on an eleven-year time series (Jan/1997-Dec/2007) obtained from the NWW3 operational model hindcast reanalysis. Information about wave climate in Brazilian waters is very scarce and mainly based on occasional short-term observations, the present analysis being the first covering such temporal and spatial scales. To define the wave climate, six sectors were defined and analyzed along the Brazilian shelf-break: South (W1), Southeast (W2), Central (W3), East (W4), Northeast (W5) and North (W6). W1, W2 and W3 wave regimes are determined by the South Atlantic High (SAH) and the passage of synoptic cold fronts; W4, W5 and W6 are controlled by the Intertropical Convergence Zone (ITCZ) and its meridional oscillation. The most energetic waves are from the S, generated by the strong winds associated to the passage of cold fronts, which mainly affect the southern region. Wave power presents a decrease in energy levels from south to north, with its annual variation showing that the winter months are the most energetic in W1 to W4, while in W5 and W6 the most energetic conditions occur during the austral summer. The information presented here provides boundary conditions for studies related to coastal processes, fundamental for a better understanding of the Brazilian coastal zone.

PARTITIONED ANALYSIS FOR DIMENSIONALLY-HETEROGENEOUS HYDRAULIC NETWORKS

In this work an iterative strategy is developed to tackle the problem of coupling dimensionally-heterogeneous models in the context of fluid mechanics. The procedure proposed here makes use of a reinterpretation of the original problem as a nonlinear interface problem for which classical nonlinear solvers can be applied. Strong coupling of the partitions is achieved while dealing with different codes for each partition, each code in black-box mode. The main application for which this procedure is envisaged arises when modeling hydraulic networks in which complex and simple subsystems are treated using detailed and simplified models, correspondingly. The potentialities and the performance of the strategy are assessed through several examples involving transient flows and complex network configurations.