Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “singular characteristics”

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 181 records · Page 10

Inflection point caustic problems and solutions for high-gain dual-shaped reflectors

The singular nature of the uniform geometrical theory of diffraction (UTD) subreflector scattered field at the vicinity of the main reflector edge (for a high-gain antenna design) is investigated. It is shown that the singularity in the UTD edge-diffracted and slope-diffracted fields is due to the reflection distance parameter approaching infinity in the transition functions. While the geometrical optics (GO) and UTD edge-diffracted fields exhibit singularities of the same order, the edge slope-diffracted field singularity is more significant and is substantial for greater subreflector edge tapers. The diffraction analysis of such a subreflector in the vicinity of the main reflector edge has been carried out efficiently and accurately by a stationary phase evaluation of the phi-integral, whereas the theta-integral is carried out numerically. Computational results from UTD and physical optics (PO) analysis of a 34-m ground station dual-shaped reflector confirm the analytical formulations for both circularly symmetric and offset asymmetric subreflectors. It is concluded that the proposed PO(theta)GO(phi) technique can be used to study the spillover or noise temperature characteristics of a high-gain reflector antenna efficiently and accurately.

Galindo-Israel, Victor↗

Chapter 12 - Flight Envelope

The term "flight envelope" is used to refer to the boundaries of aircraft loading and flight conditions within which operation of the aircraft is satisfactory, and beyond which some aspect becomes unacceptable. This flight envelope represents, in fact, the limiting conditions arising from a matrix of inter-related flight envelopes covering the appropriate variables. Thus, for each loading (i.e., external stores configuration and its associated range of weight and center of gravity (c.g.) position) and aircraft configuration (i.e., position of undercarriage (u/c), flaps, slats, etc.), the envelopes of airspeed versus altitude, airspeed versus load factor, angle of attack versus angle of sideslip, etc., must be investigated to establish the limits within which all aspects such as handling qualities, engine behavior, structural loads, etc., remain acceptable. Flight testing of new or derivative aircraft models is carried out with the initial purpose of defining a flight envelope which is, first and foremost, safe and secondarily, which enables the effective use of the vehicle for its intended purpose. Flight testing occurs only after numerous reviews of the design and review of results from ground tests and predictions of flight characteristics in such areas as structures, aerodynamics, stability and control, flight controls (particularly fly-by-wire control systems, propulsion, etc.). Accordingly, opening and expanding the envelope is a task that must be approached cautiously, systematically, and with coordination and cooperation of the many disciplines involved in the design and test of an airplane. (Sections 8 and 10 cover test planning and safety of flight considerations, respectively). The fundamental tenet in establishing a flight envelope via flight test is risk reduction. This is reflected in the typical sequence of events leading to initial flight test - design reviews (both hardware and software), then ground test involving singular disciplines (windtunnel tests for aerodynamics, structurally loading the wing/fuselage/nacelle on a ground test article with loads anticipated to occur in flight, flight control system control law checkout, propulsion test cell runs and/or flying test bed tests, etc.), and then ground tests involving multi-disciplines (See Section 9). Only after these have been accomplished will an initial, limited, low-risk, flight envelope be established. The limited envelope will typically be in the middle of the projected final flight envelope. Subsequent flight tests will then be devoted to expanding the initial envelope by operating the airplane at increasing ranges - representing increasing risk - of engine operation, airspeeds both fast and slow, altitude, load factor both above and below 1g, centers of gravity (fore and aft), and with system/subsystem failures. Whether flight tests are to define a flight envelope on a new model airplane with the attendant new airframe, new engine(s), and new subsystems (hydraulics, pressurization, etc.), or on an airplane involving only a few of these areas such as new engines in an old airframe, the fundamental approach to establishing an envelope is the same.

H Walgemoed↗

Refinement of a limit cycle oscillator model of the effects of light on the human circadian pacemaker

