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

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.↗

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.↗

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.↗

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↗

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.↗

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.↗

Convection-permitting hindcasting of diurnal variation of Mei-yu rainfall over East China with a global variable-resolution model

Two sets of 48-hr hindcast experiments for Mei-yu season using non-hydrostatic global variable-resolution model (MPAS-Atmosphere) with regional refinement at 4km resolution (V4km) and uniform 60km resolution (U60km) are first-time conducted to attempt at simulating the Mei-yu rainfall in 2015 over East China. The analysis focuses on the seasonal average of diurnal variation of precipitation amount (PA), intensity (PI), and frequency (PF). Both simulations reasonably reproduce the spatial distribution of PA. V4km is more skillful in simulating the spatial distributions and magnitudes of PF and PI. The diurnal cycle of Mei-yu rainfall shows a major early morning peak and a minor afternoon peak, contributed by precipitation during two sub-periods with distinct synoptic circulations. With strong Mei-yu in the first sub-period, the diurnal variation of PA is controlled by nocturnal southwesterly jet. V4km overestimates the morning peak mainly due to its bias in simulating boundary layer inertial oscillation. Although differences in PA between the simulations are small, U60km overestimates PF and underestimates PI. With weak Mei-yu in the second sub-period, the diurnal variation of PA is controlled by both synoptic circulation and local convection. At both resolutions, deviations in the large-scale circulation modulated by a few typhoons lead to positive biases in the morning peak of PA. After removing the typhoon impacts, V4km captures the observed diurnal cycle of PA well, while U60km significantly underestimates PA and PI particularly in the afternoon. Future studies focusing on advancing modeling of southwesterly jet and typhoons may further improve convection permitting simulation of Mei-yu rainfall.

Xu, Mingyue↗

Scaling Uintah on the Aurora Exascale System up to 122,880 Intel Ponte Vecchio Xe Stacks

The challenge of being able to scale application codes based on the Asynchronous Many-Task (AMT) Uintah framework on the Department of Energy (DOE) Aurora exascale system is addressed in this work by considering a challenging Reverse Monte Carlo Ray Tracing radiation benchmark calculation. This benchmark involves potentially global all-to-all communication and uses adaptive mesh refinement and ray tracing to achieve scalability. This benchmark has been used as part of previous scalability studies on a number of pre-exascale systems and on the DOE Frontier exascale system. This paper describes steps taken to enable this benchmark to run successfully on up to 10,240 nodes and 122,880 Intel® Ponte Vecchio Xe stacks on the DOE Aurora exascale system. This scalability was achieved through a limited number of experiments on Aurora, given machine loads and its uniqueness. These experiments constitute valuable lessons learned to achieve scalability at this level. The resulting scalability runs, while few in number, demonstrate relatively good strong-scaling characteristics. A detailed analysis of these results provides important indications about the path to scalability on Aurora for future work. Overall, these results continue the remarkable ability of this AMT approach to produce scalable solutions for challenging problems at extreme scale on heterogeneous architectures.

Garcia, Marta [Argonne National Laboratory (ANL)]↗

Finite element substructuring methods for composite mechanics

Finite element substructuring strategies are presented to obtain numerical solutions for three typical problems of interest to the composites community: (1) impact and toughness characterization of composites using Charpy's impact test specimen; (2) free-edge stress analysis of composite laminates; and (3) fracture toughness predictions of composites for individual and combined fracture of modes I, II, and III. The key issue common to these problems is the presence of singular or near singular stress fields. The regions prone to see steep stress gradients are substructured with progressively refined meshes to study the local response simultaneously with the global response. The results from the select examples indicate that finite element substructuring methods are computationally effective for composite singularity mechanics.

Murthy, Pappu L. N.↗

Elastic-plastic crack analysis using a global-local approach on a parallel computer

A global-local finite-element analysis procedure is employed to solve elastoplastic crack problems. This procedure involves a global analysis using a coarse mesh and the subsequent locally refined analysis of a number of subregions. Parallel computation is involved in the local analyses. Numerical examples are performed on a MIMD parallel computer to demonstrate the time-saving capability of the proposed procedure.

Sun, C. T.↗

Finite element substructuring methods for composite mechanics

