Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Volumetric constraint”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fourier method for 3-dimensional data fusion of X-ray Computed Tomography and ultrasound

X-ray Computed Tomography (CT) is essential for nondestructive inspection of many manufactured components but is susceptible to various forms of artifacts and noise. In particular, multi-detector row cone-beam CT systems can suffer from so called “cone-beam” artifacts and partial volume effects, particularly on planar edges at the periphery of the field of view in the cone angle dimension. Uni-directional ultrasonic testing methods generally have extremely accurate in-plane depth resolution but poorer lateral resolution due to physical and geometric constraints. This paper presents a novel technique called Computed Tomographic Fusion (CT-F) which uses three-dimensional Fourier filtering to combine x-ray cone-beam CT reconstructions and ultrasound data in the frequency space. The result is a single image with improved accuracy and quality. CT-F mitigates artifacts while allowing for rapid, accurate characterization of large three-dimensional CT volumes. Finally, using simulations and experiments, we demonstrate artifact reduction and edge contrast improvement in volumetric reconstructions.

42 ENGINEERING↗
DOE OSTI2022

Multiscale plasticity of geomaterials predicted via constrained optimization‐based granular micromechanics

Abstract A general framework to derive nonlinear elastic and elastoplastic material models from granular micromechanics is proposed, where a constraint‐based variational structure is introduced to classical grain contact‐based homogenization methods of hyperelasticity. Like the classical hyperelastic methods, reference solutions for closed‐form hyperelastic material models are analytically derived from the grain‐scale contact mechanics. However, unlike prior methods, the proposed homogenization framework defines closed‐form hyperelastoplastic material models that extend multiscale variational methods to granular plasticity. The proposed framework is used to develop novel granular micromechanics‐based macroscopic models for a Mises type solid, Drucker–Prager type plasticity, and grain‐contact cohesive‐debonding with a deviatorically and volumetrically coupled nonlinearly elastic response. Macroscopic plastic parameters and yield criteria are explicitly related to their microscale counterparts, for example, the friction coefficient governing intergranular slip. Numerical examples and comparison to measurements from the literature, including triaxial compaction of concrete, are provided to investigate model predictions and demonstrate calibration to experimental data.

Bryant, E. C.↗
DOE OSTI2025

Radiation damages the silicates present in polyphasic mineral aggregates causing concrete’s degradation

While many U.S. nuclear power plants have submitted Subsequent License Renewal Application to operate beyond 60 years, others are already considering Operations Beyond Eighty years. In such cases, concrete biological shields are exposed to neutron and gamma radiation exceeding prescribed thresholds. Radiation-induced volumetric expansion (RIVE), extensively studied in single crystals, may also contribute to the degradation of polycrystalline aggregates. Since minerals differ in atomic structure and chemistry, radiation can affect them in distinct ways. This study examines quartzite, marble, and limestone to evaluate how irradiation affects their physical attributes and chemical reactivity. Results show crystalline silicates experience significant RIVE damage and enhanced reactivity in alkaline solutions compared to non-irradiated phases. Enhanced intra- and inter-granular dissolution could compromise aggregate integrity. An empirical correlation links silicate dissolution rate to atomic constraints (density, rigidity) and radiation dose, providing a predictive framework for estimating changes in silicate aggregate properties within radiation-exposed concrete.