In 1990, Kronauer proposed a mathematical model of the effects of light on the human circadian pacemaker. Although this model predicted many general features of the response of the human circadian pacemaker to light exposure, additional data now available enable us to refine the original model. We first refined the original model by incorporating the results of a dose response curve to light into the model's predicted relationship between light intensity and the strength of the drive onto the pacemaker. Data from three bright light phase resetting experiments were then used to refine the amplitude recovery characteristics of the model. Finally, the model was tested and further refined using data from an extensive phase resetting experiment in which a 3-cycle bright light stimulus was presented against a background of dim light. In order to describe the results of the four resetting experiments, the following major refinements to the original model were necessary: (i) the relationship between light intensity (I) and drive onto the pacemaker was reduced from I1/3 to I0.23 for light levels between 150 and 10,000 lux; (ii) the van der Pol oscillator from the original model was replaced with a higher-order limit cycle oscillator so that amplitude recovery is slower near the singularity and faster near the limit cycle; (iii) a direct effect of light on circadian period (tau x) was incorporated into the model such that as I increases, tau x decreases, which is in accordance with "Aschoff's rule". This refined model generates the following testable predictions: it should be difficult to enhance normal circadian amplitude via bright light; near the critical point of a type 0 phase response curve (PRC) the slope should be steeper than it is in a type 1 PRC; and circadian period measured during forced desynchrony should be directly affected by ambient light intensity.

NASA Discipline Regulatory Physiology↗

A flexible class of priors for orthonormal matrices with basis function-specific structure

Statistical modeling of high-dimensional matrix-valued data motivates the use of a low-rank representation that simultaneously summarizes key characteristics of the data and enables dimension reduction. Low-rank representations commonly factor the original data into the product of orthonormal basis functions and weights, where each basis function represents an independent feature of the data. However, the basis functions in these factorizations are typically computed using algorithmic methods that cannot quantify uncertainty or account for basis function correlation structure a priori. While there exist Bayesian methods that allow for a common correlation structure across basis functions, empirical examples motivate the need for basis function-specific dependence structure. We propose a prior distribution for orthonormal matrices that can explicitly model basis function-specific structure. The prior is used within a general probabilistic model for singular value decomposition to conduct posterior inference on the basis functions while accounting for measurement error and fixed effects. We discuss how the prior specification can be used for various scenarios and demonstrate favorable model properties through synthetic data examples. Finally, we apply our method to two-meter air temperature data from the Pacific Northwest, enhancing our understanding of the Earth system’s internal variability.

97 MATHEMATICS AND COMPUTING↗

Greenhouse gas emissions reduction strategies that maximize portfolio-wide life cycle cost reduction, resilience, and environmental justice benefits

While strategies to achieve net-zero emissions at an individual site are well understood, new analysis methods are required for organizations seeking to achieve net-zero across multiple facilities, each with concurrent priority goals. At a portfolio level, distinct locations present varied challenges that cannot be addressed through singular solutions, and competing goals can take precedence with the assumption that net-zero emissions strategies deter from energy resilience and cost savings, therefore negatively impacting nearby communities. This study tests these assumptions by analyzing 16 diverse sites (varying in size, climate, and energy use) to identify strategies that reduce emissions and assess the impact these strategies have on life cycle costs, resilience, and communities with environmental justice concerns. Methods were developed to approximate missing information essential to net-zero evaluation. Established methods were augmented to evaluate life cycle costs, resilience, and environmental justice impacts across a set of strategies and accommodate the multi-criteria analyses. Potential benefits from identified strategies were quantified using site characteristics and a set of corresponding metrics. The net-zero analysis found that 11 sites could use on-site strategies to eliminate all but 2% of emissions generated. The remaining emissions can be offset, for instance through sequestration, executed at the portfolio scale. On-site carbon-free energy was found to reduce 51% of emissions across all sites; efficiency reduced 19% of emissions; sequestration 16%; procured carbon-free energy 15%; fuel switching 1.6%; and fleet electrification 1.3%. Building electrification, however, increased emissions by 4.4%. Different strategies also provide cost, resilience, and/or environmental justice benefits—the degree to which varies with individual site conditions. The findings indicate an advantage to considering the strategies as a comprehensive set, which leads to co-benefits, both in the ability to achieve net-zero goals and in advancing other goals. The results present the case for comprehensive advanced planning at the portfolio level to prioritize investments that will balance the minimization of emissions and life cycle cost with the maximization of resilience and environmental justice benefits. The novel methods for evaluation and integration, valuation of benefits, and consideration at the portfolio scale allow organizations to select investments that simultaneously address multiple key priorities.

Net-Zero Emissions↗

The radial action from probe amplitudes to all orders

