Discrete Second Order Adjoints in Atmospheric Chemical Transport Modeling

Atmospheric chemical transport models (CTMs) are essential tools for the study of air pollution, for environmental policy decisions, for the interpretation of observational data, and for producing air quality forecasts. Many air quality studies require sensitivity analyses, i.e., the computation of derivatives of the model output with respect to model parameters. The derivatives of a cost functional (defined on the model output) with respect to a large number of model parameters can be calculated efficiently through adjoint sensitivity analysis. While the traditional (first order) adjoint models give the gradient of the cost functional with respect to parameters, second order adjoint models give second derivative information in the form of products between the Hessian of the cost functional and a user defined vector. In this paper we discuss the mathematical foundations of the discrete second order adjoint sensitivity method and present a complete set of computational tools for performing second order sensitivity studies in three-dimensional atmospheric CTMs. The tools include discrete second order adjoints of Runge Kutta and of Rosenbrock time stepping methods for stiff equations together with efficient implementation strategies. Numerical examples illustrate the use of these computational tools in important applications like sensitivity analysis, optimization, uncertainty quantification, and the calculation of directions of maximal error growth in three-dimensional atmospheric CTMs

Research and Technology, 1989

Selected research and technology activities at Ames Research Center, including the Moffett Field site and the Dryden Flight Research Facility, are summarized. These accomplishments exemplify the Center's varied and highly productive research efforts for 1989

Large eddy simulation of buoyant plumes

A 3D parallel CFD code is written to investigate the characteristics of and differences between Large Eddy Simulation (LES) models in the context of simulating a thermal buoyant plume. An efficient multigrid scheme is incorporated to solve the Poisson equation, resulting from the fractional step, projection method used to solve the Low Mach Number (LMN) Navier-Stokes equations. A wide range of LES models are implemented, including a variety of eddy models, structure models, mixed models and dynamic models, for both the momentum stresses and the temperature fluxes. Generalised gradient flux models are adapted from their RANS counterparts, and also tested. A number of characteristics are observed in the LES models relating to the thermal plume simulation in particular and turbulence in general. Effects on transition, dissipation, backscatter, equation balances, intermittency and energy spectra are all considered, as are the impact of the governing equations, the discretisation scheme, and the effect of grid coarsening. Also characteristics to particular models are considered, including the subgrid kinetic energy for the one-equation models, and constant histories for dynamic models. The argument that choice of LES model is unimportant is shown to be incorrect as a general statement, and a recommendation for when the models are best used is given.EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Large-scale 3D environmental modelling and visualisation for flood hazard warning.

3D environment reconstruction has received great interest in recent years in areas such as city planning, virtual tourism and flood hazard warning. With the rapid development of computer technologies, it has become possible and necessary to develop new methodologies and techniques for real time simulation for virtual environments applications. This thesis proposes a novel dynamic simulation scheme for flood hazard warning. The work consists of three main parts: digital terrain modelling; 3D environmental reconstruction and system development; flood simulation models. The digital terrain model is constructed using real world measurement data of GIS, in terms of digital elevation data and satellite image data. An NTSP algorithm is proposed for very large data assessing, terrain modelling and visualisation. A pyramidal data arrangement structure is used for dealing with the requirements of terrain details with different resolutions. The 3D environmental reconstruction system is made up of environmental image segmentation for object identification, a new shape match method and an intelligent reconstruction system. The active contours-based multi-resolution vector-valued framework and the multi-seed region growing method are both used for extracting necessary objects from images. The shape match method is used with a template in the spatial domain for a 3D detailed small scale urban environment reconstruction. The intelligent reconstruction system is designed to recreate the whole model based on specific features of objects for large scale environment reconstruction. This study then proposes a new flood simulation scheme which is an important application of the 3D environmental reconstruction system. Two new flooding models have been developed. The first one is flood spreading model which is useful for large scale flood simulation. It consists of flooding image spatial segmentation, a water level calculation process, a standard gradient descent method for energy minimization, a flood region search and a merge process. The finite volume hydrodynamic model is built from shallow water equations which is useful for urban area flood simulation. The proposed 3D urban environment reconstruction system was tested on our simulation platform. The experiment results indicate that this method is capable of dealing with complicated and high resolution region reconstruction which is useful for many applications. When testing the 3D flood simulation system, the simulation results are very close to the real flood situation, and this method has faster speed and greater accuracy of simulating the inundation area in comparison to the conventional flood simulation model

Research in Natural Laminar Flow and Laminar-Flow Control, part 2

Part 2 of the Symposium proceedings includes papers addressing various topics in basic wind tunnel research/techniques and computational transitional research. Specific topics include: advanced measurement techniques; laminar flow control; Tollmien-Schlichting wave characteristics; boundary layer transition; flow visualization; wind tunnel tests; flight tests; boundary layer equations; swept wings; and skin friction

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described