Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

A long-term space astrophysics research program: The evolution of the quasar continuum

The research program supported by this grant now has great momentum. Numerous papers are in progress, and a strong multi-wavelength observing program is rapidly accumulating data on samples of high redshift quasars across the spectrum. ROSAT spectra of quasars continue to yield surprises. Of four z = 3 quasars with X-ray spectra, three show strong absorption. This contrasts strongly with the situation for luminous AGN at low redshifts where fewer than 1 in 20 show X-ray absorption. A new site for this absorption is probably needed, either around the quasar (e.g. in a cluster cooling flow) or along the line of sight (e.g. in a Damped Lyman-alpha system). The unabsorbed quasar allows limits on the physical conditions in a damped Lyman-alpha cloud to be calculated, and will allow a X-ray Gunn-Peterson test to be applied that will limit the fraction of the closure mass in an intergalactic medium. The X-ray spectral indices of these z = 3 quasars show no change from those of similar objects at low z, suggesting that 'short-lifetime' models apply. Eight other z = 3-4 quasars have been detected and their energy distributions from X-rays to Infrared (using new infrared spectrographs) have been compiled. These are now being compared with the low z continua from the 'Atlas of Quasar Energy Distributions' to search for evolutionary changes. The discovery of a likely warm absorber in 3C351 made recognition of another example simple. Also, modeling of the conditions in the absorber in 3C351 using the OVI absorption line from HST and the high ionization emission lines, suggests that the broad line region is indeed the origin of the warm absorber in this quasar, and by extension, others. Warm absorbers can now be used as a new diagnostic of this region. The X-ray spectrum of a 'Red Quasar', 3C212, has a cut-off spectrum, which could be fitted by an absorbed power-law, or more remarkably, by an unabsorbed black body. Using our quasi-simultaneous optical data and photoionization modeling we are able to rule out this latter possibility. The long-sought key to understanding red quasars has now been found. (Reddening was earlier ruled out because no 2175A absorption feature was seen. This is now seen to have been misleading.)

Elvis, M.↗

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗

The XXL Survey I. Scientific Motivations - Xmm-Newton Observing Plan - Follow-up Observations and Simulation Programme

The quest for the cosmological parameters that describe our universe continues to motivate the scientific community to undertake very large survey initiatives across the electromagnetic spectrum. Over the past two decades, the Chandra and XMM-Newton observatories have supported numerous studies of X-ray-selected clusters of galaxies, active galactic nuclei (AGNs), and the X-ray background. The present paper is the first in a series reporting results of the XXL-XMM survey; it comes at a time when the Planck mission results are being finalized. Aims. We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of approx. 5 x 10(exp 15) erg/s/sq cm in the [0.5-2] keV band (completeness limit). The surveys main goals are to provide constraints on the dark energy equation of state from the space-time-distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programs. Methods. We describe the 542 XMM observations along with the associated multi- and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. Results. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. Conclusions. The XXL multi- data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve asa calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps.

X-rays: general – large-scale structure of Unive↗

A modeling study of the time-averaged electric currents in the vicinity of isolated thunderstorms

A thorough examination of the results of a time-dependent computer model of a dipole thunderstorm revealed that there are numerous similarities between the time-averaged electrical properties and the steady-state properties of an active thunderstorm. Thus, the electrical behavior of the atmosphere in the vicinity of a thunderstorm can be determined with a formulation similar to what was first described by Holzer and Saxon (1952). From the Maxwell continuity equation of electric current, a simple analytical equation was derived that expresses a thunderstorm's average current contribution to the global electric circuit in terms of the generator current within the thundercloud, the intracloud lightning current, the cloud-to-ground lightning current, the altitudes of the charge centers, and the conductivity profile of the atmosphere. This equation was found to be nearly as accurate as the more computationally expensive numerical model, even when it is applied to a thunderstorm with a reduced conductivity thundercloud, a time-varying generator current, a varying flash rate, and a changing lightning mix.

Driscoll, Kevin T.↗

Theory of photoelectron thermalization and transport in the ionosphere

A new theoretical approach for calculating the equilibrium photoelectron flux energy and pitch angle distribution in the ionosphere is presented. Photoelectron transport, secondary electron production, and energy degradation by the excitation of the discrete energy states of the neutral atmospheric constituents and by continuous energy transfer to the ambient thermal electron gas are included in a manner consistent with the Boltzmann equation which constitutes the foundation of the theory. A difference equation, suitable for numerical solution, is given, and a numerical method for the solution of this equation is discussed in detail.

Mantas, G. P.↗

Particle erosion mechanisms and mass redistribution in Saturn's rings