Finite element substructuring strategies are presented to obtain numerical solutions for three typical problems of interest to the composites community: (1) impact and toughness characterization of composites using Charpy's impact test specimen; (2) free-edge stress analysis of composite laminates; and (3) fracture toughness predictions of composites for individual and combined fracture of modes I, II, and III. The key issue common to these problems is the presence of singular or near singular stress fields. The regions prone to see steep stress gradients are substructured with progressively refined meshes to study the local response simultaneously with the global response. The results from the select examples indicate that finite element substructuring methods are computationally effective for composite singularity mechanics.

Murthy, P. L. N.↗

IMERG V06: Changes to the Morphing Algorithm

As the US Science Team’s globally gridded precipitation product from the NASA/JAXA Global Precipitation Measurement (GPM) mission, the Integrated Multi-satellitE Retrievals for GPM (IMERG) estimates the surface precipitation rates at 0.1° every half-hour using spaceborne sensors for various scientific and societal applications. One key component of IMERG is the morphing algorithm, which uses motion vectors to perform quasi-Lagrangian interpolation to fill in gaps in the passive microwave precipitation field using motion vectors. Up to IMERG V05, the motion vectors were derived from the large-scale motions of infrared observations of cloud tops. This study details the changes introduced in IMERG V06 to derive motion vectors from large-scale motions of selected atmospheric variables in numerical models, which allow IMERG estimates to be extended from the 60°N/S latitude band to the entire globe. Evaluation against both instantaneous passive microwave retrievals and ground measurements demonstrates the general improvement in the precipitation field of the new approach. Most of the model variables tested exhibited similar performance, but total precipitable water vapor was chosen as the source of the motion vectors for IMERG V06 due to its competitive performance and global completeness. Continuing assessments will provide further insights into possible refinements of this revised morphing scheme in future versions of IMERG.

Jackson Tan↗

Enabling Real-time Scattering Data Analysis with Scalable Optimization [Slides]

Diffraction experiments produce datasets with rich multidimensional physics information such as microstructure, equations of state, crystal structure, elastoplastic properties, and other key inputs to LANL mission-essential multiphysics models. This information is typically extracted through a process called Rietveld refinement, which involves selecting appropriate models of the instrument, crystal structure, and microstructure, identifying suitable starting values, and then fitting often hundreds of model parameters using a sequence of empirical parameter turnon/off sequences within a non-global gradient-based optimization. Extensive user expertise is required to properly setup a refinement, identify appropriate models, and select initial parameter values close to truth, such that the refinement will yield parameter values that are optimally predictive. This is a very tedious manual process performed far after the beamline campaign has ended. As facilities have become capable of generating larger volumes of data, the limitation in throughput due to Rietveld refinement has led to a dramatic increase in unanalyzed data as opposed to an intended increase in new science. In our FY22 TED, we demonstrated an integrated toolset providing near real-time automated Rietveld analysis. If this toolset can be optimized to provide automated Rietveld analysis in real-time, this could alleviate the bottleneck in unanalyzed diffraction data, aid in decision-making during experiments, and increase efficiency of the facility.

74 ATOMIC AND MOLECULAR PHYSICS↗

Leveraging regional mesh refinement to simulate future climate projections for California using the Simplified Convection-Permitting E3SM Atmosphere Model Version 0

Abstract. The spatial heterogeneity related to complex topography in California demands high-resolution (< 5 km) modeling, but global convection-permitting climate models are computationally too expensive to run multi-decadal simulations. We developed a 3.25 km California climate modeling framework by leveraging regional mesh refinement (CARRM) using the U.S. Department of Energy (DOE)'s global Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) version 0. Four 5-year time periods (2015–2020, 2029–2034, 2044–2049, and 2094–2099) were simulated by nudging CARRM outside California to 1° coupled simulation of E3SMv1 under the Shared Socioeconomic Pathways (SSP)5-8.5 future scenario. The 3.25 km grid spacing adds considerable value to the prediction of the California climate changes, including more realistic high temperatures in the Central Valley and much improved spatial distributions of precipitation and snowpack in the Sierra Nevada and coastal stratocumulus. Under the SSP5-8.5 scenario, CARRM simulation predicts widespread warming of 6–10 °C over most of California, a 38 % increase in statewide average 30 d winter–spring precipitation, a near-complete loss of the alpine snowpack, and a sharp reduction in shortwave cloud radiative forcing associated with marine stratocumulus by the end of the 21st century. We note a climatological wet precipitation bias for the CARRM and discuss possible reasons. We conclude that SCREAM RRM is a technically feasible and scientifically valid tool for climate simulations in regions of interest, providing an excellent bridge to global convection-permitting simulations.

