Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global refinement”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Newmark local time stepping on high-performance computing architectures

In multi-scale complex media, finite element meshes often require areas of local refinement, creating small elements that can dramatically reduce the global time-step for wave-propagation problems due to the CFL condition. Local time stepping (LTS) algorithms allow an explicit time-stepping scheme to adapt the time-step to the element size, allowing near-optimal time-steps everywhere in the mesh. We develop an efficient multilevel LTS-Newmark scheme and implement it in a widely used continuous finite element seismic wave-propagation package. In particular, we extend the standard LTS formulation with adaptations to continuous finite element methods that can be implemented very efficiently with very strong element-size contrasts (more than 100x). Capable of running on large CPU and GPU clusters, we present both synthetic validation examples and large scale, realistic application examples to demonstrate the performance and applicability of the method and implementation on thousands of CPU cores and hundreds of GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Algorithms in Diffraction Profile Analysis

This chapter introduces several approaches to crystallographic structure refinement from diffraction profile data including the Rietveld method, global optimization methods, and Bayesian modeling. First, the Rietveld method and related nonlinear least squares approaches to analyzing diffraction profiles are briefly discussed. We then review the most common algorithms within global optimization as alternative methods of structure solution: grid search, simulated annealing, and genetic algorithms. As the Rietveld method is the standard tool for analyzing diffraction profile data, we catalogue its challenges and limitations. A Bayesian approach to diffraction profile analysis is presented as an alternative, along with several sampling algorithms including Gibbs sampling, variants of Metropolis sampling, Hamiltonian Monte Carlo, and approximate Bayesian computing. These sampling algorithms are applied to the fitting of a neutron diffraction profile from a National Institute of Standards and Technology silicon standard reference material, and their results are compared.

Singh, Susheela↗
DOE OSTI2019

CMDV-RRM: Representation of Clouds and Convection across scales in E3SM

In order to address mission-relevant questions, DOE’s Energy Exascale Earth System Model (E3SM) is pushing the limits of horizontal resolution among global Earth System Models (ESMs). The atmospheric resolution in the recently released E3SMv1 is 110 km in the standard configuration, and 25 km in the high-resolution configuration. Furthermore, E3SM is fairly unique among global ESMs with a capability for multiresolution modeling using refined grids in most of its components (atmosphere, land, ocean, sea-ice). This Regionally Refined Model (RRM) capability is critical for future E3SM simulation campaigns that have a strong regional focus to meet DOE's needs in support of energy-sector decisions. RRM grids are also particularly well suited for model evaluation and process-level studies taking advantage of the vast amount of ASR site data collected over the past several decades. The overarching goal of this project was to understand and improve the behavior of clouds and convection in E3SM across these scales using RRM. Key deliverables that resulted from the project are summarized below, with more details provided in the subsequent sections of the report.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Improved Time-Stepping Methods in Global to Regional Ocean Modeling

Improvements in time-stepping schemes can provide substantial improvements in the throughput of global ocean-climate models, but is challenging due to large variations in time scales in the dynamics and spatial scales in the global, variable resolution mesh. Two funded projects aim to speed up the Model for Prediction Across Scales-Ocean (MPAS-Ocean). The first will add a variable-timestepping formulation to MPAS-Ocean, which is a critical step towards running efficiently on global variable-resolution meshes with regional and coastal refinement. The second optimizes the barotropic (depth-integrated) mode, which is currently a performance bottleneck at high core counts. The proposed compute time will be used for extensive model validation and parameter studies for both methods.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1980

Grid and metric generation on the assembly of locally bi-quadratic coordinate transformations

The generation of metric coefficients of the coordinate transformation from a generally curved-sided domain boundary to the unit square (cube) is required for efficient solution algorithms in computational fluid mechanics. An algebraic procedure is presented for establishment of these data on the union of arbitrarily selected sub-domains of the global solution domain. A uniformly smooth progression of grid refinement is readily generated, including multiple specification of refined grids for a given macro-element domain discretization. The procedure is illustrated as generally applicable to non-simply connected domains in two and three dimensions.

Baker, A. J.↗
NASA NTRS1991

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗
NASA NTRS1991

Analysis and representation of complex structures in separated flows

