A finite element method for second order nonvariational elliptic problems

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of “finite element Hessian” and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasi-linear PDE, all in nonvariational form.

Synthesis and characterization of thermally stable second-order nonlinear optical side-chain polyimides containing thiazole and benzothiazole push-pull chromophores

Push-pull nonlinear optical (NLO) chromophores containing thiazole and benzothiazole acceptors were synthesized and characterized. Using these chromophores a series of second-order NLO polyimides were Successfully prepared from 4,4'-(hexafluoroisopropylidene) diphthalic anhydride (6FDA), pyromellitic dianhydride (PMDA) and 3,3'4,4'-benzophenone tetracarboxylic dianhydride (BTDA) by a standard condensation polymerization technique. These polyimides exhibit high glass transition temperatures ranging from 160 to 188 degrees C. UV-vis spectrum of polyimide exhibited a slight blue shift and decreases in absorption due to birefringence. From the order parameters, it was found that chromophores were aligned effectively. Using in situ poling and temperature ramping technique, the optical temperatures for corona poling were obtained. It was found that the optimal temperatures of polyimides approach their glass transition temperatures. These polyimides demonstrate relatively large d(33) values range between 35.15 and 45.20 pm/V at 532 nm. (C) 2008 Elsevier B.V. All rights reserved.

Second order Coriolis resonance in the v2, v5 Raman bands of CH3F

The vibration-rotation Raman spectrum of the ν2 and ν5 fundamentals of CH3F is reported, from 1320 to 1640 cm−1, with a resolution of about 0.3 cm−1. The Coriolis resonance between the two bands leads to many perturbation-allowed transitions. Where the resonance is still sufficiently weak that the quantum number K′ retains its meaning, perturbation-allowed transitions are observed for all values of ΔK from +4 to −4; in regions of strong resonance, however, we can only say that the observed transitions obey the selection rule Δ(k−l) = 0 or ±3. The spectrum has been analyzed by band contour simulation using a computer program based on exact diagonalization of the Hamiltonian within the ν2, ν5 vibrational levels, and improved vibration-rotation constants for these bands are reported. The relative magnitudes and relative sings of polarizability derivatives involved in these vibrations are also reported.

Second-order cybernetics and enactive perception

Purpose - To present an account of cognition integrating second-order cybernetics (SOC) together with enactive perception and dynamic systems theory. Design/methodology/approach - The paper presents a brief critique of classical models of cognition then outlines how integration of SOC, enactive perception and dynamic systems theory can overcome some weaknesses of the classical paradigm. Findings - Presents the critique of evolutionary robotics showing how the issues of teleology and autonomy are left unresolved by this paradigm although their solution fits within the proposed framework. Research limitations/implications - The paper highlights the importance of genuine autonomy in the development of artificial cognitive systems. It sets out a framework within which the robofic research of cognitive systems could succeed. Practical implications - There are no immediate practical implications but see research implications. Originality/value - It joins the discussion on the fundamental nature of cognitive systems and emphasise the importance of autonomy and embodiment.

Second order accurate upwind difference schemes for scalar conservation laws with source terms

A second order accurate, characteristic-based, finite difference scheme is developed for scalar conservation laws with source terms. The scheme is an extension of well-known second order scalar schemes for homogeneous conservation laws. Such schemes have proved immensely powerful when applied to homogeneous systems of conservation laws using flux-difference splitting. Many application areas, however, involve inhomogeneous systems of conservation laws with source terms, and the scheme presented here is applied to such systems in a subsequent paper.

Enhanced second order non-linear optical properties of silica sol-gel thin films through photo-induced poling

Sol-gel derived inorganic materials are of interest as hosts for non-linear optically active guest molecules and they offer particular advantages in the field of non-linear optics. Orientationally ordered glasses have been prepared using a sol-gel system based on tetramethoxysilane, methyltrimethoxysilane and a non-linear optical chromophore Disperse Red 1. The novel technique of photo-induced poling was used to generate enhanced levels of polar order. The level of enhancement is strongly dependent on the extent of gelation and an optimum preparation time of ∼100 h led to an enhancement factor of ∼5. Films prepared in this manner exhibited a high stability of the polar order.

