首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 922 毫秒
1.
We present a theoretical study of the resonances of a fluid–structureproblem, an elastic plate placed in a duct in the presence ofa compressible fluid. The case of a rigid plate has been largelystudied. Acoustic resonances are then associated to resonantmodes trapped by the plate. Due to the elasticity of the plate,we need to solve a quadratic eigenvalue problem in which theresonance frequencies k solve the equations (k) = k2, where are the eigenvalues of a self-adjoint operator of the formA + kB. First, we show how to study the eigenvalues locatedbelow the essential spectrum by using the min–max principle.Then, we study the fixed-point equations. We establish sufficientconditions on the characteristics of the plate and of the fluidto ensure the existence of resonances. Such conditions are validatednumerically.  相似文献   

2.
Bhupen Deka Department of Mathematics, Assam University, Silchar-788011, India A finite-element discretization, independent of the locationof the interface, is proposed and analysed for linear ellipticand parabolic interface problems. We establish error estimatesof optimal order in the H1-norm and almost optimal order inthe L2-norm for elliptic interface problems. An extension toparabolic interface problems is also discussed and an optimalerror estimate in the L2(0, T;H1())-norm and an almost optimalorder estimate in the L2(0, T;L2())-norm are derived for thespatially discrete scheme. A fully discrete scheme based onthe backward Euler method is analysed and an optimal order errorestimate in the L2(0, T;H1())-norm is derived. The interfacesare assumed to be of arbitrary shape and smooth for our purpose.  相似文献   

3.
Massimo Fornasier Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate, Università "La Sapienza" in Roma, Via Antonio Scarpa, 16/B, I-00161 Roma, Italy Rob Stevenson|| Department of Mathematics, Utrecht University, PO Box 80.010, NL-3508 TA Utrecht, The Netherlands This paper is concerned with the development of adaptive numericalmethods for elliptic operator equations. We are particularlyinterested in discretization schemes based on wavelet frames.We show that by using three basic subroutines an implementable,convergent scheme can be derived, which, moreover, has optimalcomputational complexity. The scheme is based on adaptive steepestdescent iterations. We illustrate our findings by numericalresults for the computation of solutions of the Poisson equationwith limited Sobolev smoothness on intervals in 1D and L-shapeddomains in 2D.  相似文献   

4.
A new solution is obtained for Lamé's equation in the form of a perturbation seriesabout k = 0; numerical accuracy appears to be high when v ismoderate and k not too close to 1, so that the solution is ofpractical value in problems involving elliptic cones or infinitesectors.  相似文献   

5.
An Rm-valued sequence (xk): = (xk : k = 1, 2, ...), e.g. generatedrecursively by xk = fk (xkk, Uk), is called ‘averagepth power bounded’ if (1/K) is bounded uniformly in K= 1, 2,.... (The case p = 2 may correspond to ‘power’in the physical sense.) This is a notion of stability. Givenestimates of the form: fk (x, u) < a x + ¶ k conditionsare obtained on the coefficient sequence (ak) and the inputestimates ek:=¶k (uk) which ensure this form of stabilityfor the output (xk). In particular, a condition (utilized inan application to adaptive control) is obtained which imposes(i) a bound b on (ak) and a ‘sparsity measure’ m(K) on #{kK: ak>} as K ( >1) (ii) average pth power boundednesson (ek), and (iii) a growth condition on (ek) related to b andm (•). This condition is sharp.  相似文献   

6.
Finite element approximations for the Dirichlet problem associated to a second-order elliptic differential equation are studied. The purpose of this paper is to discuss domain embedding preconditioners for discrete systems. The essential boundary condition on the interior interface is removed by introducing Lagrange multipliers. The associated discrete system, with a saddle point structure, is preconditioned by a block diagonal preconditioner. The main contribution of this paper is to propose a new operator, constructed from the -inner product, for the block of the preconditioner corresponding to the multipliers.

  相似文献   


