Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Averages”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An automated and portable method for selecting an optimal GPU frequency

Power consumption poses a significant challenge in current and emerging graphics processing unit (GPU) enabled high-performance computing systems. In modern GPUs, dynamic voltage frequency scaling (DVFS) appears to be a reliable control to regulate power consumption and performance. However, the DVFS design space is large - hence, brute-force approaches are infeasible to select the optimal frequency. Furthermore, no single frequency can be universally optimal for applications with varying computational intensities. Thus, the application's complexity and the availability of a wide range of frequency settings are a challenge in selecting the optimal frequency configuration for a given GPU workload. To that end, this paper proposes a systematic approach that consists of three steps. The feature characterization study identifies the fine-grain GPU utilization metrics that influence the power consumption and execution time of a given workload. To understand the performance, power, and energy consumption behaviors of a workload across GPU's DVFS design space, we derived analytical power and performance models using the identified fine-grain features. Here, it is shown that the same set of GPU utilization metrics can estimate both the power consumption and execution time while being agnostic of changes to frequency and input sizes. Applying a power control with the single objective of reducing power may cause performance degradation, leading to more energy consumption. A multi-objective approach is proposed to select the optimal GPU DVFS configuration for a workload that reduces power consumption with negligible degradation in performance. The evaluation was conducted using SPEC ACCEL benchmarks and three real applications - NAMD LAMMPS, and LSTM on NVIDIA GV100, GA100, and AMD MI210 GPUs. On average, real applications showed 29.6% energy savings with a performance loss of 5.2% on GA100 and 22.6% energy savings with a performance loss of 4.7% on GV100. Moreover, the proposed models are portable to real applications, GPU architectures, and vendors, and require metric collection at only the default frequency rather than all supported DVFS configurations. Additionally, we conducted a comparison between our models and the GPU assembly instructions (PTX)-based static models. The results revealed a significant reduction in the average error rates, with a decrease from 19.7% to 3.1% for power models and from 29.4% to 5.2% for performance models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Chemical characteristics of iron meteorite parent bodies

The projected relative abundances of the highly siderophile elements (HSE; Re, Os, Ir, Ru, Pt, and Pd) for bulk parent bodies of 10 magmatic iron meteorite groups/grouplet (IC, IIAB, IIC, IID, IIF, IIIAB, IIIF, IVA, IVB, and South Byron Trio) are broadly similar and show no resolvable differences between noncarbonaceous (NC) and carbonaceous (CC) genetic heritage. The processes driving genetic isotopic heterogeneity in the early Solar System, therefore, evidently did not leave discernable chemical fingerprints with respect to HSE relative abundances on the bulk planetesimal scale. By contrast, the absolute abundances of HSE projected for parent body cores, which reflect core size, are more variable and, on average, higher in CC bodies compared to NC bodies. Overall, bulk core chemical compositions, as well as core size, are linked to the distribution of Fe within a parent body, which is controlled by its oxidation state. The CC parent bodies are constrained to have formed under heterogeneous oxidizing conditions which were, on average, more oxidizing than those of the NC environment.

Hilton, Connor D.↗
DOE OSTI2020

Zeolite-templated carbons as effective sorbents to remove methylsiloxanes and derivatives: A computational screening

Though widely used in our daily lives, volatile methylsiloxanes and derivatives are emerging contaminants and becoming a high-priority environment and public health concern. Developing effective sorbent materials can remove siloxanes in a cost-effective manner. Herein, by means of Grand Canonical Monte Carlo (GCMC) simulations, we evaluated the potentials of the recently proposed 68 stable zeolite-templated carbons (ZTCs) for the removal of four linear methylsiloxanes and derivatives as well as two cyclic methylsiloxanes by the calculated average loading and average adsorption energy values. Four ZTCs, namely ISV, FAU1, FAU3, and H8326836, were identified with the top 50% adsorption performance toward all the six targeted contaminants, which outperform activated carbons. Further first principles computations revealed that steric hindrance, electrostatic interactions (further enhanced by charge transfer), and CH-π interactions account for the outstanding adsorption performance of these ZTCs. This work provides a quick procedure to computationally screen promising ZTCs for siloxane removal, and help guide future experimental and theoretical investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Quantifying the influence of different biocrust community states and their responses to warming temperatures on soil biogeochemistry in field and mesocosm studies

