Essentially optimal explicit Runge–Kutta methods with application to hyperbolic–parabolic equations

Optimal explicit Runge–Kutta methods consider more stages in order to include a particular spectrum in their stability domain and thus reduce time-step restrictions. This idea, so far used mostly for real-line spectra, is generalized to more general spectra in the form of a thin region. In thin regions the eigenvalues may extend away from the real axis into the imaginary plane. We give a direct characterization of optimal stability polynomials containing a maximal thin region and calculate these polynomials for various cases. Semi-discretizations of hyperbolic–parabolic equations are a relevant application which exhibit a thin region spectrum. As a model, linear, scalar advection–diffusion is investigated. The second-order-stabilized explicit Runge–Kutta methods derived from the stability polynomials are applied to advection–diffusion and compressible, viscous fluid dynamics in numerical experiments. Due to the stabilization the time step can be controlled solely from the hyperbolic CFL condition even in the presence of viscous fluxes.     

An algorithm for decomposing a polynomial system into normal ascending sets

We present an algorithm to decompose a polynomial system into a finite set of normal ascending sets such that the set of the zeros of the polynomial system is the union of the sets of the regular zeros of the normal ascending sets.If the polynomial system is zero dimensional,the set of the zeros of the polynomials is the union of the sets of the zeros of the normal ascending sets.     

Terminal control in flat systems

The article considers the application of dynamical feedback control to a special class of nonlinear dynamical systems — so-called flat systems. Flat system, linearizing output, and Lie-Backlund isomorphism are among the concepts reviewed. A man-machine algorithm is proposed for solving the terminal control problem for an arbitrary flat system with known linearizing outputs. __________ Translated from Nelineinaya Dinamika i Upravlenie, No. 3, pp. 191–200, 2003.     

Estimates from Below for the Spectral Function and for the Remainder in Local Weyl’s Law

We obtain asymptotic lower bounds for the spectral function of the Laplacian and for the remainder in local Weyl’s law on manifolds. In the negatively curved case, thermodynamic formalism is applied to improve the estimates. Key ingredients of the proof include the wave equation parametrix, a pretrace formula and the Dirichlet box principle. Our results develop and extend the unpublished thesis of A. Karnaukh [Ka]. The first author was supported by NSERC, FQRNT, Alfred P. Sloan Foundation Fellowship and Dawson Fellowship. The second author was supported by NSERC and FQRNT. Received: April 2006 Revision: October 2006 Accepted: October 2006     

Change in the magnetic state of a surface upon formation of hydrogen bonds during adsorption

It is established that reversible adsorption of water and methyl alcohol molecules, occurring via formation of hydrogen bonds, changes the dynamic properties of domain walls in the surface region of soft ferromagnets, as well as their initial static magnetic susceptibility. A mechanism is proposed for the effects revealed.     

A model for the interaction of near-infrared laser pulses with metal powders in selective laser sintering

A thermal model of the interaction of pulsed near-infrared laser radiation from a Nd:YAG laser was made, taking the measured powder properties such as reflectance, optical penetration depth and thermal conductivity into account. It allows an estimation of the evolution of two different temperatures: the average temperature of the powder (taken over the grains in a volume given by the laser beam diameter and the optical penetration depth) and the temperature distinction within a single grain. It showed that in pulsed mode consolidation can be achieved at much lower average power as the surface of the powder particles are molten but their cores remain at nearly room temperature. This leads to a much lower average temperature and therefore a dramatic decrease in residual thermal stresses in the finished piece. The results of the model were experimentally tested and confirmed. Received: 26 July 2001 / Accepted: 23 November 2001 / Published online: 23 January 2002     

On the Use of the Restitution Condition in Flexible Body Dynamics

The difference between the classical treatment offlexible body impact and the treatment of impact in flexiblemultibody dynamics is due to several fundamental reasons. Inthe classical impact theory, simple structures such as beamsand plates are used. Infinite dimensional models can bedeveloped for these simple structural elements to study theimpact dynamics and the wave propagation problem. Flexiblemultibody impact problems, on the other hand, involve bodieswith complex geometry that cannot be modeled using infinitenumber of degrees of freedom. Furthermore, the classicalimpact theory has been mainly concerned with the impactbetween a rigid mass that moves without constraints beforeit impacts a simple flexible structure. This is not amultibody simulation scenario in which the impact occursbetween kinematically constrained bodies that are subjectedto impulsive constraint forces in addition to the impactforces. These constraint forces can influence the motion ofthe two bodies immediately after impact, and as aconsequence, the simple classical theory scenario of impactdoes not apply. It is the objective of this paper to discussthe use of the restitution condition in flexible multibodyimpact problems and demonstrate that the use of thisapproach does not exclude the classical formulation.Nonetheless, the impulse momentum balance approach can serveas an effective and efficient procedure for solving theimpact problem in finite dimensional models that do not obeythe classical wave theory. Energy results of simplestructural elements are presented in order to demonstratethe consistency of using the impulse momentum balanceapproach in solving impact problems in finite dimensionalflexible body applications.     

Geometrically qualified ESPI vibration analysis of an engine

This paper presents an application in the automotive industry where a combination of electronic speckle-pattern interferometry and laser doppler velocimetry were used at a critical stage in the design process of an internal combustion engine. Combined deformation and surface relief measurements were used to study the phase and amplitude of deformation of a vibrating engine. The relief data was combined with the interferometer geometry and used to geometrically correct the deformation data, in an effort to improve accuracy. The measurements allowed rapid identification and quantification of design weaknesses, particularly those causing undesirable resonances. This led to a significant reduction in the design time and lowering of costs, when compared with existing design optimisation methods.     

Gas solubility of carbon dioxide in poly(lactic acid) at high pressures: Thermal treatment effect

The sorption of carbon dioxide in glassy Poly(lactic acid) (PLA) films was studied by quartz crystal microbalance (QCM) at high pressures. Two thermal treatments, melted and quenched, were performed in PLA with two different L:D contents, 80:20 and 98:2, films and compared with a third thermal protocol, annealed, and used in a previous work. The results obtained show that for pressures higher than 2 MPa, the carbon dioxide solubility is larger in PLA 80:20 than in PLA 98:2, indicating that the L:D plays a dominant role on this property. The thermal treatments only affect the gas solubility in PLA 98:2. Sorption isotherms at temperatures 303, 313, and 323 K, below the glass transition temperature of the polymer, and pressures up to 5 MPa were measured and analyzed with three different models, the dual‐mode sorption model, the Flory–Huggins equation, and a modified dual‐mode sorption model where the Henry's law term was substituted by the Flory–Huggins equation. This last model performs especially well for CO₂ in PLA 80:20, due to the convex upward curvature of the solubility isotherms for that system. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 616–625, 2007