A variety of physical processes can erode the surfaces of planetary ring particles. According to current estimates, the most efficient of these over the bulk of Saturn's rings is hypervelocity impact by 100 microns to one centimeter radius meteoroids. The atoms, molecules, and fragments ejected from ring particles by erosion arc across the rings along elliptical orbits to produce a tenuous halo of solid ejecta and an extensive gaseous atmosphere. Continuous exchange of ejecta between different ring regions can lead to net radial transport of mass and angular momentum. The equations governing this ballistic transport process are presented and discussed. Both numerical and analytic studies of idealized ring systems illustrate that ballistic transport can cause significant mass redistribution in the rings, especially near regions of high density contrast, such as the inner edges of the A and B rings. Ejecta exchanges can also alter local particle sizes and compositions and may produce pulverized regoliths at least several centimeters deep. The meteoroid erosion rate is so high that significant global torques and mass loss are possible on times shorter than a solar system life time.

Durisen, R. H.↗

A multidimensional model for aerosols - Description of computational analogs

The numerical algorithms used to simulate the advection, diffusion, sedimentation, coagulation, and condensational growth of atmospheric aerosols are described. The model can be used in one, two, or three spatial dimensions. The continuity equation in a generalized horizontal and vertical coordinate system is developed, which allows the model to be quickly adapted to a wide variety of dynamical models of global or regional scale. Algorithms are developed to treat the various physical processes, and the results of simulations are presented, which show the strengths and weaknesses of these algorithms. Although the emphasis is on the modeling of aerosols, the work is also applicable to the simulations of the transport of gases.

Toon, O. B.↗

Probabilistic Structural Analysis Methods (PSAM) for select space propulsion system components

This annual report summarizes the work completed during the third year of technical effort on the referenced contract. Principal developments continue to focus on the Probabilistic Finite Element Method (PFEM) which has been under development for three years. Essentially all of the linear capabilities within the PFEM code are in place. Major progress in the application or verifications phase was achieved. An EXPERT module architecture was designed and partially implemented. EXPERT is a user interface module which incorporates an expert system shell for the implementation of a rule-based interface utilizing the experience and expertise of the user community. The Fast Probability Integration (FPI) Algorithm continues to demonstrate outstanding performance characteristics for the integration of probability density functions for multiple variables. Additionally, an enhanced Monte Carlo simulation algorithm was developed and demonstrated for a variety of numerical strategies.

Source record↗

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

Oceanic Observations

For many years, merchant ships and the naval fleets of various countries have been the major source of data over and in the open ocean. Oceanographic research experiments and process studies in the field have also contributed to the climatological data bases for the global ocean, but, for the most part, these have been limited in duration and extent. However, over the last 10 years under the auspices of the World Climate Research Program and the International Geosphere Biosphere Program the role of the oceans in global and climate change has taken on increased significance. This has created a need for a considerably improved understanding of the seasonal, interannual, decadal and longer time-scale variability of the physical and biogeochemical attributes of the global ocean. As a result, over the past 10 years several major international field programs have been implemented and have had a tremendous impact on the number of in situ observations obtained for the global ocean. The Tropical Ocean Global Atmosphere (TOGA) program, the World Ocean Circulation Experiment (WOCE), and the Joint Global Ocean Flux Study (JGOFS) were designed with observational, modelling, and process study components aimed at analyzing different aspects of the ocean's role in the coupled climate system. In parallel with the field programs, continuous space-based observations of sea surface temperature, sea surface topography, and sea surface winds spanning nearly a decade or longer have become a reality. During this same time period, numerical ocean models and computational power have advanced to the point where the oceanographic observations, both in situ and remotely sensed, can be assimilated into numerical ocean models in order to provide a four-dimensional (x-y-z-t) depiction of the evolving state of the global ocean.

Busalacchi, Antonio J.↗

Impact into Coarse Grained Spheres

Several experimental studies [1,2,3] indicate that differences in the grain size of the target relative to the projectile could influence the cratering process. Impacts into coarse sand grains of size comparable to the projectile show some discrepancies with existing relationships for crater growth [e.g. 4]. Similarly, targets of ne grained, uniform in diameter glass spheres show differences in crater depth, transient crater diameter, and volume of ejecta excavated as a function of grain size [2,3]. The purpose of this work is to continue investigating how the relative grain size may influence early time coupling between a projectile and target, with implications for subsequent ejecta excavation and crater growth. In previous efforts we used numerical techniques to focus on the propagation of shock waves in coarse, granular media emphasizing the influence of relative grain size on crater growth, ejecta production, cratering efficiency, target strength, and crater shape [5,6,7]. In this study, we use experimental techniques - in part as a reality check for the numerical studies - to report on how coarse grained targets might influence ejecta excavation and crater shape. This body of work possesses important implications for ejecta excavation and cratering efficiency on asteroids that may possess rubble pile-like structures, and on planets that may possess either pre-fractured surfaces or large-scale heterogeneities in shock impedance.