Biocrusts influence soil biogeochemistry by fixing carbon (C) and nitrogen (N) and through leachate inputs to soils. Functional rates can vary among biocrust community states and in response to edaphic properties, heterotrophic microbial activity, and global change. Using soils and biocrusts from the Colorado Plateau, Utah, USA, we aimed to quantify the influence of early-successional (ES) and late-successional (LS) biocrust community states on soil biogeochemistry. In a field setting, we found soil was less “fertile” under ES than LS biocrust, but ES biocrust had a relative influence 1.3 times greater than LS biocrust on soil fertility. Leachate collected from LS biocrust had, on average, 6 times more organic C and 1.7 times more dissolved N than ES biocrust, but concentrations of phosphorus (P) and inorganic N did not differ among the two biocrust types. Here, to disentangle influences of biocrusts and soil properties on biogeochemical pools, we constructed mesocosms from homogenized soil and left the surface bare or covered with ES or LS biocrust, before assignment to ambient or warmed (+5 °C) temperature treatments for 3 months. Multivariate biogeochemical properties differed among cover types, yet all exhibited losses of P, N, and organic C and nearly half of the biogeochemical variables considered did not differ among cover types. Mesocosms with LS biocrust retained more dissolved N, supported 8 additional, significant correlations among biogeochemical pools of the biocrust and mineral soil layer on average, and lost fewer of these correlations under warming. Overall, while soils under LS biocrusts were more fertile (i.e., had higher nutrient concentrations) than under ES, we did not find evidence implicating leachate as the primary driver of this difference. Biocrust influences on soil fertility were greater when mineral soil nutrients were in lower concentrations, highlighting the value of even incipient biocrusts for dryland functioning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

How deep is your soil? Quantifying and spatially analyzing understudied deep soil in the United States

Deep soil is largely understudied and important in understanding biogeochemical processes in soil. Here, understudied soil is defined as the difference between soil studied to a known depth and the estimated bedrock depth. To understand more about deep soil, the understudied soil in the US was quantified and spatially analyzed using soil survey data and model estimates of bedrock depth. An equation was derived to find understudied soil using the dataset parameters “max lower depth studied”, “depth to bedrock”, and “likelihood of bedrock in the top 200 cm”. The survey data and bedrock model revealed that soil has been studied to an average depth of 1-2 meters, and the average depth to bedrock is 20 meters. Soil data density in the soil surveys was greatest in the West Coast, Midwest, and areas historically managed for agricultural, while the non-contiguous US and interior West were underrepresented. The soil had been studied deeper than the estimated soil depth in 455 out of 56,889 observation points concentrated in Alaska, California, Texas, Florida, Puerto Rico, and the US Virgin Islands. To understand the diversity and any taxonomic bias of the global soil data available, soil order was compared to US-based National Resource Conservation Service percentages and it was found that Oxisols, Alfisols, Ultisols, Andisols, and Histosols were overrepresented while Gelisols, Aridisols, Vertisols, Entisols, and Spodosols are underrepresented. Soil depth is important in exploring the complexity of biogeochemical processes that take place in soil.

Bedrock↗
DOE OSTI2025

Evolving physical and mechanical rock properties during exposure at Earth's surface

Rocks at Earth's surface experience progressive damage due to exposure to environmental conditions. Although environmental damage to rocks is recognized as being time-dependent, our quantitative understanding of its effects on physical and mechanical properties remains limited by field and laboratory studies constrained to short timescales. In this study, we aim to rigorously quantify how these properties evolve under longer natural environmental exposure spanning from the Holocene to the late Pleistocene. We conducted a ‘time series’ analysis of key physical and mechanical rock properties using multiple interior samples from 10 cogenetic granitoid boulders exposed for different amounts of time at two sites in Eastern California, USA. We find that over the late Pleistocene (~30 ka), rocks exhibit porosity increases of 40 ± 4 % and seismic velocity decreases of 7 ± 1 % at both sites. Similarly, over the same time period, rock tensile and compressive strengths declined by an average of 11 ± 4 %, and rock permeability increased on average by 104 ± 50 %. Microscopic analysis of thin sections of the same rocks indicates that microfracture density and intensity also progressively increased over the same timescales. These results suggest that time-dependent fracture propagation due to environmental exposure is likely a central contributor to the observed changes in rock properties. This study provides quantitative documentation that pronounced changes in the physical and mechanical properties of rock progressively occur within thousands of years of exposure.

