Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations
|
|
- Reynard Hunter
- 6 years ago
- Views:
Transcription
1 Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We numerically simulate the following two nonlinear evolution equations with a fourth order (biharmonic) leading term: 2 u ɛ 2 ( u 2 )u = u t in R 2 or R 3 and 2 u + ɛ 2 (( u 2 ) u) = u t in R 2 or R 3 with an initial value and a Dirichlet boundary conditions. We use a bivariate spline space like finite element method to solve these equations. We discuss the convergence of our numerical scheme and present several numerical experiments under different boundary conditions and different domains in the bivariate setting. AMS 2 Mathematics Subject Classifications: 35J35, 35K55, 65M2, 65M6 Keywords and Phrases: Nonlinear Biharmonic evolution equations, Ginzburg-Landau type equations, Bivariate Spline Method. Introduction In this paper, we are interested in the numerical behavior of the following two biharmonic equations 2 u ɛ 2 ( u 2 )u = u t, in R 2 or R 3 (.) u(x, ) = u (x), x u(x, t) = u (x) x u n (x, t) = u 2(x) x Department of Mathematics, University of Georgia, Athens, GA 362. This author is supported by the National Science Foundation under Grant #DMS His address is mjlai@math.uga.edu. Department of Mathematics, Pennsylvania State University, University Park, PA 682. His address is liu@math.psu.edu. Department of Mathematics, University of Georgia, Athens, GA 362. His address is paul@math.uga.edu.
2 and (.2) 2 u + ɛ 2 (( u 2 ) u) = u t, in R 2 or R 3 u(x, ) = u (x), x u(x, t) = u (x) x u n (x, t) = u 2(x) x, where is a bounded domain in R 2 or R 3 and n is the outward normal of the boundary. Equation (.) may be regarded as the gradient flow of the following energy functional: (.3) E(u) := 2 in the class ( u 2 + 2ɛ 2 ( u 2 ) 2 (.4) C = {u H 2 (), u = u, ) dx n u = u 2 }. Similarly, equation (.2) is the gradient flow of the energy functional (.5) E(u) := ( u 2 + ) 2 2ɛ 2 ( u 2 ) 2 dx in the class (.4). The understanding of the above two equations were originally motivated by our earlier study in the smectic liquid crystals (cf. [LLW3]). In that paper, the existence, uniqueness and stability of the solutions of both equations were established. In particular, for both equations, one of the important issues is to study the asymptotic behavior as ɛ. However, little is known about the variational problem involving the second order derivative term as those in (.3) and (.5). From the analysis point of view, one main difficulty is the lack of the maximum principle. In [LLW3], we used a special bivariate spline space in our proofs of the global existence of the weak solutions. We recall the weak formulation for (.): find u L 2 (, T, H 2 ()) such that u t vdx+ u vdx + ɛ (.6) 2 (u 2 )uvdx = f(u)vdx for all v H 2 (). The weak formulation for (.2) is as follows: Find u L 2 (, T, H 2 ()) such that [ u t vdx+ u v + ] ε (.7) 2 ( u 2 ) u v dx = f(u)vdx 2
3 for all v H 2 (). Here, we have already substracted a biharmonic function with the original boundary conditions in (.) and (.2). See [LLW3] for a detailed treatment and the compatiblity of the boundary conditions. We used the bivariate spline space S3r( +) r as finite dimensional approximation space in our Galerkin s procedure. The spline space mentioned above is defined as S r 3r( +) = {S C r () : s t P 3r, t +}, where P 3r is the space of all polynomials of degree 3r and + is a triangulated quadrangulation of the domain R 2 which is obtained from a nondegenerate quadrangulation of by adding two diagonals of each quadrilateral in (cf. [Lai and Schumaker99]). Such a special triangulation is general enough to partition any polygonal domains. We will show the convergence of the spline solution to the weak solution of the equations (.) and (.2) in 2 and 3, respectively. We explore the asymptotic behavior of the numerical solutions as time t and ɛ. We encountered many interesting phenomena indicating the possibility of the nonuniqueness of the solutions to the steady state equations (.) and (.2) as ɛ very small. That is, there could be many critical points for (.3) and (.5) or there are many paths leading to the minimizer. 2. Numerical Analysis of Equation (.) In this section, we will present numerical analysis and simulation of the spline solutions to both the steady state equation (2.) below and the time dependent equation (.). We will Newton s method to solve the nonlinear equations (2.) and Crank-Nicolson s method in time steps for (.). Find u V H 2 () such that (2.) a 2 (u, v) + b(u, u, v) f(u), v g, v = for all v V H 2 (), where and a 2 (φ, ψ) = b(θ, φ, ψ) = ɛ 2 φ ψdx ( θ 2 )φψdx. The stationary problem of equation (.) can be converted into (2.) for a certain f(u) (cf. [LLW3]). Let V = S r 3r( +) H 2 () be a spline subspace in H 2 (). We may absorb g, v into f(u), v and follow the same arguments of Theorem 2.3 in [LLW3] to have Theorem 2.. For any ε >, there exists a weak solution S u V satisfying (2.). Alternatively we may use the Brouwer fixed point theorem to prove Theorem 2.. (cf. [Lai and Wenston 98] for such kind of proof for the nonlinear biharmonic equation associated with the Navier-Stokes equations.) Note that the proof of Lemma 2.5 in [LLW3] can be used to prove the boundedness of S u. That is, 3
4 Lemma 2.. Let S u be a weak solution satisfying (2.) with V = S3r( +) r H 2 (). Then S u 2 dx C(ɛ). Also, the proof of Theorem 2.4 in [LLW3] yields the following Theorem 2.2. If ɛ is not very small, then the weak solution S u V is unique. We now see how close S u is to the weak solution u. Consider ɛ 2 a 2 (u S u, u S u ) = ɛ 2 a 2 (u S u, u S) + ɛ 2 a 2 (u S u, S S u ) where S is a spline function in S r 3r( +) H 2 (). Note that we have a 2 (u S u, S S u ) =b(s u, S u, S S u ) b(u, u, S S u ) + f(s u ) f(u), S S u =b(s u, S u, u S u ) + b(s u, S u, S u) b(u, u, u S u ) + b(u, u, u S) 3 ɛ 2 φ(s u + u)(s u u), S S u 3 ɛ 2 φ(s u u), S S u For b(s u, S u, S u) + b(u, u, u S), we have, by using a well-known inequality f 4 dx 2 f 2 dx f 2 dx and by Poincaré s inequality, b(s u, S u, u S) b(u, u, u S) b(s u, u S u, u s) + b(u, u, u S) b(s u, u, u S) ɛ 2 C ( S u 2, + ) S u u, u S, + ɛ 2 C 2 S u u, ( S u, + u, ) u, u S, C(ɛ, K, ) S u u 2, u S,, where k, is the usual semi norm of H k () and C(ɛ, K, ) is a constant. Here, we have used Lemma 2. and Theorem 2.5 in [LLW3]. For b(s u, S u, u S u ) b(u, u, u S u ), we have ɛ 2 (b(s u, S u, u S u ) b(u, u, u S u )) = S u u 4 dx 3 u 2 S u 2 dx + 4 S u u 2 dx.
5 For ɛ 2 φ(s u + u)(s u u), S S u + ɛ 2 φ(s u u), S S u, we have, φ(s u + u)(s u u), S S u + φ(s u u), S S u φ(s u + u)(s u u), u S u + φ(s u u), u S u + φ(s u + u)(s u u), S u + φ(s u u), S u 2Km( S u, + u ) u S u 2, + 2mK 2 u S u 2, + 2K 2 m( S u, + u, ) u S u, u S, + 2K 2 m u S u, u S, C 3 (m, K, ɛ, ) u S u 2 2, + C 3 (m, K, ɛ, ) u S u 2, u S 2, for a constant C 3 which is dependent on m, K, ɛ and, where K is the Poincaré constant. We now summarize the above discussion to have ɛ 2 a 2 (u S u, u S u ) + S u u 4 dx + 3 u 2 S 4 4 u 2 2 dx = ɛ 2 a 2 (u S u, u S) + S u u 2 dx Thus, we have + ɛ 2 b(s u, S u, S u) + ɛ 2 b(u, u, u S) 3 φ(s 2 u u 2 ), S S u 3 φ(s u u), S S u ɛ 2 u S u 2, u S 2, + K 4 u S u 2 2, + C(ɛ, K, ) u S u 2, u S 2, + C 3 (m, K, ɛ ) u S u 2, + C 3 (m, K, ɛ ) u S u 2, u S 2,. (ɛ 2 K 4 C 3 (m, ɛ, K, )) u S u 2, (ɛ 2 + C(ɛ, K, ) + C 3 (m, K, ɛ, )) u S 2,. That is, u S u 2, α u S 2, for any S S3( ) H 2 (), where constant α = (ɛ 2 + C(ɛ, K, )+C 3 (m, K, ɛ ))/(ɛ 2 K 4 C(m, ɛ, K, )). By using the spline approximation property (c.f. Theorem. in [Lai and Schumaker 99]), we complete the proof of the following Theorem 2.3. Suppose that ɛ is not very small. Suppose that φ H k () with k 3. Then the bivariate spline solution S u approximates the weak solution u and satisfies S u u 2, C + k 2 φ k,. Next we discuss the implementation of bivariate spline spaces to numerically solve both the equation (.) and its steady state equation (2.). We apply bivariate spline space 5
6 S 3( +) to discretize the space variables of the equations (.) using Galerkins s method and use Crank-Nicolson s method to discretize the time variable of the equations (.). For each step we employ the Newton method to iterate the nonlinear term in (.). We set the error between iterations to be 8 for a stopping criterion. For the time dependent problem (.), we use the solution at t = t i to be an initial guess for the solution for t = t i+ and then apply the Newton method. More precisely, our numerical scheme for (2.) is ) Start with S u, S 3( +) which is the numerical solution of a biharmonic equation: S u, vdx = gvdx, u = u (x), n u = u 2 (x) for all v S 3( +) H 2 () 2) For k =, 2, 3,, use Newton s method to solve S u,k S 3( +) satisfying 3) Letting S u,k = i= S u,k vdx + f(s u,k)vdx = gvdx u = u (x), n u = u 2 (x). N c i,k φ i, if {c i,k } {c i,k } l < 8, we stop the Newton iteration in 2) and use S u,k as the spline approximation of the weak solution of (2.). Finally, we present some numerical experiments with the equations (.) and (2.). We did the following two sets of numerical experiments. First, we test if our programs work correctly by inputing the programs the right-hand side and boundary conditions computed from a known test function u, computing a numerical approximation S u, and comparing the maximum errors between u and S u for several refinements of a underlying quadrangulation of a given domain. We did such a test for several functions over different domains. We tabulate some experimental results below in Tables In the following Tables , we numerically solve equation (2.) with known exact solutions u(x, y) = x 4 + y 4 and u(x, y) = sin(xy). We list the maximum errors of the numerical solutions against the exact solutions evaluated over points equally-spaced over [, ] [, ] for different values ɛ. Table 2.. The maximum errors of numerical solutions of (2.) against the exact solution u(x, y) = x 4 + y 4. matrix size ɛ = ɛ =. ɛ =. ɛ =. ɛ =
7 Table 2.2. The maximum errors of numerical solutions of (2.) against the exact solution u(x, y) = sin(xy). ²= ² =. ² =. ² =. ² = When ² =.5, our numerical experiments confirm the nonuniqueness of the solution of (2.). That is, starting from the spline solution of the biharmonic equation, the iterations converge. However, they converge to a different solution than the exact solution. That is why the errors in Tables 4. 2 do not decrease as the refinement level increases. In Figures , we show graphs of two numerical solutions of (2.) associated with the same boundary conditions. matrix size Fig. 2. Two numerical solutions of the equation (2.) with the same boundary conditions and ² = Fig. 4.2 Another two numerical solutions of the equation (2.) with the same boundary conditions. Here, ² =.5. 7
8 Secondly, we test if the numerical solutions of (.) converge to the solution of the steady state equation (2.). We input the programs the right-hand sides and boundary conditions computed from a known test function u, input the programs with various initial values, compute the numerical solutions S u at t =, 2, 5 and, and then compare the maximum errors between u and S u. We list the maximum errors for various initial values, different times t =, 2, 5 and and various values ɛ. Table 2.3. The maximum errors of numerical solutions of (.) with boundary conditions obtained from u(x, y) = x 4 + y 4 and initial value u(x, y) = (x 4 + y 4 )( sin(πx) sin(πy)) against the exact solution u(x, y) = x 4 + y 4. time ɛ = ɛ =. ɛ =. ɛ =. ɛ =.5 t = t = t = t = Table 2.4. The maximum errors of numerical solutions of (.) with boundary conditions obtained from u(x, y) = x 4 + y 4 and initial value u(x, y) = against the exact solution u(x, y) = x 4 + y 4. time ɛ = ɛ =. ɛ =. ɛ =. ɛ =.5 t = t = t = t = Table 2.5. The maximum errors of numerical solutions of (.) with boundary conditions obtained from u(x, y) = sin(xy) and initial value u(x, y) = sin(xy)( sin(πx) sin(πy)) against the exact solution u(x, y) = sin(xy). time ɛ = ɛ =. ɛ =. ɛ =. ɛ =.5 t = t = t = t =
9 Table 2.6. The maximum errors of numerical solutions of (.) with boundary conditions obtained from u(x, y) = sin(xy) and initial value u(x, y) = against the exact solution u(x, y) = sin(xy). time ɛ = ɛ =. ɛ =. ɛ =. ɛ =.5 t = t = t = t = From Tables , we see that for ɛ >.5, the numerical solutions of (.) converges to the exact solution of the steady state equation (2.) as t increases. This verifies Theorem 2.4. However, when ɛ =.5, the numerical solution of (.) in Table does not converge to x 4 + y 4 and sin(xy). Instead, it converges to the other numerical solution of (2.) as shown in the right graph of Fig. 2. and Fig We have evaluated their energy functional at the solutions. The energy of the solutions in the right graph of Figs 2. and 2.2 is smaller than the energy of the exact solutions. 9
10 3. Numerical Simulation of the Equation (.2) In this section, we study the numerical solution to both the steady state equation (3.) below and time dependent equation corresponding to (3.). (3.) 2 u + ɛ ( u 2 ) u) = g, 2 u = u (x), n u = u 2 (x). The study is similar to 2 although it is more complicated. Let us outline the main steps and ingredients. We first express the equation (3.) in a weak formulation. Find u V H 2 () such that in (3.) a 2 (u, v) + b(u, u, v) f(u), v g, v = for all v V H 2 (), where b(θ, φ, ψ) = ɛ 2 ( θ 2 ) φ ψdx and ɛ 2 f(u), v = u 2 φ vdx+ 2 u φ u vdx φ 2 u vdx u φ φ vdx φ vdx + φ 2 φ vdx with φ being the biharmonic function with the original boundary conditions. Let V = S r 3r( +) H 2 () be a spline subspace in H 2 () as before. We may use the Brouwer fixed point theorem to prove Theorem 3. (cf. [Lai and Wenston ] for such kind of proof for the nonlinear biharmonic equation associated with Navier-Stokes equations). Theorem 3.. For any ε >, there exists a weak solution S u V satisfying (5.). Note that the proof of Lemma 3.4 in [LLW3] can be used to prove the boundedness of S u. That is, Lemma 3.. Let S u be a weak solution satisfying (3.) with V = S r 3r( +) H 2 (). Then S u 2 dx C(ɛ). Similar to Theorem 3.4 in [LLW3], we have
11 Theorem 3.2. If ɛ is not very small, then the weak solution S u is unique. We next study the approximation of S u to u. Consider ɛ 2 a 2 (u S u, u S u ) = ɛ 2 a 2 (u S u, u S) + ɛ 2 a 2 (u S u, S S u ) for any spline function S V. Note that we have a 2 (u S u, S S u ) =b(s u, S u, S S u ) b(u, u, S S u ) + f(s u ) f(u), S S u =b(s u, S u, u S u ) + b(s u, S u, S u) b(u, u, u S u ) + b(u, u, u S) + ɛ 2 ( S u 2 u 2 ), φ (S S u ) + ɛ 2 φ 2, (S u u) (S S u ) + 2 ɛ 2 (S u u) φ, φ (S S u ) + 2 ɛ 2 (S u u) φ, S u (S S u ) + 2 ɛ 2 (S u u) (S S u ), φ u. The last four terms in the above can be bounded by C(m, ɛ, K, ) S S u 2, S u u 2, as in 4, where C(m, ɛ, K, ) is a positive constant. Also, the first four terms on the righthand side in the above can be treated exactly as in 4 with φ in place of φ. Therefore we should get the same result. That is, Theorem 3.3. Suppose that ɛ is not very small. Suppose that u H k () with k 3. Then the bivariate spline solution S u satisfying (3.) approximates the weak solution u such that S u u 2, C + k 2 u k,. Finally we discuss the implementation of S3( +) to numerically solve both the equation (.2) and its steady state equation (3.). We use S3( +) to discretize the space variables of the equations (.2) and (3.) by Galerkins s method and use Crank-Nicolson s method to discretize the time variable of the equations (.2). For each step we employ the Newton method to iterate the nonlinear term in (.2) and (3.). We set the error between iterations to be 8 for a stopping criterion. For the time dependent problem (.2), we use the solution at t = t i to be an initial guess for the solution for t = t i+. Our numerical scheme for (.2) is ) Start with S u, S3( +) which is the numerical solution of a biharmonic equation: S u, vdx = gvdx, for all v S 3( +) H 2 () u = u (x), n u = u 2 (x) 2) For k =,,, solve S u,k S3( +) satisfying S u,k vdx = gvdx f(s u,k ) vdx, u = u (x), n u = u 2 (x)
12 3) Letting S u,k = N c i,k φ i, if {c i,k } {c i,k } < 8 we stop the above iteration. i= Let us explain the computation of the nonlinear terms in (3.) in detail. Let S u be an approximation of u in the spline space S3( +). Note that f(u, ɛ) is a function of x u and y u. Thus, we write f(u, ɛ) = [f ( x u, x u)), f 2( x u, y u)] and we approximate f and f 2 by quadratic finite elements associated with the underlying triangulation +. That is, letting W = {w i, i =,, M} be the collection of all vertices of + and the midpoints of all edges of +, f ( x u, M y u) = i= f ( x u(w i), y u(w i))h i (x, y) and similar for f 2, where h i is the continuous quadratic finite element basis function satisfying h i (w j ) = δ ij for all i, j =,, M. We numerically solve equation (3.) with known exact solutions u(x, y) = x 4 + y 4 and u(x, y) = sin(xy). In the following Tables , we list the maximum errors of the numerical solutions against the exact solutions evaluated over points equallyspaced over [, ] [, ] for different values ɛ. Table 3.. The maximum errors of numerical solutions of (3.) against the exact solution u(x, y) = x 4 + y 4. matrix size ɛ =. ɛ =. ɛ =. ɛ = Table 3.2. The maximum errors of numerical solutions of (.9) against the exact solution u(x, y) = sin(xy). matrix size ɛ =. ɛ =. ɛ =. ɛ = From Tables 3. and 3.2, we see that for ɛ >., the numerical solutions of (3.) converges to the exact solution as the refinement level increases. However, when 2
13 Fig Numerical solution of equation (.9) with ² =.5.2. and ² = Fig Numerical solution of equation (.9) with ² =.2. and ² =. ²., the numerical solution of (3.) in Tables 3. and 3.2 does not converge to x4 + y 4 and sin(xy). Instead, it converges to the other numerical solution of (3.) as shown in Fig. 3. and Fig We have computed their energy functional. The energy of the solutions in the Figs are almost the same as the energy of the exact solutions x4 +y 4 and sin(xy), respectively. We have also applied the continuation or homotopy technique to solve the equation. That is, we started with the numerical solution associated with ² = as an initial value to compute the solution associated with ² =.. Then using the solution with ² =. as an initial value, we compute the numerical solution associated with ² =.. And so on. We were able to find a very good approximation of the exact solution for ² =. and ² =. for all refinements of underlying quadrangulations. References [] S. Allen and J. Cahn, A microscopic theory for the antiphase boundary motion and its application to antiphase domain coarsening, Acta Metallurgica, 27 (979), pp [2] S. Angenet and M. Gurtin, Multiphase thermo-mechanics with interfacial structure. 3
14 Evolution of an isothermal interface, Acta Rational Mech. Math. Anal. 8(989), pp [3] L. Bronsard and R. Kohn. Motion by mean curvature as the singular limit of Ginsburg- Landau dynamics, J.D.H. Eqs. 2(99), pp [4] M. Carme Calderer, C. Liu and K. Voss, Radial configurations of smectic A materials and focal conics, Submitted [5] M. Carme Calderer, C. Liu and Karl Voss, Smectic A Liquid Crystal Configurations with Interface Defects Submitted to Physica D (998). [6] Du, Q., M. D. Gunzburger, J. S. Peterson, Analysis and approximation of the Ginsburg Landau model of superconductivity, SIAM Review, 34(992) pp [7] W. E, Nonlinear Continuum Theory of Smectic-A Liquid Crystals, Preprint, 997. [8] P. Grisvard, Elliptic Problems on Nonsmooth Domains, Pitman, Boston, 985. [9] W.M. Jin, Singular Perturbations, Fold Energies and Micromagnetics, Ph. D. thesis, Courant Institute, 997. [] D. Kinderlehrer and C. Liu, Revisiting the Focal Conic Structures in Smectic A. Proc. Symposium on Elasticity to honor Professor J. L. Ericksen 996, ASME Mechanics and Materials Conference. [] R. Kohn, A singular perturbed variational problem describing thin film blisters, IMA workshop, Feb. 5 9, 996. [2] M. J. Lai and L. L. Schumaker, On the approximation power of splines on triangulated quadrangulations, SIAM J. Num. Anal., 36(999), pp [3] M. J. Lai, C. Liu and P. Wenston, On two nonlinear biharmonic evolution equations: existence, uniqueness and stability, to appear in Applicable Analysis, 23. [4] M. J. Lai and P. Wenston, Bivariate spline method for the steady state Navier-Stokes equations, Numerical Methods for PDE, 6(2), pp [5] M. Ortiz and G. Gioia, 994, The morphology and folding patterns of buckling-driven thin-film blisters, J. Mech. Phys. Solids, 42(994), pp [6] J. Rubinstein, P. Sternberg and J. Keller. Fast reaction, slow diffusion and curve shortening. SIAM J. Appl. Math., 49(989), pp [7] R. Kohn and P. Sternberg, Local minimizers and singular perturbations, Proc. Roy. Soc. Edinburgh Sect. A (989), pp [8] P. Sternberg, The effect of a singular perturbation on nonconvex variational problems, Arch. Rational. Mech. Anal, (988), pp
On Two Nonlinear Biharmonic Evolution Equations: Existence, Uniqueness and Stability
On Two Nonlinear Biharmonic Evolution Equations: Existence, Uniqueness and Stability Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We study the following two nonlinear evolution equations with a fourth
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationMaximum Norm Estimate for Bivariate Spline Solutions to Second Order Elliptic Partial Differential Equations in Non-divergence Form
Maximum Norm Estimate for Bivariate Spline Solutions to Second Order Elliptic Partial Differential Equations in Non-divergence Form Ming-Jun Lai January 2, 2017 Abstract The convergence of the bivariate
More informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationCoupled second order singular perturbations for phase transitions
Coupled second order singular perturbations for phase transitions CMU 06/09/11 Ana Cristina Barroso, Margarida Baía, Milena Chermisi, JM Introduction Let Ω R d with Lipschitz boundary ( container ) and
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationSpectral theory for magnetic Schrödinger operators and applicatio. (after Bauman-Calderer-Liu-Phillips, Pan, Helffer-Pan)
Spectral theory for magnetic Schrödinger operators and applications to liquid crystals (after Bauman-Calderer-Liu-Phillips, Pan, Helffer-Pan) Ryukoku (June 2008) In [P2], based on the de Gennes analogy
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 6 () 77 745 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Multiphase
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationNonlinear Schrödinger problems: symmetries of some variational solutions
Nonlinear Differ. Equ. Appl. (3), 5 5 c Springer Basel AG -97/3/35- published online April 3, DOI.7/s3--3- Nonlinear Differential Equations and Applications NoDEA Nonlinear Schrödinger problems: symmetries
More informationClassical solutions for the quasi-stationary Stefan problem with surface tension
Classical solutions for the quasi-stationary Stefan problem with surface tension Joachim Escher, Gieri Simonett We show that the quasi-stationary two-phase Stefan problem with surface tension has a unique
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationAnalysis of a non-isothermal model for nematic liquid crystals
Analysis of a non-isothermal model for nematic liquid crystals E. Rocca Università degli Studi di Milano 25th IFIP TC 7 Conference 2011 - System Modeling and Optimization Berlin, September 12-16, 2011
More informationLOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S PARABOLIC-HYPERBOLIC COMPRESSIBLE NON-ISOTHERMAL MODEL FOR LIQUID CRYSTALS
Electronic Journal of Differential Equations, Vol. 017 (017), No. 3, pp. 1 8. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu LOCAL WELL-POSEDNESS FOR AN ERICKSEN-LESLIE S
More informationIMA Preprint Series # 2143
A BASIC INEQUALITY FOR THE STOKES OPERATOR RELATED TO THE NAVIER BOUNDARY CONDITION By Luan Thach Hoang IMA Preprint Series # 2143 ( November 2006 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationGAKUTO International Series
GAKUTO International Series Mathematical Sciences and Applications, Vol.28(2008) Proceedings of Fourth JSIAM-SIMMAI Seminar on Industrial and Applied Mathematics, pp.139-148 A COMPUTATIONAL APPROACH TO
More informationUniform Convergence of a Multilevel Energy-based Quantization Scheme
Uniform Convergence of a Multilevel Energy-based Quantization Scheme Maria Emelianenko 1 and Qiang Du 1 Pennsylvania State University, University Park, PA 16803 emeliane@math.psu.edu and qdu@math.psu.edu
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationTHE EXISTENCE OF GLOBAL ATTRACTOR FOR A SIXTH-ORDER PHASE-FIELD EQUATION IN H k SPACE
U.P.B. Sci. Bull., Series A, Vol. 8, Iss. 2, 218 ISSN 1223-727 THE EXISTENCE OF GLOBAL ATTRACTOR FOR A SIXTH-ORDER PHASE-FIELD EQUATION IN H k SPACE Xi Bao 1, Ning Duan 2, Xiaopeng Zhao 3 In this paper,
More informationNUMERICAL APPROXIMATIONS OF ALLEN-CAHN AND CAHN-HILLIARD EQUATIONS
UMERICAL APPROXIMATIOS OF ALLE-CAH AD CAH-HILLIARD EQUATIOS JIE SHE AD XIAOFEG YAG Abstract. Stability analyses and error estimates are carried out for a number of commonly used numerical schemes for the
More informationRESEARCH ARTICLE. A strategy of finding an initial active set for inequality constrained quadratic programming problems
Optimization Methods and Software Vol. 00, No. 00, July 200, 8 RESEARCH ARTICLE A strategy of finding an initial active set for inequality constrained quadratic programming problems Jungho Lee Computer
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationHeat equations with singular potentials: Hardy & Carleman inequalities, well-posedness & control
Outline Heat equations with singular potentials: Hardy & Carleman inequalities, well-posedness & control IMDEA-Matemáticas & Universidad Autónoma de Madrid Spain enrique.zuazua@uam.es Analysis and control
More informationA Product Property of Sobolev Spaces with Application to Elliptic Estimates
Rend. Sem. Mat. Univ. Padova Manoscritto in corso di stampa pervenuto il 23 luglio 2012 accettato l 1 ottobre 2012 A Product Property of Sobolev Spaces with Application to Elliptic Estimates by Henry C.
More informationOn the Midpoint Method for Solving Generalized Equations
Punjab University Journal of Mathematics (ISSN 1016-56) Vol. 40 (008) pp. 63-70 On the Midpoint Method for Solving Generalized Equations Ioannis K. Argyros Cameron University Department of Mathematics
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationwith deterministic and noise terms for a general non-homogeneous Cahn-Hilliard equation Modeling and Asymptotics
12-3-2009 Modeling and Asymptotics for a general non-homogeneous Cahn-Hilliard equation with deterministic and noise terms D.C. Antonopoulou (Joint with G. Karali and G. Kossioris) Department of Applied
More informationA numerical solution of a Kawahara equation by using Multiquadric radial basis function
Mathematics Scientific Journal Vol. 9, No. 1, (013), 115-15 A numerical solution of a Kawahara equation by using Multiquadric radial basis function M. Zarebnia a, M. Takhti a a Department of Mathematics,
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationPattern formation in compressed elastic films on compliant substrates: an explanation of the herringbone structure
Pattern formation in compressed elastic films on compliant substrates: an explanation of the herringbone structure Robert V. Kohn Courant Institute, NYU SIAM Annual Meeting, July 2012 Joint work with Hoai-Minh
More informationC 0 IPG Method for Biharmonic Eigenvalue Problems
C IPG Method for Biharmonic Eigenvalue Problems Susanne C. Brenner, Peter Monk 2, and Jiguang Sun 3 Department of Mathematics and Center for Computation & Technology, Louisiana State University, Baton
More informationSharp blow-up criteria for the Davey-Stewartson system in R 3
Dynamics of PDE, Vol.8, No., 9-60, 011 Sharp blow-up criteria for the Davey-Stewartson system in R Jian Zhang Shihui Zhu Communicated by Y. Charles Li, received October 7, 010. Abstract. In this paper,
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationContinuum Methods 1. John Lowengrub Dept Math UC Irvine
Continuum Methods 1 John Lowengrub Dept Math UC Irvine Membranes, Cells, Tissues complex micro-structured soft matter Micro scale Cell: ~10 micron Sub cell (genes, large proteins): nanometer Macro scale
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
More informationCURRICULUM VITAE SHIBIN DAI
CURRICULUM VITAE SHIBIN DAI Address Department of Mathematics, The University of Alabama, Tuscaloosa, AL 35487 Email: sdai4@ua.edu Web: http://sdai.people.ua.edu Education Ph.D. in Applied Mathematics,
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationHIGHLY EFFICIENT AND ACCURATE NUMERICAL SCHEMES FOR THE EPITAXIAL THIN FILM GROWTH MODELS BY USING THE SAV APPROACH
HIGHLY EFFICIENT AND ACCURATE NUMERICAL SCHEMES FOR THE EPITAXIAL THIN FILM GROWTH MODELS BY USING THE SAV APPROACH QING CHENG, JIE SHEN, XIAOFENG YANG Abstract. We develop in this paper highly efficient,
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
More informationTwo-dimensional examples of rank-one convex functions that are not quasiconvex
ANNALES POLONICI MATHEMATICI LXXIII.3(2) Two-dimensional examples of rank-one convex functions that are not quasiconvex bym.k.benaouda (Lille) andj.j.telega (Warszawa) Abstract. The aim of this note is
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationAn a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element
Calcolo manuscript No. (will be inserted by the editor) An a posteriori error estimate and a Comparison Theorem for the nonconforming P 1 element Dietrich Braess Faculty of Mathematics, Ruhr-University
More informationSymmetry of entire solutions for a class of semilinear elliptic equations
Symmetry of entire solutions for a class of semilinear elliptic equations Ovidiu Savin Abstract. We discuss a conjecture of De Giorgi concerning the one dimensional symmetry of bounded, monotone in one
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationPolynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes
Polynomial Preserving Recovery for Quadratic Elements on Anisotropic Meshes Can Huang, 1 Zhimin Zhang 1, 1 Department of Mathematics, Wayne State University, Detroit, Michigan 480 College of Mathematics
More informationCurriculum Vitae Mingji Zhang
Curriculum Vitae Mingji Zhang Contact Information Department of Mathematics Michigan State University 619 Red Cedar Road, East Lansing, MI 48824 Phone: (517)884-1455 Fax: (517)432-1562 Email: mzhang0129@gmail.com
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationNONCLASSICAL SHOCK WAVES OF CONSERVATION LAWS: FLUX FUNCTION HAVING TWO INFLECTION POINTS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 149, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONCLASSICAL
More informationCURRICULUM VITAE SHIBIN DAI
CURRICULUM VITAE SHIBIN DAI Address Department of Mathematics, University of Alabama, Tuscaloosa, AL 35487 Phone: (office) 205-348-5156 Email: sdai4@ua.edu Web: http://sdai.people.ua.edu Education Ph.D.
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationA Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term
A Ginzburg-Landau Type Problem for Nematics with Highly Anisotropic Elastic Term Peter Sternberg In collaboration with Dmitry Golovaty (Akron) and Raghav Venkatraman (Indiana) Department of Mathematics
More informationHONGJIE DONG. Assistant Professor of Applied Mathematics Division of Applied Mathematics Brown University
HONGJIE DONG Assistant Professor of Applied Mathematics Division of Applied Mathematics Brown University Address Division of Applied Mathematics 182 George Street Providence, RI 02912 Phone: (401) 863-7297
More informationTRAVELLING WAVES FOR A NONLOCAL DOUBLE OBSTACLE PROBLEM
TRAVELLING WAVES FOR A NONLOCAL DOUBLE OBSTACLE PROBLEM PAUL C. FIFE, Mathematics Department, University of Utah, Salt Lake City, Utah 84112, USA Abstract Existence, uniqueness, and regularity properties
More informationarxiv: v1 [math.ap] 28 Mar 2014
GROUNDSTATES OF NONLINEAR CHOQUARD EQUATIONS: HARDY-LITTLEWOOD-SOBOLEV CRITICAL EXPONENT VITALY MOROZ AND JEAN VAN SCHAFTINGEN arxiv:1403.7414v1 [math.ap] 28 Mar 2014 Abstract. We consider nonlinear Choquard
More informationNew exact multiplicity results with an application to a population model
New exact multiplicity results with an application to a population model Philip Korman Department of Mathematical Sciences University of Cincinnati Cincinnati, Ohio 45221-0025 Junping Shi Department of
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationEXISTENCE OF WEAK SOLUTIONS FOR A NONUNIFORMLY ELLIPTIC NONLINEAR SYSTEM IN R N. 1. Introduction We study the nonuniformly elliptic, nonlinear system
Electronic Journal of Differential Equations, Vol. 20082008), No. 119, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu login: ftp) EXISTENCE
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationUniversität Regensburg Mathematik
Universität Regensburg Mathematik On the stable discretization of strongly anisotropic phase field models with applications to crystal growth John W. Barrett, Harald Garcke and Robert Nürnberg Preprint
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationSolution of fractional oxygen diffusion problem having without singular kernel
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 1 (17), 99 37 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa Solution of fractional oxygen diffusion
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationSEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT
Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC
More informationIMA Preprint Series # 2260
LONG-TIME BEHAVIOR OF HYDRODYNAMIC SYSTEMS MODELING THE NEMATIC LIUID CRYSTALS By Hao Wu Xiang Xu and Chun Liu IMA Preprint Series # 2260 ( June 2009 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY
More informationA Numerical Algorithm for Finding Positive Solutions of Superlinear Dirichlet Problem
ISSN 1479-3889 (print), 1479-3897 (online) International Journal of Nonlinear Science Vol.3(2007) No.1,pp.27-31 A Numerical Algorithm for Finding Positive Solutions of Superlinear Dirichlet Problem G.A.Afrouzi
More informationWeak solutions for the Cahn-Hilliard equation with degenerate mobility
Archive for Rational Mechanics and Analysis manuscript No. (will be inserted by the editor) Shibin Dai Qiang Du Weak solutions for the Cahn-Hilliard equation with degenerate mobility Abstract In this paper,
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationOn Moving Ginzburg-Landau Vortices
communications in analysis and geometry Volume, Number 5, 85-99, 004 On Moving Ginzburg-Landau Vortices Changyou Wang In this note, we establish a quantization property for the heat equation of Ginzburg-Landau
More informationPattern Formation in a Compressed Elastic Film on a Compliant Substrate
Pattern Formation in a Compressed Elastic Film on a Compliant Substrate Robert V. Kohn Courant Institute, NYU ICERM Workshop, Providence, May 30 - June 1, 2012 Joint work with Hoai-Minh Nguyen Funding
More informationFormulation of the problem
TOPICAL PROBLEMS OF FLUID MECHANICS DOI: https://doi.org/.43/tpfm.27. NOTE ON THE PROBLEM OF DISSIPATIVE MEASURE-VALUED SOLUTIONS TO THE COMPRESSIBLE NON-NEWTONIAN SYSTEM H. Al Baba, 2, M. Caggio, B. Ducomet
More informationTHE CAUCHY PROBLEM AND STEADY STATE SOLUTIONS FOR A NONLOCAL CAHN-HILLIARD EQUATION
Electronic Journal of Differential Equations, Vol. 24(24), No. 113, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) THE CAUCHY
More informationMEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW
MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW GREGORY DRUGAN AND XUAN HIEN NGUYEN Abstract. We present two initial graphs over the entire R n, n 2 for which the mean curvature flow
More informationDissipative solutions for a hyperbolic system arising in liquid crystals modeling
Dissipative solutions for a hyperbolic system arising in liquid crystals modeling E. Rocca Università degli Studi di Pavia Workshop on Differential Equations Central European University, Budapest, April
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More information