7.
We observe that the Bargmann transform arises from the polardecomposition of a simple restriction operator. By virtue ofthe Bargmann transform, the Wigner transform is imbedded inthe geometric symmetry group of phase space. From this imbedding,we easily obtain that the spectrum of the Wigner transform consistsof {eik/12:k=0, 1, ..., 23}. 1991 Mathematics Subject Classification81R30, 47B38, 81S30.  相似文献   

8.
Sorin Micu This paper studies the numerical approximation of the boundarycontrol for the wave equation in a square domain. It is knownthat the discrete and semi-discrete models obtained by discretizingthe wave equation with the usual finite-difference or finite-elementmethods do not provide convergent sequences of approximationsto the boundary control of the continuous wave equation as themesh size goes to zero. Here, we introduce and analyse a newsemi-discrete model based on the space discretization of thewave equation using a mixed finite-element method with two differentbasis functions for the position and velocity. The main theoreticalresult is a uniform observability inequality which allows usto construct a sequence of approximations converging to theminimal L2-norm control of the continuous wave equation. Wealso introduce a fully discrete system, obtained from our semi-discretescheme, for which we conjecture that it provides a convergentsequence of discrete approximations as both h and t, the timediscretization parameter, go to zero. We illustrate this factwith several numerical experiments.  相似文献   

9.
Let be a discrete group acting on a compact manifold X, letV be a -equivalent Hermitian vector bundle over X, and let Dbe a first-order elliptic self-adjoint -equivalent differentialoperator acting on sections of V. This data is used to defineToeplitz operators with symbols in the transformation groupC*-algebra C(X), and it is shown that if the symbol of sucha Toeplitz operator is invertible, then the operator is Fredholm.In the case where is finite and acts freely on X, a geometric-topologicalformula for the index is stated that involves an explicitlyconstructed differential form associated to the symbol. 2000Mathematics Subject Classification 47A53 (primary), 19K56, 47B35,46L87 (secondary).  相似文献   

10.
Let A be a unital maximal full algebra of operator fields withbase space [0, 1]k and fibre algebras . It is shown in this paper that the stable rankof A is bounded above by the quantity , where ‘sr’ means stable rank. Usingthe above estimate, the stable ranks of the C*-algebras of the(possibly higher rank) discrete Heisenberg groups are computed.2000 Mathematics Subject Classification 47L99.  相似文献   

11.
This paper makes contributions to the general theory of wavepropagation in conservative systems under conditions when theproportional change in amplitude or wavenumber over a distanceof one wavelength is very small. For linear systems, such propagationis governed by the well-known theory of group velocity; thereis "frequency dispersion", in the sense that energy in componentsof different frequency is propagated at different group velocities.For non-linear systems without frequency dispersion, e.g. acousticsystems, a different, but also well-known, modification of thewaveform occurs. It may be called "amplitude dispersion", inthat different values of an amplitude variable like the pressureare propagated at different speeds. A much more general theory of non-linear systems, where frequencydispersion and amplitude dispersion would be expected to bein competition, has been given by Whitham (1965b). Energy doesnot play a key role in the theory, because it is easily transferredbetween components of different frequencies. The fundamentalequation follows from Hamilton's principle in an averaged form. In examples given by Whitham, changes in, for example, wavenumber(or amplitude) are propagated at two different velocities, becausethe fundamental equation is hyperbolic. However, in the limitingcase of infinitesimal amplitude, the equation is parabolic andonly one velocity of propagation (the group velocity) occurs.Thus, Whitham showed that non-linearity can "split" the groupvelocity. This paper is concerned with the inference of detailed conclusionsfrom Whitham's theory, to enable comparisons with experimentthat will show the range of applicability of the theory. Itattempts to obtain these in the simplest case, namely, thatof one-dimensional propagation when Whitham's "pseudo-frequencies"are absent. If the relationship between frequency and wavenumber k forinfinitesimal amplitude is = f(k), then for finite amplitudethe equation is shown to be hyperbolic or elliptic respectively,according as [–f(k)]f*(k) takes positive or negative values.For gravity waves on deep water this product is negative andthese, it is inferred, may be good for comparison of theorywith experiment in the elliptic case. A new non-linear non-perturbationaltheory of waves under the combined action of gravity and surfacetension is used to indicate that waves at 9.6 c/s on mercurymay be suitable for comparison with experiment in the hyperboliccase. When non-linear effects are only moderate, approximate transformationsof Whitham's equation to the axisymmetric Laplace and wave equationsrespectively, in the elliptic and hyperbolic cases, are usedto obtain particular solutions for comparison with experiment.A feature of these solutions is the appearance of discontinuitiesin wavelength. For example, when a wavemaker creates gravity waves of fixedfrequency whose amplitude first increases and then decreases,the theory predicts that the length of waves in the group decreasesahead of the point of maximum amplitude and increases behindit. This produces in turn a concentration of energy towardsthe centre of the group, which continues during the whole periodbefore a discontinuity in wavelength actually forms. This solutionin the elliptic case is obtained with the aid of the theoryof imaginary characteristics.  相似文献   

