SEARCH · Engineering Papers
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.
Well-balanced reconstruction techniques have been developed for stellar hydrodynamics to address the challenges of maintaining hydrostatic equilibrium during evolution. I show how to adapt a simple well-balanced method to the piecewise parabolic method for hydrodynamics. A python implementation of the method is provided.
We present joint South Pole Telescope and XMM-Newton observations of eight massive galaxy clusters (0.8–2 × 10$^{15}$ M$_{⊙}$) spanning a redshift range of 0.16–0.35. Employing a novel Sunyaev–Zel’dovich + X-ray fitting technique, we effectively constrain the thermodynamic properties of these clusters out to the virial radius. The resulting best-fit electron density, deprojected temperature, and deprojected pressure profiles are in good agreement with previous observations of massive clusters. For the majority of the cluster sample (five out of eight clusters), the entropy profiles exhibit a self-similar behavior near the virial radius. We further derive hydrostatic mass, gas mass, and gas fraction profiles for all clusters up to the virial radius. Comparing the enclosed gas fraction profiles with the universal gas fraction profile, we obtain nonthermal pressure fraction profiles for our cluster sample at >0.5R$_{500}$, demonstrating a steeper increase between R$_{500}$ and R$_{200}$ that is consistent with the hydrodynamical simulations. Our analysis yields nonthermal pressure fraction ranges of 8%–28% (median: 15% ± 11%) at R$_{500}$ and 21%–35% (median: 27% ± 12%) at R$_{200}$. Notably, weak-lensing mass measurements are available for only four clusters in our sample, and our recovered total cluster masses, after accounting for nonthermal pressure, are consistent with these measurements.
Abstract. Stellar occultation observations from space can probe the stratosphere and mesosphere at a fine vertical scale around the globe. Unlike other measurement techniques like radiosondes and aircraft, stellar occultation has the potential to observe the atmosphere above 30 km, and unlike radio occultation, stellar occultation probes fine-scale phenomena with potential to observe atmospheric turbulence. We imaged the refractive bending angle of a star centroid for a series of occultations by the atmosphere. Atmospheric refractivity, density, and then temperature are retrieved from the bending observations with the Abel transformation and Edlén's law, the hydrostatic equation, and the ideal gas law. The retrieval technique is applied to data collected by two nanosatellites operated by Terran Orbital. Measurements were primarily taken by the GEOStare SV2 mission, with a dedicated imaging telescope, supplemented with images captured by spacecraft bus sensors, namely the star trackers on other Terran Orbital missions. The bending angle noise floor is 10 and 30 arcsec for the star tracker and GEOStare SV2 data, respectively. The most significant sources of uncertainty are due to centroiding errors due to the fairly low-resolution stellar images and telescope pointing knowledge derived from noisy satellite attitude sensors. The former mainly affects the star tracker data, while the latter limits the GEOStare SV2 accuracy, with both providing low vertical resolution. This translates to a temperature profile retrieval up to roughly 20 km for both star tracker and GEOStare SV2 datasets. In preparation of an upcoming 2023 mission designed to correct these deficiencies, SOHIP, we simulated bending angle measurements with varying magnitudes of error. The expected maximum altitude of retrieved temperature is 41 km on average for these simulated measurements with a noise floor of 0.39 arcsec. Our work highlights the capabilities of stellar occultation observations from nanosatellites for atmospheric sounding. Future work will investigate high-frequency observations of atmospheric gravity waves and turbulence, mitigating the major uncertainties observed in these datasets.
Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.
Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.
We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.
This data collection was contributed to the Visualisation Hackathon 2022 (#VisMetHack2022), in conjunction with the Using ECMWF's Forecasts (UEF2022) workshop. The European Center for Medium-Range Weather Forecasts (ECMWF) and the Oak Ridge National Laboratory (ORNL) are pleased to announce access to the data collection from global 1-km nature run (NR) simulations using the Integrated Forecast System (IFS) with explicit convection. We invite you to join us in exploring this precursor to a digital twin of the earth! The NR simulations reveal unprecedented detail of the earth’s atmosphere, and the then outgoing Editor-in-Chief of AGU JAMES commended the project as one of “stunning ambitions,” enabled by computational capacity at scale. The project also won the 2020 HPCwire Readers Choice Award for Best Use of HPC in Physical Sciences. A set of two NR seasonal simulations have been completed, one corresponding to the northern hemispheric winter months (NDJF) and the other for the North Atlantic tropical cyclone season (ASO). The project used the Summit supercomputer at the Oak Ridge Leadership Computing Facility (OLCF). The simulations were facilitated with an INCITE award from the US Department of Energy Office of Science. For the first seasonal run of four months (NDJF), the hydrostatic IFS model was initialized at 00Z on 1 November 2018. The NR for the TC season (AS) was initialized at 00Z on 1 August 2019. The NR simulations were constrained only by sea surface temperatures (SST) at the lower boundary. The IFS output was saved every 3 hours. After feedback and interest from the scientific community, the simulations were rerun for four specific extreme events, with output every 15 minutes. The special cases include a tropical cycle and three severe storm events over the continental USA.
Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.
In this summary paper, the effects of initial deformations on the buckling and postbuckling characteristics of circular cylindrical shells under hydrostatic pressure is determined in an approximate manner. The influence of initial axisymmetric deformations is stressed. Also, the classical buckling of an axially compressed, noncircular (oval) cylindrical shell is studied. The results show that the major-minor axis ratio of the cross section has a marked effect on the critical load, and that use of the maximum radius of curvature in the formula for the classical buckling stress of a circular cylindrical shell leads to good results for moderate eccentricities.
Large deflection theory is used to compute buckling loads of simply supported initially perfect cylinders under axial compression, external hydrostatic pressure, and combinations of axial compression and internal or external pressure. Important results are obtained by taking into account prebuckling deformations and stresses induced by edge support. For example, the presence of these deformations and stresses can decrease the axial-compression buckling load of an unpressurized perfect cylinder by 50 percent or more.
The effect of out-of-roundness on the yield collapse strength of ring stiffened circular cylindrical shells under hydrostatic pressure has been determined experimentally. On the basis of the experimental data and theoretical stress distributions, empirical relations have been developed which, in conjunction with the theoretically predictable yield collapse strength of round ring-stiffened circular cylindrical shells, permits the prediction of the yield collapse strength of shells with out-of-roundness.
Lunar chemical composition, hydrostatic equilibrium, thermal history and material evaporation are considered in a study of the internal structure of the moon
Photodisintegration of iron into helium and neutrino annihilation process cause hydrostatic and dynamic instabilities that produce star collapse
An essentially frictionless four-point support system which utilizes bearings and pistons allows for determination of vibration frequencies of large structures. Retardation of vertical or horizontal motion is due to the viscous damping by the hydrostatic pressure of the oil or by adjustment of the gas volume in the accumulator.
Test method detects leaks that occur during hydrostatic pressure testing of welded joints in metal tanks. A strip of aluminum foil and a strip of water-soluble paper are placed over the weld. A voltage applied between the tank wall and the foil strip is monitored to detect a decrease in ohmic resistance caused by water leakage into the paper layer.
Buckling characteristics of eccentrically stiffened cylindrical shells under hydrostatic pressure loading and vibration, buckling and flutter characteristics of stiffened flat plates
Leaks in thin wall tubing joints are sealed by cleaning and silver plating the hot gas side of the joint in the leakage area. The pressure differential across the silver during hydrostatic test and subsequent use forces the ductile silver into the leak area and seals it.
Design of precision multimode telescope mounting structure with hydrostatic bearings