We discuss our recent work on extraction and visualization of topological information in separated fluid flow data sets. As with scene analysis, an abstract representation of a large data set can greatly facilitate the understanding of complex, high-level structures. When studying flow topology, such a representation can be produced by locating and characterizing critical points in the velocity field and generating the associated stream surfaces. In 3D flows, the surface topology serves as the starting point. The 2D tangential velocity field near the surface of the body is examined for critical points. The tangential velocity field is integrated out along the principal directions of certain classes of critical points to produce curves depicting the topology of the flow near the body. The points and curves are linked to form a skeleton representing the 2D vector field topology. This skeleton provides a basis for analyzing the 3D structures associated with the flow separation. The points along the separation curves in the skeleton are used to start tangent curve integrations. Integration origins are successively refined to produce stream surfaces. The map of the global topology is completed by generating those stream surfaces associated with 3D critical points.

Helman, James↗
NASA NTRS1999

Studies for the Loss of Atomic and Molecular Species from IO

The general objective of this project has been to advance our theoretical understanding of lo's atmosphere and how various atomic and molecular species are lost from this atmosphere and are distributed in the circumplanetary environment of Jupiter. This grant has provided support for the activities of Dr. Michael Combi at the University of Michigan to serve as a small part in collaboration with a larger project awarded to Atmospheric & Environmental Research, Inc., with primary principal investigator Dr. William H. Smyth. Dr. Combi is the Principal Investigator and Project Manager for the Michigan grant NAG5-6187. This Michigan grant has provided for a continuation of a collaboration between Drs. Smyth and Combi in related efforts beginning in 1981, and with the object to develop and apply sophisticated theoretical models to interpret and to relate a number of new and exciting observations for the atmospheric gases of the satellite. The ability to interpret and then to relate through the theoretical fabric a number of these otherwise independent observations are a central strength of this program. This comprehensive approach provides a collective power, extracting more from the sum of the parts and seeing beyond various limitations that are inherent in any one observation. Although the approach is designed to unify, the program is divided into well-defined studies for the likely dominant atmospheric gases involving species of the SO2 family (SO2, SO, O2, S and O) and for the trace atmospheric gas atomic sodium and a likely escaping molecular ion NaX(+) (where Na(X) is the atmospheric molecule and X represents one or more atoms).Attachments: IO's sodium corona and spatially cloud: a consistent flux speed distribution. and Io's plasma environment during the Galileo flyby: global three-dimensional MHD modeling with adaptive mesh refinement.

Combi, Michael R.↗
NASA NTRS1993

Venus small volcano classification and description

The high resolution and global coverage of the Magellan radar image data set allows detailed study of the smallest volcanoes on the planet. A modified classification scheme for volcanoes less than 20 km in diameter is shown and described. It is based on observations of all members of the 556 significant clusters or fields of small volcanoes located and described by this author during data collection for the Magellan Volcanic and Magmatic Feature Catalog. This global study of approximately 10 exp 4 volcanoes provides new information for refining small volcano classification based on individual characteristics. Total number of these volcanoes was estimated to be 10 exp 5 to 10 exp 6 planetwide based on pre-Magellan analysis of Venera 15/16, and during preparation of the global catalog, small volcanoes were identified individually or in clusters in every C1-MIDR mosaic of the Magellan data set. Basal diameter (based on 1000 measured edifices) generally ranges from 2 to 12 km with a mode of 34 km, and follows an exponential distribution similar to the size frequency distribution of seamounts as measured from GLORIA sonar images. This is a typical distribution for most size-limited natural phenomena unlike impact craters which follow a power law distribution and continue to infinitely increase in number with decreasing size. Using an exponential distribution calculated from measured small volcanoes selected globally at random, we can calculate total number possible given a minimum size. The paucity of edifice diameters less than 2 km may be due to inability to identify very small volcanic edifices in this data set; however, summit pits are recognizable at smaller diameters, and 2 km may represent a significant minimum diameter related to style of volcanic eruption. Guest, et al, discussed four general types of small volcanic edifices on Venus: (1) small lava shields; (2) small volcanic cones; (3) small volcanic domes; and (4) scalloped margin domes ('ticks'). Steep-sided domes or 'pancake domes', larger than 20 km in diameter, were included with the small volcanic domes. For the purposes of this study, only volcanic edifices less than 20 km in diameter are discussed. This forms a convenient cutoff since most of the steep-sided domes ('pancake domes') and scalloped margin domes ('ticks') are 20 to 100 km in diameter, are much less numerous globally than are the smaller diameter volcanic edifices (2 to 3 orders of magnitude lower in total global number), and do not commonly occur in large clusters or fields of large numbers of edifices.