12.
Charalambos Makridakis In this paper, we derive a posteriori error estimates for space-discreteapproximations of the time-dependent Stokes equations. By usingan appropriate Stokes reconstruction operator, we are able towrite an auxiliary error equation, in pointwise form, that satisfiesthe exact divergence-free condition. Thus, standard energy estimatesfrom partial differential equation theory can be applied directly,and yield a posteriori estimates that rely on available correspondingestimates for the stationary Stokes equation. Estimates of optimalorder in L(L2) and L(H1) for the velocity are derived for finite-elementand finite-volume approximations.  相似文献   

13.
Wolfgang Hackbusch In this paper, we discuss the application of hierarchical matrixtechniques to the solution of Helmholtz problems with largewave number in 2D. We consider the Brakhage–Werner integralformulation of the problem discretized by the Galerkin boundary-elementmethod. The dense n x n Galerkin matrix arising from this approachis represented by a sum of an -matrix and an 2-matrix, two different hierarchical matrix formats.A well-known multipole expansion is used to construct the 2-matrix. We present a new approach to dealingwith the numerical instability problems of this expansion: theparts of the matrix that can cause problems are approximatedin a stable way by an -matrix. Algebraic recompression methods are used to reducethe storage and the complexity of arithmetical operations ofthe -matrix.Further, an approximate LU decomposition of such a recompressed-matrix is aneffective preconditioner. We prove that the construction ofthe matrices as well as the matrix-vector product can be performedin almost linear time in the number of unknowns. Numerical experimentsfor scattering problems in 2D are presented, where the linearsystems are solved by a preconditioned iterative method.  相似文献   

14.
A metric space X has the unique midset property if there isa topology-preserving metric d on X such that for every pairof distinct points x, y there is one and only one point p suchthat d(x, p) = d(y, p). The following are proved. (1) The discretespace with cardinality n has the unique midset property if andonly if n 2, 4 and n c, where c is the cardinality of thecontinuum. (2) If D is a discrete space with cardinality notgreater than c, then the countable power DN of D has the uniquemidset property. In particular, the Cantor set and the spaceof irrational numbers have the unique midset property. A finite discrete space with n points has the unique midsetproperty if and only if there is an edge colouring of the completegraph Kn such that for every pair of distinct vertices x, ythere is one and only one vertex p such that (xp) = (yp). Letump(Kn) be the smallest number of colours necessary for sucha colouring of Kn. The following are proved. (3) For each k 0, ump(K2k+1) = k. (4) For each k 3, k ump(K2k) 2k–1.  相似文献   