We extract the relativistic classical radial action from scattering amplitudes, to all orders in perturbation theory, in the probe limit. Our sources include point charges and monopoles, as well as the Schwarzschild and pure-NUT gravitational backgrounds. A characteristic relativistic effect, that scattering trajectories may wind around these sources any number of times, can be recovered when all-order amplitudes are available. We show that the amplitude for scattering a probe off a pure NUT is given by the solution of a transcendental equation involving continued fractions, and explain how to solve this equation to any desired loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Experimental Investigation on the Static Equilibria and Dynamics of Liquid Bridges

A liquid bridge is a volume of liquid held between two or more solid supports. In the case of small disk supports with a sharp edge, the contact line between the bridge and the support disk will be anchored along the edge of the disk. For these cases the solid presents a geometrical singularity and the contact angle is indeterminate within a given range. This dissertation presents research conducted on liquid bridges with anchored contact lines. The three major topics covered are: determining the role of support geometry on static equilibria, liquid bridge dynamical behavior, and forces exerted by a liquid bridge on a support structure. The work was primarily experimental and conducted in a "Plateau tank" that allowed for the simulation of equivalent low-gravity conditions. The main thrust of the experimental work involved the use of a high resolution optical measurement system for imaging the dynamic zone shape, measurement of the static and dynamic contact angles and non-invasive analysis of excited surface modes. The liquid bridge was manipulated by computer controlled linear actuators which allowed precise control over the physical characteristics of the bridge. Experiments have been carried out to locate a bifurcation point along the maximum volume axisymmetric stability margin. Below the critical slenderness the bifurcation from an axisymmetric to a stable nonaxisymmetric configuration is supercritical. However, above this critical slenderness, the bifurcation is subcritical. A series of experiments analyzed the effect on axisymmetric bridge stability by using support disks of different radii, The shape behavior as transition points were approached, as well as the limiting case of a vanishing support radius was investigated. Experiments were performed to determine the resonant frequencies of axisymmetric bridges subject to lateral vibrations. Anomolous results led to a series of experiments to characterize nonlinearities present in the dynamic bridge shape. Finally, an attempt was made to experimentally measure the force exerted by the bridge on the lower support disk. This was done through use of a force balance apparatus. Particular attention was paid to the behavior of the bridge as the minimum volume stability limit was approached.

Resnick, Andrew Howard↗

One-shot gas detection with transformer paired neural networks in Mako collected longwave infrared hyperspectral imagery

To date, careful data treatment workflows and statistical detectors are used to perform hyperspectral image (HSI) detection of any gas contained in a spectral library, which is often expanded with physics models to incorporate different spectral characteristics. In general, surrounding evidence or known gas-release parameters are used to provide confidence in or confirm detection capability, respectively. This makes quantifying detection performance difficult as it is nearly impossible to develop an absolute ground truth for gas target pixel presence in collected HSI. Consequently, developing and comparing new detection methods, especially machine learning (ML)-based methods, is susceptible to subjectivity in derived detection map quality. Here, in this work, we demonstrate the first use of transformer-based paired neural networks (PNNs) for one-shot gas target detection for multiple gases while providing quantitative classification and detection metrics for their use on labeled data. Terabytes of training data are generated from a database of long-wave infrared HSI obtained from historical Mako sensor campaigns over Los Angeles. By incorporating labels, singular signature representations, and a model development pipeline, we can tune and select PNNs to detect multiple gas targets that are not seen in training on a quantitative basis. We additionally assess our test set detections using interpretability techniques widely employed with ML-based predictors, but less common with detection methods relying on learned latent spaces.

Hyperspectral imaging↗

Analytical method for the ditching analysis of an airborne vehicle

A simple analytical method has been introduced for aerohydrodynamic load analysis of an airborne configuration during water ditching. The method employs an aerodynamic panel code, based on linear potential flow theory, to simulate the flow of air and water around an aircraft configuration. The free surface separating the air and water region is represented by doublet" singularities. Although all the theoretical load distributions are computed for air, provisions are made to correct the pressure coefficients obtained on the configuration-wetted surfaces to account for the water density. As an analytical tool, a vortex aerodynamic code is chosen to carry out the present investigation. After assessing· the validity of the method, an application is presented for the water ditching of the Space Shuttle configuration at a 12-deg attitude.

Pressure coefficient↗

An analytical method for the ditching analysis of an airborne vehicle