Critical zone↗
DOE OSTI2020

Consensus Weighting of a Multi-Agent System for Load Shedding

An agent-based scheme is proposed for distributed underfrequency load shedding (UFLS). The key concept is a consensus weighting protocol (CWP) for agents to reach an agreement. Unlike the well adopted average-consensus protocol (ACP), the proposed CWP enables each agent to converge to its weighted portion of the sum of all initial values. Thus, the proposed CWP is more suitable for the UFLS application, as it allows load buses with a higher level of loading to shed more, rather than reducing by the same average value over all load buses as required by the ACP. Proof of convergence for the proposed CWP is derived and presented in this paper, together with the constraints. For implementation of the multi-agent system (MAS), the monitoring, estimation, and distribution steps are developed. Two study cases and the simulation results are provided to validate the performance of the proposed agent-based UFLS scheme.

Xie, Jing↗
DOE OSTI2021

Energy use of piperazine with the advanced stripper from pilot plant testing

Here, piperazine with the Advanced Stripper (PZAS™) was tested at the National Carbon Capture Center (NCCC) at coal and NGCC conditions. Energy performance was measured with accounting for heat loss. The net heat duty was found to be identical at 2.45 GJ/tonne for both 4 and 12% CO 2 at 90% CO 2 removal. The average measured heat loss was 0.38 GJ/tonne and represented 30% of the gross heat rate. For PZAS, the warm rich piping before the control valve is the largest in surface area. Heat loss at this location increased heat duty by 80% of the heat loss value. Therefore, the total heat duty was decreased by 80% of the heat loss. The Independence™ model was validated using 112 steady-state pilot plant runs and predicted average net heat duty with an error of 9% and lean loading with an error of 10%. The model was used to simulate process impacts of heat loss for the PZAS and simple stripper configurations. For both strippers, heat loss was found to impact heat duty at locations other than just the steam heater. While heat loss in the steam heater was the most important for PZAS, heat loss in the column was the most important for the simple stripper due to a large top-side pinch and CO 2 reabsorption at the same location. CO 2 reabsorption in PZAS occurred only at high heat loss values (~0.14 GJ/hr), compared to low heat loss values (~ 0.03 GJ/hr) in the PZ-simple stripper. Excessive packing in the PZ-simple stripper caused a large top-side pinch at 0.23 mol/mol lean loading, making the column very sensitive to heat loss. Heat loss in the PZ-simple stripper increased heat duty by 120–150% at typical heat loss values and lean loading of 0.23 mol/mol. Heat loss in the MEA-simple stripper only increased heat duty by 1–13% in a typical lean loading range of 0.1–0.2 mol/mol. A large bottom side pinch at low lean loading for 7 m MEA made the effect of heat loss less pronounced.

42 ENGINEERING↗
DOE OSTI2022

CO 2 Pipeline risk assessment and comparison for the midcontinent United States

For this work, a comprehensive quantitative risk assessment for the construction and operation of CO 2 transportation networks considered for the Midcontinent United States was conducted. The results showed risks associated with CO 2 pipelines were significantly less than those of other pipeline types. The assessment used four conceptual pipelines of different lengths to discuss risks operators may see. The assessment evaluated the risk associated with construction and operation using data from the US Occupational Safety Health Administration to determine the risk of injury or death for pipeline workers and data from the US Pipeline and Hazardous Materials Safety Administration for CO 2 , natural gas distribution, natural gas transmission/gathering, and non-CO 2 hazardous liquid pipelines to develop quantitative likelihood and severity values leading to risk values. The data for the assessment covered incidents from 2010 to 2017 for CO 2 pipelines. The average risk for construction and 30 years of operation for four CO 2 pipeline configurations ranging between 79 and 1,546 miles in length was found. The construction and operational risk averaged between $\$1,400,521$ (approximately $\$0.02$/tonne of CO 2 ) for the shorter pipeline (79 miles) and $\$27,481,939$ (approximately $\$0.10$/tonne of CO 2 ) for a longer pipeline (1,546 miles). The largest risks of fatality for CO 2 pipelines comes from vehicle transport. The largest operational risk to the pipeline was due to leakage. Public pipeline opposition is also a significant risk; it was not quantified but is addressed.