Bouissonnié, Arnaud [Univ. of California, Los Ange↗
DOE OSTI2024

A Bulk versus Nanoscale Hydrogen Storage Paradox Revealed by Material-System Co-Design

Metal hydrides are serious contenders for materials-based hydrogen storage to overcome constraints associated with compressed or liquefied H 2 . Their ultimate performance is usually evaluated using intrinsic material properties without considering a systems design perspective. An illustrative case with startling implications is (LiNH 2 +2LiH). Using models that simulate the storage system and associated fuel cell of a light-duty vehicle (LDV), the performance of the bulk hydrides is compared with a nanoscaled version in porous carbon (PC), (LiNH 2 +2LiH)@(6-nm PC). Using experimental material properties, the simulations show that (LiNH 2 +2LiH)@(6-nm PC) counterintuitively has higher usable gravimetric and volumetric capacities than the bulk counterpart on a system basis despite having lower capacities on a materials-only basis. Nanoscaling increases the thermal conductivity and lowers the desorption enthalpy, which consequently increases heat management efficiency. In a simulated drive cycle for fuel cell-powered LDV, the fuel cell is inoperable using bulk (LiNH 2 +2LiH) as the storage material but completes the drive cycle using the nanoscale material. Further, these results challenge the notion that nanoscaling incurs mass and volume penalties. Instead, the synergistic nanoporous host-hydride interaction can favorably modulate chemical and heat transfer properties. Moreover, a co-design approach considering application-specific tradeoffs is essential to accurately assess a material's potential for real-world hydrogen storage.

08 HYDROGEN↗
NASA NTRS2008

Adaptation of the Skeletal System during Long-duration Spaceflight

This review will highlight evidence from crew members flown on space missions greater than 90 days to suggest that the adaptations of the skeletal system to mechanical unloading may predispose crew members to an accelerated onset of osteoporosis after return to Earth. By definition, osteoporosis is a skeletal disorder - characterized by low bone mineral density and structural deterioration - that reduces the ability of bones to resist fracture under the loading of normal daily activities. Involutional or agerelated osteoporosis is readily recognized as a syndrome afflicting the elderly population because of the insipid and asymptomatic nature of bone loss that does not typically manifest as fractures until after age approximately 60. It is not the thesis of this review to suggest that spaceflight-induced bone loss is similar to bone loss induced by metabolic bone disease; rather this review draws parallels between the rapid and earlier loss in females that occurs with menopause and the rapid bone loss in middle-aged crew members that occurs with spaceflight unloading and how the cumulative effects of spaceflight and ageing could be detrimental, particularly if skeletal effects are totally or partially irreversible. In brief, this report will provide detailed evidence that long-duration crew members, exposed to the weightlessness of space for the typical long-duration (4-6 months) mission on Mir or the International Space Station -- 1. Display bone resorption that is aggressive, that targets normally weight-bearing skeletal sites, that is uncoupled to bone formation and that results in areal BMD deficits that can range between 6-20% of preflight BMD; 2. Display compartment-specific declines in volumetric BMD in the proximal femur (a skeletal site of clinical interest) that significantly reduces its compressive and bending strength and which may account for the loss in hip bone strength (i.e., force to failure); 3. Recover BMD over a post-flight time period that exceeds spaceflight exposure but for which the restoration of whole bone strength remains an open issue and may involve structural alteration; and 4. Display risk factors for bone loss -- such as the negative calcium balance and down-regulated calcium-regulating hormones in response to bone atrophy -- that can be compounded by the constraints of conducting mission operations (inability to provide essential nutrients and vitamins). The full characterization of the skeletal response to mechanical unloading in space is not complete. In particular, countermeasures used to date have been inadequate and it is not yet known whether more appropriate countermeasures can prevent the changes in bone that have been found in previous flights, knowledge gaps related to the effects of prolonged (greater than or equal to 6 months) space exposure and to partial gravity environments are substantial, and longitudinal measurements on crew members after spaceflight are required to assess the full impact on skeletal recovery.

Sibonga, Jean D.↗
NASA NTRS2005

High-Power Radar Sounders for the Investigation of Jupiter Icy Moons

The high power and high data rate capability made available by a Prometheus class spacecraft could significantly enhance our ability to probe the subsurface of the planets/moons and asteroid/comets. The main technology development driver for our radar is the proposed Jupiter Icy Moon Orbiter (or JIMO) mission due to its harsh radiation environment. We plan to develop a dual-band radar at 5 and 50 MHz in response to the two major science requirements identified by the JIMO Science Definition Team: studying the near subsurface (less than 2 km) at high resolution and detection of the ice/ocean interface for Europa (depth up to 30 km). The 50-MHz band is necessary to provide high spatial resolution (footprint and depth) as required by the JIMO mission science requirements as currently defined. Our preliminary assessment indicates that the 50-MHz system is not required to be as high-power as the 5-MHz system since it will be more limited by the surface clutter than the Jupiter or galactic background noise. The low frequency band (e.g. 5 MHz), which is the focus of this effort, would be necessary to mitigate the performance risks posed by the unknown subsurface structure both in terms of unknown attenuation due to volumetric scattering and also the detection of the interface through the attenuative transition region at the ice/ocean interface. Additionally, the 5-MHz band is less affected by the surface roughness that can cause loss of coherence and clutter noise. However, since the Signal-to-Noise-Ratio (SNR) of the 5-MHz radar band is reduced due to Jupiter noise when operating in the Jupiter side of the moon, it is necessary to increase the radiated power. Our challenge is to design a high-power HF radar that can hnction in Jupiter's high radiation environment, yet be able to fit into spacecraft resource constraints such as mass and thermal limits. Our effort to develop the JIMO radar sounder will rely on our team's experience with planetary radar sounder design gained during our participation in the MARSIS radar sounder implementation.

Safaeinili, A.↗
NASA NTRS1991

Pictorial communication in virtual and real environments

Papers about the communication between human users and machines in real and synthetic environments are presented. Individual topics addressed include: pictorial communication, distortions in memory for visual displays, cartography and map displays, efficiency of graphical perception, volumetric visualization of 3D data, spatial displays to increase pilot situational awareness, teleoperation of land vehicles, computer graphics system for visualizing spacecraft in orbit, visual display aid for orbital maneuvering, multiaxis control in telemanipulation and vehicle guidance, visual enhancements in pick-and-place tasks, target axis effects under transformed visual-motor mappings, adapting to variable prismatic displacement. Also discussed are: spatial vision within egocentric and exocentric frames of reference, sensory conflict in motion sickness, interactions of form and orientation, perception of geometrical structure from congruence, prediction of three-dimensionality across continuous surfaces, effects of viewpoint in the virtual space of pictures, visual slant underestimation, spatial constraints of stereopsis in video displays, stereoscopic stance perception, paradoxical monocular stereopsis and perspective vergence. (No individual items are abstracted in this volume)

Ellis, Stephen R.↗
NASA NTRS1982

Effects of subsurface volume scattering on the lunar microwave brightness temperature spectrum

The effects of volumetric scattering on the lunar microwave brightness temperature are examined for a broad range of feasible lunar rock population distributions. Mie-scattering phase functions and the radiative transfer method are utilized. Surveyor and Apollo data relevant to lunar rock size distributions are discussed, and parameters are chosen for nine scattering models which liberally cover the range of studied rock population distributions. Scattering model brightness temperature predictions are analyzed in terms of the lunar disk center emission averaged over a lunation for wavelengths of 3-30 cm. The effects of scattering on the amplitude of disk center brightness temperature variations and resultant deductions of regolith electrical loss are examined. Constraints on the global scale variability of subsurface scatterers imposed by microwave brightness temperature maps are considered.

Keihm, S. J.↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2021

Beyond Taylor’s hypothesis: a novel volumetric reconstruction of velocity and density fields for variable‑density and shear flows

Abstract This work presents a novel numerical procedure for reconstructing volumetric density and velocity fields from planar laser-induced fluorescence (PLIF) and stereoscopic particle image velocimetry (SPIV) data. This new method is theoretically and practically demonstrated to provide more accurate 3D vortical structures and density fields in high shear flows than reconstruction methods based on the mean convective velocity. While Taylor’s hypothesis of frozen turbulence is commonly applied by using the local mean streamwise velocity, the proposed algorithm uses the measured local instantaneous velocity for data convection. It consists of a step-by-step reconstruction based on a mixed Lagrangian–Eulerian solver that includes the 3D interpolation of scattered flow data and that relaxes the Taylor’s hypothesis by iterative enforcement of the incompressibility constraint on the velocity field. This methodology provides 3D fields with temporal resolution, spatial resolution, and accuracy comparable to that of real 3D snapshots, thus providing a practical alternative to tomographic measurements. The procedure is validated using numerical data of the constant-density channel flow available on the Johns Hopkins University Turbulence Database (JHTDB), showing the accurate reconstruction of the 3D velocity field. The algorithm is applied to an experimental dataset of PLIF and SPIV measurements of a variable-density jet flow, demonstrating its capability to provide 3D velocity and density fields that are more consistent with the Navier–Stokes equations compared to the mean flow convective method. Graphic abstract

42 ENGINEERING↗
DOE OSTI2026

Memory-Aware External Facelist Calculation: A Data-Parallel Atomic Hash Counting Approach

Unstructured volumetric meshes serve as fundamental data representations in various scientific simulations and analyses. They play a crucial role in representing complex computational domains and are essential for important numerical techniques, such as finite element analysis. Whenever such a mesh is read from a file, streamed in-situ, or generated by algorithms, scientific visualization libraries rely on calculating the external surface of a geometry, named “external facelist”, to produce a polygonal mesh for rendering. Consequently, external facelist calculation has become one of the most widely used algorithms in the scientific visualization domain, necessitating optimal performance. In this paper, we explore relevant work on external facelist calculation algorithms in two common visualization libraries, VTK and Viskores, assess their performance and memory constraints, and introduce a novel memory-aware external facelist calculation algorithm employing an atomic hash counting approach. This algorithm fully leverages Viskores' data-parallel primitive operations, facilitating its execution across diverse many-core architectures. Our algorithm features the lowest memory footprint on the GPU and the second-lowest on the CPU among all evaluated methods, and it also delivers the fastest performance on both CPU and GPU. It has been made available under an open-source license in the VTK and Viskores visualization systems.

Tsalikis, Spiros [Kitware] (ORCID:0000000151137195↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
NASA NTRS1991

Sources of volcanic aerosols: Petrologic and volcanological constraints

Global climatic effects brought about by volcanism are related to the impact of volcanic gases and their derivative aerosols on the atmosphere, rather than the effects of volcanic ash. Evidence from both historic eruptions and polar ice cores indicate that volcanic sulfur gases are the dominant aerosol-forming component, resulting in produciton of a sulfuric acid-rich stratosphere aerosol that can have profound effects on the earth radiation budget over periods of a few years. Due to highly variable sulfur content of different magma types, the climatic effects do not relate simply to total erupted mass. There is a close relationship between volcanic sulfur yield to the atmospheric and hemispheric surface temperature decrease following an eruption, with up to 1 C surface temperature decrease indicated following a major volcanic event such as the 1815 Tambora eruption. While the erupted mass of HCl and HF is equal to or greater than that of sulfur gases in some volcanic events, the halogens do not form known aerosols nor are they abundant in ice core acidity layers. The early removal of halogens from eruption columns occurs by rain flushing and adsorption onto tephra particles, but the fate of halogens in the atmosphere following very large explosive eruptions is unknown. The CO2 flux to the atmosphere from volcanic eruptions is volumetrically one of the most important of the gas species, but owing to the huge size of the atmospheric reservoir of this gas, the volcanic contribution is likely to have negligible effects.

Sigurdsson, Haraldur↗
NASA NTRS1980

A numerical investigation of thermal convection in a heat-generating fluid layer

Finite difference solutions of the equations governing thermal convection driven by uniform volumetric energy sources are presented for two-dimensional flows in a rectangular domain. The boundary conditions are a rigid (i.e, zero slip), zero heat-flux lower surface, rigid adiabatic sides, and either a rigid or free (i.e., zero shear) isothermal upper surface. Computations are carried out for Prandtl numbers from 0.05 to 20 and Rayleigh numbers from 5 x 10 to the 4th to 5 x 10 to the 8th. Nusselt numbers and average temperature profiles within the layer are in good agreement with experimental data for rigid-rigid boundaries. For rigid-free boundaries, Nusselt numbers are larger than in the former case. The structure of the flow and temperature fields in both cases is dominated by rolls, except at larger Rayleigh numbers where large-scale eddy transport occurs. Generally, low velocity upflows over broad regions of the layer are balanced by higher velocity downflows when the flow exhibits a cellular structure. The hydrodynamic constraint at the upper surface and the Prandtl number are found to influence only the detailed nature of flow and temperature fields. No truly steady velocity and temperature fields are found despite the fact that average Nusselt numbers reach steady values.

Emara, A. A.↗
DOE OSTI2020

3D Volumetric Structural Hierarchy Induced by Colloidal Polymerization of a Quantum-Dot Ionic Liquid Monomer Conjugate

Hierarchically structured emissive quantum dot (CdS/CdSe) poly(ionic liquid) composites are synthesized through colloidal polymerization within a lyotropic liquid crystal, yielding a solvent-responsive nanoparticle (NP) composite possessing four ordered structural motifs. The composite is prepared by coassembling an amphiphilic ionic liquid (IL) monomer, 3-decyl-1-vinylimidazolium chloride, [C 10 VIm + ][Cl - ], and IL monomer, 3-(10-mercaptodecyl)-1-vinylimidazolium bromide surface-functionalized quantum dots, QD-[SC 10 VIm + ][Br - ], in 30 wt % water. Photo (UV)-initiated free radical polymerization serves to form a colloidal copolymer that collimates pairs of QDs within the core of large cylinders (average diameter = 21 nm and length = 422 nm), as evidenced by ultrasmall-angle X-ray scattering (USAXS). As revealed by continuous wave terahertz multispectral reconstructive 3-D imaging, the subsurface structure along the cylinder axis undergoes dynamic rearrangement between 2D hexagonal and cubic packing. The polymerized ionic liquid sheath surrounding the QDs is composed of concentric hexagonally perforated lamellae (HPL) with an interlamellar repeat distance of 3 nm, as determined by small-angle X-ray scattering (SAXS). The encapsulated QDs retain crystallinity (TEM and WAXS) and emission characteristics (photoluminescence spectroscopy, lambda(EM) = 638 nm). Phenomenological time-dependent Ginzburg-Landau calculations verify the mesoscale cylinders arise from the topological constraints imposed by the in-plane perforations of the polymerized ionic liquid as it wraps around the hard (unmalleable) spherical QDs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

Emplacement of Volcanic Domes on Venus and Europa

Placing firmer constraints on the emplacement timescales of visible volcanic features is essential to obtaining a better understanding of the resurfacing history of Venus. Fig. 1 shows a Magellan radar image and topography for a putative venusian lava dome. 175 such domes have been identified, having diameters that range from 19 - 94 km, and estimated thicknesses as great as 4 km [1-2]. These domes are thought to be volcanic in origin [3], having formed by the flow of a viscous fluid (i.e., lava) onto the surface. Among the unanswered questions surrounding the formation of Venus steep-sided domes are their emplacement duration, composition, and the rheology of the lava. Rheologically speaking, maintenance of extremely thick, 1-4 km flows necessitates higher viscosity lavas, while the domes' smooth upper surfaces imply the presence of lower viscosity lavas [2-3]. Further, numerous quantitative issues, such as the nature and duration of lava supply, how long the conduit remained open and capable of supplying lava, the volumetric flow rate, and the role of rigid crust in influencing flow and final morphology all have implications for subsurface magma ascent and local surface stress conditions. The surface of Jupiter's icy moon Europa exhibits many putative cryovolcanic constructs [5-7], and previous workers have suggested that domical positive relief features imaged by the Galileo spacecraft may be volcanic in origin [5,7-8] (Fig. 2). Though often smaller than Venus domes, if emplaced as a viscous fluid, formation mechanisms for europan domes may be similar to those of venusian domes [7]. Models for the emplacement of venusian lava domes (e.g. [9-10]) have been previously applied to the formation of putative cryolava domes on Europa [7].

timescales↗
NASA NTRS2008

Structural Definition and Mass Estimation of Lunar Surface Habitats for the Lunar Architecture Team Phase 2 (LAT-2) Study

The Lunar Architecture Team Phase 2 study defined and assessed architecture options for a Lunar Outpost at the Moon's South Pole. The Habitation Focus Element Team was responsible for developing concepts for all of the Habitats and pressurized logistics modules particular to each of the architectures, and defined the shapes, volumes and internal layouts considering human factors, surface operations and safety requirements, as well as Lander mass and volume constraints. The Structures Subsystem Team developed structural concepts, sizing estimates and mass estimates for the primary Habitat structure. In these studies, the primary structure was decomposed into a more detailed list of components to be sized to gain greater insight into concept mass contributors. Structural mass estimates were developed that captured the effect of major design parameters such as internal pressure load. Analytical and empirical equations were developed for each structural component identified. Over 20 different hard-shell, hybrid expandable and inflatable soft-shell Habitat and pressurized logistics module concepts were sized and compared to assess structural performance and efficiency during the study. Habitats were developed in three categories; Mini Habs that are removed from the Lander and placed on the Lunar surface, Monolithic habitats that remain on the Lander, and Habitats that are part of the Mobile Lander system. Each category of Habitat resulted in structural concepts with advantages and disadvantages. The same modular shell components could be used for the Mini Hab concept, maximizing commonality and minimizing development costs. Larger Habitats had higher volumetric mass efficiency and floor area than smaller Habitats (whose mass was dominated by fixed items such as domes and frames). Hybrid and pure expandable Habitat structures were very mass-efficient, but the structures technology is less mature, and the ability to efficiently package and deploy internal subsystems remains an open issue.

Dorsey, John T.↗