Home About us Contact
|
Home About us Contact | ||
|
|
||
Smooth Solutions (smooth + solution)
Selected AbstractsLocal existence for the one-dimensional Vlasov,Poisson system with infinite massMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 5 2007Stephen Pankavich Abstract A collisionless plasma is modelled by the Vlasov,Poisson system in one dimension. We consider the situation in which mobile negative ions balance a fixed background of positive charge, which is independent of space and time, as ,x, , ,. Thus, the total positive charge and the total negative charge are both infinite. Smooth solutions with appropriate asymptotic behaviour are shown to exist locally in time, and criteria for the continuation of these solutions are established. Copyright © 2006 John Wiley & Sons, Ltd. [source] Radial basis functions for solving near singular Poisson problemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 5 2003C. S. Chen Abstract In this paper, we investigate the use of radial basis functions for solving Poisson problems with a near-singular inhomogeneous source term. The solution of the Poisson problem is first split into two parts: near-singular solution and smooth solution. A method for evaluating the near-singular particular solution is examined. The smooth solution is further split into a particular solution and a homogeneous solution. The MPS-DRM approach is adopted to evaluate the smooth solution. Copyright © 2003 John Wiley & Sons, Ltd. [source] State-space time integration with energy control and fourth-order accuracy for linear dynamic systemsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2006Steen Krenk Abstract A fourth-order accurate time integration algorithm with exact energy conservation for linear structural dynamics is presented. It is derived by integrating the phase-space representation and evaluating the resulting displacement and velocity integrals via integration by parts, substituting the time derivatives from the original differential equations. The resulting algorithm has an exact energy equation, in which the change of energy is equal to the work of the external forces minus a quadratic form of the damping matrix. This implies unconditional stability of the algorithm, and the relative phase error is of fourth-order. An optional high-frequency algorithmic damping is constructed by optimal combination of three different damping matrices, each proportional to either the mass or the stiffness matrix. This leads to a modified form of the undamped algorithm with scalar weights on some of the matrices introducing damping of fourth-order in the frequency. Thus, the low-frequency response is virtually undamped, and the algorithm remains third-order accurate even when algorithmic damping is included. The accuracy of the algorithm is illustrated by an application to pulse propagation in an elastic medium, where the algorithmic damping is used to reduce dispersion due to the spatial discretization, leading to a smooth solution with a clearly defined wave front. Copyright © 2005 John Wiley & Sons, Ltd. [source] Limiting Amplitude principle in the scattering by wedgesMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 10 2006A. I. Komech Abstract We consider a nonstationary scattering of plane waves by a wedge. We prove that the Sommerfeld-type integral, constructed in (Math. Meth. Appl. Sci. 2005; 28:147,183; Proc. Int. Seminar ,Day on Diffraction-2003', University of St. Petersburg, 2003; 151,162), is a classical smooth solution from a functional space, and prove the Limiting Amplitude principle. Copyright © 2006 John Wiley & Sons, Ltd. [source] Higher-order analogues of the Tracy-Widom distribution and the Painlevé II hierarchyCOMMUNICATIONS ON PURE & APPLIED MATHEMATICS, Issue 3 2010Tom Claeys We study Fredholm determinants related to a family of kernels that describe the edge eigenvalue behavior in unitary random matrix models with critical edge points. The kernels are natural higher-order analogues of the Airy kernel and are built out of functions associated with the Painlevé I hierarchy. The Fredholm determinants related to those kernels are higher-order generalizations of the Tracy-Widom distribution. We give an explicit expression for the determinants in terms of a distinguished smooth solution to the Painlevé II hierarchy. In addition, we compute large gap asymptotics for the Fredholm determinants. © 2009 Wiley Periodicals, Inc. [source] The design of an optimal filter for monthly GRACE gravity modelsGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 2 2008R. Klees SUMMARY Most applications of the publicly released Gravity Recovery and Climate Experiment monthly gravity field models require the application of a spatial filter to help suppressing noise and other systematic errors present in the data. The most common approach makes use of a simple Gaussian averaging process, which is often combined with a ,destriping' technique in which coefficient correlations within a given degree are removed. As brute force methods, neither of these techniques takes into consideration the statistical information from the gravity solution itself and, while they perform well overall, they can often end up removing more signal than necessary. Other optimal filters have been proposed in the literature; however, none have attempted to make full use of all information available from the monthly solutions. By examining the underlying principles of filter design, a filter has been developed that incorporates the noise and full signal variance,covariance matrix to tailor the filter to the error characteristics of a particular monthly solution. The filter is both anisotropic and non-symmetric, meaning it can accommodate noise of an arbitrary shape, such as the characteristic stripes. The filter minimizes the mean-square error and, in this sense, can be considered as the most optimal filter possible. Through both simulated and real data scenarios, this improved filter will be shown to preserve the highest amount of gravity signal when compared to other standard techniques, while simultaneously minimizing leakage effects and producing smooth solutions in areas of low signal. [source] A hybridizable discontinuous Galerkin method for linear elasticityINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2009S.-C. Soon Abstract This paper describes the application of the so-called hybridizable discontinuous Galerkin (HDG) method to linear elasticity problems. The method has three significant features. The first is that the only globally coupled degrees of freedom are those of an approximation of the displacement defined solely on the faces of the elements. The corresponding stiffness matrix is symmetric, positive definite, and possesses a block-wise sparse structure that allows for a very efficient implementation of the method. The second feature is that, when polynomials of degree k are used to approximate the displacement and the stress, both variables converge with the optimal order of k+1 for any k,0. The third feature is that, by using an element-by-element post-processing, a new approximate displacement can be obtained that converges at the order of k+2, whenever k,2. Numerical experiments are provided to compare the performance of the HDG method with that of the continuous Galerkin (CG) method for problems with smooth solutions, and to assess its performance in situations where the CG method is not adequate, that is, when the material is nearly incompressible and when there is a crack. Copyright © 2009 John Wiley & Sons, Ltd. [source] Smooth finite element methods: Convergence, accuracy and propertiesINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 2 2008Hung Nguyen-Xuan Abstract A stabilized conforming nodal integration finite element method based on strain smoothing stabilization is presented. The integration of the stiffness matrix is performed on the boundaries of the finite elements. A rigorous variational framework based on the Hu,Washizu assumed strain variational form is developed. We prove that solutions yielded by the proposed method are in a space bounded by the standard, finite element solution (infinite number of subcells) and a quasi-equilibrium finite element solution (a single subcell). We show elsewhere the equivalence of the one-subcell element with a quasi-equilibrium finite element, leading to a global a posteriori error estimate. We apply the method to compressible and incompressible linear elasticity problems. The method can always achieve higher accuracy and convergence rates than the standard finite element method, especially in the presence of incompressibility, singularities or distorted meshes, for a slightly smaller computational cost. It is shown numerically that the one-cell smoothed four-noded quadrilateral finite element has a convergence rate of 2.0 in the energy norm for problems with smooth solutions, which is remarkable. For problems with rough solutions, this element always converges faster than the standard finite element and is free of volumetric locking without any modification of integration scheme. Copyright © 2007 John Wiley & Sons, Ltd. [source] Local maximum-entropy approximation schemes: a seamless bridge between finite elements and meshfree methodsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 13 2006M. Arroyo Abstract We present a one-parameter family of approximation schemes, which we refer to as local maximum-entropy approximation schemes, that bridges continuously two important limits: Delaunay triangulation and maximum-entropy (max-ent) statistical inference. Local max-ent approximation schemes represent a compromise,in the sense of Pareto optimality,between the competing objectives of unbiased statistical inference from the nodal data and the definition of local shape functions of least width. Local max-ent approximation schemes are entirely defined by the node set and the domain of analysis, and the shape functions are positive, interpolate affine functions exactly, and have a weak Kronecker-delta property at the boundary. Local max-ent approximation may be regarded as a regularization, or thermalization, of Delaunay triangulation which effectively resolves the degenerate cases resulting from the lack or uniqueness of the triangulation. Local max-ent approximation schemes can be taken as a convenient basis for the numerical solution of PDEs in the style of meshfree Galerkin methods. In test cases characterized by smooth solutions we find that the accuracy of local max-ent approximation schemes is vastly superior to that of finite elements. Copyright © 2005 John Wiley & Sons, Ltd. [source] Some results on the accuracy of an edge-based finite volume formulation for the solution of elliptic problems in non-homogeneous and non-isotropic mediaINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2009Darlan Karlo Elisiário de Carvalho Abstract The numerical simulation of elliptic type problems in strongly heterogeneous and anisotropic media represents a great challenge from mathematical and numerical point of views. The simulation of flows in non-homogeneous and non-isotropic porous media with full tensor diffusion coefficients, which is a common situation associated with the miscible displacement of contaminants in aquifers and the immiscible and incompressible two-phase flow of oil and water in petroleum reservoirs, involves the numerical solution of an elliptic type equation in which the diffusion coefficient can be discontinuous, varying orders of magnitude within short distances. In the present work, we present a vertex-centered edge-based finite volume method (EBFV) with median dual control volumes built over a primal mesh. This formulation is capable of handling the heterogeneous and anisotropic media using structured or unstructured, triangular or quadrilateral meshes. In the EBFV method, the discretization of the diffusion term is performed using a node-centered discretization implemented in two loops over the edges of the primary mesh. This formulation guarantees local conservation for problems with discontinuous coefficients, keeping second-order accuracy for smooth solutions on general triangular and orthogonal quadrilateral meshes. In order to show the convergence behavior of the proposed EBFV procedure, we solve three benchmark problems including full tensor, material heterogeneity and distributed source terms. For these three examples, numerical results compare favorably with others found in literature. A fourth problem, with highly non-smooth solution, has been included showing that the EBFV needs further improvement to formally guarantee monotonic solutions in such cases. Copyright © 2008 John Wiley & Sons, Ltd. [source] On the hyperbolic system of a mixture of Eulerian fluids: a comparison between single- and multi-temperature modelsMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 7 2007Tommaso Ruggeri Abstract The first rational model of homogeneous mixtures of fluids was proposed by Truesdell in the context of rational thermodynamics. Afterwards, two different theories were developed: one with a single-temperature (ST) field of the mixture and the other one with several temperatures. The two systems are from the mathematical point of view completely different and the relationship between their solutions was not clarified. In this paper, the hyperbolic multi-temperature (MT) system of a mixture of Eulerian fluids will be explained and it will be shown that the corresponding single-temperature differential system is a principal subsystem of the MT one. As a consequence, the subcharacteristic conditions for characteristic speeds hold and this gives an upper-bound esteem for pulse speeds in an ST model. Global behaviour of smooth solutions for large time for both systems will also be discussed through the application of the Shizuta,Kawashima condition. Finally, as an application, the particular case of a binary mixture is considered. Copyright © 2006 John Wiley & Sons, Ltd. [source] The asymptotic behaviour of global smooth solutions to the multi-dimensional hydrodynamic model for semiconductorsMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 14 2003Ling Hsiao Abstract We establish the global existence of smooth solutions to the Cauchy problem for the multi-dimensional hydrodynamic model for semiconductors, provided that the initial data are perturbations of a given stationary solutions, and prove that the resulting evolutionary solution converges asymptotically in time to the stationary solution exponentially fast. Copyright © 2003 John Wiley & Sons, Ltd. [source] Asymptotic behaviour of global smooth solutionsto the multidimensional hydrodynamic model for semiconductorsMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 8 2002Ling Hsiao Abstract In this paper, we study asymptotic behaviour of the global smooth solutions to the multidimensional hydrodynamic model for semiconductors. We prove that the solution of the problem converges to a stationary solution time asymptotically exponentially fast. Copyright © 2002 John Wiley & Sons, Ltd. [source] Thermoelasticity with second sound,exponential stability in linear and non-linear 1-dMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 5 2002Reinhard Racke We consider linear and non-linear thermoelastic systems in one space dimension where thermal disturbances are modelled propagating as wave-like pulses travelling at finite speed. This removal of the physical paradox of infinite propagation speed in the classical theory of thermoelasticity within Fourier's law is achieved using Cattaneo's law for heat conduction. For different boundary conditions, in particular for those arising in pulsed laser heating of solids, the exponential stability of the now purely, but slightly damped, hyperbolic linear system is proved. A comparison with classical hyperbolic,parabolic thermoelasticity is given. For Dirichlet type boundary conditions,rigidly clamped, constant temperature,the global existence of small, smooth solutions and the exponential stability are proved for a non-linear system. Copyright © 2002 John Wiley & Sons, Ltd. [source] Exchange of conserved quantities, shock loci and Riemann problemsMATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 13 2001Michael Sever Systems of conservation laws admitting extensions, such as entropy density/flux functions, generate related systems obtained by exchanging the extension with one of the constituent equations. Often if not always, the smooth solutions of the two systems coincide, and weak solutions of one system containing only small discontinuities are approximate weak solutions of the other. The adiabatic approximation for the Euler system illustrates the utility of this procedure. Such an exchange of conserved quantities preserves hyperbolicity and genuine non-linearity in the sense of Lax. On the other hand, the topological structure of the shock locus of a point in phase space and the solvability of Riemann problems in the large can be strongly affected. A discussion of when and how this occurs is given here. In this paper the exchange of conserved quantities is conveniently described by a simple homotopy in an extended version of the usual ,symmetric variables'. A dynamical system in phase space is constructed, the trajectories of which describe the Hugoniot locus of a fixed point in phase space at each state of the homotopy. The appearance of critical points for this dynamical system is identified with the alteration of the topological structure of the Hugoniot locus by the exchange of conserved quantities. Copyright © 2001 John Wiley & Sons, Ltd. [source] Hydrodynamic limits with shock waves of the Boltzmann equationCOMMUNICATIONS ON PURE & APPLIED MATHEMATICS, Issue 3 2005Shi-Hsien Yu We show that piecewise smooth solutions with shocks of the Euler equations in gas dynamics can be obtained as the zero Knudsen number limit of solutions of the Boltzmann equation for hard sphere collision model. The construction of the Boltzmann solutions is done in two steps. First we introduce a generalized Hilbert expansion with shock layer correction to construct approximations to the solutions of the Boltzmann equations with small Knudsen numbers. We then apply the recently developed macro-micro decomposition and energy method for Boltzmann shock layers to construct the exact Boltzmann solutions through the stability analysis. © 2004 Wiley Periodicals, Inc. [source] Low-curvature image simplifiers: Global regularity of smooth solutions and Laplacian limiting schemesCOMMUNICATIONS ON PURE & APPLIED MATHEMATICS, Issue 6 2004Andrea L. Bertozzi We consider a class of fourth-order nonlinear diffusion equations motivated by Tumblin and Turk's "low-curvature image simplifiers" for image denoising and segmentation. The PDE for the image intensity u is of the form where g(s) = k2/(k2 + s2) is a "curvature" threshold and , denotes a fidelity-matching parameter. We derive a priori bounds for ,u that allow us to prove global regularity of smooth solutions in one space dimension, and a geometric constraint for finite-time singularities from smooth initial data in two space dimensions. This is in sharp contrast to the second-order Perona-Malik equation (an ill-posed problem), on which the original LCIS method is modeled. The estimates also allow us to design a finite difference scheme that satisfies discrete versions of the estimates, in particular, a priori bounds on the smoothness estimator in both one and two space dimensions. We present computational results that show the effectiveness of such algorithms. Our results are connected to recent results for fourth-order lubrication-type equations and the design of positivity-preserving schemes for such equations. This connection also has relevance for other related fourth-order imaging equations. © 2004 Wiley Periodicals, Inc. [source] | |||||||||||||