Barnouin-Jha, O. S.↗

Transport and Use of a Centaur Second Stage in Space

As nations continue to explore space, the desire to reduce costs will continue to grow. As a method of cost reduction, transporting and/or use of launch system components as integral components of missions may become more commonplace in the future. There have been numerous scenarios written for using launch vehicle components (primarily space shuttle used external tanks) as part of flight missions or future habitats. Future studies for possible uses of launch vehicle upper stages might include asteroid diverter using gravity orbital perturbation, orbiting station component, raw material at an outpost, and kinetic impactor. The LCROSS (Lunar CRater Observation and Sensing Satellite) mission was conceived as a low-cost means of determining whether water exists at the polar regions of the moon. Manifested as a secondary payload with the LRO (Lunar Reconnaissance Orbiter) spacecraft aboard an Atlas V launch vehicle, LCROSS guided its spent Centaur Earth Departure Upper Stage (EDUS) into the lunar crater Cabeu's, as a kinetic impactor. This paper describes some of the challenges that the LCROSS project encountered in planning, designing, launching with and carrying the Centaur upper stage to the moon.

Strong, James M.↗

Transmission Electron Microscopy of Itokawa Regolith Grains

Introduction: In a remarkable engineering achievement, the JAXA space agency successfully recovered the Hayabusa space-craft in June 2010, following a non-optimal encounter and sur-face sampling mission to asteroid 25143 Itokawa. These are the first direct samples ever obtained and returned from the surface of an asteroid. The Hayabusa samples thus present a special op-portunity to directly investigate the evolution of asteroidal sur-faces, from the development of the regolith to the study of the effects of space weathering. Here we report on our preliminary TEM measurements on two Itokawa samples. Methods: We were allocated particles RA-QD02-0125 and RA-QD02-0211. Both particles were embedded in low viscosity epoxy and thin sections were prepared using ultramicrotomy. High resolution images and electron diffraction data were ob-tained using a JEOL 2500SE 200 kV field-emission scanning-transmission electron microscope. Quantitative maps and anal-yses were obtained using a Thermo thin-window energy-dispersive x-ray (EDX) spectrometer. Results: Both particles are olivine-rich (Fo70) with μm-sized inclusions of FeS and have microstructurally complex rims. Par-ticle RA-QD02-0125 is rounded and has numerous sub-μm grains attached to its surface including FeS, albite, olivine, and rare melt droplets. Solar flare tracks have not been observed, but the particle is surrounded by a continuous ~50 nm thick, stuctur-ally disordered rim that is compositionally similar to the core of the grain. One of the surface adhering grains is pyrrhotite show-ing a S-depleted rim (~8-10 nm thick) with nanophase Fe metal grains (<5 nm) decorating the outermost surface. The pyrrhotite displays a complex superstructure in its core that is absent in the S-depleted rim. Particle RA-QD02-0211 contains solar flare particle tracks (~2x109 cm-2) and shows a structurally disordered rim ~100 nm thick. The track density corresponds to a surface exposure of ~103-104 years based on the track production rate of [1]. The dis-ordered rim is nanocrystalline with minor amorphous material between crystalline domains. Quantitative element maps show the outermost ~10 nm of the disordered rim is Si-rich. Discussion and Conclusions: Both particles record the ef-fects of space weathering processes on Itokawa. Noguchi et al. [2] proposed that the disordered rims they observed on Itokawa particles largely result from solar wind radiation damage and we arrive at a similar conclusion for the two particles we analyzed. The microstructure of the S-depleted layer on the pyrrhotite grain in RA-QD02-0125 is similar to that observed in troilite irradiated with ~1018 4 kV He+ [3, 4]. Prolonged irradiation has also been shown to disorder pyrrhotite such that the superstructure reflec-tions are lost [5].

Keller, Lindsay P.↗

The Aqua-planet Experiment (APE): Response to Changed Meridional SST Profile