A comparison of photoinduced poling and thermal poling of azo‐dye‐doped polymer films for second order nonlinear optical applications

Photoinduced poling (PIP) is a new technique which allows the room‐temperature preparation of guest/host polymer films exhibiting significant polar order for nonlinear optical applications. We report a comparison of this novel technique with the conventional electrode poling procedure performed at the glass transition temperature of the polymer using disperse red 1/poly(methylmethacrylate) films. In particular, in situ second harmonic generation measurements show that levels of polar order achieved using these two techniques are similar. In contrast, the stability of the polar order is reduced by up to 20 times in terms of the decay time constant in films prepared using PIP although the stability is very dependent upon the temperature at which the poling was performed.

Boundary value problems for third-order linear PDEs in time-dependent domains

We present the extension of a methodology to solve moving boundary value problems from the second-order case to the case of the third-order linear evolution PDE qt + qxxx = 0. This extension is the crucial step needed to generalize this methodology to PDEs of arbitrary order. The methodology is based on the derivation of inversion formulae for a class of integral transforms that generalize the Fourier transform and on the analysis of the global relation associated with the PDE. The study of this relation and its inversion using the appropriate generalized transform are the main elements of the proof of our results.

A finite element method for nonlinear elliptic problems

We present a Galerkin method with piecewise polynomial continuous elements for fully nonlinear elliptic equations. A key tool is the discretization proposed in Lakkis and Pryer, 2011, allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretization method is that a recovered (finite element) Hessian is a byproduct of the solution process. We build on the linear method and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems as well as the Monge–Amp`ere equation and the Pucci equation.

Approximate Gauss-Newton methods for nonlinear least squares problems

The Gauss–Newton algorithm is an iterative method regularly used for solving nonlinear least squares problems. It is particularly well suited to the treatment of very large scale variational data assimilation problems that arise in atmosphere and ocean forecasting. The procedure consists of a sequence of linear least squares approximations to the nonlinear problem, each of which is solved by an “inner” direct or iterative process. In comparison with Newton’s method and its variants, the algorithm is attractive because it does not require the evaluation of second-order derivatives in the Hessian of the objective function. In practice the exact Gauss–Newton method is too expensive to apply operationally in meteorological forecasting, and various approximations are made in order to reduce computational costs and to solve the problems in real time. Here we investigate the effects on the convergence of the Gauss–Newton method of two types of approximation used commonly in data assimilation. First, we examine “truncated” Gauss–Newton methods where the inner linear least squares problem is not solved exactly, and second, we examine “perturbed” Gauss–Newton methods where the true linearized inner problem is approximated by a simplified, or perturbed, linear least squares problem. We give conditions ensuring that the truncated and perturbed Gauss–Newton methods converge and also derive rates of convergence for the iterations. The results are illustrated by a simple numerical example. A practical application to the problem of data assimilation in a typical meteorological system is presented.

Application of the Lighthill-Ford theory of spontaneous imbalance to clear-air turbulence forecasting

A new method of clear-air turbulence (CAT) forecasting based on the Lighthill–Ford theory of spontaneous imbalance and emission of inertia–gravity waves has been derived and applied on episodic and seasonal time scales. A scale analysis of this shallow-water theory for midlatitude synoptic-scale flows identifies advection of relative vorticity as the leading-order source term. Examination of leading- and second-order terms elucidates previous, more empirically inspired CAT forecast diagnostics. Application of the Lighthill–Ford theory to the Upper Mississippi and Ohio Valleys CAT outbreak of 9 March 2006 results in good agreement with pilot reports of turbulence. Application of Lighthill–Ford theory to CAT forecasting for the 3 November 2005–26 March 2006 period using 1-h forecasts of the Rapid Update Cycle (RUC) 2 1500 UTC model run leads to superior forecasts compared to the current operational version of the Graphical Turbulence Guidance (GTG1) algorithm, the most skillful operational CAT forecasting method in existence. The results suggest that major improvements in CAT forecasting could result if the methods presented herein become operational.