A Variational Assimilation Method for Satellite and Conventional Data: a Revised Basic Model 2B

A variational objective analysis technique that modifies observations of temperature, height, and wind on the cyclone scale to satisfy the five 'primitive' model forecast equations is presented. This analysis method overcomes all of the problems that hindered previous versions, such as over-determination, time consistency, solution method, and constraint decoupling. A preliminary evaluation of the method shows that it converges rapidly, the divergent part of the wind is strongly coupled in the solution, fields of height and temperature are well-preserved, and derivative quantities such as vorticity and divergence are improved. Problem areas are systematic increases in the horizontal velocity components, and large magnitudes of the local tendencies of the horizontal velocity components. The preliminary evaluation makes note of these problems but detailed evaluations required to determine the origin of these problems await future research

Macro-element interpolation on tensor product meshes

A general theory for obtaining anisotropic interpolation error estimates for macro-element interpolation is developed revealing general construction principles. We apply this theory to interpolation operators on a macro type of biquadratic $C^1$ finite elements on rectangle grids which can be viewed as a rectangular version of the $C^1$ Powell-Sabin element. This theory also shows how interpolation on the Bogner-Fox-Schmidt finite element space (or higher order generalizations) can be analyzed in a unified framework. Moreover we discuss a modification of Scott-Zhang type giving optimal error estimates under the regularity required without imposing quasi uniformity on the family of macro-element meshes used. We introduce and analyze an anisotropic macro-element interpolation operator, which is the tensor product of one-dimensional $C^1-P_2$ macro interpolation and $P_2$ Lagrange interpolation. These results are used to approximate the solution of a singularly perturbed reaction-diffusion problem on a Shishkin mesh that features highly anisotropic elements. Hereby we obtain an approximation whose normal derivative is continuous along certain edges of the mesh, enabling a more sophisticated analysis of a continuous interior penalty method in another paper

A multivariate variational objective analysis-assimilation method. Part 1: Development of the basic model

The variational method of undetermined multipliers is used to derive a multivariate model for objective analysis. The model is intended for the assimilation of 3-D fields of rawinsonde height, temperature and wind, and mean level temperature observed by satellite into a dynamically consistent data set. Relative measurement errors are taken into account. The dynamic equations are the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation. The model Euler-Lagrange equations are eleven linear and/or nonlinear partial differential and/or algebraic equations. A cyclical solution sequence is described. Other model features include a nonlinear terrain-following vertical coordinate that eliminates truncation error in the pressure gradient terms of the horizontal momentum equations and easily accommodates satellite observed mean layer temperatures in the middle and upper troposphere. A projection of the pressure gradient onto equivalent pressure surfaces removes most of the adverse impacts of the lower coordinate surface on the variational adjustment

Locking-free two-layer Timoshenko beam element with interlayer slip

A new locking-free strain-based finite element formulation for the numerical treatment of linear static analysis of two-layer planar composite beams with interlayer slip is proposed. In this formulation, the modified principle of virtual work is introduced as a basis for the finite element discretization. The linear kinematic equations are included into the principle by the procedure, similar to that of Lagrangian multipliers. A strain field vector remains the only unknown function to be interpolated in the finite element implementation of the principle. In contrast with some of the displacement-based and mixed finite element formulations of the composite beams with interlayer slip, the present formulation is completely locking-free. Hence, there are no shear and slip locking, poor convergence and stress oscillations in these finite elements. The generalization of the composite beam theory with the consideration of the Timoshenko beam theory for the individual component of a composite beam represents a substantial contribution in the field of analysis of non-slender composite beams with an interlayer slip. An extension of the present formulation to the non-linear material problems is straightforward. As only a few finite elements are needed to describe a composite beam with great precision, the new finite element formulations is perfectly suited for practical calculations. (c) 2007 Elsevier B.V. All rights reserved

Spectral-Element and Adjoint Methods in Seismology

We provide an introduction to the use of the spectral-element method (SEM) in seismology. Following a brief review of the basic equations that govern seismic wave propagation, we discuss in some detail how these equations may be solved numerically based upon the SEM to address the forward problem in seismology. Examples of synthetic seismograms calculated based upon the SEM are compared to data recorded by the Global Seismographic Network. Finally, we discuss the challenge of using the remaining differences between the data and the synthetic seismograms to constrain better Earth models and source descriptions. This leads naturally to adjoint methods, which provide a practical approach to this formidable computational challenge and enables seismologists to tackle the inverse problem

Nonlinear nonlocal multicontinua upscaling framework and its applications

In this paper, we discuss multiscale methods for nonlinear problems. The main idea of these approaches is to use local constraints and solve problems in oversampled regions for constructing macroscopic equations. These techniques are intended for problems without scale separation and high contrast, which often occur in applications. For linear problems, the local solutions with constraints are used as basis functions. This technique is called Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM). GMsFEM identifies macroscopic quantities based on rigorous analysis. In corresponding upscaling methods, the multiscale basis functions are selected such that the degrees of freedom have physical meanings, such as averages of the solution on each continuum. This paper extends the linear concepts to nonlinear problems, where the local problems are nonlinear. The main concept consists of: (1) identifying macroscopic quantities; (2) constructing appropriate oversampled local problems with coarse-grid constraints; (3) formulating macroscopic equations. We consider two types of approaches. In the first approach, the solutions of local problems are used as basis functions (in a linear fashion) to solve nonlinear problems. This approach is simple to implement; however, it lacks the nonlinear interpolation, which we present in our second approach. In this approach, the local solutions are used as a nonlinear forward map from local averages (constraints) of the solution in oversampling region. This local fine-grid solution is further used to formulate the coarse-grid problem. Both approaches are discussed on several examples and applied to single-phase and two-phase flow problems, which are challenging because of convection-dominated nature of the concentration equation

Phase-field approach to polycrystalline solidification including heterogeneous and homogeneous nucleation

Advanced phase-field techniques have been applied to address various aspects of polycrystalline solidification including different modes of crystal nucleation. The height of the nucleation barrier has been determined by solving the appropriate Euler-Lagrange equations. The examples shown include the comparison of various models of homogeneous crystal nucleation with atomistic simulations for the single component hard-sphere fluid. Extending previous work for pure systems (Gránásy L, Pusztai T, Saylor D and Warren J A 2007 Phys. Rev. Lett. 98 art no 035703), heterogeneous nucleation in unary and binary systems is described via introducing boundary conditions that realize the desired contact angle. A quaternion representation of crystallographic orientation of the individual particles (outlined in Pusztai T, Bortel G and Gránásy L 2005 Europhys. Lett. 71 131) has been applied for modeling a broad variety of polycrystalline structures including crystal sheaves, spherulites and those built of crystals with dendritic, cubic, rhombododecahedral, truncated octahedral growth morphologies. Finally, we present illustrative results for dendritic polycrystalline solidification obtained using an atomistic phase-field model

A Variational Assimilation Method for Satellite and Conventional Data: Development of Basic Model for Diagnosis of Cyclone Systems

A summary is presented of the progress toward the completion of a comprehensive diagnostic objective analysis system based upon the calculus of variations. The approach was to first develop the objective analysis subject to the constraints that the final product satisfies the five basic primitive equations for a dry inviscid atmosphere: the two nonlinear horizontal momentum equations, the continuity equation, the hydrostatic equation, and the thermodynamic equation. Then, having derived the basic model, there would be added to it the equations for moist atmospheric processes and the radiative transfer equation