This paper explores the sensitivity of Atmospheric General Circulation Model (AGCM) simulations to changes in the meridional distribution of sea surface temperature (SST). The simulations are for an aqua-planet, a water covered Earth with no land, orography or sea- ice and with specified zonally symmetric SST. Simulations from 14 AGCMs developed for Numerical Weather Prediction and climate applications are compared. Four experiments are performed to study the sensitivity to the meridional SST profile. These profiles range from one in which the SST gradient continues to the equator to one which is flat approaching the equator, all with the same maximum SST at the equator. The zonal mean circulation of all models shows strong sensitivity to latitudinal distribution of SST. The Hadley circulation weakens and shifts poleward as the SST profile flattens in the tropics. One question of interest is the formation of a double versus a single ITCZ. There is a large variation between models of the strength of the ITCZ and where in the SST experiment sequence they transition from a single to double ITCZ. The SST profiles are defined such that as the equatorial SST gradient flattens, the maximum gradient increases and moves poleward. This leads to a weakening of the mid-latitude jet accompanied by a poleward shift of the jet core. Also considered are tropical wave activity and tropical precipitation frequency distributions. The details of each vary greatly between models, both with a given SST and in the response to the change in SST. One additional experiment is included to examine the sensitivity to an off-equatorial SST maximum. The upward branch of the Hadley circulation follows the SST maximum off the equator. The models that form a single precipitation maximum when the maximum SST is on the equator shift the precipitation maximum off equator and keep it centered over the SST maximum. Those that form a double with minimum on the equatorial maximum SST shift the double structure off the equator, keeping the minimum over the maximum SST. In both situations only modest changes appear in the shifted profile of zonal average precipitation. When the upward branch of the Hadley circulation moves into the hemisphere with SST maximum, the zonal average zonal, meridional and vertical winds all indicate that the Hadley cell in the other hemisphere dominates.

AGCM↗

On Promoting the Use of Lidar Systems in Forest Ecosystem Research

Forest structure is an important driver of ecosystem dynamics, including the exchange of carbon, water and energy between canopies and the atmosphere. Structural descriptors are also used in numerous studies of eco-logical processes and ecosystem services. Over the last 20+ years, lidar technology has fundamentally changed the way we observe and describe forest structure, and it will continue to impact the ways in which we investigate and monitor the relations between forest structure and functions. Here we present the currently available lidar system types (ground, air, and space-based), we highlight opportunities and challenges associated with each system, as well as challenges associated with a wider use of lidar technology and wider availability of lidar derived products. We also suggest pathways for lidar to further contribute to addressing questions in forest ecosystem science and increase benefits to a wider community of researchers.

Lidar systems↗

The solution of three-variable duct-flow equations

This paper establishes a numerical method for the solution of three-variable problems and is applied here to rotational flows through ducts of various cross sections. An iterative scheme is developed, the main feature of which is the addition of a duplicate variable to the forward component of velocity. Two forward components of velocity result from integrating two sets of first order ordinary differential equations for the streamline curvatures, in intersecting directions across the duct. Two pseudo-continuity equations are introduced with source/sink terms, whose strengths are dependent on the difference between the forward components of velocity. When convergence is obtained, the two forward components of velocity are identical, the source/sink terms are zero, and the original equations are satisfied. A computer program solves the exact equations and boundary conditions numerically. The method is economical and compares successfully with experiments on bent ducts of circular and rectangular cross section where secondary flows are caused by gradients of total pressure upstream.

Stuart, A. R.↗

2008 NASA Range Safety Annual Report

Welcome to the 2008 edition of the NASA Range Safety Annual Report. Funded by NASA Headquarters, this report provides a NASA Range Safety overview for current and potential range users. This year, along with full length articles concerning various subject areas, we have provided updates to standard subjects with links back to the 2007 original article. Additionally, we present summaries from the various NASA Range Safety Program activities that took place throughout the year, as well as information on several special projects that may have a profound impact on the way we will do business in the future. The sections include a program overview and 2008 highlights of Range Safety Training; Range Safety Policy; Independent Assessments and Common Risk Analysis Tools Development; Support to Program Operations at all ranges conducting NASA launch operations; a continuing overview of emerging Range Safety-related technologies; Special Interests Items that include recent changes in the ELV Payload Safety Program and the VAS explosive siting study; and status reports from all of the NASA Centers that have Range Safety responsibilities. As is the case each year, contributors to this report are too numerous to mention, but we thank individuals from the NASA Centers, the Department of Defense, and civilian organizations for their contributions. We have made a great effort to include the most current information available. We recommend that this report be used only for guidance and that the validity and accuracy of all articles be verified for updates. This is the third year we have utilized this web-based format for the annual report. We continually receive positive feedback on the web-based edition, and we hope you enjoy this year's product as well. It has been a very busy and productive year on many fronts as you will note as you review this report. Thank you to everyone who contributed to make this year a successful one, and I look forward to working with all of you in the years to come.

Lamoreaux, Richard W.↗