Aubele, J. C.↗
NASA NTRS1990

Evidence from carbon isotope measurements for diverse origins of sedimentary hydrocarbons

The organic matter found in sedimentary rocks must derive from many sources; not only from ancient primary producers but also from consumers and secondary producers. In all of these organisms, isotope effects can affect the abundance and distribution of 13C in metabolites. Here, by using an improved form of a previously described technique in which the effluent of a gas chromatograph is continuously analysed isotopically, we report evidence of the diverse origins of sedimentary organic matter. The record of 13C abundances in sedimentary carbonate and total organic carbon can be interpreted in terms of variations in the global carbon cycle. Our results demonstrate, however, that isotope variations within sedimentary organic mixtures substantially exceed those observed between samples of total organic carbon. Resolution of isotope variations at the molecular level offers a new and convenient means of refining views both of localized palaeoenvironments and of control mechanisms within the global carbon cycle.

Non-NASA Center↗
NASA NTRS2007

Structural Analysis of the Redesigned Ice/Frost Ramp Bracket

This paper describes the interim structural analysis of a redesigned Ice/Frost Ramp bracket for the Space Shuttle External Tank (ET). The proposed redesigned bracket consists of mounts for attachment to the ET wall, supports for the electronic/instrument cables and propellant repressurization lines that run along the ET, an upper plate, a lower plate, and complex bolted connections. The eight nominal bolted connections are considered critical in the summarized structural analysis. Each bolted connection contains a bolt, a nut, four washers, and a non-metallic spacer and block that are designed for thermal insulation. A three-dimensional (3D) finite element model of the bracket is developed using solid 10-node tetrahedral elements. The loading provided by the ET Project is used in the analysis. Because of the complexities associated with accurately modeling the bolted connections in the bracket, the analysis is performed using a global/local analysis procedure. The finite element analysis of the bracket identifies one of the eight bolted connections as having high stress concentrations. A local area of the bracket surrounding this bolted connection is extracted from the global model and used as a local model. Within the local model, the various components of the bolted connection are refined, and contact is introduced along the appropriate interfaces determined by the analysts. The deformations from the global model are applied as boundary conditions to the local model. The results from the global/local analysis show that while the stresses in the bolts are well within yield, the spacers fail due to compression. The primary objective of the interim structural analysis is to show concept viability for static thermal testing. The proposed design concept would undergo continued design optimization to address the identified analytical assumptions and concept shortcomings, assuming successful thermal testing.

Phillips, D. R.↗
NASA NTRS2008

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗
NASA NTRS2017

Improving Estimation of Ground Casualty Risk From Reentering Space Objects

A recent improvement to the long-term estimation of ground casualties from reentering space debris is the further refinement and update to the human population distribution. Previous human population distributions were based on global totals with simple scaling factors for future years, or a coarse grid of population counts in a subset of the world's countries, each cell having its own projected growth rate. The newest population model includes a 5-fold refinement in both latitude and longitude resolution. All areas along a single latitude are combined to form a global population distribution as a function of latitude, creating a more accurate population estimation based on non-uniform growth at the country and area levels. Previous risk probability calculations used simplifying assumptions that did not account for the ellipsoidal nature of the Earth. The new method uses first, a simple analytical method to estimate the amount of time spent above each latitude band for a debris object with a given orbit inclination and second, a more complex numerical method that incorporates the effects of a non-spherical Earth. These new results are compared with the prior models to assess the magnitude of the effects on reentry casualty risk.

Ostrom, Chris L.↗
NASA NTRS2017

Improving Estimation of Ground Casualty Risk from Reentering Space Objects

A recent improvement to the long-term estimation of ground casualties from reentering space debris is the further refinement and update to the human population distribution. Previous human population distributions were based on global totals with simple scaling factors for future years, or a coarse grid of population counts in a subset of the world's countries, each cell having its own projected growth rate. The newest population model includes a 5-fold refinement in both latitude and longitude resolution. All areas along a single latitude are combined to form a global population distribution as a function of latitude, creating a more accurate population estimation based on non-uniform growth at the country and area levels. Previous risk probability calculations used simplifying assumptions that did not account for the ellipsoidal nature of the earth. The new method uses first, a simple analytical method to estimate the amount of time spent above each latitude band for a debris object with a given orbit inclination, and second, a more complex numerical method that incorporates the effects of a non-spherical Earth. These new results are compared with the prior models to assess the magnitude of the effects on reentry casualty risk.

