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.
This dataset provides an estimation of the surface turbulent fluxes (eddy covariance) calculated with EddyPro v7.0.9 software, 10-min averaging, from raw 20Hz sonic anemometer data.
ABSTRACT Separating earthquake source spectra from propagation effects is challenging. The propagation effect contains a site-dependent term related to the high attenuation of shallow sediments. Neglecting the site-dependent attenuation can cause large biases and scattering in the corner-frequency (fc) estimates, resulting in significant stress-drop deviations. In this study, we investigate shallow attenuation at the LArge-n Seismic Survey in Oklahoma (LASSO) and site-related biases and scattering in source parameter measurements due to simplified attenuation models. We measure the high-frequency spectral decay parameter kappa on the vertical acceleration spectra of regional earthquakes (125 km away). The site-dependent kappa (κ0,acc) suggests that attenuation increases rapidly at shallow depth and is highly site-dependent. 10%–75% of the attenuation is site-dependent for S waves and even larger for P waves. The quality factor for S waves (QS) ranges from 10 to 100 in the upper 400 m. QP for P waves is mostly below 10 within the same depth. The Quaternary sediments tend to be more attenuating (QS<30), but the Permian rocks also can have high attenuation. We demonstrate that using a non-site-dependent attenuation model in single-spectra fitting leads to large scattering in fc estimates among stations with apparent good fits. The apparent fc can significantly deviate when the range of site-dependent kappa is large or with a higher assumed source spectral fall-off rate n. The biases in apparent fc depend on site condition and distance; however, the correlation between fc and these factors might not be obvious, depending on model assumptions. An apparent increase of stress drop with magnitude in a previous study for local microearthquakes (1.3
Abstract Investigations into the role of anthropogenic emissions in the occurrence of extreme weather often use a method that compares simulations of atmospheric climate models run under a factual scenario of historical boundary conditions observed during the period of the event against simulations run under a counterfactual scenario of what those boundary conditions might naturally have been over that same period in the absence of anthropogenic emissions. A particular requirement for this experiment design is an accurate estimation of ocean surface boundary conditions for use by the counterfactual natural simulations. Here we use output from the CMIP5 multi‐climate‐model archive to develop a robust estimate of sea surface temperatures and sea ice conditions for use in counterfactual natural simulations, intended as a benchmark estimate to facilitate comparison across climate models and across studies. This development includes tests to ensure that the final estimate is stable from year‐to‐year and stable against other perturbations to the methodology, as well as consideration of the strengths and weaknesses in comparison to other available attributable warming estimates. While this estimate is tailored specifically for the International CLIVAR C20C+ Detection and Attribution Project, it can be used by related projects as well.
Measuring surface water temperature spatial variability is needed to estimate the interaction between surface water and groundwater, evaluate fish habitat and thermal inertia, and to estimate streamflow frequency and duration. Fiber optic distributed temperature sensing (FO-DTS) has been used in rivers and lakes, providing high-resolution and sensitive temperature monitoring over large temporal and spatial scales. However, in streams with cobbly or bedrock-lined streambeds and variable bathymetry, use of FO-DTS to measure temperature close to the surface water and groundwater interface can be challenging if even feasible. FO-DTS can also be costly, involve difficult installations, and require an advanced understanding of the technology, calibration, and data processing. In this study, we compared FO-DTS stream temperature survey results to an alternative temperature survey method employing a towed transect of high-resolution temperature loggers spaced at 1-m and transported in the stream along the study reach, to measure the spatial distribution of stream-water temperature in East Fork Poplar Creek near Oak Ridge, Tennessee, USA. We assessed the applicability and limitations of the two methods, and quantitatively compared in-situ temperature survey results measured simultaneously with each method. Regression results showed strong temporal and spatial correlation between the two methods. Differences were only elevated near the stream banks in areas that were coincident with correlation slope deviations from unity, which was attributed to shallower water and lower data density. Kriging standard errors were also low at channel center with minor increases near the stream banks. Furthermore, the results suggested that the array of the individual temperature sensors can provide a practical alternative to FO-DTS for thermal characterization of surface water, providing slightly lower spatial and temporal resolution, but with higher accuracy of temperature measurement, with greater simplicity, and with a broader range of conditions where it may be applied.
The last decade has seen a large number of severe heatwaves that were unprecedented in the observational record, highlighting challenges associated with observationally-based statistical quantification of the likelihood and magnitude of future extreme temperatures. An alternative to such probabilistic assessments is identification of upper bounds that quantify the hottest surface air temperatures that can possibly be achieved by the end of the 21st century. Theory, simulations, and observational analyses support the existence of a finite upper bound for surface air temperature; however, estimates for future upper-bound values that are realistic and usable for planning remain unavailable. Here, we combine atmospheric theory with large ensembles of dynamically downscaled projections to estimate historical and end-of-century upper bounds for surface air temperatures. A number of physical mechanisms can influence upper bounds, and at the end of the 21st century, estimates based on mechanisms that yield more moderate upper-bounds produce values around 60∘C for much of the Western United States and in excess of 80∘C for the hottest parts of the domain. Even cooler high-altitude locations have end-of-century upper bounds over 50∘C. Although these upper-bound estimates might seem implausibly large, increases in the upper bounds over the 21st century are similar to increases in dynamically downscaled peak surface temperatures after adjusting those downscaled temperatures to eliminate the possibly biased model trends in surface specific humidity. While upper bound estimates are high relative to historical observations, they nonetheless suggest that heatwave intensity risk is bounded, with uncertainty dominated by projections of surface and upper-level humidity.
Abstract The determination of basic plasma parameters in atmospheric pressure discharges is critical to advancing their use in applications. Atmospheric pressure plasma jets have found use in the fields of medicine, agriculture, material modification and others. Atmospheric pressure plasma jets often generate plasma surface ionization waves (SIW) which interact with and propagate over surfaces. Electrical diagnostics are challenging in SIW due to high collision frequencies and small scale of the plasma discharge. This work employs a passive optical emission line ratio technique to estimate the mean electron energy in SIW over planar dielectric substrates. The method uses an intensity ratio of two helium triplet lines: He( 3 3 S ) at 706.5 nm and He( 3 3 D ) at 587.56 nm. A collisional-radiative model is used to correlate line ratio to mean electron energy and determine dependencies on electron density and He/air gas mixture. Mean electron energies ranging from 3–8 eV are determined in He/air mixtures and are found to remain constant as the surface wave propagates radially. This work provides a 2D, time-resolved, mean electron energy diagnostic for surface ionization wave propagation and validation of numerical modeling in atmospheric pressure systems with spatially varying He/air gas mixtures. The model in question is designed for use with any He line ratio in the n = 3 excitation level.
We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.
In-situ monitoring technologies for laser powder bed fusion (LPBF) additive manufacturing often face one key challenge, extracting the ultrafast melt pool (MP) signatures for understanding the localized part properties. Further, the spatial information of each monitored MP signature is essential for correlating the MP – part property. This spatial information is often unavailable especially from commercial LPBF printers. Many MP monitoring methods have been reported and utilized. However, very few of these have the MP’s spatial information. To overcome this challenge, in this work we report a method for spatially registering the key MP signatures (MP intensity, temperature, and area) to the monitored print parts. The MP signatures are obtained from our coaxial high-speed single-camera based two-wavelength imaging pyrometry (STWIP) system and the MP spatial information is obtained from an off-axis camera system. A machine learning aided image analysis method is employed to retrieve the spatial distribution of MPs within the corresponding part’s coordinates system. Then, the MP signature maps (MPSMs) are reconstructed by mapping the STWIP measured MP signatures to the registered MP coordinates. Further, a long short-term memory (LSTM) neural network is developed for estimating the layer surface topography from the registered MPSMs. The obtained results indicate that the layer surface topography can be more accurately estimated by using MP temperature signature rather than MP intensity and/or area signatures as in common practice. Finally, our developed methods for MP monitoring, registration, and MP-surface topography prediction offer advanced capabilities for the online detection of process anomalies and part defects.
Considering the current booming interest for the large-scale deployment of bifacial photovoltaic modules, the solar industry now requires accurate estimates of broadband surface albedo at high spatial resolution. In this context, the present study evaluates the adequacy and performance over North America of two Moderate Resolution Imaging Spectroradiometer (MODIS) white-sky albedo products (at 500-m and 1-km resolution) and the National Solar Radiation Database (NSRDB) product at 4-km resolution. Two variations of the 500-m MODIS product are also considered: black-sky albedo and all-sky albedo. Albedo observations from 36 radiometric stations during 2011–2015 are analyzed while considering the apparent homogeneity of the surface characteristics over the 4x4 km NSRDB pixel in which they are located. Even at sites where the albedo around the station has been found “homogeneous” in the literature, marked differences are found between the daily observations and the gridded estimates at any spatial resolution. Furthermore, differences in seasonal behavior and between the three different albedo types also impact the accuracy of the albedo estimates, with overtones caused by local specificities and inhomogeneities. All this precludes the desirable evaluation of the local albedo at a specific site of relatively small size compared to its corresponding 4x4-km pixel if only the mean albedo over that pixel is known. Significant discrepancies are also found at snow-impacted sites, most importantly in the case of the NSRDB albedo estimates, which are typically much too high.
Abstract. Permafrost currently stores more than a fourth of global soil carbon. A warming climate makes this carbon increasingly vulnerable to decomposition and release into the atmosphere in the form of greenhouse gases. The resulting climate feedback can be estimated using land surface models, but the high complexity and computational cost of these models make it challenging to use them for estimating uncertainty, exploring novel scenarios, and coupling with other models. We have added a representation of permafrost to the simple, open-source global carbon–climate model Hector, calibrated to be consistent with both historical data and 21st century Earth system model projections of permafrost thaw. We include permafrost as a separate land carbon pool that becomes available for decomposition into both methane (CH4) and carbon dioxide (CO2) once thawed; the thaw rate is controlled by region-specific air temperature increases from a preindustrial baseline. We found that by 2100 thawed permafrost carbon emissions increased Hector’s atmospheric CO2 concentration by 5 %–7 % and the atmospheric CH4 concentration by 7 %–12 %, depending on the future scenario, resulting in 0.2–0.25 ∘C of additional warming over the 21st century. The fraction of thawed permafrost carbon available for decomposition was the most significant parameter controlling the end-of-century temperature change in the model, explaining around 70 % of the temperature variance, and was distantly followed by the initial stock of permafrost carbon, which contributed to about 10 % of the temperature variance. The addition of permafrost in Hector provides a basis for the exploration of a suite of science questions, as Hector can be cheaply run over a wide range of parameter values to explore uncertainty and can be easily coupled with integrated assessment and other human system models to explore the economic consequences of warming from this feedback.
The Surface Atmosphere Integrated Field Laboratory (SAIL) campaign generated an unprecedented, multisensory record of cold-season precipitation in the Upper Colorado River Basin (UCRB). To complement the Surface Quantitative Precipitation Estimation (SQUIRE) snow water equivalent (SWE) product derived from the Colorado State University (CSU) X-band radar, a hydrometeor phase classification data set was developed to systematically identify phase of the near-surface precipitation. The algorithm uses Corrected Moments in Antenna Coordinates (CMAC) fields with semi-supervised classification (SSC) and fuzzy-logic hydrometeor classification (FHC) methods and maps hydrometeor classes to a compact set of phase classes suitable for Quantitative Precipitation Estimation (QPE) and hydrologic applications. We used SQUIRE grid projections to mitigate beam-blockage issues by selecting the near-surface gates for each beam as well as to make it directly usable for the QPE products.
Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.
The structure and evolution of the atmospheric boundary layer (ABL) under clear-sky fair weather conditions over mountainous terrain is dominated by the diurnal cycle of the surface energy balance and thus strongly depends on surface snow cover. We use data from three passive ground-based infrared spectrometers deployed in the East River Valley in Colorado's Rocky Mountains to investigate the response of the thermal ABL structure to changes in surface energy balance during the seasonal transition from low to high snow cover. Temperature profiles were retrieved from the infrared radiances using the optimal estimation physical retrieval Tropospheric Remotely Observed Profiling via Optimal Estimation. A nocturnal surface inversion formed in the valley during clear-sky days, which was subsequently mixed out during daytime with the development of a convective boundary layer when snow cover was low. Over high snow cover, a very shallow convective boundary layer formed, above which the inversion persisted through the daytime hours. We compare these observations to NOAA's operational High-Resolution-Rapid-Refresh model and find large warm biases on clear-sky days resulting from the model's inability to form strong nocturnal inversions and to maintain the stable stratification in the valley during daytime when there was snow on the ground. We suggest several factors contributing to the large model errors. These are (a) the inability of the model to represent well-developed thermally driven flows likely due to the too coarse horizontal grid spacing (3 km), (b) too much convective mixing during daytime, and (c) too strong vertical coupling between the valley atmosphere and the free troposphere.
Abstract Among the subsurface geophysical methods used in the critical zone investigations, induced polarization (IP) shows great vitality thanks to its unique ability to assess porosity via bulk conduction and estimate permeability through surface conduction and/or polarization. However, such an advantageous separation between bulk and surface is mostly implemented by multi‐salinity experiments in the laboratory, which is incredibly difficult to realize in the field. One promising approach to address such an obstinate issue is to gauge the surface conductivity ( σ s ) from the quadrature conductivity ( σ ″) or normalized chargeability ( M n ) with the ratios between ( l = σ ″/ σ s , l mn = M n / σ s ). While these ratios are known not to be universal, the underlying principles are not fully understood and relevant theoretical studies are rare, which makes quantitative IP applications difficult. Here we scrutinize the conduction and polarization mechanisms of geomaterials and pinpoint that the two ratios are inherently functions of salinities and frequencies rather than only determined by the properties of the electrical double layer (EDL), hence representative samples from the investigated field must be calibrated in the laboratory and a characteristic frequency should be chosen for their usage. Besides the macro‐scale ratios l and l mn , we define two micro‐scale ratios χ and χ mn directly from the EDL, such that the new ratios exclude the effect of salinity and frequency and offer the opportunity to characterize and monitor changes of the EDL. Our study demonstrates that the existing macro‐scale ratios converge toward the values of novel micro‐scale ratios at high water salinity.
(Section 5.1) Composites made from Arkema’s Elium® thermoplastic resin and Johns Manville fiberglass were researched during this project for applications in wind blade manufacturing. A techno-economic model was developed to model this wind blade manufacturing process using these materials in place of traditional composites made with thermoset resin. This model was based on manufacturing a 61.5-meter wind blade, which showed a 4.7% reduction in wind blade cost as compared traditional thermoset materials. These cost savings were not from the thermoplastic material costing less than traditional thermoset materials, but rather from decreased capital costs, faster cycle times and reduced energy requirements and labor costs. (Section 5.2) An infusion and curing model was developed for thermoplastic composite wind blades using PAM-RTM. The primary goal was to demonstrate the infusion simulation for the Elium® resin system on a 13-meter wind blade. Additionally, the exotherm temperature was predicted and compared to measurements, which showed model results within 10% of actual measurements. (Section 5.3) Composite laminate panels and composite sandwich panels with a balsa core were produced; specimens were cut and characterized. Similar composite specimens were made with Elium® thermoplastic resin and Hexion thermoset epoxy (RIMR135/RIMH1366) to enable comparisons between these resin systems. The static test methods included: tensile, compression, in-plane shear, interlaminar shear, flexural, sandwich core shear flexure, and single cantilever beam tests for sandwich beams. Fatigue testing at room temperature was completed to composite laminate panels at a stress ratio of R=0.1 and R=10. In addition, fatigue testing to laminate panels was completed at -30°C, and at room temperature after conditioning specimens at 70°C and 90% relative humidity. Overall, mechanical test results from Elium® composites are similar to epoxy composites. (Section 5.4) Elium composite panels were produced with intentional defects such as voids and nonwetting of fibers to begin to understand performance sensitivity to defects. A thermal digital image correlation (TDIC) method provides high spatial resolution strain field at elevated temperatures and can be used to identify defective regions within composite panels. Flexural modulus differences of 21% were seen between defect and non-defect panels. Other Elium® composite panels were forced to be defective by boiling the resin after infusion, which created voids throughout the composite laminate. X-ray computed tomography scanning was used to view the internal structure of the defect panels. Defect panels had a significant reduction in fatigue life as compared to baseline panels produced without intentional defects. (Section 5.5) Lap shear specimens were fabricated to compare the lap shear strength of an off-the-shelf adhesive (Plexus MA590) and two new adhesives developed by Arkema (Bostik SAF30 90 and Bostik SAF30 120). ISO standard 4587:2003 was used to standardize the testing method and sample fabrication. Lap shear specimens were made at 1mm, 3mm, and 10mm thicknesses. The Bostik adhesive lap shear test results were similar to Plexus for all thicknesses. (Section 5.6) Fiber-reinforced polymer (FRP) composites are typically used in high-performance applications (e.g., aerospace), and their expansion into high-volume industries (e.g. consumer automotive and wind turbine blade manufacturer or similar) is hindered by their cost and a lack of efficient manufacturing techniques. Monitoring the curing process of these composites during manufacturing can improve the efficiency of the process, and therefore reduce the manufacturing cost. Cure monitoring techniques were developed that use probabilistic estimation methods and surface temperature measurements made using infrared cameras. These techniques enable real-time monitoring of the infusion process to locate manufacturing flaws, and they can, potentially, estimate residual stresses in the part. Their commercialization will help facilitate expansion of FRP composites in high-volume industries. (Section 5.7) A 13-meter composite wind blade was produced with Elium® resin and Johns Manville fiberglass; this blade was made with VARTM processing similar to how megawatt-scale wind blades are currently manufactured, but no post-mold heating was used for this thermoplastic composite blade. The wind blade underwent full-scale validation for static loading (4-different load orientations) and flapwise fatigue loading to simulate 20-years of operational loads. The thermoplastic composite wind blade withstood the loading without any noted issues and performed similar to results from a previous full-scale validation to an equivalent epoxy composite wind blade produced with the same blade molds. (Section 5.8) A study was conducted to determine the feasibility of recycling composite wind turbine blade components fabricated with glass fiber reinforced Elium® thermoplastic resin. Dissolution, which is a process unique to thermoplastic matrices, allows recovery of both the polymer matrix and full-length glass fibers, while maintaining their stiffness and strength throughout the recovery process. The economics of recycling is favorable if 50% of the glass fiber is recovered and resold for a process of $\$$ 0.28/kg, and 90% of the resin is recovered and resold at a price of $\$$ 2.50/kg.(Section 10) Recommendations are outlined for commercializing thermoplastic resin for composite wind blade production, in addition to recommended areas for future research.
Accurate and fine-scale estimates of biogenic carbon fluxes are critical for measuring and monitoring the biosphere's responses and feedback to the climate system. Currently available data products from flux towers and model-intercomparison projects struggle to adequately represent spatiotemporal dynamics of surface biogenic carbon fluxes, and to quantify their uncertainties, which also are crucial to atmospheric inversion systems. To address these gaps, we introduce a new perturbed-parameter model ensemble with the CASA model to estimate surface biogenic carbon fluxes at monthly and 3-hourly scales for North America at ~500-m and 5-km resolutions. We first use the Extended Fourier Amplitude Sensitivity Testing to choose the three most sensitive parameters to be perturbed, maximum light use efficiency (E max ), optimal temperature of photosynthesis (T opt ), and temperature response of respiration (Q 10 ). The initial range for each parameter is broadly sampled for the L1 ensemble, but then we pruned E max with site-level primary productivity to derive an L2 ensemble with narrower uncertainty ranges. Ensembles are strongly correlated with site-level results at both monthly and 3-hourly scales, and the spread across L1/L2 ensemble members encompasses the range of AmeriFlux observations. Monthly variability in the L2 ensemble mean is 85% of the observed variability. The L2 ensemble outperforms diverse data products with the highest Taylor skill scores at diurnal to annual scales. Finally, the ensemble's seasonality agrees well with other models for most biome types and in high and middle latitudes, but inconsistencies are found in subtropical and tropical ecoregions and for annual totals over North America.
Quantitatively predicting the corrosion resistance of alloys is vital to the design and application of the next generation of superior corrosion resistant alloys (CRAs). Here, the Chloride Susceptibility Index (CSI), a scientifically based quantitative descriptor for repassivation tendency of CRAs, is proposed as a metric constructed through combining atomistically resolved information regarding adsorption and alloy surface composition with environmental conditions such as applied potential, temperature, pH, and Cl – concentration. Using this method, the response of CRAs to relevant environments can be predicted. CSI is determined by: (1) estimation of thermoequilibrium surface coverages of O and Cl by a Langmuir isotherm model using adsorption energies of different species (O, Cl, OH, H 2 O) obtained from DFT. (2) determination of CSI by integrating Cl surface coverage over a realistic window of applied potentials, thereby incorporating the effect of Cl surface coverage and electrochemical environment simultaneously. A series of Ni-Cr-X alloys was used as an example to develop and validate the method through studying the effect of different alloy solutes X on chloride resistance. Here, the trends and relations predicted by CSI are in qualitative agreements with experimental observations. Moreover, a quantitative correlation is found between CSI and the repassivation potential.
he Areal Averaged Albedo VAP yields areal averaged surface spectral albedo estimates from MFRSR measurements collected under fully overcast conditions via a simple one-line equation (Barnard et al., 2008), which links cloud optical depth, normalized cloud transmittance, asymmetry parameter, and areal averaged surface albedo under fully overcast conditions.