58 GEOSCIENCES↗

Site selection for Mars exobiology

The selection of sites on Mars that have a high priority for exobiological research is fundamental for planning future exploration. The most immediate need is to identify targets for high resolution orbital imaging during the Mars Observer and Mars '94/'96 missions that can be used to refined site priorities for surface exploration. We present an objective approach to site selection whereby individual sites are selected and scored, based on the presence of key geological features which indicate high priority environments. Prime sites are those that show evidence for the prolonged activity of liquid water and which have sedimentary deposits that are likely to have accumulated in environments favorable for life. High priority areas include fluvio-lacustrine (stream-fed lake systems), springs, and periglacial environments. Sites where mineralization may have occurred in the presence of organisms (e.g. springs) are given high priority in the search for a fossil record on Mars. A systematic review of Viking data for 83 sites in the Mars Landing Site Catalog (MLSC) resulted in the selection of 13 as being of exobiological interest. The descriptions of these sites were expanded to address exobiological concerns. An additional five sites were identified for inclusion in the second edition of the MLSC. We plan to broaden our site selection activities to include a systematic global reconnaissance of Mars using Viking data, and will continue to refine site priorities for exobiological research based on data from future missions in order to define strategies for surface exploration.

Farmer, J.↗

Global stochastic optimization of stellarator coil configurations

In the construction of a stellarator, the manufacturing and assembling of the coil system is a dominant cost. These coils need to satisfy strict engineering tolerances, and if those are not met the project could be cancelled as in the case of the National Compact Stellarator Experiment (NCSX) project. Therefore, our goal is to find coil configurations that increase construction tolerances without compromising the performance of the magnetic field. In this paper, we develop a gradient-based stochastic optimization model which seeks robust stellarator coil configurations in high dimensions. In particular, we design a two-step method: first, we perform an approximate global search by a sample efficient trust-region Bayesian optimization; second, we refine the minima found in step one with a stochastic local optimizer. To this end, we introduce two stochastic local optimizers: BFGS applied to the sample average approximation; and Adam, equipped with a control variate for variance reduction. Numerical simulations performed on a W7-X-like coil configuration demonstrate that our global optimization approach finds a variety of promising local solutions at less than 0.1% of the cost of previous work, which considered solely local stochastic optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spaceborne Remote Sensing of Aerosol Type: Global Distribution, Model Evaluation and Translation into Chemical Speciation

It is essential to evaluate and refine aerosol classification methods applied to passive satellite remote sensing. We have developed an aerosol classification algorithm (called Specified Clustering and Mahalanobis Classification, SCMC) that assigns an aerosol type to multi-parameter retrievals by spaceborne, airborne or ground-based passive remote sensing instruments [1]. The aerosol types identified by our scheme are pure dust, polluted dust, urban-industrialdeveloped economy, urban-industrialdeveloping economy, dark biomass smoke, light biomass smoke and pure marine. We apply the SCMC method to inversions from the ground-based AErosol RObotic NETwork (AERONET [2]) and retrievals from the space-borne Polarization and Directionality of Earths Reflectances instrument (POLDER, [3]). The POLDER retrievals that we use differ from the standard POLDER retrievals [4] as they make full use of multi-angle, multispectral polarimetric data [5]. We analyze agreement in the aerosol types inferred from both AERONET and POLDER and evaluate GEOS-Chem [6] simulations over the globe. Finally, we use in-situ observations from the SEAC4RS airborne field experiment to bridge the gap between remote sensing-inferred qualitative SCMC aerosol types and their corresponding quantitative chemical speciation. We apply the SCMC method to airborne in-situ observations from the NASA Langley Aerosol Research Group Experiment (LARGE, [7]) and the Differential Aerosol Sizing and Hygroscopicity Spectrometer Probe (DASH-SP, [8]) instruments; we then relate each coarsely defined SCMC type to a sum of percentage of individual aerosol species, using in-situ observations from the Particle Analysis by Laser Mass Spectrometry (PALMS, [9]), the Soluble Acidic Gases and Aerosol (SAGA, [10]), and the High - Resolution Time - of - Flight Aerosol Mass Spectrometer (HR ToF AMS, [11]).

airborne↗