Ostrom, C.↗
NASA NTRS2012

Twenty Years of Progress on Global Ocean Tides: The Impact of Satellite Altimetry

At the dawn of the era of high-precision altimetry, before the launch of TOPEX/Poseidon, ocean tides were properly viewed as a source of noise--tidal variations in ocean height would represent a very substantial fraction of what the altimeter measures, and would have to be accurately predicted and subtracted if altimetry were to achieve its potential for ocean and climate studies. But to the extent that the altimetry could be severely contaminated by tides, it also represented an unprecedented global-scale tidal data set. These new data, together with research stimulated by the need for accurate tidal corrections, led to a renaissance in tidal studies in the oceanographic community. In this paper we review contributions of altimetry to tidal science over the past 20 years, emphasizing recent progress. Mapping of tides has now been extended from the early focus on major constituents in the open ocean to include minor constituents, (e.g., long-period tides; non-linear tides in shelf waters, and in the open ocean), and into shallow and coastal waters. Global and spatially local estimates of tidal energy balance have been refined, and the role of internal tide conversion in dissipating barotropic tidal energy is now well established through modeling, altimetry, and in situ observations. However, energy budgets for internal tides, and the role of tidal dissipation in vertical ocean mixing remain controversial topics. Altimetry may contribute to resolving some of these important questions through improved mapping of low-mode internal tides. This area has advanced significantly in recent years, with several global maps now available, and progress on constraining temporally incoherent components. For the future, new applications of altimetry (e.g., in the coastal ocean, where barotropic tidal models remain inadequate), and new mission concepts (studies of the submesoscale with SWOT, which will require correction for internal tides) may bring us full circle, again pushing further development of tidal models as corrections.

Egbert, Gary↗
NASA NTRS1991

Monitoring, analyzing, and modeling global climate

Diabatic heating rate estimates as residuals of the dry thermodynamic equation were generated for May 1, 1985 to December 1989 in pentad resolution. Published results show moderate correlations (approx. .6) between heating rate and outgoing long wave radiation for periods under 90 days in the tropics and many extratropical locations. Nine years of simulation with the Community Climate Model 1 (CCM1) using R15 and observed sea surface temperatures shows that the model retains significantly more heat at the surface and in the free atmosphere than does the actual earth system. A post-processor for the CCM1, with capabilities to produce simulated microwave sounding unit (MSU) brightness temperatures was written. Techniques were refined considerably and validation studies were carried out to verify the globally distributed free atmosphere temperature anomalies derived from MSU data. The precision is such that detailed, long-term climate monitoring is well within the capability of these data.

Christy, John R.↗
NASA NTRS1988

Geological remote sensing signatures of terrestrial impact craters

Geological remote sensing techniques can be used to investigate structural, depositional, and shock metamorphic effects associated with hypervelocity impact structures, some of which may be linked to global Earth system catastrophies. Although detailed laboratory and field investigations are necessary to establish conclusive evidence of an impact origin for suspected crater landforms, the synoptic perspective provided by various remote sensing systems can often serve as a pathfinder to key deposits which can then be targetted for intensive field study. In addition, remote sensing imagery can be used as a tool in the search for impact and other catastrophic explosion landforms on the basis of localized disruption and anomaly patterns. In order to reconstruct original dimensions of large, complex impact features in isolated, inaccessible regions, remote sensing imagery can be used to make preliminary estimates in the absence of field geophysical surveys. The experienced gained from two decades of planetary remote sensing of impact craters on the terrestrial planets, as well as the techniques developed for recognizing stages of degradation and initial crater morphology, can now be applied to the problem of discovering and studying eroded impact landforms on Earth. Preliminary results of remote sensing analyses of a set of terrestrial impact features in various states of degradation, geologic settings, and for a broad range of diameters and hence energies of formation are summarized. The intention is to develop a database of remote sensing signatures for catastrophic impact landforms which can then be used in EOS-era global surveys as the basis for locating the possibly hundreds of missing impact structures. In addition, refinement of initial dimensions of extremely recent structures such as Zhamanshin and Bosumtwi is an important objective in order to permit re-evaluation of global Earth system responses associated with these types of events.

Garvin, J. B.↗