A simple analytical method has been introduced for aerohydrodynamic load analysis of an airborne configuration during water ditching. The method employs an aerodynamic panel code, based on linear potential flow theory, to simulate the flow of air and water around an aircraft configuration. The free surface separating the air and water region is represented by doublet sheet singularities. Although all the theoretical load distributions are computed for air, provisions are made to correct the pressure coefficients obtained on the configuration wetted surfaces to account for the water density. As an analytical tool, the Vortex Separation Aerodynamic (VSAERO) code is chosen to carry out the present investigation. After assessing the validity of the method, its first application is to analyze the water ditching of the Space Shuttle configuration at a 12 degree attitude.

Pressure coefficient↗

Possible Evidence for an Event Horizon in Cyg XR-1

The X-ray emitting component in the Cyg XR-1/HDE226868 system is a leading candidate for identification as a stellar-mass sized black hole. The positive identification of a black hole as predicted by general relativity requires the detection of an event horizon surrounding the point singularity. One signature of such an event horizon would be the existence of dying pulse trains emitted by material spiraling into the event horizon from the last stable orbit around the black hole. We observed the Cyg XR-1 system at three different epochs in a 1400 - 3000 A bandpass with 0.1 ms time resolution using the Hubble Space Telescope's High Speed Photometer. Repeated excursions of the detected flux by more than three standard deviations above the mean are present in the UV flux with FWHM 1 - 10 ms. If any of these excursions are pulses of radiation produced in the system (and not just stochastic variability associated with the Poisson distribution of detected photon arrival times), then this short a timescale requires that the pulses originate in the accretion disk around Cyg XR-1. Two series of pulses with characteristics similar to those expected from dying pulse trains were detected in three hours of observation.

Dolan, Joseph F.↗

Alkali-Metal Interlocking of 2D V 4 O 10 Sheets Defines Discretized Interlayer Shear Relationships

Low-dimensional materials manifest structural anisotropy, quantum confinement, and tightly bound excitonic states, which make them attractive building blocks that can be assembled within three-dimensional laterally stitched heterostructures, stacked van der Waals solids, and complex moiré superlattices. Ion intercalation in the galleries between layered materials provides a means of modifying interlayer separation and coupling, but it is also known to drive the shearing of the layers. In this article, we explore the distinct ligand coordination environments afforded by vanadyl oxygens of singular [V 4 O 10 ] sheets and examine how the size, polarizability, and stoichiometry of Group I cations sandwiched between such layers determine the interlocking of the sheets in stacked structures. Based on the topochemical insertion of alkali-metal ions into the layered λ-V 2 O 5 , we identify seven types of guest ion coordination sites discretized into four distinct regimes of interlayer shear in units of half octahedral widths. The coordination preferences of intercalated cations govern how they interlock 2D [V 4 O 10 ] sheets and engender specific shear conformations. We present evidence that static and dynamic disorder in guest ion arrangement modulate the magnetic structure of the intercalated compounds based on electrostatic polarization, localization of charge and spin density, and lattice distortion. The results illustrate the use of topochemical ion insertion to modulate stacking relationships and magnetic transition characteristics.

36 MATERIALS SCIENCE↗

Characteristics of Control Laws Tested on the Semi-Span Super-Sonic Transport (S4T) Wind-Tunnel Model

The Semi-Span Supersonic Transport (S4T) is an aeroelastically scaled wind-tunnel model built to test active controls concepts for large flexible supersonic aircraft in the transonic flight regime. It is one of several models constructed in the 1990's as part of the High Speed Research (HSR) Program. Control laws were developed for the S4T by M4 Engineering, Inc. and by Zona Technologies, Inc. under NASA Research Announcement (NRA) contracts. The model was tested in the NASA-Langley Transonic Dynamics Tunnel (TDT) four times from 2007 to 2010. The first two tests were primarily for plant identification. The third entry was used for testing control laws for Ride Quality Enhancement, Gust Load Alleviation, and Flutter Suppression. Whereas the third entry only tested FS subcritically, the fourth test demonstrated closed-loop operation above the open-loop flutter boundary. The results of the third entry are reported elsewhere. This paper reports on flutter suppression results from the fourth wind-tunnel test. Flutter suppression is seen as a way to provide stability margins while flying at transonic flight conditions without penalizing the primary supersonic cruise design condition. An account is given for how Controller Performance Evaluation (CPE) singular value plots were interpreted with regard to progressing open- or closed-loop to higher dynamic pressures during testing.

Christhilf, David M.↗

Mode properties of low-frequency waves: Kinetic theory versus Hall-MHD