03 NATURAL GAS↗
DOE OSTI2023

Measured CO 2 sorption isotherms with 25 Bakken Petroleum System rock samples from the Lower and Upper Shales, Middle Bakken, and Three Forks formations

Isotherms were measured with 25 Bakken Petroleum System (BPS) rock samples using a magnetic suspension balance at reservoir conditions of 110 °C and pressures up to 345 bar. Samples from four wells represented the Middle Bakken (MB) and Three Forks (TF) production zones, and the Upper Bakken Shale (UBS) and Lower Bakken Shale (LBS) source shales. UBS and LBS CO 2 sorption capacities at 345 bar averaged 12–14 g/kg rock compared to 1.3 to 2.5 g/kg rock for MB and TF. Increasing pressure from 150 to 345 bar resulted in only small increases in CO 2 sorption for all lithofacies. For the UBS and LBS source shales, the higher average total organic carbon (TOC) values (14.7 and 11.4 wt.%) and clay values (25 and 28 wt.%) corresponded to their higher CO 2 sorption values while for the MB and TF samples, lower TOCs (0.3 and 0.5 wt.%) and clay values (8.2 and 18.2 wt.%) corresponded to lower CO 2 sorption values. Here, when all 25 samples (from all lithofacies) were evaluated as a group, sorption capacities were highly correlated with TOC (Spearman's rho values +0.84 to +0.85), with rock density (rho values -0.75 to -0.76), and clay content (rho values +0.66). However, correlations with porosity and thermal maturity (Tmax) were weak. Measured sorption capacities and rock porosities showed that sorption accounted for about 10-times as much CO 2 storage as the available pore space for the UBS and LBS lithofacies. In contrast, the pore space dominated the storage resource of the TF and MB lithofacies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗
DOE OSTI2020

Elucidating the Impact of Flow on Material-Sensitive Critical Heat Flux and Boiling Heat Transfer Coefficients: An Experimental Study with Various Materials

In this work, steady-state internal flow boiling critical heat flux (CHF) experiments were carried out at 10°C inlet subcooling and atmospheric pressure (84 kPa at an altitude of about 5300 feet) with mass flux ranging from 200 to 2,000 kg/m 2 -s on seven tube materials, including the accident-tolerant fuel cladding candidate FeCrAl (C26M) alloy. At a relatively low flow rate, the sensitivity of CHF to different materials was appreciable (the average difference among the tested materials was 19% at 200 kg/m 2 -s). As the absolute CHF differences stayed constant for higher flow rates, the increased flow rate led to reduced relative CHF differences among the tested materials (average relative difference of 15% at 1,000 kg/m 2 -s). FeCrAl alloy claddings that were oxidized under simulated Light Water Reactor (LWR) conditions yielded CHF comparable to as-received FeCrAl alloy cladding, with less than 1.8% difference. Experimental CHF results imply marginal material sensitivity impacts on the departure from nucleate boiling ratio in steady-state nuclear reactor cores. Yet, for accident progression analysis, material sensitivity may hold importance, as boiling behavior variation at low mass flux is appreciable. Macroscopic material parameters (surface wettability, surface roughness, and thermal properties) do not provide a suitable explanation for observed material-sensitive CHFs. Yet, surface roughness and thermal effusivity revealed power and linear relations, respectively, to the nucleate boiling heat transfer coefficients at low mass flux. At high mass flux, material sensitivity on nucleate boiling heat transfer significantly decreased due to the increased shear stress on boiling surface overwhelming material-related boiling implications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Numerical Simulations of High Prandtl Number Liquid Jets Impinging on a Flat Plate

In this work, 3D simulations of oil jets impinging on a flat, heated wall are presented. The numerical setup uses the Volume of Fluid (VoF) method to model the two-phase flow. A careful grid definition across the liquid film, along with the use of the Conjugate Heat Transfer (CHT) approach allowed local heat transfer to be solved with fine resolution at the wall. Variations of liquid flow rate, liquid temperature and surface temperature allow to cover a wide range of local Reynolds and Prandtl numbers (226 < Re < 2850, 77 < Pr < 161). Resulting surface-averaged heat transfer compares very well with experimental measurements conducted in a previous study. In-depth analysis of the flow has identified expected features from the literature. In particular, the impact of jet axial velocity profiles on the heat transfer distribution in the stagnation zone was clearly stated. The increase in heat transfer when warming the liquid film was also reproduced and explained by a decrease in oil viscosity and an increase in film velocity. All those effects were taken into account in correlations for stagnation and local values of Nusselt number. A grid sensitivity study was also conducted, showing that if the grid solving the thermal boundary layer in the stagnation zone can be coarsened without impacting local and surface-averaged predictions of heat transfer, a minimum resolution (2 to 3 cells) within the thermal boundary layer is however required for an accurate prediction of heat transfer.

conjugate heat transfer↗
DOE OSTI2021

Shear assisted processing and extrusion of enhanced strength aluminum alloy tubing

Hollow extrusions are used in many industries to make strong, lightweight, and affordable structures. Shear Assisted Processing and Extrusion (ShAPE) enables the extrusion of many alloys, often with noteworthy properties. In this investigation, ShAPE was used to extrude 1 and 2 mm wall tubes of aluminum alloy 6063 measuring 12 mm in diameter at extrusion speeds up to 3.8 m/min. This speed is 10 and 15 times faster than has been reported for any alloy using ShAPE and friction stir back extrusion, respectively. Increasing the extrusion speed from 0.7 to 3.8 m/s resulted in using 68% less process energy at steady state. Microstructural analysis and tensile testing were performed on as-extruded tubes and tubes with a direct T5 heat treatment. Grain refinement from 65 µm to below 20 µm was observed. Growth of macroscale Mg 2 Si strengthening precipitates was not observed. ShAPE processing resulted in recrystallization, and extrusion speed was shown to influence final extrudate texture. Here, transmission electron microscopy analysis revealed that the as-extruded microstructure was free of nanoscale ß" strengthening precipitates. Nanoscale ß" was prominent in direct-aged (T5) ShAPE tubes, similar to a slightly over-aged conventional T6 microstructure. As-extruded tubes had ultimate tensile strengths on par with conventional T5 extrusions and approximately double the total elongation. Tubes that underwent a T5 heat treatment had yield and ultimate strengths averaging 198 and 234 MPa, respectively, with an average total elongation of 11%. The tensile test results are comparable to a conventional T6 heat treatment, but without the solutionization heat treatment step. Mechanical properties either remained flat or improved with increased extrusion speed, which suggests that ShAPE can be scaled up to an industrial process to create energy- and economically efficient high-strength extrusions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermal power plant upgrade via a rotating detonation combustor and retrofitted turbine with optimized endwalls

Over the past decade, pressure gain combustion research has promised over 10 percentage-points of increase in power plant thermal efficiency. Alas, to realize such potential gain, one must effectively couple the turbine with the detonation combustor, whose exhaust conditions differ substantially from current state of the art gas turbines. This paper presents a modeling approach that enables a superior thermodynamic cycle with a rotating detonation combustor and a retrofitted gas turbine by including a diffuser downstream of the combustor and by contouring the turbine endwall, while preserving the airfoil geometry. We propose a multi-step optimization strategy, parametrizing the endwall geometry with a few control points, without altering the airfoil geometry, and with the stage turbine efficiency as objective function. In a first step the turbine performance is assessed with steady inlet conditions by solving the steady three-dimensional Reynolds-Averaged Navier-Stokes equations. In a second step, the inlet conditions are unsteady, as predicted from a detonation combustor and a diffuser, and three-dimensional full unsteady simulations are performed with an unsteady Reynolds-Averaged Navier-Stokes solver. By altering the vane endwall, the steady optimization yielded an efficiency increase of 12% relative to the baseline, while the unsteady optimization resulted in 21% increase compared to the datum turbine. Finally, a full engine analysis demonstrated the superiority of pressure gain combustion which included a realistic thermodynamic cycle of the combustor, the diffuser and the optimized turbine components.

42 ENGINEERING↗
DOE OSTI2022

Automated bubble analysis of high-speed subcooled flow boiling images using U-net transfer learning and global optical flow

Capturing and analyzing the bubble dynamics is crucial to improving the understanding of boiling heat transfer mechanisms and predicting boiling heat transfer coefficient and boiling crisis. High speed video (HSV) imaging has been used for decades towards this end. Still, there is no universal approach to quantitatively analyze bubble dynamics from HSV images. In this study, we propose a data-driven post-processing approach to segment, track, and identify wall-attached vapor bubbles from HSV images of the boiling process in subcooled flow conditions. Firstly, we employ a transfer learning framework with a U-Net-based convolution neural network (CNN) architecture to detect and segment bubbles in HSV images of diverse contrast and surface texture using very little data (e.g., 10 images) for training. Then, we evaluate the trained CNN model with 100 ground-truth images, and the validation results show that the model accuracy and precision in detecting the optical footprint of bubbles are higher than 90%. Finally, we suggest a criterion to identify a condensing bubble based on the divergence of the bubble displacement, which is calculated from sequential segmented bubble images using a global optical flow code. Using this combination of machine learning and optical flow, we can identify nucleation sites and track the growth of bubbles nucleating at each site to quantify nucleation site density, nucleation frequency, and other fundamental boiling parameters. The proposed system is validated using results obtained on a special heater, which enables both infrared (IR) thermometry and HSV imaging on a metallic surface. We compare the fundamental boiling parameters obtained by the two different diagnostics. The results show good agreement. In conclusion, the difference between the measurements of nucleation site density, averaged nucleation frequency, and averaged growth time performed with the two techniques is always within ± 20% and mostly ± 10% of the values measured with IR thermometry.

42 ENGINEERING↗
DOE OSTI2020

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Detection of surface water temperature variations of Mongolian lakes benefiting from the spatially and temporally gap-filled MODIS data

Lakes provide critical water resources for human activities and ecosystems, particularly in the Mongolian Plateau (MP), which is characterized by a dry climate and a harsh environment. As a region that is sensitive to anthropogenic warming, tracking lake surface water temperature (LSWT) changes in Mongolian lakes is crucial for understanding the consequences of a warming climate on lake ecosystems. However, the long-term monitoring of LSWT is restricted by the spatiotemporal gaps in the raw imagery of remote sensing-based land surface temperature (LST), e.g., the commonly used Moderate Resolution Imaging Spectroradiometer (MODIS) LST products. This study applied an improved gap-filling method by utilizing the discrete cosine transform-based penalized least squares (DCT-PLS) strategy in the spatial domain combined with the linear interpolation (LI) algorithm in the temporal domain. The method was applied to fill gaps in the LSWT imagery of 12 representative lakes across MP. The randomly sampled high-quality MODIS LSWT values in the spatial and temporal domains were excavated as false data gaps and considered “virtual true” validation datasets. The spatial validation results showed that the estimated LSWT for all the lake cases were comparable with the “virtual true” LSWT values, with the average values of the coefficient of determination, mean absolute error, mean square error, and root mean square error being 0.98, 0.38 °C, 0.45 °C, and 0.59 °C, respectively. Meanwhile, the error of nighttime LSWT results was relatively lower than that of daytime LSWT. For temporal interpolation validation, the LI algorithm exhibited relatively better performance and could more objectively indicate the variation in LSWT. Benefiting from the spatially and temporally well-constrained data, we analyzed the interannual and intra-annual change characteristics of the LSWTs of the 12 lakes. The long-term variations of annual and seasonal mean LSWTs in the 12 selected lakes exhibited no evident trends in 2000–2020, while presented apparent interannual fluctuations. The slight changes in the average LSWTs of the 12 selected lakes were in excellent synchronization with the surrounding LST derived from the reanalysis datasets, confirming the widely reported phenomenon of “global warming hiatus” that occurred in the early 21st century. This study improves the understanding of the LSWT variations in Mongolian lakes in response to global climate change. It has the potential to provide an effective approach for monitoring LSWT changes in other large-scale studies.

54 ENVIRONMENTAL SCIENCES↗