15.
W. Kotarski Institute of Informatics, Silesian University, Bedzinska 60, 41-200 Sosnowiec, Poland Email: bahaa_gm{at}hotmail.com Email: kotarski{at}gate.math.us.edu.pl Received on March 14, 2006; Accepted on December 20, 2006 A distributed control problem for n x n parabolic coupled systemsinvolving operators with infinite order is considered. The performanceindex is more general than the quadratic one and has an integralform. Constraints on controls are imposed. Making use of theDubovitskii–Milyutin theorem, the necessary and sufficientconditions of optimality are derived for the Dirichlet problem.Yet, the problem considered here is more general than the problemsin El-Saify & Bahaa (2002, Optimal control for n x n hyperbolicsystems involving operators of infinite order. Math. Slovaca,52, 409–424), El-Zahaby (2002, Optimal control of systemsgoverned by infinite order operators. Proceeding (Abstracts)of the International Conference of Mathematics (Trends and Developments)of the Egyptian Mathematical Society, Cairo, Egypt, 28–31December 2002. J. Egypt. Math. Soc. (submitted)), Gali &El-Saify (1983, Control of system governed by infinite orderequation of hyperbolic type. Proceeding of the InternationalConference on Functional-Differential Systems and Related Topics,vol. III. Poland, pp. 99–103), Gali et al. (1983, Distributedcontrol of a system governed by Dirichlet and Neumann problemsfor elliptic equations of infinite order. Proceeding of theInternational Conference on Functional-Differential Systemsand Related Topics, vol. III. Poland, pp. 83–87) and Kotarskiet al. (200b, Optimal control problem for a hyperbolic systemwith mixed control-state constraints involving operator of infiniteorder. Int. J. Pure Appl. Math., 1, 241–254).  相似文献   

16.
Mingzhu Liu Department of Mathematics, Harbin Institute of Technology, Harbin 150001, China Email: ghu{at}hit.edu.cn, ghuca{at}yahoo.ca Email: mzliu{at}hit.edu.cn Received on January 31, 2006; Accepted on December 9, 2006 In this paper, we are concerned with the properties of the weightedlogarithmic matrix norms. A relation between the elliptic logarithmicmatrix norm and the weighted logarithmic matrix norm is given.Based on Lyapunov equations, two weighted logarithmic matrixnorms are constructed which are less than 1-logarithmic matrixnorm and -logarithmic matrix norm, respectively. Then, an iterativescheme is presented to obtain the logarithmically -efficientmatrix norm. Numerical examples are given to illustrate theresults.  相似文献   

17.
We prove that, if 2 k1 k2, then there is no infinite sequence of positive integers such that the representation functionr(n) = #{(a, a'): n = k1a + k2a', a, a' } is constant for nlarge enough. This result completes the previous work of Diracand Moser for the special case k1 = 1 and answers a questionposed by Sárkozy and Sós.  相似文献   

18.
Recall that an infinite Toeplitz matrix, , say, is one for which the values tjk dependonly on jk, so we may write tjk = ajk; thusthe entries are constant on diagonals sloping down from topleft to bottom right:

Sucha matrix is said to be banded if there are only finite manynon-zero aj, say . The propertiesof such a matrix, when defining a linear operator T on a sequencespace, are often closely related to the properties of the Laurentpolynomial

for zlying on the unit circle in the complex plane. For example,suppose that the matrix acts on the most familiar sequence space,the Hilbert space 2; then it has been known for many years thatthe operator norm of T is the maximum  相似文献   

19.
Using Z/k-manifolds we give a geometric interpretation of thering homomorphism from cobordismwith Z/k coefficients to elliptic cohomology with Z/k coefficients,induced by the map of spectra MSO Ell associated to the naturaltransformation of cohomology theories MSO* Ell*.  相似文献   

20.
In this paper we examine computation of optimal control u* ofthe exact controllability problem (referred to as the constraintproblem) governed by the following type of linear parabolicdifferential equations: (y/t) + Ay = u in Q y = 0 on y(0) = y0 on where A is the second-order elliptic differential operator, is a bounded domain in k with smooth boundary , Q = (0, T)x , = (0, T) x and T > 0. This is achieved by approximatingu* through a sequence {un} of controls corresponding to unconstrainedproblems involving a penalty function arising from the controllabilityconstraint.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号