In fluid theory, the ordering of low-frequency modes in a homogeneous plasma is based on the phase velocity, since modes do not intersect each other in dispersion diagrams as a function of wavenumber or other parameters. In linear kinetic theory, modes cross each other. Thus a consistent and useful classification should be based on the physical properties of the modes instead. This paper attempts such a classification by documeting the dispersion and general mode properties of the low-frequency waves (omega much less than (OMEGA(sub ci) OMEGA(sub ce) (exp 1/2)), where OMEGA(sub ci), OMEGA(sub ce) are the cyclotron frequencies of the ions and electrons, respectively) in kinetic theory, and by comparing them to the results of two-fluid theory. Kinetic theory gives a seperate Alfven/ion-cyclotron (A/IC) wave with phase speed Omega/k approximately = v(sub A) cos theta for omega much less than OMEGA(sub ci), where v(sub A) is the Alfven velocity and theta the angle of propagation between wave vector k and background magnetic field B(sub o). For a given wavenumber, the magnetosonic mode is a double-valued solution with a singular point in theta, beta parameter space, where beta is the ratio of thermal pressure to magnetic pressure. It is shown that a branch cut starting at the singular point theta approximately 30 deg, beta approximately 3 and leading to larger beta gives a practical and consitent seperation of this double-valued magnetosonic solution. Selection of this branch cut results in a moderately damped fast/magnetos onic and a heavily damped slow/sound wave. A comprehensive review of the polarization, compressibility and other mode properties is given and shown to be consistent with the selected branch cut. At small wavenumbers, the kinetic mode properties typically start to deviate significantly from their fluid counterparts at beta approximately 0.5. At larger beta, there is no longer a consistent correspondence between the fluid and kinetic modes. Kinetic theory also dictates the use of different mode properties to distinguish between them in observational data. For example, the phase between the density and magnetic field perturbation may become useless at high beta, whereas the direction of the magnetic field perturbations with respect to k and B(sub o) remains a useful characteristic. Two quantities based on this characteristic are suggested and are shown to be useful also to distinguish between the mirror mode and A/IC waves in a plasma with temperature anisotropy.

Krauss-Varban, D.↗

Powering Data Centers with Clean Energy: A Techno-Economic Case Study of Nuclear and Renewable Energy Dependability

Rising data demands from artificial intelligence (AI) and large language models (LLMs) generating images, videos, and text have prompted increased need for larger and more robust data centers in the United States. Major companies interested in these larger data centers face the choice of linking them to existing regional grids, building stand-alone power supplies onsite, or a combination of both. The request, review, and approval process for new transmission lines to grids in the United States, however, has grown in recent years to times spans rivaling those of new construction for nuclear power plants. Building an islanded power supply for each data center is therefore becoming a prominent option. In this case study, several technologies are modeled in techno-economic simulations for long-term system costs subject to fixed electricity demand from a singular data center. A 250 MWe data center is assumed with additional 50 MWe for resiliency. Techno-economic simulations are conducted using the Holistic Energy Resource Optimization Network (HERON) software, which is a part of the Framework for Optimization of Resources and Economics (FORCE) tool suite. Technologies considered include solar, wind, lithium-ion batteries, and several types of nuclear reactors: large-scale reactors, small modular reactors, and microreactors. A low- and high-cost estimate for each technology is assumed to develop a range of expected economic performance. Low-cost estimates included several clean energy production tax credits. Different combinations of renewable energy generators with nuclear reactors are considered, ranging from a fully renewable-powered data center to a fully nuclear-powered data center. Historic time series of wind and solar availability from the Texas grid are used to train a reduced order model; this model then generates unique time series with similar characteristics of the training dataset. Multiple scenarios of weather and subsequent operations are simulated for each renewable-nuclear combination to determine total costs throughout the project lifetime. Fully renewable-powered configurations required large amounts of installed capacity (GW scale) in the simulations to meet the fixed demand of the data center. This is due to some scenarios in the historical dataset which captured low-wind and low-solar days, requiring over-building of these technologies as well as batteries to compensate for the low amounts of electricity generation. Fully nuclear-powered configurations outperformed the fully renewable and mixed renewable-nuclear configurations in terms of cost, with ranges between $1B and $10B in 2023 USDs compared to $40B+ for fully renewable configurations. Of the nuclear technologies, small modular reactors performed better economically than large-scale nuclear models due to lower projected capital costs, and both performed better than the microreactor models. These results demonstrate the applicability of firm, dispatchable electricity resources from baseload generators like nuclear power plants for operating facilities that run at constant power without daily variability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Closing the Loops on Solar Photovoltaics Modules: An Agent-Based Modeling Approach for the Study of Circular Economy Strategies

Solar photovoltaics (PV) installed capacity have grown exponentially since the early 2000s (average annual growth rate of 50%). With 4,700 GW projected installed capacity by 2050, the volume of PV panels waste is also expected to become substantial. Though renewables are a sine qua non to the establishment of a truly circular economy (CE), the issue arising from their end-of-life management needs to be resolved. The management of PV end-of-life also represent a singular opportunity to create value, from recovered valuable, rare or critical materials (e.g. silver, tellurium, indium). Moreover, the PV case illustrates some of the current barriers to CE; for instance, the need for a common definition of waste (PV are defined as e-waste in the European Union but as general waste in the United States) and potential loss of innovations' advantages (PV average efficiency has continuously grown the past 10 years). In this study, an agent-based modeling (ABM) approach is proposed to simulate PV end-of-life management in the United States. The model explores how the decisions of the various actors involved in handling PV waste affect the quantities of PV that are reused, recycled, or land-filled. In the model, manufacturers, residential, and nonresidential PV owners, installers, and recyclers are represented by agents while governmental policies and regulations are treated as exogenous variables. PV-market data are used to define agents' characteristics (e.g., installed PV capacities or producing costs). Agents' decision rules draw on the literature related to industrial symbiosis (IS) and peoples' waste behaviors as they appear as the most widely used model to implement CE principles at a meso level. Specifically, the concepts of mutual trust between IS actors and knowledge about the IS philosophy are yielded to define the agents' decision process. The primary outputs of the ABM are the costs and volume of waste associated with each end-of-life pathway. Preliminaries results indicate recycling of residential PV is affected by the cost, perceived difficulty of recycling behaviors, and social norms. Moreover, the type of network connecting agents and the initial recycling rate strongly influence results. The model also highlights the crucial role of recycling behavior adoption: in case of early failure of PV, recycled volumes increase by about 40%. This represents an additional 5.7 billion USD from potential recovered silver. Finally, although results for the PV case are presented here, the developed ABM aims at being a general tool for the study of CE strategies and the CE transition. Furthermore, steps of this research include extending the model to encompass circular economy strategies (e.g., design for recycling or lifetime extension) and validating its general architecture from various case studies.

14 SOLAR ENERGY↗

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗

The Cretaceous-Tertiary Impact Crater and the Cosmic Projectile that Produced it

Evidence gathered to date from topographic data, geophysical data, well logs, and drill-core samples indicates that the buried Chicxulub basin, the source crater for the approximately 65 Ma Cretaceous-Tertiary (K/T) boundary deposits, is approximately 300 km in diameter. A prominent topographic ridge and a ring of gravity anomalies mark the position of the basin rim at approximately 150 km from the center. Wells in this region recovered thick sequences of impact-generated breccias at 200-300 m below present sea level. Inside the rim, which has been severely modified by erosion following impact, the subsurface basin continues to deepen until near the center it is approximately 1 km deep. The best planetary analog for this crater appears to be the 270 km-diameter Mead basin on Venus. Seismic reflection data indicate that the central zone of downward displacement and excavation (the transient crater is approximately 130 km in diameter, consistent with previous studies of gravity anomaly data). Our analysis of projectile characteristics utilizes this information, coupled with conventional scaling relationships, and geochemical constraints on the mass of extraterrestrial material deposited within the K/T boundary layer. Results indicate that the Chicxulub crater would most likely be formed by a long-period comet composed primarily of nonsilicate materials (ice, hydrocarbons, etc.) and subordinate amounts (less than or equal to 50 percent) primitive chondritic material. This collision would have released the energy equivalent to between 4 x 10(exp 8) and 4 x 10(exp 9) megatons of TNT. Studies of terrestrial impact rates suggest that such an event would have a mean production rate of approximately 1.25 x 10(exp -9) y(exp -1). This rate is considerably lower than that of the major mass extinctions over the last 250 million years (approximately 5 x 10(exp -7) y(exp -1). Consequently, while there is substantial circumstantial evidence establishing the cause-effect link between the Chicxulub basin forming event and the K/T biological extinctions, the results of our analysis do not support models of impact as a common or singular causative agent of mass extinctions on Earth.

Sharpton, Virgil L.↗